《热对流方程》课件

热对流方程欢迎大家学习热对流方程课程热对流是自然界和工程应用中最常见的传热方式之一，它结合了流体流动和热传导两种物理过程通过本课程，我们将深入探讨热对流的基本理论、数学描述以及工程应用热对流方程是描述流体流动与热量传递耦合现象的数学模型，对于理解从大气环流到工业换热器等各种现象至关重要在接下来的课程中，我们将系统地学习这一领域的核心概念和应用方法课程目标掌握热对流基本理论理解热对流的物理本质、基本概念和分类，掌握描述热对流过程的基本物理量及其关系熟悉热对流方程推导学习质量、动量和能量守恒定律，理解如何从这些基本定律推导出热对流的控制方程组应用热对流计算方法掌握热对流方程的数值解法，了解实际工程中的应用技巧，能够运用所学知识解决实际工程问题了解前沿研究方向了解热对流领域的最新研究进展和未来发展趋势，为进一步深入学习和研究打下基础热对流的基本概念定义机理热对流是指流体流动过程中伴流体宏观流动携带热量对流随的热量传递现象它结合了输运与流体微观分子运动传热传导和流体流动两种物理过递热量传导共同作用的结果程，是一种更为复杂的传热方流体温度梯度导致密度差异式，进而影响流动状态特点热对流传热效率通常高于纯导热；热对流过程通常伴随复杂的流体动力学现象；温度场和速度场相互耦合，形成复杂的非线性系统热对流的类型自然对流强制对流由于温度差异引起流体密度变化，在重力场作用下产生的流在外部动力（如泵、风机等）驱动下，流体强制流动过程中体运动和热量传递发生的热量传递例如热空气上升、冷空气下沉；热水上浮、冷水下沉；自例如风扇冷却计算机CPU；暖气片加热室内空气；换热器然风形成等现象中的流体换热等特点无需外部动力驱动；流速较低；受布斯涅斯克近似控特点需要外部动力驱动；流速较高；雷诺数是关键参数；制；格拉晓夫数是关键参数换热效率通常高于自然对流热对流过程中的物理量温度T速度v压力p表征物体热状态的物理量，流体质点的运动速率和方向流体单位面积上的法向力，是热对流研究的核心物理量，通常表示为三维矢量场单位为帕斯卡Pa压力梯之一通常采用开尔文K、单位为米/秒m/s速度场度是驱动流体流动的主要因摄氏度℃或华氏度℉为单描述了流体的运动状态，是素之一，与速度场密切相关位温度梯度是导致热量传研究热对流的另一个核心物递的根本驱动力理量热流密度q单位时间内通过单位面积的热量，单位为瓦特/平方米W/m²是表征热量传递强度的重要物理量，与温度梯度和对流换热系数有关流体力学基础回顾纳维-斯托克斯方程1描述流体运动的基础方程欧拉方程2忽略粘性的流体运动方程连续性方程3质量守恒状态方程4流体热力学性质方程热对流方程的推导基于流体力学的基本理论要全面理解热对流现象，必须先掌握流体力学的基础知识，包括各种守恒定律和控制方程这些方程描述了流体的运动规律和特性，是研究热对流问题的理论基础流体力学与热传递的结合形成了热对流的复杂过程流体动力学方程组与热传递方程组耦合，构成了完整的热对流方程体系连续性方程1物理意义2微分形式连续性方程表达了流体质量守∂ρ/∂t+∇·ρv=0，其中ρ为恒的原理，即在流动过程中，流体密度，v为流体速度矢量流体的质量既不会凭空产生，，t为时间，∇·表示散度算子也不会凭空消失对于任意控这个方程适用于任何流体（制体积，进入的流体质量等于可压缩或不可压缩）的流动情流出的流体质量加上控制体积况内流体质量的变化率3不可压缩流体简化形式当流体可视为不可压缩时（如液体或低速气体流动），密度为常数，ρ连续性方程简化为∇·v=0，这表明不可压缩流体的速度场是无散度的动量方程（纳维斯托克斯方程）-物理意义微分形式纳维-斯托克斯方程表达了流体动ρ∂v/∂t+v·∇v=-∇p+μ∇²v+量守恒的原理，即牛顿第二定律ρg，其中ρ为流体密度，v为流体在流体系统中的应用它描述了速度矢量，p为压力，μ为动力粘流体质点在各种力（压力梯度力度，g为重力加速度左侧表示流、粘性力、重力等）作用下的运体质点的加速度（局部加速度和动规律对流加速度），右侧表示作用在流体上的各种力热对流中的特点在热对流问题中，温度变化会引起流体密度变化，进而影响动量方程通常采用布西内斯克近似，即仅在浮力项中考虑密度变化，表示为ρg≈ρ₀[1-βT-T₀]g，其中β为体积膨胀系数能量方程物理意义1能量方程表达了热力学第一定律，即能量守恒原理在流体系统中的应用它描述了流体内能的变化与热量传递和功的关系，是热对流微分形式2问题的核心方程之一ρcp∂T/∂t+v·∇T=k∇²T+Φ，其中ρ为流体密度，cp为定压比热容，T为温度，k为热导率，Φ为粘性耗散项左侧表示流体单位体积边界条件3内能的变化率，右侧表示热传导和粘性耗散产生的热量求解能量方程需要适当的边界条件，包括温度边界条件（第一类边界条件）；热流密度边界条件（第二类边界条件）；对流换热边界条件（第三类边界条件）；辐射换热边界条件等边界层理论简介边界层定义边界层特性1流体在固体表面附近形成的薄层区域，其中厚度随雷诺数增大而减小，内部流动状态可2速度和温度梯度较大为层流或湍流边界层影响边界层方程4决定了表面阻力和热传递效率，是热对流研简化的纳维-斯托克斯方程，省略了垂直于3究的核心问题壁面方向的项边界层理论由普朗特于1904年提出，它极大地简化了流体力学问题的求解在边界层内，流体的粘性影响显著，而在边界层外，流体可以视为无粘理想流体这一理论为解决实际工程中的流动和传热问题提供了有效的理论框架在热对流问题中，边界层理论尤为重要，因为大部分热量传递都发生在边界层内准确理解和计算边界层特性是解决热对流问题的关键速度边界层边界层厚度δ速度边界层厚度定义为流体速度达到主流速度的99%处的垂直距离平板上的层流边界层厚度近似公式δ≈

