高等流体力学-第五讲
- 格式:ppt
- 大小:2.43 MB
- 文档页数:51
概念第一章绪论连续介质:但流体力学研究的是流体的宏观运动,不以分子作为流动的基本单元,而是以流体质点为基本单元,把流场看做是由无数流体质点组成的连续体。
流体质点:流场中一个体积很小并可以忽略其几何尺寸,但与分子相比,这个体积可容纳足够多的分子数目的流体元,有一个稳定的平均特性,即满足大数定律理想流体:忽略流体黏性的流体,即μ=0.可压缩流体与不可压缩流体:简单地讲,密度为常数的流体为不可压缩流体,如水、石油及低速流动的气体。
反之,密度不为常数的流体为可压缩流体。
牛顿流体与非牛顿流体:根据流体流动时切应力与流速梯度之间的关系,即牛顿内摩擦定律。
凡是符合牛顿内摩擦定律的成为牛顿流体,如水、空气、石油等。
否则为非牛顿流体,如污泥、泥石流、生物流体、高分子溶液等动力粘度与运动粘度:动力粘度又成为动力黏度系数,动力黏度是流体固有的属性。
运动粘度又称为运动粘性系数,运动黏性系数则取决于流体的运动状态体积力与表面力:体积力亦称质量力,是一种非接触力,即外立场对流体的作用,且外立场作用于流体每一质点上,如重力、惯性力、离心力。
表面力是一种表面接触力,指流体与流体之间或流体与物体之间的相互作用,主要指压力、切应力、阻力等定常流与非定常流:又称恒定流与非恒定流。
若流场中流体质点的所有运动要素均不随时间变化,则这种流动称为定常流;反之只要有一个运动要素随时间变化则为非定常流大气层分为5层:对流层、同温层、中间层、电离层及外逸层第二章流体运动学描述流体质点的位置、速度及加速度的两种方法,即拉格朗日法和欧拉法质点导数:亦称随体导数,表示流体质点的物理量对时间的变化率,亦即跟随流体质点求导数那布拉P9流体质点的运动轨迹称为迹线流线:此曲线上任一点的切线方向就是该点流速方向依照一定次序经过流场中某一固定点的各个质点连线称为脉线,也叫序线。
流体线:在流场中任意指定的一段线,该段线在运动过程中始终保持由原来那些规定的质点所组成。
高等流体力学第一章 流体力学的基本概念连续介质:流体是由一个紧挨着一个的连续的质点所组成的,没有任何空隙的连续体,即所谓的连续介质。
流体质点:是指微小体积内所有流体分子的总和。
欧拉法质点加速度:zuu y u u x u u t u dt du a x z x y x x x x x ∂∂+∂∂+∂∂+∂∂==zu u yu u xu u tu dtdu a y zy yy xy y y ∂∂+∂∂+∂∂+∂∂==zu u y u u x u u t u dt du a z z z y z x z z z ∂∂+∂∂+∂∂+∂∂==质点的随体导数:质点携带的物理量随时间的变化率称为质点的随体导数,用dtd表示。
在欧拉法描述中的任意物理量Q 的质点随体导数表述如下:kQu t Q dt dQ k ∂∂+∂∂= 式中Q 可以是标量、矢量、张量。
质点的随体导数公式对任意物理量都成立,故将质点的随体导数的运算符号表示如下:ku t dt d k ∂∂+∂∂= 其中t ∂∂称为局部随体导数,ku k ∂∂称为对流随体导数,即在欧拉法描述的流动中,物理量的质点随体导数等于局部随体导数与对流随体导数之和。
体积分的随体导数:()dV divv dt d dV v div t dS u dV t dV dt d v v n s v v ⎥⎦⎤⎢⎣⎡Φ+Φ=⎥⎦⎤⎢⎣⎡Φ+∂Φ∂=Φ+∂Φ∂=Φ⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰()dV adivv dt da dV av div t a dS au dV t a adV dt d v v n s v v ⎥⎦⎤⎢⎣⎡+=⎥⎦⎤⎢⎣⎡+∂∂=+∂∂=⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰ 变形率张量: 11ε 12ε13εD ij = 21ε 22ε 23ε31ε 32ε 33ε其中ii ε表示所在方向的线性变形率,其余ij ε(j i ≠)为角变形率。
D ij 为变形张量。
⎪⎪⎭⎫⎝⎛∂∂+∂∂=i j j i ij x u x u 21ε 旋转角速度: 0 z ω- y ωR ij =z ω 0 x ω-y ω- x ω 0z ω=⎪⎪⎭⎫ ⎝⎛∂∂-∂∂y u x u x y 21y ω=⎪⎭⎫⎝⎛∂∂-∂∂x u z u z x 21x ω=⎪⎪⎭⎫ ⎝⎛∂∂-∂∂z u y u y z 21 判断有旋流和无旋流:x ω=y ω=z ω=0,z ω=⎪⎪⎭⎫ ⎝⎛∂∂-∂∂y u x u x y 21=0,y ω=⎪⎭⎫⎝⎛∂∂-∂∂x u z u z x 21=0 x ω=⎪⎪⎭⎫⎝⎛∂∂-∂∂z u y u y z 21=0 ,y u x u x y ∂∂=∂∂x u z u z x ∂∂=∂∂,z u y u yz ∂∂=∂∂ 涡量与速度环量的关系:涡量,流体力学中多用涡量来表示流体微团的旋转。
高等计算流体力学讲义(5)§7. TVD 格式一、背景1.求解线性波动方程0t x u au a const +==的经典差分格式 (1)一阶迎风格式(First order )11100n nj j n nj j nnj j u u a u u a u u a λ++-⎧-<⎪=-⋅⎨->⎪⎩, 其中txλ∆=∆。
上式也可以写为: 11/21/21/2111/211()()()22()()22n n n n j j j j n nn n nj j j j j n n n n n j j j j j u u f f a a f u u u u a a f u u u u λ++-+++---=--=+--=+--(2)Lax -Friedrichs 格式 (First order )()111111202n nn n n jj j j j u u u u u a t x +-++--+-+=∆∆或11/21/21/2111/211()1()()22()()22n n n nj j j j n n n nn j j j j j n n n n n j j j j j u u f f a f u u u u a f u u u u λλλ++-+++---=--=+--=+-- (3)Lax -Wendroff 格式 (second order )()()1211112222n n n n j j j jn n n j j j uu u u a t au u u txx++-+---∆+=-+∆∆∆。
或11/21/221/21121/211()()()22()()22n n n n j j j j n n nn n j j j j j n n nn n j j j j j u u f f a a fu u u u a a f u u u u λλλ++-+++---=--=+--=+--(4)Warming -Beam 格式 (Second order )()()()()12122121212212342022342022n n nn n jj j j jn n n j j j n n n n n j j j j jn n nj j j uu u u u a t au u u a txxu u u u u a t au u u a txx+----+++++-+-∆+=-+>∆∆∆-+--∆+=-+<∆∆∆2.二阶以上的差分或有限体积格式在间断附近的解可能会出现振荡。
第五讲 气动函数及压力波一、 气流参数(一)滞止参数如果按照一定的过程将气流速度滞止到零,此时气流的参数就叫做滞止参数。
滞止状态的概念可以很形象地用图5-1来表示。
它是假想把某一点处的气流引入一个容积很大的贮气箱,使其速度滞止到零。
根据一元稳定绝能流动的能量方程式2211221122h v h v +=+ 可知气体的焓值随气流速度的减小而增大。
如果把气流由速度v 1=v (焓h 1=h )绝能地滞止到v 2=0,此时所对应的焓值h 2就称为滞止焓,用符号h *表示,则*212h h v =+如果研究的是定热比容的完全气体,h =c p T ,则式(9一22)可改p c v T T /212*+= (5-1) 式中 T *称为滞止温度,它是把气流速度绝能滞止到零时的温度。
将式(5—1)两边同除以T ,则有2*2221111/1/()12212p kR k v T T v c T v T k c -=+=+=+- 所以*211Ma 2k T T -=+(5-2) 前面得到了滞止温度与温度的比与Ma 数的关系式,下面我们来推导一下其它滞止参数的表达式。
完全气体的状态方程和滞止状态的状态方程可表示为p =ρRT 和p *=ρRT * ,两者相除则有 ***()()p p ρρT T =。
(a ) 对等熵流动有p */ ρ*k =常数,p / ρk =常数,两者相比,则有**()k p p ρρ=。
(b ) 由式(a )和(b )可得**2111()(1Ma )2k kk k k p p T T ---==+ (5-3)图5-1 滞止参数模型11**2111()(1Ma )2k k k ρρT T ---==+ (5-4)由式(9-2、3、4)可知,气流参数与其滞止参数的比值只是气流Ma 数的函数。
这种函数关系是分析和计算气体流动的基础,在气体动力学中占有非常重要地位。
这里应强调的是,在气体动力学中,引进滞止状态的概念是把它作为一个参考状态。