fluent流--固耦合传热-图文

fluent流--固耦合传热-图文
一两端带法兰弯管置于大空间内，管外壁与空气发生自然对流换热；内通烟气并与管内壁发生强制对流换热。

结构和尺寸及其它条件如图。

计算任务为用计算流体力学/计算传热学软件Fluent求解包括管内流体和管壁固体在内的温度分布，其中管壁分别采用薄壁和实体壁两种方法处理。

所需的边界条件采用对流换热实验关联式计算。

要求在发动机数值仿真实验室的计算机上完成建立几何模型、生成计算网格、建立计算模型、提交求解、和结果后处理等步骤，并分别撰写计算任务的报告，计算报告用计算机打印。

计算报告包括以下与计算任务相关的项目和内容：
(1)...............................传热过程简要描述包括传热方式、流动类型等；(2)计算方案分析包括所求解的控制方程及其简化、边界条件及其确定方法和主要计算过程；
(3)计算网格简报包括网格划分方案、单元拓扑、单元和节点数量、网格质量等；
(4)计算模型描述包括流体物性、边界条件、湍流模型、辐射模型及近壁处理等；
(5)求解过程简报包括求解方法、离散格式、迭代过程监控、收敛准则等；(6)................................................计算结果及分析给出下列图表和数据：纵剖面和中间弯管45°方向横剖面上的温度、温度梯度、速度分布图，以
及法兰和中间弯管处的局部放大图。

管内壁面上的温度、热流密度和表面传热系数分布，包括三维分布和沿管长度方向上的分布。

............................................................ .........................................总热流量。

由2种数值计算方法求得管内外烟气和空气之间换热的平均传热系数和烟气出口温度，并与工程算法得到的数值对比。

1、传热过程简述
计算任务为用计算流体力学/计算传热学软件Fluent求解通有烟气的法兰弯管包括管内流体和管壁固体在内的温度分布，其中管壁分别采用薄壁和实体壁两种方法处理。

在进行分析时要同时考虑导热，对流，热传递三种传热方式。

（1）、直角弯管内外壁间的热传导。

注意：如果是按薄壁方法处理时不用考虑该项，此时管壁厚度可以忽略，内壁和外壁温差几乎为零。

（2）、管子外壁与环境发生的自然对流换热。

由于流体浮生力与粘性力对自然对流的影响，横管与竖管对流换热系数略有不同的。

计算公式也不一样。

同时，还有内壁同烟气发生的强制对流换热。

（3）、管子外壁和大空间（环境）的辐射换热
通过进气温度和流量，我们可以推断出管内烟气为湍流流动。

这在随后的计算中可以得到证实。

2、计算方案分析
（1）、控制方程及简化
质量守恒方程：
uvw0t某yz它表述的是单位时间内流体微元体中质量的增加，等于
同一时间间隔内流入该微元体的净质量。

对定常不可压缩流动该方程可简化为
uvw0某yz动量守恒方程：
updivUudivgraduSut某vpdivUvdivgradvSvtywpdivUwdivgradwSwtz
该方程的依据是，微元体中的流体的动量对时间的变化率等于外界作用在
该微元体上的各种力之和。

式中Su、Sv、Sw是动量方程中的广义源项。

和前面一样上式可以简化为：
u1pdivUudivgradut某
v1pdivUvdivgradvtyw1pdivUwdivgradwtz
能量守恒方程
其依据四能量守恒定律：微元体中能量的增加率等于进入微元体的净
热流量加上体力与面力对微元体所做的功。

huhvhwhpdivUdivgradTSht某yz：流体导热系数；Sh：流体内热源
：由于粘性作用，机械能转换为热能的部分；
STTdivUTdivgradT对本题可化为：
ctp（2）边界条件的定义
1）发生在管内、外壁上的对流换热属于第三类边界条件
2）弯管两端法兰截面是绝热的，属于第二类边界条件（热流量为零）3）外表面的辐射换热边界条件
1）内壁上强制对流换热系数
h1
0.509kg/m3；5.9某102W/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.4，冷却流体时n0.3
式中:定性温度采用流体平均温度，特征长度为管内径。

将数据代入上式得Nuf0.0231.111050.80.6370.3218.39
5.9102257.7W2得到强制对流换热系数h1Nu218.39mkd0.052）外壁自然流换热系数h2、h3
自然对流试验关联式：NuCGrPr，NuGr3.51108
nhl，
其中C0.59、n1/4、Pr0.701【查表】
，注意：对于竖管段特征长度取管长，横管特征长度取管直径。

定性温度采用tm(twt)/2700300/2500K，
ttwtf700300400K
查表知500K时干空气热物性参数：
0.705kg/m3、4.1某102W/m某K、37.83某106m2/得竖管对流换热系数h2NulCGrPrngtlCPr2l314l，代数数据得
9.84000.4h20.590.70125003144.11024.47W2
m某k0.4同理得横管段Gr7.11107，则c=0.48,n1/4
9.84000.073h30.480.7012500140.0723.6Wm2k注：横管的C=0.48（必须
知道h2Nu
l中l的含义，以及Grgtl32中l的含义）
3、网格简报
本模型用UG进行建模，采用GAMBIT来划分网格。

一般说来，网格越密越能保证精确
的计算结果。

但我们也要具体问题具体分析，针对不同模型采用不同的处理方法。

在硬件配置允许的条件下，画出合理的高质量的网格。

（1）网格划分方案
1）管体是狭长规则几何体，优先采用四边形（quadrilateral）、六面体网格（he某ahedralcell），它们相比较其他类型网格允许较大的纵
横比。

几何体弯曲部分要给以加密，修正由弯曲造成的网格放大。

2）由
传热及流动过程界定网格布局：
首先，我们想到的是弯管的内壁，此处烟气和管壁发生强制对流换热。

