ASME V&V 20: CFD的验证与确认
ASME V&V 20的理论基础
概述
老师！今天要讲的是ASME V&V 20: CFD的验证与确认吗？这是什么东西呢？
基于ASME V&V 20规格的CFD结果不确定性评估。数值不确定性、输入不确定性、实验不确定性的整合。
等等，等等，基于规格这样的话，意思是这种情况也能用吗？
控制方程
嗯……的验证与确认看似简单，但实际上涵盖很深啊。
离散化手法
这个方程在计算机上具体怎么求解呢？
使用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，构建全局刚度方程。
矩阵求解算法
矩阵求解算法具体是什么意思呢？
直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程。大规模问题中预处理迭代法有效。
| 求解方法 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小~中规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小~中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模、非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说，在有限元法这个阶段出问题的话，后面会很吃苦头啊。我牢记在心！
商用工具中的实现
那么，要进行ASME V&V 20: CFD的验证与确认，需要什么样的软件呢？
| 工具名称 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| MSC Nastran / NX Nastran | MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software） | .bdf, .dat, .f06, .op2, .pch |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| Ansys Mechanical (原ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
| OpenFOAM | 开源（OpenCFD/ESI、OpenFOAM Foundation） | 字典文件（blockMeshDict等）, .foam |
供应商系统与产品整合的发展过程
各个软件的成长历程是不是，有点戏剧性呢？
MSC Nastran / NX Nastran
接下来是MSC Nastran的话题吗。什么内容呢？
NASA结构分析（NASTRAN）在1960年代开发。MSC Software进行商业化，后来UGS（现在是Siemens）派生出NX Nastran。MSC在2017年被Hexagon AB收购。
现在的所属: MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）
Abaqus FEA (SIMULIA)
Abaqus FEA是什么意思呢？
1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购并整合到SIMULIA品牌。
现在的所属: Dassault Systèmes SIMULIA
等等，等等，结构解析的话，也就是说这样的情况也能用吗？
Ansys Mechanical (原ANSYS Structural)
关于"Ansys Mechanical"，能讲讲吗！
1970年由Swanson Analysis Systems Inc. (SASI) 开发。以APDL（Ansys Parametric Design Language）为基础。
现在的所属: Ansys Inc.
哦哦~，结构解析的话，超有意思啦！ 再多讲讲吧。
文件格式与相互操作性
不同软件间要交换数据的时候，有什么需要注意的吗？
| 格式 | 扩展名 | 类型 | 概述 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | ISO 10303标准的3D CAD数据交换格式。形状+PMI支持。 |
| IGES | .igs/.iges | 中立CAD | 早期CAD数据交换规格。曲面数据兼容存在问题。向STEP迁移进行中。 |
| VTK | .vtk/.vtu | 可视化 | Visualization Toolkit格式。用于ParaView等。 |
在不同求解器间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载·边界条件的表现差异。特别是高阶单元或特殊单元（内聚单元、用户定义单元等）在求解器间往往无法直接转换。
原来……格式看似简单，实际涵盖很深啊。
实务注意事项
教科书里没有的"现场智慧"有吗？
网格收敛性确认、边界条件妥当性验证、材料参数的灵敏度分析特别重要啊。
哎呀，ASME V&V 20: CFD的验证与确认涵盖这么深……但听了老师的讲解，整理得差不多了！
嗯，进展不错啊！ 实际操作是最好的学习，有什么不懂的随时问。
ASME V&V 20的数值计算手法
数值手法详情
具体用什么算法来求解ASME V&V 20: CFD的验证与确认呢？
离散化的表述
利用形状函数 $N_i$ 逼近未知量:
用公式表示是这样的。
基础方程的离散形式
用公式表示是这样的。
嗯……只是公式的话，不太能接受……什么意思呢？
连续体的支配方程离散化后，得到以下代数方程系:
其中$[K]$是全局刚度矩阵（或等价的系统矩阵），$\{u\}$是未知节点变量向量，$\{F\}$是外力向量啊。
啊，是这样！ 连续体的支配方程离散化是这样的原理啊。
单元技术
"单元技术"听过，但可能理解不到位……
| 单元类型 | 阶数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中~高 | 中 |
积分体系
积分体系具体是什么意思呢？
这样听了，总算明白为什么单元类型这么重要了！
收敛性和稳定性
没有收敛的话，首先要检查什么？
收敛速度: 二阶单元以 $O(h^2)$ 阶数误差减少（光滑解情况）
原来……网格细化看似简单，其实涵盖很深啊。
求解器设置建议
具体用什么算法来求解ASME V&V 20: CFD的验证与确认呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理手法 | ILU(0) or AMG | 根据问题规模 |
| 最大反复次数 | 1000 | 非收敛时需见设置 |
| 内存模式 | In-core | 尽可能 |
低阶单元
计算成本低且实现简洁，但精度有限。粗网格时可能产生大误差。
高阶单元
相同网格下达到更高精度。计算成本增加，但所需单元数往往减少。
牛顿-拉夫森法
非线性问题的标准方法。收敛半径内有二次收敛。$||R|| < \epsilon$ 时判定收敛。
时间积分
显式法: 条件稳定（CFL条件）。隐式法: 无条件稳定，但每步需求解联立方程。
验证数据的可视化
理论值与计算值的比对定量呈现。误差5%以内为合格基准。
| 评估项目 | 理论值/参考值 | 计算值 | 相对误差 [%] | 判定 |
|---|---|---|---|---|
| 最大位移 | 1.000 | 0.998 | 0.20 | PASS |
| 最大应力 | 1.000 | 1.015 | 1.50 | PASS |
| 固有振动数（1阶） | 1.000 | 0.997 | 0.30 | PASS |
| 反力合计 | 1.000 | 1.001 | 0.10 | PASS |
| 能量守恒 | 1.000 | 0.999 | 0.10 | PASS |
判定基准: 相对误差 < 1%: ■ 优良、1~5%: ■ 允许、> 5%: ■ 需检讨
ASME V&V 20的实务应用
实践指南
老师，关于"实践指南"请讲讲！
ASME V&V 20: CFD的验证与确认的实际解析流程与注意点来讲解。
老师讲解通俗易懂！ 的验证与确认的模糊不清的地方解开了。
解析流程
从第一步开始讲讲吧！ 怎么开始呢？
1. 预处理 (Pre-processing)
- CAD数据的输入和形状简化
- 材料特性的定义
- 网格生成（单元类型、大小的决定）
- 边界条件和荷载条件的设置
2. 求解 (Solving)
- 求解器设置（解法、收敛基准、输出控制）
- 作业投入和计算执行
- 收敛监视
3. 后处理 (Post-processing)
- 结果的可视化（位移、应力、其他物理量）
- 结果的检验和妥当性确认
- 报告作成
网格生成的最佳实践
网格的好坏怎么判断呢？
单元品质指标
关于"单元品质指标"请讲讲！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 长宽比 | 1.0 | < 5.0 | 精度降低 |
| 雅可比比 | 1.0 | > 0.3 | 单元退化 |
| 翘曲 | 0° | < 15° | 精度降低 |
| 倾斜度 | 0° | < 45° | 收敛性恶化 |
| 锥形比 | 0 | < 0.5 | 精度降低 |
网格密度的决定
网格密度的决定具体是什么意思呢？
边界条件设置指南
边界条件，听说这里搞错全完蛋啊……
啊，是这样！ 过拘束注意原来是这样的原理啊。
商用工具分别实施步骤
有这么多软件，对吧？ 各自的特点讲讲吧！
| 工具名称 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| MSC Nastran / NX Nastran | MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software） | .bdf, .dat, .f06, .op2, .pch |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| Ansys Mechanical (原ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
| OpenFOAM | 开源（OpenCFD/ESI、OpenFOAM Foundation） | 字典文件（blockMeshDict等）, .foam |
MSC Nastran / NX Nastran
接下来是MSC Nastran的话题吗。什么内容呢？
NASA结构分析（NASTRAN）在1960年代开发。MSC Software进行商业化，后来UGS（现在是Siemens）派生出NX Nastran。MSC在2017年被Hexagon AB收购。
现在的所属: MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）
Abaqus FEA (SIMULIA)
Abaqus FEA是什么意思呢？
1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购并整合到SIMULIA品牌。
现在的所属: Dassault Systèmes SIMULIA
老师讲解通俗易懂！ 工具名的模糊不清的地方解开了。
常见失败与对策
新手容易犯的失误类型有吗？ 事先知道想！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、不当边界条件 | 网格改善、拘束条件见设置 |
| 应力异常大 | 应力特异点、网格依赖 | 特异点回避、局部网格细化 |
| 位移非现实的 | 材料常数错误、单位系统不一致 | 输入数据确认 |
| 计算时间过长 | 不必要的细化、非高效解法 | 网格最优化、并行计算 |
品质保证检查清单
教科书里没有的"现场智慧"有吗？
哎呀，ASME V&V 20: CFD的验证与确认涵盖这么深……但听了老师的讲解，整理得差不多了！
嗯，进展不错啊！ 实际操作是最好的学习，有什么不懂的随时问。
ASME V&V 20的ASME V&V 20软件比较
商用工具比较
有这么多软件，对吧？ 各自的特点讲讲吧！
ASME V&V 20: CFD的验证与确认支持的主要商业CAE工具的功能比对与各产品的历史背景详述。
支持工具列表
那么，要进行ASME V&V 20: CFD的验证与确认，需要什么样的软件呢？
| 工具名称 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| MSC Nastran / NX Nastran | MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software） | .bdf, .dat, .f06, .op2, .pch |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| Ansys Mechanical (原ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
| OpenFOAM | 开源（OpenCFD/ESI、OpenFOAM Foundation） | 字典文件（blockMeshDict等）, .foam |
MSC Nastran / NX Nastran
接下来是MSC Nastran的话题吗。什么内容呢？
NASA结构分析（NASTRAN）在1960年代开发。MSC Software进行商业化，后来UGS（现在是Siemens）派生出NX Nastran。MSC在2017年被Hexagon AB收购。
现在的所属: MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）
Abaqus FEA (SIMULIA)
Abaqus FEA是什么意思呢？
1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购并整合到SIMULIA品牌。
现在的所属: Dassault Systèmes SIMULIA
等等，等等，结构解析的话，也就是说这样的情况也能用吗？
Ansys Mechanical (原ANSYS Structural)
关于"Ansys Mechanical"，能讲讲吗！
1970年由Swanson Analysis Systems Inc. (SASI) 开发。以APDL（Ansys Parametric Design Language）为基础。
现在的所属: Ansys Inc.
Ansys Fluent
接下来是Ansys Fluent的话题吗。什么内容呢？
Fluent Inc.开发。2006年被Ansys收购。非结构化网格基础的通用CFD求解器。
现在的所属: Ansys Inc.
哦哦~，结构解析的话，超有意思啦！ 再多讲讲吧。
功能比对矩阵
预算和时间有限，成本效益最强的是哪个？
| 功能 | Nastran | Abaqus | Ansys Mechanical | Fluent |
|---|---|---|---|---|
| 基本功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU支持 | △ | △ | △ | ○ |
转换时的风险
转换时的风险具体是什么意思呢？
啊，是这样！ 不同工具间的数据转换就是这样的原理啊。
许可形式
"许可形式"听过，但可能理解不到位……
| 工具 | 许可 | 特点 |
|---|---|---|
| 商业FEA | 节点锁定/浮动 | 高额但附官方支持 |