雙曲守恒律方程的Lax—Wendroff時間離散WENO格式

李興華+孫陽+艾曉輝

摘 要：雙曲守恒型方程的高精度、高分辨率計算格式的研究一直是計算流體力學的熱點問題。針對原WENOJS格式分辨率較低和計算量偏大的不足問題，提出利用簡單的重構數(shù)值通量的方法以提高計算效率，構造了新的簡單限制器的5階迎風型WENO格式。通過MATLAB軟件的仿真對LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式的實驗結果進行了分析，并比較了這四種計算格式的計算效率和計算精度。數(shù)值實驗表明：新格式LaxWendroff簡單限制器WENO格式在保持原WENO分辨率的前提下，計算速度有明顯提高，減少了20%的計算時間。

關鍵詞：高精度；WENO；RungeKutta；LaxWendroff；時間離散

DOI：10.15938/j.jhust.2017.06.026

中圖分類號： O175

文獻標志碼： A

文章編號： 1007-2683（2017）06-0134-06

Abstract：The research of high accuracy and high resolution schemes have been a hot topic in computational mathematics. According to low resolution and large amount of calculation of the original WENOJS scheme， we propose a simple new limiter fifth order upwind WENO scheme to reconstruct the numerical flux of the simple structure to improve the computational efficiency. Compared with other efficient high accuracy schemes such as ENO and WENO， it is shown that the computational cost of this scheme is less than that of WENOJS in the same accuracy. By use of MATLAB software， we compared and analyzed computational efficiencies and computational accuracies of LaxWendroff WENOJS scheme， LaxWendroff simple limiter WENO scheme， RungeKutta simple limiter WENO scheme and RungeKutta WENOJS scheme. The numerical results show that the new LaxWendroff simple limiter WENO scheme can improve the computing speed and reduce the computing time by 20% while maintaining the original WENO resolution.

Keywords：high accuracy； WENO； RungeKutta； LaxWendroff； time discretization

0 引 言

雙曲守恒律方程（組）為科學理論和工程應用研究中一類非常重要的偏微分方程（組）。空氣動力學、爆炸力學、流體力學等許多力學問題的求解都與其密切相關。由于即使在初始條件充分光滑的條件下，雙曲型守恒律方程的解仍可能出現(xiàn)間斷。因此，為了能夠更加高效地捕捉到間斷，并且避免間斷附近出現(xiàn)數(shù)值振蕩，計算流體力學及計算數(shù)學等領域一直致力于研究穩(wěn)定高效的數(shù)值計算方法[1-5]。

低階精度的離散方法可能會對許多流動結構造成“失真”模擬，不能真實反映真實的流動現(xiàn)象。高精度格式具有較小的耗散誤差和色散誤差，能夠更準確地模擬雙曲守恒律方程的時空多尺度特性，在同樣網格條件下比低階格式能分辨出更加精細的流場、捕捉到其更細微的結構變化。然而，高階精度格式也存在不足：求解在間斷處可能出現(xiàn)Gibbs現(xiàn)象即產生偽物理振蕩，從而導致非線性不穩(wěn)定現(xiàn)象。因此，發(fā)展高分辨率、高精度離散格式是計算流體力學發(fā)展中的迫切需要，成為計算流體力學工作者的一大研究方向[6-10]。

TVD格式的出現(xiàn)，標志著計算流體力學步入了高精度計算格式階段，大量高精度計算格式不斷出現(xiàn)，而WENO格式由于其良好的計算穩(wěn)定性，已經成為一類極為重要的計算格式。眾所周知，經典的WENOJS格式存在分辨率較低和計算量偏大兩方面問題，因此計算數(shù)學工作者嘗試對其改進，構造了諸如WENOZ、WENOM等很多改進格式。根據(jù)文獻[11]的基本思想構造的簡化限制器的5階WENO格式是對原WENO格式的很好的補充和修正。