5.0x/√Re_x，其中x为沿流动方向的距离，Re_x为局部雷诺数边界层内速度分布层流边界层内的速度分布可用三次多项式近似u/U∞=2y/δ-y/δ²，其中u为x方向速度，U∞为主流速度，y为垂直于壁面的距离湍流边界层内的速度分布更为复杂，通常采用对数律描述壁面剪切应力壁面剪切应力τw=μ∂u/∂yy=0，它与边界层内速度梯度直接相关，决定了流体对固体表面的摩擦阻力壁面剪切应力系数Cf=τw/½ρU∞²是表征边界层特性的重要无量纲参数温度边界层温度边界层定义与速度边界层的关系壁面热流密度当流体流过温度与主流不同的物体表面温度边界层与速度边界层同时发展，但壁面热流密度qw=-k∂T/∂yy=0，它与时，在表面附近形成的温度梯度较大的二者厚度通常不同，其比值δt/δ受普朗温度边界层内温度梯度直接相关，决定区域其中温度从壁面温度Tw变化到主特数Pr控制Pr1时，δtδ；Pr1了固体表面与流体之间的热量传递通流温度T∞的区域温度边界层厚度δt时，δtδ；Pr=1时，δt=δ水的Pr过引入对流换热系数h，可以表示为qw定义为温度变化达到主流温度99%处的约为7，空气的Pr约为

