2 d盖子空泡流驱动的

这对lid-driven回购MATLAB示例代码提供了一个空泡流,数值求解不可压缩纳维斯托克斯方程使用一个简单的二阶有限差分格式在交错网格系统。

(左:Re = 100,右:Re = 10000)

箭头表示速度场,等高线表示它的大小。

第1部分:空泡流的开始

• 点击在这里详细的英文文档。
• 日本語のドキュメントはこちらから

保存原始数值方案;粘性项的明确的治疗(解决方案在低雷诺数发散),集成是欧拉和时间。

まずは単純な手法でキャビティ流れのシミュレーションを実施します。

第2部分:隐格式的粘性

• 点击在这里详细的英文文档。
• 日本語のドキュメントはこちらから

隐式实现治疗粘性项,即Crank-Nicolson方法。更好的稳定性的非线性项,Adams-Bashforth 3 steps-Runge-Kutta也实现。

拡散項に対して陰解法を実装しました。対流項へアダムス・バッシュフォースを使用したもの,3段階のルンゲクッタ法の2つの時間発展を実装しています。

第3部分:隐式方法的性能比较

• 点击在这里详细的英文文档。
• 日本語のドキュメントはこちらから

隐式治疗粘性项结果在解决离散亥姆霍兹方程在每一个时间步。我们比较五种方法的性能。

拡散項に対して陰解法を実装すると離散ヘルムホルツ方程式を毎時間ステップ解く必要が出てきます。直接法,反復法,離散サイン変換を使用した方法などの処理速度を検証します。

第4部分:数值计算的验证方案

• 点击[这]用英语详细文档。(没准备好)
• 日本語のドキュメントはこちらから

空间和时间收敛性测试的结果。融合运行测试使用的方法制造解决方案,navier - stokes方程被迫这样规定的解决方案将是一个时间函数。万博 尤文图斯

ある外力項を加えたn - s方程式の数値解と解析解と比較することで,時間積分の精度(1次精度)と空間微分の精度(2次精度)を確認します。

下一个来

这个计划是让任意边界条件模拟更多的乐趣。

环境

• MATLAB R2019b
• 信号处理工具箱求解泊松eqn如果使用dct。

待办事项

1. 实现隐式处理粘性项
2. 实现crank-Nicolson非线性项
3. 域内允许的障碍
4. 允许流入从墙上
5. 3 d

- - -

版权(c) 2020年,MathWorks公司。