在時間離散方面，本文研究了LaxWendroff型時間離散格式，其具有較RungeKutta方法更為優(yōu)秀的計算效率。本文將WENO格式與LaxWendroff型時間離散相結合，得到了一種時空同步離散的新數(shù)值方法——基于LaxWendroff型時間離散的WENO格式。本文研究可以說是對高精度數(shù)值方法研究方面的有益補充和探索。

1 WENO5格式

針對以上兩種WENO格構造不同的時間離散格式，能夠得到不同的數(shù)值方法。針對雙曲守恒律方程給出了具體的格式構造過程，主要采用如下介紹的兩種時間離散方法LaxWendroff時間離散和TVD RungeKutta時間離散。下面，我們將介紹這兩種常用的時間離散方法。endprint

2 時間離散格式

目前雙曲守恒律方程最為常見的時間離散格式為TVD RungeKutta時間離散格式，也稱為SSP（strong stability preserving）RungeKutta時間離散。當一階向前Euler迭代下的空間離散格式為TVD時，所構造的全離散格式也TVD的。本文采用的RungeKutta時間離散格式為三階RungeKutta格式（簡記為RK3）[12-16]如下：

3 結果與分析

考慮線性對流方程和非線性方程兩種類型的方程的數(shù)值求解[17-20]，利用LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式四種格式進行數(shù)值模擬，研究數(shù)值方法的穩(wěn)定性、收斂性和計算效率等特性，所采用多的通量分裂為LaxFredrichs通量分裂。

表4和表5分別給出LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式計算結果。表4為四種格式的L1誤差和精度表，而表5為四種格式的L1誤差和精度表。

具有如下初值u（x，0）=0.5+sin（πx），計算到t=32π時刻，邊界條件為周期條件，N=400個計算點。此時在1.2～1.3之間的位置產生了一強激波，圖1和圖2為其計算結果。

計算結果中簡單限制器的WENO格式僅僅用到了2-3個計算點就可以識別激波，因此可以斷定簡單限制器的WENO格式的分辨率要高于WENOJS格式，而計算時間上，LaxWendroff時間離散具有一定的優(yōu)勢，WENO LW3和SWENO LW3兩種格式的計算時間為最少。

4 結 論

針對線性和非線性標量守恒律方程進行了數(shù)值求解。研究了不同初值條件下的，LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式的計算效率和計算精度。LaxWendroff簡單限制器WENO格式不但形式簡單且便于應用到CFD代碼中。算例表明，LaxWendroff WENO格式的分辨率與RungeKutta型WENO格式較為接近，而LaxWendroff簡單限制器WENO格式具有更好的激波及小尺度波分辨率，特別計算速度提高了WENO5的20%以上。本文從分析標量雙曲守恒律方程解的性質入手，分析了雙曲守恒律方程自身的特點和原有的WENOJS差分格式的優(yōu)缺點后，提出一種簡化限制器的WENO格式，新格式較原格式的分辨率有所提高。LaxWendroff時間離散格式由于其計算效率方面的優(yōu)勢，可以作為RungeKutta方法的一種時間離散改進格式。本文研究了四種不同的數(shù)值方法：LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式計算了一維標量守恒型方程，并給出了數(shù)值實驗的結果。數(shù)值試驗表明這四種方法均出了高精度高分辨的特性，都具有更好的分辨激波捕捉的能力，在激波間斷處不但保持了陡峭的的形狀，而且沒有出現(xiàn)偽物理振蕩。我們發(fā)現(xiàn)，對于同樣的問題，LaxWendroff簡單限制器WENO方法能在更短的時間和更小的儲存量。當然，從數(shù)值實驗的結果來看，方法仍有改進的空間。未來需要關注的主要方面在于：在精度相同的情況下， LaxWendroff型的計算量要小于RungeKutta，因而更具有效率上的優(yōu)勢價值。然而，在推廣到雙曲守恒律方程組時，需要進行張量運算這給編程帶來了一定的困難。下一步作者將試著將推廣到高維守恒律方程方程組問題，現(xiàn)公式推到已經完成，正在進行數(shù)值實驗，以得到較優(yōu)結果。

