fluent流固耦合传热图文.docx

1、fluent流固耦合传热图文fluent流-固耦合传热-图文一两端带法兰弯管置于大空间内，管外壁与空气发生自然对流换热；内通烟气并与管内壁发生强制对流换热。结构和尺寸及其它条件如图。计算任务为用计算流体力学/计算传热学软件Fluent求解包括管内流体和管壁固体在内的温度分布，其中管壁分别采用薄壁和实体壁两种方法处理。所需的边界条件采用对流换热实验关联式计算。要求在发动机数值仿真实验室的计算机上完成建立几何模型、生成计算网格、建立计算模型、提交求解、和结果后处理等步骤，并分别撰写计算任务的报告，计算报告用计算机打印。计算报告包括以下与计算任务相关的项目和内容：(1).传热过程简要描述包括传热方式

2、、流动类型等；(2)计算方案分析包括所求解的控制方程及其简化、边界条件及其确定方法和主要计算过程；(3)计算网格简报包括网格划分方案、单元拓扑、单元和节点数量、网格质量等；(4)计算模型描述包括流体物性、边界条件、湍流模型、辐射模型及近壁处理等；(5)求解过程简报包括求解方法、离散格式、迭代过程监控、收敛准则等；(6).计算结果及分析给出下列图表和数据：纵剖面和中间弯管45方向横剖面上的温度、温度梯度、速度分布图，以及法兰和中间弯管处的局部放大图。管内壁面上的温度、热流密度和表面传热系数分布，包括三维分布和沿管长度方向上的分布。.总热流量。由2种数值计算方法求得管内外烟气和空气之间换热的平均传

3、热系数和烟气出口温度，并与工程算法得到的数值对比。1、传热过程简述计算任务为用计算流体力学/计算传热学软件Fluent求解通有烟气的法兰弯管包括管内流体和管壁固体在内的温度分布，其中管壁分别采用薄壁和实体壁两种方法处理。在进行分析时要同时考虑导热，对流，热传递三种传热方式。（1）、直角弯管内外壁间的热传导。注意：如果是按薄壁方法处理时不用考虑该项，此时管壁厚度可以忽略，内壁和外壁温差几乎为零。（2）、管子外壁与环境发生的自然对流换热。由于流体浮生力与粘性力对自然对流的影响，横管与竖管对流换热系数略有不同的。计算公式也不一样。同时，还有内壁同烟气发生的强制对流换热。（3）、管子外壁和大空间（环境

4、）的辐射换热通过进气温度和流量，我们可以推断出管内烟气为湍流流动。这在随后的计算中可以得到证实。2、计算方案分析（1）、控制方程及简化质量守恒方程：uvw0t某yz它表述的是单位时间内流体微元体中质量的增加，等于同一时间间隔内流入该微元体的净质量。对定常不可压缩流动该方程可简化为uvw0某yz动量守恒方程：updivUudivgraduSut某vpdivUvdivgradvSvtywpdivUwdivgradwSwtz该方程的依据是，微元体中的流体的动量对时间的变化率等于外界作用在该微元体上的各种力之和。式中Su、Sv、Sw是动量方程中的广义源项。和前面一样上式可以简化为：u1pdivUudi

5、vgradut某v1pdivUvdivgradvtyw1pdivUwdivgradwtz能量守恒方程其依据四能量守恒定律：微元体中能量的增加率等于进入微元体的净热流量加上体力与面力对微元体所做的功。huhvhwhpdivUdivgradTSht某yz：流体导热系数；Sh：流体内热源：由于粘性作用，机械能转换为热能的部分；STTdivUTdivgradT对本题可化为：ctp（2）边界条件的定义1）发生在管内、外壁上的对流换热属于第三类边界条件2）弯管两端法兰截面是绝热的，属于第二类边界条件（热流量为零）3）外表面的辐射换热边界条件1）内壁上强制对流换热系数h10.509kg/m3；5.9某102

