大气湍流N-S方程
前面复习
什么是湍流？ 湍流与层流有什么区别？ 雷诺数Re的表达式和物理意义？ 湍流有哪些理论？ 流体运动的稳定性指的是什么？ 处理流体运动的稳定性问题时，什么是 小扰动法和能量法？
流体力学和N-S方程
流体力学是力学的一个分支，它是研究 流体 ( 包括液体及气体 ) 这样一个连续介质 的宏观运动规律以及它与其他运动形态之 间的相互作用。
通常所说的流体力学就是 指建立在连续介质假设基础上的流体力学。
连续介质假设认为真实流体所占有的空 间可近似地看做是由“流体质点”连续无 空隙地充满着的。
所谓流体质点指的是微 观上充分大，宏观上充分小的分子团．
流体运动的描述
欧拉方法着眼于流场空间的固定点， 拉格朗日着眼于确定的流体质点。
两种方法可以互换。
K qi = qi (r , t )
qi = qi (ξ , t )
物理量的物质导数和当地导数
在欧拉方法的表达式中，专门引进了一 个运算符号d/dt，它表示某确定流体质点的 物理量随时间的变化率，称为该物理量的 物质导数；同时，将欧拉表述下物理量函 数对时间的偏导数，即空间固定点上物理 量的时间变化率，称为当地导数，记作э/эt。
dq ∂q K = + (v ⋅ ∇ ) q dt ∂t
M 1m/s M 2m/s
M’ 2m/s (t=0) M’ 3m/s (t=1s)
应力张量
流体质点所受的力需要用二阶张量来描 述，σji。
在通过某点并具有任意方向n的面 元上，应力矢量 T(n) 为二阶张量和该面元 的法向单位矢n唯一确定。
K Ti (n ) = σ ji n j
流体力学方程组
连续方程、动量方程、能量方程
dVi ∂Vi ∂Vi 1 ∂σ ij = +Vj = + fi dt ∂t ∂x j ρ ∂x j
∂ρ ∂ + ( ρvi ) = 0 ∂t ∂xi
σ ij = − pδ ij + λekk δ ij + 2 μeij
1 ∂vi ∂v j eij = ( + ) 2 ∂v j ∂vi
∂u i 1 ∂qi dε 1 = σ ij − ∂x j ρ ∂xi dt ρ
流体的分类
粘性流体Re 无粘流体Re→∞ 不可压缩流体 可压缩流体
流动的分类
有旋运动 无旋运动 定常运动 不定常运动 一维运动 二维运动和三维运动
利用不可压缩流的连续方程对动量方程 进行简化即可得到N-S方程：
粘性流体的不可压缩流动 N－S方程
K K K dv ∂v K 1 K K = + ( v ⋅ ∇ ) v = − ∇ p + υ∇ 2 v + f dt ∂t ρ
。
第二章大气运动的基本特征
第二章 大气运动的基本特征地球大气的各种天气现象和天气变化都与大气运动有关。
大气运动在时间和空间上具有很宽的尺度谱，天气学所研究的是那些与天气和气候有关的大气运动。
对这些运动，可忽略离散的分子特性，可以视大气为连续的流体介质，表征大气状态的物理变量（如气压、密度、温度）在大气这具有单一的值，这些场变量和它们的导数是空间和时间的连续函数，控制大气运动的流体力学和热力学基本定律可以用场变量作为因变量和空间、时间变量作为自变量的偏微分方程表示。
大气运动受质量守恒、动量守恒和能量守恒等基本物理定律所支配。
§2—1 影响大气运动的作用力一、 基本作用力影响大气运动的基本作用力：是指大气与地球或大气之间的相互作用而产生的真实力，它的存在与参考系无关。
1、 气压梯度力：作用于单位质量的气块上的净压力，称为气压梯度力。
当气压分布不均匀时，气块会受到净压力的作用。
P G ∇-=ρ1（1）其中，ρ为气块密度， k zp j y p i x p P∂∂+∂∂+∂∂=∇称为气压梯度力。
P ∇是由于气压分布不均匀而造成的。
气压梯度力与气压梯度成正比，与密度成反比。
方向指向P ∇-的方向，即由高压指向低压的方向。
2、 地心引力由牛顿万有引力定理说明，宇宙间任何两个物体之间都具有引力：⎪⎭⎫ ⎝⎛-=r r r GMm F g2所以，地球对单位质量空气的引力（地心引力）为：⎪⎭⎫ ⎝⎛-=r r r GM m F g2 设：地球的半径为a（地心到海平面的距离），海拔高度为z，则()()⎪⎭⎫⎝⎛+-=⎪⎭⎫ ⎝⎛+-=r r az a Gm r r z a GM g222*112*01⎪⎭⎫⎝⎛+=a z g 在气象学范围内，z 的值一般为数十公里，而地球半径a 竟达6000多公里，故**g g ≈可作为常数。
地心引力始终是作用于大气的真实的力。
3、 摩擦力大气是一种粘性流体，它同任何实际流体一样都受内摩擦的影响。
可压缩湍流化学反应流的数值模拟研究
可压缩湍流化学反应流的数值模拟研究可压缩湍流化学反应流的数值模拟研究在化学工程领域具有重要意义，能够帮助工程师和研究人员更好地理解和优化工业过程中的化学反应流动。
本文将介绍湍流化学反应流的数值模拟研究的基本原理、方法和应用。
在可压缩湍流化学反应流的数值模拟研究中，首先需要建立数学模型来描述流体运动和化学反应。
通常采用的是基于Navier-Stokes方程和化学动力学方程的耦合模型。
Navier-Stokes方程描述了流体的运动，包括惯性力、粘性力和压力力等因素；化学动力学方程则描述了化学物质在流动中的传输和反应。
通过求解这些方程组，可以得到流体流动和化学反应的详细信息。
数值模拟研究中，对于湍流流动的模拟通常采用雷诺平均N-S方程（RANS）和大涡模拟（LES）等方法。
RANS方法适用于大规模湍流结构较为简单的情况，而LES方法则能够模拟湍流中的小尺度涡旋，更适用于复杂流动的模拟。
在化学反应方面，通常采用有限元法或有限体积法求解化学动力学方程，考虑各种化学反应过程的影响。
可压缩湍流化学反应流的数值模拟研究在石油化工、环境保护、新能源开发等领域具有广泛的应用。
例如，在石油炼制过程中，通过模拟反应流动的变化，可以提高反应器的效率和产率；在环境保护方面，可以模拟大气中的化学反应过程，预测有害物质的扩散和传输；在新能源开发中，可以模拟燃烧过程，优化能源利用效率。
总之，可压缩湍流化学反应流的数值模拟研究是一项复杂而重要的工作，需要综合考虑流体力学、化学动力学和数值方法等多方面知识。
通过不断改进模型和算法，将能够更好地应用于工程实践，推动工业化学反应流动的进步和发展。
大气湍流基础汇总
• （7）流动特性。湍流不是流体的物理属性，而是流动的运 动性质，所以不同的流体其湍流特征往往也不一样，例如 边界层湍流与尾迹湍流，正因为如此（湍流依赖于外部条 件，如边界条件），所以工程上很难对湍流进行统一的模 式处理，但是湍流的一些本质特征是普适的，寻找这些普 遍规律正是湍流理论研究的中心任务；
• 适用条件: 湍强不太大 均匀湍流 平稳湍流
大气边界层湍流风速时间序列(一)
大气边界层湍流风速时间序列(二)
湍流是随机的，复杂的 但可以通过统计的方法 来研究
• 大气边界层内的湍流总是包括很多大小不 同、相互叠加的湍涡，这些不同尺度湍涡 的相对强度定义为湍流谱。
• 最大的边界层湍涡接近边界层的厚度（100 －3000km），最小湍涡尺度只有几毫米， 由于分子粘性的耗散作用，其Βιβλιοθήκη 度非常微 弱。小湍涡以大湍涡为能源。
大气湍流基础
王成刚 大气物理系
平流层 对流层 边界层
~ 10 km 1~2 km
什么是湍流？——湍流现象
达芬奇描绘的湍流
火山爆发
杂乱、随机、无序
什么是湍流——湍流的定义
• Von.Karman和I.G Taylor对湍流的定义：湍流是流体和气 体中出现的一种无规则流动现象，当流体流过固体边界或 相固流体相互流过时会产生湍流。
湍流的主要特征（三）
• （8）记忆特性（相关性）。湍流运动在不同的时刻或空间不同 点上并不是独立的，而是有相互关联，但这种关联随着时间间隔 或空间距离的增大而变小，最后趋近于零；
• （9）间歇性。内间歇：充分发展的湍流场中某些物理量（特别 是高阶统计量）并不是在空间（或时间）的没一点上都存在的， 即有奇异性。外间歇：指湍流区与非湍流区边界的时空不确定性， 例如积云与蓝天之间的界面。间歇现象是近代湍流研究的重大发 现之一，目前是湍流理论研究的前沿课题;
大涡模拟滤波网格尺度研究及其应用
大涡模拟滤波网格尺度研究及其应用一、本文概述本文旨在深入探讨大涡模拟（Large Eddy Simulation, LES）中的滤波网格尺度问题，以及其在流体动力学领域的应用。
大涡模拟作为一种重要的湍流模拟方法，能够捕捉到湍流中的大尺度结构，并通过模型描述小尺度运动对大尺度的影响。
滤波网格尺度作为大涡模拟中的关键参数，其选择直接影响到模拟的精度和效率。
因此，研究滤波网格尺度对于提高大涡模拟的准确性和适用性具有重要意义。
本文首先将对大涡模拟的基本理论和方法进行概述，介绍滤波网格尺度在大涡模拟中的作用和影响。
然后，通过对不同滤波网格尺度下的模拟结果进行比较分析，探讨滤波网格尺度对模拟精度和计算效率的影响机制。
在此基础上，本文将提出一种优化的滤波网格尺度选择方法，以提高大涡模拟的准确性和效率。
本文还将探讨大涡模拟在流体动力学领域的应用，特别是在复杂流动和工程实际问题中的应用。
通过具体案例的分析和讨论，展示大涡模拟在解决实际问题中的潜力和优势。
本文将全面系统地研究大涡模拟中的滤波网格尺度问题及其应用，为大涡模拟在流体动力学领域的应用提供理论支持和实践指导。
二、大涡模拟理论基础大涡模拟（Large Eddy Simulation，简称LES）是一种介于直接数值模拟（DNS）和雷诺平均N-S方程（RANS）之间的湍流数值模拟方法。
它的主要思想是将湍流运动通过某种滤波函数分解为大尺度运动和小尺度运动两部分，大尺度运动通过直接求解滤波后的N-S方程得到，而小尺度运动对大尺度运动的影响则通过模型来模拟。
在LES中，滤波函数的选择至关重要。
常用的滤波函数包括盒式滤波、高斯滤波等。
滤波后的N-S方程会包含一个新的未知量，即亚格子应力张量。
为了封闭这个方程，需要引入亚格子尺度模型（Subgrid-Scale Model，简称SGS模型）。
SGS模型的作用是模拟小尺度湍流对大尺度湍流的影响，从而使方程封闭可解。
在大涡模拟中，网格尺度是一个关键参数。
第九章大气湍流扩散
∂Xv (xv,t)
=
v V(
xv,
t)
=
uv(
v X
,
t
)
=
uv(
v X
( xv,
t ), t )
∂t
∫ v
X
( xv, t )
=
xv
+
t
v V
(
xv
,
t
′)dt
′
t0
粒子在τ时间内的位移
令Y(τ)为粒子在时间间隔（t0，t0+τ）内
的位移矢量，τ为扩散时间，即
∫ Yv(τ ) = Xv (xv, t0 + τ ) − Xv (xv, t0 ) =
t
′)dt
′
t0
Y2′(τ ) x 1
t0 +τ
横向位移脉动均方差
不难求出位移脉动均方值为
t0 +τ
∫ ∫ σ
2 Y2′
(τ
)
=
V2′(x, t1 )V2′(x, t2 )dt1dt2
t0
Lagrange横向相关系数RL
定义Lagrange横向相关系数RL：
RL
(t1
−
t2
)
=
V2′(
x,
t1 )V2′(x, V2′2
复习
动力分子粘滞系数
μ=1.789×10-5kg·m-1·s-1，
运动分子粘滞系数
ν=1.461×10-5m2·s-1。
分子热传导率
κθ=2.53×10-2W·m-1·K-1， 分子热扩散率νθ=2.06×10-5m2·s-1。 定压干空气比热Cpd=1004.67m2·s-2·K-1。
瞬时n-s方程
瞬时n-s方程
瞬时n-s方程描述了流体力学中的速度和压力之间的关系，它是求解流体运动问题的重要工具。
这个方程是基于质量守恒和牛顿第二定律的基础上推导出来的，它可以描述流体在不同位置和时间的速度和压力分布。
在瞬时n-s方程中，速度是描述流体运动状态的重要参数。
它可以用来描述流体的流动性质，比如流速、流量等。
速度的大小和方向可以随着时间和位置的变化而变化。
而压力则是描述流体内部分子之间相互作用力的参数。
压力的大小和方向可以随着时间和位置的变化而变化。
在瞬时n-s方程中，速度和压力是通过一系列偏微分方程来描述的。
瞬时n-s方程的求解可以通过数值模拟和实验方法来进行。
数值模拟方法是通过将流体运动问题转化为一系列离散的时间和空间网格上的方程，然后利用数值方法求解这些方程得到速度和压力的分布。
实验方法则是通过实验测量来获取速度和压力的分布。
瞬时n-s方程的应用非常广泛，可以用来研究各种流体运动问题。
比如在航空航天领域，可以用瞬时n-s方程来研究飞机的空气动力学性能；在能源领域，可以用瞬时n-s方程来研究燃烧和传热过程；在环境工程领域，可以用瞬时n-s方程来研究水污染和水资源管理等问题。
瞬时n-s方程是流体力学中非常重要的方程，它可以用来描述流体的运动状态，求解流体运动问题。
通过研究和应用瞬时n-s方程，可以深入理解流体运动的规律，为解决各种工程和科学问题提供重要的参考依据。
空气动力学部分知识要点
空气动力学及飞行原理课程空气动力学部分知识要点一、流体属性与静动力学基础1、流体与固体在力学特性上最本质的区别在于：二者承受剪应力和产生剪切变形能力上的不同。
2、静止流体在剪应力作用下（不论所加剪切应力T多么小，只要不等于零）将产生持续不断的变形运动（流动），换句话说，静止流体不能承受剪切应力，将这种特性称为流体的易流性。
3、流体受压时其体积发生改变的性质称为流体的压缩性，而抵抗压缩变形的能力和特性称为弹性。
4、当马赫数小于0.3 时，气体的压缩性影响可以忽略不计。
5、流层间阻碍流体相对错动（变形）趋势的能力称为流体的粘性，相对错动流层间的一对摩擦力即粘性剪切力。
6、流体的剪切变形是指流体质点之间出现相对运动（例如流体层间的相对运动）流体的粘性是指流体抵抗剪切变形或质点之间的相对运动的能力。
流体的粘性力是抵抗流体质点之间相对运动（例如流体层间的相对运动）的剪应力或摩擦力。
在静止状态下流体不能承受剪力；但是在运动状态下，流体可以承受剪力，剪切力大小与流体变形速度梯度有关，而且与流体种类有7、按照作用力的性质和作用方式，可分为彻体力和表面力（面力）两类。
例如重力，惯性力和磁流体具有的电磁力等都属于彻体力，彻体力也称为体积力或质量力。
8、表面力：相邻流体或物体作用于所研究流体团块外表面，大小与流体团块表面积成正比的接触力。
由于按面积分布，故用接触应力表示，并可将其分解为法向应力和切向应力：9、理想和静止流体中的法向应力称为压强，其指向沿着表面的内法线方向，压强的量纲是［力］/［长度］210、标准大气规定在海平面上，大气温度为15 C 或T o =288.15K，压强p o = 760毫米汞柱二101325牛/米2，密度p二1.225 千克/米311 、从基准面到11 km 的高空称为对流层，在对流层内大气密度和温度随高度有明显变化，温度随高度增加而下降，高度每增加1km，温度下降6.5 K。
从11 km到21km的高空大气温度基本不变，称为同温层或平流层，在同温层内温度保持为216.5 K。
高等大气物理：第六章 湍流的发生与描述
8
湍流运动的发生机制和运动结构
？
湍流运动的形成首先与流体运动的 稳定性密切相关
9
9
二、流体运动的稳定性
1、Re数：1883年， O.Reynolds在圆管水流试验中系统地 研究了湍流现象，引进了Re数
Re
> <
Rec
湍流 层流
不可压缩流体的稳定性可以用Re这样的无因次特征量 来表达
10
10
2、线性小扰动方法
• 奇怪吸引子的出现是由于高维相空间中的耗散系统，在演 化过程中要耗损掉快弛豫参量，剩下决定系统长时间行为 的慢弛豫参量； • 在这过程中，系统的相体积要不断地收缩,并趋向一个维 数比原来相空间维数低的有限区域——吸引子上；方程的 非线性，使得某些方向上的运动是不稳定的，局部看来呈 指数分离。为了在有限的区域里进行指数分离，空间运动 轨道只能采取无穷次折迭起来的办法
讨论基本流场的稳定条件 平均场＋扰动场 求解流体扰动衰减条件
11
不可压缩流体的navier-stokes方程
1 ∂p ∂ui ∂ui + uj =− + ν∆ui ∂x j ρ ∂x i ∂t ∂u j ∂x j =0
（3）
（2）
描述液体力学和等离子体的基本方程是纳维尔-期斯托克 斯方程，在2000年5月14日被巴黎国际数学大会列为21世纪世 界数学的七大难题之一
第六章 湍流的发生与描述
提出问题
湍流是怎样产生的？
利用实验方法来研究：动力扰动、热力扰动 利用数学方法来研究：Lorenz非周期
寻找解决 方法和思 路
得出结论
湍流形成过程流体运动的特征 混沌运动的特征
湍流的混沌特征
2
2
湍流运动方程
湍流运动方程
湍流运动方程是描述湍流现象的基本方程之一，它在流体力学中具有重要的意义。
湍流是指流体在运动过程中产生的无规则、混乱的流动状态。
它具有高度的不可预测性和复杂性，因此对湍流的研究一直是流体力学领域的重点之一。
湍流运动方程能够描述湍流的运动规律和特性，为湍流现象的研究和应用提供了理论基础。
湍流运动方程的基本形式可以用一维不可压缩Navier-Stokes方程组来表示，即连续性方程和动量方程。
连续性方程描述了流体质点的质量守恒，它表明在流体运动过程中，单位体积内的质量保持不变。
动量方程则描述了流体质点的动量守恒，它包含了流体的惯性力、压力力和粘性力等因素的影响。
湍流运动方程的求解是一个极为复杂的问题，因为湍流本身的非线性和不可预测性使得方程的解析解难以得到。
目前，研究者主要采用数值模拟方法来求解湍流运动方程，通过计算机模拟湍流现象的发展和演化过程，从而获得湍流的统计特性和动力学行为。
湍流运动方程的研究不仅对于科学研究具有重要意义，也在工程技术领域有着广泛的应用。
例如，湍流运动方程的研究可以用于优化飞机、汽车和船舶等交通工具的设计，改善其流体力学性能；同时，湍流运动方程的研究也可以应用于气象预报、水力学和环境工程等领域，提高相关问题的解决效率和准确性。
湍流运动方程是描述湍流现象的基本方程，它具有重要的理论和应用价值。
通过研究湍流运动方程，可以深入了解湍流的运动规律和特性，为湍流现象的研究和应用提供理论支持。
湍流运动方程的研究不仅对于科学研究有着重要意义，也在工程技术领域具有广泛的应用前景。
通过进一步深入研究和探索，相信湍流运动方程的应用将会得到进一步的拓展和发展。
第二章 光在湍流大气中传输的理论概述
2.1 大气折射率在光学频率范围内，对流层（高度<17km）中的地球大气的空气折射率表示如下：n=1+77.6（1+7.52×10-3λ-2）（p/T）×10-6 （2.1）式中，p是以mbar为单位的大气气压，T是热力学温度，λ是以μm为单位的光波波长，由于地面上温度对n1（r）的贡献＜1%，故（2.1）式中忽略了与水汽压相关的项，当然这一项对水上传播光路是不可忽略的。
2. 2 大气湍流描述自然界中的流体运动存在着二种不同的形式：一种是层流，看上去平顺、清晰，没有掺混现象；另一种是湍流，看上去毫无规则，显得杂乱无章。
例如，如果流体以一定的速度流过一个管子，我们可以用带颜色的染料对它进行观察，在流体速度低的时候，流线光滑面清晰，流体处于层流状态；不断增加流体速度，当流速达到一定值时，流线就不再是光滑的了，整个流体开始作不规则的随机运动，流体处于湍流状态。
自从1883 年Reynolds 做了著名的湍流实验以来，以Monin-Obukhov 提出的相似理论、Deardorff 提出的大涡模拟、美国Kansas 州观测实验等为代表，大气湍流的研究已经取得了很大的进展和丰硕的成果，并在天气、气候研究和工程实际中获得成功地应用。
湍流对大气中声、光和其它电磁波的传播具有极为重要的影响，例如湍流风速、温度和湿度的脉动都会引起声音散射和减弱，大气小尺度光折射率的起伏(称为光学湍流)，会严重影响光的传播和光学成像的质量等等。
长期以来，以Tatarskii 的工作为代表，声光电传播的湍流效应大都是按照Kolmogorov 的均匀、平稳和各向同性假设处理的，而实际的湍流经常不满足这些假设，要建立更加完善的波动传播模型就必须考虑湍流的各向异性、以及间歇性的影响。
2. 3 折射率湍流模型在湍流大气中，折射率在不同地点、不同时刻都是变化的。
一方面，我们还不可能对这些变化作出预测；另一方面，即使已知这些变化，要对所有时刻、所有地点的值作出描述实际上也是不可能的。