高等流体力学湍流模型
流体力学复习要点(计算公式)
D D y S x e P gh2 gh1 h2 h1 b L y C C D D y x P hc 第一章 绪论 单位质量力： m F f B m = 密度值： 3 m kg 1000=水ρ， 3 m kg 13600=水银ρ， 3 m kg 29.1=空气ρ 牛顿内摩擦定律：剪切力： dy du μ τ=， 内摩擦力：dy du A T μ= 动力粘度： ρυ μ= 完全气体状态方程：RT P =ρ 压缩系数： dp d 1dp dV 1ρρκ= -=V （N m 2 ） 膨胀系数：T T V V V d d 1d d 1ρρα - == (1/C ?或1/K) 第二章 流体静力学+ 流体平衡微分方程： 01;01;01=??-=??-=??- z p z y p Y x p X ρρρ 液体平衡全微分方程：)(zdz ydy xdx dp ++=ρ 液体静力学基本方程：C =+ +=g p z gh p p 0ρρ或 绝对压强、相对压强与真空度：a abs P P P +=；v a abs P P P P -=-= 压强单位换算：水银柱水柱mm 73610/9800012 ===m m N at 2/101325 1m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D +=α sin 1 ) () 2(32121h h h h L e ++= ρ 若01 =h ，则压强为三角形分布，3 2L e y D == ρ 注：①图算法适合于矩形平面；②计算静水压力首先绘制压强分布图， α 且用相对压强绘制。 （2）解析法 A gh A p P c c ρ== 作用点A y I y y C xc C D + = 矩形12 3 bL I xc = 圆形 64 4 d I xc π= 曲面上的静水总压力： x c x c x A gh A p P ρ==；gV P z ρ= 总压力z x P P P += 与水平面的夹角 x z P P arct an =θ 潜体和浮体的总压力： 0=x P 排浮gV F P z ρ== 第三章 流体动力学基础 质点加速度的表达式??? ? ? ? ??? ??+??+??+??=??+??+??+??=??+??+??+??=z u u y u u x u u t u a z u u y u u x u u t u a z u u y u u x u u t u a z z z y z x z z y z y y y x y y x z x y x x x x A Q V Q Q Q Q Q G A = === ? 断面平均流速重量流量质量流量体积流量g udA m ρρ 流体的运动微分方程： t z t y t x d du z p z d du y p Y d du x p X = ??-=??-=??- ρρρ1;1;1 不可压缩流体的连续性微分方程 ： 0z u y u x u z y x =??+??+?? 恒定元流的连续性方程： dQ A A ==2211d u d u 恒定总流的连续性方程：Q A A ==2211νν 无粘性流体元流伯努利方程：g 2u g p z g 2u g p z 2 2 222 111++=++ρρ 粘性流体元流伯努利方程： w 2 2222111'h g 2u g p z g 2u g p z +++=++ρρ
湍流模型概述
大多数飞行器都是在高Re数下飞行，表面的流态是湍流。为了准确地确定湍流流态下的摩阻、热流，湍流成为一个重要而困难的研究课题。 （一）DNS 目前处理湍流数值计算问题有三种方法，第一种方法即所谓直接数值模拟方法（DNS方法），直接求解湍流运动的N-S方程，得到湍流的瞬时流场，即各种尺度的随机运动，可以获得湍流的全部信息。随着现代计算机的发展和先进的数值方法的研究，DNS方法已经成为解决湍流的一种实际的方法。但由于计算机条件的约束，目前只能限于一些低Re数的简单流动，不能用于工程应用。目前国际上正在做的湍流直接数值模拟还只限于较低的需诺数(Re～200)和非常简单的流动外形，如平板边界层、完全发展的槽道流，以及后台阶流动等。用直接数值模拟方法处理工程中的复杂流动问题，即使是当前最先进的计算机也还差三个量级。 （二）LES 另一种方法称做大涡模拟方法（LES方法）。这是一种折衷的方法，即对湍流脉动部分直接地模拟，将N-S方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。小涡对大涡的影响会出现在大涡方程中，再通过建立模型（亚格子尺度模型）来模拟小涡的影响。由于湍流的大涡结构强烈地依赖于流场的边界形状和边界条件，难以找出普遍的湍流模型来描述具有不同的边界特征的大涡结构，宜做直接模拟。相反地，小尺度涡对边界条件不存在直接依赖关系，而且一般具有各向同性性质。所以亚格子模型具有更大的普适性，比较容易构造，这是它比雷诺平均方法要优越的地方。自从1970年Deardorff第一次给出具有工程意义的LES计算以来，LES方法已经成为计算湍流的最强有力的工具之一，应用的方向也在逐步扩展，但是仍然受计算机条件等的限制，使之成为解决大量工程问题的成熟方法仍有很长的路要走。 （三）RANS 目前能够用于工程计算的方法就是模式理论。所谓湍流模式理论，就是依据湍流的理论知识、实验数据或直接数值模拟结果，对Reynolds应力做出各种假设，即假设各种经验的和半经验的本构关系，从而使湍流的平均Reynolds方程封闭。随着计算流体力学的发展，湍流模式理论也有了很大的进步，有了非常丰硕的成果。从对模式处理的出发点不同，可以将湍流模式理论分类成两大类：一类称为二阶矩封闭模式，另一类称涡粘性封闭模式。 （1）雷诺应力模式 所谓二阶矩封闭模式，是从Reynolds应力满足的方程出发，将方程右端未知的项（生成项，扩散项，耗散项等）用平均流动的物理量和湍流的特征尺度表示出来。典型的平均流动的变量是平均速度和平均温度的空间导数。这种模式理论，由于保留了Reynolds应力所满足的方程，如果模拟的好，可以较好地反映Reynolds应力随空间和时间的变化规律，因而可以较好地反映湍流运动规律。因此，二阶矩模式是一种较高级的模式，但是，由于保留了Reynolds应力的方程，加上平均运动的方程整个方程组总计15个方程，是一个庞大的方程组，应用这样一个庞大的方程组来解决实际工程问题，计算量很大，这就极大地限制了二阶矩模式在工程问题中的应用。 （2）涡粘性模式
流体力学基本公式
1流体中稳定流动和均匀流动的区别 （1）①根据当地加速度是否为0，即流体运动要素是否随时间变化，流体分为 稳定流动和不稳定流动。 ②根据迁移加速度是否为0，即流体运动要素是否随空间参数变化，流体 分为均匀流和非均匀流。（非均匀流又分为缓变流和急变流） （2）稳定流动是流场中流体质点通过空间点时所有的运动要素都不随时间改变 的流动。 （3）均匀流动是指流场中同一直线上的各流体质点的运动要素沿程不变（不随 空间参数变化）的流动。 （4）稳定流的流线可以为曲线。均匀流的流线不能为曲线，只能是一元流动。 2迹线方程最后是写成多个还是整合成一个？ 答：如果迹线方程可以合并为一个，尽量合并为一个，并且尽量消掉参数t 。如果不能合并，就不用合并。理论上说都是可以的，但是从考试的答案来说，基本上都是合并的。 流体力学基本公式 1.牛顿内摩擦定律 （1）表达式： dy du μτ±=。 （2）内摩擦定律与三个因素相关，粘性切应力与流体粘度和速度梯度有关，与 压力的大小关系不大。 （3）适用条件：牛顿流体的层流运动。 2.欧拉平衡微分方程 （1）01=??-x p X ρ，01=??-y p Y ρ，01=??-z p Z ρ （2）适用于绝对静止状态和相对静止状态，可压缩流体和不可压缩流体。 3.静力学基本方程式 （1） g p z g p z ρρ2 211+=+ （2）适用条件：重力作用下、静止的、连通的、均质流体。 （3）几何意义：静止流体中，各点的测压管水头为常数。 （4）物理意义：静止流体中，各点的总比能为常数。 4.连续性方程
（1）适用于系统的质量守恒定律在控制体上的应用。 （2）三种形式：一般形式，恒定流，不可压缩流。 ①一般形式：0)()()(=??+??+??+??z u y u x u t z y x ρρρρ ②恒定流：0)()()(=??+??+??z u y u x u z y x ρρρ ③不可压缩流体：0=??+??+??z u y u x u z y x 5.欧拉运动方程 （1） dt du z p Z dt du y p Y dt du x p X z y x =??-=??-=??-ρρρ1,1,1 （2）适用条件：所有理想流体。 6．理想流体的伯努利方程 （1）2211221222p u p u z z g g g g ρρ++=++ （2）适用条件：理想流体；不可压缩流体；质量力只有重力；沿稳定流的流线 或微小流束。 （3）几何意义：沿流线总水头为常数。 （4）物理意义：沿流线总比能为常数。 7.实际流体总流的伯努利方程 （1）221112221222w p v p v z z h g g g g ααρρ++=+++ （2）适用条件：实际流体稳定流；不可压缩流体；质量力只有重力；所取断面 为缓变流断面。 （3）动能修正系数α：总流有效断面上的实际动能与按平均流速算出的假想动 能的比值。1α>，由断面上的速度分布不均匀引起，不均匀性越大，α越大。 8.动量方程 （1）() 21=Q F v v ρ-∑
计算流体力学
1、数值的耗散与频散： 在数值解中出现的振幅衰减波长加宽的现象叫数值耗散，与高阶偶次空间偏导数有关；在数值解中出现解得主波后有一系列频及传播速度不等的尾波的现象叫数值频散，与高阶奇次偏导数有关。 2、湍流模型理论：湍流模式理论或简称湍流模型，就是以雷诺平均运动方程与脉动运动方程为基础，依靠理论与经验的结合，引进一系列模型假设，而建立起得一组 描写湍流平均量的封闭方程组。 3、修正的偏微分方程：与差分方程相等价的微分方程称之为修正的微分方程。 4、自适应网格:为了计算具有高雷诺数的流场，必须将流场内的网格加密，但是实际计算中并不需要对全流场的网格所有部分同样加密，只需在某些部分，如物面附近、尾流区等得网格加密即可。因此需要事先估计一些变化较快的区域，但这种估计又是是正确的。有时则不正确。特别是不定常流动，流动过程本身就是变化的，所以需要不断的调整网格的位置和疏密，这样就产生了自适应网格。 5、CFL 条件：定义t C x μ ?=? ，不等式1C ≤ 称为CFL 条件，此条件一般应用于双曲线偏微分方程的显式格式。物理意义：即在时间步长内，波的位移应小于空间步长。 数学意义：差分方程解的依赖区域包含微分方程解得依赖区域。 1、简答CFD 方法求解流动问题的基本步骤 答：①确定流动模型；②计算区域离散化；③用离散节点变量代替场；④将控制方程中偏导数进行离散，得到线性方程组；⑤边界条件和初值条件离散化；⑥离散的线性方程组求解，得到离散值；⑦计算结果数据处理。 2、简述离散偏微分方程的三个原则及LAX 定理 三原则；相容性、稳定性、收敛性。 LAX 定理：对于一个选定的线性偏微分方程的初值问题，对应的差分方法是相容的，则差分方程解得收敛性和稳定性事等价的或者说稳定性是收敛性的充要条件。 3、简述差分格构造的基本规律，并应用规律方程 0t x μμ λ??+=?? 利用网格点（） （）（）构造方程的差分格式，并验证其离散格式的精度等级。 答：构造的基本规律 ：①为保证均匀流场，差分的分子各项系数之和为零 ②分母向量级与微分的阶数一致 ③构造差分级指明针对哪点构造 ④差分格式的精度 由网格点（）（）（）规律方程（）构造得 1 11 1 0n n n n j j j j x x μμμμλ +++---+=?? 令112j j x k k x μμμ-+= ? 用泰勒公 式展开的23 126j j x xx x x x x μμμμμ-??=-?+- 所以12101k k k +=?? -=? 得12 11k k =-??=? 所以1j j x x μμμ--+=? 所以具有一阶精度 4、简要概括流动的数值计算对网格的基本要求 答：①计算域边界上的网格节点都应在边界上 ②物理域上的特点与计算域上的节点要求一一对应 ③网格应尽量尺寸匀称，相邻网格长度比应小于2 ④物理域网格夹角不宜太小（≥45°） ⑤流动参数梯度大的地方网格要加密，否则稀疏。 5、简述人工压缩方法（时间相关法）的基本思想 答：用非定常流动方程来求解定常流动问题，用其稳态求解定常流动的解，将不可压缩的粘性流动的连续方程，添加到可压缩项。则与动量方程构成定常粘性流动时间相关方程，可把非定常流动的稳态解作为非定常流动的解。
第1章流体力学的基本概念
第1章 流体力学的基本概念 流体力学是研究流体的运动规律及其与物体相互作用的机理的一门专门学科。本章叙述在以后章节中经常用到的一些基础知识，对于其它基础内容在本科的流体力学或水力学中已作介绍，这里不再叙述。 连续介质与流体物理量 连续介质 流体和任何物质一样，都是由分子组成的，分子与分子之间是不连续而有空隙的。例如，常温下每立方厘米水中约含有3×1022 个水分子，相邻分子间距离约为3×10－8 厘米。因而，从微观结构上说，流体是有空隙的、不连续的介质。 但是，详细研究分子的微观运动不是流体力学的任务，我们所关心的不是个别分子的微观运动，而是大量分子“集体”所显示的特性，也就是所谓的宏观特性或宏观量，这是因为分子间的孔隙与实际所研究的流体尺度相比是极其微小的。因此，可以设想把所讨论的流体分割成为无数无限小的基元个体，相当于微小的分子集团，称之为流体的“质点”。从而认为，流体就是由这样的一个紧挨着一个的连续的质点所组成的，没有任何空隙的连续体，即所谓的“连续介质”。同时认为，流体的物理力学性质，例如密度、速度、压强和能量等，具有随同位置而连续变化的特性，即视为空间坐标和时间的连续函数。因此，不再从那些永远运动的分子出发，而是在宏观上从质点出发来研究流体的运动规律，从而可以利用连续函数的分析方法。长期的实践和科学实验证明，利用连续介质假定所得出的有关流体运动规律的基本理论与客观实际是符合的。 所谓流体质点，是指微小体积内所有流体分子的总体，而该微小体积是几何尺寸很小（但远大于分子平均自由行程）但包含足够多分子的特征体积，其宏观特性就是大量分子的统计平均特性，且具有确定性。 流体物理量 根据流体连续介质模型，任一时刻流体所在空间的每一点都为相应的流体质点所占据。流体的物理量是指反映流体宏观特性的物理量，如密度、速度、压强、温度和能量等。对于流体物理量，如流体质点的密度，可以地定义为微小特征体积内大量数目分子的统计质量除以该特征体积所得的平均值，即 V M V V ??=?→?'lim ρ （1-1） 式中，M ?表示体积V ?中所含流体的质量。 按数学的定义，空间一点的流体密度为 V M V ??=→?0 lim ρ （1-2）
湍流模型发展综述
湍流模型发展综述 摘要：在概述了湍流问题的基础上，本文简要介绍了湍流的四种模型，对湍流模型在不同情况下的模拟能力进行了对比，最后简述了湍流模型的发展方向。 关键词：湍流模型；Navier-Stokes方程组；J-K模型 Abstract：On the basis of introducing the problems of turbulence, this paper briefly analyzed four kinds of turbulence models and compared their ability of simulation in different situations. At last, the paper expounded the development direction of the turbulence model. Key words：Turbulence model; Navier-Stokes equations; J-K model 一、引言 湍流又称紊流,是自然界中常见的一种很不规则的流动现象。当粘性阻尼无法消除惯性的影响时,自然界中的绝大部分流动都是湍流。 湍流运动的实验研究表明，虽然湍流结构十分复杂，但它仍然遵循连续介质的一般动力学规律，湍流流动的各物理量的瞬时值也应该服从一般的N-S方程。对粘性流体服从的N-S方程进行时均化，就可以得到雷诺平均方程。与定常的N-S方程相比，不同之处是在该式右边多了九项与脉动量有关的项，这脉动量的乘积的平均值与密度的乘积是湍流流动中的一种应力，称为湍流应力或雷诺应力。其中，法向雷诺应力和切向雷诺应力各有三个。 湍流问题就是在给定的边界条件下解雷诺方程。由于雷诺平均方程中未知数个数远多于方程个数而出现了方程不封闭的问题，这就需要依据各种半经验理论提出相应的补充方程式，即各种湍流模型。一般按照所用湍流量偏微分方程的物理含义或者数量进行区分，分别称为梅罗尔—赫林方法和雷诺方法。而后者又将湍流模型分成四类。（1）零方程模型；（2）一方程模型；（3）二方程模型；（4）应力方程模型。下面就对这些模型进行简单的描述。 二、湍流模型简介 1、零方程模型 最初的湍流模型只考虑了一阶湍流计算统计量的动力学微分方程，即平均方程，没有引进高阶统计量的微分方程，因而称之为一阶封闭模式或零方程模型。零方程模型又称为代数模型，代数模型又可以分成以下几种模型：（1）Cebeci —Smith 模型，（2）Baldwin—Lomax 模型，（3）Johnson—King 模型。 其中，B-L与C-S模型的不同之处在于外层湍流粘性系数取法不同。后者适用于湍流边界层，而前者则可用于 N-S方程的计算。此两模型已在工程计算中
流体力学公式总结
工程流体力学公式总结 第二章流体得主要物理性质 ?流体得可压缩性计算、牛顿内摩擦定律得计算、粘度得三种表示方法。１.密度ρ= m／V 2.重度γ= G /V 3．流体得密度与重度有以下得关系:γ= ρg或ρ= γ/ g ４.密度得倒数称为比体积,以υ表示υ= １/ ρ= V/m 5.流体得相对密度:d = γ流/γ水= ρ流/ρ水 6．热膨胀性 ７.压缩性、体积压缩率κ 8．体积模量 ９.流体层接触面上得内摩擦力 10．单位面积上得内摩擦力(切应力)(牛顿内摩擦定律) 1１.、动力粘度μ： 12.运动粘度ν:ν=μ/ρ 1３．恩氏粘度°E:°E = t 1 /t 2 第三章流体静力学 ?重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、压强关系换算、相对静止状态流体得压强计算、流体静压力得计算（压力体)。 1.常见得质量力: 重力ΔW = Δmg、 直线运动惯性力ΔFI ＝Δm·ａ 离心惯性力ΔＦR ＝Δm·rω２、 ２.质量力为F。:F= ｍ·aｍ= m(ｆｘi+f yj＋ｆzk) am ＝F/m = f xi＋f yj＋fｚk为单位质量力,在数值上就等于加速度 实例:重力场中得流体只受到地球引力得作用，取z轴铅垂向上,xoy为水平面，则单位质量力在x、y、z轴上得分量为 fｘ= ０,ｆy＝0 , fz=-mg/ｍ= －ｇ式中负号表示重力加速度g与坐标轴z方向相反 ３流体静压强不就是矢量,而就是标量,仅就是坐标得连续函数。即:p=p(x,y,z),由此得静压强得全微分为： 4.欧拉平衡微分方程式 单位质量流体得力平衡方程为:
流体力学计算公式
1、单位质量力：m F f B B = 2、流体的运动粘度：ρ μ=v （μ[动力]粘度，ρ密度） 3、压缩系数：dp d dp dV V ρρκ?=?-=11（κ的单位是N m 2）体积模量为压缩系数的倒数 4、体积膨胀系数：dT d dT dV V v ρρα?-=?=11（v α的单位是C K ?1,1） 5、牛顿内摩擦定律：为液体厚）为运动速度，以应力表示为y u dy du dy du A T (,μτμ== 6、静止液体某点压强：为该点到液面的距离）h gh p z z g p p ()(000ρρ+=-+= 7、静水总压力： )h (为受压面积，为受压面形心淹没深度为静水总压力，A p ghA A p p c ρ== 8、元流伯努利方程;'2221112w h g p z g u g p z ++=++ρρ('w h 为粘性流体元流单位重量流体由过流断面1-1运动至过流断面2-2的机械能损失，z 为某点的位置高度或位置水头，g p ρ为测压管高度或压强水头，g u ρ2是单位流体具有的动能，u gh g p p g u 22'=-=ρ，u gh C g p p g C u 22'=-=ρC 是修正系数，数值接近于1） 9、总流伯努利方程:w h g v g p z g v g p z +++=++222 221221111αραρ(α为修正系数通常取1） 10、文丘里流量计测管道流量：)21)(41()()(42 122211g d d d k h k g p z g p z k Q -=?=+-+=πμρρμ 11、沿程水头损失一般表达式：g v d l h f 22 λ=（l 为管长，d 为管径，v 为断面平均流速，g 为重力加速度，λ为沿程阻力系数）
计算流体力学实例
汽车外部气体流动模拟 振动和噪声控制研究所 1.模型概述 在汽车外部建立一个较大的长方体几何空间，长度约为30m，宽度和高度约为5m，在空间内部挖出汽车形状的空腔，汽车尺寸参照本田CRV为4550mm*1820mm*1685mm。由于汽车向前开进，气体从车头流向车尾，因此将汽车前方空间设为气体入口，后方空间设为气体出口，模拟气体在车外的流动。另外为了节省计算成本将整个模型按1:100的比例缩小，考虑到模型和流体均是对称的，因此仅画出几何模型的一半区域，建立对称面以考虑生成包含理想气体的流体域。在Catia中建立的模型如图1.1所示。 图1.1几何模型 2.利用ICEM CFD进行网格划分 a)导入有Catia生成的stp格式的模型； b)模型修复，删除多余的点、线、面，允许公差设为0.1； c)生成体，由于本模型仅为流体区域，因此将全部区域划分为一个体，选取方法可以 使用整体模型选取； d)为了后面的设置边界方便，因此将具有相同特性的面设为一个part,共设置了in， out，FreeWalls，Symmetry和Body； e)网格划分，设置Max element=2，共划分了1333817个单元，有225390个节点； f)网格输出，设置求解器为ANSYS CFX，输出cfx5文件。 3.利用ANSYS CFX求解 a)生成域，物质选定Air Ideal Gas，参考压强设为1atm，浮力选项为无浮力模型，
域运动选项为静止，网格变形为无；流体模型设定中的热量传输设定为Isothermal，流体温度设定为288k，湍流模型设定为Shear Stress Transport模型，壁面函数 选择Automatic。 b)入口边界设定，类型为Inlet，位置选定在in，质量与栋梁选定Normal Speed，设 定为15m/s，湍流模型设定类型为Intensity and Length Scale=0.05，Eddy Len.Scale=0.1m。 c)出口边界设定，边界类型为Outlet，位置选out。质量与动量选项为Static Pressure，相对压强为0pa。 d)壁面边界设定，边界类型为Wall，位置选在FreeWalls。壁面边界详细信息中指定 WallInfluence On Flow为Free Slip。 e)对称边界设定，边界类型为Symmetry，位置选在Symmetry。 f)汽车外壁面设定，边界类型为Wall，位置设在Body，壁面详细信息选项中指定Wall Influence On Flow为No Slip，即汽车壁面为无滑移壁面。 g)初始条件设定，初始速度分量设为U方向为15m/s，其他两个方向的速度为零。 h)求解设置，残差类型选为RMS，残差目标设定为1e-5，当求解达到此目标时，求解 自动终止。求解之前的模型如图3.1所示。 图3.1求解之前的模型 4.结果后处理 从图4.1中可以看出计算收敛。
流体力学计算公式
C3.6.2 达西摩擦因子 为了确定λ与Re 的关系,人们作了大量实验和理论研究，下面介绍有代表性的结果。 1．尼古拉兹实验 尼古拉兹（J.Nikuradse，1932）分析了达西的圆管沿程阻力实验数据后，发现壁面粗糙度对λ的影响很大，决定用人工粗糙度方法实现对粗糙度的控制。他用当地黄砂砂粒经筛选后分类均匀粘贴在管内壁上，相对粗糙度ε/d 从1/30—1/1014分6种，测得λ与Re 的关系，得到尼古拉兹图（图C3.6.1）。 2． 常用计算公式 从尼古拉兹图中看到在不同Re 数和ε/d 值的区域，λ有不同的变化规律。 图C3.6.1
（1）层流区 由泊肃叶定律推导的沿程水头损失（C3.4.10）式可得 代入达西公式（C3.6.3）式，可得层流区λ的解析式 上式表明层流区λ与管壁粗糙度无关，写成常用对数形式为 上式在双对数坐标系中是一条直线，与尼古拉兹图吻合。 （2）过渡区 该区是层流向湍流的转捩区（2000 【最新整理，下载后即可编辑】 工程流体力学公式总结 第二章 流体的主要物理性质 ? 流体的可压缩性计算、牛顿内摩擦定律的计算、粘度的三种表示方法。 1．密度 ρ = m /V 2．重度 γ = G /V 3．流体的密度和重度有以下的关系：γ = ρ g 或 ρ = γ/ g 4．密度的倒数称为比体积，以υ表示υ = 1/ ρ = V/m 5．流体的相对密度：d = γ流 /γ水 = ρ流 /ρ水 6．热膨胀性 7．压缩性. 体积压缩率κ 8．体积模量 9．流体层接触面上的内摩擦力 10．单位面积上的内摩擦力（切应力）（牛顿内摩擦定律） 11．.动力粘度μ： T V V ??=1αp V V ??-=1κV P V K ??- =κ1n A F d d υμ=dn d v μτ±=n v d /d τμ= 12．运动粘度ν ：ν = μ/ρ 13．恩氏粘度°E ：°E = t 1 / t 2 第三章 流体静力学 ? 重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、压强关系换算、相对静止状态流体的压强计算、流体静压力的计算（压力体）。 1．常见的质量力： 重力ΔW = Δmg 、 直线运动惯性力ΔFI = Δm ·a 离心惯性力ΔFR = Δm ·rω2 . 2．质量力为F 。：F = m ·am = m (f xi+f yj+f zk) am = F /m = f xi+f yj+f zk 为单位质量力，在数值上就等于加速度 实例：重力场中的流体只受到地球引力的作用，取z 轴铅垂向上，xoy 为水平面，则单位质量力在x 、y 、 z 轴上的分量为 fx = 0 , fy = 0 , fz = -mg /m = -g 式中负号表示重力加速度g 与坐标轴z 方向相反 3流体静压强不是矢量，而是标量，仅是坐标的连续函数。即： p = p (x ,y ,z )，由此得静压强的全微分为: 4．欧拉平衡微分方程式 z z p y y p x x p p d d d d ??????++=d d d d d d 0x p f x y z x y z x ??-=ρd d d d d d 0y p f x y z x y z y ??-=ρd d d d d d 0z p f x y z x y z z ??- =ρ 工程流体力学公式总结 第二章 流体的主要物理性质 流体的可压缩性计算、牛顿内摩擦定律的计算、粘度的三种表示方法。 1．密度 ρ = m /V 2．重度 γ = G /V 3．流体的密度和重度有以下的关系：γ = ρ g 或 ρ = γ/ g 4．密度的倒数称为比体积，以υ表示υ = 1/ ρ = V/m 5．流体的相对密度：d = γ流 /γ水 = ρ流 /ρ水 6．热膨胀性 7．压缩性. 体积压缩率κ 8．体积模量 9．流体层接触面上的内摩擦力 10．单位面积上的内摩擦力（切应力）（牛顿内摩擦定律） 11．.动力粘度μ： 12．运动粘度ν ：ν = μ/ρ 13．恩氏粘度°E ：°E = t 1 / t 2 第三章 流体静力学 重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、压强关系换算、相对静止状态流体的压强计算、流体静压力的计算（压力体）。 1．常见的质量力： 重力ΔW = Δmg 、 直线运动惯性力ΔFI = Δm·a 离心惯性力ΔFR = Δm·r ω2 . T V V ??=1αp V V ??-=1κV P V K ??-=κ1n A F d d υμ=dn d v μτ±=n v d /d τμ= 2．质量力为F 。：F = m ·am = m (f xi+f yj+f zk) am = F /m = f xi+f yj+f zk 为单位质量力，在数值上就等于加速度 实例：重力场中的流体只受到地球引力的作用，取z 轴铅垂向上，xoy 为水平面，则单位质量力在x 、y 、 z 轴上的分量为 fx = 0 , fy = 0 , fz = -mg /m = -g 式中负号表示重力加速度g 与坐标轴z 方向相反 3流体静压强不是矢量，而是标量，仅是坐标的连续函数。即：p = p (x ,y ,z )，由此得静压强的全微分为: 4．欧拉平衡微分方程式 单位质量流体的力平衡方程为： 5．压强差公式（欧拉平衡微分方程式综合形式） 6．质量力的势函数 7．重力场中平衡流体的质量力势函数 z z p y y p x x p p d d d d ??????++=d d d d d d 0x p f x y z x y z x ??-=ρd d d d d d 0y p f x y z x y z y ??-=ρd d d d d d 0z p f x y z x y z z ??-=ρ0 1=??-x p f x ρ10y p f y ??-=ρ01=??-z p f z ρz z p y y p x x p z f y f x f z y x d d d )d d d (??+??+??=++ρ) d d d (d z f y f x f p z y x ++=ρd (d d d )x y z p f x f y f z dU ρ=++=ρd d d d x y z U U U U x y z =f dx f dy f dz x y z gdz ??????=++++=- 《流体力学与流体机械》(上)主要公式及方程式 1．流体的体积压缩系数计算式：β1dρ p=-1dV Vdp=ρdp 流体的体积弹性系数计算式：E=-Vdpdp dV=ρdρ 流体的体积膨胀系数计算式：βdV T=1 VdT=-1dρ ρdT 2．等压条件下气体密度与温度的关系式：ρ0 t=ρ 1+βt，其中β=1 273。 3T=±μAdu dy 或τ=Tdu A=±μdy 恩氏粘度与运动粘度的转换式：ν=(0.0731E-0.0631 E)?10-4 f1?p? x-ρ?x=0?fr-1?p=0? ?ρ?r?? 4．欧拉平衡微分方程式： f? y-1?p ρ?y=0??和fθ-1?p ρ=0? f1?p?r?θ ρ?z=0?? ??f1?p? z-z-ρ?z=0?? 欧拉平衡微分方程的全微分式：dp=ρ(fxdx+fydy+fzdz) dp=ρ(frdr+fθrdθ+fzdz) 5 fxdx+fydy+fzdz=0 frdr+fθrdθ+fzdz=0 6p γ+z=C 或 p1 γ+zp21=γ+z2 或p1+ρgz1=p2+ρgz2 相对于大气时：pm+(ρ-ρa)gz=C 或pm1+(ρ-ρa)gz1=pm2+(ρ-ρa)gz2 7p=p0+γh，其中p0为自由液面上的压力。 8．水平等加速运动液体静压力分布式：p=p0-ρ(ax+gz)；等压面方程式： ax+gz=C；自由液面方程式：ax+gz=0。注意：p0为自由液面上的压力。 1 9．等角速度旋转液体静压力分布式：p=p0+γ(ω2r2 2g-z)；等压面方程式：ω2r2 2-gz=C；自由液面方程式：ω2r2 2-gz=0。注意：p0为自由液面上的压力。 10．静止液体作用在平面上的总压力计算式：P=(p0+γhc)A=pcA，其中p0为自由液面上的相对压力。压力中心计算式：yD=yc+γsinαIxc (p0+γycsinα)A Ixc ycA或yD-yc=Ixc ycA。当自由液面上的压力为大气压时：yD=yc+ 矩形截面的惯性矩Ixc计算式：Ixc= 圆形截面的惯性矩Ixc计算式：Ixc11bh3；三角形截面的惯性矩Ixc计算式：Ixc=bh3 1236π4=d 64 11．静止液体作用在曲面上的总压力的垂直分力计算式：Pz=p0Az+γVP，注意：式中p0应为自由液面上的相对压力。 12 ?ux?ux?ux?ux?+ux+uy+uz?τ?x?y?z???uy?uy?uy?uy?+ux+uy+uz直角坐标系：ay=? ?τ?x?y?z??u?uz?uz?uz?az=z+ux+uy+uz?τ?x?y?z??ax= ?ur?ur?ur?uruθ2ar=+ur+uθ+uz-?τ?rr?θ?zr ?u?u?u?uuu圆柱坐标系：aθ=θ+urθ+uθθ+uzθ+rθ ?τ?rr?θ?zr ?u?uz?uz?uzaz=z+ur+uθ+uz?τ?rr?θ?z????????? 流体质点的压力、密度等流动参量对时间的变化率计算式： dp?p?p?p?p=+ux+uy+uzdτ?τ?x?y?z dρ?ρ?ρ?ρ?ρ=+ux+uy+uz?τ?x?y?z dτ 13 drrdθdzdxdydz==== 及uxuyuzuruθuz2 ?ρ?(ρux)?(ρuy)?(ρuz)14．三维连续性方程式的一般式：+++=0 ?τ?x?y?z ?ρρur?(ρur)?(ρuθ)?(ρuz)++++=0 ?τr?rr?θ?z ?ux?uy?uz15．不可压缩流体的三维连续性方程式：++=0 ?x?y?z ur?ur?uθ?uz+++=0?rr?θ?z r 16M=ρ11A1=ρ22A2 对于不可压缩流体： Q=1A1=2A2 随体倒数 ()D u D t t ααα?= +??? ()() u u i v j w k i j k u v w x y z x y z ??????????=++?++=++ ????????? 雷诺输运定理：对系统的随体倒数求法 [()][ )]V V k V V k D dv u dv D t t D dv u dv D t t x φφφφφφ?=+????? = +?????? （ ij i j e e δ=? ()i j k i jkl l jkl il jki ijk e e e e e εεδεε??=?=== i j ijk k e e e ε?= ()()()() i j i j i j i j i i e e e e x x x x x x φ φφφ?????????=?=?=?????? ()i i i i e e x x φφφ???==?? ()i i j j i i a a e a e x x ??????=?= ????? ()()j j k i j j i j ijk k ijk i i i i j a a a a e a e e e e e x x x x εε??????=?=?==???? 1、i j u x ?? ?? ?????? ：速度梯度张量 应变率张量：表示微团的变形运动 1122112211 22ij u u v u w x y x z x v u v v w s x y y z y w u w v w x z y z z ?? ?? ???????++ ? ? ?????????? ? ? ?? ?? ????? ?=++ ? ?????? ? ???? ? ?? ??????? ? ++ ? ? ???????? ?? ? ? 旋转张量：表示旋转 3231 210 0 0ij a ωωωωωω-?? ?= ? ?-?? － 质量守恒： ()0k k u t x ρρ??+=?? 0k k u D D t x ρρ ?+=? 第二那诺雷诺输运定律： V V D D dv dv D t D t αραρ= ? ? 动量守恒定律：() u u u f t ρ ρρ?+??=??+?σ ij i i j D u f D t x σρ ρ?= +? ij i i j i j j u u u f t x x σρ ρρ???+= +??? D u f D t ρ ρ=??+σ 能量守恒定律：()1 2i i i j ij i i i i q D e u u u u f D t x x ρ σρ???? +=+- ????? 231a ω=-312 a ω=-123a ω=-ij ijk k a εω=- 第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1 计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学"。 从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解析解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解析解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力 第一章流体力学基本知识 学习本章的目的和意义：流体力学基础知识是讲授建筑给排水的专业基础知识，只有掌握了该部分知识才能更好的理解建筑给排水课程中的相关内容。 §1-1 流体的主要物理性质 1．本节教学内容和要求： 1.1本节教学内容： 流体的4个主要物理性质。 1.2教学要求： （1）掌握并理解流体的几个主要物理性质 （2）应用流体的几个物理性质解决工程实践中的一些问题。 1.3教学难点和重点： 难点：流体的粘滞性和粘滞力 重点：牛顿运动定律的理解。 2．教学内容和知识要点： 2.1 易流动性 （1）基本概念：易流动性——流体在静止时不能承受切力抵抗剪切变形的性质称易流动性。 流体也被认为是只能抵抗压力而不能抵抗拉力。 易流动性为流体区别与固体的特性 2.2密度和重度 （1）基本概念：密度——单位体积的质量，称为流体的密度即： M ρ= V M——流体的质量，kg ; V——流体的体积，m3。 常温，一个标准大气压下Ρ水=1×103kg/ m3 Ρ水银=13.6×103kg/ m3 基本概念：重度：单位体积的重量，称为流体的重度。重度也称为容重。 G γ= V G——流体的重量，N ; V——流体的体积，m3。 ∵G=mg ∴γ=ρg 常温，一个标准大气压下γ水=9.8×103kg/ m3 γ水银=133.28×103kg/ m3密度和重度随外界压强和温度的变化而变化 液体的密度随压强和温度变化很小，可视为常数，而气体的密度随温度压强变化较大。 2..3 粘滞性 （1）粘滞性的表象 基本概念：流体在运动时抵抗剪切变形的性质称为粘滞性。当某一流层对相邻流层发生位移而引起体积变形时，在流体中产生的切力就是这一性质的表 现。 为了说明粘滞性由流体在管道中的运动速度实验加以分析说明。用流速仪测出管道中某一断面的流速分布如图一所示 设某一流层的速度为u,则与其相邻的流层为u+du,du为相邻流层的速度增值，设相邻流层的厚度为dy，则du/dy叫速度梯度。 由于各流层之间的速度不同，相邻流层间有相对运动，便在接触面上产生一种相互作用的剪切力，这个力叫做流体的内摩擦力，或粘滞力。 平板实验 （2）牛顿内摩擦定律 基本概念：牛顿在平板实验的基础上于1867年在所著的《自然哲学的数学原理》中提出了流体内摩擦力的假说——牛顿内摩擦定律： 当切应力一定时，粘性越大，剪切变形的速度越小，所以粘性又可定义为流体 流体力学 一、流体的主要物性与流体静力学 1、静止状态下的流体不能承受剪应力，不能抵抗剪切变形。 2、粘性：内摩擦力的特性就是粘性，也是运动流体抵抗剪切变形的能力，是运动流体产生机械能损失的根源；主要与流体的种类和温度有关，温度上升粘性减小，与压强没关系。 3、牛顿内摩擦定律： du F A dy μ= F d u A d y τμ== 相关因素：粘性系数、面积、速度、距离；与接触面的压力没有关系。 例1：如图6-1所示，平板与固体壁面间间距为1mm,流体的动力黏滞系数为0.1Pa.S, 以50N 的力拖动，速度为1m/s,平板的面积是（ ）m 2。 解：F F A du dy δ μνμ= ==0.5 例2：如图6-2所示，已知活塞直径d=100mm,长l=100mm 气缸直径D=100.4mm,其间充满黏滞系数为0.1Pa·s 的油，活塞以2m/s 的速度运动时，需要的拉力F 为（ ）N 。 解：332 0.1[(10010)0.1]31.40.210 du F A N dy μπ--==?????=? 4、记忆个参数，常温下空气的密度3 1.205/m kg ρ=。 5、表面力作用在流体隔离体表面上，起大小和作用面积成正比，如正压力、剪切力；质量 力作用在流体隔离体内每个流体微团上，其大小与流体质量成正比，如重力、惯性力，单位质量力的单位与加速度相同，是2 /m s 。 6、流体静压强的特征： A 、垂直指向作用面，即静压强的方向与作用面的内法线方向相同； B 、任一点的静压强与作用面的方位无关，与该点为位置、流体的种类、当地重力加速度等因素有关。 7、流体静力学基本方程 0p p gh ρ=+ 2198/98at kN m kPa ==流体力学公式总结(完整资料).doc
流体力学公式总结
《流体力学》Ⅰ主要公式及方程式讲解
流体力学-公式
第二章计算流体力学的基本知识
流体力学基础学习知识知识
流体力学总结笔记