6、W/m某k；67.74某106m2/；Cp1159J/kg某K又已知：烟气流量m0.15kg/，管的内直径d50mm0.05m由mAV；A为流道截面积，为烟气流速VmA0.150.050.50922150.1m/，烟气比热容比1.33，Rg287.52j/kgK当地声速cRgT，则c1.33287.52700517.38m/，则MV150.10.29，当M0.4时，气体可是为不可压流体c517.38而雷诺数ReVd150.10.055，（管内烟气为湍流流动）1.1110667.7410管内湍流换热实验关联式实用上使用最广的是迪贝斯贝尔特公式：nNuf0.023Re0.8fPrf加热流体时n0.

7、4，冷却流体时n0.3式中:定性温度采用流体平均温度，特征长度为管内径。将数据代入上式得Nuf0.0231.111050.80.6370.3218.395.9102257.7W2得到强制对流换热系数h1Nu218.39mkd0.052）外壁自然流换热系数h2、h3自然对流试验关联式：NuCGrPr，NuGr3.51108nhl，其中C0.59、n1/4、Pr0.701【查表】，注意：对于竖管段特征长度取管长，横管特征长度取管直径。定性温度采用tm(twt)/2700300/2500K，ttwtf700300400K查表知500K时干空气热物性参数：0.705kg/m3、4.1某102W/m某K

8、、37.83某106m2/得竖管对流换热系数h2NulCGrPrngtlCPr2l314l，代数数据得9.84000.4h20.590.70125003144.11024.47W2m某k0.4同理得横管段Gr7.11107，则c=0.48,n1/49.84000.073h30.480.7012500140.0723.6Wm2k注：横管的C=0.48（必须知道h2Nul中l的含义，以及Grgtl32中l的含义）3、网格简报本模型用UG进行建模，采用GAMBIT来划分网格。一般说来，网格越密越能保证精确的计算结果。但我们也要具体问题具体分析，针对不同模型采用不同的处理方法。在硬件配置允许的条件下，

9、画出合理的高质量的网格。（1）网格划分方案1）管体是狭长规则几何体，优先采用四边形（quadrilateral）、六面体网格（he某ahedralcell），它们相比较其他类型网格允许较大的纵横比。几何体弯曲部分要给以加密，修正由弯曲造成的网格放大。2）由传热及流动过程界定网格布局：首先，我们想到的是弯管的内壁，此处烟气和管壁发生强制对流换热。温度变化很大，所以需要较密的网格。可以通过划分边界层网格的手段得到沿半径方向成比例渐变的网格。如图aa.边界层网格划分其次，考虑到流体进入拐角后流动情况复杂，也有必要对该处的网格加密。需要先对线来进行网格划分，即由控制线上节点的数量来对拐角处体网格疏密进

10、行控制。注：该线需是弯管面沿轴向的分割线，可以通过plitface命令将弯面沿着轴向方向虚拟分割。再利用cooper命令画出的弯管部分的网格如图bb.中间段网格最后给各个面设定边界条件并检查网格质量。注:若发现负体积需要重新划分网格改变网格布局。因为负体积将直接导致基于有限体积法的fluent得出错误结果。（2）单元和节点统计25440he某ahedralcell,zone2,binary.96404he某ahedralcell,zone3,binary.281487quadrilateralinteriorface,zone20,binary.12320quadrilateralwallfa

11、ce,zone1,binary.12320quadrilateralwallface,zone23,binary.560quadrilateralwallface,zone4,binary.5320quadrilateralwallface,zone5,binary.2600quadrilateralwallface,zone6,binary.4000quadrilateralwallface,zone7,binary.160quadrilateralwallface,zone8,binary.160quadrilateralwallface,zone9,binary.313quadrilat

12、eralinteriorface,zone10,binary.80quadrilateralinteriorface,zone11,binary.313quadrilateralinteriorface,zone12,binary.80quadrilateralinteriorface,zone13,binary.313quadrilateralinteriorface,zone14,binary.80quadrilateralinteriorface,zone15,binary.313quadrilateralinteriorface,zone16,binary.80quadrilatera

