流体力学基本方程组总结

流体力学根本方程组总结流体力学根本方程组包括连续性方程、运动方程、组分质量守恒方程、能量方程、本构方程、状态方程及通用形式守恒方程。虽各相关文献都有介绍这些根本方程组，但多采用作者熟悉的方式表示，往往同一方程具有多种形式，而难于直观比照。以下内容是对文献报道的各种形式的总结和比照，并分析了它们之间的转化关系，以期彻底理解〔切实掌握微分方程中每一项的物理意义〕流体力学根本方程组的数学物理意义，为离散计算该方程组打下根底。1连续性方程根据文献ADDINEN.CITE<EndNote><Cite><Author>吴望一</Author><Year>1982</Year><RecNum>6</RecNum><record><rec-number>6</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">吴望一</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">流体力学〔上册〕</style></title></titles><dates><year>1982</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">北京大学出版社</style></publisher><urls></urls></record></Cite></EndNote>[1]连续性方程可由四种方法得到，分别为拉格朗日法下对有限体积和体积元应用质量守恒定律、在欧拉法下对有限体积应用质量守恒定律及在直角坐标系中直接应用质量守恒定律。L法有限体积分析取体积为，质量为的一定流体质点团，那么有：〔1〕因为速度散度的物理意义是相对体积膨胀率及密度随体导数，即：〔2〕〔3〕代入式〔1〕得〔4〕运用奥高定理〔5〕得〔6〕上式即是连续性方程的积分形式。假定被积函数连续，而且体积是任意选取的，由此可知被积函数必须等于零，即：〔7〕或〔8〕在直角坐标系中连续性方程为：〔9〕或〔10〕连续性方程〔10〕说明，密度变化〔随时间和位置〕等于密度和体积变形的乘积ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陈敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流体力学〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化学工业出版社</style></publisher><urls></urls></record></Cite></EndNote>[2]。L法体积元分析考虑质量为的体积元，对其用拉格朗日观点，根据质量守恒定律有：〔11〕〔12〕两边同除以，得〔13〕或写成〔14〕上式说明要维持质量守恒定律，相对体积变化率必须等于负的相对密度变化率。1.3E法有限体积分析着眼坐标空间，取空间中以面为界的有限体积，那么称面为控制面，为控制体。取外法线方向为法线的正方向，为外法线方向的单位矢量。考虑该体积内流体质量的变化，该变化主要以下两方面原因引起。第一，通过外表有流体流出或流入，单位时间内流出流入变化的总和为：〔15〕第二，由于密度场的不定常性〔注意，欧拉观点下空间点是固定的，密度的变化只由场的不定常性刻画〕，单位时间内体积的质量将变化，变化量为：〔16〕上述两者应相等，即〔17〕由于体积是任意的，且被积函数连续，那么〔18〕1.4E法直角坐标系分析单位时间内通过外表EFGH的通量为：通过外表ABCD的通量为：其他三对外表类似，另外，该控制体内质量的变化率为：那么〔19〕特殊情况下的连续性方程：定常态：不可压缩流体：2动量方程任取一体积为的流体，它的边界为。根据动量定理，体积中流体动量的变化率等于作用在该体积上的质量力和面力〔应力〕之和。单位面积上的面力，其中是二阶对称应力张量，所以不是通常指的在〔单位体积面元的法线方向〕方向的分量。单位质量上的质量力为。那么作用在该体积上的质量力和面力分别为〔20〕及〔21〕动量变化率为〔22〕上述动量变化率的表达式可有两种处理方法ADDINEN.CITE<EndNote><Cite><Author>吴望一</Author><Year>1982</Year><RecNum>6</RecNum><record><rec-number>6</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">吴望一</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">流体力学〔上册〕</style></title></titles><dates><year>1982</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">北京大学出版社</style></publisher><urls></urls></record></Cite></EndNote>[1]，如下〔1〕求解上式右边第二项内对体积元的随体导数，那么〔23〕〔2〕对动量变化率表达式右边第二项应用质量守恒定律〔24〕由上可得两种积分形式的动量方程，即〔25〕或〔26〕由上，动量方程的微分形式为：〔27〕或〔28〕微分方程中各项的物理意义为，表示单位体积上惯性力，为单位体积上的质量力，为单位体积上应力张量的散度，它是与面力等效的体力分布函数〔由奥高公式转化而来〕。在直角坐标系下以应力表示的运动方程可采取以下形式〔29〕或〔30〕这两种表达方式的等号左边实际只差了一个连续性方程，由根本微分公式〔31〕得〔32〕由连续性方程知〔33〕所以有〔34〕上述运动方程是以应力表示的粘性流体的运动方程，它们对任何粘性流体，任何运动状态都是适用的。但它没有反映出不同属性的流体受力后的不同表现。另外，方程数和未知量之数不等，运动方程有三个，加上连续性方程共四个，但未知量却有九个〔六个应力张量分量〔九个张量分量因对称关系减少为六个〕和三个速度分量〕，所以该方程组不封闭。为使该方程组可解，必须考虑应力张量和变形速度张量之间的关系〔将应力张量用速度分量表示出来〕，补足所需的方程ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陈敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流体力学〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化学工业出版社</style></publisher><urls></urls></record></Cite></EndNote>[2]。3本构方程本构方程是表征流体宏观性质的一种微分方程，它是表达流体粘性定律的应力张量和变形速度张量之间的关系。最简单的应力与应变之间的关系是牛顿流体作一维运动，即牛顿剪切定律：〔35〕要得到普遍意义上的广义牛顿定律需作一定假设，而首先应理解流体速度分解定理和变形速度张量。文献ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陈敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流体力学〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化学工业出版社</style></publisher><urls></urls></record></Cite></EndNote>[2]对速度分解定理虽作了较直观的描述和推导但不严格，而文献ADDINEN.CITE<EndNote><Cite><Author>吴望一</Author><Year>1982</Year><RecNum>6</RecNum><record><rec-number>6</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">吴望一</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">流体力学〔上册〕</style></title></titles><dates><year>1982</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">北京大学出版社</style></publisher><urls></urls></record></Cite></EndNote>[1]对该局部内表达较详细。3.1速度分解定理刚体运动包括平动和转动两局部，一般可表为〔36〕其中是刚体中选定一点上的平动速度，是刚体绕点转动的瞬时角速度矢量，就是要确定速度那一点到点的矢径。转动角速度可用表示〔37〕故〔38〕流体运动除平动、转动外还有变形运动。设微团内点的速度为，邻域内任一点的速度为。将在点泰勒展开并略去二阶无穷小项，得〔39〕显然，是一个二阶张量〔局部速度梯度张量〕，由张量分解定理可将该张量分解成对称张量和反对称张量之和，于是〔40〕所以〔41〕上式右边第二、三项可具体表示为〔42〕及〔43〕其中〔44〕另外；所以〔45〕上式说明流体运动可分为平动、转动和变形三种形式组成，称为变形速度张量，该定理称为亥姆霍兹〔Helmholtz〕速度分解定理。另外，流变学中常用应变速率张量来表示流体的变形和拉伸〔或压缩〕，而用转动张量表示转动，它们与流体力学中的变形速度张量和转动张量的关系是：3.2变形速度张量的物理意义写出变形速度的表达式〔46〕经分析可得；其中，及是角变形速率，亦称剪切应变速率〔称拉伸应变速率〕。文献ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陈敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流体力学〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化学工业出版社</style></publisher><urls></urls></record></Cite></EndNote>[2]中定义的剪切速率的值是这里的一半，这是有问题的，因为剪切速率本身的值应以这里为准，但变形速度张量内剪切变形的量值为该剪切速率的一半。由上可知，变形速度张量的对角线分量，，的物理意义分别是轴线上线段元的相对拉伸速度或相对压缩速度。而非对角线分量的物理意义分别是与轴、与轴、与轴之间夹角的剪切速率的负值。3.3广义牛顿定律及根本假设〔1〕运动流体的应力张量在运动停止后应趋于静止流体的应力张量。据此将应力张量写成各向同性局部和各向异性局部是方便的。〔47〕是除去后得到的张量，称为偏应力张量。当运动消失时它趋于零。可见，偏应力张量和应力张量一样也是对称张量。〔2〕偏应力张量的各分量是局部速度梯度张量各分量的线性齐次函数。当速度在空间均匀分布时，偏应力张量为零；当速度偏离均匀分布时，在粘性流体中产生了偏应力，它力图使速度回复到均匀分布情形。〔3〕流体是各向同性的，即流体性质不依赖于方向或坐标系的转换。根据假设〔2〕，有〔48〕显然是一四阶张量，它是表征流体粘性的常数，共个。根据假设〔3〕，是各向同性张量且对，对称，故〔49〕观察上式可知，对也是对称的，物性常数减少至只有2个即第二粘度和粘度，证明见下。将上式代入偏应力表达式〔反对称项为零〕得〔50〕那么应力张量为〔51〕引进〔52〕那么〔53〕根据上式〔54a〕〔54b〕〔54c〕将上三式等号两边相加，得〔55〕对可压缩流体，流体的体积在运动过程中发生膨胀或收缩，它将引起平均法应力〔由奥高公式可证某固定点处所有方向上法应力的平均值等于三个方向上法应力的平均值，这是一个不随坐标系改变的不变量〕的值发生的改变，称为第二粘性系数亦称膨胀粘性系数。应用斯托克斯假定，即，那么本构方程为〔56〕〔57〕〔58〕一般处理的是不可压缩流体，那么〔59〕〔60〕〔61〕在直角坐标系下有〔62〕这里，有的文献中将应力张量用表示。将上述应力张量与变形速度张量的关系式代入运动方程，得即〔63〕写成直角坐标系下的形式〔64〕在数值传热学中ADDINEN.CITE<EndNote><Cite><Author>陶文铨</Author><Year>2001</Year><RecNum>8</RecNum><record><rec-number>8</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">陶文铨</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">数值传热学〔第二版〕</style></title></titles><dates><year>2001</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">西安</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">西安交通大学出版社</style></publisher><urls></urls></record></Cite></EndNote>[3]或CFD计算ADDINEN.CITE<EndNote><Cite><Author>王福军</Author><Year>2004</Year><RecNum>9</RecNum><record><rec-number>9</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">王福军</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">计算流体动力学——</style><styleface="normal"font="Arial"size="100%">CFD</style><styleface="normal"font="default"charset="134"size="100%">软件原理与应用</style></title></titles><dates><year>2004</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">清华大学出版社</style></publisher><urls></urls></record></Cite></EndNote>[4]中常把上式等号右边表示分子粘性作用的三项做如下变化，以第一式为例其中据此，有及〔65〕其中广义源项定义为当流体粘度不变且不可压缩时〔牛顿流体〕，有所以运动方程简化为〔66〕其中是运动粘度，亦是动量扩散系数，单位。本构方程和运动方程是紧密联系在一起的，通过本构方程可将应力张量用变形速度张量表示出来，即应力可用应变速率表示，而应变速率实际由速度分量决定，故使运动方程和连续性方程原那么上封闭可解。需要注意的是，这里讨论的本构方程仅局限于牛顿流体，符合广义牛顿定律的流体称为牛顿流体，否那么称为非牛顿流体。非牛顿流体的本构方程不能用广义牛顿定律描述，如对聚合物溶液等流体应该参考相关文献。4能量方程由能量守恒知，体积内流体的动能和内能的变换率等于单位时间内质量力和外表力所作的功加上单位时间内给予体积的热量。体积内流体的动能和内能的总和为：〔67〕其中是单位体积内的流体内能。质量力对体积内流体所作的功为〔单位时间内移动距离，点积求做功〕：〔68〕外表力对体积内流体所作的功为〔单位时间内移动距离，点积求做功〕：〔69〕单位时间内以热传导方式通过外表传给体积的热量为：〔70〕上式被积函数实际就是傅里叶热传导定律，即热流密度矢量正比于传热面法向温度梯度。单位时间内由于辐射或其他原因〔反响、蒸发等〕传入内的总热量为〔为单位时间内传入单位质量的热量分布函数〕：〔71〕将上述各式进行守恒计算，得〔72〕这是积分形式的能量守恒方程。求解体积分的随体导数并运用奥高公式把面积分转化为体积分可得微分形式的能量守恒方程，即因为质量守恒定律所以〔73〕另外〔74〕〔75〕那么能量方程的微分形式为：〔76〕或〔77〕或〔78〕上式各项的物理意义如下，左边第一、二项代表内能和动能的随体导数，右边第一项为哪一项单位体积内的质量力做功，第二项是单位体积内面力所作的功，第三项是单位体积内热传导输入的热量，最后一项表示由于辐射或其他物理或化学原因的热量奉献。能量守恒方程的另一种形式为〔79〕此式的物理意义为：单位体积内由于流体变形面力所作的功加上热传导及辐射等其他原因传入的热量恰好等于单位体积内的内能在单位时间内增加。将该式进一步简化，有设〔80〕为由于粘性作用机械能转化为热能的局部，称为耗散函数〔dissipationfunction〕。另外，在考虑液体流体时，比焓与内能值可看作相等，即，压力不作功。那么所以有〔81〕其中是单位体积内热源或由于辐射或其他物理或化学原因的热量奉献。一般较小可以忽略。对液体及固体可以取，进一步取为常数，并把耗散函数纳入到源项，于是〔82〕对不压缩流体〔83〕对于可以忽略粘性耗散作用的稳态低速流，能量方程可以简化为〔84〕及〔85〕取速度为零，那么可得到稳态的热传导方程〔对流项消失〕：〔86〕5状态方程由连续性方程、运动方程、能量方程确定的未知量有六个，但方程数只有五个，为使方程组封闭需补充一个联系的状态方程：〔87〕6组分质量守恒方程在一个特定的系统中，可能存在质的交换，或者存在多种化学组分，每一种组分都需要遵守组分质量守恒定律，即系统内某种化学组分对时间的变化率，等于通过系统界面的净扩散流量与由反响产生的生成率之和，可表为ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陈敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流体力学〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化学工业出版社</style></publisher><urls></urls></record></Cite><Cite><Author>帕坦卡</Author><Year>1980</Year><RecNum>10</RecNum><record><rec-number>10</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"size="100%">S</style><styleface="normal"font="default"charset="134"size="100%">.</style><styleface="normal"font="default"size="100%">V</style><styleface="normal"font="default"charset="134"size="100%">帕坦卡</style></author></authors><subsidiary-authors><author><styleface="normal"font="default"charset="134"size="100%">张政</style></author></subsidiary-authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">传热与流体流动的数值计算</style></title></titles><dates><year>1980</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">科学出版社</style></publisher><urls></urls></record></Cite></EndNote>[2,5]〔88〕其中代表单位体积内组分的质量变化率，是组分的对流流量密度。代表扩散流量密度，它由Fick定律给出。是单位体积内组分的生成率。费克定律：〔89〕其中为扩散系数。将扩散定律代入守恒方程，得〔90〕7控制方程的通用形式前面在牛顿流体的根底上，即在采用牛顿流体本构方程的根底上推导分析了运动方程和能量守恒方程，获得了较全面的流体力学方程组，同时也采用了张量不变性记法、张量分量记法及直角坐标记法三种不同方式来表示这些根本方程组，可以说各方程之间到达了初