Невязка. Как ее уменьшить?

[TDL 0]

При решении задачи, невязка установилась на величинах 0,1 до 0,01 изменяясь периодически в зависимости от номера итерации. Т.е. решение и не разваливается, но и не сходится.

Что можно сделать в таком случае, для получения решения ?

[Iatsouk 3]

В какой системе работаете ?

По какому параметру невязка ?

В качестве примера действий:

Во Fluente при расчете течения я поставил на боковых границах условие pressure outlet. У меня через эти границы начала зачем-то ( 8))) ) втекать жидкость (предполагалось, что она будет только понемножку вытекать). Невязка в уравнении неразрывности не скакала, но стояла на отметке 0.1, что в моём представлении очень много.

Плюнул, задал на границах условие wall со скольжением (нулевое трение) - стало получше.

[TDL 0]

В моей понимании невязка как бы ее не определили должна зависит от параметров сетки и других параметров дискретизации (какого порядка дискретизации используется схема). На пример, можно ли утверждать, что с увеличением числа ячеек невязка уменьшается и стремится к 0 если число ячеек стремится к бесконечности? А как, влияет на невязку точность расчета (8 машинных знаков или 16). С какой точностью стоит считать, если минимальный объем ячейки 10e-17 при расчете методом контрольных объемов - (Real 8 или 16 а может 32 )?

Для определенности можно взять понятие невязки как во Fluent для Residuals for the Segregated Solver.

А у меня ситуация вообще не понятная для меня. Ведется расчет нестационарной задачи. При одних временах решение сходится, при других временах невязка стабилизируется и колеблется в течение всего времени итерации для данного времени. Не сходится критерий по уравнению неразрывности (самый первый критерий) по другим параметрам сходится. (шаг по времени 1e-5, использовался адаптивный шаг по времени по критерию который использует Fluent, правда данный шаг Fluent не увеличивал, реагирующая смесь ).

[Miki 3]

так временной шаг является тоже одним из параметров кроме дискретизации и сетки, влияющем на невязки. тем более в нестационаре.

Вас чем то не устроил тот шаг при котором все сошлось ? значит в процессе есть что то, что требует такого шага, выведите невязки в результаты и посмотрите где они, обычно в разных тупиковых областях и вихрях, и после этого уже думайте что с этим делать

[Miki 3]

имхо, не стоит придавать обьему элемента столь важную роль, уменьшение обьема далеко не всегда двигает решение в сторону сходимости, так же как и очень маленькие временные шаги.

И уж тем более обьем элемента стараться тащить до 1е-24 м^3 если гекса то грань 1/100000 миллиметра !! или Вы там DNS занимаетесь ?? но даже если не занимаетесь то интересно, что у Вас там за железо стоит, для столь мелких сеток ??

вобщем универсальный метод-- выводите невязки в результаты и смотрите

[TDL 0]

Ну вывел я в fail. Допустим невязка по уравнению неразрывности скачет 0,1 да 0,01 стабильно с некоторым периодом в зависимости от итерации. Так что сделать, чтобы уменьшить невязку. Вы можете посоветовать, кроме того как

вобщем универсальный метод-- выводите невязки в результаты и смотрите

Какие нужно сделать стандарные (т.е. какакие народ сделал шаги) шаги

1. Увеличили число ячеек.

2.

3.

4.

5.

..........

[TBI 0]

На сколько я понимаю, при измельчении сетки по координате алгоритм обычно хуже себя ведет в плане сходимости... Из моего опыта (работаю с Phoenics) могу сказать, что на мелкой сетке невязки обычно больше, чем на исходной. Измельчение шага по времени положительно влияет на сходимость - в частности в phoenics предусмотрены параметры релаксации, т.н. псевдо-временной шаг (false DT) и стационарные задачи решаются в "псевдо нестационаре". Таким образом, попробуйте уменьшить (на сколько это возможно) количество ячеек сетки по координате, шаг по времени можете попробовать также уменьшить (опять же на сколько это возможно).

[TDL 0]

1) А можно ли утверждать, что если невязка стабилизируется на одном значении, то ошибка расчета не превосходит exp(значение невязки*время)?

2) сказал, что он экспериментально наблюдал

]

Из моего опыта (работаю с Phoenics) могу сказать, что на мелкой сетке невязки обычно больше, чем на исходной.

<{POST_SNAPBACK}>

интересно аналитически это доказано (где можно посмотреть методику). И с чем это связано (с точность расчетов или с методом реализации дискретизации).

[TBI 0]

интересно аналитически это доказано (где можно посмотреть методику). И с чем это связано (с точность расчетов или с методом реализации дискретизации).

<{POST_SNAPBACK}>

Честно говоря, не встречал статей относительно этого факта. Как объяснить - пока тоже не знаю.