參 考 文 獻：

[1] 陸金甫. 偏微分方程數(shù)值解法[M]. 北京：清華大學出版社， 2004.

[2] 袁光偉，杭旭登，盛志強，等. 福射擴散計算方法若干研究進展[M]. 計算物理，2009，26（4）： 475-500.

[3] 歐陽頎. 非線性科學與斑圖動力學導論[M]. 北京：北京大學出版社，2010.

[4] 劉儒勛，舒其望. 計算流體力學的若干新方法[M]. 北京：科學出版社，2003.

[5] WARMING R F. Flux Vector Splitting of the Inviscid GasDynamic Equation with Application of Finite Difference Method [J]. Journal of Computational， 1981，40： 263-293.

[6] GODONOV S K. A Difference Schemes for Numerical Computation of Discontinuous Solution of Hydrodynamics Equations [J]. Math. Sbornik， 1959，47： 271-306.

[7] ENGQUIST，OSHER S. Stable and Entropy Satisfying Approximation for Transonic Flow Calculation [J]. Mathematics of Computation， 1980，34： 45-75.

[8] COCKBURN B，LIN S Y，SHU C W. TVB RungeKutta Local Projection Discontinuous Galerkin Finite Element Method for Conservation Laws III： OneDimensional Systems [J]. J. Comput. Phys.，1989，84： 90-113.endprint

[9] COCKBURN B，HOU S， SHU C W. TVB RungeKutta Local Projection Discontinuous Galerkin Finite element method for conservation laws IV [J]. Math. Comput.，1990，54： 545-581.

[10]LIU H，XU H，GONG H. Modeling the Asymmetry in Traffic Flow （b）：Macroscopic Approach [J]. Appl. Math. Model，2013，37（22）： 9441-9450.

[11]ZHU J，QIU J X. A New Fifth Order Finite Difference WENO Scheme for Solving Hyperbolic Conservation Laws [J]. J. Comput. Phys.，2016，318： 110-121.

[12]LIU X， ZHANG S，ZHANG H，SHU C W. A New Class of Central Compact Schemes with Spectrallike Resolution II： Hybrid Weighted Nonlinear Schemes [J]. J. Comput. Phys.，2015，284： 133-154.

[13]PIROZZOLI S. Conservative Hybrid CompactWENO Schemes for ShockTurbulence Interaction [J]. J. Comput. Phys.，2002，178 （1）： 81-117.

[14]QIU J X， DUMBSER M， SHU C W. The Discontinuous Galerkin Method with LaxWendroff type Time Discretizations [J]. Comput. Methods Appl. Mech. Engrg. 2005， 194： 4528-4543.

[15]LAX P D. Hyperbolic System of Conservation Laws II [J]. Comm. Pure and Applied Mathematics. 1957，10： 537-566.

[16]ZHANG X，SHU C W. On MaximumPrincipleSatisfying High Order Schemes for Scalar Conservation Laws [J]. J. Comput. Phys.，2010，229 （9）： 3091-3120.

[17]ZHANG X，SHU C W. MaximumPrincipleSatisfying and PositivityPreserving HighOrder Schemes for Conservation Laws： Survey and New Developments Proc. R. Soc. A [J]. Math. Phys. Eng. Sci.，2011，467（2）： 2752-2776.

[18]LIU Y J，SHU C W，TADMOR E，ZHANG M P. Central Discontinuous Galerkin Methods on Overlapping Cells with A NonOscillatory Hierarchical Reconstruction[J]. SIAM J. Numer. Anal.，2007， 45：2442-2467.

[19]姚慧麗，宋曉秋，李興華. 一類半線性微分方程的漸進概自守溫和解 [J].哈爾濱理工大學學報，2012，17（1）： 72-78.

[20]姚慧麗，卜憲江，宋曉秋. 一類微分方程的指數(shù)增長的溫和漸進概自守解[J]. 哈爾濱理工大學學報，2014，19（5）：23-26.

（編輯：王 萍）endprint