计算流体力学2

计算流体力学的求解步骤
计算流体力学（Computational Fluid Dynamics，简称 CFD）是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。

其求解步骤通常包括以下几个方面：
1. 建立物理模型：根据实际问题建立相应的物理模型，包括流动区域、边界条件、流体性质等。

2. 数学模型：将物理模型转化为数学模型，通常使用 Navier-Stokes 方程等流体动力学基本方程来描述流体的运动和行为。

3. 网格生成：将计算区域划分为离散的网格单元，以便在每个网格点上进行数值计算。

4. 数值方法：选择合适的数值方法，如有限差分法、有限体积法或有限元法等，对数学模型进行离散化，将其转化为代数方程组。

5. 求解算法：使用适当的求解算法，如迭代法或直接解法，求解代数方程组，得到各个网格点上的流体变量的值。

6. 结果可视化：将计算得到的结果以图形或图表的形式展示出来，以便对流体的流动情况进行分析和评估。

7. 结果验证：将计算结果与实验数据或其他可靠的参考数据进行比较，验证计算结果的准确性和可靠性。

8. 优化与改进：根据结果验证的情况，对物理模型、数学模型、网格生成、数值方法或求解算法等进行优化和改进，以提高计算精度和效率。

需要注意的是，计算流体力学的求解步骤可能因具体问题和应用领域的不同而有所差异。

在实际应用中，还需要根据具体情况选择合适的软件工具和计算平台来执行上述步骤。

高等计算流体力学讲义（2）第二章 可压缩流动的数值方法§1. Euler 方程的基本理论 0 概述在计算流体力学中，传统上，针对可压缩Navier －Stokes 方程的无粘部分和粘性部分分别构造数值方法。

其中最为困难和复杂的是无粘部分的离散方法；而粘性项的离散相对简单，一般采用中心差分离散。

所以，本章主要研究无粘的Euler 方程的解法。

在推广到Navier －Stokes 方程时，只需在Euler 方程的基础上，加上粘性项的离散即可。

Euler 方程是一种典型的非线性守恒系统。

下面我们将讨论一般的非线性守恒系统以及Euler 方程的一些数学理论，作为研究数值方法的基础。