0.7=hTw-T∞，这就是牛顿冷却定律垂直距离无量纲参数介绍普朗特数Pr雷诺数Re2表征动量扩散与热扩散的比值，影响速度与温表征惯性力与粘性力的比值，决定流动状态（度边界层的相对厚度1层流或湍流）努塞尔数Nu3表征对流换热与纯传导换热的比值，是评价换热效果的关键参数瑞利数Ra5格拉晓夫数Gr格拉晓夫数与普朗特数的乘积，表征自然对流的强度表征浮力与粘性力的比值，是自然对流的重要4参数无量纲参数在热对流问题中具有重要意义，它们可以简化分析过程，使研究结果具有普遍适用性通过相似性原理，可以将实验室小尺度模型的结果推广到实际工程大尺度应用中不同的无量纲参数反映了流体流动和热传递过程中不同物理机制的相对重要性，是热对流研究和工程应用中的重要工具雷诺数（）Re1定义2物理意义雷诺数是惯性力与粘性力的比值雷诺数表征了流体流动中惯性效，定义为Re=ρvL/μ=vL/ν，应与粘性效应的相对重要性雷其中ρ为流体密度，v为特征速度诺数较小时，粘性力起主导作用，L为特征长度，μ为动力粘度，，流动呈层流状态；雷诺数较大为运动粘度雷诺数由英国物时，惯性力占主导，流动趋于湍ν理学家奥斯本·雷诺（Osborne流状态雷诺数是判断流动状态Reynolds）于1883年首次提出的重要依据3临界值不同流动情况有不同的临界雷诺数管道内流动的临界雷诺数约为2300；平板上的临界雷诺数约为5×10⁵当实际雷诺数超过临界值时，流动从层流转变为湍流湍流状态下的换热效率通常高于层流普朗特数（）Pr1定义2物理意义普朗特数是动量扩散与热扩散普朗特数反映了流体传递动量的比值，定义为Pr=ν/α=能力与传递热量能力的相对强μcp/k，其中ν为运动粘度，α弱Pr1表示热扩散快于动量为热扩散率，μ为动力粘度，扩散，Pr1表示动量扩散快于cp为定压比热容，k为热导率热扩散普朗特数决定了速度普朗特数由德国物理学家路边界层与温度边界层的相对厚德维希·普朗特（Ludwig度关系Prandtl）提出3常见流体的Pr值液态金属Pr≈

0.01-

0.05（热传导很强）；气体Pr≈

0.7-

1.0；水（20℃）Pr≈7；油类Pr≈50-2000不同普朗特数的流体在换热器设计中有不同的适用场景和优化方向努塞尔数（）Nu定义经验关联式努塞尔数是对流换热与纯传导换热的比值，定义为Nu=强制对流Nu=C·Re^m·Pr^nhL/k，其中h为对流换热系数，L为特征长度，k为流体热导自然对流Nu=C·Gr·Pr^m=C·Ra^m率努塞尔数由德国工程师威廉·努塞尔（Wilhelm Nusselt）提出其中C、m、n为经验常数，与具体几何形状和流动状态有关这些关联式是工程计算中最常用的工具，通过实验数据拟努塞尔数本质上表示了实际对流换热效果相当于纯传导的倍合得到数，反映了对流强化热传递的程度Nu=1表示纯传导，无对流效应；Nu1表示对流增强了热传递例如，平板层流强制对流的关联式Nu_x=

0.332·Re_x^1/2·Pr^1/3，适用于Re_x5×10⁵,Pr

0.6格拉晓夫数（）Gr定义1浮力与粘性力的比值计算公式2Gr=gβTs-T∞L³/ν²参数含义3g-重力加速度，β-体积膨胀系数，L-特征长度，ν-运动粘度格拉晓夫数（Grashof number）是自然对流中的关键无量纲参数，由德国工程师弗兰兹·格拉晓夫（Franz Grashof）提出它表征了流体因温度差异引起的浮力效应与粘性阻力的相对大小，类似于强制对流中雷诺数的作用格拉晓夫数越大，浮力效应越显著，自然对流越强烈当Gr较小时，流动呈层流状态；当Gr超过特定临界值时，流动转变为湍流对于垂直平板，临界格拉晓夫数约为10⁹在自然对流换热的关联式中，通常使用格拉晓夫数与普朗特数的乘积（即瑞利数Ra=Gr·Pr）来表征换热效果热对流方程的推导概述基本假设牛顿流体、连续介质、热物性参数恒定（除密度外）、忽略粘性耗散和压缩功、布西内斯克近似等这些假设简化了问题，使方程更易处理，同时保留了主要物理特性守恒定律质量守恒（连续性方程）、动量守恒（牛顿第二定律）、能量守恒（热力学第一定律）这三个基本守恒定律是热对流方程组的理论基础推导方法控制体积分析法对微小控制体积应用守恒定律，导出微分方程也可采用场论方法，直接建立守恒量的输运方程两种方法在数学上等价，但物理解释略有不同完整方程组包括连续性方程、动量方程（纳维-斯托克斯方程）和能量方程，再加上状态方程和边界条件，形成完整的热对流方程组这是求解热对流问题的数学基础质量守恒方程1推导过程2向量形式考虑一个微小控制体积，在连续性方程的向量形式为时间间隔dt内，进入控制体∂ρ/∂t+∇·ρV=0，其中积的质量减去流出控制体积V=u,v,w是速度向量，∇·的质量，等于控制体积内质表示散度算子这种形式更量的增加量对于三维流动加简洁，适用于任何坐标系，表达为∂ρ/∂t+统中的表达∂ρu/∂x+∂ρv/∂y+∂ρw/∂z=03不可压缩流体简化对于不可压缩流体，密度ρ为常数，方程简化为∇·V=0，或分量形式∂u/∂x+∂v/∂y+∂w/∂z=0这一简化形式在液体流动和低速气体流动（马赫数Ma

