下载文档原格式
实习报告实习单位：XX大学气象研究所实习时间：20XX年XX月XX日—20XX年XX月XX日实习内容：中尺度数值模拟一、实习背景及目的随着科学技术的不断发展，数值模拟技术在气象领域得到了广泛的应用。
中尺度数值模拟作为气象研究中的一种重要方法，可以有效地模拟和研究气象现象及其演变过程。
本次实习旨在通过实际操作，掌握中尺度数值模拟的基本原理和方法，提高自己在气象领域的实际工作能力。
二、实习内容及过程1. 实习前的准备工作在实习开始前，我参加了由导师组织的中尺度数值模拟培训课程，学习了中尺度数值模拟的基本原理、方法和实际应用。
此外，我还通过查阅相关文献，对中尺度数值模拟的研究现状和发展趋势有了更深入的了解。
2. 实习内容实习期间，我主要进行了以下几个方面的工作：（1）了解并熟悉中尺度数值模拟所使用的软件和硬件设备，掌握相关操作技巧。
（2）根据实际气象资料，设置数值模拟的初始条件和边界条件。
（3）运行中尺度数值模拟软件，观察模拟过程中气象参数的变化，并对模拟结果进行分析。
（4）根据模拟结果，撰写实习报告，总结中尺度数值模拟的特点和不足。
3. 实习过程在实习过程中，我首先选择了某一地区的一次降水过程作为研究对象。
根据实测气象资料，我设置了数值模拟的初始条件和边界条件，并使用中尺度数值模拟软件进行模拟。
在模拟过程中，我密切关注气象参数的变化，并对模拟结果进行了详细的分析。
通过实习，我了解到中尺度数值模拟在气象研究中的应用广泛，可以模拟各种气象现象，如降水、台风、沙尘暴等。
同时，中尺度数值模拟还能够为气象预报提供重要的参考依据。
三、实习收获及体会通过本次实习，我对中尺度数值模拟的基本原理和方法有了更加深入的了解，提高了自己的实际操作能力。
同时，实习过程中我对气象研究的重要性有了更深刻的认识，为自己的未来研究方向和方法的选择提供了宝贵的经验。
实习期间，我也发现中尺度数值模拟存在一些不足之处，如模拟结果受初始条件和边界条件的影响较大，有时难以准确反映实际气象状况。
中小尺度中尺度带状对流系统由对流单体侧向排列而成的中尺度对流系统一般称为带状对流系统。
结构：飑线作为一个中尺度系统，应包括对流区和非对流(层状云)区两部分。
对流区包含强烈的、垂直延伸的强回波核，而层状区域由一些降水构成均匀（不是绝对均匀）纹理。
概念：为—种带(或线)状中尺度系统,是非锋面性狭窄的活跃的雷暴带(或不稳定线)。
其中有许多雷暴单体（其中包括若干超级单体）侧向排列而形成的，是风向、风速气压、温度等突变的狭窄强对流云带。
为破坏力严重的灾害性天气。
飑线的一般特征○!发生地点：出现在中纬度的某些大陆地区以及主要的热带大陆和热带海洋地区。
温带地区的飑线常发生在春夏之交的过渡季节，有的出现在冷锋前或气旋波的暖区，有的在冷锋后的冷区里，还有在冷暖锋上或切变线(辐合线)附近生成的。
并大致与锋面相平行。
○2尺度：长约几百千米，宽度约50～100km。
飑线由若干“飑段”组成。
每个飑段包含若干大而孤立的相互分离的风暴。
○3时间尺度：几小时至十几小时。
○4飑线的地面要素场的结构：飑线由雷暴单体侧向排列而成，每个单体成熟期都有地面冷丘及水平外流和阵风锋；小冷丘和阵风锋结合起来形成小尺度雷暴高压和阵风锋。
阵风锋又称为飑锋：处在雷暴高压边缘。
具有很强的温度梯度、气压梯度和风速和风向水平切变，它也叫气压涌升线或跳跃线。
○5过境特征：由于飑锋附近是各种气象要素水平梯度很大的地带，因此当飑锋过境时，气象要素将发生急剧的变化。
通常表现为气压涌升、气温急降、风向突变、风速剧增以及强烈降水等。
飑线前低压：飑锋前方一般有中尺度低压。
它的形成可能与飑线前方高层的补偿下沉气流引起的绝热增温有关。
尾流低压：雷暴高压后方的中尺度低压，它的形成与雷暴高压后部的尾流效应相联系。
飑中系统：包括飑线、飑线前低压、雷暴高压以及尾流低压统称为飑中系统。
※飑中系统的全部系统一般只在成熟阶段才同时出现。
不同阶段系统的强度和结构是不同的。
两类比较常见飑线：1、具有前导对流线和尾随层状云区以及具有由前向后和由后向前两支入流的飑线发生在风垂直切变相对小的环境中的飑线飑线的前方有一支由前向后的入流迎着飑锋上升，到高层分裂成向前和向后的两支气流，其后部中层则另有一支由后向前的入流。
1.背风波：当风速随高度增大时，则可在背风坡出现波动气流，这种波动成为背风波。
2.多单体雷暴：由一些处于不同发展阶段的生命期短暂的对流单体组成，是具有统一环流的雷暴系统3.龙卷风暴：产生龙卷的强风暴系统称为龙卷风暴。
4.温带飑线：为—种带(或线)状中尺度系统,是非锋面性狭窄的活跃的雷暴带(或不稳定线)。
其中有许多雷暴单体（其中包括若干超级单体）侧向排列而形成的，是风向、风速、气压、温度等突变的狭窄强对流云带。
为破坏力严重的灾害性天气。
5.对流复合体（MCC）：指由若干对流单体或孤立对流系统及其衍生的层状云系所组成的对流系统，其空间尺度和时间尺度具有幅度很广的谱。
最简单的是二维的线状对流系统，最大而复杂的是一种具有近于圆形团状结构的MCC这两种系统位于对流复合体波谱的两端。
6. 对称不稳定：在流体静力、地砖平衡且具有水平切变的情况下，浮力和旋转会共同起作用，这两种效应会导致一种新的浮力惯性不稳定，即对称不稳定，对称不稳定是中尺度雨带与雨团形成的主要不稳定机制。
7.条件性不稳定:对干空气是静力稳定的，而对饱和湿空气静力不稳定的情况。
8.对流性不稳定：不论气层原先的层结性如何，在其被抬升达到饱和后，如果是不稳定的则称对流性不稳定。
9.第二类条件性不稳定：大尺度流场通过摩擦边界层的抽吸作用，为积云对流提供了必须的水汽辐合与上升运动，反过来积云对流释放凝结潜热又成为驱动大尺度扰动所需要的能量，于是小尺度积云对流和大尺度流场通过相互作用，相辅相成的都得到了发展。
这种通过不同尺度运动的相互作用使对流和大尺度流场不稳定增长的物理机制称为第二类条件性不稳定。
10. 超级单体风暴：直径达20~40km 以上，生命期达数小时以上，即比普通的成熟单体雷暴更巨大、更持久、天气更猛烈的单体强雷暴系统。
它具有近于稳定的、高度有组织的内部环流，并且连续地向前传播可达数百公里。
11. .暖输送带：在槽前辐合区的边界上通常可以看到一支狭长的云带。
中尺度大气数值模拟及其进展中尺度大气数值模拟及其进展中尺度大气数值模拟是指对中尺度大气运动、湍流、边界层、云微物理、辐射传输等过程进行数值模拟的一种方法。
近年来，随着计算机技术的快速发展和观测技术的不断进步，中尺度大气数值模拟的研究已经取得了许多重要的进展，对于气象预报、气候变化研究和环境污染预测等方面都起到了重要的作用。
中尺度大气数值模拟的目标是通过计算空间和时间上的大量物理量，来模拟和预测中尺度大气运动过程。
中尺度大气运动是指介于大尺度天气系统和小尺度湍流系统之间的系统，其典型特征是空间尺度在几十公里到几百公里之间，时间尺度在几分钟到几小时之间。
中尺度大气运动包括了许多重要的现象，如大气锋面、对流云团、飑线等，对于气象预报和气候变化研究具有重要的意义。
中尺度大气数值模拟的基本原理是通过数值方法将大气方程离散化，并通过数值解算得到大气运动的演化过程。
其中，最常用的模型是基于Navier-Stokes方程的大气动力学模型，通过有限差分、谱方法等数值技术对方程进行求解。
此外，为了更好地模拟大气过程，中尺度大气数值模拟还必须考虑到湍流的影响，湍流参数化是其中的关键技术之一。
近年来，随着计算机技术的不断进步，中尺度大气数值模拟的能力也得到了极大的提高。
传统的数值模拟方法需要通过将整个大气划分成若干个网格，然后分别对每个网格进行计算，这种方法在计算量和存储空间上都有较大的挑战。
为了克服这些问题，新型的数值模拟方法应运而生，如有限元方法、有限体积方法和伪谱法等。
这些方法可以更好地处理复杂的地形、不均匀的边界条件和非线性问题，提高了数值模拟的计算效率和精度。
除了数值方法的发展，观测技术的进步也为中尺度大气数值模拟提供了更多的观测数据，从而提高了数值模拟的准确性和可靠性。
现代大气观测技术，如雷达、卫星和飞机观测等，可以提供高时空分辨率的大气观测数据，在验证和改进数值模拟模型方面发挥重要作用。
此外，数据同化技术的应用也为中尺度大气数值模拟提供了新的思路和方法，通过将观测数据与数值模拟结果进行融合，可以进一步提高数值模拟的准确性和预报能力。
中尺度大气数值模拟及其进展中尺度大气数值模拟及其进展一、引言大气数值模拟是一种使用数学方程和计算机算法来模拟大气运动和气象现象的方法，它不仅能够帮助预测和研究天气、气候变化等现象，还可为决策提供重要参考。
在气象学研究领域，中尺度大气数值模拟被广泛应用，具有重要的意义。
本文将介绍中尺度大气数值模拟的基础理论和方法，并探讨其在气象学领域中的进展。
二、中尺度大气数值模拟的基础理论和方法中尺度指大气运动的空间尺度在几十到几百公里之间，时间尺度在几小时到几天之间。
中尺度大气数值模拟的基础理论是对大气运动和物理过程的基本方程进行数学化处理，建立相应的模型。
其中，最常用的模型是基于质量守恒、动量守恒、热量守恒和状态方程的Navier-Stokes方程。
为了简化计算，通常还采用了一些物理参数化方案，如湍流参数化、云微物理参数化等。
中尺度大气数值模拟的方法可以分为欧拉法和拉格朗日法。
欧拉法是在空间网格上离散化基本方程，通过数值迭代求解得到大气场的时空分布。
拉格朗日法则是跟踪气体的运动轨迹，通过将大气分成许多气团来模拟大气运动。
三、中尺度大气数值模拟在气象学领域的应用中尺度大气数值模拟在气象学领域有着广泛的应用。
首先，它可以用于天气预报，通过模拟大气运动，结合实时观测数据，可以提供准确的天气预报结果。
其次，中尺度大气数值模拟还可以用于研究气象灾害，如暴雨、台风等的形成机制和前后过程，从而为灾害预防和减轻提供科学依据。
此外，中尺度大气数值模拟还可以用于研究气候变化，如模拟气候系统中的能量和水分交换，探索气候变化的内在机制。
四、中尺度大气数值模拟的进展随着计算机技术的不断发展和模型改进，中尺度大气数值模拟在气象学领域取得了许多重要的进展。
首先，模拟精度显著提高，模型对大气物理过程的描述更加准确。
其次，模拟时间和空间分辨率不断增加，模拟结果更加细致。
此外，数据同化技术的应用使得模拟结果与实况数据更加吻合，提高了模式的可信度。
中小尺度中尺度带状对流系统由对流单体侧向排列而成的中尺度对流系统一般称为带状对流系统。
结构：飑线作为一个中尺度系统，应包括对流区和非对流(层状云)区两部分。
对流区包含强烈的、垂直延伸的强回波核，而层状区域由一些降水构成均匀（不是绝对均匀）纹理。
概念：为—种带(或线)状中尺度系统,是非锋面性狭窄的活跃的雷暴带(或不稳定线)。
其中有许多雷暴单体（其中包括若干超级单体）侧向排列而形成的，是风向、风速气压、温度等突变的狭窄强对流云带。
为破坏力严重的灾害性天气。
飑线的一般特征○!发生地点：出现在中纬度的某些大陆地区以及主要的热带大陆和热带海洋地区。
温带地区的飑线常发生在春夏之交的过渡季节，有的出现在冷锋前或气旋波的暖区，有的在冷锋后的冷区里，还有在冷暖锋上或切变线(辐合线)附近生成的。
并大致与锋面相平行。
○2尺度：长约几百千米，宽度约50～100km。
飑线由若干“飑段”组成。
每个飑段包含若干大而孤立的相互分离的风暴。
○3时间尺度：几小时至十几小时。
○4飑线的地面要素场的结构：飑线由雷暴单体侧向排列而成，每个单体成熟期都有地面冷丘及水平外流和阵风锋；小冷丘和阵风锋结合起来形成小尺度雷暴高压和阵风锋。
阵风锋又称为飑锋：处在雷暴高压边缘。
具有很强的温度梯度、气压梯度和风速和风向水平切变，它也叫气压涌升线或跳跃线。
○5过境特征：由于飑锋附近是各种气象要素水平梯度很大的地带，因此当飑锋过境时，气象要素将发生急剧的变化。
通常表现为气压涌升、气温急降、风向突变、风速剧增以及强烈降水等。
飑线前低压：飑锋前方一般有中尺度低压。
它的形成可能与飑线前方高层的补偿下沉气流引起的绝热增温有关。
尾流低压：雷暴高压后方的中尺度低压，它的形成与雷暴高压后部的尾流效应相联系。
飑中系统：包括飑线、飑线前低压、雷暴高压以及尾流低压统称为飑中系统。
※飑中系统的全部系统一般只在成熟阶段才同时出现。
不同阶段系统的强度和结构是不同的。
两类比较常见飑线：1、具有前导对流线和尾随层状云区以及具有由前向后和由后向前两支入流的飑线发生在风垂直切变相对小的环境中的飑线飑线的前方有一支由前向后的入流迎着飑锋上升，到高层分裂成向前和向后的两支气流，其后部中层则另有一支由后向前的入流。
中尺度数值模拟报告中尺度数值模拟是一种重要的气象预报手段，可以对天气过程进行较准确的预测和分析，尤其在短期天气预报中具有很高的实用价值。
以下是一份中尺度数值模拟报告的范例。
报告名称：2021年8月21日北京市短期天气预报预报时间：2021年8月20日15时一、天气概况北京市区今天（8月20日）自早晨以来开始阴雨天气，气温明显下降。
预计明天（8月21日）北京市有小到中雨，其中西南部地区部分地方有暴雨，受降雨影响，气温下降较大，最高气温不超过27℃。
二、气象预报1. 降水预报北京市区明天上午有小到中雨，中午时段转为零散小雨。
西南部地区降水较强，局地有暴雨，建议做好防御准备。
预计24小时内，北京市区累计降水量为10-25毫米，局部西南部地区可能达到30-50毫米。
2. 温度预报明天北京市气温将继续下降，最高气温不超过27℃，最低气温为20℃左右。
各区气温预计变化范围为：东城区、西城区、朝阳区、海淀区、石景山区、丰台区、通州区、房山区、顺义区、门头沟区、昌平区、大兴区、平谷区最高气温均在27℃以下。
3. 风力预报明天北京市区气流较强，东部地区有6-7级偏东北大风，其他区域风力为4-5级偏东北风。
三、预警提示根据气象预报，预计明天北京市西南部地区降水较强，局地有暴雨，建议留意山区洪水和滑坡灾害的可能性，及时采取措施，确保人身安全。
四、评估分析此次天气系统来袭，与强冷空气和副高相互作用使得北京市气温下降，降水增多的趋势很明显。
目前各项数据稳定，预报准确度较高。
综合分析，明天北京市仍有较强的降水和大风天气，需要做好防护措施。
五、预报措施依据气象预报，明天初始化观测方案包括增加对西南部地区的降水监测和洪水及滑坡等风险评估，及时调整预警方案，避免因天气带来的自然灾害。
同时，加强监测台站、拓展网络、科学管理，不断提高短期天气预报的准确率和精度。
以下均为老师的重点，按顺序整理1.简述Orlanski分类法对中尺度的分类？2.简述中尺度天气系统的基本特征？（1）空间尺度小，生命期短。
β中尺度系统的水平尺度(L)一般为20-200 km，垂直尺度(H)为10 km左右，因而形态比H／L为10-1-100β中尺度系统的生命史一般在几小时到十几小时，（2）气象要素梯度大。
中尺度天气系统气象要素的梯度很大，气压达1~3 hPa／10 km，温度3℃以上／10 km，中尺度系统如飑线过境时，变温为10℃／15 min左右，变压为6 hPa／15 min 左右。
中尺度天气系统产生的天气现象一般比较激烈（3）非地转平衡和非静力平衡及强的垂直运动。
尺度分析表明，中尺度系统的动量方程中，加速度项与地转偏向力和气压梯度力具有相同的量级，不能满足地转平衡关系对于β中尺度，尤其是γ中尺度运动，准静力平衡假定对所描述的中尺度系统有明显的歪曲，垂直运动速度也明显大于大尺度运动（4）小概率和频谱宽、大振幅事件。
中尺度系统在统计的意义上是小概率的它的空间尺度跨越的范围宽（频谱宽）且中尺度系统影响时的要素变化激烈，表明它是频谱宽的大振幅事件。
3.滞弹性近似、包辛内斯克近似及适用条件？1.包辛内斯克近似：(1)在连续性方程中不考虑密度的个别变化，是完全非弹性，因此是速度无辐散的；(2)与重力相联系的方程中要部分考虑密度的影响，(3)状态方程或热流量方程中要考虑密度变化的影响，(主要是由受热不均匀即温度变化引起的密度变化)，不考虑压力效应对密度变化的影响；它主要适用于浅对流的中尺度运动2.滞弹性近似比包辛内斯克流体在弹性近似方面更接近实际流体，其主要特征是：（1）连续性方程中虽然不考虑密度的个别变化，但保留了平均密度的垂直变化，因而是滞弹性的，或称之为质量无辐散，（2）与重力相联系的方程及状态方程和热流量方程中要同时考虑压缩效应和热膨胀效应引起的密度变化。
因而滞弹性流体主要适用于研究深对流的中尺度运动，它是另一种形式的包辛内斯克近似4.什么是“对称不稳定”?所谓对称不稳定，从物理上看就是大气运动在垂直方向上是对流稳定的和水平方向上是惯性稳定的情况下，作倾斜上升运动时仍然可能发生的一种不稳定大气现象5.简述强风暴发生的天气学必要条件？①位势不稳定层结，并常有逆温层存在；②低层有水汽辐合；③有不稳定的释放的机制；④强的风垂直切变；⑤低空急流；⑥中空干冷空气等6.什么是条件不稳定、对流不稳定？其适用条件各是什么？对流性天气一般发生在条件性不稳定层结的情况下，但有时在上干下湿的条件性稳定层结下，如果有较大的抬升运动，特别是发生整层大气得到抬升时，原先的条件性稳定层结变成不稳定的了，这种不稳定层结称为对流性不稳定，其适用于气层。
1、 控制方程：描述某一物理现象的（材料成形过程中的物理现象）的微分方程（组）称为~（也称场方程）。
它分为解析方法，数值方法。
2、 数值模拟：根据工程问题的基本规律（控制方程和边界条件），利用计算机程序求出满足工程要求的数值解。
分为有限元法，有限差分法，边界元法。
3、 CAE ：计算机辅助工程分析，使用计算机工程分析软件来辅助工程师进行设计后的分析、预测或进行同步工程。
4、 有限差商：通常把一个连续函数f （x ）的增量与自变量的增量的比值。
5、 步长：区域离散化时剖分单元的尺寸。
可分为等步长、变步长；时间步长、空间步长。
6、 有限元分析：利用数学近似的方法对真实物理系统（集合和载荷工况）进行模拟，利用简单而又相互作用的元素（单元），用有限数量的未知量去逼近无限未知量的真实系统。
7、 自由度：（DOFs ）用于描述一个物理场的响应特性。
结构---位移；热---温度；流体---压力；电---电位；磁---磁位。
8、 节点：有限元模型中空间的坐标位置，具有一定自由度和存在相互物理作用。
9、 单元：有限元模型中节点间相互联系的抽象物体，具有节点自由度间相互作用的数值、矩阵描述（称为刚度或系数矩阵）。
10、弹性有限元分析分为位移法、力法、混合法。
11、位移法：选择节点位移作为基本未知量。
应用范围最广。
12、力学：选择节点力作为基本未知量。
13、混合法：取一部分节点力和一部分节点位移作为基本未知力。
14、刚度法：通过建立节点力和节点位移的关系来求解结构问题的一种有限元方法。
15、单元刚度矩阵：假设单元的节点位移为de ，该单元的节点力为fe ，两者的关系可表示为：是一个方阵）称为单元刚度矩阵。
则(e e e e K d K f =。
16、整体刚度矩阵：将单元刚度矩阵扩展到研究区域的所有单元（N 个单元），即称为整体刚度矩阵。
则K K K N e e ∑==117、将位移表示为坐标变量的简单函数，这种函数称为~或位移函数、最常用的位移模式是多项式。
中小尺度数值模拟报告姓名：巴桑扎西学号：20121301002学院：大气科学专业：大气科学12级气科1班一、实习要求：1、wps的参数配置，包括namelist.wps和pltgrid.ncl所画的图。
2、WRF的参数配置，即WRF的namelist.input。
3、台风，飑线每三小时，超级单体每1小时的组合回波及地面降水图。
4、根据预报结果的适量分析。
二、实习内容：（1）、模拟个例的简要介绍：（2）、本次模拟时间段：初始时间：'2005-08-03_00:00:00','2005-08-03_00:00:00', 结束时间：'2005-08-03_12:00:00','2005-08-03_12:00:00',三、实习结果：（1）、WPS部分：在namelist.wps的内容以及ncl脚本生成的模拟区域图&sharewrf_core = 'ARW',max_dom = 2,start_date = '2005-08-03_00:00:00','2005-08-03_00:00:00',end_date = '2005-08-03_12:00:00','2005-08-03_12:00:00',interval_seconds = 21600io_form_geogrid = 2,/&geogridparent_id = 1, 1,parent_grid_ratio = 1, 3,i_parent_start = 1, 10,j_parent_start = 1, 20,e_we = 74, 112,e_sn = 61, 97,geog_data_res = '10m','2m',dx = 45000,dy = 45000,map_proj = 'mercator',ref_lat = 24.00,ref_lon = 130.00,truelat1 = 30.0,truelat2 = 30.0,stand_lon = 130.0,geog_data_path =’/opt/GEOG’/&ungribout_format = 'WPS',prefix = 'FILE',/&metgridfg_name = 'FILE'io_form_metgrid = 2,/（2）、WRFV3部分：在WRFV3目录里的namelist.input的&physics部分&physicsmp_physics = 3, 3, 3, ra_lw_physics = 1, 1, 1, ra_sw_physics = 1, 1, 1,radt = 45, 45, 45, sf_sfclay_physics = 1, 1, 1, sf_surface_physics = 1, 1, 1, bl_pbl_physics = 1, 1, 1, bldt = 0, 0, 0, cu_physics = 1, 1, 0, cudt = 5, 5, 5, isfflx = 1,ifsnow = 0,icloud = 1,surface_input_source = 1,num_soil_layers = 5,sf_urban_physics = 0, 0, 0, /（3）、模拟结果：台风路径三、WRFV3部分。
中小尺度数值模拟第二章 控制大气的基本方程组和地图投影1、数值模式下的平均运动方程组⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧===∙∇+=∙∇+∂∂+⨯Ω-+∇-=RT p S dt dq Q V p dt dT C V t F V g p dtV d v ρρρρρρ,0213333333 2、次网格过程用格点值不能直接描述的过程。
如湍流，积云对流→对流扰动（空间）。
平均值对应网格点值，扰动值对应次网格过程。
3、气象上常用的地图投影有哪些？在数值模拟过程中如何确定使用何种地图投影？（1）极射赤面投影：高纬地区变形小，多用作极地天气图和北半球天气图的底图。
（k=1）（2）兰伯特投影：中纬地区变形小，适用于中纬度地区天气图的底图，如亚欧天气图。
（3）麦卡托投影：低纬地区变形小，适用于作低纬地区天气图的底图。
（k=0）第三章 方程组的离散—差分法4、中小尺度现象的模拟方法（1） 实验室模型：动力相似性（2） 解析解：数学方法求解方程（连续）（3） 数值模式：离散数值求解方法5、数值求解的主要方法及其适用范围有限差分法：利用Taylor 级数展开，使用最广泛。
有限元法：根据能量最小原理，将积分区域划分为有限的、不重叠但互相连接的单元，每个单元选择基函数，用单元基函数的线性组合逼近单元中的真解，整体区域的基函数可以看作由每个单元基函数组成，整个区域的解可以看作由所有单元的近似解，适用于不规则区域。
谱方法：Fourier 变换（谱展开），主变换要用于全球模式，计算精度高，现已经在中尺度模式中应用。
6、什么是差分方程的一致性、收敛性和线性稳定性？一致性（相容性）：即差分方程是否逼近微分方程（当步长→0）收敛性：在一定的定解条件下，差分方程的解U 是否逼近微分方程的解u 。
稳定性：在时间积分过程中，当时间步长趋向于0时，在整个求解区域内，舍入误差U U -保持有界则是稳定的。
即差分解的误差不随时间增长。
7、什么是差分方程的非线性计算不稳定？产生的原因是什么？答：在满足线性稳定性条件下，由于非线性作用而产生的不稳定，称为非线性不稳定。
数值天气预报第十章几种数值模式及模拟试验举例兰州大学大气科学学院中小尺度天气系统常与暴雨、冰雹、雷雨大风等剧烈天气过程联系在一起。
随着探测手段的进步，监测和跟踪能力的提高，对中小尺度天气系统发生、发展机制的探讨及预报方法的研究，近年来取得迅速的进展。
1986——1990 年期间，中国建立了京津冀地区的中尺度天气系统的监测和预报基地。
数值模拟是用试验的手段分析中小尺度天气系统的理论工具。
用其研究中小尺度天气过程，可以避免研究中求解非线性方程组的困难，但却较真实的揭示出影响中小尺度天气过程的物理因子以及演变的细节。
本章内容描述中小尺度天气系统的基本方程组描述中小尺度天气系统的线性方程组的动力学特征模拟中小尺度天气系统应考虑的物理因子中尺度天气系统数值模拟实例3.其他方程的简化热流量方程、水分方程及其他气体和气溶胶方程，一般较少做简化，采用原有方程形式。
只有当源汇项较小时，可将它们略去。
小结：简化的中小尺度系统的方程组，其连续方程根据系统深厚或浅薄分别取(10.20)式或(10.22)式；其垂直运动方程根据系统的水平尺度分别取(10.27)或(10.28)式；水平运动方程采用(10.32)式和(10.33)式；而其他方程一般仍取原有形式，即(10.2)及(10.4)至(10.8)式。
这些方程则构成了描述中小尺度系统的基本方程组。
(10.75)(10.76)引入这两个参数的目的是为了在分析中追踪它们所代表的项；以便很清楚地了解采用流体静力假设或滞弹性假设应被忽略的项。
二、形式解与频率方程上述方程组中含有6个未知数及，有6个方程，则方程组闭合，可以求解。
设各未知数有下列形式的解：11,0,λ⎧=⎨⎩21,0,λ⎧=⎨⎩表示流体是非静力的表示流体是静力平衡的表示流体是可压缩的表示流体是滞弹性的u v w p ρ′′′′′、、、、θ′(10.77)将上述形式解代入方程组(10.69)～(10.74)中，可得含有6个未知数的线性齐次代数方程：%()()()() ()() ()()()()%()(),,,,,,,,,,,,x z x z x z x z x z x z i k x k z t x zi k x k z t x z i k x k z t x z i k x k z t x z i k x k z t x z i k x k z t x z u u k k e v v k k ew w k k e p p k k e k k e k k eωωωωωωωωωωρρωθθω+−+−+−+−+−+−⎧′=⎪⎪′=⎪′⎪=⎪⎨′=⎪⎪′=⎪⎪′=⎪⎩%力假定或滞弹性流体假定去除。
中小尺度数值模拟复习整理第一章绪论第二章控制大气运动的基本方程组和地图投影1、中尺度动力学控制方程（1）为什么需要将方程组重新写成平均意义上的方程形式答：①气象实际观测资料是在一定意义下的时空平均值②数值模式由于时空离散化，变量值均为间隔为x, y, z, t点上的值，即网格点上值，也可理解为网格区域的平均值2、数值模拟：给定初值和边值条件，数值求解控制大气运动方程组。
其中，主要用格点法（或谱方法）把连续的时间空间变为离散的时间空间，用变量在离散点上（网格点上）的值来代替连续变量场，用差商代替微商，将微分方程变为差分方程（谱方程）3、次网格过程：用格点值不能直接描述的过程4、湍流与对流有区别：湍流：发生在边界层，耗散动量、输送热量和水汽对流：发生在对流层，对流凝结释放潜热加热大气，输送热量和水汽到中高层，也具有耗散作用5、地图投影：把地球投影到一个投影曲面上，然后把投影曲面展开成平面，这就是投影平面，最后将投影平面缩小一定的倍数，就成为地图。
6、比例尺：投影平面上某一线段长度和地球上相应长度的比值。
7、为什么要进行地图投影？答：①进行有限区域数值模拟时，常常使用平面直角坐标系，因此常需将球面投影于平面上，重新划分网格。
②应用于天气图分析（天气图底图是用来填写各地气象站观测记录的特种空白地图）和和数值预报（数值模拟初值来自客观分析，结果也表示在天气图上）8、投影面：把地球表面的地形，经纬线等投影到一个有简单图形的面上，如：平面，圆锥面，圆柱面（后二者可展开成平面），这些面称为投影面。
9、投影变形：地球表面经投影的方法描绘在平面上，地球上地理区域的距离、方向、面积、形状等特征会有变形。
10、等距投影：距离不变的投影11、正形投影：两条交线交角保持不变、且各方向放大(或缩小倍数)相等的投影（即地图上某一个点各个方向的比例尺保持不变）12、透视：按几何透视来分，如一平面与地球北极相切，球心有一光源，由球心发出的光线将经纬线投影到平面。
一、实习背景中尺度数值模拟是气象、海洋、大气科学等领域中研究自然现象的重要手段。
通过对中尺度大气和海洋过程进行数值模拟，可以揭示其内在规律，为天气预报、气候预测和灾害预警提供科学依据。
为了提高自己的专业技能，我于20xx年暑假期间参加了中尺度数值模拟实习。
二、实习目的1. 掌握中尺度数值模拟的基本原理和方法；2. 学习数值模拟软件的使用，提高数值模拟技能；3. 通过实际案例分析，加深对中尺度现象的理解；4. 培养团队协作精神和沟通能力。
三、实习内容1. 中尺度数值模拟基本原理实习期间，我学习了中尺度数值模拟的基本原理，包括动力方程、热力学方程、微物理过程等。
通过对这些基本原理的掌握，我能够更好地理解中尺度数值模拟的物理过程。
2. 数值模拟软件使用实习过程中，我学习了WRF（Weather Research and Forecasting）数值模拟软件的使用。
WRF是一种广泛应用于中尺度数值模拟的软件，具有强大的功能和灵活的配置。
通过实际操作，我掌握了WRF的基本操作流程，包括模式设置、数据准备、模拟运行和结果分析等。
3. 实际案例分析为了加深对中尺度现象的理解，实习期间，我选取了以下几个实际案例进行分析：（1）2005年山东半岛大雪过程数值模拟与中尺度特征分析通过WRF模型模拟2005年山东半岛大雪过程，分析了冷空气爆发、地形影响和海陆风环流等因素对大雪过程的影响。
模拟结果表明，冷空气爆发是导致大雪的主要原因，而地形和海陆风环流也起到了关键作用。
（2）基于LAPS雷达资料变分分析技术的暴雨数值模拟及其中尺度结构特征研究利用LAPS雷达资料变分分析技术，对暴雨数值模拟进行了研究。
分析了雷达资料同化对暴雨预报能力的影响，以及暴雨中尺度结构特征。
结果表明，LAPS雷达资料变分分析技术能够有效提高暴雨预报的准确率。
（3）南海北部以及吕宋海峡次级中尺度动力过程数值模拟研究通过使用单向嵌套方案AGRIF和原始方程海洋模式ROMS，对南海北部陆架海区的次级中尺度过程进行了模拟。
中小尺度数值模拟
第二章 控制大气的基本方程组和地图投影
1、数值模式下的平均运动方程组
⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧===∙∇+=∙∇+∂∂+⨯Ω-+∇-=RT p S dt dq Q V p dt dT C V t F V g p dt
V d v ρρ
ρρρρ,0213
333333 2、次网格过程
用格点值不能直接描述的过程。
如湍流，积云对流→对流扰动（空间）。
平均值对应网格点值，扰动值对应次网格过程。
3、气象上常用的地图投影有哪些？在数值模拟过程中如何确定使用何种地图投影？
（1）极射赤面投影：高纬地区变形小，多用作极地天气图和北半球天气图的底图。
（k=1）
（2）兰伯特投影：中纬地区变形小，适用于中纬度地区天气图的底图，如亚欧天气图。
（3）麦卡托投影：低纬地区变形小，适用于作低纬地区天气图的底图。
（k=0）
第三章 方程组的离散—差分法
4、中小尺度现象的模拟方法
（1） 实验室模型：动力相似性
（2） 解析解：数学方法求解方程（连续）
（3） 数值模式：离散数值求解方法
5、数值求解的主要方法及其适用范围
有限差分法：利用Taylor 级数展开，使用最广泛。
有限元法：根据能量最小原理，将积分区域划分为有限的、不重叠但互相连接的单元，每个单元选择基函数，用单元基函数的线性组合逼近单元中的真解，整体区域的基函数可以看作由每个单元基函数组成，整个区域的解可以看作由所有单元的近似解，适用于不规则区域。
谱方法：Fourier 变换（谱展开），主变换要用于全球模式，计算精度高，现已经在中尺度模式中应用。
6、什么是差分方程的一致性、收敛性和线性稳定性？
一致性（相容性）：即差分方程是否逼近微分方程（当步长→0）
收敛性：在一定的定解条件下，差分方程的解U 是否逼近微分方程的解u 。
稳定性：在时间积分过程中，当时间步长趋向于0时，在整个求解区域内，舍入误差U U -保持有界则是稳定的。
即差分解的误差不随时间增长。
7、什么是差分方程的非线性计算不稳定？产生的原因是什么？
答：在满足线性稳定性条件下，由于非线性作用而产生的不稳定，称为非线性不稳定。
非线性计算不稳定主要表现为短波能量的虚假增长。
产生的原因：①可能是混淆误差 ②可分辨波能量之间的相互转换中的误差 ③还没有完全弄清楚
8、混淆误差
有限网格系统能分辨的波的最短波长为x ∆2，对于非线性作用产生的波长小于x ∆2的波动，网格系统不能正确地分辨，而把它错误地表示成为某种波长大于x ∆2的波，从而造成了这种波的误差——混淆误差。
9、如何消除差分方程的非线性不稳定？
（1）空间、时间平滑，滤去短波分量，滤去波长< x ∆4的波可消去alias 解。
（2）加水平扩散项A r 2
∇，A 为物理量，r 为扩散系数，可根据作用大小选择。
（3）构造具有隐式平滑式某种选择性衰减作用的差分格式。
以上方法优点：实用有效；缺点：克服短波同时，也会削弱长波，影响预报效果。
（4）构造差分格式：尽量保持连续系统的一些积分约束的差分格式。
能量守恒差分格式得总能量约束。
涡度拟能守恒差分格式得能谱约束，抑制能量向短波转移。
10、为何能量只会在短波上发生虚假的自反馈？
第四章 离散网格、资料处理与初始化
11、数值模拟中模拟区域和格距大小的设置要求有哪些
区域设置：（1）中小尺度过程引起的变量的有意义变化必须被包括在模式空间中，即模式范围要足够大，从而边界上基本不受中小尺度过程影响。
（2）外强迫的尺度和范围（如地形、加热等）
（3）对强迫的响应尺度
格距设置：（1）中小尺度系统的特征必须能足够被分辨，即中小尺度系统至少要占4个格距。
（2）设置的格距要适合计算方案的稳定性要求。
（3）计算代价最小，必须考虑计算资源的占用问题。
12、水平网格结构
①等距网格。
y x ∆∆,都是常数，优点：简单
②拉伸网格。
y x ∆∆,不是常数，视需要调整。
优点：节省计算资源，同样的网格数可以有更大的模拟区域。
③嵌套网格。
分粗网格和细网格，优点：节省计算资源，通过多层嵌套，可以实现高分辨率的模拟。
④移动网格。
可跟随着中尺度系统的发展来移动，优点：简单，计算代价小，适用移动范围较大的系统的模拟。
13、何为交错网格？为什么垂直方向要采用拉伸网格？
交错网格：将不同的气象要素有规律的写在不同网格上。
大部分中尺度天气现象的强迫主要位于对流层低层大气，绝大多数的数值模式在垂直方向上都采用拉伸网格结构。
14、为什么由客观分析得到的格点值并不能直接作为初始场
观测资料或分析资料的误差会破坏风场和气压场之间的不平衡。
初始资料和模式之间可能不协调，如果不进行初始化处理，在模式积分的一段容易出现虚假的高频波。
15、初始化方法
客观分析方法：有限元方法，逐步订正法，九点拉格朗日插值，最小二乘法等。
另外，由于非定时探空资料日益增多，客观分析还包括了卫星资料、雷达资料等多种数据的使用。
16、边界设置
海绵边界条件，即在边界附近设一缓冲带，一般三至五圈，设为模式中变量，有模式计算输出的局地变化值为经海绵边界处理后实际采用的值。
海绵边界条件的实质是在边界上对固定的外边界和内部的预报场之间起一种平滑协调的作用，它可以对向外传播的惯性重力波较好地吸收，减小边界上波的反射。
对于深对流系统的模拟，一般取上边界在对流层顶或平流层下层。
对于浅对流系统的模拟，上边界可以取在对流层中层，模式层顶的选取主要依赖于实际模拟系统的系统性质。
17、大气边界层的作用
①厚度：1~2km （分三层）
②主要过程是湍流摩擦（分子摩擦较小，可略）
③边界层内动量、热量和水汽的湍流垂直通量最大
④是大气重要的动量汇，热量和水汽的源。
因此，大气通过边界层的湍流输送而和下垫面发生热力和动力作用。
第五章数值模拟参数化
18、参数化——用大尺度变量表征次网格或小尺度作用总体效应的方法。
积云对流参数化：将大尺度模式不能显式分辨的对流凝结和对流引起的热量、水分和动量的输送与模式的预报变量联系起来。
实际做法是用大尺度网格点上的变量值来表示次网格尺度积云的凝结加热以及垂直输送效应的作用。
积云参数化方案（CPSs）通常要求在网格尺度未饱和时产生能垂直输送热量、水汽和其他物理量的次网格尺度隐式云。
闭合假设需要确定隐式云和网格尺度变量之间的关系。
一般来说，<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设置垂直层数。
