大气湍流N-S方程

第二章 大气运动的基本特征地球大气的各种天气现象和天气变化都与大气运动有关。

大气运动在时间和空间上具有很宽的尺度谱，天气学所研究的是那些与天气和气候有关的大气运动。

对这些运动，可忽略离散的分子特性，可以视大气为连续的流体介质，表征大气状态的物理变量（如气压、密度、温度）在大气这具有单一的值，这些场变量和它们的导数是空间和时间的连续函数，控制大气运动的流体力学和热力学基本定律可以用场变量作为因变量和空间、时间变量作为自变量的偏微分方程表示。

大气运动受质量守恒、动量守恒和能量守恒等基本物理定律所支配。

§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方法则能够模拟湍流中的小尺度涡旋，更适用于复杂流动的模拟。

在化学反应方面，通常采用有限元法或有限体积法求解化学动力学方程，考虑各种化学反应过程的影响。

可压缩湍流化学反应流的数值模拟研究在石油化工、环境保护、新能源开发等领域具有广泛的应用。

例如，在石油炼制过程中，通过模拟反应流动的变化，可以提高反应器的效率和产率；在环境保护方面，可以模拟大气中的化学反应过程，预测有害物质的扩散和传输；在新能源开发中，可以模拟燃烧过程，优化能源利用效率。

总之，可压缩湍流化学反应流的数值模拟研究是一项复杂而重要的工作，需要综合考虑流体力学、化学动力学和数值方法等多方面知识。

通过不断改进模型和算法，将能够更好地应用于工程实践，推动工业化学反应流动的进步和发展。

• （6）连续性。湍流是一种连续介质的运动现象，因此满足 连续介质力学的基本规律，例如N-S方程；
• （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方程，可以深入理解流体运动的规律，为解决各种工程和科学问题提供重要的参考依据。

空气动力学及飞行原理课程空气动力学部分知识要点一、流体属性与静动力学基础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
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 折射率湍流模型在湍流大气中，折射率在不同地点、不同时刻都是变化的。

一方面，我们还不可能对这些变化作出预测；另一方面，即使已知这些变化，要对所有时刻、所有地点的值作出描述实际上也是不可能的。

基于混沌理论对湍流现象的一些探讨江苏大学京江学院 张炳乾摘要：本文简要介绍混沌理论，并在前辈研究的基础上，用混沌理论对湍流现象展开思考。

正反两极的系统造就整体稳定局部不稳定的运动状态，在内部不断变化的机制下，用矛盾对立统一的辩证思维，从事物演化方向进行探讨。

关键词：混沌理论；湍流；随机性；奇怪吸引子；分形中图分类号：O415.5 文献标识码：A 文章编号：2096-4595（2020）30-0047-0002湍流现象在日常生活中十分常见，例如海洋里的波涛汹涌、烟囱上的袅袅炊烟、大气中的乱云飞渡等，但是湍流的问题至今还未能解决。

长期以来，湍流都被认为是经典物理留下的世纪难题[1]，令许多科学家们望而生畏。

湍流最基本的特征就是流体微团运动的随机性，正是这种随机运动引起的变化使得湍流异常复杂。

而混沌理论与这种随机性极为相似，为此，从混沌理论角度上讨论湍流现象，希望从中得到科学的启发。

一、混沌理论“混沌”一词最早出现在中国和希腊的神话故事里，指代宇宙最开始的自然状态，是古代的朴素宇宙观。

随着人类文明的发展，“混沌”一词在文学、科学和宗教等著作中广泛使用。

混沌理论是系统演化对初始条件的异常敏感，差异的发展导致未来路径的演化发生变化。

美国物理学家约克（James Yorke）提出：“混沌是宇宙的自然状态，在混乱中，复杂系统可以被不断完善，秩序形成于混沌之中。

混沌不是混乱的、随机分散的，相反，其中的模式是非常有序的，只是较为复杂，混沌指的就是这种复杂的秩序化。

”[2]我们常常认为混沌是一种虚假空幻的状态，世界本没有混沌，而是人们对客观世界的思考意识受到世俗文化影响被潜移默化。

但实际上，混沌的确存在，不过不是我们常规思维所意识到的直接观念。

