搜索
一两端带法兰弯管置于大空间内，管外壁与空气发生自然对流换热；内通烟气并与管内壁发生强制对流换热。结构和尺寸及其它条件如图。计算任务为用计算流体力学/计算传热学软件Fluent求解包括管内流体和管壁固体在内的温度分布，其中管壁分别采用薄壁和实体壁两种方法处理。所需的边界条件采用对流换热实验关联式计算。要求在发动机数值仿真实验室的计算机上完成建立几何模型、生成计算网格、建立计算模型、提交求解、和结果后处理等步骤，并分别撰写计算任务的报告，计算报告用计算机打印。计算报告包括以下与计算任务相关的项目和内容：(1)...............................传热过程简要描述包括传热方式、流动类型等；(2)计算方案分析包括所求解的控制方程及其简化、边界条件及其确定方法和主要计算过程；(3)计算网格简报包括网格划分方案、单元拓扑、单元和节点数量、网格质量等；(4)计算模型描述包括流体物性、边界条件、湍流模型、辐射模型及近壁处理等；(5)求解过程简报包括求解方法、离散格式、迭代过程监控、收敛准则等；(6)...............................................计算结果及分析给出下列图表和数据：w纵剖面和中间弯管45°方向横剖面上的温度、温度梯度、速度分布图，以及法兰和中间弯管处的局部放大图。w管内壁面上的温度、热流密度和表面传热系数分布，包括三维分布和沿管长度方向上的分布。w..................................................................................................总热流量。w由2种数值计算方法求得管内外烟气和空气之间换热的平均传热系数和烟气出口温度，并与工程算法得到的数值对比。1、传热过程简述计算任务为用计算流体力学/计算传热学软件Fluent求解通有烟气的法兰弯管包括管内流体和管壁固体在内的温度分布，其中管壁分别采用薄壁和实体壁两种方法处理。在进行分析时要同时考虑导热，对流，热传递三种传热方式。（1）、直角弯管内外壁间的热传导。注意：如果是按薄壁方法处理时不用考虑该项，此时管壁厚度可以忽略，内壁和外壁温差几乎为零。（2）、管子外壁与环境发生的自然对流换热。由于流体浮生力与粘性力对自然对流的影响，横管与竖管对流换热系数略有不同的。计算公式也不一样。同时，还有内壁同烟气发生的强制对流换热。（3）、管子外壁和大空间（环境）的辐射换热通过进气温度和流量，我们可以推断出管内烟气为湍流流动。这在随后的计算中可以得到证实。2、计算方案分析（1）、控制方程及简化质量守恒方程：()()()0uvwtxyzρρρρ∂∂∂∂+++=∂∂∂∂它表述的是单位时间内流体微元体中质量的增加，等于同一时间间隔内流入该微元体的净质量。对定常不可压缩流动该方程可简化为0uvwxyz∂∂∂++=∂∂∂动量守恒方程：()()()uupdivUudivgraduStxρρµ∂∂+=+−∂∂()()()vvpdivUvdivgradvStyρρµ∂∂+=+−∂∂()()()wwpdivUwdivgradwStzρρµ∂∂+=+−∂∂该方程的依据是，微元体中的流体的动量对时间的变化率等于外界作用在该微元体上的各种力之和。式中uS、vS、wS是动量方程中的广义源项。和前面一样上式可以简化为：()()1updivUudivgradutxνρ∂∂+=−∂∂()()1vpdivUvdivgradvtyνρ∂∂+=−∂∂()()1wpdivUwdivgradwtzνρ∂∂+=−∂∂能量守恒方程其依据四能量守恒定律：微元体中能量的增加率等于进入微元体的净热流量加上体力与面力对微元体所做的功。()()()()()hhuhvhwhpdivUdivgradTStxyzρρρρλ∂∂∂∂+++=−++Φ+∂∂∂∂λ：流体导热系数；hS：流体内热源Φ：由于粘性作用，机械能转换为热能的部分；对本题可化为：()TpSTdivUTdivgradTtcλρρ⎛⎞∂+=+⎜⎟⎜⎟∂⎝⎠（2）边界条件的定义1）发生在管内、外壁上的对流换热属于第三类边界条件2）弯管两端法兰截面是绝热的，属于第二类边界条件（热流量为零）3）外表面的辐射换热边界条件1）内壁上强制对流换热系数1h查阅资料[]确定700K烟气的热物性参数如下30.509/kgmρ=；()25.9*10/*Wmkλ−=；6267.74*10/msν−=；1159/*CpJkgK=又已知：烟气流量0.15/mkgs=，管的内直径500.05dmmm==由mAVρ=；A为流道截面积，Ｖ为烟气流速20.15150.1/0.050.5092mVmsAρπ===⎛⎞⎜⎟⎝⎠，烟气比热容比1.33γ=，287.52j/kgKgR=当地声速gcRTγ=，则1.33287.52700517.38/cms=××=，则150.10.29517.38VMc===，当0.4M时，气体可是为不可压流体而雷诺数56150.10.051.111067.7410eVdRν−×===××，（管内烟气为湍流流动）管内湍流换热实验关联式实用上使用最广的是迪贝斯－贝尔特公式：0.80.023RePrnfffNu=加热流体时=0.4n，冷却流体时=0.3n。式中:定性温度采用流体平均温度，特征长度为管内径。将数据代入上式得()0.850.30.0231.11100.637218.39fNu=×=得到强制对流换热系数2215.910218.39257.70.05WhNumkdλ−×==×=2）外壁自然流换热系数2h、3h自然对流试验关联式：()PrnNuCGr=，hlNuλ=，83.5110Gr=×其中0.59C=、1/4n=、Pr0.701=【查表】，注意：对于竖管段特征长度取管长，横管特征长度取管直径。定性温度采用()()/2700300/2500mwtttK∞=+=+=，700300400wftttK∆=−=−=查表知500K时干空气热物性参数：30.705/kgmρ=、()24.1*10/*WmKλ−=、6237.83*10/msν−=得竖管对流换热系数()13422PrPrngtlhNuCGrClllλλαλν⎛⎞∆===⎜⎟⎝⎠，代数数据得()13422229.84000.44.1100.590.7014.47*5000.4Whmkν−⎛⎞×××==⎜⎟⎜⎟×⎝⎠同理得横管段77.1110Gr=×，则c=0.48,1/4n=()1342329.84000.070.480.70123.65000.07Whmkλν⎛⎞××==⎜⎟⎜⎟×⎝⎠注：横管的C=0.48（必须知道2hNulλ=中l的含义，以及32gtlGrαν∆=中l的含义）3、网格简报本模型用UG进行建模，采用GAMBIT来划分网格。一般说来，网格越密越能保证精确的计算结果。但我们也要具体问题具体分析，针对不同模型采用不同的处理方法。在硬件配置允许的条件下，画出合理的高质量的网格。（1）网格划分方案1）管体是狭长规则几何体，优先采用四边形（quadrilateral）、六面体网格（hexahedralcells），它们相比较其他类型网格允许较大的纵横比。几何体弯曲部分要给以加密，修正由弯曲造成的网格放大。2）由传热及流动过程界定网格布局：首先，我们想到的是弯管的内壁，此处烟气和管壁发生强制对流换热。温度变化很大，所以需要较密的网格。可以通过划分边界层网格的手段得到沿半径方向成比例渐变的网格。如图aaaaa....边界层网格划分其次，考虑到流体进入拐角后流动情况复杂，也有必要对该处的网格加密。需要先对线来进行网格划分，即由控制线上节点的数量来对拐角处体网格疏密进行控制。注：该线需是弯管面沿轴向的分割线，可以通过splitface命令将弯面沿着轴向方向虚拟分割。再利用cooper命令画出的弯管部分的网格如图bb.b.b.b.中间段网格最后给各个面设定边界条件并检查网格质量。注:若发现负体积需要重新划分网格改变网格布局。因为负体积将直接导致基于有限体积法的fluent得出错误结果。（2）单元和节点统计25440hexahedralcells,zone2,binary.96404hexahedralcells,zone3,binary.281487quadrilateralinteriorfaces,zone20,binary.12320quadrilateralwallfaces,zone1,binary.12320quadrilateralwallfaces,zone23,binary.560quadrilateralwallfaces,zone4,binary.5320quadrilateralwallfaces,zone5,binary.2600quadrilateralwallfaces,zone6,binary.4000quadrilateralwallfaces,zone7,binary.160quadrilateralwallfaces,zone8,binary.160quadrilateralwallfaces,zone9,binary.313quadrilateralinteriorfaces,zone10,binary.80quadrilateralinteriorfaces,zone11,binary.313quadrilateralinteriorfaces,zone12,binary.80quadrilateralinteriorfaces,zone13,binary.313quadrilateralinteriorfaces,zone14,binary.80quadrilateralinteriorfaces,zone15,binary.313quadrilateralinteriorfaces,zone16,binary.80quadrilateralinteriorfaces,zone17,binary.313quadrilateraloutflowfaces,zone18,binary.313quadrilateralvelocity-inletfaces,zone19,binary.63440quadrilateralinteriorfaces,zone21,binary.12320shadowfacepairs,binary.141166nodes,binary.141166nodeflags,binary.（3）网格单元质量Equisizeskew质量类型时最差网格质量值为0.40125EquiAngleskew质量类型时最差网格质量值也为0.40125可见网格质量很好。４、计算模型描述高质量的网格为随后的计算提供了有力保障，但仅这样还不够。针对不同模拟对象我们还需要选择不同的计算方式，这就需要对计算模型进行描述。其中主要是边界条件、湍流模型、近壁处理。(1)流体物性查资料得700K烟气的热物性如下30.509/kgmρ=；()25.9*10/*Wmkλ−=；6267.74*10/msν−=；1158.8/*CpJkgK=(2)边界条件1）法兰管端面为绝热，热流量为零。定义名称为flux_1、flux_22）由于横管竖管自然对流换热系数不同，需要分别定义：竖管对流换热系数命名为up_conv，值为24.47Wmk⋅；横管对流换热系数命名为down_conv，值为223.6Wmk⋅；弯曲段对流换热系数elbow_conv取前两者的平均214.04Wmk⋅；入口命名为inlet，取速度入口V=150.1m/s内壁命名为inner_wall，对流换热系数值为2257.7Wmk⋅管子外表面辐射率和吸收率均取0.8(3)湍流模型目前湍流模型有四种：Spart-Alpla单方程模型；K-e双方程模型；雷诺应力模型；大涡模拟模型。单方程模型在这几种模型的中的计算量最小，它是一种刚刚发展起来的湍流模型，主要针对于航空流体机械的数值模拟，对于