高 云， 樂勵華

(東華理工大學(xué)數(shù)學(xué)與信息科學(xué)學(xué)院，江西撫州 344000)

新的格子Bhatnagar-Gross-Krook模型求解修正的Burgers方程

高 云， 樂勵華

(東華理工大學(xué)數(shù)學(xué)與信息科學(xué)學(xué)院，江西撫州 344000)

隨著計算機技術(shù)的發(fā)展，數(shù)值模擬方法求解偏微分方程得到越來越廣泛的應(yīng)用。格子Boltzmann方法是一種新型的模擬方法，由于該方法具有計算效率高、邊界條件容易處理、完全并行性等獨特的優(yōu)點，使得它具有廣泛的應(yīng)用領(lǐng)域。利用格子Bhatnagar-Gross-Krook模型來求解修正的Burgers方程，首先用該方法正確的恢復(fù)了宏觀方程，然后數(shù)值模擬了兩個具有解析解的修正Burgers方程。把模擬解與解析解進行對比，發(fā)現(xiàn)數(shù)值解與解析解和前人研究中的數(shù)值解都吻合很好。

格子BGK模型;修正的burgers方程;數(shù)值解

高云，樂勵華.2012.新的格子Bhatnagar-Gross-Krook模型求解修正的Burgers方程［J］.東華理工大學(xué)學(xué)報:自然科學(xué)版，35(1):89-93.

GaoYun，Le Li-hua.2012.A New Lattice Bhatnagar-Gross-Krook Model for the Modified Burgers’Equation［J］.Journal of East China Institute of Technology(Natural Science Edition)，35(1):89-93.

本研究考慮的修正burgers方程為

其中Re是描述粘性系數(shù)v(=1/Re)的Reynolds數(shù)。

修正的burgers方程是重要的非線性對流擴散方程。Bateman(1915)第一次提出burgers方程，之后，很多學(xué)者作了進一步的研究和發(fā)展。近年來，修正的burgers方程被應(yīng)用于實際的輸運問題中，例如湍流、波在熱彈性介質(zhì)中的傳播、河流和沉淀物中污染物的運輸和散射、低頻吸收的媒介中非線性波的傳播等。許多學(xué)者一直在構(gòu)造修正burgers方程的高效數(shù)值方法。Irk(2009)采用六次樣條函數(shù)的配置方法、Bratsos(2010)用有限差分方法、Ramadan(2005a)采用五次樣條插值的配置方法來求解修正burgers方程。

格子Boltzmann方法在模擬流體復(fù)雜的物理問題時是一種有效的方法。與傳統(tǒng)的計算方法相比，該方法具有許多獨特的優(yōu)勢，如計算效率高、邊界條件容易處理、具有完全并行性等。格子Bhatnagar-Gross-Krook(BGK)模型是目前格子Boltzmann方法研究和應(yīng)用的主要模型之一，已經(jīng)在多孔介質(zhì)流(Xu et al.，2004)、血液流(Fang et al.，2004)、多相流(Hou，1997)等領(lǐng)域取得了很大的成功。由于格子Boltzmann方法的獨特優(yōu)勢，目前，該方法正逐步應(yīng)用于模擬非線性系統(tǒng)，如對流擴散方程(Shi et al.，2009，2010)，反應(yīng)擴散方程(Chen et al.，1995)，泊松方程(Chai et al.，2008)，KdV方程(Chai et al.，2008)等。

本文在基于格子BGK模型和Chapman-Enskog展開的基礎(chǔ)上來求解方程(1)，發(fā)現(xiàn)數(shù)值解與解析解和其它模型的數(shù)值解都吻合很好。

1 格子BGK模型

使用Qian等(1992)提出的D1Q3模型，該格子BGK模型分布函數(shù)的演化方程為:

盡管，閻廣武等(1999)應(yīng)用多尺度展開技術(shù)，通過選擇平衡態(tài)分布函數(shù)的高階矩的格子Boltzmann方法研究了Burgers’方程，對于一維模型，將速度離散成4個方向，每一個節(jié)點與相鄰近的4個節(jié)點相連，即D1Q4模型。Duan等(2008)針對該方程提出2-bit模型，即D1Q2模型，將速度離散成2個方向，每一個節(jié)點與相鄰近的2個節(jié)點相連，通過Taylor展開和多尺度分析，從格子Boltzmann方程恢復(fù)了Burgers’方程。但是，對于Burgers’方程本研究選取D1Q3模型，介于即D1Q4模型與D1Q2模型之間，本研究模型更簡單，精度也更高。

2 數(shù)值模擬

為了驗證D1Q3格子BGK模型的有效性，模擬了修正 Burgers方程。在模擬中，使用了 Guo等(2002)提出的非平衡態(tài)外推方法來處理邊界條件。為了檢驗數(shù)值方法的有效性，誤差定義為:

在模擬中，選取 α=0.9，Δx=0.005，Δt=0.01。表1給出數(shù)值模擬結(jié)果與解析解的比較，從中可以看出數(shù)值解與解析解和 Ramadand等(2005b)模擬的解差別很小。表2顯示了誤差的比較，可以看出本文的誤差較小。表3給出了不同粘性系數(shù)v的誤差，v越小的，誤差越小。因此，格子BGK模型更適合小的粘性系數(shù)v。另外，還模擬了v=0.001和v=0.000 5兩種情況(圖1，圖2)，從圖中可以看出數(shù)值解和解析解非常吻合。

