单调半拉格朗日方案在对流云模拟上的应用.doc

PAGE PAGE 8 單調半拉格朗日方案在對流雲模擬上的應用 曾忠一12 王旌熒2 1.中央研究院物理研究所 2.國立台灣大學大氣科學研究所 摘要 由於半拉格朗日法相對於尤拉方案可增加時間步長而且不影響穩定性，因此模式積分較有效率但不失準確性。在使用半拉格朗日法時需要用到內插法來預測下一個時刻的值，但這些內插法常會造成預報值的過量（overshoot）或不足（undershoot），尤其是在梯度很大或不連續的地方更容易如此。以水汽場來說，內插後得到的預報值可能會過飽和或是負值，以致造成微物理過程參數化的困擾。解決這種問題的許多方法被提出來，最重要的是形狀保持方案（shape preserving）或單調方案（monotone）。 本研究使用的方案有準單調方案（QMSL; quasi-monotone semi- Lagrangian）、線性約束方案（LCSL; linear constraint semi-Lagrangian）和五次內插多項式法。QMSL利用高階和低階內插法的組合以建立的一個單調方案，LCSL則是利用出發點鄰近的格點來判斷是否需要加入線性約束條件來達到單調性。使用五次內插多項式配合形狀保持性質的導數估計值，以做為半拉格朗日法的高階內插公式，同樣的也具有單調的性質。 本研究探討二維無輻散旋轉流場、一維水汽平流凝結過程和二維暖雲模式中上述各種單調方案的特性。 1.前言 由於大氣問題的數值模擬中數值解與真解有所不同，多年來學者們致力於改善數值誤差，如Boris and Book（1973, 1976）等。大氣模擬中的一個重要部份是處理平流項，也就是對平流項做數值近似，其中一個常用的方法是半拉格朗日（semi- Lagrangian）平流法。這是一種有效率而準確的平流方案（Bates 1985; Rasch and Williamson 1990; Smolarkiewicz and Pudykiewicz 1992; Priestley 1993）。半拉格朗日法是由某一格點沿風場特徵線回溯到上一時刻的出發點，然後以內插法求出預報值。這種平流方案的優點，除了可在規則格點上進行計算外，還可使用較大的時間步長，因此在計算上比尤拉（Euler）方案更有效率。以往在尤拉方案中，令人困擾的是時間步長受到穩定性的限制。半拉格朗日法，則不受穩定性的限制，只須考慮空間的截斷誤差。Robert（1981）使用三時間層半隱式（semi-implicit）半拉格朗日法的淺水方程模式，研究結果顯示，穩定的時間步長比尤拉方案大6倍。Bermejo（1990）指出，在無輻散流場中使用半拉格朗日法，並以三次樣條函數（cubic spline）進行內插時，質量是守恆的，這是到目前為止具有質量守恆性質的半拉格朗日平流法。 使用半拉格朗日法時需用到內插法來求得預報值，但是在內插出發點處物理場量，會因內插方案的特性而產生非物理性的負值或振盪，尤其在強梯度或不連續處更容易發生。當這種情形發生在水汽時，會產生負值的混合比和過大的水汽混合比，因而導致微物理參數化上的困擾。為了使內插值正確的表達初始物理場，內插方案必須是正定的（positive definite）甚至是單調的。 Boris and Book（1973, 1976）, Zalesak（1979）提出一種在有限差分方案中對通量項做非線性修正的方法，稱為FCT（flux-corrected transport）方案。Carpenter et al.（1990）的PPM（piecewise parabolic method）方案也是一種保持物理場單調性的方法。Williamson and Rasch（1989）指出，使用不同的內插法也可達到物理場的單調性。Bermejo and Staniforth（1992）在內插方案中加入簡單的約束條件，大量的減少數值振盪，防止不合理的負值物理場發生，來維持物理場的單調性。 當方程中沒有強迫項時，單調方案可適當的維持物理場的單調性。Grabowski and Smolarkiewicz（1990）指出，在多個方程間有強迫項的交互作用時，使用單調平流方案還是不足以防止非物理性的數值振盪發生，這是由於強迫項與守恆方程之間的非線性相互作用所造成的影響。Pellerin et al.（1995）指出，即使不使用單調方案的半拉格朗日法，在不連續處的數值振盪也比使用單調方案的尤拉法要小。 本研究的內容分成6個部分，第2節中將介紹本研究使用的三種單調方案，分別是QMSL方案、LCSL方案和五次內插多項式法。第3節是單調半拉格朗日方案的數值實驗，將一缺口圓柱體放在二維無輻散旋轉流場內，將它平流六圈後，來看缺口圓柱體的形狀保持性。第4節是一維平流凝結實驗，使用一組有相互作用