廣義 Rosenau-Kawahara 方程的有效譜方法

文賢 王中慶

摘要：針對廣義Rosenau-Kawahara方程提出了Legendre dual-Petrov-Galerkin 譜方法，并基于對 角化技巧，構建了快速有效算法。在此基礎上研究了單個孤立波的傳播、守恒律及波的生成等物理現象。數值結果驗證了所提算法的有效性。

關鍵詞：Legendre dual-Petrov-Galerkin譜方法；廣義Rosenau-Kawahara方程；孤立波；守恒律

中圖分類號：O 241.82 ?文獻標志碼：A

An efficient spectral method for the generalized Rosenau-Kawahara equation

WEN Xian， WANG Zhongqing

（College of Science， University of Shanghai for Science and Technology， Shanghai 200093， China）

Abstract： The Legendre dual-Petrov-Galerkin spectral method was proposed to the generalized Rosenau-Kawahara equation. A fast and efficient algorithm was constructed based on the diagonalization technique. The motion of single solitary wave solution， conservation laws and the phenomena of wave generation were also studied. Numerical results illustrate the effectiveness of the suggested approach.

Keywords： Legendre dual-Petrov-Galerkin spectral method; generalized Rosenau-Kawahara equation; solitary wave; conservation laws

1 問題提出及預備知識

非線性波的傳播是自然界中普遍存在的物理現象，對其進行數值研究具有十分重要的理論意義和應用價值。在過去的幾十年里， Korteweg-de Vries （ KdV ）方程主要用來描述波動行為。為了克服傳統 KdV 方程無法處理波與波及波與墻相互作用的缺點，在離散動力學系統的研究中經常使用 Rosenau 方程。作為非線性波的進一步考慮，在Rosenau 方程中添加兩個黏性項uxxx和?uxxxxx， 即可 得到 Rosenau-Kawahara 方程[1]。不失一般性，本文 將主要考慮如下廣義 Rosenau-Kawahara 方程的初 邊值問題[2]：

式中： p 為正整數； β， γ， ?，η 為常數； λ 為非負常數。

作為廣義 Rosenau-Kawahara方程的特殊情況， Kawahara 方程、 Rosenau-Kawahara 方程等在流體力學、電磁學、等離子體、非線性光學等領域都有廣泛的應用[1-5]。

文獻[1-2]運用 sine–cosine 和 tanh 等方法構造了廣義 Rosenau-Kawahara方程的漸進解。目前僅有為數不多的文獻討論了廣義 Rosenau-Kawahara 方程的數值方法。如文獻[3]對 Rosenau-Kawahara 方程考慮了兩種守恒有限差分格式，文獻[4]提出了 Rosenau-Kawahara 方程的一種緊致保結構有限差分格式，文獻[2]則發展了廣義 Rosenau-Kawahara 方程的一個守恒有限差分格式。

作為微分方程的主要數值解法之一，譜方法以其高精度優勢漸受關注，并被廣泛應用于流體力學、海洋科學和大氣科學等問題的數值模擬中[5-8]。目前尚未發現使用譜方法求解廣義 Rosenau- Kawahara 方程的相關文獻。若使用傳統譜方法進行求解，則離散后最終得到的系數矩陣是一個十一對角矩陣，因而極大地影響了計算效率。近年來，文獻[9-12]提出了橢圓型問題和 Camassa- Holm 方程等問題的對角化譜方法，其主要優點是離散后的線性系統為對角矩陣，進而從本質上提高了計算效率。本文針對廣義 Rosenau-Kawahara 方程構造了一組 Sobolev 雙正交基函數，由此提出一個有效的 Legendre dual-Petrov-Galerkin 譜方法，并對單個孤立波的傳播、守恒律及波的生成等物理現象進行了數值模擬。

下面回顧 Legendre 多項式的基本知識。記Ln（ξ）為n階 Legendre 多項式，滿足以下三遞推關系[8]：

其正交性為

式中，δmn為 Kronecker 符號。

進一步地，

令PN為不超過N次代數多項式的集合。正如文獻[5]所指出的，針對五階微分方程，由于算子的非對稱性，應使用 dual-Petrov-Galerkin譜方法，即分別考慮如下的測試函數空間和檢驗函數空間：

相應地，兩組基函數選擇如下[5]：

顯然

記?（2k+ m）n+1：=（2k+ m）（2k+2+ m）···（2k+2n+ m），al，k =（?k， ψl）， bl，k =（?k ， ?xψl），cl，k =（?x?k ， ?x（2）ψl）， dl，k =（?x（2）?k ， ?x（3）ψl）。由?Legendre 多項式的性質式（2）～（4）和兩組基函數的定義可得，矩陣A =（alk）k>0，l0，l0， l0， l0，l