表1 數(shù)值解與解析解的比較在v=0.001，t=4.0Table 1 Comparison of numerical solution and analytic solution while v=0.001，t=4.0

表2 v=0.001，Δx=0.005，Δt=0.01時誤差比較Table 2 Error comparison while v=0.001，Δx=0.005，Δt=0.01

表3 不同粘性系數(shù)v的誤差比較Table 3 Error comparison of different v of viscosity coefficient

選取Δx=0.01，Δt=0.01，α=0.1，v=0.005進行模擬。然后將所得的數(shù)值結(jié)果與其它文獻結(jié)果進行L2，L∞誤差比較(表4)。發(fā)現(xiàn)格子BGK模型的誤差更小。圖3顯示了粘性系數(shù)時的數(shù)值結(jié)果。圖4顯示了數(shù)值解與解析解的誤差情況。

表4 v=0.004在不同時間的誤差比較Table 4 Error comparison with different time while v=0.004

4 結(jié)語

本研究在基于格子BGK模型的基礎(chǔ)上，對修正的Burgers方程進行了數(shù)值模擬。發(fā)現(xiàn)本文模擬的數(shù)值結(jié)果和方程的解析解吻合很好。由于采用D1Q3模型，將速度離散成3個方向，每一個節(jié)點與相鄰近的3個節(jié)點相連，簡化了算法的同時，也減少了誤差。同時，也與以前文獻的結(jié)果也進行了比較，發(fā)現(xiàn)本研究的誤差更小。對模型穩(wěn)定性分析是下一步工作的目標(biāo)。

閻廣武.1999.用格子Boltzmann方法研究Burgers方程［J］，力學(xué)學(xué)報，31:143-151.

Bateman H.1915.Some recent researches on the motion of fluids［J］.Monthly Weather Rev，43:163-170.

Bratsos A G.2010.A fourth-order numerical scheme for solving the modified Burgers equation［J］.Comput.Math.Appl 60:1393-1400.

Bratsos A G，Petrakis L A.2011.An explicit numerical scheme for the modified Burgers’equation［J］.Int.J.Numer.Meth.Biomed.Engng 27:232-237.

Chai Z H，Shi B C，Zheng L.2008.A unified lattice Boltzmann model for some nonlinear partial differential equations［J］.Chaos Solitons Fractals，36:874-882

Chai Z H，Shi B C.2008.A novel lattice Boltzmann model for the Poisson equation［J］.Appl.Math.Model 32:2050-2058.

Chen S，Dawson S P，Doolen G D，et al.1995.Lattice methods and their applications to reacting systems［J］.Comput.Chem.Eng.，19:617-646.

Duan Y L，Liu R X，Jiang Y Q，2008.Lattice Boltzmann model for the modified Burgers’equation，Appl.Math.Comput 202:489-497.

Fang H P，Chen S Y.2004.Boltzmann method for three-dimensional moving particles in a Newtonian fluid［J］.Chin.Phys.，13:47-53.

Guo Z L，Zheng C G，Shi B C.2002.Non-equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method［J］.Chin.Phys.，11:366-374.

Hou S，Shan X，Zou Q，et al.1997.Evaluation of two lattice Boltzmann models for multiphase flows［J］.J.Comput.Phys.，138:695-713.

Irk D.2009.Sextic B-spline collocation method for the modified Burgers equation［J］.Kybernetes，38:1599-1620.

Qian Y，d’Humiéres D，Lallemand P.1992.Lattice BGK models for Navier-Stokes equation［J］.Europhysics Letters.，17:479-484.

Ramadan M A，El-Danaf T S.2005a.Numerical treatment for the modified Burgers equation［J］，Math.Comput.Simul，70:90-98.

Ramadan M A，El-Danaf T S，Alaal F.2005b.A numerical solution of the Burgers’equation using septic B-splines［J］.Chaos Solitons Fractals，26:1249-1258.

Saka B，Da 2008.A numerical study of the Burgers’equation［J］.J.Franklin Inst 345:328-348.

Shi B C，Guo Z L.2009.Lattice Boltzmann model for nonlinear convection-diffusion equations［J］.Phys.Rev.E，79:016701.

Shi B C，Guo Z L.2011.Lattice Boltzmann simulation of some nonlinear convection-diffusion equations［J］.Comput.Math.Appl.，61:3443-3452.

Xu Y S，Zhong Y J，Huang G X.2004.Lattice Boltzmann method for diffusion-reaction-transport processes in heterogeneous porous media［J］.Chin.Phys.Lett.，21:1298-1301.

A New Lattice Bhatnagar-Gross-Krook Model for the Modified Burgers’Equation

GAO Yun， LE Li-hua
(Dept.of Mathematics＆Information Science，East China Institute of Technology，F(xiàn)uzhou，JX 344000，China)

The methods of numerical simulation are applied widely with the development of computer technology.The lattice Boltzmann method has a more wide application owing to its high computing efficiency，easier processed boundary conditions，and parallel computation.In this study，the Bhatnagar-Gross-Krook model is used to solve the modified Burgers’equation.First，the macroscopic equation is recovered with this model.Then，the two modified Burgers’equations with analytical solutions are numerical simulated.Results show numerical simulation and analytical solutions match well with numerical simulation of previous publication.

LBGK model;modified burgers’equation;numerical solution

O24

A

1674-3504(2012)01-089-05

10.3969/j.issn.1674-3504.2012.01.013

2011-09-18 責(zé)任編輯:張國慶

高 云(1986—)，女，碩士生，研究方向:微分方程。E-mail:gaoyun1986@126.com