第2章 一维定常流动的基本方程(Part1.四个方程)
- 格式:ppt
- 大小:874.00 KB
- 文档页数:36
第二章 基本流动方程本项目研究需要计算进气道内外流耦合的复杂三维粘性流动,采用的流场控制方程为三维的Navier-Stokes 方程组,它是连续方程、动量方程、能量方程的联立方程组。
本章给出了积分形式和微分形式的N-S 方程组,以及在有限体积法离散中需使用的坐标转换后N-S 方程组,另外还给出了湍流计算使用的平均化后的湍流N-S 方程。
§2.1 积分形式N-S 方程组在直角坐标系下,忽略重力做功和辐射传热的积分形式N-S 方程组可写为如下的矢量形式:质量守恒方程⎰⎰⎰⎰⎰=⋅+∂∂sds n v d t 0ρνρν (2-1) 动量守恒方程⎰⎰⎰⎰⎰⎰⎰⎰⎰=⋅-+⋅+∂∂s ss ds n ds n p ds n v v p d t v 0ˆ)()(τνρν(2-2) 能量守恒方程⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰=⋅+⋅-⋅+⋅+∂∂s ss s ds q n ds n v ds n v p ds n v E d t E 0)ˆ()()()(τρνρν (2-3) 这是一个关于时间的双曲型方程组,其中E 代表总比能(内能和动能之和);q为能量流矢量(energy flux vector ),这里假定能量流矢量仅仅表达分子能量的运输,它可以用Fourier 定律来描述,即T grad K q ⋅-=(2-4)K 为热传导系数,根据Prandtl 数:KC p μ⋅=Pr 为常数,及Sutherland 公式:TT TK K f +⨯=0230来确定,其中0K 和 0f T 是与流体有关的常数,详见表2-1。
对于理想气体,总比能可表述为下述张量缩并的形式:()212i i i i u u pu u e E +-=+=γρ 并且总焓为:ρpE H +=这里e代表内能,H为总焓,γ为比热比。
方程(2-2)中τ是粘性应力张量(viscous stress tensor ),对于牛顿流体,假定应力张量τˆ随着变形率张量Sˆ连续变化(Stokes 假设) )(32ˆ2ˆV S⋅∇-=μμτ (2-5)变形率张量Sˆ表为 ⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂+∂∂∂∂+∂∂∂∂+∂∂∂∂∂∂+∂∂∂∂+∂∂∂∂+∂∂∂∂=z wyw z v x w z u y w z v y v x v y u x w z u xv y u x u S)(21)(21)(21)(21)(21)(21ˆ μ为粘性系数,层流的粘性系数与压力的关系可以忽略不计,于是可采用Sutherland 公式 TT T+⨯=0230μμ,其中0μ和 0T 是与流体有关的常数,详见表2-1。