这为不可压缩的Navier Stokes方程在交错网格系统上的数值求解提供了一个MATLAB例子代码。(日本語ドキュメントもあります)
第1部分:从空腔流动开始
数值格式保持原始;粘性项的显式处理(解在低雷诺数时发散),时间积分为欧拉。
まずは単純な手法でキャビティ流れのシミュレーションを実施します。
第二部分:粘性项的隐式格式
粘性方面的暗示疗法的实施,即曲柄尼科尔森方法。对于用于非线性项,亚当斯-Bashforth,和3个更好的稳定性步骤-龙格 - 库塔也被实现。
拡散項に対して陰解法を実装しました。対流項へアダムス・バッシュフォースを使用したもの,3段階のルンゲクッタ法の2つの時間発展を実装しています。
第3部分:隐式方法的性能比较
对粘性项的隐式处理使得离散化后的亥姆霍兹方程在每一时刻都可以求解。我们比较了五种方法的性能。
拡散項に対して陰解法を実装すると離散ヘルムホルツ方程式を毎時間ステップ解く必要が出てきます。直接法,反復法,離散サイン変換を使用した方法などの処理速度を検証します。
下一个来:
该计划允许任意边界条件,以进行更有趣的模拟。
道夫(2020)。CFD101: 2D Lid驱动空腔流GitHub (https://www.github.com/mathworks/2D-Lid-Driven-Cavity-Flow-Incompressible-Navier-Stokes-Solver)。检索。
1.2.0 | 增加了第3部分:隐式方法的性能比较 |
|
1.1.1 | 更新描述 |
|
1.1.0 | 增加了第2部分文档。Crank-Nicolson和Adams-Bashforth,三步龙格-库塔实现。 |
|
1.0.4 | 更新solvePoissonEquation_direct。m,通过使用持久变量消除重新创建矩阵 |
|
1.0.3 | 更新GitHub URL |
|
1.0.2中 | 更新缩略图GIF |
|
1.0.1 | 改变了标题 |
健太(查看配置文件)
E. Cheynet(查看配置文件)