我如何写出一个简单的刚性物体的二维 Navier-Stokes 方程,该物体浸入流动中并由于重力而自由下落?我试图查看由流体-固体相互作用引起的涡度。
我正在使用 Matlab 进行模拟。
我如何写出一个简单的刚性物体的二维 Navier-Stokes 方程,该物体浸入流动中并由于重力而自由下落?我试图查看由流体-固体相互作用引起的涡度。
我正在使用 Matlab 进行模拟。
挑战不是可视化,而是编写 Navier-Stokes 求解器。这是因为这不仅是一个难以解决的方程,而且您现在还使用了一个随时间变化的域,最重要的是,它的变化方式(即椭圆如何下落)取决于流动。
处理这个问题的科学计算领域被称为“流固耦合”。对于您的特定情况,人们使用的典型方法称为“浸没边界法”,您会发现相当多的出版物或多或少地涉及您关心的问题。一个合理的起点是 Luca Heltai ( https://scholar.google.it/citations?user=9TUDvmgAAAAJ&hl=en ) 或他本人 Charles Peskin ( https://scholar.google.com/scholar?hl ) 的论文=en&as_sdt=0%2C6&q=charles+peskin&btnG=&oq=charles+pesk )。您可能还对实现示例感兴趣,例如本教程程序:https ://www.dealii.org/developer/doxygen/deal.II/step_60.html