现代混沌理论强调有序与无序、确定性与不确定性、周期性与非周期性，看似矛盾，实则对立统一。

混沌是一种联系的状态，以整体的形式存在，是一个不能分割、相互联系、不断运动的实体[3]。

2、风速随高度的变化
地表附近的气流运动受下垫面（地表面的树林、田野、建筑物 等）摩擦力的阻滞，使气流在垂直方向上产生的风速梯度。
表示风速随高度变化的曲线，称为“风速廓线”。
风速随高度的增加而增加。 风速大小直接决定了该地区的大气污染稀释能力（环境容量）。
对数律
u*
z
u ln
k
z0
u——高度为z处的风速，m/s ； u*——摩擦速度，m/s ； k—— Karman常数，一般取0.4； z0 ——地面粗糙度，m. z0定义为离地面高度z0处，平均风速u＝0.
注意：=>
m
u2

u 1

Z2 Z1

Z 200 m
Z 2 200 m ，取 Z 2 200 m
稳定度等级 城市 乡村 城乡平均
不同稳定度下的风廓线指数m的取值
A 0.10 0.07 0.085
B 0.15 0.07 0.11
C 0.20 0.10 0.15
D 0.25 0.15 0.20
对数律较适用于描述近地层中性层结时的风廓线。
指数律（幂律）
m
z u u 10
10
u——高度为z处的平均风速，m/s ； u10——10m处的平均风速，m/s ； m——风廓线指数，也称风速高度指数。 与大气稳定度和地面粗糙度有关。
气象台站的观测风速是在距地面10m高度10min内的平均风速。
由各种热力因子诱生热力作用，也称为对流。如太阳加热地表 导致热对流泡向上运动，地表受热不均匀，气层不稳定等。
3、尺度
湍流是由大大小小的湍涡组成。 ► 最 大 的 湍 涡 尺 度 约 和 边 界 层 厚 度 相 当 ， 1000 ～ 10000km ， 生 命 周期为半天到几天； ►最小的湍涡仅几个mm，生命周期仅百分之几秒。

大气模型目标距离：t Z （m ） 目标天顶角：θ（°） 横向目标速度：(/sec)t U m 风速：(/sec)U m ω 湍流结构常数：223()nC m -吸收系数：1()m α- 散射系数：1()m σ- 激光出射半径：()m R m 平均功率：((/sec))P W J 波长：()m λ 光束质量：M 抖动：2J σ2()rad漂移：22()rad ωσ激光到达目标表面时的光斑半径为：2222222222()()L B D T J B R Z Z ωσσσσσσσ=+=++++吸收数：A t N Z α= 旋转数：/s t w N U U =菲涅尔数：2(2)/()F m t N R Z πλ=热畸变数：93(2.33310)[()/()]D t A m w N PZ N R U -≅⨯线性扩展：0220.5[/()]Dm R σλπ=， 2'20.5[()/(2)]p m mM R σλ=湍流均方束散半角22025223000.182()()()[() 1.18()e D T e e D D M r D D Mr r σσσ⎛⎫⎪⎪= ⎪- ⎪⎝⎭ 有效光束质量：0122222''()DTJQ D m N m σσσσ++=， 012222'()/D TJ QD N Mσσσσ++=1大气衰减模型假设地面有一个功率为P 的连续高能激光发射系统，目标距离为Z ，则激光平均功率密度：20.63exp()/I P Z R μπ=-exp()Z μ-为大气透过率若激光射出几条谱线，各谱线的功率权重因子为i ε ，对应的功率为i P则i i P P ε=。

11nii ε==∑1exp[]exp[]ni i i P h P h αεα--∆=-∆∑平均衰减系数为11ln[exp[]]ni i i h h αεα-=--∆∆∑斜程光学厚度1sec ns i j j h ταθ==∆∑斜程透过率1exp[()]exp[sec ]exp[]Lns j j s j T Z dZ h ααθτ==-=-∆=-∑⎰2热晕模型：热晕畸变常数0)T ZtDP w mn R P N n C Cl R αρ-'=,T nT d N d =为压强一定时大气折射率随温度变化率,m R 为出射半径，P 为激光功率，Cp 为大气定压比热容，ρ、0n 分别为未受扰动大气密度和折射率0D D F P m nN N N '==3湍流模型 3-1 B-L 模型湍流粘性系数为0min[(),()]t t i t μμμ=，（i 、0分别指内外层的值）其中，2()Re t i L μρ=Ω，且Re Re/)∞=，Ω为涡矢量大小。

