FLUENT-计算步骤.ppt

第二章FLUENT计算步骤 曹双华主讲07 03 本章通过计算一个二维弯管中的湍流流动和热传导过程再次演示FLUENT的求解过程 在这个算例中可以看到FLUENT计算的标准流程 其中包括 1 如何读入网格文件 2 如何使用混合的单位制定义几何模型和物质属性 3 如何设定边界条件和和物质属性 4 如何初始化计算并用残差曲线监视计算进程 5 如何用分离求解器计算流场 6 如何用FLUENT的图形显示功能检查流场 7 如何用二阶精度离散格式获得更高精度的流场 8 以温度梯度为基准调整网格以提高对温度场的计算精度 2 1问题概述 如图2 1所示温度为26 的冷流体流过弯管温度为40 的热流体从转弯处流入 并与主流中的冷流体混合 管道的尺寸如图2 1所示 单位为英寸 而边界条件和流体材料性质则采用国际单位制 入口处的雷诺数为2 03x105 因此必须使用湍流模型 图2 1弯管流动图示 2 2处理网格 网格处理包括网格的输入 检查 光顺 比例转换和显示等操作 2 2 1读入网格文件首先启动FLUENT的2D版 然后读入网格文件elbow msh File Read Case 2 2 2检查网格执行下列菜单操作 进行网格检查 Grid Check此时控制台窗口中会显示与网格有关的信息 包括网格空间范围 体积信息 表面积信息 节点信息等等 其中最需要检查的是网格单元的体积不能为负值 否则计算将无法继续下去 2 2 3光顺并转换网格执行下列菜单操作 打开Smooth Swap 网格光顺和转换 面板 Grid Smooth Swap 光顺网格可以提高网格质量 提高计算精度 在启动光顺 转换网格面板后 点击Smooth 光顺 按钮 然后再点击Swap 转换 按钮 控制台上将报告被转换网格的数量 反复点击Smooth 光顺 按钮和Swap 转换 按钮 直到报告中被转换 Swap 的网格数量降低为零 则网格光顺处理过程结束 图2 2Smooth SwapGrid 光顺 转换网格 面板 2 2 4按比例调整网格执行下列菜单操作 打开ScaleGrid 按比例调整网格 面板 Grid Scale 图2 3ScaleGrid 按比例转换网格 面板 在UnitsConversion 单位转换 下面的GridWasCreatedIn列表中选择in 英寸 作为网格创建数据单位 点击Scale 比例 按钮进行比例计算 注意在计算域范围 DomainExtents 中 网格的单位显示仍然是国际单位制中的m 米 点击ChangeLengthUnits 改变长度单位 按钮 可以将英寸变为长度单位 同时可以从面板下方计算域范围中看到x y坐标的最大值改变为64英寸 2 2 5显示网格执行下列菜单操作 打开GridDisplay 网格显示 面板 如图2 4所示 Display Grid 选择所有表面 然后点击Display 显示 按钮 FLUENT随即打开一个新的显示窗口显示网格 如图2 5所示 用右键点击图形显示窗口中的一个边界 则与这个边界相关的信息就会显示在控制台窗口中 在网格区域很多时 用这种方法可以比较方便地分辨各个区域 图2 4GridDisplay 网格显示 面板 图2 5网格显示 2 3计算模型 图2 6Solver 求解器 面板 2 3 1设置求解器参数执行下列菜单操作 打开Solver 求解器 面板 Define Models Solver 求解器的缺省设置如图2 6所示 本算例中可以保持这些缺省参数设置 2 3 2选择湍流模型执行下列菜单操作 打开Vicous 粘性 面板 如图2 7所示 Define Models Viscous 在Model列表中选择k epsilon 2eqn 即k 两方程模型 点击OK按钮接受系统对k 模型的缺省设置 包括标准模型设置 壁面函数设置和湍流模型常数的设置等等 图2 7粘性模型面板 2 3 3启用能量方程计算执行下列菜单操作 打开Energy 能量 面板 并激活能量计算 如图2 8所示 Define Models Energy 图2 8能量方程面板 2 4定义材料性质 创建一种新的物质 命名为water 水 并设置相关参数 首先执行下列菜单操作 打开Materials 材料 面板 如图2 9所示 Define Materials 在Name中键入材料名称water 然后设定水的各种物理性质 其中Density为1000kg m3 Cp为4216J kg K thermalconductivity为0 677W m K viscosity为8x10 4kg m s 最后点击Change Create按钮完成创建过程 在FLUENT提示是否用新的材料数据覆盖air的数据时 点击NO按钮予以否认 实际上FLUENT的材料性质数据库中已经存在水的数据 因此也可以简单地将相关数据直接从Database 数据库 中调入 如果个别数据与前面的数据有出入 可以直接在面板上进行修改 设置完毕后 关闭Materials 材料 面板 图2 9Materials 材料 面板 2 5定义边界条件 打开边界条件面板 如图2 10所示 Define BoundaryConditions 2 5 1设定流体条件在Zone 区域 列表中选择fluid 9 Type 类型 列表中相应地反白显示为fluid 流体 点击Set 设置 按钮打开Fluid 流体 面板 如图2 11所示 在MaterialName 材料名称 列表中选择water 水 点击OK按钮结束设置过程 图2 10边界条件面板 图2 11Fluid 流体 面板 2 5 2设置速度入口边界条件 1 设置velocity inlet 5上的边界条件 在Zone列表中选择velocity inlet 5 然后点击Set 按钮 进入速度入口设置面板 如图2 12所示 在VelocitySpecificationMethod 速度定义方法 中 选择Components 分量形式 定义速度 然后将X Velocity X向速度分量 设为0 2m s 将Temperature设为293K 在TurbulenceSpecificationMethod中 选择IntensityandHydraulicDiameter 强度和水力直径 定义湍流 然后将TurbulenceIntensity 湍流强度 设为5 将HydraulicDiameter 水力直径 设为32in 点击OK按钮结束对velocity inlet 5的设置 图2 12VelocityInlet 速度入口 面板 2 设置velocity inlet 6上的边界条件 方法同上 相关参数如下 Y Velocity Y向速度分量 为1 0m s temperature 温度 为313K TurbulenceIntensity 湍流强度 为5 HydraulicDiameter 水力直径 为8in 英寸 2 5 3设置压强出口边界条件pressure outlet 7上的压强出口边界条件的相关参数为 图2 13 表压0pascal 回流总温为300KBackflowDirectionSpecificationMethod为NormaltoBoundaryTurbulenceSpecificationMethod为IntensityandHydraulicDiameterBackflowTurbulenceIntensity为5 BackflowHydraulicDiameter为32in 图2 13PressureOutlet 压强出口 面板 2 5 4壁面条件设置壁面wall 4的边界条件保持缺省设置 其中HeatFlux 热流通量 的值为0 壁面wall 8也保持缺省设置 2 6求解过程 参数设置完毕并检查无误后 即可以开始流场计算 2 6 1初始化执行下列菜单操作打开SolutionIntialization 求解初始化 面板 如图2 14所示 Solve Initialize Initialize 在ComputeFrom 计算起始于 列表中选择velocity inlet 5 然后在YVelocity中加入0 2m sec 点击Init初始化 最后点击Close关闭 图2 14SolutionInitialization面板 2 6 2设置残差监视器执行下列菜单操作 打开ResidualMonitors 残差监视器 面板 Solve Monitors Residual 在Options 选项 下选择Plot 绘图 点击OK按钮关闭面板 图2 15ResidualMonitors面板 2 6 3保存算例文件执行下列菜单操作保存算例文件 File Write Case 2 6 4打开迭代面板执行下列菜单操作 打开Iterate 迭代 面板 Solve Iterate 将NumberofIterations 设为100 点击Iterate开始计算 在开始计算时残差监视器将同时被打开 2 6 5监视收敛过程图2 16是计算进行到第60步时残差曲线的走势 因为没有普适的收敛判断标准 所以在观察残差曲线时 不要仅仅监视残差曲线下降的数量级 最好同时能够监视相关流场变量的变化情况 图2 16残差监视器 简单地说 可以用三种方法判断计算是否已经收敛 1 观察残差曲线 可以在残差监视器面板中设置ConvergenceCriterion 收敛判据 比如设为10 3 则残差下降到小于10 3时 系统既认为计算已经收敛并同时终止计算 2 流场变量不再变化 有时候不论怎样计算 残差都不能降到收敛判据以下 此时可以用具有代表性的流场变量来判断计算是否已经收敛 如果流场变量在经过很多次迭代后不再发生变化 就可以认为计算已经收敛 3 总体质量 动量 能量达到平衡 在FluxReports 通量报告 面板中检查质量 动量 能量和其他变量的总体平衡情况 通过计算域的净通量应该小于0 1 FluxReports 通量报告 面板如图2 17所示 其启动方法为 Report Fluxes 图2 17FluxReports 通量报告 面板 2 6 6保存数据文件执行下列菜单操作保存数据文件 File Write Data 2 7显示计算结果 2 7 1等值线显示执行下列菜单操作 启动Contours 等值线 面板 如图2 18所示 Display Contours 首先显示速度值的填充等值线 在Contoursof下的列表中选择Velocity 和VelocityMagnitude 速度值 然后在Options下面选择Filled最后点击Display显示速度的等值线 用同样的方法可以显示温度 压强等变量的等值线图 图2 18Contours 等值线 面板 图2 19速度值的等值线图 2 7 2显示矢量图除了等值线图 还可以用矢量图观察流场中的速度矢量 首先执行下列菜单命令 启动Vectors 矢量 面板 Display Vectors 在Vectorsof下选择Velocity 然后点击Display就可以在不改变缺省设置的情况下显示速度矢量图 如果对显示效果不满意 可以在Vectors面板中修改显示设置 包括改变Scale 比例 和Style 风格 如果计算的是三维问题 则需要先在Surfaces下选择一个显示面 才能进行显示 显示的结果是速度矢量在这个面上的投影 图2 20Vectors 矢量 面板 图2 21速度矢量图 2 7 3XY曲线显示出口处的温度变化可以用XY曲线形式予以显示 首先启动SolutionXYPlot面板 Plot XYPlot 在YAxisFunction下面的列表中选择Temperature 和StaticTemperature作为Y轴变量 在Surfaces列表中选择pressure outlet 7 点击Plot 则静温在出口处的分布曲线就显示在图形窗口中 用同样方法还可以显示压强等变量在出口 入口等处的分布曲线 2 7 4定义场函数定义动压头的场函数 首先启动CustomFieldFunctionCalculator 场函数计算器 面板 Define CustomFieldFunctions 将新函数的名称加入NewFunctionName栏 然后根据动压头的定义公式进行公式编辑 动压头的输入内容如图2 24所示 在公式中需要输入流场变量 比如密度 速度 时 就在FieldFunctions 场函数 列表下进行选择 在使用数字和运算操作符的时候 则直接从面板左部的键盘中点取 输入完毕后 点击Define 定义 按钮完成函数定义 场函数的显示与其他场变量的显示方法是完全相同的 动压头的显示与速度值的显示方法相同 1 启动等值线面板 2 然后从Contoursof列表中选择CustomFieldFunction 修改场函数 和dynam head 动压头 3 设置显示选项 点击Display 显示 按钮