温度变化很大，所以需要较密的网格。

可以通过划分边界层网格的手段得
到沿半径方向成比例渐变的网格。

如图a
a.边界层网格划分
其次，考虑到流体进入拐角后流动情况复杂，也有必要对该处的网格
加密。

需要先对线来进行网格划分，即由控制线上节点的数量来对拐角处
体网格疏密进行控制。

注：该线需是弯管面沿轴向的分割线，可以通过plitface命令将弯面沿着轴向方向虚拟分割。

再利用cooper命令画出的
弯管部分的网格如图b
b.中间段网格
最后给各个面设定边界条件并检查网格质量。

注:若发现负体积需要
重新划分网格改变网格布局。

因为负体积将直接导致基于有限体积法的fluent得出错误结果。

（2）单元和节点统计
25440he某ahedralcell,zone2,binary.96404he某
ahedralcell,zone3,binary.
281487quadrilateralinteriorface,zone20,binary.12320quadrilat eralwallface,zone1,binary.12320quadrilateralwallface,zone23,bina ry.560quadrilateralwallface,zone4,binary.5320quadrilateralwallfa ce,zone5,binary.2600quadrilateralwallface,zone6,binary.4000quadr ilateralwallface,zone7,binary.160quadrilateralwallface,zone8,bin ary.160quadrilateralwallface,zone9,binary.313quadrilateralinteri
orface,zone10,binary.80quadrilateralinteriorface,zone11,binary.3 13quadrilateralinteriorface,zone12,binary.80quadrilateralinterio rface,zone13,binary.313quadrilateralinteriorface,zone14,binary.8 0quadrilateralinteriorface,zone15,binary.313quadrilateralinterio rface,zone16,binary.80quadrilateralinteriorface,zone17,binary.31 3quadrilateraloutflowface,zone18,binary.
313quadrilateralvelocity-inletface,zone19,binary.
63440quadrilateralinteriorface,zone21,binary.12320hadowfacep air,binary.141166node,binary.
141166nodeflag,binary.（3）网格单元质量
Equiizekew质量类型时最差网格质量值为0.40125EquiAnglekew质量类型时最差网格质量值也为0.40125可见网格质量很好。

４、计算模型描述
高质量的网格为随后的计算提供了有力保障，但仅这样还不够。

针对不同模拟对象我们还需要选择不同的计算方式，这就需要对计算模型进行描述。

其中主要是边界条件、湍流模型、近壁处理。

(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弯曲段对流换热系
数elbow_conv取前两者的平均入口命名为inlet，取速度入口V=150.1m/ m2k；
m2k；
14.04W内壁命名为inner_wall，对流换热系数值为257.7W管子外表面辐射率和吸收率均取0.8
m2k
(3)湍流模型
目前湍流模型有四种：Spart-Alpla单方程模型；K-e双方程模型；
雷诺应力模型；大涡模拟模型。

单方程模型在这几种模型的中的计算量最小，它是一种刚刚发展起来
的湍流模型，主要针对于航空流体机械的数值模拟，对于其他复杂流动的
计算还没有经过验证双方程湍流模型能够比较准确地模拟各种复杂流动，
而且计算量也在工程可以接受的范围；标准k-ε模型解决一般的流动问题，RNGk-ε模型主要应用于旋转坐标系下的流动问题（旋转机械），Realizablek-ε模型主要用于射流、大分离、回流等问题。

雷诺应力模型和大涡模拟主要应用于湍流运动的机理研究中，由于计算量非
常大，因此目前还很少用于有复杂几何形状的工程问题中网格类型的选择。

本题采用工程中较常用的双方程k-ε模型来进行湍流模拟。

（4）近
壁处理
对于粘性流动问题，FLUENT默认设置是壁面无滑移条件，但你也可
以指定壁面切向速度分量（壁面平移或者旋转运动时），也可以给出壁面
切应力从而模拟壁面滑移。

根据当地流动情况，可以计算壁面切应力和与
流体换热情况。

湍流流动在近壁区域受壁面的影响很大。

５、求解过程简报（1）导入网格
１）启动fluent，选择3D模式；
２）Grid=>Check读入网格数据并检查，确保网格最小体积不能为负，同时留意其他“警告”提示。

必要时重新划分网格。

３）Grid=>Scale根据题目要求确定长度单位为mm；
４）Diplay=>Grid显示网格，再次确认之前在gambit里定义的边界
条件。

（2）定义模型及边界
１）Define=>Model=>Solve定义求解器
在gradientoption项里选择LeatSquareCellBaed，其他项默认不变。

２）Define=>Model=>Energy
选择能量方程并设置湍流模型，在Model项里选择k-epilon(2eqn)，点OK３）Define=>Material定义流体属性
依次填入700K烟气的Denity（密度）、Cp（定压比热容）、ThermalConductiviy（导热系数）、Vicoity（动力粘度）。

点
Change/Create完成设置。

４）Define=>BoundaryCondition定义边界条
件
依据前面的计算结果输入相应边界条件。

注意，入口面速度默认是沿
某轴正向的，随几何模型不同入口朝向也各不相同，应根据实际情况定义。

（3）求解
２）设置监视窗口
Solver=>Monitor=>Reidual定义残差，在Option栏选取Plot，点
OK
Solver=>Monitor=>SurfaceMonitor定义监视面，这里取出口为监视面，监视出口的Turbulence（湍流度）以及Ma-Weighted-Average（平均
比重）File=>Write=>Autoave设置自动保存，这里取每20次迭代自动保
存一次Solver=>Iterate...设置迭代次数500开始迭代，当迭代至199
次时达到收敛（这里未对残差进一步设置，所以较快就收敛了）残差迭代的图
出口流量监视图。