[2018年最新整理]5求解器的设置.pptVIP

大纲 使用求解器（大致求解过程） 设置求解器参数 收敛 定义 监视 稳定性 加速收敛 准确度 网格独立性 网格自适应 非定常流动模型 非定常流动问题设置 非定常流动模型选项 摘要 附录 非定常流动模型 解算在每个时间步内都迭代到收敛，然后前进到下一步 解的初始化定义了初始条件而且这些条件必须是真实的 非迭代时间前进法对更快的计算机时间适用（看附录中的细节） 对基于压力的解算器： 时间步长, ?t, 在 Iterate 面板中设置 ?t 必须足够小 用来解决依赖于时间的特征; 保证在每个时间步的最大迭代步数中能够达到收敛 合适的时间步长可以由下估计得出： 时间步长估算也能够被选择为了能够解算得出非定常流动的特征（比如流动在一个已知的波动周期内） 使用零时间步数来定义及时不前进的迭代 PISO格式可能对许多非定常流动的加速收敛都有帮助 非定常流动模型选项 自适应时间步 通过分析局部舍入误差来自动调节时间步长 可以通过UDF来定制 时间统计平均 对LES紊流模型特别适用 如果值得的话，在迭代之前应该设置运动（对于可以看到的流动） 对于基于密度的求解器，柯朗数定义了： 对于密度显式求解器的全局时间步长 对于密度隐式求解器的虚的时间步长 真实的时间步长依然在Iterate 面板中定义 小结 压力模式和密度模式的求解器的求解过程是相同的 一直算到一个收敛的结果 获得一个二阶的解（建议） 优化网格并且重新计算直到获得一个独立网格的解 所有的求解器提供工具用来判定和改善收敛和确保稳定 所有的求解器提供工具来检查和改善精确度 解的精确性依赖于你所选择的物理模型和定义的边界条件是否合适 附录 背景 有限体积法 显式对比隐式 独立对比耦合 瞬态解 NITA和ITA格式的流程图 The Finite Volume Method有限体积法 FLUENT 求解器基于有限体积法. 计算域离散为控制体的有限集合 全部的质量，动量，能量等输运方程离散后应用于每一个单元体 解出所有方程用来显示流场 被一系列控制体积离散的管道流体区域 控制体积 Unsteady Convection Diffusion Generation 方程 变量 连续性 1 X 动量 u Y动量 v Z动量 w 能量 h 有限体积法 每个输运方程离散为代数形式，对单元体 P来说, 离散后的方程需要单元体中心和表面上的信息 流场数据（材料属性，速度，等等）存储在单元体中心 面上信息由当地或者相邻的单元体信息插值得到 离散精确度依赖于 截断误差. 离散后的方程可以简单表示为 对计算域中的每一个控制体适用的方程构成一个方程组 face f adjacent cells, nb cell p 线性化 方程组反复迭代求解 系数 ap 和 anb 是解的变量 系数矩阵 (非线性耦合). 系数从之前的迭代计算中的解算变量中写入 线性化除去对系数?的相关性. 退耦过程除去其它求解变量的系数相关性 系数随着每次外迭代实时更新 对于所给的内迭代，系数是常数（冻结） ?p 既能够显式解出也能够隐式解出. 显式、隐式解对比 关于 ?nb 的假设 显式线性化 从每一个单元体上的已知数和它们之间的关系计算出未知数(?nb 假定从之前的迭代中已知). ?p 然后用龙格－库塔法显式解出 隐式线性化 ?p 和?nb 假定为未知的，用线性方程方法来解 方程隐式线性化趋向于对稳定性要求较低 方程组在解的同时要使用二次迭代循环（比如，指向高斯－赛德尔） 基于压力的对比基于密度的求解器 基于压力的求解器 如果假设一个方程中只有一个未知数，那么方程组就能够不依赖其它变量的解而解出了 简单输入，每个控制方程都可以独立解出 在这种情况下，系数ap 和anb 是标量值 基于密度的求解器 如果在每个方程中不止一个变量是未知的，而且每个变量是由它自身的输运方程所定义的，那么这个方程组就是耦合在一起的 在这种情况下, 系数 ap 和 anb 是 Neq × Neq 矩阵. ? 是一个因变量矢量, {p, u, v, w, T, Y}T 基于压力的求解器 在基于压力的求解器中，每个方程都能够独立解出 连续性方程把压力修正方程的形式作为SIMPLE算法的一部分 欠松弛因子包含在离散化的方程中 包括它为了改善迭代过程的稳定性 显式欠松弛因子α ，限制了变量从一个迭代步到下个迭代步之间的改变 更新属性 解动量方程 (u, v, w 速度) 解压力修正 (连续性)方程更新压力场和面的质量流率 求解能量、组分、湍能和其他标量方程 是 否 收敛? 停 基于密度的求解器 连续性，动量，能量和多物质在密度求解器中同时解出 修正后的方程可解出可压流和不可压流 衰减项经常包括其中 随着时间推进和瞬时现象趋近于零，稳态解就