本文選題：非靜力模式 + 譜方法　； 參考：《國防科學(xué)技術(shù)大學(xué)》2015年博士論文

【摘要】：數(shù)值天氣預(yù)報(bào)在過去幾十年的發(fā)展中一直采用一個(gè)重要的近似:靜力近似。隨著計(jì)算機(jī)集群計(jì)算能力不斷提升,數(shù)值預(yù)報(bào)模式的分辨率也越來越高,能夠?yàn)橛脩籼峁└泳_的預(yù)報(bào)服務(wù)。一般而言,當(dāng)水平分辨率優(yōu)于10km時(shí),靜力近似已經(jīng)不適用了。因此開發(fā)非靜力模式是未來數(shù)值預(yù)報(bào)系統(tǒng)發(fā)展的迫切需求。相比靜力模式,一方面非靜力模式要增加預(yù)報(bào)變量,方程的形式更加復(fù)雜;另一方面非靜力模式水平網(wǎng)格尺度更細(xì),為了能夠準(zhǔn)確地顯示高分辨率下小尺度地形引起的局部氣象要素,需要很多特定的數(shù)值處理技術(shù)。本文針對現(xiàn)有非靜力模式缺乏高精度垂直離散格式、模式誤差增長估計(jì)手段不足等問題,在已有的數(shù)值方法基礎(chǔ)上,對全球非靜力譜模式絕熱動(dòng)力框架若干關(guān)鍵問題進(jìn)行深入研究,主要工作包括:1.實(shí)現(xiàn)有限差分二維X-Z平面譜模式。二維X-Z平面模式是數(shù)值預(yù)報(bào)模式的研究模式,其方程組和數(shù)值離散要點(diǎn)同三維有限區(qū)域預(yù)報(bào)模式保持一致。對垂直有限差分離散格式進(jìn)行了精巧設(shè)計(jì),使得其滿足非靜力模式的垂直約束。模式采用半隱式半拉格朗日時(shí)間平流離散格式,同時(shí)也對側(cè)邊界、上下邊界條件分別進(jìn)行了處理以保證計(jì)算結(jié)果的正確。在二維模式下,測試了四種山峰波測試樣例以及重力波樣例,驗(yàn)證了模式采用的數(shù)值離散格式的準(zhǔn)確性和穩(wěn)定性。2.實(shí)現(xiàn)垂直混合有限元有限差分二維X-Z平面模式。基于質(zhì)量坐標(biāo)的非靜力模式垂直離散存在一系列的約束條件,應(yīng)用有限元垂直離散被證明是難以實(shí)現(xiàn)的。本文提出了一個(gè)高精度混合離散格式,對半隱式格式的線性部分進(jìn)行有限差分離散,對非線性部分進(jìn)行有限元離散。同時(shí)為提高有限差分離散的精度,增加線性部分計(jì)算的分層數(shù)量。有限元采用四階B-spline樣條基函數(shù)。為方便處理垂直邊界,先把垂直場投射到三階樣條插值函數(shù)空間形成多項(xiàng)式分段函數(shù),在處理完邊界后再投射到B-spline空間。如此避免了B-spline樣條基函數(shù)直接處理邊界,使得有限元計(jì)算更加簡潔高效。對二維有限元模式也進(jìn)行了多種測試樣例,并同有限差分進(jìn)行了定性和定量的對比,結(jié)果表明混合有限元格式在精度上要優(yōu)于有限差分格式。3.實(shí)現(xiàn)三維全球非靜力譜模式。三維球面模式和二維模式的動(dòng)力框架基本一致,但采用的譜方法基函數(shù)有區(qū)別,球面譜模式采用球面諧波函數(shù)為基函數(shù)。球面模式水平網(wǎng)格采用精簡高斯網(wǎng)格,模式方程的曲率項(xiàng)通過半拉格朗日插值時(shí)方向轉(zhuǎn)換矩陣引入�？剖狭�(xiàng)采用的是隱式處理形式。對三維模式進(jìn)行了一系列的測試,包括穩(wěn)定狀態(tài)測試、斜壓不穩(wěn)定波測試、Rossby-Haurwitz波測試、山峰波測試等。測試結(jié)果驗(yàn)證了三維模式的準(zhǔn)確性和穩(wěn)定性。4.研究混沌模型的初始誤差增長特性,發(fā)現(xiàn)混沌模型的平均相對誤差飽和值同初始誤差值之間存在一個(gè)簡單的對數(shù)線性關(guān)系:二者的對數(shù)和為常數(shù)。利用該關(guān)系,提出平均絕對誤差的概念,混沌系統(tǒng)的平均絕對誤差飽和值為與初始誤差無關(guān)的常數(shù)。利用這一特性,給出一個(gè)定量計(jì)算可預(yù)報(bào)期限的數(shù)學(xué)估計(jì)模型。研究了混沌系統(tǒng)模型誤差增長特性,發(fā)現(xiàn)模型平均絕對誤差增長在短期內(nèi)呈指數(shù)增長,一段時(shí)間后達(dá)到飽和�？衫媚Ｐ驼`差增長指數(shù)大小和達(dá)到飽和值的時(shí)間作為衡量模型優(yōu)劣的評價(jià)手段。
[Abstract]:The numerical weather forecast has been using an important approximation in the development of the past several decades: the static approximation. With the increasing computing power of the computer cluster, the resolution of the numerical prediction model is also getting higher and higher, which can provide more accurate prediction service for the users. Generally, when the horizontal resolution is better than 10km, the static approximation has been approximated. It is not applicable. Therefore, the development of non static model is an urgent need for the development of the future numerical forecast system. Compared with the static model, the non hydrostatic model should increase the prediction variable and the form of the equation more complex. On the other hand, the scale of the non static model is finer, in order to accurately display the small scale terrain under high resolution. The local meteorological elements need a lot of specific numerical processing techniques. In this paper, in view of the lack of high precision vertical discrete schemes and the lack of model error growth estimation methods, the main problems of the key problems of the global non static spectral model thermal dynamic force framework are studied in this paper. The work includes: 1. to realize the finite difference two-dimensional X-Z plane spectrum pattern. The two-dimensional X-Z plane model is the model of the numerical prediction model. The equations and the numerical discrete points are consistent with the three-dimensional finite area prediction model. The vertical finite difference dispersion scheme is designed so that it satisfies the vertical constraint of the non static model. The semi implicit semi Lagrange time advection scheme is adopted in the formula, while the side boundary and the upper and lower boundary conditions are treated respectively to ensure the correctness of the calculation results. In the two-dimensional model, four peak wave samples and gravity wave examples are tested. The accuracy and stability of the numerical discrete scheme used in the model are verified.2. real. The vertical mixed finite difference two-dimensional X-Z plane model is presented. There is a series of constraints on the vertical discrete model of the non static mode based on the mass coordinates. The application of the finite element vertical dispersion is proved to be difficult to realize. A high precision mixed discrete scheme is proposed in this paper, and the finite difference dispersion is carried out for the linear part of the semi implicit scheme. In order to improve the precision of the finite difference separation and increase the accuracy of the finite difference separation dispersion and increase the number of layers calculated by the linear part, the finite element uses the four order B-spline spline basis function. In order to handle the vertical boundary, the vertical field is first projected to the three order spline interpolation function space to form a polynomial piecewise function, after the processing of the boundary. Then it is projected into the B-spline space, so that the B-spline spline basis function is avoided directly and the finite element calculation is more concise and efficient. A variety of test examples are also carried out for the two-dimensional finite element model, and the qualitative and quantitative comparison with the finite difference is made. The results show that the hybrid finite element scheme is better than the finite difference. The three-dimensional global non static spectral model is realized by format.3.. The three-dimensional spherical model is basically the same as the dynamic frame of the two-dimensional model, but the basis function of the spectral method is different. The spherical spectral model uses the spherical harmonic function as the base function. The spherical pattern horizontal grid uses the simplification of the Gauss grid, the curvature of the mode equation is interpolated by half Lagrange interpolation The time direction transformation matrix is introduced. The Coriolis force term is used in implicit treatment. A series of tests are carried out on the three-dimensional model, including stable state test, baroclinic unstable wave test, Rossby-Haurwitz wave test, peak wave test and so on. The test results verify the accuracy and stability of the three-dimensional model and the initial error of the chaotic model of the.4. model. It is found that there is a simple logarithmic linear relationship between the average relative error saturation value of the chaotic model and the initial error value: the logarithm and the constant of the two is a constant. Using this relation, the concept of the mean absolute error is proposed. The mean absolute error saturation value of the chaotic system is a constant independent of the initial error. A mathematical estimation model is given for the quantitative calculation of the predictable term. The error growth characteristic of the chaotic system model is studied. It is found that the average absolute error growth of the model increases exponentially in the short term and reaches the saturation after a period of time. The time of the model error growth index to the saturation value can be used as the evaluation of the model. Price means.
【學(xué)位授予單位】：國防科學(xué)技術(shù)大學(xué)
【學(xué)位級別】：博士
【學(xué)位授予年份】：2015
【分類號】：P456.7;O241.8

【相似文獻(xiàn)】

相關(guān)期刊論文 前10條

1 張朝林,鄭慶林;一個(gè)改進(jìn)的譜模式非線性項(xiàng)的計(jì)算方法及其應(yīng)用[J];科學(xué)通報(bào);2002年19期

2 曹鴻興,蔣維東,威廉·哥道夫斯基;自憶譜模式及其初步應(yīng)用[J];大氣科學(xué);1998年01期

3 呂克利，王柏強(qiáng);大氣平衡態(tài)研究中高截?cái)嗟妥V模式的問題[J];熱帶氣象學(xué)報(bào);1996年01期

4 曹鴻興,谷湘潛;自憶譜模式對中期環(huán)流預(yù)報(bào)的改進(jìn)[J];自然科學(xué)進(jìn)展;2001年03期

5 張道民,紀(jì)立人,吳萬里;全球和半球譜模式的比較[J];大氣科學(xué);1990年03期

6 M·Kanamitsu;毛賢敏;王雪晶;陳艷秋;;日本氣象廳業(yè)務(wù)譜模式評述[J];遼寧氣象;1989年02期

7 谷湘潛;一個(gè)基于大氣自憶原理的譜模式[J];科學(xué)通報(bào);1998年09期

8 曹杰;尤亞磊;;低階譜模式的一種新方案及其應(yīng)用[J];應(yīng)用數(shù)學(xué)和力學(xué);2006年04期

9 鄭峰;張小偉;苗長明;LIQIANG SUN;;RSM區(qū)域譜模式特征分析[J];安徽農(nóng)業(yè)科學(xué);2009年12期

10 鄭慶林,宋青麗,蔣平;一個(gè)改進(jìn)的T42L10譜模式及其在月預(yù)報(bào)中的應(yīng)用[J];應(yīng)用氣象學(xué)報(bào);1993年S1期

相關(guān)會議論文 前3條

1 張朝林;鄭慶林;;一個(gè)改進(jìn)的譜模式非線性項(xiàng)的計(jì)算方法及其應(yīng)用[A];第七屆北京青年科技論文評選獲獎(jiǎng)?wù)撐募痆C];2003年

2 張朝林;鄭慶林;宋青麗;;一個(gè)改進(jìn)的譜模式半隱式時(shí)間差分方案及其應(yīng)用研究[A];大氣科學(xué)發(fā)展戰(zhàn)略——中國氣象學(xué)會第25次全國會員代表大會暨學(xué)術(shù)年會論文集[C];2002年

3 張朝林;鄭慶林;;一個(gè)改進(jìn)的譜模式非線性計(jì)算方法及其在國際通用動(dòng)力學(xué)框架中的應(yīng)用[A];第五屆全國優(yōu)秀青年氣象科技工作者學(xué)術(shù)研討會學(xué)術(shù)論文集[C];2002年

相關(guān)博士學(xué)位論文 前2條

1 楊錦輝;半隱式半拉格朗日非靜力數(shù)值天氣預(yù)報(bào)譜模式動(dòng)力框架若干關(guān)鍵技術(shù)研究[D];國防科學(xué)技術(shù)大學(xué);2015年

2 張朝林;譜模式基本動(dòng)力學(xué)框架的研究[D];中國氣象科學(xué)研究院;2000年

，

本文編號：2084468