共查询到18条相似文献，搜索用时 875 毫秒
1.
本文用差分模拟绝对涡度守恒定律,构成正压涡度方程的两类差分格式:积分守恒型格式和相当于准拉格朗日型格式。前一类型格式的显式部分是条件稳定的,隐式部分是无条件稳定的。后者至为无条件稳定的显式格式,可以依据实际需要和可能,来放大无条件稳定格式的时间步长。 本文还根据格式的计算要求,讨论了“风”场的计算和平滑问题,以及差分方程的边、初条件(计算边、初条件)的给法。并通过理想场计算和实例预报,检验了几种显式格式的计算效果。 相似文献
2.
本文以带耗散项的非线性方程——Burger方程为模型,用截谱形式的差分解代入相应的差分方程,分析了数值天气预报中常用的一类三时间层差分格式的计算稳定性问题,特别讨论了计算初值、时间步长和耗散因子对计算稳定性的影响.利用定性分析和数值试验的方法证实这类三时间层格式有时容易出现非线性计算不稳定,尤其是两步计算初值反号和耗散系数过小的时候是如此.但是耗散系数也不能取得过大,否则也会出现计算不稳定. 相似文献
3.
天气数值预报工作中所用的初始方程显式差分格式,时间步长取得很短,实际应用很不方便。比如,空间步长取300公里,时间步长只能用10分钟。本文提出一种新型显式差分格式,其计算稳定性不受时间步长(τ)和空间步长(d)大小的限制,由|λ_e|τ/d(e=1,2,3)的整数部分m_e=[|λ_e|τ/d],控制执行“浮动”计算,保持计算稳定。这种差分格式曾用104电子计算机作检验计算,空间步长用25公里,时间步长取20分钟,计算过程稳定,结果精度较高,可与文献[1]中的结果相比较。这表明,初始方程的差分计算,适当延长时间步长是有可能的。 相似文献
4.
在边界层动力学中，涡动粘性系数是影响边界层风场结构的一个重要参数。本文利用边界层动力学中的Ekman动量近似理论，给出了涡动粘性系数随高度缓变条件下的Ekman动量近似边界层模式解，着重讨论了边界层的风场结构、水平散度、垂直涡度以及边界层顶部的垂直速度。结果分析表明：与常值涡动粘性系数情况相比，在边界层低层随高度增加的涡动粘性系数可以导致低层边界层风速随高度迅速增加，即风速垂直切变增加，同时风速矢与地转风之间的夹角减小。惯性项作用可以导致上述作用在气旋性区域减小、而在反气旋性区域增大。随高度增加的涡动粘性系数导致水平散度绝对值、垂直涡度绝对值以及边界层顶部的垂直速度绝对值在气旋性区域减小，而在反气性旋区域增大。涡动粘性系数与惯性之间的非线性相互作用是边界层动力学中重要过程。 相似文献
5.
采用与实测较接近的二次函数来表达Ｅｋｍａｎ层中的湍流粘性系数Ｋ，在圆形气压场条件下，求得了山地上空边界层中的风速，进而求得散度、涡度和垂直速度等场变量随高度的分布。并作图分析了这些场变量的一些动力学特征。改进了以往在求解析解时，略去运动方程中湍流粘性力项中的关于高度的一阶导数项，以及取山坡面上风速为零作下边界条件等欠合理欠精确的做法。所求得的风速、散度、涡度和垂直速度均用简单的初等函数表示出来，有助于边界层参数化和深化对边界层动力学的认识。 相似文献
6.
针对■ζ/■_t+u■ζ/■_x=0方程,本文首先讨论了在一般情形下所谓“先向前差、后中央差”差分格式的稳定性,并给出了稳定的条件。其结果和一般只用中央差分的情形有所不同。利用三种差分格式进行比较,发现目前广泛应用的用向前差计算第一个时间步长,以后全用中央差的格式是最差的;当λ(=uΔt/Δs)→1—0时,计算并不稳定。如果用另两种格式中任一种代替向前差,则可以避免这个困难。其中又以“逆2—3格式”较好。 相似文献
7.
运用判定非线性发展方程差分格式计算稳定性的Hirt启发性分析方法,对一类非线性Schrdinger方程差分格式的计算稳定性进行分析,得到了保证差分格式计算稳定的必要条件.数值试验结果进一步表明,得到的稳定性判据不仅是保证差分格式计算稳定的必要条件,而且在实际中也是非常有效的. 相似文献
8.
变步长显式完全平方守恒差分格式 总被引：8，自引：0，他引：8
综合隐式完全平方守恒差分格式和显式瞬时平方守恒差分格式的优点，针对一类非线性发展方程构造了一种通过自动调节时间步长来保持平方守恒性的显式差分格式。它基于加小耗散的思想，但又与小耗散法有所不同。本文取的耗散项不是一般的人工耗散，而是取能够弥补由于时间离散所产生的截断误差的所谓（时间）协调耗散。因此，该格式具有较高的时间精度。在数值试验中，该类格式可取得满意的效果。 相似文献
9.
将包含水汽潜热和摩擦因子的平衡模式ω方程适当简化,形式上与准地转方程相似,在简单谐波分布的假定下,由此求得各标准气压层上垂直运动和涡度倾向的近似解析解,它是包含所有各气压层上涡度源变和温度源变(或称“广义平流”)以及地形强迫作用的一组线性多项式。各项的系数分别表现出该项源变或地形强迫作用对某一气压层上垂直运动或涡度倾向的影响。针对几种典型情况,计算出这些影响系数,这有助于定性地诊断天气尺度和次天气尺度系统短期演变的物理过程,适用于副热带情况。 相似文献
10.
进一步讨论了有关非线性不稳定的一些问题,其主要内容有： １．考察了有代表性的三类发展方程,指出其对应的差分格式是否出现非线性计算不稳定,与原微分方程解的性质密切相关。 ２．进一步讨论了带周期边条件的守恒型差分格式的非线性计算稳定性问题,总结了克服非线性不稳定的有效措施。 ３．以非线性平流方程为例,着重分析了带非周期边条件的非守恒差分格式的非线性计算稳定性问题,给出了判别其计算稳定性的“综合分析判别法”。 相似文献
11.
在数值预报和数值模拟中,描述空间微分项的最主要的方法是有限差分法,但使用差分方法会引入截断误差.伍荣生1979年指出,通过在原物理场的基础上构造一个新的物理场,替代原物理场进行差分计算,可以达到减小误差的目的.该文是伍荣生1979年工作的继续,目的在于解释伍荣生1979年所构造的差分格式并得到更为一般化的差分格式.文中给出新的差分格式结合了经典有限差分方法的快速计算和谱方法的高精度的优点.如果在一个给定的网格上对气象要素场进行离散傅里叶级数展开,则基函数(正弦或余弦)的频谱是事免已知的.作者将伍荣生1979年构造物理场的方法视为对物理场的一次平滑,探讨了获取二次平滑场、多次平滑的一般化方法.获取平滑场的基奉原理是使得在固定频谱上的差分逼近程度达到最优.通过对频谱上的累计误差的下降速度分析表明,平滑次数的上限为3次.数值分析的结果表明,二次平滑的最大误差是未作任何平滑的最大误差的0.04倍,在使用相同计算代价的情况下,二次平滑的最大误差是经典的差分格式的0.3倍.平流试验的结果也表明,新的差分格式即一次平滑、二次平滑方案的结果远远优于经典的差分格式.新的差分格式意义在于,在不加密网格的情况下提供了一条提高数值计算精度的途径. 相似文献
12.
AN IMPROVED SEMI-IMPLICIT TIME DIFFERENCE SCHEME OF SPECTRAL MODEL AND NUMERICAL APPLICATIONS 总被引：1，自引：0，他引：1
In fact,the popular semi-implicit time difference scheme of spectral model still includes someimportant linear terms using time explicit difference scheme,and the major terms are directlyrelated to fast internal-and external-gravity waves in the atmospheric forecasting equation.Additionally,due to using time difference on two terms at different time.the popular schemeartificially introduces unbalance between pressure gradient force and Coriolis force terms whilenumerically computing their small difference between large quantities.According to thecomputational stability analysis conducted to the linear term time difference scheme in simpleharmonic motion equation,one improved semi-implicit time difference scheme is also designed inour study.By adopting a kind of revised time-explicit-difference scheme to these linear terms thatstill included in spectral model governing equations,the defect of spectral model which only partlyusing semi-implicit integrating scheme can be overcome effectively.Moreover,besides all spectralcoefficients of prognostic equations,especially of Helmholtz divergence equation,can be workedout without any numerical iteration,the time-step (computation stability) can also be enlarged(enhanced) by properly introducing an adjustable coefficient. 相似文献
13.
AN IMPROVED SEMI-IMPLICIT TIME DIFFERENCE SCHEME OF SPECTRAL MODEL AND NUMERICAL APPLICATIONS* 
下载免费PDF全文
下载免费PDF全文 In fact,the popular semi-implicit time difference scheme of spectral model still includes some important linear terms using time explicit difference scheme,and the major terms are directly related to fast internal-and external-gravity waves in the atmospheric forecasting equation.Additionally,due to using time difference on two terms at different time.the popular schemear tificially introduces unbalance between pressure gradient force and Coriolis force terms while numerically computing their small difference between large quantities.According to the computational stability analysis conducted to the linear term time difference scheme in simple harmonic motion equation,one improved semi-implicit time difference scheme is also designed in our study.By adopting a kind of revised time-explicit-difference scheme to these linear terms that still included in spectral model governing equations,the defect of spectral model which only partlyusing semi-implicit integrating scheme can be overcome effectively.Moreover,besides all spectral coefficients of prognostic equations,especially of Helmholtz divergence equation,can be worked out without any numerical iteration,the time-step (computation stability) can also be enlarged (enhanced) by properly introducing an adjustable coefficient. 相似文献
14.
15.
利用WRF中尺度数值模式,NCEP/NCAR分析资料,多普勒雷达观测资料等,对2016年7月25日一次东北冷涡下的飑线过程进行数值模拟,研究了飑线形成和维持与水平涡度的关系及飑线过程中中尺度对流涡旋(MCV)的形成机制,分析发现,高低层水平涡度逆时针旋转对本次飑线的形成和维持有很好的指示意义。(1)飑线发生前,高层渤海湾西侧出现水平涡度的逆时针旋转中心,并有较强的辐散配合,低层水平涡度为逆时针弯曲,为飑线产生提供了有利的上升运动条件。随后高层多个对流单体的水平涡度气旋式涡旋合并形成较大范围的气旋式涡旋结构,触发低层的上升运动,同时低层对流区前部形成一致的气旋式弯曲使得对流单体组织成带状结构,形成飑线。(2)飑线成熟时期高层水平涡度表现为统一大范围气旋式涡旋结构,低层则呈现典型的S型弯曲结构,水平涡度x方向的分量沿对流带从南至北表现为正负正,y方向的分量始终为正,并由对流带的中心向两侧减小,显示出水平涡度矢量旋转的方向对飑线影响的重要性。(3)由垂直涡度方程的分析得出,在飑线发展中期,MCV形成前,雷达反射率回波在500 hPa左右表现出明显的旋转,此时主要与500 hPa以上强的正涡度水平平流项及中层倾侧项和水平散度项有关,之后,在这几项的作用下使得中层风场产生气旋式旋转,形成MCV。 相似文献
16.
完全非内插半拉格朗日格式在一维
Burgers方程及二维浅水波方程上的应用 总被引：3，自引：2，他引：1
在作者过去提出的完全非内插半拉格朗日格式的基础上,针对半拉格朗日格式由于内插带来预报场人为的光滑性问题,进一步发展了这种计算格式,证明了此格式的计算稳定性。为检验这种新的计算格式的性能,在一维和二维问题上进行了应用。在一维问题中采用了一维无粘Burgers方程（方程中有突变点）;二维问题采用了浅水波方程,同时将这些计算结果与Ritchie方案及欧拉方案或一般半拉格朗日内插方案的计算结果进行了比较,发现新格式消除了内插和预报场的人为光滑,并且计算精度有一定程度的提高,这为以后将此格式推广到全球谱模式打下了基础。 相似文献
17.
湿位涡和倾斜涡度发展 总被引：266，自引：22，他引：266
从完整的原始方程出发，在导出精确形式的湿位涡方程的基础上，证得绝热无摩擦的饱和湿空气具有湿位涡守恒的特性。并由此去研究湿斜压过程中涡旋垂直涡度的发展。结果表明，在湿等熵坐标中，涡旋的发展与对流稳定度的减少，等熵面上的辐合和潜热的释放有关。由于等熵位涡分析的应用受等熵面倾斜的限制，又进而发展了Ｚ坐标及Ｐ坐标中的倾斜涡度发展理论。指出无论是湿对称不稳定或对流不稳定大气，还是湿对称稳定或对流稳定大气，除对流稳定度的影响外，风的垂直切变的增加或水平湿斜压的增加均能因湿等熵面的倾斜而引起垂直涡度的增长。湿等熵面的倾斜越大，这种由干湿斜压性加强所引起的涡旋发展更激烈。在梅雨锋附近及其南侧暖湿区的北端，湿等熵面十分陡立，是涡旋发展及暴雨发生的重要地区。对１９９１年６月１２—１５日江淮流域暴雨过程的湿位涡分析表明，湿位涡分析，尤其是等压面上湿位涡量ＭＰＶ１和ＭＰＶ２的分析不仅在中高纬有效，在低纬度及低对流层也十分有效，是暴雨诊断和预报的有力工具。 相似文献
18.
The complete form of the vertical vorticity tendency equation (the complete-form vorticity equation) is derived from the Ertel potential vorticity equation to contain thermodynamic factors. In this study, a new complete-form vorticity equation, which has the same form as the original complete-form vorticity equation, is deduced from the absolute vorticity vector equation combined with the continuity equation and the expression of three-dimensional (3D) entropy gradient. By comparing the complete-form vorticity equation with the classical vertical vorticity equation, it is found that regardless of whether or not the isentropic surface is tilting, the two vorticity equations are in essence the same. The “baroclinic term” of the complete-form vorticity equation is exactly equal to the solenoidal term of the classical one, and there is a significant amount of cancellation between the two baroclinic items (the “slantwise term” and the horizontal vorticity change term) in the complete-form vorticity equation. In operational weather analysis, the tilt of the isentropic surface can be diagnosed according to the density of the isotherm on the upper-level isobaric map. For synoptic-scale motion, the vertical vorticity produced by the tilt of the isentropic surface is due to the contribution of atmospheric baroclinicity, which is measured by the solenoid. The 3D solenoid is parallel to the isentropic surface, so the more tilted the isentropic surface, the bigger the projection of the 3D solenoid in the vertical direction. The baroclinic contribution can be interpreted based on the PV thinking theory, but the relationship between the vorticity field and the potential vorticity field is not immediate. 相似文献