1非线性守恒系统和Euler 方程一维一阶非线性守恒系统(守恒律)可写为下列一般形式=∂∂+∂∂xF tU ，0,>∈t R x（1）其中U 称为守恒变量，是有m 个分量的列向量，即T m u u u U ),...,(21=。

T m f f f F ),...,(21=称为通量函数，是U 的充分光滑的函数，且满足归零条件，即：0)(lim=→U F U即通量是对守恒变量的输运，守恒变量为零时，通量也为零。

守恒律的物理意义设U 的初始值为：0(,0)(),U x U x x =∈R 。

如果0()U x 在x ∈R 中有紧支集（即0U 在有限区域以外恒为零），则0(,)()U x t dx U x dx =⎰⎰RR。

即此时虽然(,)U x t 的分布可以随时间变化，但其总量保持守恒。

多维守恒律可以写为)(=++∙∇+∂∂k H j G i F tU（2）守恒律的空间导数项可以写为散度形式。

守恒系统（1）可以展开成所谓拟线性形式)(=∂∂+∂∂xU U A tU （3）A 是m m ⨯矩阵，称为系数矩阵或Jacobi 矩阵，其具体形式为111122221212.........m m m m mm f f f u u u f f f u u u A f f f u u u ∂∂∂⎡⎤⎢⎥∂∂∂⎢⎥∂∂∂⎢⎥⎢⎥∂∂∂=⎢⎥⎢⎥⎢⎥∂∂∂⎢⎥∂∂∂⎢⎥⎣⎦（4），容易验证：F U Axx∂∂=∂∂，通常也记F A U∂=∂。

第二章计算流体力学的基本知识流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。

这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。

2.1计算流体力学简介2.1.1计算流体力学的发展流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。

20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。

数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学"。

从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。

数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。

数值计算方法最近发展很快，其重要性与日俱增。

自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。

最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。

航空技术的发展强烈推动了流体力学的迅速发展。

流体运动的规律由一组控制方程描述。

计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。

但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解读解。

计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力学这门交叉学科。

计算流体力学是一门用数值计算方法直接求解流动主控方程(Euler或Navier-Stokes方程)以发现各种流动现象规律的学科。

流体力学中的流体流量与流速计算流体力学是研究流体在运动过程中的性质和行为的学科。

其中，流体流量和流速是流体力学中的重要概念，用于描述流体运动的特征和量度。

本文将介绍流体流量与流速的概念及计算方法。

一、流体流量的概念及计算方法流体流量是指单位时间内通过某一截面的流体体积。

按照定义，流体流量的计算公式为：Q = A * v其中，Q表示流体流量，A表示截面面积，v表示流速。

二、流速的概念及计算方法流速是指单位时间内流体通过一个截面的体积。

流速的计算公式可以根据具体情况而定，以下是常见的几种计算方法：1. 定常流的流速计算在定常流动情况下，流体的质量流率和体积流率保持不变。

流速的计算公式为：v = Q / A其中，v表示流速，Q表示流体流量，A表示截面面积。

2. 非定常流的流速计算在非定常流动情况下，流体的流速可能随时间和空间的变化而变化。

针对不同的情况，可以采用不同的方法计算流速，如通过流速图、针对特定位置的流速计算等。

三、流体流量与流速的应用流体流量和流速是流体力学中的基本概念，广泛应用于各个领域，包括但不限于以下几个方面：1. 水泵和液压系统的设计在水泵和液压系统的设计中，流体流量和流速是重要的设计参数。

通过合理计算流体流量和流速，可以确定水泵和液压系统的工作参数，确保其正常运行。

2. 水流和气流的测量与控制在环境监测、水利工程、能源利用等领域，对水流和气流的测量与控制是常见需求。

通过准确计算流体流量和流速，可以帮助实现对水流和气流的精确测量和控制。

3. 管道流量的计算与优化对于管道流动问题，合理计算流体流量和流速有助于分析和优化管道系统的性能。

通过调整管道直径、流速等参数，可以实现管道系统的节能、减压等目标。

四、总结流体流量和流速是流体力学中的重要概念，用于描述流体运动的特征和量度。

在实际应用中，合理计算流体流量和流速，可以帮助我们设计、控制和优化各类流体系统。

因此，对于流体力学中的流体流量与流速的计算方法和应用有深入的了解，对于工程实践具有重要意义。

流体力学计算公式流体力学是研究流体的运动规律和性质的一门学科，广泛应用于工程和科学领域中。

在流体力学的研究过程中，有许多重要的计算公式和方程被提出和应用。

下面是一些重要的流体力学计算公式。

1.压力力学方程:压力力学方程是描述流体力学中流体静压力分布和变化的方程。

对于稳定的欧拉流体，方程为:∇P=-ρ∇φ其中，P是压力，ρ是流体的密度，φ是流体的势函数。

2.欧拉方程:欧拉方程用于描述流体的运动,它是流体运动的基本方程之一:∂v/∂t+v·∇v=-1/ρ∇P+g其中，v是流体的速度，P是压力，ρ是流体的密度，g是重力加速度。

3.奇异体流动方程:奇异体流动是流体与孤立涡流动的一种类型，其方程为:D(D/u)/Dt=0其中，D/Dt是对时间的全导数，u是速度向量。

4.麦克斯韦方程:5.纳维-斯托克斯方程:纳维-斯托克斯方程是描述流体的动力学行为的方程，它是流体力学中最重要的方程之一:∂v/∂t+v·∇v=-1/ρ∇P+μ∇²v其中，v是速度矢量，P是压力，ρ是密度，μ是动力黏度。

6.贝努利方程:贝努利方程描述了在不可压缩流体中流体静力学的变化。

贝努利方程给出了伯努利定律，即沿着一条流线上的速度增加，压力将降低，反之亦然。

贝努利方程的公式为:P + 1/2ρv^2 + ρgh = const.其中，P是压力，ρ是密度，v是流体速度，g是重力加速度，h是流体高度。

7.流量方程:流量方程用于描述流体在管道或通道中的流动。

Q=A·v其中，Q是流量，A是截面积，v是流速。

8.弗朗脱方程:弗朗脱方程是描述管道中流体流动的方程，其中考虑了摩擦阻力的影响:hL=f(L/D)(v^2/2g)其中，hL是管道摩擦阻力头损失，f是阻力系数，L是管道长度，D 是管道直径，v是流速，g是重力加速度。

以上是一些重要的流体力学计算公式。

这些公式和方程在流体力学中具有广泛的应用，是工程和科学领域中进行流体流动分析和计算的基础。

流体力学流速计算公式一、伯努利方程推导流速公式（理想不可压缩流体定常流动）1. 伯努利方程。

- 对于理想不可压缩流体作定常流动时，在同一条流线上有p+(1)/(2)ρ v^2+ρ gh = C（p是流体压强，ρ是流体密度，v是流速，h是高度，C是常量）。

- 假设水平流动（h_1 = h_2），则方程变为p_1+(1)/(2)ρ v_1^2=p_2+(1)/(2)ρ v_2^2。

- 由此可推导出流速公式v_2=√(v_1^2)+(2(p_1 - p_2))/(ρ)。

2. 适用条件。

- 理想流体（无粘性），实际流体在粘性较小时可近似使用。

- 不可压缩流体，像水在大多数情况下可视为不可压缩流体，气体在低速流动时也可近似为不可压缩流体。

- 定常流动，即流场中各点的流速等物理量不随时间变化。

3. 示例。

- 已知水管中某点1处的压强p_1 = 2×10^5Pa，流速v_1 = 1m/s，另一点2处的压强p_2 = 1.5×10^5Pa，水的密度ρ = 1000kg/m^3。

- 根据v_2=√(v_1^2)+(2(p_1 - p_2))/(ρ)，将数值代入可得：- v_2=√(1^2)+frac{2×(2×10^{5-1.5×10^5)}{1000}}- 先计算括号内的值：2×(2×10^5-1.5×10^5)=2×5×10^4=10^5。

- 则v_2=√(1 + 100)= √(101)≈10.05m/s。

二、连续性方程推导流速公式（不可压缩流体定常流动）1. 连续性方程。

- 对于不可压缩流体的定常流动，有S_1v_1 = S_2v_2（S_1、S_2分别是流管中两个截面的面积，v_1、v_2是相应截面处的流速）。

- 由此可推导出流速公式v_2=(S_1)/(S_2)v_1。

2. 适用条件。

- 不可压缩流体，如液体或低速流动的气体。

计算流体力学的数学模型与方法计算流体力学（Computational Fluid Dynamics，简称CFD）是研究流体运动的力学现象而采用的计算方法。

它结合了数学模型和计算方法，通过数值计算和模拟的手段，来解决流体问题。

本文将从数学模型和计算方法两个方面，探讨计算流体力学的基本原理与应用。

一、数学模型数学模型是计算流体力学的基础，它描述了流体运动的基本方程和边界条件。

常用的数学模型包括Navier-Stokes方程、动量守恒方程、质量守恒方程和能量守恒方程等。

1. Navier-Stokes方程Navier-Stokes方程是描述流体的速度和压力随时间和空间变化的方程。

其一般形式为：\[\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0\]其中，$\rho$表示流体的密度，$\mathbf{v}$表示流体的速度。

2. 动量守恒方程动量守恒方程描述了流体运动中动量的变化。

它可以表示为：\[\frac{\partial (\rho \mathbf{v})}{\partial t} + \nabla \cdot (\rho\mathbf{v} \mathbf{v}) = -\nabla p + \nabla \cdot \mathbf{\tau}\]其中，$p$表示压力，$\mathbf{\tau}$表示粘性应力张量。

3. 质量守恒方程质量守恒方程描述了流体质量的守恒。

它可以表示为：\[\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0\]4. 能量守恒方程能量守恒方程描述了流体能量的守恒。

它可以表示为：\[\frac{\partial (\rho e)}{\partial t} + \nabla \cdot (\rho e \mathbf{v}) =\nabla \cdot (\lambda \nabla T) + \nabla \cdot (\mathbf{\tau \cdot v}) + \rho \mathbf{v} \cdot \mathbf{g}\]其中，$e$表示单位质量流体的总能量，$T$表示温度，$\lambda$表示热导率。

经常用到的给排水流体力学计算公式：
1、h f=(λL/d)*(v2/2g)
h f ——流段的沿程水头损失（m液柱或气柱）
L——流段的长度（m）
d——管段的直径（m）
v——流体的流动速度（m/s）
λ——沿程阻力系数（或摩擦阻力系数），在层流运动中，该值可根据λ=64/Re求出。

给水工程经常采用钢管和铸铁管，由于管内壁容易锈蚀和积垢，所以管壁的粗糙度按旧钢管和铸铁管考虑，并为一个常数。

管内水流温度一般为10℃左右，运动粘度也可以为一个常数。

这样是的沿程阻力系数λ的经验计算公式比较简单，在紊流区内：
v＜1.2 m/s时，λ=(0.0179/d0.3)*(1+0.867/ v)0.3
v≥1.2 m/s时，λ=0.021/ d0.3
上式中，d为管道的内径（m），不是公称直径；v为流速（m/s）。

2、v=(1/n)R2/3i1/2
n——粗糙系数
R——过流断面的水利半径（m）
i——渠底或管底的坡度
常用材料的粗糙系数n值。

《流体力学》Ⅱ主要公式及方程式(总8页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--《流体力学与流体机械》(下)主要公式及方程式1．流体力学常用准数： (1) 雷诺准数 μρlu =Re (2) 欧拉准数 2Eu u p ρ= (3) 牛顿准数 22Ne l u F ρ=(4) 付鲁德准数 lg u 2Fr = (5) 马赫准数 a u=M (6) 斯特罗哈准数 l u τ=St(7) 阿基米德准数 TTu l g ∆=2Ar (8) 格拉晓夫准数23G r νβt l g ∆= (9) 韦伯准数 σρl u 2We =2．气体等压比热和等容比热计算式：1p -=k Rk C ； 1v -=k R C 3．完全气体比焓定义式：T C RT e pe i p =+=+=ρ4．完全气体状态方程式：T R p ρ= 状态方程微分式：TT p p d d d +=ρρ 5．完全气体等熵过程方程式：C p=kρ等熵过程方程微分式：ρρd d kp p = 气体压力p 、密度ρ和温度T 之间的等熵关系：1k k12k 1212)()(-==T Tp p ρρ6．气体熵增计算式：)]()ln[(ln ln 211k k121212p 12p pT T R p p R T T C s s -=-=-7．热力学第一定律的能量方程式：we u z g p q e u z g p ++++=++++22222212111122ρρ 可压缩理想流体绝热流动能量方程式： 022221122i u i u i =+=+ 以温度和流速表述： 0p 222p 211p 22T C u T C u T C =+=+ 以温度和流速表述：022221112121T R k k u T R k k u T R k k -=+-=+-以压力、密度和流速表述： 002222211112121ρρρp k k u p k k u p k k -=+-=+- 以音速和流速表述： 121212022222121-=+-=+-k a u k a u k a 8．完全气体的音速公式：T R k pk pa ===ρρd d9．理想流体一维稳定流动连续性方程式：C uA Q ==ρ 连续性方程微分式：0d d d =++AA u u ρρ10．欧拉运动方程的积分式：C u z g p=++⎰2d 2ρ 或简化为 C u p =+⎰2d 2ρ 欧拉运动方程的微分式：0d d d =++u u z g pρ或简化为0d d =+u u pρ11．理想流体稳定流动的动量方程式： ⎪⎭⎪⎬⎫-=∑-=∑-=∑)()()(z1z2z y 1y 2y x1x2x u u Q F u u Q F u u Q F ρρρ一维稳定流动动量方程微分式：0d d x=++AR u u pρδρ12．气体极限速度及临界速度计算式：120m ax -=k T kR u ； 120*+=k kRT u 13．流动参量与滞止参量间的关系：20211M k T T -+=； 1k k20)211(--+=M k p p 1k 120)211(--+=M k ρρ； 2120)211(M k a a -+= 14．无因次速度Λ与马赫数M 间的关系： 222)1()1(2ΛΛ--+=k k M15．流速的计算式： ])(1[12k1k 00---=p pRT k k u ； 或 ])(1[12k1k 000---=p pp k k u ρ无因次速度计算式：k1k 00max)(11--=-=p p T T u u16．质量流量的计算式： ])()[(12k1k 0k 2000+--=p pp p p k k AG ρ1)2(k 1k 200)211(-+--+=M k M p k A G ρ 最大质量流量计算式：00*1)2(k 1k max )12(ρp k A k G -++= 或 00*1)2(k 1k max )12(T P A k R k G -++= 17．喷管出口马赫数计算式： ]1)[(12k 1k e0e --=-p p k M 18．正激波在静止气体中传播速度计算式： 121212w ρρρρ⋅--=p p u 19．正激波后气流速度计算式： 211212)()(ρρρρ--=p p u20．正激波前后速度关系式： 2*21a u u =21．正激波前后马赫数间的关系式： )1(2)1(2212112121212122---+=---+=k M k M k k M k M k M 22．正激波前后气流参量比与波前M 1数的关系式：2121212112)1(2)1(21121M k M k M k M k -++=-++=ρρ11122112+--+=k k M k k p p]1)1(2)[112()11(2121212+---+-=M k M k k k k T T11)1(22112+-++=k k M k u u1k 1211k k21210102)1112(])1(2)1([--+--+-++=k k M k k M k M k p p23．范诺流极限管长计算式： ])1(2)1(ln 211[21212121max M k M k k k M k M D L -++++-=λ24．范诺流参量变化关系式：2*)1(21M k k T T -++=； 2122*])1(2)1([Mk M k u u -++= 2122*])1()1(2[Mk M k +-+=ρρ； 212*])1(21[1M k k M p p -++= 1)2(k 1k 2*00)1112(1-++-++=M k k k M p p 25．瑞利流参量变化关系式：2*11M k k p p ++=； 222*)11(M k k M T T ++= )11(122*kM k M ++=ρρ； )11(22*M k k M u u ++= ]1)1(2[)11(2222*00+-+++=k M k Mk k M T T 1k k 22*00]1)1(2[11-+-+++=k M k M k k p p 26．瑞利流能量方程式： 22222211u i q u i +=++ 27．等温流能量方程式： 0201i q i =+ 或 222221u q u =+ 28．等温流压降计算式：)ln2(1212112221Dl u u p u p p λρ+=- 等温流压降近似计算式：211211211M k Dlp T R u D l p p λλ-=-= 29．等温流质量流量计算式：)(16222152p p TR l D G -=λπ 30．等温流极限管长计算式： )]ln(1[212121maxM k M k M k D L +-=λ 31．等温流参量变化关系式：M k u u =∆； Mk p p 1==∆∆ρρ； T R u =∆32．等温流可能的最小压力： 11min M p k p p ==∆ 33．紊流射流主要参量计算式：35．阿基米德准数：对圆截面射流a 0200Ar T T u R g ∆=，对平面射流a200Ar T T u B g ∆=。