Комментарии 18
Студенческая работа?
Ага
Но меня интересует мнение профессионалов в этой области...
Профи в этой области на хабре что забыли ;) Они варятся в инженерии связанной с аэро/гидро динамикой. Это военка (этих сразу посадят), и производители авто, самолётов, f1 и всякие автоспорты, этих просто уволят, как минимум. ;)
Например, интересует мнение @kbtsiberkin
За что посадят/уволят? За мнение о стандартной студенческой задаче? :)
Знаете, на одной конференции слушал сообщение об «обтекании конического тела гиперзвуковым потоком», и ничего. Была рассказана красивая и интересная фундаментальная научная задача.
А автор живет в рф, и всё ещё на свободе? Или там была чисто гипотетическая "практически не применимая в инженерии дальних скоростных полётов конусных тел" задача ;) Если это было не публично на кухне, тогда не считается.
Среди прочего речь скорее шла про применимость пакета OpenFOAM при моделировании обтекания указанного предмета. Конференция «Вихри и волны в сложных средах» (Москва, Институт проблем механики РАН), год в районе или 2014–2015, или 2019, когда я сам туда ездил.
Докладчик, отвечая на вопросы, в один момент вместо «тело» неожиданно сказал «изделие», после чего дискуссия как‑то сама собой затухла. Но точно всё прошло спокойно и беспроблемно. Не нужно преувеличивать масштабы.
Правильно ли я посчитал и визуализировал магнитные поля соленоида и постоянных магнитов?
У меня только что закончился отпуск, а впереди работа на полную ставку профессора с руководством образовательной программой (это локальный аналог заведования кафедрой), трое аспирантов, один магистрант, и ожидается по одному студенту на 3 и 4 курсах бакалавриата, может ещё приклеятся.
Думаете, я буду вот так всё подряд просматривать и консультировать на разных ресурсах?
Если автору интересно развиваться в этой области, то дальше можно ознакомиться с другими группами методов вычислительной гидродинамики. Часть из них не решает уравнение Навье-Стокса напрямую (например, метод решеточных уравнений Больцмана), позволяют работать напрямую с воксельной геометрией и дают возможность иначе взглянуть на уже привычные гидродинамические процессы. И с точки зрения реализации в виде кода будет интереснее)
Зачем‑то меня в комментарии призвали. Я последнюю пока свою гидродинамическую работу опубликовал лет пять назад, переключившись на магнетизм. А работы с конечно‑разностными моделями и того давнее.
Имею сказать, что центральные конечные разности — плохая идея. Они при некоторых условиях вообще безусловно неустойчивые, здесь, полагаю, чисто за счёт нелинейности задачи какие‑то решения выживают. Устойчивые схемы реализуются разного рода разностями против потока, которые, правда, отягощены высокой схемной (ненастоящей) вязкостью — поэтому и стабильны. Либо иными хитрыми подходами.
Сетка 200×5 узлов — это трындец. Отношение размера ячеек больше 2:1 вообще не рекомендуется, погрешность счёта же пропорциональна шагу или квадрату шага. В итоге, когда по одной оси погрешность в 40 раз больше, чем по другой (при первой степени аппроксимации; у вас вторая степень - итого в 1600 раз больше!) расчёт вообще смысл теряет. В рассматриваемой задаче (так понимаю из кода — увы, не приведены в тексте граничные условия) движется правая стенка? Т.е. стандартный тест «течение в каверне». Результат даже близко на него не похож.
Течение в каверне с ResearchGate. Подвижная верхняя стенка, расчёт в OpenFOAM.

Советую погрузиться в упомянутую монографию Роуча (третья глава — что ни параграф, то та или иная рабочая разностная схема), труды А.А. Самарского (и его ученика П.Н. Вабищевича) или одного из моих учителей Е.Л. Тарунина.
Новизны же в здешней работе, достойной публикации (об этом меня автор спросил в личке), увы, нет. Как сказано выше — студенческая работа. Причём не самая толковая.
А правильно ли здесь примененять проекционный метод?

Гидродинамика на python: Пишем CFD симуляцию плоского течения несжимаемой жидкости