數值天氣預報常用計算方法
正文
數值天氣預報中所用的方程大多是非線性的,迄今還沒有一種解析求解方法,常用的是數值求解方法。其中最常用的是差分法,其次是譜方法。差分法 即用差商代替微商的方法。考慮任意函式f(x,у,t),其偏微商 дf/дx可以用幾種不同的形式來近似表示。如
L.F.理查孫最早將這種方法套用於天氣預報問題。他用中央差代替空間微商,用向前差代替時間微商,認為這樣一步步地計算,就可以作出預報。如對於平流方程:
不過,對於非線性方程,即使線性穩定性判據得到滿足,計算也不一定是穩定的。1959年,N.A.菲利普斯在實際計算中發現,存在一種無論怎樣縮小 Δt 也不能排除的不穩定現象。他認為是由非線性作用產生的短波被虛假地表示了的所謂“混淆現象”引起的,稱為“非線性不穩定性”。這種不穩定性,並非大氣運動中由於物理原因產生的,而是由於在構造差分格式時,破壞了原微分方程的性質而造成的,因此也稱為計算不穩定性。為了克服這個困難,可以考慮使差分方程保留原微分方程所具有的在計算域內總的守恆性質〔如總能量守恆,總位渦度(見大氣動力方程)平方守恆等〕的格式。這時,在時間積分過程中預報量保持有界,因而這種格式對非線性不穩定性有抑制作用。
由於實際需要,有時還須要製作較小範圍的區域性預報。根據實踐,這種預報以用“套格線”製作效果較好。所謂“套格線”,是指在計算域中置兩種以上格線,其中一部分格線距較大,另一部分較小,而細格線計算域又包含在粗格線計算域中。用這種套格線法可以提高局部地區的解析度,從而提高預報的準確率。
譜方法 將微分方程組中函式的空間變化用正交函式的級數的前有限項展開,通過一系列積分運算,使微分方程組變換成以展開係數和其對時間的微商的常微分方程組,以求得近似數值解的方法。正交函式的選擇,依賴於區域的幾何形狀。譜方法通常用來解半球或全球問題,並多用球面調和函式。1954年,有人提出用球面調和函式解正壓渦度方程的譜方法,時間外推方法和差分法相同。以後,雖然有不少人繼續研究,但由於這種方法的計算量很大,特別是非線性項,更是如此。長期以來,這種方法一直停留在研究試驗階段。1970年,A.埃利亞森等利用當時剛發展起來的快速富氏變換計算非線性項,空間微商用譜方法進行,乘積運算在格線點上完成,迴避了直接計算非線性項中相互作用係數的大量工作,使計算量大為減少。
一般說來,譜方法的優點是:①空間微商的計算精確,有利於減小位相誤差;②可以避免非線性混淆現象,使非線性不穩定性不易產生;③便於解泊松方程;④能自動並徹底地濾去短波,比一般差分法中用平滑算符好;⑤解全球或半球問題可以沒有奇異點。缺點是:①表示非線性項所需的計算量和存儲量均較大,計算量隨波數增加得太快;②對分布不太連續的物理量,容易發生跳躍現象,必須用較多的波才能表示;③象降水那樣的局部地區天氣現象和其伴隨的潛熱加熱作用,必須知道整個場的預報量,這就經常要把所有的諧波分量重新組合起來。總之,在研究局地現象時還是採用粗細格線相套的差分法更加靈活方便。