0.3）的热对流问题中广泛适用动量守恒方程推导基础纳维-斯托克斯方程动量守恒方程基于牛顿第二定律F=ma，即控制体积内流体x方向ρ∂u/∂t+u∂u/∂x+v∂u/∂y+w∂u/∂z=-∂p/∂x+所受的合外力等于其动量变化率考虑各个方向上的力平衡μ∂²u/∂x²+∂²u/∂y²+∂²u/∂z²+ρgₓ，包括表面力（压力和粘性力）和体积力（如重力）y方向ρ∂v/∂t+u∂v/∂x+v∂v/∂y+w∂v/∂z=-∂p/∂y+对于牛顿流体，粘性应力与变形速率成正比，比例系数为动μ∂²v/∂x²+∂²v/∂y²+∂²v/∂z²+ρgᵧ力粘度μ这一构成关系是推导动量方程的关键步骤z方向ρ∂w/∂t+u∂w/∂x+v∂w/∂y+w∂w/∂z=-∂p/∂z+μ∂²w/∂x²+∂²w/∂y²+∂²w/∂z²+ρgz能量守恒方程热力学第一定律流体内能方程能量守恒方程基于热力学第一ρ∂e/∂t+u∂e/∂x+v∂e/∂y+定律，即系统内能的增加等于w∂e/∂z=-p∂u/∂x+∂v/∂y+外界对系统做功加上系统吸收∂w/∂z+k∂²T/∂x²+∂²T/∂y²+的热量对于流体微元，需考∂²T/∂z²+Φ，其中e为单位质虑对流热传递、热传导、内能量内能，为粘性耗散函数Φ变化以及各种功的作用流体温度方程对于不可压缩流体，假设比热容cp为常数，引入T=e/cp，可得ρcp∂T/∂t+u∂T/∂x+v∂T/∂y+w∂T/∂z=k∂²T/∂x²+∂²T/∂y²+∂²T/∂z²+Φ，这是热对流问题中常用的能量方程形式热对流方程的基本形式完整方程组布西内斯克近似热对流的完整控制方程组包括在大多数热对流问题中，密度变化仅在浮力项中考虑，即•连续性方程∂ρ/∂t+∇·ρV=0ρ=ρ₀[1-βT-T₀]，其中β为体积膨胀系数•动量方程ρ∂V/∂t+V·∇V=-∇p+μ∇²V+ρg•能量方程ρcp∂T/∂t+V·∇T=k∇²T+Φ这一近似使方程组显著简化，同时保留了温度引起浮力效应•状态方程ρ=ρp,T的关键特性浮力项表示为ρg≈ρ₀g[1-βT-T₀]布西内斯克近似适用于小温差问题βT-T₀≪1，如大多数自然对流问题热对流方程的简化二维情况对于二维稳态热对流问题，控制方程可大幅简化假设流动在x-y平面内，且物性参数恒定，应用布西内斯克近似，方程组可表示为连续性方程∂u/∂x+∂v/∂y=0x方向动量方程ρu∂u/∂x+v∂u/∂y=-∂p/∂x+μ∂²u/∂x²+∂²u/∂y²y方向动量方程ρu∂v/∂x+v∂v/∂y=-∂p/∂y+μ∂²v/∂x²+∂²v/∂y²-ρ₀gβT-T₀能量方程ρcpu∂T/∂x+v∂T/∂y=k∂²T/∂x²+∂²T/∂y²热对流方程的简化一维情况1一维流动情形2平行平板间稳态层流在某些特殊情况下，流动和对于平行平板间的充分发展温度分布可简化为一维问题层流，速度仅是y的函数，，如充分发展的管道流动、温度是x和y的函数动量方层状流动等这种简化大大程简化为0=-dp/dx+降低了问题的复杂性，使得μd²u/dy²，能量方程简化为解析解成为可能ρcpudT/dx=kd²T/dy²3垂直平板自然对流对于垂直平板附近的自然对流边界层，假设边界层很薄，压力与自由流体静压相同，则动量方程简化为u∂u/∂x+v∂u/∂y=ν∂²u/∂y²+gβT-T∞，能量方程简化为u∂T/∂x+v∂T/∂y=α∂²T/∂y²边界条件的设定速度边界条件温度边界条件壁面无滑移条件固体表面处流第一类边界条件（Dirichlet条件）体速度为零，即u=v=w=0自由指定表面温度T=Ts第二类边表面条件法向速度为零，切向界条件（Neumann条件）指定剪切应力为零入口条件指定表面热流密度-k∂T/∂n=qs第三流体流入速度分布出口条件类边界条件（Robin条件）指定通常假设流体充分发展，速度梯表面对流换热-k∂T/∂n=hTs-T∞度为零特殊边界条件周期性边界条件计算域边界上的物理量在空间上周期性重复对称性边界条件利用问题的对称性，简化计算域界面条件流体-流体或流体-固体界面上的速度连续、温度连续、热流连续等条件初始条件的设定时变问题的初始条件稳态问题的初始条件常见的初始条件类型对于非稳态问题，需要指定计算域内各点在虽然稳态问题理论上不需要初始条件，但在均匀初始场整个计算域内速度和温度取同初始时刻t=0的速度场Vx,y,z,0和温度场数值求解过程中常采用时间推进法逐步逼近一值分段初始场不同区域设定不同的初Tx,y,z,0初始条件的选择应尽量接近问题稳态解此时需要设置初始猜测值，这实际始值渐进初始场物理量按空间位置线性的物理实际，以减少计算时间和提高数值稳上就是一种初始条件初始猜测值的选择或非线性变化也可以使用简化模型的解作定性会影响计算的收敛速度为初始条件，或使用前一个工况的计算结果作为新工况的初始条件层流与湍流的区别层流特性湍流特性层流是指流体质点沿着流线有序、规则地流动，流体层与层湍流是指流体质点做无规则随机运动，流体层之间发生强烈之间不发生混合的流动状态层流的特点包括流动规律、混合的流动状态湍流的特点包括流动无规则、紊乱；流有序；流线光滑、稳定；流体质点轨迹可以准确预测；动量线扭曲、不稳定；流体质点轨迹不可精确预测；存在大小不、能量和质量的交换主要通过分子扩散实现；摩擦阻力与流同的湍流涡旋；动量、能量和质量的交换主要通过涡流混合速成正比实现；摩擦阻力与流速的平方成正比层流中的热传递效率相对较低，因为热量主要通过导热方式湍流中的热传递效率较高，因为涡流混合加强了热量在流体在流体中传播中的传播层流热对流方程1层流的数学描述2层流方程的特点层流状态下，流体质点运动满方程中不包含时间平均项和脉足纳维-斯托克斯方程，没有动项，物理量（速度、压力、额外的湍流效应需要考虑控温度等）直接代表实际值，没制方程组包括连续性方程、动有统计平均的概念非线性项量方程和能量方程，这些方程主要来自对流项（如u∂u/∂x）构成了一个确定性的偏微分方，这些非线性项在高雷诺数时程组会导致数值解不稳定3求解方法对于简单几何形状和边界条件，层流方程可以获得解析解例如，平行平板间的泊肃叶流动、圆管内的哈根-泊肃叶流动等复杂情况下需要采用数值方法求解，如有限差分法、有限体积法或有限元法湍流热对流方程湍流特性与挑战雷诺平均方程闭合问题湍流流动具有多尺度、随机性和强烈的在雷诺平均方法中，物理量分解为平均RANS方程组中的雷诺应力和湍流热通量三维性，直接数值模拟DNS需要极高的值和脉动值u=ū+u，T=T̄+T将需要通过湍流模型来闭合常见的闭合计算资源为了工程应用，通常采用统此代入原始方程并取时间平均，得到雷模型包括代数模型（如混合长度模型计平均方法，如雷诺平均RANS或大涡诺平均纳维-斯托克斯RANS方程和雷诺）；一方程模型（如Spalart-Allmaras模模拟LES来处理湍流问题平均能量方程这些方程中出现了雷诺型）；两方程模型（如k-ε模型、k-ω模型应力-ρuv和湍流热通量-ρcpTv等新）；雷诺应力模型RSM等项湍流模型简介代数模型涡粘性假设通过简单代数关系计算湍流粘性，如混合长2将湍流脉动引起的动量和热量传递效应等效度模型1为增强的粘性和导热性一方程模型求解湍流动能或湍流粘性的输运方程，如3Spalart-Allmaras模型雷诺应力模型5两方程模型直接求解雷诺应力的输运方程，避免涡粘性假设的局限性4同时求解两个湍流特性量的输运方程，如k-ε和k-ω模型湍流模型的选择需要根据具体问题特点而定两方程模型如标准k-ε模型在工程应用中最为广泛，它求解湍流动能k和湍流耗散率ε的输运方程，具有计算效率高、适用范围广的特点，但在强旋转流动、强曲率流动和强浮力流动等情况下精度有限湍流热通量通常通过湍流普朗特数概念与湍流动量通量联系起来-ρcpTv=μt/Prt·∂T̄/∂y，其中μt为湍流粘度，Prt为湍流普朗特数（通常取

