
FLUENT软件实际入门操作.doc
23页FLUENT软件实际入门操作 求解技术(Solve) Solve>Controls>Solution… 计算格式的选择 一阶迎风格式:适用于流淌方向与网格方向根本全都,构造化网格具有稳定性高,计算速度快的优点在网格方向与流淌方向不全都时,产生的数值误差比拟大 二阶格式:计算时间比拟长,收敛性差 适宜的计算方式:在计算开头时先用一阶格式进展计算以获得一个相对粗糙的解,在计算收敛后再用二阶格式完成计算以提高解的精度避开二阶格式收敛性差、计算时间长的问题,也避开了一阶格式在简单流场计算中数值误差大的问题 QUICK格式:对于构造网格计算旋转流淌问题时,计算精度高,但在其它状况下,QUCIK格式的精度与二阶格式相当 指数律格式:与一阶格式精度根本一样 中心差分:在LES湍流模型中使用,且应当在网格足够密集、局部Peclet数小于1的状况下使用 压强插值格式的选择 1在彻体力对流场有很大影响的状况下,应当选择彻体力加权 (body-force-weighted)格式。
2 在流场中有涡量很大的集中涡、高雷诺数自然对流、高速旋转流、多孔介质,以及流线曲率很大时,应当选择PRESTO!格式 3 对于可压流,应当使用二阶格式 4 二阶格式不能用于多孔介质计算和多相流计算中的混合物模型及VOF 模型在其他状况下,为了提高精度可以选用二阶格式 密度插值格式的选择 在用分别算法计算单相可压流时,有三种密度插值格式可供选择,即一阶迎风格式、二阶格式和 QUICK 格式一阶迎风格式具有良好的稳定性,但是在计算带激波的可压流时,会对激波解产生“抹平”作用,因此应当选用二阶格式或 QUICK 格式在用四边形网格、六面体网格或混合网格计算带激波的流淌时,最好使用 QUICK 格式计算全部变量需要留意的是,在计算可压多项流时,只能用一阶迎风格式计算可压缩相的流淌 Solve>Controls>Solution… Discretization(离散) 定义动量、能量、湍流淌能等工程,有一阶迎风格式、二阶迎风格式、指数律格式、QUICK格式和中心差分格式(在LES湍流模式计算中),也可以在使用耦合求解器时,定义湍流淌能、湍流耗散率等工程,并为这些工程选择一阶迎风格式、二阶迎风格式。
Pressure-Velocity Coupling(压强——速度关联) SIMPLE:FLUENT缺省设定格式 SIMPLEC:(1)稳定性较好,可以将亚松弛因子适当放大 (2)在层流计算中,在没有使用辐射模型等帮助方程时可以大大加 速计算收敛速度在简单流淌计算中,二者收敛速度相差不多 PISO:(1)非定常计算 (2)允许使用较大的时间步进步行计算,可以缩短计算时间(大涡模拟 LES这类网格划分较密集,时间步长小的除外) (3)处理网格畸变较大的问题,假如在 PISO 格式中使用邻近修正(对 动量进展修正)(neighbor correction),可以将亚松弛因子设为1.0或接近于1.0的值而在使用畸变修正(skewness correction)时,则应当将动量和压强的亚松弛因子之和设为1.0,比方将压强的亚松弛因子设为 0.3,将动量的亚松弛因子设为0.7假如同时采纳两种修正 形式,则应将全部松弛因子设为1.0或接近于1.0的值 Under-Relaxation Factors(亚松弛因子) 条件:简单流淌状况下,缺省设置不能满意稳定性要求,计算过程中可能消失振 荡、发散等状况,减小亚松弛因子,以保证计算收敛。
假如实际的残差曲线向上进展,则中断计算,适当调整亚松弛因子 计算发散时,可以考虑将压强、动量、湍流淌能和湍流耗散率的亚松弛因子的缺省值分别降低为 0.2、0.5、0.5、0.5 库朗数 库朗数越大,时间步长就越长,计算收敛速度就越快库朗数都应当在允许的范围内尽量取最大值 Solve>Controls>Limits…设置求解极限 意义:避开消失非物理解,比方密度或温度变成负值,或者大得远选超过真实值 Solve>Initialize>Initialize…全局初始化 开头进展计算之前,必需为流场设定一个初始值设定初始值的过程被称为“初始化” 初始化的步骤: (1) 设定初始值 ? 用某个区域上设定的初始值进展全局初始化 Compute From----选择需要定义初始值的区域名----Initial Values ? 用平均值的方法对流场进展初始化 Compute From----all-zones----Initial Values (2)对于动网格,通过“Absolute”or“Relative to Cell Zone”来打算。
(3)点击Init Solve>Initialize>Initialize…初始值进展局部修补 步骤: (1) 在Variable(变量)列表中选择需要修补的变量名 (2) 在Zones To Patch(需要修补的区域)或Registers To Patch(需要修补的 标记区)中选择需要修补变量所在的区域 (3) 假如需要将变量的值修补为常数,则直接在输入栏中输入变量的值假如 需要用一个预先设定的函数定义变量,可以在Use Field Function(使用场变量函数)中的Field Function(场函数)列表中选择适宜的场函数 (4) 假如需要修补的变量为速度,则除了定义速度的大小,还要定义速度为绝 对速度,还是相对速度 (5)点击Patch(修补)按钮更新流场数据 流场求解 Solve>Iterate…定常问题计算 时间相关问题的计算 FLUENT模拟的时间相关问题:旋涡卷动、可压缩振荡、瞬态热传导、瞬态化学反响与混合流淌等问题,比方圆柱后面卡门涡街中旋涡的生成和脱落过程等等。
求解器设置步骤: (1) 求解器面板的启动步骤为:Define->Models->Solver... 首先选中Unsteady(非定常)选项,然后再在Unsteady Formulation(非定常算法)中选择所需要的算法在选择耦合显式求解器时,可以使用 Explicit(显式)格式显式格式主要用于激波等波动解的捕获问题 在使用分别求解器时,还可以选择 Use Frozen Flux Formulation(使用冻结通量格式)选项,但是这个选项只能在用分别求解器计算单相流问题时使用,同时计算中还不能使用动网格技术 (2) 定义相关物理模型和边界条件 (3) 假如选用了分别求解器,则应当在压强-速度关联方法中选择 PISO格式 在用LES湍流模型计算湍流问题时,不应当选择 PISO格式,那样会增加系统资源消耗 (4) 假如使用的是显式非定常计算或适应性时间推动法,应当在每步迭代中打 开当前时间或当前时间步长的print选项这项操作在 Statistic Monitors(统计数据监视器)面板中完成,通过以下菜单操作可以进入这个面板: Solve->Monitors->Statistic... 在 Statistics(统计数据)列表中选择time(时间)或delta_time(时间增量),并翻开Print选项。
当FLUENT在掌握台窗口打印残差信息时, 就会消失当前时间或当前时间步长 (5) 可以用Force Monitors(受力监视器)或 Surface Monitors(外表监视器) 面板监视随时间变化的受力系数的值或平均值、质量平均值、积分或场变量通量,或者任意一个函数在外表上随时间变化的状况 (6) 用Solution Initialization(求解初始化)面板设定初始条件: Solve->Initialize->Initialize... 还可以读入一个定常状态的数据文件来设定初始条件: File->Read->Data... (7) 用自动保存功能定义文件名和在迭代过程中自动保存算例文件和数据文 件的频率: File->Write->Autosave... 还可以用 Execute Commands(执行命令)面板自动执行系统命令 (8) 假如想创立解的动画过程,可以使用 Solution Animation(求解动画)面 板设置动画显示参数 (9) 监视时间数据的步骤如下: 1)在 Iterate(迭代)面板中翻开 Data Sampling for Time Statistics(时间数据采样)选项: Solve->Iterate... 翻开这个选项可以同时显式平均值和均方根值。
2)初始化流场统计数据: Solve->Initialize->Reset Statistics 用这个菜单还可以在收集了一些数据后重置流场数据比方在将时间采样选项翻开的状况下先计算10步,检查计算结果,然后再计算10步,则时间统计数据中将包括前面10步的数据,除非重新对流场统计数据进展初始化 (10)指定时间相关解的参数并开头计算,对于定常、非定常计算的步骤如下: 1)假如选择一阶或二阶精度隐式计算,步骤如下: Solve->Iterate... 通常在计算开头时,时间步长取得较小,而在后续的计算中则可以逐步增加时间步长在缺省设置中,时间步长是固定的,假如需要让时间步长随计算进程发生变化,可以选择Adaptive Time Stepping(适应性时间推动)法中的 Adaptive(适应性步长)选项 2)显式非定常计算的步骤如下: 在Solution Controls(求解过程掌握)面板中保存Solver Parameter(求解器参数)缺省设置: Solve->Controls->Solution... Solve->Iterate...进展迭代 监视残差Solve->Monitors->Residual... or Plot->Residuals... 分别算法中的残差是守恒方程两端数值的差值在全部网格点上的和,耦合算法中的残差则是各变量的时间变化率。












