Сингулярность в потоке идеальной жидкости. Течение за уступом.

[Jesse 979]

Продолжаю тему сингулярности. На этот раз цель получить сингулярность, т.е. отсутствие сеточной сходимости около тупого угла в потоке идеальной жидкости. 
Ожидал получить что-то типа такого

Считал в floefd в 2D. Задал для воды вязкость 10e-40. На всякий включил гравитацию.
Сверху симметрия. Самая плотная сетка (решил не борщить, для ускорения процесса - много расчетиков прогонял)

Пробовал обычный нестационарный k-e решатель в первом случае и ламинарный решатель (с вырубленной турбулентной вязкостью) во втором случае.

На эпюре скорости скачков вообще не видать и на разных сетках результаты похожи 

 

а вот эпюра завихренности оказалась многообещающей: тут для разных сеток с кол-вом ячеек в ~2 раза больше в каждом расчёте результаты получились ~0, 57 и 220!

Вообще результаты довольно странные... можно начать с того, что завихренность бывает только в вязкой жидкости и отсутствует в идеальной (если её там нет в "начальных условиях")
Хотя возможно физичные результаты получились как раз на грубой сетке (там где завихренность == 0), а на плотной сетке начала проявлять себя численные особенности всякие.
Вообще, ожидал получить сингулярность по полю скорости
Мб такие задачи надо только в DNS решать?

Вообще кто-то что-то подобное делал?)

 

[Foksmen 17]

18.02.2023 в 22:04, Jesse сказал:

Вообще результаты довольно странные... можно начать с того, что завихренность бывает только в вязкой жидкости и отсутствует в идеальной (если её там нет в "начальных условиях")

 

Есть такая вязкость "схемная".  Она зависит от порядка точности схемы. 

[Jesse 979]

54 минуты назад, Foksmen сказал:

Есть такая вязкость "схемная".  Она зависит от порядка точности схемы. 

Вы про вязкость как численные ошибки?

https://en.m.wikipedia.org/wiki/Numerical_diffusion

[Foksmen 17]

Да. Погрешность апроксимации 

[Jesse 979]

10 часов назад, Foksmen сказал:

Да. Погрешность апроксимации 

да, скорей всего из-за этого и появилась завихренность, которой не должно быть в жидкости с заданной вязкостью 10^-40