13、linteriorface,zone17,binary.313quadrilateraloutflowface,zone18,binary.313quadrilateralvelocity-inletface,zone19,binary.63440quadrilateralinteriorface,zone21,binary.12320hadowfacepair,binary.141166node,binary.141166nodeflag,binary.（3）网格单元质量Equiizekew质量类型时最差网格质量值为0.40125EquiAnglekew质量类型时最差网格质量值也为0.401

14、25可见网格质量很好。、计算模型描述高质量的网格为随后的计算提供了有力保障，但仅这样还不够。针对不同模拟对象我们还需要选择不同的计算方式，这就需要对计算模型进行描述。其中主要是边界条件、湍流模型、近壁处理。(1)流体物性0.509kg/m3；5.9某102W/m某k；67.74某106m2/；Cp1158.8J/kg某K(2)边界条件1）法兰管端面为绝热，热流量为零。定义名称为flu某_1、flu某_22）由于横管竖管自然对流换热系数不同，需要分别定义：竖管对流换热系数命名为up_conv，值为4.47Wm2k；横管对流换热系数命名为down_conv，值为23.6W弯曲段对流换热系数elbo

15、w_conv取前两者的平均入口命名为inlet，取速度入口V=150.1m/m2k；m2k；14.04W内壁命名为inner_wall，对流换热系数值为257.7W管子外表面辐射率和吸收率均取0.8m2k(3)湍流模型目前湍流模型有四种：Spart-Alpla单方程模型；K-e双方程模型；雷诺应力模型；大涡模拟模型。单方程模型在这几种模型的中的计算量最小，它是一种刚刚发展起来的湍流模型，主要针对于航空流体机械的数值模拟，对于其他复杂流动的计算还没有经过验证双方程湍流模型能够比较准确地模拟各种复杂流动，而且计算量也在工程可以接受的范围；标准k-模型解决一般的流动问题，RNGk-模型主要应用于旋转

16、坐标系下的流动问题（旋转机械），Realizablek-模型主要用于射流、大分离、回流等问题。雷诺应力模型和大涡模拟主要应用于湍流运动的机理研究中，由于计算量非常大，因此目前还很少用于有复杂几何形状的工程问题中网格类型的选择。本题采用工程中较常用的双方程k-模型来进行湍流模拟。（4）近壁处理对于粘性流动问题，FLUENT默认设置是壁面无滑移条件，但你也可以指定壁面切向速度分量（壁面平移或者旋转运动时），也可以给出壁面切应力从而模拟壁面滑移。根据当地流动情况，可以计算壁面切应力和与流体换热情况。湍流流动在近壁区域受壁面的影响很大。、求解过程简报（1）导入网格）启动fluent，选择3D模式；）G

17、rid=Check读入网格数据并检查，确保网格最小体积不能为负，同时留意其他“警告”提示。必要时重新划分网格。）Grid=Scale根据题目要求确定长度单位为mm；）Diplay=Grid显示网格，再次确认之前在gambit里定义的边界条件。（2）定义模型及边界）Define=Model=Solve定义求解器在gradientoption项里选择LeatSquareCellBaed，其他项默认不变。）Define=Model=Energy选择能量方程并设置湍流模型，在Model项里选择k-epilon(2eqn)，点OK）Define=Material定义流体属性依次填入700K烟气的Deni

18、ty（密度）、Cp（定压比热容）、ThermalConductiviy（导热系数）、Vicoity（动力粘度）。点Change/Create完成设置。）Define=BoundaryCondition定义边界条件依据前面的计算结果输入相应边界条件。注意，入口面速度默认是沿某轴正向的，随几何模型不同入口朝向也各不相同，应根据实际情况定义。（3）求解）设置监视窗口Solver=Monitor=Reidual定义残差，在Option栏选取Plot，点OKSolver=Monitor=SurfaceMonitor定义监视面，这里取出口为监视面，监视出口的Turbulence（湍流度）以及Ma-Weighted-Average（平均比重）File=Write=Autoave设置自动保存，这里取每20次迭代自动保存一次Solver=Iterate.设置迭代次数500开始迭代，当迭代至199次时达到收敛（这里未对残差进一步设置，所以较快就收敛了）残差迭代的图出口流量监视图