2 Legendre dual-Petrov-Galerkin 譜方法

2.1 譜格式

為了使用 Legendre dual-Petrov-Galerkin 譜方法求解原方程（1），進行?（x）= x/L ，= t/L線性變換，仍用（x， t）表示（?（x）， ）?？紤]如下壓縮的廣義 Rosenau-Kawahara 方程：

記時間步長為τ， M = T/τ， u（k）= u（x， kτ）， k =0， 1， ··· ， M。上述方程的標準中心差分格式為

相應的全離散 Legendre dual-Petrov-Galerkin 譜格式如下：

求u +1）∈ VN ， 使得

其中，

2.2 對角化方法

為了更加高效地求解廣義 Rosenau-Kawahara方程，應使用對角化技巧，即重新選擇兩組基函數，使得方程組（6）中的系數矩陣為對角矩陣，從而提高在每個時間層求解方程組（6）的效率。為此，先介紹兩組修正的基函數Φn（x）：=Φn，α，β，?，λ，η（x）和Ψn（x）：=Ψn，α，β，?，λ，η（x） ， 它們關于如下雙線性算子?Sobolev 雙正交，即Bα，β，?，λ，η（u， v）：=α（u， v）?β（u， v x）+?（ux ， v xx）+λ（u xx ， v xx）+η（u xx ， v xxx）。顯然，該算子是式（6）中雙線性算子的一般化。

引理1 設Φ0（x）=?0（x）， Ψ0（x）=ψ0（x），且Φn（x）∈

Vn+5， Ψn ∈?Vn（?）+5是兩組?Sobolev 雙正交多項式，滿足

Φn ??n ∈?Vn+4 ， Ψn ?ψn ∈?Vn（?）+4 ， 且

Bα，β，?，λ，η（Φn， Ψm）=ρnδn，m ，?n， m >0

則下述遞推式成立，即

Φn =?n ? anΦn?1? bnΦn?2? cnΦn?3?

當 n <0時，Φn（x）≡Ψn（x）≡0， ρn =0；當 n <1時， an =0；當n <2時， bn =0；當n <3時， cn =0；當n <4時， dn =0；當n <5時， en =0，且

ρn = b（a） n，ρ（） ，c（n）ρ（β） ndn2（?） ρ（e）n（n）5（η）+ an（2）ρn?1?

an = （an，n?1α+ bn，n?1β+ cn，n?1?+ dn，n?1λ+

n，end（en） an?1ρn?2? cnbn?1ρn?3+ dncn?1ρn?4?

bn = n（2）2ρ（an） encn（+ bn） ）cn，n?2?+ cnan?2ρn?3?

cn =ρ3（an，n?3α+ bn，n?3β+ dnan?3ρn?4? enbn?3ρn?5）

e（d）n（n） α（an），n?4α+ bn，n?4β+ enan?4ρn?5）

限于文章篇幅，具體證明過程可參照文獻[9-12]類似可得。特別地，取兩組修正的 Sobolev 雙正交基函數Φ（?）n（x）和Ψ（?）m（x），滿足

則問題的數值解u +1）可表示為

式（6）在時間方向上是一個隱格式，實際計算包含非線性項的右端項時需要使用迭代方法求解。顯然，式（6）的系數矩陣是對角矩陣，而若使用經典的基函數式（5），則系數矩陣為十一對角矩陣。因此，本文的對角化方法極大地提升了計算效率。對于區間（a， b）上的問題，只需作線性變換

3 數值結果

針對廣義 Rosenau-Kawahara方程給出一系列數值結果，用以描述廣義 Rosenau-Kawahara 方程單個孤立波的傳播、守恒律及波的生成等物理現象。除特殊說明外，以下均取L =100， N =256，τ=10?4。

3.1 孤立波解

眾所周知，孤立波在許多自然科學領域有重要價值。首先考慮 Rosenau-Kawahara方程的孤立波解[1]：

式中，uex（x， t）= Asech4（B （x?vt））。A =35（√156（205）?13），B =1√12（205）?13和?v =分別表示孤立波的振幅、波數和波速。方程的該解表示一個初始以 x =0為中心、波速為v的自左向右傳播的單個孤立波。

采用前述譜格式進行數值求解。圖1表明孤立波在不同時刻以恒定的速度和振幅自左向右傳播，并展示了相應的三維傳播圖?？梢杂^察到數值解的振幅和波速幾乎是常數，且與精確值完全一致。圖2分別給出了數值解的L2和L∞誤差，可以看到，在t =0.5（對應實際時間t =50，下同）時的誤差量級在10?7，該結果表明本文算法提供了比已有文獻更好的數值結果。

