深入浅出Fluent5.ppt

深入浅出Fluent（五） 求解器 适应性网格 本讲概要 求解器 Fluent求解概述 求解过程 选择求解器 初始化计算 余元（残差） 收敛性 计算精度 适应性网格 作用 示例 适应性方针 适应性算法比较 进行网格适应 Fluent求解概述 求解过程 求解设置 选择求解器 离散方法 初始化 收敛性 监视收敛性 求解稳定性 亚松弛因子 Courant数设置 加速收敛 求解精度 网格独立性 网格自适应 选择求解器 求解器类型 压力基求解器(pressure-based solver) 密度基求解器(density-based solver) 共同点 通过网格划分将空间区域分解成由离散的控制体组成的集合。 在控制体上用积分形式构造离散变量的代数方程。 将离散方程线性化，然后通过求解线性方程获得变量的迭代解。 两种求解器的区别在于它们所使用的线性化方法和求解离散方程的方法不同 选择求解器 求解器类型 压力基求解器(pressure-based solver) 速度场通过动量方程求解 压力场通过对由连续方程和动量方程进行的压力修正求解 隐式求解 密度基求解器(density-based solver) 将连续性方程、动量方程、能量方程和组分输运方程作为一个方程组或者向量方程进行同时求解。其他的标量控制方程依次进行求解。 显式求解 + 隐式求解 选择求解器 若密度、能量、动量或组分相关度较高，最好选择密度基求解器 如：高速可压流体或通过化学反应动力学建模的流动 隐式密度基求解器比显式更好用 时间方面：隐式求解器求解速度大概快2倍 内存占用：隐式求解器占用显式求解器大概2倍的内存 显式密度基求解器一般只用在时间尺度特征与声波在同一个数量级的非稳态问题上 如：追踪瞬态激波 其他类型问题建议使用压力基（隐式）求解器 内存消耗远小于密度基求解器 压力基求解方法提供了更大的灵活性 方程离散方法 有限体积法中变量值存储与体心，需要进行差值得到面上的值： 对流项差值方法： 一阶迎风格式 一阶精度，收敛性最好 指数格式 Re5时精度比一阶迎风更好 二阶迎风格式 二阶精度，尤其适用于三角形或四面体网格，收敛较慢 QUICK格式 适用于四边形或六面体网格，三角形网格不适用。三阶精度 压力 - 速度耦合 压力 – 速度耦合是通过由连续性方程推导而来的通量守恒作为压力的一个条件来进行的。压力基求解器可对流动问题用分离或耦合方法求解。 分离式压力 – 速度耦合算法： SIMPLE 默认方法，鲁棒性较好 SIMPLEC 对于简单问题收敛更快，如无其它模型的层流流动 Skewness Correction PISO 属于SIMPLE系算法 对于非稳态流动问题或高斜度网格工况更有效 压力 - 速度耦合 耦合式压力 – 速度耦合算法 耦合方法在单相稳态工况中鲁棒性和效率更好，而且性能更优 耦合算法能更好地处理较差的网格或者大的时间步长。 压力基分离算法分别求解动量方程和压力修正方程，这种半隐式求解方法收敛较慢。 压力基求解器限制 压力基求解器对于下列问题不适用 Eulerian multiphase 模型 The non-iterative time advancement solver (NITA) 周期质量流量边界条件 初始化计算 迭代之前需要对所有求解变量进行初始化 Solve ? Initialize ? Initialize... 合理的初场有助于改善收敛性 某些工况，正确的初场很重要： 例如：用高温初场初始化燃烧区域 “Patch” 对于独立变量的某些区域 Solve ? Initialize ? Patch... 自由喷流(patch 出口高速) 燃烧计算(patch 高温?点火) 初始化计算 Patch窗口 余元（残差） 监视残差 Solve ? Monitors ? Residual... Plot ? Residuals... 余元（残差） Convergence criterion: absolute ?? 默认设置。对于稳态工况，只能选择absolute和none。每次迭代的残差和用户定义值比较，如果残差小于用户定义值，则认为该时间步长的计算收敛。 relative ?? 计算残差和该时间步长开始时残差比较，如果比值小于用户定义值，则认为该时间步长的计算收敛。 relative or absolute ?? Absolute或者relative中的某一个准则满足，则认为收敛 none ? 不检查收敛 收敛性 没有绝对的方法判断收敛 工况收敛特征： 所有的离散方程守恒误差在可接受范围内 流场不再随着迭代继续进行而变化 全场质量、动量、能量和其他标量平衡. 利用残差监视收敛性： 残差判断收敛性有时有效，有时也会误导 10-3 量