CFX前处理——求解器控制.pptVIP

1/25/05 ANSYS, Inc. Propietary CFX前处理求解器控制 讲座 5 求解器控制 在求解过程中控制求解器的行为 差分格式 收敛控制 收敛准则 计算时间控制 方程组设置 高级选项 如何定义求解器设置 如前所示… 差分格式 差分格式 CFX中的插分格式是以下面的形式实现的: 其中 是迎风节点上的值, 是从上一个迎风节点到ip 的矢量， 是临界节点的平均梯度 时间步 ANSYS CFX 采用了所谓的伪瞬态运算法则. 它在模拟迭代的过程中需要时间步. 在稳态模拟计算中，时间步提供了非线性方程的松弛因子 一个稳态的模拟是从初始猜测值到稳态情况下的瞬态发展 收敛的结果相对于所使用的时间步独立 时间步选择 为了获得成功的收敛性, 时间步的选择在其中起着非常重要的角色. 如果时间步过大, 收敛性就会变得有弹性，或者会导致求解器的失败. 如果时间步过小, 收敛性将会非常的慢，结果可能不会十分准确 对于水平对流为主的流动来说, 流体在计算域内停留的时间分数经常是确定时间步的好参考 时间步选择 通常选择长度尺度/速度尺度的1/3作为时间步 可能在迭代的最初几步或是复杂的物理现象，跨声速流动等模拟时，需要较小的时间步 对于旋转机械来讲, 1/? 是个好的选择 对于浮力驱动流, 时间步应该基于重力，热膨胀率，温度差和长度尺度的函数 时间尺度控制 时间尺度控制 自动时间尺度 求解器基于边界条件和初始条件或是求解问题和计算域长度尺度计算一个固定的时间步 对于计算域的长度尺度或是特定的值，使用一个保守的或是大胆的估计. 时间步每隔几个迭代步就会重新计算或是更新 可以设置最大的时间尺度来提供上限 自动时间尺度趋于保守 物理时间尺度 指定时间步. 通常为一个常数，但有时也可以通过CEL编辑一个变量 可以比自动时间尺度更快的收敛 时间尺度控制 当地时间尺度系数 时间尺度在计算域内变化 在当地时间尺度系数中输入的值是默认的当地基于单元的时间尺度(5)的乘数 可以加速收敛，特别是当当地的速度尺度变化非常剧烈时 当网格的大小变化过于剧烈时，收敛性差 对于完全收敛的结果不要使用; 总是在结束时伴随着一个常数时间步 收敛方案 收敛方案残差 ANSYS CFX 使用一个多重网格（MG）加速的不完全降低上层（ILU）因数分解技术 (迭代求解器) 来求解分离系统的线性方程形式: [A] [Φ] = [b] 其中 [A] 是系数矩阵, [Φ] 是求解矢量并且 [b] 方程右边 方程是迭代求解的: 新的结果是从一个老的结果的基础上修正得到的: Φn+1 = Φn + Φ’ 余项被定义为Rn = b – AΦn 每个变量最初的残差会被规则化，这是为了求解监控和收敛准则的要求 RΦ~ = RΦ / (ap ΔΦ) 其中 ap 是控制体积系数 ，ΔΦ代表了在计算域中一个变化的范围 收敛方案残差 残差类型 MAX: 收敛基于计算域内最大的残差 RMS: 收敛基于计算域内所有控制体积的平均残差 平均平方根 = 残差目标 对于一种合理的收敛情况来讲，最大的残差应该小于1.0E-3, 平均残差至少要小于1.0E-4. 目标残差要基于精确度的需要. 可能为了获得更高的精度需要更小的残差目标值 收敛方案守恒目标 守恒目标 一个额外的收敛指标应该在求解结束之前达到 设置这个选项有利于结果 使用1% (0.01) 或更小: 在求解其求解结束前 流入 – 流出 1% 逝去时间控制 逝去时间控制 逝去时间控制 设置最大的计算时间 固体计算域时间尺度控制 固体计算域时间尺度控制 这个选项只有当存在固体计算域是才可用 固体计算域时间尺度应该选择为远远大于流体计算域中的时间尺度 (10 或是 100倍的流体计算域中的时间尺度通常是合法的) 能量方程在固体计算域中通常是很稳定的 固体计算域中时间尺度要远远大于流体计算域中时间尺度 流体计算域中时间尺度可以有下面的公式确定: L/V 固体计算域中时间尺度自动计算为: L2/α 其中: α=λ/ρcp 高级设置 方程组设置 差分格式和时间步可以针对每个方程进行控制 不需要经常使用 高级选项 高级求解器控制选