流體力学基本方程组总结.docVIP

流体力学基本方程组总结 流体力学基本方程组包括连续性方程、运动方程、组分质量守恒方程、能量方程、本构方程、状态方程及通用形式守恒方程。虽各相关文献都有介绍这些基本方程组，但多采用作者熟悉的方式表示，往往同一方程具有多种形式，而难于直观对比。以下内容是对文献报道的各种形式的总结和对比，并分析了它们之间的转化关系，以期彻底理解（切实掌握微分方程中每一项的物理意义）流体力学基本方程组的数学物理意义，为离散计算该方程组打下基础。 1 连续性方程 根据文献[1]连续性方程可由四种方法得到，分别为拉格朗日法下对有限体积和体积元应用质量守恒定律、在欧拉法下对有限体积应用质量守恒定律及在直角坐标系中直接应用质量守恒定律。 L法有限体积分析 取体积为，质量为的一定流体质点团，则有： （1） 因为速度散度的物理意义是相对体积膨胀率及密度随体导数，即： （2） （3） 代入式（1）得 （4） 运用奥高定理 （5） 得 （6） 上式即是连续性方程的积分形式。 假定被积函数连续，而且体积是任意选取的，由此可知被积函数必须等于零，即： （7） 或 （8） 在直角坐标系中连续性方程为： （9） 或 （10） 连续性方程（10）表明，密度变化（随时间和位置）等于密度和体积变形的乘积[2]。 L法体积元分析 考虑质量为的体积元，对其用拉格朗日观点，根据质量守恒定律有： （11） （12） 两边同除以，得 （13） 或写成 （14） 上式表明要维持质量守恒定律，相对体积变化率必须等于负的相对密度变化率。 1.3 E法有限体积分析 着眼坐标空间，取空间中以面为界的有限体积，则称面为控制面，为控制体。取外法线方向为法线的正方向，为外法线方向的单位矢量。考虑该体积内流体质量的变化，该变化主要以下两方面原因引起。第一，通过表面有流体流出或流入，单位时间内流出流入变化的总和为： （15） 第二，由于密度场的不定常性（注意，欧拉观点下空间点是固定的，密度的变化只由场的不定常性刻画），单位时间内体积的质量将变化，变化量为： （16） 上述两者应相等，即 （17） 由于体积是任意的，且被积函数连续，则 （18） 1.4 E法直角坐标系分析 单位时间内通过表面EFGH的通量为： 通过表面ABCD的通量为： 其他三对表面类似，另外，该控制体内质量的变化率为： 则 （19） 特殊情况下的连续性方程： 定常态： 不可压缩流体： 2 动量方程 任取一体积为的流体，它的边界为。根据动量定理，体积中流体动量的变化率等于作用在该体积上的质量力和面力（应力）之和。单位面积上的面力，其中是二阶对称应力张量，所以不是通常指的在（单位体积面元的法线方向）方向的分量。单位质量上的质量力为。则作用在该体积上的质量力和面力分别为 （20） 及 （21） 动量变化率为 （22） 上述动量变化率的表达式可有两种处理方法[1]，如下 （1）求解上式右边第二项内对体积元的随体导数，则 （23） （2）对动量变化率表达式右边第二项应用质量守恒定律 （24） 由上可得两种积分形式的动量方程，即 （25） 或 （26） 由上，动量方程的微分形式为： （27） 或 （28） 微分方程中各项的物理意义为，表示单位体积上惯性力，为单位体积上的质量力，为单位体积上应力张量的散度，它是与面力等效的体力分布函数（由奥高公式转化而来）。 在直角坐标系下以应力表示的运动方程可采取下列形式 （29） 或 （30） 这两种表达方式的等号左边实际只差了一个连续性方程，由基本微分公式 （31） 得 （32） 由连续性方程知