海上无测风区域风速评估摘要：海上风电场风能资源评估通常需要采集风电场位置完整年数据进行分析计算。

本文提出基于风资源图谱成果的风速分析法，基于3Tier中尺度数据集的风速分析法和基于WT的风速分析法对海上无测风区域110m高度的风速大小进行评估。

三种评估方法计算得到风电场场址中心点轮毂高度出的风速分别为7.82m/s、7.87m/s和7.84m/s。

三种评估方法的评估结果基本一致，可为无实测资料地区海上风电场风资源评估分析提供参考。

关键词：风电场风资源；风速评估0 导语海上风电场风能资源评估通常需要采集风电场位置完整年数据进行分析计算。

目前，海上风电场风能资源观测和获取有两种主要途径[1]：一是海上测风塔，通过在特定高度安装机械式测风仪的方式进行观测；二是海上漂浮式激光雷达测风[2]，通过搭载在浮体上的激光雷达进行不同高度风资源观测。

近年来，由于海上风电技术的飞速发展，海上风电开发逐步从近海走向深远海。

随着海上风电机组大型化发展趋势，需要观测风资源高度逐步提高。

复杂的水深地质条件、观测高度的提升以及有限的海上施工窗口期，使得建海上测风塔设计建设难度增大、施工周期变长、建设成本增加，难以满足海上风电开发建设进度要求。

而目前国内外海上漂浮式激光雷达测风技术仍处于技术验证阶段，还没实现完全商业化，技术适用性和可靠性有待进一步验证。

当前国内外众多研究机构推出基于不同物理模型和全球不同观测站实测数据得出了不同尺度下的风资源图谱产品[3]，如中国气象局近海风资源数据集、GlobalWindAltas全球风能资源数据集、3Tier中尺度数据集等，对无实测数据区域尤其是海上风资源评估具有重要的参考意义。

法国美迪公司开发的Meteodyn WT是一款可供风电场进行风资源评估的流体力学数值模拟软件，能够实现不同结果点位置不同高度的风资源评估。

本文在此背景下，提出基于风资源图谱成果的风速分析法、基于3Tier中尺度数据集的风速分析法和基于WT的风速分析法等三种方法对海上无测风区域进行风速评估，为海上风电场风资源评估分析提供参考。

• （10）猝发与拟序结构。这也是近代湍流研究的重大发现，试验 表明，在湍流混合层和剪切湍流边界层中存在大尺度的相干结构 和猝发现象，说明湍流不是完全无秩序、无内部结构的运动。这 促使人们改变了对湍流的某些传统观念。
从层流到湍流（一）
从层流到湍流 Corssin and Karweit 1969
从层流到湍流（二）
• （6）连续性。湍流是一种连续介质的运动现象，因此满足 连续介质力学的基本规律，例如N-S方程；
• （7）流动特性。湍流不是流体的物理属性，而是流动的运 动性质，所以不同的流体其湍流特征往往也不一样，例如 边界层湍流与尾迹湍流，正因为如此（湍流依赖于外部条 件，如边界条件），所以工程上很难对湍流进行统一的模 式处理，但是湍流的一些本质特征是普适的，寻找这些普 遍规律正是湍流理论研究的中心任务；
ui ui ui

q q q

p p p
平均场方程描述长时间过程， 脉动场方程描述短时间过程。
以状态方程为例：
p Ra
p'




'
T
v

Tv
进行雷诺平均后：
p Ra


Tv
'Tv '
右边第二项很小可以略去不计：
p Ra
u3 w
3、动量守恒方程
ui t
uj
ui x j
i3g 2ijk juk

1

p xi

2ui
x
2 j
ⅠⅡ Ⅲ
Ⅳ
ⅤⅥ
存 平 直重 储 流 方力 项 传 向项
输 作， 项 用仅
在 垂