0.9）强制对流换热湍流强制对流1换热效率最高，应用最广泛过渡区强制对流2层流向湍流转变的复杂区域层流强制对流3换热效率较低但流动阻力小强制对流是在外力作用下（如泵、风机等）引起的流体流动所产生的对流换热它的特点是流速较高，对流换热系数大，换热效率高强制对流的强度主要由雷诺数Re表征，雷诺数越大，对流换热效果越好强制对流换热系数通常通过努塞尔数Nu关联式计算Nu=C·Re^m·Pr^n，其中C、m、n为与具体几何形状和流动状态有关的常数对于不同的流动条件（层流、过渡区或湍流）和几何形状（如平板、圆管、环形通道等），有不同的关联式可供选择在工程应用中，为了提高换热效率，常采用各种强化换热技术，如增加表面粗糙度、使用肋片或扰流件、改变流道几何形状等平板层流强制对流边界层形成1流体流过平板前缘时，由于粘性作用，在平板表面附近形成速度边界层和温度边界层两种边界层厚度的比值由普朗特数决定δt/δ≈Pr^-1/3边界层发展2随着流体沿平板流动，边界层厚度逐渐增加，满足δ≈

5.0x/√Re_x在某个临界位置，当局部雷诺数Re_x超过5×10⁵时，边界层由层流转变为湍流局部换热系数3层流区局部努塞尔数Nu_x=