下面考慮廣義 Rosenau-Kawahara方程的孤立波解：

式中，?uex（x， t）= Asech4/7（B （x?vt））。

分別表示孤立波的振幅、波數和波速。

采用前述譜格式進行數值求解。圖3表明孤立波在不同時刻以恒定的速度和振幅自左向右傳播，并展示了相應的三維傳播圖形?？梢杂^察到數值解的振幅和波速幾乎是常數，且與精確解完全一致。

3.2 守恒律

對廣義 Rosenau-Kawahara 方程，文獻[2-4]給出了3個基本的守恒量，即質量、動量和能量：

下面考慮對應的離散守恒量及其相對誤差。作為例子，本文僅考慮 Rosenau-Kawahara方程，其初值

在圖4中，分別畫出3個守恒量及其相對誤差隨時間的變化情況。從圖中可以看到， Rosenau-Kawahara 方程的3個守恒量保持得非常好，且滿足：

3.3 波的生成

考慮 Rosenau-Kawahara 方程和初值 u0（x）=κsech4（√√205?13/12x），對于給定的κ， 該初值可用來模擬單個孤立波生成多個波的過程， 取區間I =（?2， 2）。

首先，令κ=6。圖5描繪了波的生成過程。由圖5可以發現，最終有4個孤立波生成，各個孤立波之間的距離隨著時間推移而逐漸增加，且各自以不同的波速和振幅自左向右傳播，振幅和波速大的孤立波會逐步遠離那些振幅和波速小的波。另外，也可以清晰地看到，第2，3，4個孤立波是在傳播過程中生成的，其生成時刻分別約為t =0.03， 0.05， 0.20。特別地，通過數值實驗可以發現，對于給定的正整數κ= M（M +1），上述初值會生成M +2個波。

下面考慮負數κ=?2。圖6描繪了不同時刻波的形狀。數值結果顯示最終沒有孤立波生成，相反地，形成了一系列沿負方向傳播、振幅不斷衰減的振蕩色散波。

參考文獻：

[1]ZUO J M. Solitons and periodic solutions for the Rosenau- KdV ?and Rosenau-Kawahara ?equations[J]. Applied Mathematics and Computation， 2009， 215（2）：835–840.

[2]HE D D. New solitary solutions and a conservative numerical method for the Rosenau-Kawahara equation with?power law nonlinearity[J]. Nonlinear Dynamics， 2015，82（3）：1177–1190.

[3]HU J S， XU Y C， HU B， et al. Two conservative difference schemes for Rosenau-Kawahara equation[J]. Advances in Mathematical Physics， 2014， 2014：217393.

[4] WONGSAIJAI B， CHAROENSAWAN P， CHAOBAN- KOH T， et al. Advance in compact structure-preserving manner to the Rosenau-Kawahara model of shallow-water wave[J]. Mathematical Methods in the Applied Sciences， 2021， 44（8）：7048–7064.

[5] YUAN J M， SHEN J， WU J H. A dual-Petrov-Galerkin method for the Kawahara-type equations[J]. Journal of Scientific Computing， 2008， 34（1）：48–63.

[6] BOYD J P. Chebyshev and Fourier spectral methods[M].2nd ed. Mineola： Dover Publications， 2001.

[7] GUO B Y. Spectral methods and their applications[M]. Singapore： World Scientific， 1998.

[8] SHEN J， TANG T， WANG L L. Spectral methods： algorithms， analysis and applications[M]. Heidelberg： Springer， 2011.

[9] AI Q， LI H Y， WANG Z Q. Diagonalized Legendre?spectral methods using Sobolev orthogonal polynomials for elliptic boundary value problems[J]. Applied Numerical Mathematics， 2018， 127：196–210.

[10] YU X H， YE X Q， WANG Z Q. A fast solver of Legendre- Laguerre spectral element method for the Camassa-Holm equation[J]. Numerical Algorithms， 2021， 88（1）：1–23.

[11] LI S， YAN S M， WANG Z Q. Efficient Legendre dual- Petrov-Galerkin methods ?for ?odd-order differential equations[J]. Discrete and Continuous Dynamical Systems- B， 2020， 25（4）：1543–1563.

[12] ZHANG Y， YU X H， WANG Z Q. Efficient and accurate Legendre spectral element methods for one-dimensional higher order problems[J]. Numerical Mathematics： Theory， Methods and Applications， 2021， 14（2）：461–487.

（編輯：丁紅藝）