AdiabaticFluent辐射与自然对流模拟引言在这个算例中,将会解决二维方箱中的辐射与自然对流相结合的问题,网格采用四边形单元网格在这个算例中将会学到以下知识点:1. 应用Fluent中各种辐射模型Rosseland;2. 使用Boussinesqmodel定义密度;3. 设定辐射与自然对流传热问题的边界条件;4. 将单一的墙划分为多个墙区域;5. 对已有的流体物性进行修改;6. 用隔离求解器求解;7. 显示速度矢量和流函数等值线,以及温度等值线问题描述将被考虑的问题如图5.1所示,一个边长为L的正方形箱体,右墙温度为2000K,左墙温度为1000K,上下墙绝热,重力向下,由于热重引起密度梯度所以发展为浮力流箱体中的介质被认为是有吸收性和散射性的,因此墙壁间的辐射交换因存在吸收被减弱,同时也因为介质的散射作用而增强了所有墙壁被认为是黑体,目的在于应用有效的辐射模型计算箱体中流场和温度场分布,以及墙壁的热流量,并且对于不同光学深度aL比较所表现出的特性工质普朗特数大约为0.71,基于L的雷诺数为500000,这说明流动相当于层状流动,应用Boussinesq假设来模拟浮力流动普朗克数为0.02,用于考虑传导与辐射的相对重要性,其中,T0=(Th+Tc)/2。
在这个算例中将有三种opticalthickness的情况会被考虑到,分别是aL=0,aL=0.2,andaL=5注意:物理属性和工作条件(重力加速度)都已经给定以适合于产生的想要的普朗特数,雷诺数和普朗克数如下图所示:p=1000kg/c=11030x1Pk=15309VM=10"3kgp=10'51/g=-6,96x11』=0,0工5L=1mJKa=5x10Fr=0.71Fl=0.02第1步:网格将网格文件rad/rad.msh拷至fluent的工作目录下(就像在指南1中描述的一样),并起动fluent的二维单精度解算器1. 读取网格文件rad.msh.FileReadCase...当网格读入的时候,在Fluent控制窗口会显示相应的信息,会报告网格有2500个单元2. 检查网格质量GridCheck...Fluent会对网格进行各种各样的检查,并会在控制窗口显示信息特别注意最小体积,确保它是正数值3. 显示出网格(如图5.2)注意:此时所有的墙体为一个整体,wall-4下面需要将其分成四个独立的区域,才能对不同的区域定义不同的边界条件4. 将单一的墙分为四部分GridSeparateFaces...在Options.在Zones列定义角度值单击Sepa)b)c)d)arate。
带法向矢量的面在被89度分隔后会成为单独的区域,既然四个墙的区域是正交的,wall-4就被分为四个独立的区域5. 重新显示出网格a)在griddisplay对话框中选定所有的面后单击Display注意:现在有四个不同的墙区域而不是一个另外:你可以用鼠标右键来验证墙区域名称所对应的墙在图形窗口中用鼠标右键点击任何一个边界,区域名称,类型都在Fluent控制窗口中显示出来当同一个类型有几个不同区域,而你又需要迅速的区分开来时,这个用途的作用就特别明显有些时候,可以使内部网格不显示,以便更准确的选择边界第2步:模型选择辐射模型采用的是Rosselandmodel1. 选择默认解算器HSolverFormulationSolverVelocityFormulation<•Absolute厂RelativeGradientOption<•Green-GaussCellBased厂Green-GaussNodeBased「LeastSquaresCellBasedPorousFormulation<•SuperficialVelocityPhysicalVelocity2.打开Rosseland辐射模型。
DefineModelsRadiation...HRadiationModelModelOffRosselandP1DiscreteTransfer(DTRM]SurfacetoSurface(S2S)DiscreteOrdinates(DO)IHelp会中现一个信息提示框,告诉你新Cancel1”口十—当你在RadiationModel的材料物性被添加了,你将在后面设置物性参数,因此现在只需单击ok确认这个信息即可注意:当你激活辐射模型时,fluent会自动打开能量解算器,因此不用另外再激活Energy对话框了3.在模型中添加重力影响DefineOperatingConditions...|£J|OperatingCoinditions1)PressureOperatingPressure(pascal)|101325RefereneePressureLocationGravityOK|Zncel|Help—|开Gravity2) 设置y方向的GravitationalAcceleration值为-6.94e-5m/s2文中提到过,重力加速度调整到产生合适的无因次数(普朗特数,雷诺数和普朗克数),见图5.1及相关的说明。
3) 设置工作温度为1000k.在下一步要激活的theBoussinesqmodel中将会用到工作温度值第3步:材料物性默认的流体材料为空气,在此次要解决的问题中工作介质也正是空气但却是假想的流动,使各个物性都适合于给出特定的无因次数,因此必须对默认的空气参数进行修改应用光学深度aL=0.2进行计算假设L=l,吸收系数则应设为0.2)在后面,aL=5或者无光学介质(aL=0)时,计算得到的结果将会用来比较不同的光学介质对于的辐射模型的不同DefineMaterials...NrihaMidleinlullypeOrderMatsrialsftyairIfluid茴NameCheiiiiciJlbDrmuliiFluentFluidMaterialsrChemicalrormtils1|dlrJrinsntDatabase...口MaterialsMlxturp:nonelooperliesUKpr-DRflnfid3Density(kg/rnl)boussincsqtdiCp『kgk)COnislclBl
6. 在ScatteringCoefficient和ScatteringPhaseFunction中保持默认值,在要解决的问题中不涉及到散射问题;7. 设定热扩散系数(用boussinesq模型时)为1e-5K-18. 单击Change/Create,关闭Materials面板第4步:边界条件DefineBoundaryConditions・・・1.对底部的墙设置边界条件(wall-4.006)注意:如果不能确定wall-4.006对应的是否为底墙,可以用鼠标右键在图形窗口中单击底墙,在边界条件面板中Zone列表中会自动显示相应区域的名称在设定其它墙的边界条件时也可以这样做,一定要保证设定正确的边界条件1) .将ZoneName改名为left.(2) 在ThermalConditions下选择Temperature,并设置温度值为1000k;3•对右墙设置边界条件(wall-4:007.)(1) 将ZoneName改名为right;(2) 在ThermalConditions下选择Temperature,并设置温度值为2000k;4.对顶墙设置边界条件(wall-4:005);(1) 将ZoneName改名为top;(2) 对于绝热墙,保持默认设置的热边界条件(heatflux=0)。
第5步:Rosseland模型求解1.为控制解设置参数SolveControlsSolution...(1) 在Equations和Under-RelaxationFactors.下,保持默认值2) 在Discretization下,Pressure选择PRESTO!,MomentumandEnergy选SecondOrderUpwind2. 流场初始化SolveInitializeInitialize...HSolutionInitialization(1)将温度设置为1500k,3.计算时,激活残Residual...SolveMonitorsGaugePressure(pascal)|qXVelocity(m⑸|qYVelocity(m⑸[g单击1(k)[1500差曲线显示IJ®RelativetoCellZone厂AbsoluteComputeFromRefereneeFrameInitialValues在Options单击okRosseland(1)(2)注意:线来4. 保存case文件FileWriteCase5. 开始计算,进行SolveIterate..大约在180次时计算就会收敛。
6. 保存数据文件(rad_ross.dat).FileWriteData...第6步:Rosseland模型后处理1.显示收敛曲线SolveMonitorsResidual..plotResidualscontinuityx-velocityy-veloGityeneray020钩6QSO10012014016QIterations2.显示速度矢量FLlDMayF^e^iors...2.11e-042„01e-041.90e-04180e-04169e-04158e-041.48e-041.37e-04<4HJl8.45a-a57.40e-05034e-O55,28e-054.23e-05I3I7e-O52,11e-05。