87年国科会计划报告.docVIP

半拉格朗日格式特徵線法和雲模式中的地形效應 曾忠一 張博雄 中央研究院 中央氣象局 國立台灣大學 國立台灣大學 摘 要 在進行較小尺度大氣現象的數值模擬時，重力波的影響不能忽略，重力波也會限制時步的大小。雖然採用半隱式格式可以增長時步，但這種格式也會使重力波的相速度與事實不符。本文中使用的特徵線法考慮了重力波項的效應，如此可以採用顯式格式而不需擔心會改變重力波相速度，而且時步也不會受到限制。在一維非線性淺水方程特徵線法中使用單調格式，除了時步不受限制外，還可消除重力波附近的寄生波，也使得邊界條件的給定更為簡單。 在雲模式的測試方面，半拉格朗日法會影響雲微物理參數化格式中水物質項的生成時間，但對流的結構特徵不受到影響。本文最後研究進行地形對對流的影響，在沒有強對流發生的情況下，低層水汽扮演減弱山岳波的角色。若流場為均勻流，則通常在山岳下游會形成類似雲牆的結構。 關鍵詞：半拉格朗日法，特徵線法，雲模式 一、前言 對於中小尺度的模擬，半拉格朗日法採用或不採用單調內插格式，模擬的結果會有明顯不同(張與曾, 1998)。加入單調內插格式可以抑制數值方案所產生的計算性誤差，因而不需要額外的人為平滑。在含有水汽的大氣模式中採用單調內插格式，除了防止產生不合理的負值或水汽過飽和外，還可以避免非物理雜波產生(張與曾, 1999)。雖然許多學者採用多項式做為半拉格朗日法的內插函數，但仍然有其他的內插算符可採用。許多學者在不同坐標的實驗結果表明，構造半拉格朗日法的單調內插算符是可行的。然而要維持形狀保存性或單調性的主要困難在於如何決定內插算符所需具備的性質以及如何運用它們，這並沒有一個絕對的選擇標準。由於近年來半拉格朗日法中小尺度數值模式的飛躍發展，因此必須有合適的單調內插格式來處理模式中複雜的預報變數，以便透過不同的測試來了解這些單調內插格式在不同大氣流場以及推廣到三維數值模式時的表現。 基本方程的重力波項可以用顯式或隱式格式計算。如果採用隱式格式，則時步(time step)將會受到限制。在中尺度模式中，重力波項顯得比較不重要，此時採用半隱式格式(semi-implicit scheme)就能得到不錯的結果。但是對小尺度來說，這時重力波開始影響數值模擬的結果，不僅是由於它的快速移動限制了時步的大小，更因為採用隱式格式將使得重力波模減速，這種情形在重力波模以不同方向移動時更為顯著。 通常處理重力波項的問題可採用顯式或隱式格式，但是這兩種方法各有利弊，因此Erbes(1993)提出利用特徵線法(method of characteristics)來整合這兩種方法的優點。這個方法必須先求得基本方程組的特徵方程，然後利用顯式格式來求解，此時雖然時步可以較大，但仍然沒有解決平流場的單調性問題。本文中我們結合準單調半拉格朗日法(QMSL)(Bermejo and Staniforth, 1992)與特徵線法進行測試。 本研究的雲微物理參數化法採用Lin et al.(1983)的格式，除了一般降水粒子，也考慮冰相粒子對降水結構的影響，使降水過程較為健全。在微物理參數化格式中假設雲滴、雪和冰珠呈Marshall-Palmer分布，降水質點的拖曳力等於質點重量，飽和水汽壓採用Tetens的公式計算。雲內部允許有過飽和的現象，若產生過飽和則採用冰水調整方式修正。進行調整過程時需注意在調整之前先進行半拉格朗日法平流，考慮平流擴散過程，再進行雲微物理過程。由於以往採用的是較短時步下的雲物理過程，因此必須修正雲微物理參數化過程，另外也將探討地形所引發的對流。 半拉格朗日法的效率主要來自於平流項的積分是沿著質點運動軌跡的，因此得到的數值解較尤拉(Euler)格式準確而且引起的誤差較小。只要前一時步質點出發點的估計夠準確，則半拉格朗日法在穩定性上的限制少於尤拉格式，因而才能使用較長的時步。為了建立一個有效率且準確的中小尺度雲模式，我們將進行單調半拉格朗日特徵線法以及雲模式中地形對流的數值實驗。 二、特徵線法的數值實驗 從事較小尺度大氣現象的數值模擬時，重力波會限制時步的大小。半隱式格式雖然可以增長時步，但也會使重力波的相速度與事實不符。Erbes(1993)提出將特徵線法應用在半拉格朗日法中，這個方法考慮了重力波項的效應，如此可以對模式採用顯式格式而不需擔心會改變重力波相速度，且時步也不會有限制。在構造半拉格朗日法特徵線格式時，首先需找出基本方程的特徵方程組，然後將特徵方程組的平流項採用半拉格朗日格式，而特徵值就是平流速度。 任何線性或非線性雙曲線型方程組可以表示如下： , , q為向量，A為係數矩陣，R為A的特徵向量矩陣。我們可以找到N個特