0.332·Re_x^1/2·Pr^1/3，适用于Re_x5×10⁵，Pr

0.6这表明局部换热系数沿流动方向逐渐减小，符合h_x∝x^-1/2的规律平均换热系数4整个平板的平均努塞尔数Nu_L=

0.664·Re_L^1/2·Pr^1/3，适用于全层流区域如果包含湍流区，则需分段计算或使用复合关联式圆管内层流强制对流入口区流动充分发展区流动流体进入圆管后，首先经历流动在充分发展区，速度分布不再随发展区（水动力入口区），速度轴向位置变化，呈抛物线形ur分布从均匀分布逐渐发展为抛物=2U_m[1-r/R²]，其中U_m为平线形分布水动力入口区长度约均流速温度分布的形态取决于为L_e≈

0.05·Re·D同时存在热入边界条件，常见的是壁面恒温条口区，温度分布从初始分布逐渐件或壁面恒热流条件发展为稳定分布换热特性壁面恒温条件下，充分发展区Nu_D=

3.66；壁面恒热流条件下，充分发展区Nu_D=

4.36入口区的换热系数高于充分发展区，可通过修正因子考虑入口效应的影响考虑全管长的平均Nu数，可用关联式Nu_D=

3.66+

0.0668·D/L·Re·Pr/1+

0.04·D/L·Re·Pr^2/3。