基于AVL-FIRE的CFD分析过程及案例操作.pptVIP

分析流程 线网格：表面网格的特征边缘线 * 分析流程 划分面网格——体网格 * 设定整体网格最大最小尺寸 对需要细化的selection单独设定尺寸 设定边界层数 其他网格处理 生成体网格 网格质量检查 减小尺寸 分析流程 网格质量 * 分析流程 网格质量 * 分析流程 网格质量 * 分析流程 网格质量 * 分析流程 网格质量 * 边界条件 边界条件 所谓边界条件，是指在求解域的边界上所求解的变量或其一阶导数随地点及时间变化的规律，只有给定了合理的边界条件，才能计算得出流场的解。 FIRE中边界条件分类 * 边界条件 进/出口边界 对称边界 壁面边界 边界条件 进/出口边界 Static pressure：静压 Total pressure：总压 Stagnation pressure：用于亚音速和超音速流动问题，需输入动压和马赫数 总压=静压+动压 Normal velocity：法向速度 Velocity：速度（x/y/z任一方向速度） * 边界条件 进/出口边界 Mass flow：质量流量 Gradient=0： Gradient =0-pressure：作为静态压力边界 Gradient=0-all variables：不保证整体的连续，适用于可压缩流动 * 边界条件 1）对于风道稳态计算，一般进口用流量，出口用静压或梯度为零；若是瞬态计算，往往进出口都采用压力边界条件，如果计算不易收敛，则建议进出口都采用流量（出口流量的值要改变符号） 2）对于冷却水套一般进口定义流量，出口定义静压或梯度为零。 建议是有实测的静压值，这样有助于加快计算收敛 3）入口处湍流值的给定： 一般turb.ref.velocity处填入进口平均速度（可根据流量估算） % of mean velocity一般填1-10 % of hydraulic diameter一般填5-10 这样Turb.kin.energy及下面两项都会由程序自动算出 一般来说，TKE的值大一些对计算收敛有帮助 * 初始条件 流体物性可预设为常量(空气,流体或新物性), 也可用公式定义；所有变量在全场的初始化可通过均匀初始化和势流场初始化来实现；对于瞬态计算，初始条件必须准确给定，因为后续时间步的结果直接受其影响。 * 离散 计算边界值：Extrapolate and Mirror Mirror 对于边界网格质量不好的情况更为适用，可作为默认选项 计算导数的方法：Least Sq. Fit（最小二乘法） and Gauss（高斯法），Gauss 作为默认选项 Cell face adjustment: 限制扩散项，影响能量守衡；但对质量不好的网格计算有帮助 人工可压缩性（Artificial Compressibility)：实现了每个单元采用不同的时间步长 ，这样松弛因子可以加大 ，减少计算时间 * 算法 Simple Simple/c Simple/ Piso Simples算法是基础算法，piso算法能减少30%-50%的计算时间 SIMPLEC 对于一些压力-速度耦合起的作用较大的应用会得到更好的结果. 而在一些有其他源项的情况下, 如较强的湍流, 喷雾, 燃烧时SIMPLE 法算得会更好 * 控制方程 S-A单方程模型：一般用于高马赫数流动 涡粘性/耗散模型（k-e）：双方程模型，复杂流动的模拟不够准确。优点是计算稳定性好，对计算资源的要求和花费低 k-??f：四方程模型，精度和稳定性都较好，推荐使用，计算时间仅比k-e模型多15% RSM模型：对每个湍流应力分量的求解能准确模拟湍流应力场及其各相异性。缺点是对计算资源有很高的要求，计算稳定性较差 AVL复合湍流模型 （HTM)：结合k-e模型与RSM模型的优点 * 控制方程 * 壁面处理 壁面的层流区对于各种湍流方程不适用，与之对应的解决方法有两种 壁面函数法 低RE数法 * 收敛标准 * 最大迭代次数 残差满足收敛 收敛？ 物理判断 收敛 案例操作 案例描述：常温流体（空气）以指定的质量流量从进口面，进入计算域，从出口面流出，计算该计算域内的流动 * 入口 出口 边界条件 进口 指定质量流量：m=0.0381kg/s 指定温度：T=293K 指定湍流条件：湍动能+湍流长度尺寸 出口：静压：P=1bar * 导入STL * 创建Selections * 生成特征线 * 划分体网格 用Hybrid Assistant 划分混合六面体网格 * 求解器设置 * 边界条件 初始条件 流体属性 控制方程 收敛标准 计算监控 * 结果分析+后处理 * 速度切片图 结果分析+后处理 * 总压切片图 * 基于AVL-FIRE的CFD分析过程及案例操作