迎风差分格式
風上差分スキーム
理论与物理
概述
老师！今天要讲迎风格式对吧？那是什么东西呢？
一阶迎风、二阶迎风、TVD限制函数。数值耗散与稳定性的权衡。
控制方程
$$ \phi_f=\phi_U+\psi(r)\frac{\phi_D-\phi_U}{2} $$
$$ r=\frac{\phi_U-\phi_{UU}}{\phi_D-\phi_U} $$
离散化方法
这个方程，在计算机上实际是怎么求解的呢？
要使用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，构建整体刚度方程。
矩阵求解算法
矩阵求解算法，具体是指什么呢？
通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解线性方程组。对于大规模问题，带预处理的迭代法非常有效。
| 求解器 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小～中规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小～中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说在有限元法这部分偷懒的话，后面会吃苦头对吧。我铭记在心！
商用工具中的实现
那么，要用什么软件来做迎风格式呢？
| 工具名称 | 开发商/现状 | 主要文件格式 |
|---|---|---|
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Ansys CFX | Ansys Inc. | .cfx, .def, .res, .ccl |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| OpenFOAM | 开源（OpenCFD/ESI、OpenFOAM Foundation） | 字典文件（blockMeshDict等）, .foam |
供应商谱系与产品整合历程
各个软件的诞生过程，是不是挺有戏剧性的？
Ansys Fluent
接下来是Ansys Fluent的内容对吧。讲什么呢？
Ansys CFX
请讲讲「Ansys CFX」！
由AEA Technology (UK) 开发的CFX。2003年被Ansys收购。以耦合型求解器为特点。
当前所属: Ansys Inc.
听到这里，终于明白开发为什么重要了！
Simcenter STAR-CCM+
接下来是Simcenter STAR的内容对吧。讲什么呢？
哇～，开发的故事，超级有趣！请再多讲一些。
文件格式与互操作性
在不同求解器之间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、载荷与边界条件的表达差异。特别是高阶单元和特殊单元（如粘聚单元、用户自定义单元等），在求解器之间往往无法直接转换。
原来如此…格式看起来简单，实际上内涵很深呢。
实务注意事项
有没有教科书上不写的“现场智慧”之类的东西？
网格收敛性确认、边界条件合理性验证、材料参数敏感性分析非常重要。
- 网格依赖性验证: 至少用3个级别的网格密度确认收敛性
- 边界条件合理性: 设置物理上有意义的约束条件
- 结果验证: 与理论解、实验数据、已知基准问题进行比较
哎呀，迎风格式真是深奥啊… 不过多亏了老师的讲解
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告