亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        復(fù)Ginzburg-Landau方程的格子Boltzmann模型分析解

        2019-11-28 11:39:50張建影閆廣武

        張建影,閆廣武,李 婷,3

        (1.長春工業(yè)大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,長春 130012;2.吉林大學(xué) 數(shù)學(xué)學(xué)院,長春 130012;3.吉林化工學(xué)院 理學(xué)院,吉林 吉林 132022)

        復(fù)Ginzburg-Landau方程在許多物理領(lǐng)域中應(yīng)用廣泛[1-2].文獻(xiàn)[3-4]研究表明,用格子Boltzmann方法(LBM)求解復(fù)Ginzburg-Landau方程具有簡單有效、數(shù)值精度高等優(yōu)點(diǎn).該方法的思想是通過構(gòu)造格子Boltzmann方程中的平衡態(tài)分布函數(shù)得到所模擬的宏觀方程,通過尋找平衡態(tài)分布函數(shù)的矩函數(shù)給出平衡態(tài)分布函數(shù)與宏觀量之間的變換,再進(jìn)一步通過格子Boltzmann方程給出下一時刻的分布函數(shù)[5].在某種條件下,上述分布函數(shù)存在分析解[6-7].本文利用Chapman分析方法,給出系列偏微分方程以及Chapman多項(xiàng)式的一般形式,通過求解復(fù)Ginzburg-Landau方程的平衡態(tài)分布函數(shù),給出不同時間尺度上的分布函數(shù)表達(dá)式,從而不需要格子Boltzmann方程迭代可以直接給出分布函數(shù),進(jìn)而得到復(fù)Ginzburg-Landau方程的格子Boltzmann分析解.

        1 系列偏微分方程及Chapman多項(xiàng)式

        考慮二維的FHP(Frisch-Hasslacher-Pomeau)網(wǎng)格,在位置x、時刻t定義具有粒子速度eα的復(fù)分布函數(shù)Fα(x,t),其實(shí)部和虛部可視為兩種單粒子分布.復(fù)變量A(x,t)定義為

        (1)

        為得到穩(wěn)定的宏觀量A(x,t),假設(shè)分布函數(shù)Fα(x,t)具有平衡態(tài),即

        (2)

        定義Knudsen數(shù)ε=l/L,其中:l為粒子平均自由程;L為特征長度.選擇時間步長Δt與Knudsen數(shù)相等[3],則復(fù)格子Boltzmann方程為

        (3)

        其中: 實(shí)常數(shù)τ為單松弛時間因子;變量ωα(x,t)為非碰撞項(xiàng),表示分布函數(shù)的增量, 假設(shè)其具有多尺度形式:

        (4)

        在小Knudsen數(shù)的假設(shè)下,對Fα(x,t)做Chapman-Enskog展開,得

        (5)

        (6)

        引入時間多尺度t0,t1,…,滿足

        tn=εnt,n=0,1,…,

        (7)

        將式(5),(7)代入式(3),并做Taylor展開,可得前6個不同時間尺度上的系列偏微分方程[3]:

        在方程(8)~(13)中,

        (14)

        ωα(x,t)=ε2θα(x,t),

        (15)

        則由式(8),(9)可得到系列偏微分方程的一般表達(dá)式:

        (16)

        其中:

        (17)

        表達(dá)式(16)對于i>6的情形仍然成立.

        2 復(fù)Ginzburg-Landau方程的格子Boltzmann模型

        其中:δkj表示Kronecker符號;λ為由Knudsen數(shù)和松弛因子確定的實(shí)參數(shù).對式(8)兩端關(guān)于α求和,得到一階宏觀方程,即t0時間尺度上的守恒律:

        (21)

        將式(8)+式(9)×ε,再關(guān)于α求和,可得二階宏觀方程:

        (22)

        式(22)是具有二階截?cái)嗾`差O(ε2)的復(fù)Ginzburg-Landau方程,其中

        λε(τ-1/2)=1,

        (23)

        (24)

        (25)

        (26)

        在式(25),(26)中,D(=2)為空間維度,b(=6)為連接到相鄰節(jié)點(diǎn)的方向數(shù),c=|eα|為粒子運(yùn)動速度,參數(shù)β可以是實(shí)數(shù)或復(fù)數(shù).由于方程中不包含源項(xiàng)的對流項(xiàng),故假設(shè)θα和α(α=1,2,…,b)無關(guān)[3],則可得

        θα(x,t)=ξH(A),α=1,2,…,b,

        (27)

        θ0(x,t)=(1/ε-bξ)H(A),

        (28)

        其中ξ是待定參數(shù).

        3 任意階分布函數(shù)的計(jì)算

        式(16)也可寫成

        將式(9)+式(10)×ε+…+式(29)×εi-1,再關(guān)于α求和,得

        (30)

        在式(30)中,誤差項(xiàng)Ei-1為

        (31)

        (35)

        作用于θα的算子Δ的指數(shù)部分滿足如下關(guān)系:

        (36)

        p+s-Nsj,s-1=i-Nij,i-1,

        (37)

        p+s-2-Nsj,s-3=i-2-Nij,i-3.

        (38)

        因此,式(34)可以寫成

        (41)

        從而求得分布函數(shù)Fα解析解的級數(shù)形式.用上述方法可以得到任意階精度的解析解.例如,若取到前三階,則Fα為

        將平衡態(tài)分布(25),(26)及附加分布(27),(28)代入式(42),再由式(21)可得

        (44)

        將式(9),(10)關(guān)于α求和,并由式(21)可得

        (45)

        (46)

        將其代入式(43),(44)可得Fα的解析解形式為

        綜上可見,當(dāng)構(gòu)造復(fù)Ginzburg-Landau方程的格子Boltzmann 模型時,由于使用的附加項(xiàng)是小Knudsen數(shù)的二階假設(shè),使得任意階的系列偏微分方程變成顯式,從而可得到不同尺度上的分布函數(shù).該方法對其他非線性偏微分方程的格子Boltzmann模型解析解的構(gòu)造有參考作用.

        青青草视频在线免费视频 | 超碰97人人射妻| 亚洲国产精品情侣视频| 久久精品国产av麻豆五月丁| 亚洲av精二区三区日韩| 国色天香精品一卡2卡3卡4| 97久久超碰国产精品2021| 人妻被猛烈进入中文字幕| 久久亚洲一区二区三区四区五| 后入丝袜美腿在线观看| 乱老年女人伦免费视频| 免费人成视频x8x8入口| 久久99精品国产99久久| 国产人成无码中文字幕| 日本精品人妻在线观看| 少妇高潮呻吟求饶视频网站| 人妻少妇被粗大爽视频| 亚洲综合av大全色婷婷| 国产精品成人aaaaa网站| 日本高清www午色夜高清视频 | 乱人伦人妻中文字幕不卡| 美利坚合众国亚洲视频| 国产成人亚洲精品无码青| 色狠狠色噜噜av天堂一区| 欧美最猛性xxxxx免费| 国产真人无遮挡作爱免费视频| 国产啪精品视频网给免丝袜| 亚洲愉拍自拍视频一区| 国产免费观看久久黄av麻豆| 欧美牲交a欧美牲交aⅴ免费真| 少妇的丰满3中文字幕| 一二三四中文字幕日韩乱码| 国产在线一区二区三区不卡| 精品亚洲一区二区三区四区五| 亚洲精品无码久久久久av老牛| 亚洲精品成人专区在线观看| 国产在线精品福利大全| 亚洲人妻御姐中文字幕| 国产免费久久精品99久久| 亚洲精品suv精品一区二区| 成人免费网站视频www|