《中小尺度数值模拟》复习.docx
- 格式:docx
- 大小:73.38 KB
- 文档页数:4
目录Chapter1 Introduction (2)1.Water wave phenomena(水波现象) (2)Chapter2 Theories of typical wave models (2)2.Wave deformations in coastal region (2)3.Main types of wave model (2)4.Modeling for different engineering interests (2)5.Models based on solution to the N-S equation (2)6.Nonhydrostatic model (7)7.Models based on potential flow solver (8)8.Boussinesq type models (9)9.Shallow water wave (=long wave) equation (13)10.Models based on mild-slope equation (MSE) (16)11.Phase-averaged models (17)12.Summary of various kinds of wave models (17)Chapter3 Fundamentals of finite difference method (18)13.Typical forms of partial differential equation (PDE) for fluid flow (18)14.Basic concepts of finite difference method (FDM) (18)15.Modified Partial Differential Equation (MPDE) (21)16.Stability(稳定性) and Convergence(收敛性) (22)17.Discussions on the numerical stability of various FDE schemes for specific schemesbased on calculation (25)18.Fourier analysis for numerical stability (26)19.Upwind concept (27)20.Monotonicity(单调性) (28)21.Conservation (28)22.Further discussions on the properties of various FDE schem (28)23.Well-posed problem (30)Chapter4 Fundamentals of FVM (Finite V olume Method) (30)24.Differences among FDM, FVM and FEM (30)25.FVM scheme (30)26.Gridding system (33)27.An example of FVM (34)Chapter5 Solution method to some typical wave models (35)28.Model based on solution to N-S equations (35)29.Models based on Boussinesq equation (36)30.Models based on Parabolic equation (37)Chapter1 Introduction1. Water wave phenomena (水波现象)Concept: water wave is a kind of water flow oscillation with free surface.driven by various kind of forces; sustained by gravity.Chapter2 Theories of typical wave models2. Wave deformations in coastal regionWave diffraction(绕射)Wave refraction (折射):波峰线与海岸线成一定夹角入射的波浪,波峰线会逐渐旋转,使得波峰线逐渐平行于海岸线。
第一章中尺度天气系统的特征1、中尺度气象学:水平尺度: 10-1000km对象:中尺度环流系统内容:中尺度环流系统的结构、形成和发展演变规律、机制及其分析预报方法意义:①许多灾害性天气(如暴雨、大风、冰雹、龙卷等)都是由中小尺度系统造成的。
②中尺度气象学是甚短期预报和临近预报的理论基础。
(长期>10天,中期3-10天,短期1-3天,甚短期0-12h,临近0-2h)③中尺度环流系统是大气环流重要成员(大尺度背景场依存条件)2、天气系统的尺度划分:(一)经验分类法(经典方法)小尺度系统(雷暴、龙卷)和大尺度系统(锋面、气旋)中尺度系统(飑线、中气旋等)(二)动力学定义可利用罗斯贝数(Ro)和弗劳德(Froude)数(Fr)来描述大气的时空尺度。
Ro = U/fL (惯性力/柯氏力);Fr=U2/gL(△ρ/ρ)(惯性力/浮力)(三)实用(几何)分类3、中尺度大气运动的基本特征(1)尺度:水平尺度在2-2000km之间,时间尺度在几十分钟至几天之间。
范围很宽。
性质不同。
(2)散度、涡度、垂直速度:取V~10m/s,H~10km,对α,β,γ中尺度W分别为10-1m/s, 100m/s和 101m/s,垂直速度、散度、涡度都比大尺度运动大1到几个量级。
(3)地转偏向力和浮力的作用:中尺度运动中,地转偏向力和浮力的作用都必须考虑。
大尺度运动:地转偏向力重要,浮力可略小尺度运动:浮力重要,地转偏向力可略中尺度运动:地转偏向力和浮力都考虑(4)质量场和风场的适应关系:质量场(气压场)适应风场。
大尺度运动: 风场适应质量场(气压场)。
中尺度运动: 质量场(气压场)适应风场。
第二章地形性中尺度环流1、中尺度大气环流系统分为:地形性环流系统、自由大气环流系统2、地形波:一般把气流过山所引起的气流称为地形波。
3、地形波的基本类型:层状气流(山脉波):山脉上空的平滑浅波 ,风小。
驻涡气流(驻涡):山脉背风面的半永久性涡旋,山顶以上风速大。
数值天气预报第一章1、名词解释数值天气预报:所谓数值天气预报,就是在给定初始条件和边界条件的情况下,数值求解大气运动基本方程组,由已知的初始时刻的大气状态预报未来时刻的大气状态。
因此,大气运动基本方程组是制作数值天气预报的基础。
初始条件:初始条件就是初始时刻各因变量(即气象要素)的空间分布。
其一般形式为u=u(x,y,z,O) v=v(x,y,z,O)t=0, w=w(x,y,z,O)p=p(x,y,z,O)T=T(x,y,z,O)边界条件:边界条件就是所研究区域的大气边界上气象要素应满足的条件。
研究全球范围的大气运动,如果大气内部各气象要素都是连续的,则只需给出大气的下边界条件和上边界条件;如果大气内部存在不连续面,则还需给出内边界条件。
尺度分析:所谓尺度分析就是根据某种类型运动的特征尺度来估计基本方程组中各项数量级的大小,从而使方程组得到简化的一种方法。
特征尺度:物理变量的特征尺度是指某种类型运动所占据的空间范围、维持的时间、各场变量及时空变化的典型值。
大气模式:2、问答四种坐标系的优缺点?第二章1、名词解释地图放大因子:映像比例尺m是映像平面上的距离与地球表面上相对应距离的比值,称之为地图放大系数或地图放大因子。
正形投影:投影光源位于球心,映像面为圆锥面,映像面圆锥角为a (0<a<180),标准纬度为<Po 2、问答特点% —- *第二早1、名词解释平滑:所谓平滑就是用某点周围若干点的值进行加权平均来代替该点的值,经过这样处理的物理量场可以衰减或者滤掉短波分量。
2、问答什么是差分格式的相容性、收敛性和稳定性?它们之间的关系如何?答:截断误差是否随着网格距和时间步长趋于零而趋于零,称为解的收敛性问题。
舍入误差是否随着网格距和时间步长趋于零而在整个求解区域内保持有界,称为解的稳定性。
相容性、收敛性和稳定性之间的关系,即为克拉斯等价定理:对于一个线性微分方程的适定初值问题,若其差分方程和微分方程是相容的,则稳定性是收敛性的充分必要条件。
水锤的数值模拟
一GAMBIT部分
第一步:打开;
第二步:;
第三步:依次打开,,,
;点击apply,得
第四步:划分网格,依次点击,,得
,点击apply,得
依次点击,,得下图
第六步:保存msh文件。
,得
,点击。
二fluent计算
第一步:启动
第二步:读取shuichui.msh文件。
第三步:网格检查,
第四步:算法,define---solver
第五步:选取湍流模型
第六步:选择材料
第七步:设置工作环境
,默认设置。
第八步:UDF
第九步:设置边界条件
A:进口条件
B:出口条件
第十步:求解控制参数设置,依次点击solve—controls---solution,
保持默认设置就可以了。
第十一步:打开残差图,点击plot,其他默认。
第十二步:初始化。
依次点击solve—initialize—initialize,在compute form 下拉表中选择inlet,依次点击init,close就可以了!
第十三步:动画设置
第十四步:迭代
三后处理。
中微尺度耦合
【原创版】
目录
1.中微尺度耦合的定义和背景
2.中微尺度耦合的研究方法和应用
3.中微尺度耦合的挑战和未来发展
正文
中微尺度耦合是一种研究气象和环境过程的方法,它关注的是天气和气候系统中的中小尺度过程。
这些过程对于天气和气候系统的演变有着重要的影响,但由于其尺度较小,传统的观测和模拟手段难以准确捕捉。
近年来,随着计算机技术和观测手段的进步,中微尺度耦合的研究得到了广泛的关注和发展。
中微尺度耦合的研究方法主要包括观测、数值模拟和统计分析。
观测是通过各种手段,如气象站、雷达、卫星等,获取大气和环境参数的实时数据。
数值模拟则是使用数值模型,如数值天气预报模型和气候模型,模拟大气和环境过程。
统计分析则是通过统计方法,如相关分析、回归分析等,分析观测数据和模拟结果之间的关系。
中微尺度耦合的应用主要包括天气预报、气候预测和环境监测。
天气预报是指通过分析中小尺度过程,提高天气预报的准确性。
气候预测是指通过研究中小尺度过程对气候系统的影响,提高气候预测的准确性。
环境监测是指通过观测和分析中小尺度过程,监测环境质量的变化。
尽管中微尺度耦合的研究已经取得了一定的进展,但仍然面临着一些挑战。
首先,中小尺度过程的观测和模拟难度大,需要进一步提高观测和模拟能力。
其次,中微尺度耦合的研究方法和应用还需要进一步完善和优化。
最后,中微尺度耦合的研究还需要更好地与其他气象和环境研究领域相结合,以提高研究的全面性和深入性。
未来,中微尺度耦合的研究将会继续发展,以更好地理解和预测天气和气候系统的演变。
航空数值天气预报模式下风场及风切变小尺度数值模拟蒋立辉;秦宇焘;杨依莹;王智【摘要】针对风切变这一影响航空飞行的主要天气现象,进行小尺度下的风场数值模拟研究.实验模拟时段为协调世界时2014年6月8日06时至2014年6月8日24时(9日00时),对研究区域内生成并过境的强对流天气系统进行中小尺度的绝对涡度、水平风场和垂直风场的模拟;并提取小尺度下的风速风向值,尝试对风切变提取.结果表明:①数值模式能够清晰立体地模拟出这一次强对流天气系统的生成和运动的整体状况;②小尺度下能够对该强对流天气系统的局部进行精细化的立体的模拟;③通过分析小尺度下的数值模拟输出产品,能够判别出风切变;④对于风向数据,采用余弦平滑法能够进行有效的质量控制观察分析.【期刊名称】《科学技术与工程》【年(卷),期】2016(016)013【总页数】6页(P156-161)【关键词】数值模式;风切变;强对流天气系统;小尺度【作者】蒋立辉;秦宇焘;杨依莹;王智【作者单位】中国民航大学空中交通管理学院,中国民航大学民用航空气象技术研究所,天津市智能信号与图像处理重点实验室,天津300300;中国民航大学空中交通管理学院,中国民航大学民用航空气象技术研究所,天津市智能信号与图像处理重点实验室,天津300300;中国民航大学空中交通管理学院,中国民航大学民用航空气象技术研究所,天津市智能信号与图像处理重点实验室,天津300300;中国民航大学空中交通管理学院,中国民航大学民用航空气象技术研究所,天津市智能信号与图像处理重点实验室,天津300300【正文语种】中文【中图分类】P456.7风切变是指逆风或顺风出现超过几秒钟的突然转变,导致飞机升力的变化。
它不仅可以使飞机跟踪偏差;也可能让飞机失去稳定。
风切变对飞机起飞或降落来讲,是一个潜在的危险因素;也是国际航空领域和气象界一直认可的事实;这是一个很难解决气象问题。
1976年以后,国际民航组织和世界气象组织积极开展对风切变的研究工作,起到了重要的组织作用,明确表示把风切变列为重要研究课题。
隐式( Implicit ):对于给定变量,单元内的未知值用邻近单元的已知和未知值计算得出。
因此,每一个未知值会在不止一个方程中出现,这些方程必须同时解来给出未知量。
隐式求解法将冲压成型过程的计算作为动态问题来处理后,就涉及到时间域的数值积分方法问题。
在80年代中期以前,人们基本上使用牛曼法进行时间域的积分。
根据牛曼法,位移、速度和加速度有着如下的关系:上面式子中,分别为当前时刻和前一时刻的位移,和为当前时刻和前一时刻的速度,和为当前时刻和前一时刻的加速度,β和γ为两个待定参数。
由上式可知,在牛曼法中任一时刻的位移、速度和加速度都相互关联,这就使得运动方程的求解变成一系列相互关联的非线性方程的求解。
这个求解过程必须通过迭代和求解联立方程组才能实现。
这就是通常所说的隐式求解法。
隐式求解法可能遇到两个问题。
一是迭代过程不一定收敛;二是联立方程组可能出现病态而无确定的解。
隐式求解法的最大优点是它具有无条件稳定性,即时间步长可以任意大。
显式求解法如果采用中心差分法来进行动态问题的时域积分,则有如下位移、速度和加速度关系由上式可以看出,当前时刻的位移只与前一时刻的加速度和位移有关,这就意味着当前时刻的位移求解无需迭代过程。
另外,只要将运动方程中的质量矩阵和阻尼矩阵对角化,前一时刻的加速度求解无需解联立方程组,从而使问题大大简化,这就是所谓的显式求解法。
显式求解法的优点是它即没有收敛性问题,也不需求解联立方程组,其缺点是时间步长受到数值积分稳定性的限制,不能超过系统的临界时间步长。
由于冲压成型过程具有很强的非线性,从解的精度考虑,时间步长也不能太大,这就在很大程度上弥补了显式求解法的缺陷。
在80年代中期以前显式算法主要用于高速碰撞的仿真计算,效果很好。
自80年代后期被越来越广泛地用于冲压成型过程的仿真,目前在这方面的应用效果已超过隐式算法。
显式算法在冲压成型过程的仿真中获得成功应用的关键,在于它不像隐式算法那样有解的收敛性问题。
第一章中尺度天气系统的特征1、中尺度气象学:水平尺度: 10-1000km对象:中尺度环流系统内容:中尺度环流系统的结构、形成和发展演变规律、机制及其分析预报方法意义:①许多灾害性天气(如暴雨、大风、冰雹、龙卷等)都是由中小尺度系统造成的。
②中尺度气象学是甚短期预报和临近预报的理论基础。
(长期>10天,中期3-10天,短期1-3天,甚短期0-12h,临近0-2h)③中尺度环流系统是大气环流重要成员(大尺度背景场依存条件)2、天气系统的尺度划分:(一)经验分类法(经典方法)小尺度系统(雷暴、龙卷)和大尺度系统(锋面、气旋)中尺度系统(飑线、中气旋等)(二)动力学定义可利用罗斯贝数(Ro)和弗劳德(Froude)数(Fr)来描述大气的时空尺度。
Ro = U/fL (惯性力/柯氏力);Fr=U2/gL(△ρ/ρ)(惯性力/浮力)(三)实用(几何)分类3、中尺度大气运动的基本特征(1)尺度:水平尺度在2-2000km之间,时间尺度在几十分钟至几天之间。
范围很宽。
性质不同。
(2)散度、涡度、垂直速度:取V~10m/s,H~10km,对α,β,γ中尺度W分别为10-1m/s, 100m/s和101m/s,垂直速度、散度、涡度都比大尺度运动大1到几个量级。
(3)地转偏向力和浮力的作用:中尺度运动中,地转偏向力和浮力的作用都必须考虑。
大尺度运动:地转偏向力重要,浮力可略小尺度运动:浮力重要,地转偏向力可略中尺度运动:地转偏向力和浮力都考虑(4)质量场和风场的适应关系:质量场(气压场)适应风场。
大尺度运动: 风场适应质量场(气压场)。
中尺度运动: 质量场(气压场)适应风场。
第二章地形性中尺度环流1、中尺度大气环流系统分为:地形性环流系统、自由大气环流系统2、地形波:一般把气流过山所引起的气流称为地形波。
3、地形波的基本类型:层状气流(山脉波):山脉上空的平滑浅波,风小。
驻涡气流(驻涡):山脉背风面的半永久性涡旋,山顶以上风速大。
中小尺度数值模拟
第二章控制大气的基本方程组和地图投影
1、数值模式下的平均运动方程组
——= gP + g — ZQx% + F dt p 沦皿0
吧于心Q 肆亠
二
PRT
dt p
2、次网格过程
用格点值不能直接描述的过程。如湍流,积云对流一对流扰动(空间)。平均值对应网格点
值,扰动值对应次网格过程。
3、 气象上常用的地图投影有哪些?在数值模拟过程中如何确定使用何种地图投影?
(1) 极射赤面投影:高纬地区变形小,多用作极地天气图和北半球天气图的底图。(k=l)
(2) 兰伯特投影:中纬地区变形小,适用于中纬度地区天气图的底图,如亚欧天气图。
(3) 麦卡托投影:低纬地区变形小,适用于作低纬地区天气图的底图。(k=0)
第三章方程组的离散一差分法
4、 中小尺度现象的模拟方法
(1) 实验室模型:动力相似性
(2) 解析解:数学方法求解方程(连续)
(3) 数值模式:离散数值求解方法
5、 数值求解的主要方法及其适用范闱
有限差分法:利用Taylor级数展开,使用最广泛。
有限元法:根据能量最小原理,将积分区域划分为有限的、不重叠但互相连接的单元,每个 单
元选择基函数,用单元基函数的线性组合逼近单元屮的真解,整体区域的基函数可以看作 由每
个单元基函数组成,整个区域的解可以看作由所有单元的近似解,适用于不规则区域。 谱方
法:Fourier变换(谱展开),主变换要用于全球模式,计算精度髙,现已经在屮尺度模 式中应
用。
6、 什么是差分方程的一致性、收敛性和线性稳定性?
一致性(相容性):即差分方程是否逼近微分方程(当步长一0)
收敛性:在一定的定解条件下,差分方程的解U是否逼近微分方程的解u。
稳定性:在时间积分过程中,当时间步长趋向于0时,在整个求解区域内,舍入误差U-U
保持有界则是稳定的。即差分解的误差不随时间增长。
7、 什么是差分方程的非线性计算不稳定?产生的原因是什么?
答:在满足线性稳定性条件下,由于非线性作用而产生的不稳定,称为非线性不稳定。非线 性
计算不稳定主要表现为短波能量的虚假增长。
产生的原因:①可能是混淆误差②可分辨波能量Z间的相互转换屮的误差③还没有完全弄 清楚
8、 混淆误差
有限网格系统能分辨的波的最短波反为2心,对于非线性作用产生的波长小于2心的波动, 网
格系统不能正确地分辨,而把它错误地表示成为某种波长大于2心的波,从而造成了这 种波的
误差一一混淆误差。
9、 如何消除差分方程的非线性不稳定?
(1) 空间、时间平滑,滤去短波分量,滤去波长v 4心的波可消去alias解。
(2) 加水平扩散项rV2A, A为物理量,i•为扩散系数,可根据作用大小选择。
(3) 构造具有隐式平滑式某种选择性衰减作用的差分格式。
以上方法优点:实用有效;缺点:克服短波同时,也会削弱长波,影响预报效果。
(4) 构造差分格式:尽量保持连续系统的一些积分约束的差分格式。能量守恒差分格式得
总能量约束。涡度拟能守恒差分格式得能谱约束,抑制能量向短波转移。
10、为何能量只会在短波上发生虚假的自反馈?
S
波能量发生变化是虚假的,不反映真实过程。 由于差分近似和非线性相互作
用引起的。
不断重复过程,4心波能量可能会迅速增长。
=>非线性计算不稳定
可证:能量只有在短波上发生虚假的自反馈:
第四章离散网格、资料处理与初始化
it、数值模拟中模拟区域和格距大小的设置要求有哪些
区域设置:(1)中小尺度过程引起的变量的有意义变化必须被包括在模式空间中,即模式范
围要足够大,从而边界上基本不受中小尺度过程影响。
(2) 外强迫的尺度和范围(如地形、加热等)
(3) 对强迫的响应尺度
格距设置:(1)中小尺度系统的特征必须能足够被分辨,即中小尺度系统至少要占4个格距。
若宀W
(2 心)*(4
心)
非线性相互作用
则
S = —
4
(4
心)
(2) 设置的格距要适合计算方案的稳定性要求。
(3) 计算代价最小,必须考虑计算资源的占用问题。
12、 水平网格结构
① 等距网格。心,都是常数,优点:简单
② 拉伸网格。心,不是常数,视需要调整。优点:节省计算资源,同样的网格数可以有 更大的
模拟区域。
③ 嵌套网格。分粗网格和细网格,优点:节省计算资源,通过多层嵌套,可以实现高分辨率
的模拟。
④ 移动网格。可跟随着中尺度系统的发展來移动,优点:简单,计算代价小,适用移动范围
较大的系统的模拟。
13、 何为交错网格?为什么垂直方向要采用拉伸网格?
交错网格:将不同的气象要素有规律的写在不同网格上。
大部分中尺度天气现彖的强迫主要位于对流层低层大气,绝大多数的数值模式在垂直方向上 都
采用拉伸网格结构。
14、 为什么由客观分析得到的格点值并不能直接作为初始场
观测资料或分析资料的误差会破坏风场和气压场之间的不平衡。初始资料和模式之间可能不
协调,如果不进行初始化处理,在模式积分的一段容易出现虚假的高频波。
15、 初始化方法
客观分析方法:有限元方法,逐步订正法,九点拉格朗tl插值,最小二乘法等。
另外,由于非定时探空资料日益增多,客观分析还包括了卫星资料、雷达资料等多种数据的
使用。
16>边界设置
海绵边界条件,即在边界附近设一缓冲帯,一般三至五圈,设为模式中变量,有模式计算输
出的局地变化值为经海绵边界处理后实际采用的值。
海绵边界条件的实质是在边界上对固定的外边界和内部的预报场之间起一种平滑协调的作
用,它可以对向外传播的惯性重力波较好地吸收,减小边界上波的反射。
对于深对流系统的模拟,一般取上边界在对流层顶或平流层下层。
对于浅对流系统的模拟,上边界可以取在对流层中层,模式层顶的选取主要依赖于实际模拟
系统的系统性质。
17、大气边界层的作用
① 厚度:1〜2km (分三层)
② 主要过程是湍流摩擦(分子摩擦较小,可略)
③ 边界层内动量、热量和水汽的湍流垂直通量最大
④ 是大气重要的动量汇,热量和水汽的源。
因此,大气通过边界层的湍流输送而和下垫面发生热力和动力作用。
第五章数值模拟参数化
18、 参数化一一用大尺度变量表征次网格或小尺度作用总体效应的方法。
积云对流参数化:将人尺度模式不能显式分辨的对流凝结和对流引起的热量、水分和动量的 输
送与模式的预报变量联系起来。实际做法是用大尺度网格点上的变量值来表示次网格尺度 积云
的凝结加热以及垂直输送效应的作用。
枳云参数化方案(CPSs)通常要求在网格尺度未饱和时产生能垂直输送热量、水汽和其他 物理量
的次网格尺度隐式云。闭合假设需要确定隐式云和网格尺度变量Z间的关系。
一般来说,<3km,显式对流方案;>50km,隐式对流方案;3〜50km,显式+隐式。
19、 湍流和对流的区别
湍流一一发生在边界层,耗散动量、输送热量和水汽。
对流一一发生在对流层,对流凝结释放潜热加热大气,输送热量和水汽到中高层,也具有耗 散
作用。
20、 资料同化
定义:将不同时间、不同时间间隔和不同地区通过不同方式观测収得的资料,在数理方程的 约
束下组合成为统一的资料系统。
方法:三维变分、四维变分、卡尔曼滤波法
实习内容
1、 使用WRF模式进行数值模拟的具体步骤
(1) 修改 namelist.wps 中的&share 和&geogrid 部分,运行 geogrid.exe
(2) 链接Vtable,链接外部数据fnl,运行ungrib.exe
(3) 运行metgrid:将气象要素场水平插值到地形资料上
(4) 链接met数据,运行real.exe,运彳亍wrf.exe
2、 WRF模式中每个进行步骤的作用是什么
(1) 用geogrid处理地形资料和网格数据
(2) 解码GRIB资料,得到气彖要素场
(3) 将抽取出来的气象场水平插值到模式格点
(4) 将气象场垂直插值到WRF的eta层中
(5) 生成模式预报
3、 WRF模式中如何设置水平网格和垂直分层
(1) 在namelist.wps中的&geogrid部分部分先设置母区域的中心经纬度、参考经度、网格 数和
格距,如果嵌套,再设置子区域在母区域中的位置和网格数
(2) 在namelist.input中通过e_vert设置垂直层数