大气湍流传输特性研究及模拟方法摘要：大气湍流传输特性是气象学和大气科学研究的重要课题之一，对于空气质量、气候变化以及环境污染等问题有着重要的影响。

本文将探讨大气湍流传输特性的研究方法和模拟方法，介绍湍流的形成机制和传输过程，并且分析现有的模拟方法及其应用领域。

通过深入研究大气湍流传输特性，有助于我们更好地了解和预测大气环境的变化。

1. 引言大气湍流传输特性是指大气中的湍流现象对物质传输的影响。

湍流是指流动介质中的无规则运动，具有剧烈的速度波动和能量传递。

大气湍流传输特性的研究对于理解大气环境变化以及对空气质量的影响具有重要意义。

2. 湍流的形成机制湍流的形成机制主要由湍流能量输入和湍流能量耗散两个过程共同决定。

湍流能量输入包括大气层之间的能量交换以及地表和大气之间的能量交换。

湍流能量的耗散主要通过黏性阻尼来实现。

3. 湍流传输过程湍流传输过程包括湍流扩散、湍流对流和湍流输运等几个方面。

湍流扩散主要指物质在湍流的作用下的横向扩散过程，湍流对流指因湍流引起的物体内部和外部的无规则运动，湍流输运是指物质在湍流流场中由于速度非均匀性而发生的输运。

4. 大气湍流传输特性的模拟方法现代大气科学研究中，模拟方法是一种常用的研究手段。

常见的大气湍流传输的模拟方法包括数值模拟和实验模拟两种。

4.1 数值模拟方法数值模拟方法是通过数学模型来模拟大气湍流传输特性。

常用的数值模拟方法有雷诺平均N-S方程模型（RANS）、大涡模拟（LES）和直接数值模拟（DNS）三种。

RANS方法是常用的湍流模拟方法，Les方法相对而言更加精确，DNS方法则是最为精确的模拟方法。

4.2 实验模拟方法实验模拟方法是通过实验设备来模拟大气中的湍流传输过程。

常见的实验模拟方法包括风洞实验、水洞实验和人工湍流实验等。

这些实验模拟方法可以精确测量大气湍流传输特性的相关参数，为湍流传输特性的研究提供了可靠的数据。

5. 模拟方法的应用领域大气湍流传输特性的模拟方法在空气质量、气候变化和环境污染等研究领域有着广泛的应用。

第二章 大气边界层湍流基础
大气边界层湍流基础
第一节 平均场与湍流场 第二节 湍流特征量及基本统计学方法 第三节 大气湍流谱 第四节 大气湍流通量及输送 第五节 大气湍流动能(TKE)
湍流运动特征
三维，非线性，涡旋运动——耗散性， 即湍流运动能量以非线性方式由大湍 涡向小湍涡传递，最后耗散于分子热 能运动
(1)采样时段泄露效应 (2)采样间隔混迭效应
(1)在有限时段上采样，谱S(n) 被修改成ST(n)：
ST(n) S(f)s in (T n( n f)f)df
修改后的谱ST(n)中带有虚假的高频成份 T 越大， ST(n) 越接近真实谱S(n) T 越小，泄露影响越大, ST(n)与S(n)之间的差别越大
应用于空间某一特定点，用变量A在某一时间 周期τ上的积分或总和表示：
A=A(t,s), t : 时间; s: 空间
1
At(s)
A(t,s)dt
t0
1 N1
At(s)Ni0 A(i,s)
离散情况下：t/N
2 空间平均
对某一固定时间t，用变量A在空间区域S上的 积分或总和来表示：
如果两个变量完全相关(即变化方向一致)，则 r =+1 如果完全负的相关(即反方向变化)，则 r =-1 如果两变量变化不相关，则 r =0
自相关测量某一波动在某一时间序列或空 间序列总体上的持续性。因为规则变化可 能与诸如涡动等物理现象有关，因此在序 列中确定持续波或振荡的可能性是特别有 用的。
（2）消除趋势项（Trend removed，去倾），周期大 于样本记录长度的频率成分称为趋势项。在湍流分析中， 通常需要消除这种大尺度的影响，因为去倾与否得出的相 关函数和谱在低频部分有明显区别。