• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    捕捉激波的群速度控制方法

    2014-04-30 07:24:44李新亮傅德薰馬延文
    空氣動力學(xué)學(xué)報 2014年5期
    關(guān)鍵詞:群速度波數(shù)激波

    李新亮,傅德薰,馬延文

    (中國科學(xué)院力學(xué)研究所 高溫氣體動力學(xué)國家重點實驗室,北京 100190)

    捕捉激波的群速度控制方法

    李新亮,傅德薰,馬延文

    (中國科學(xué)院力學(xué)研究所 高溫氣體動力學(xué)國家重點實驗室,北京 100190)

    回顧及綜述了群速度控制方法基本原理、發(fā)展過程及最新進展。群速度控制方法是通過控制色散誤差以捕捉間斷的一種數(shù)值方法。該方法通過在間斷兩側(cè)采用不同色散類型的格式,將數(shù)值波的傳播壓制在間斷附近,避免數(shù)值振蕩向周圍擴散,從而達到抑制數(shù)值振蕩的效果。該方法主體上采用數(shù)值色散抑制振蕩,實際使用中僅需補充較低水平的數(shù)值耗散就可將振蕩壓制在可接受的范圍內(nèi)。因而群速度控制格式的總體耗散水平較低,易于實現(xiàn)小尺度波的高分辨率捕捉,是復(fù)雜流動高分辨率數(shù)值模擬的有效方法。

    群速度控制格式;激波捕捉格式;高分辨率;色散;耗散

    0 引 言

    高精度激波捕捉方法在計算流體力學(xué)領(lǐng)域占有重要地位。該方法要求具有較高的間斷分辨率(即,數(shù)值解間斷盡量銳利),且間斷周圍數(shù)值振蕩盡量小,而在非間斷區(qū)(光滑區(qū))需具有較高的數(shù)值精度及分辨率[1]。傳統(tǒng)意義上講,不低于二階精度的激波捕捉方法均可稱為高精度激波捕捉方法。代表性的高精度激波捕捉格式有TVD格式[2],NND格式[3],WENO格式[4]以及高階保單調(diào)格式[5]等。這些格式在光滑區(qū)(非間斷區(qū))采用較高數(shù)值精度及較低數(shù)值耗散的方法,以保證較高的數(shù)值分辨率,而在間斷附近則采用較低數(shù)值精度及較高數(shù)值耗散的格式以抑制數(shù)值振蕩。這些差分格式大多是通過數(shù)值粘性壓制數(shù)值波的傳播,從而起到抑制振蕩的效果。

    群速度控制方法(Group Velocity Control Method,GVC)是傅德薰等人提出的一種新型激波捕捉方法[6-7],與TVD及WENO等格式的機制不同,該格式并不直接通過耗散機制抑制激波,而是通過色散控制捕捉激波。該格式利用色散誤差控制數(shù)值波傳播的群速度,使得數(shù)值傳播的方向指向間斷,以避免從間斷發(fā)出的數(shù)值振蕩污染流場。由于采用數(shù)值色散抑制振蕩,因而群速度控制格式總體上數(shù)值耗散較小。且該方法實現(xiàn)簡潔,其計算量明顯低于WENO等高精度激波捕捉方法。利用該思想,先后構(gòu)造了緊致型群速度控制格式[8],優(yōu)化型群速度控制格式[9],加權(quán)型群速度控制格式[10]以及助波器型群速度控制格式[11,7]等多種群速度控制格式,為含間斷復(fù)雜流動的數(shù)值模擬提供了有效的高分辨率計算方法。

    1 群速度控制方法的基本思想

    1.1 數(shù)值波的群速度:快格式、慢格式和混合格式

    以線性對流方程

    為例,討論數(shù)值波的群速度。

    對于方程(1),假設(shè)初值為如下單波:

    u(x,0)=eikx

    其中

    為修正波數(shù)。對于給定的線性格式,該修正波數(shù)可以通過理論分析得到[12]。而對于非線性格式,該修正波數(shù)的特性也可通過數(shù)值方法近似獲得[13-14]。

    該修正波數(shù)的虛部ki=ki(α)反應(yīng)了格式的色散誤差,即數(shù)值波傳播的相位(或波速)誤差。當ki(α)=α時,數(shù)值波傳播速度與精確解相同;當ki(α)>α時,數(shù)值波傳播快于精確解;反之,則慢于精確解。而對于一般情況,數(shù)值解中包含了多個波的傳播,因而觀測到波傳播的速度通常為波包的速度,即波的群速度,該群速度定義為[6-7]:

    在整個波數(shù)范圍內(nèi)(0<α≤π),D0(α)>0的格式稱為“快格式”,對于快格式,數(shù)值波傳播的群速度超過精確解;在整個波數(shù)范圍內(nèi)D0(α)<0的格式稱為“慢格式”,對于慢格式,數(shù)值波傳播的速度低于精確解;在部分波數(shù)范圍內(nèi)D0(α)>0,在另外部分波數(shù)范圍內(nèi)D0(α)<0的格式稱為“混合格式”。圖1給出了某些差分格式的群速度D0(α)=dki/dα,從中可以看出2階中心格式以及6階對稱緊致格式為“慢格式”;2階迎風格式以及5階迎風緊致格式為“混合格式”,這兩種格式在低波數(shù)范圍內(nèi)具有“快格式”特征。

    圖1 某些差分格式的群速度[6]Fig.1 Group velocities for four finite difference schemes[6]

    1.2 數(shù)值振蕩產(chǎn)生原因及群速度控制的基本思想

    群速度控制理論將間斷周圍數(shù)值振蕩的原因歸結(jié)為由色散誤差導(dǎo)致的群速度誤差。假設(shè)方程(1)的初值具有間斷,則根據(jù)Fourier分解,該間斷可分解為不同波數(shù)波的疊加,如圖2所示。如果利用快格式進行數(shù)值計算,這些數(shù)值波傳播的速度均快于理論速度a,而且不同波傳播的群速度不同,通常情況下高波數(shù)成分的傳播群速度更快。經(jīng)過一段時間后高波數(shù)成分傳播到了主波(間斷)的前方,在間斷前方形成了數(shù)值振蕩,而波后振蕩不明顯。同樣,利用慢格式計算時,所有數(shù)值波傳播的速度均低于理論速度a,但高波數(shù)成分傳播得更慢,因而數(shù)值解解在波后有明顯振蕩,而波前振蕩不明顯。

    因而,如果混合使用“快格式”和“慢格式”,在波前使用“慢格式”,在波后使用“快格式”,則可以盡量壓縮數(shù)值振蕩的傳播范圍,從而起到抑制數(shù)值振蕩的作用。這就是群速度控制的基本思想[6-7]。

    圖2 間斷分解為若干波的疊加Fig.2 Diagrammatic sketch for the decomposition of a discontinuity

    實際情況中,純粹的“快格式”很難構(gòu)造,而多數(shù)格式為混合型格式(如圖1所示)。很多情況下利用低波數(shù)區(qū)為“快格式”的“混合格式”來代替純“快格式”。這種情況下,由于這些“混合格式”在高波數(shù)范圍具有“慢格式”特征,會影響群速度控制的效果。但通常這些格式在高波數(shù)區(qū)具有較強的數(shù)值粘性,具有抑制振蕩效果,因而數(shù)值振蕩通常會被壓制在允許范圍內(nèi)。

    2 群速度控制方法的具體實現(xiàn)及應(yīng)用

    2.1 二階精度的群速度控制格式

    在文獻[6]中,傅德薰等除了介紹了群速度控制的基本思想,還給出了一種具有二階精度的群速度控制格式,這是最早的群速度控制格式。該格式利用2階迎風型格式作為“快格式/混合格式(MXD)”,利用二階中心格式作為“慢格式(SLW)”。此外,該方法利用激波形狀函數(shù)(Shock-Structure function,簡稱SS函數(shù))判斷激波前與激波后。SS函數(shù)的定義為:

    對于向右傳播的間斷波,SS(u)=-1為波前;SS(u)=1為波后。圖3為數(shù)值解中間斷波的示意圖,圖中標示了SS函數(shù)如何區(qū)分波前和波后。

    對于正通量,該二階精度群速度控制格式的數(shù)值通量為:

    其中:

    分別為作為“快格式/混合格式”使用的二階迎風格式以及作為“慢格式”使用的二階中心格式的數(shù)值通量。SS函數(shù)的值可由式(4)數(shù)值離散給出。

    圖3 數(shù)值解中傳播的間斷波示意圖[8]Fig.3 Diagrammatic sketch for a N-Sshock[8]

    2.2 緊致型群速度控制格式

    在文獻[8]中,馬延文和傅德薰利用文獻[6]的群速度思想,構(gòu)造了緊致型群速度控制格式。該方法利用文獻[6]的三階迎風緊致格式作為“快格式”(實際上,該格式為混合格式,但在低波數(shù)范圍內(nèi)呈現(xiàn)快格式特征);利用四階中心緊致格式[12]作為“慢格式”。利用以上兩種基礎(chǔ)格式,構(gòu)造出一種含自由參數(shù)的緊致型差分格式[7-8,15]:

    其中Fj/Δx是函數(shù)?f/?x的逼近式,

    分別為一階及二階導(dǎo)數(shù)中心差分算子。σ為自由參數(shù)。顯然,當σ=0時,格式為四階對稱緊致差分(“慢格式”);當σ=±1/6時,格式為三階迎風緊致差分(“快格式/混合格式”)。通過調(diào)節(jié)該自由參數(shù)的數(shù)值,在激波(間斷)前實現(xiàn)慢格式,在激波(間斷)后實現(xiàn)快格式,從而實現(xiàn)群速度控制。

    對于正通量,緊致型群速度控制格式的具體表達式如下[8]:

    其中

    1≤σ0≤2,0.8≤γ0≤1是格式中的常數(shù)。δ-fj=fjfj-1,δ+fj=fj+1-fj是向后、向前差分算子。P和c分別為當?shù)氐膲毫奥曀?,利用這兩個參數(shù)的目的是探測激波的強度,并由此強化或弱化σ的變化范圍。

    對于負通量的緊致型群速度控制格式同理構(gòu)造,具體表達式見文獻[8]。

    圖4為采用該群速度控制格式計算Sod激波管問題,得到i=0.14時刻的流場分布圖(圖片引自文獻[8])。從中可以看出,該計算具有很高的激波分辨率及較低的數(shù)值振蕩水平。從該格式的表達式可以看出,格式并沒有人為添加額外的數(shù)值耗散。計算結(jié)果顯示,格式的整體耗散保持在較低的水平,大致與迎風緊致格式持平,但激波捕捉能力得到了很大提升。需要指出的是,在文獻[8]中,并沒有利用到局部特征分解,而是直接基于Steger-Warming分裂后的通量進行差分計算。如果利用了局部特征分解,則激波附近的數(shù)值振蕩還可以進一步降低[16]。

    圖4 Sod激波管問題t=0.14時刻的流場分布[8]Fig.4 Distributions of flow at t=0.14 of Sod problem[8]

    2.3 八階精度的優(yōu)化型群速度控制格式

    在文獻[9]中,李新亮等利用群速度控制思想,構(gòu)造了八階精度的普通型(非緊致)群速度控制格式(GVC8)。在格式的構(gòu)造過程中還采用了基于算例的優(yōu)化技術(shù),針對具體算例(Sod激波管問題),以魯棒性為目標,對格式的系數(shù)進行優(yōu)化,從而在保證計算精度及激波分辨率的情況下,提高了方法的魯棒性。

    針對方程(1),構(gòu)造線性差分格式:

    其中:

    格式使用9個網(wǎng)格點構(gòu)造數(shù)值通量,理論上最高可具有9階精度。GVC8要求格式具有8階精度,因而,有一個參數(shù)可以自由調(diào)節(jié)。通過優(yōu)化調(diào)節(jié)該系數(shù),得到兩組不同的線性差分格式,使得一組為“快格式”(實際上是混合格式),另外一組為“慢格式”。在格式的構(gòu)造過程中,系數(shù)的優(yōu)化除了考慮格式的群速度特性,還根據(jù)算例進行了系數(shù)優(yōu)化,使得格式具有更好的魯棒性。

    在該格式采用更為簡潔的判據(jù)以判斷某點位于激波的前方還是激波的后方。判據(jù)為:對于向右傳播的波,如果|uj-uj-1|<|uj+1-uj|,則認為該點(j點)位于波后,否則,則認為該點位于波前。如示意圖5所示,通常情況下,越靠近激波(間斷),數(shù)值振蕩越激烈,因而振蕩(變差)大的方向為激波所在的位置。

    圖5 激波附近數(shù)值解示意圖Fig.5 Diagrammatic sketch for numerical solution near a shock

    對于正通量(方程(1),a>0),GVC8格式為式(6),格式的系數(shù)為:

    系數(shù)的具體值如表1[9]。

    表1 的值[9]Table 1 Values of[9]

    文獻[9]通過Shu-Osher問題以及可壓縮各向同性湍流DNS對該格式進行了測試,結(jié)果顯示,該格式具有較低的數(shù)值耗散及較高的尺度分辨率。此外,該格式的另外一個優(yōu)點是構(gòu)造簡潔,計算量小。數(shù)值測試顯示,其計算量接近9階精度的線性格式,遠低于5階精度的WENO格式。該格式的不足是魯棒性略弱,尤其是對于高Mach數(shù)的算例。

    2.4 加權(quán)型群速度控制格式(WGVC)

    利用群速度控制思想,何志偉等人構(gòu)造了加權(quán)群速度控制格式(WGVC,Weighted Group Velocity Control Scheme)[10],并進一步構(gòu)建了WGVC-WENO混合格式,提升了格式的魯棒性。

    WGVC格式離散的正通量表達式為:

    其中fj+1/2(σ)是含有自由參數(shù)σ,且經(jīng)過優(yōu)化的線性格式。隨著參數(shù)σ的變化,該格式可以表現(xiàn)出“慢格式”以及“混合格式/快格式”的特性。兩個自由參數(shù)σm和σs以及權(quán)重ωm、ωs根據(jù)流場特性數(shù)值計算產(chǎn)生[10]。在波前區(qū),σs≈1,σm≈0,格式(7)表現(xiàn)出“慢格式”特征;反之在波后區(qū)格式表現(xiàn)出“混合格式/快格式”特征。格式的群速度特性由參數(shù)σm和σs根據(jù)流場的光滑性質(zhì)自適應(yīng)給出。在此基礎(chǔ)上,何志偉等人構(gòu)造了WGVC-WENO混合格式,在間斷區(qū)借助WENO格式提升格式的魯棒性。圖6為Shu-Osher問題的計算結(jié)果(201網(wǎng)格點),從中可以看出,與7階WENO格式相比,7階精度的WGVC-WENO混合格式具有更好的數(shù)值分辨率。

    2.5 群速度助推法

    在文獻[11]中,馬延文等提出了一種群速度助波器(Group Velocity Booster,GVB),可將常規(guī)線性差分格式改造成群速度控制格式。

    對于微分方程:

    圖6 Shu-Osher問題計算結(jié)果:WGVC-WENO7 vs.WENO7,201網(wǎng)格點[10]Fig.6 Distributions of density for Shu-Osher problem:WGVC-WENO7 vs.WENO7,201 points[10]

    其中Dh為“慢格式”型差分算子。實際上,很多迎風型差分算子為慢格式(例如三、五、七階迎風偏心型差分)[7]。通過添加“助波器”,可將式(8)改造成群速度控制格式。改造后的格式為:

    其中,算子WB為群速度助波器,通過添加該助波器,可以將原先的“慢格式”改造成群速度控制格式。常見的助波器形式為[7,11]:

    2.6 其他群速度控制格式

    除了前文介紹的方法,還有很多研究者利用群速度控制思想進行了激波捕捉格式的構(gòu)建及應(yīng)用。

    Zhu等人[17]運用群速度控制思想并結(jié)合系數(shù)優(yōu)化,構(gòu)建了緊致型群速度控制格式,并通過飛船返回艙高速繞流算例進行了測試。結(jié)果顯示,該方法具有較好的魯棒性及激波分辨率。周超紅等人[18]利用群速度控制思想,構(gòu)造了三階及四階精度的群速度控制格式,并將該方法推廣到求解淺水波方法,數(shù)值實驗表明,該方法可以很好地控制非物理振蕩并能保證計算精度。田保林等人[19]根據(jù)數(shù)值解的群速度特性,在三階迎風緊致格式的基礎(chǔ)上引入了群速度控制機制,使其能夠正確模擬含有強激波和強接觸間斷的復(fù)雜流動。計算結(jié)果表明,該方法對激波和接觸間斷的分辨效果是令人滿意的。Cheng等人[20]在原群速度控制格式的基礎(chǔ)上,進一步減小了數(shù)值耗散,構(gòu)造了一種低耗散的群速度控制格式。高慧等人構(gòu)造具有六階精度的緊致型群速度控制差分格式[21],并通過Sod激波管問題、定常激波問題以及二維激波反射等問題進行了測試,結(jié)果顯示,該方法具有網(wǎng)格基架點小和捕捉激波能力較強的優(yōu)點。

    3 結(jié)論及展望

    低耗散、高分辨率激波捕捉方法是近年來計算流體力學(xué)領(lǐng)域研究的熱點。群速度控制方法則是近年來提出的一種新型的激波捕捉方法。不同于大多數(shù)通過耗散壓制數(shù)值振蕩的方法,群速度控制方法主張通過控制色散誤差抑制數(shù)值振蕩。該方法根據(jù)數(shù)值波傳播的特性,將格式區(qū)分為“快格式”和“慢格式”,并通過在間斷前使用“慢格式”,在間斷后使用“快格式”的方法,限制數(shù)值振蕩的傳播,從而起到抑制間斷附近數(shù)值振蕩的效果。由于該方法主體上采用數(shù)值耗散抑制振蕩,實際使用過程中僅需再添加少量數(shù)值耗散就可將振蕩抑制在可接受的范圍內(nèi)。因而,該格式的總體耗散水平較低,適用于多尺度復(fù)雜流動的高分辨率數(shù)值模擬。此外,與WENO等格式相比,群速度控制格式實現(xiàn)簡單,計算效率較高。

    群速度控制格式也有需要改進之處。通常情況下純粹的“快格式”比較難于構(gòu)造,因而經(jīng)常使用“混合格式”代替“快格式”,這種情況下高波數(shù)振蕩仍需要一定的耗散來壓制。但由于大多數(shù)“混合格式”本身的耗散特性具有壓制高波數(shù)振蕩能力,因而通常狀況下無需額外添加粘性(或僅需添加較弱的格式粘性即可)。此外,群速度控制格式的魯棒性與WENO等格式相比仍略顯不足,這可以通過構(gòu)建混合型格式等方法來彌補。相信隨著數(shù)值計算技術(shù)的完善和發(fā)展,群速度控制格式也會逐漸發(fā)展并會獲得較好的應(yīng)用。

    致謝:感謝國家超級計算天津中心,中國科學(xué)院網(wǎng)絡(luò)中心超級計算中心以及國家超級計算濟南中心提供計算機時。

    [1]PIROZZOLIS.Numerical methods for high-speed flows[J].Annual Review Fluid Mechanics,2011,43:94-163.

    [2]HARTEN A.High resolution schemes for hyperbolic systems of conservation laws[J].Journal of Computational Physics,1983,49:357-393.

    [3]張涵信.無波動、無自由參數(shù)的耗散差分格式[J].空氣動力學(xué)報,1988,6(2):143-164.

    [4]JIANG G S,SHU C W.Efficient implementation of weighted ENO schemes[J].Journal of Computational Physics,1996,126:202-228.

    [5]SURESH A,HUYNH H T.Accurate monotonicity-preserving schemes with runge-kutta time stepping[J].Journal of Computational Physics,1997,136:83-99.

    [6]FU D X,MA Y W.A high order accurate difference scheme for complex flow fields[J].Journal of Computational Physics,1997,134:1-15.

    [7]傅德薰,馬延文,李新亮,等.可壓縮湍流直接數(shù)值模擬[M].科學(xué)出版社,2010.

    [8]MA Y W,F(xiàn)U D X.Fourth order accurate compact scheme with group velocity control(GVC)[J].Science in China A,2001,44(9):1197-1204.

    [9]LI X L,F(xiàn)U D X,MA Y W.Optimized group velocity control scheme and DNS of decaying compressible turbulence of relativehigh turbulent Mach number[J].International Journal for Numerical Methods in Fluids,2005,48:835-852.

    [10]HE Z W,LI X L,LIANG X.Nonlinear spectral-like schemes for hybrid schemes[J].Science China:Physics,Mechanics& Astronomy,2014,57(4):753-763.

    [11]MA Y W,F(xiàn)U D X,LI X L.A new high order accurate shock capture method with wave booster[C].Proceeding of 2nd International Symposium on Computational Mechanics,Hong Kong and Macau,Nov 30-Dec 3,2009.

    [12]LELE S K.Compact finite difference schemes with spectral-like resolution[J].Journal of Computational Physics,1992,103:16-42.

    [13]PIROZZOLI S.On the spectral properties of shock-capturing schemes[J].J.Comput.Phys.,2006,219:489-497.

    [14]LI X L,LENG Y,HE Z W.Optimized sixth-order monotonicity-preserving scheme by nonlinear spectral analysis international[J].Journal for Numerical Methods in Fluids,2013,73:560-577.

    [15]MA Y W,F(xiàn)U D X,Scheme construction with numerical flux residual correction(NFRC)and group velocity control(GVC)[J].Progress in Natural Science,2006,16(12):1252-1259.

    [16]HE Z W,LI X L,F(xiàn)U D X,et al.Monotonicity-preserving upwind compact difference schemes[J].Science in China G,2011,54(3):511-522.

    [17]ZHU Q Y,LI Y.An upwind compact approach with group velocity control for compressible flow fields[J].International Journal for Numerical Methods in Fluids,2004,44(5):463-482.

    [18]ZHOU C H,YAO Q H,ZHU Q Y.The modified finite difference schemes with group velocity control for solving 2D shallow water equations[J].ACTA Scientiarum Naturalium Universitatis Sungatseni,2012,51(1):19-24.(in Chinese)

    周超紅,姚清河,朱慶勇.求解二維淺水波方程的群速度修正法[J].中山大學(xué)學(xué)報(自然科學(xué)版),2012,51(1):19-24.

    [19]TIAN B L,F(xiàn)U D X,MA Y W.Group velocity control scheme and two-dimensional Riemann slover[J].Chinese Journal of Computational Mechanics,2005,22(1):104-108.(in Chinese)

    田保林,傅德薰,馬延文.群速度控制格式及二維Riemann解[J].計算力學(xué)學(xué)報,

    [20]CHENG J B,F(xiàn)U D X,MA Y W.Group velocity control scheme with low dissipation[J].Chinese Journal of Aeronautics,2000,13(3):138-145.

    [21]GAO H,MA Y W,F(xiàn)U D X.Sixth order accurate compact scheme with group velocity control and application[J].Journal of Aerospace Power,2003,18(1):24-31.(in Chinese)

    高慧,馬延文,傅德薰.六階精度的群速度直接控制緊致格式及其應(yīng)用[J].航空動力學(xué)報,2003,18(1):24-31.

    Review for group velocity control schemes

    LI Xinliang,F(xiàn)U Dexun,MA Yanwen
    (LHD,Institute of Mechanics,Chinese Academy of Sciences,Beijing 100190,China)

    The basic principles and the recent developments of group velocity control(GVC)method are reviewed.GVC method is a new shock capture method by controlling the dispersion error of numerical schemes.In GVC method,the“fast”scheme is used in front of the shock and the“slow”one is used behind the shock,and in this way,the spread of numerical oscillations is limited in a very small region near the shock,and thus the oscillations are suppressed.The GVC method capture shock mainly by using dispersion technique,therefore,the dissipation of GVC schemes is relatively low.The recently developments of GVC schemes,such as compact GVC scheme,optimized GVC scheme,hybrid weighted GVC-WENO scheme and GVC scheme by group velocity booster,are also reviewed.GVC schemes are useful methods for detail simulation of multi-scale flows.

    group velocity control;shock capture scheme;high resolution,dispersion;dissipation

    V211.3

    Adoi:10.7638/kqdlxxb-2014.0104

    0258-1825(2014)05-0575-06

    2014-08-20;

    2014-09-10

    國家自然科學(xué)基金(1372330,11472010,11472278);863項目(2012AA01A304);中國科學(xué)院知識創(chuàng)新工程項目(KJCX2-EWJ01,XXH12503-02-02-04)

    李新亮(1972-),男,研究員,主要研究方向:計算流體力學(xué),湍流.E-mail:lixl@imech.a(chǎn)c.cn

    李新亮,傅德薰,馬延文.捕捉激波的群速度控制方法[J].空氣動力學(xué)學(xué)報,2014,32(5):575-580.

    10.7638/kqdlxxb-2014.0104. LI X L,F(xiàn)U D X,MA Y W.Review for group velocity control schemes[J].ACTA Aerodynamica Sinica,2014,32(5):575-580.

    猜你喜歡
    群速度波數(shù)激波
    聲場波數(shù)積分截斷波數(shù)自適應(yīng)選取方法
    一種基于SOM神經(jīng)網(wǎng)絡(luò)中藥材分類識別系統(tǒng)
    電子測試(2022年16期)2022-10-17 09:32:26
    VTI介質(zhì)地震波群速度三維擴展各向異性線性近似表征
    激光調(diào)制摻鉺光纖中光速及其在高靈敏光纖傳感領(lǐng)域的潛在應(yīng)用研究
    一種基于聚類分析的二維激波模式識別算法
    基于HIFiRE-2超燃發(fā)動機內(nèi)流道的激波邊界層干擾分析
    物質(zhì)波的波速與頻率公式中的能量
    斜激波入射V形鈍前緣溢流口激波干擾研究
    小議超光速
    適于可壓縮多尺度流動的緊致型激波捕捉格式
    这个男人来自地球电影免费观看| 麻豆国产av国片精品| 嫁个100分男人电影在线观看 | 热99国产精品久久久久久7| 日本欧美国产在线视频| 另类亚洲欧美激情| 国产不卡av网站在线观看| 19禁男女啪啪无遮挡网站| 老司机影院毛片| 免费黄频网站在线观看国产| 国产xxxxx性猛交| 看十八女毛片水多多多| 免费在线观看黄色视频的| 9色porny在线观看| 老汉色∧v一级毛片| 久久毛片免费看一区二区三区| 亚洲av日韩在线播放| 桃花免费在线播放| 一本综合久久免费| 巨乳人妻的诱惑在线观看| 欧美黄色片欧美黄色片| 大香蕉久久网| 精品国产乱码久久久久久男人| 久久99一区二区三区| 欧美精品一区二区免费开放| 18禁黄网站禁片午夜丰满| 久久久久国产一级毛片高清牌| www.精华液| 丁香六月欧美| 女性被躁到高潮视频| 精品少妇一区二区三区视频日本电影| 午夜91福利影院| 欧美 亚洲 国产 日韩一| 欧美精品啪啪一区二区三区 | 黄色怎么调成土黄色| 热re99久久国产66热| 99国产精品一区二区三区| 一级黄片播放器| 久久久久久久国产电影| 亚洲国产精品国产精品| 日本五十路高清| 日韩电影二区| 黄色视频不卡| 51午夜福利影视在线观看| 一区在线观看完整版| 91老司机精品| 亚洲国产成人一精品久久久| 人人妻人人澡人人爽人人夜夜| 国产极品粉嫩免费观看在线| 国产在视频线精品| av欧美777| 伊人亚洲综合成人网| 高清av免费在线| 老司机影院成人| 国产免费又黄又爽又色| 1024视频免费在线观看| 久久精品久久精品一区二区三区| 亚洲伊人色综图| 亚洲人成77777在线视频| 亚洲欧美成人综合另类久久久| 国精品久久久久久国模美| 丝袜在线中文字幕| 国产男女超爽视频在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美乱码精品一区二区三区| 国产成人系列免费观看| 国产女主播在线喷水免费视频网站| 婷婷丁香在线五月| 又紧又爽又黄一区二区| 另类亚洲欧美激情| 成年人黄色毛片网站| 高潮久久久久久久久久久不卡| 飞空精品影院首页| 国产精品av久久久久免费| 99热网站在线观看| 国产精品一区二区精品视频观看| 中文字幕人妻丝袜一区二区| 99精品久久久久人妻精品| av国产精品久久久久影院| 悠悠久久av| 女性生殖器流出的白浆| 中文字幕色久视频| 亚洲精品一卡2卡三卡4卡5卡 | 国产一区二区激情短视频 | 婷婷色麻豆天堂久久| 这个男人来自地球电影免费观看| 久9热在线精品视频| 狂野欧美激情性xxxx| 久久久久精品国产欧美久久久 | 国产男女内射视频| 亚洲国产精品一区三区| 国产成人免费无遮挡视频| 亚洲欧美中文字幕日韩二区| 999精品在线视频| 国产在线观看jvid| 99re6热这里在线精品视频| 国产免费福利视频在线观看| 亚洲欧美成人综合另类久久久| 女人爽到高潮嗷嗷叫在线视频| 国产亚洲av片在线观看秒播厂| 欧美激情高清一区二区三区| 欧美精品av麻豆av| 亚洲精品中文字幕在线视频| 久久鲁丝午夜福利片| 少妇人妻 视频| 国产免费一区二区三区四区乱码| 亚洲欧美色中文字幕在线| 男女下面插进去视频免费观看| 亚洲精品国产一区二区精华液| 狂野欧美激情性bbbbbb| www.av在线官网国产| 国产麻豆69| av电影中文网址| 国产免费又黄又爽又色| 国产一区二区三区av在线| 国产av国产精品国产| 日本五十路高清| 男女免费视频国产| 日本av手机在线免费观看| 中文字幕亚洲精品专区| 男女免费视频国产| 亚洲精品久久午夜乱码| 国产亚洲av片在线观看秒播厂| 男女下面插进去视频免费观看| av一本久久久久| 精品亚洲成国产av| 黄频高清免费视频| 黑人欧美特级aaaaaa片| 自线自在国产av| 免费观看人在逋| 日本猛色少妇xxxxx猛交久久| 欧美成人午夜精品| 波野结衣二区三区在线| 亚洲欧洲精品一区二区精品久久久| 男人爽女人下面视频在线观看| 亚洲精品一区蜜桃| 国产午夜精品一二区理论片| 汤姆久久久久久久影院中文字幕| 一区二区三区四区激情视频| 国产在线观看jvid| 一区二区三区精品91| av在线播放精品| 精品亚洲成a人片在线观看| 满18在线观看网站| 中文字幕av电影在线播放| 久久久国产欧美日韩av| 欧美亚洲 丝袜 人妻 在线| xxxhd国产人妻xxx| 少妇精品久久久久久久| 国产色视频综合| av国产久精品久网站免费入址| 国产免费又黄又爽又色| 欧美精品一区二区免费开放| 夜夜骑夜夜射夜夜干| 日韩欧美一区视频在线观看| 久9热在线精品视频| 亚洲国产最新在线播放| 国产精品久久久人人做人人爽| 制服诱惑二区| 欧美国产精品va在线观看不卡| 久久亚洲国产成人精品v| av福利片在线| 欧美日韩精品网址| 成人亚洲精品一区在线观看| 中文字幕最新亚洲高清| 日本91视频免费播放| 日韩人妻精品一区2区三区| 欧美大码av| 婷婷色av中文字幕| 99国产综合亚洲精品| 熟女少妇亚洲综合色aaa.| 欧美亚洲日本最大视频资源| 亚洲黑人精品在线| 黄色毛片三级朝国网站| 久久久精品免费免费高清| 青草久久国产| 日本欧美国产在线视频| 五月开心婷婷网| 亚洲国产欧美在线一区| 在线亚洲精品国产二区图片欧美| 国产免费福利视频在线观看| 真人做人爱边吃奶动态| 18禁国产床啪视频网站| 久久久久国产一级毛片高清牌| 黄网站色视频无遮挡免费观看| 国产成人精品在线电影| 国产一区二区 视频在线| 亚洲精品一区蜜桃| 久久天躁狠狠躁夜夜2o2o | 国产精品av久久久久免费| 男女高潮啪啪啪动态图| av在线老鸭窝| 婷婷色综合大香蕉| 五月天丁香电影| 久久亚洲国产成人精品v| 丝袜人妻中文字幕| 伦理电影免费视频| 国产精品久久久久久精品古装| 91麻豆精品激情在线观看国产 | 国产av精品麻豆| 桃花免费在线播放| 亚洲精品国产av成人精品| av线在线观看网站| 久久久国产精品麻豆| 大片电影免费在线观看免费| 国产深夜福利视频在线观看| 一区二区三区精品91| 一级毛片我不卡| 80岁老熟妇乱子伦牲交| 一本色道久久久久久精品综合| 一边摸一边做爽爽视频免费| 久久99一区二区三区| av在线播放精品| 国产精品香港三级国产av潘金莲 | 国产精品麻豆人妻色哟哟久久| 嫁个100分男人电影在线观看 | 丝袜人妻中文字幕| 亚洲自偷自拍图片 自拍| 高清视频免费观看一区二区| 色综合欧美亚洲国产小说| 欧美日韩黄片免| 久久人妻熟女aⅴ| 亚洲精品久久午夜乱码| 1024香蕉在线观看| 久久久国产一区二区| 中文字幕亚洲精品专区| 亚洲国产欧美一区二区综合| 不卡av一区二区三区| 国产欧美日韩精品亚洲av| 啦啦啦在线观看免费高清www| 亚洲精品中文字幕在线视频| 这个男人来自地球电影免费观看| 男人爽女人下面视频在线观看| 五月开心婷婷网| 国产高清videossex| 国产亚洲精品第一综合不卡| 在线观看免费午夜福利视频| 建设人人有责人人尽责人人享有的| 狠狠婷婷综合久久久久久88av| 久久久久久久久免费视频了| 另类亚洲欧美激情| 久久久国产一区二区| 久久人人97超碰香蕉20202| 久久天躁狠狠躁夜夜2o2o | 亚洲精品日韩在线中文字幕| 免费日韩欧美在线观看| 亚洲人成77777在线视频| 高清av免费在线| 午夜两性在线视频| 在线 av 中文字幕| 久热爱精品视频在线9| 亚洲国产精品国产精品| 少妇裸体淫交视频免费看高清 | 久久久国产精品麻豆| 国产亚洲午夜精品一区二区久久| 三上悠亚av全集在线观看| 高清黄色对白视频在线免费看| 天堂8中文在线网| 丝袜喷水一区| 青草久久国产| 免费日韩欧美在线观看| 日韩,欧美,国产一区二区三区| 精品久久蜜臀av无| 亚洲av综合色区一区| 性高湖久久久久久久久免费观看| 亚洲国产av新网站| 亚洲 欧美一区二区三区| 在线亚洲精品国产二区图片欧美| 亚洲成色77777| 国产极品粉嫩免费观看在线| 久久精品久久精品一区二区三区| 水蜜桃什么品种好| 精品一品国产午夜福利视频| 亚洲国产精品一区二区三区在线| xxx大片免费视频| 久久青草综合色| 男女边吃奶边做爰视频| 99热全是精品| 午夜免费观看性视频| 欧美亚洲 丝袜 人妻 在线| www.熟女人妻精品国产| 日本欧美国产在线视频| 一个人免费看片子| 一二三四社区在线视频社区8| 亚洲av欧美aⅴ国产| 久久久精品94久久精品| 国产黄色视频一区二区在线观看| 欧美变态另类bdsm刘玥| 日韩人妻精品一区2区三区| 国产一区亚洲一区在线观看| 亚洲av片天天在线观看| 在线亚洲精品国产二区图片欧美| 亚洲国产日韩一区二区| 亚洲自偷自拍图片 自拍| 久久精品aⅴ一区二区三区四区| 最新在线观看一区二区三区 | 国产成人一区二区在线| 国产日韩一区二区三区精品不卡| 蜜桃在线观看..| 一边亲一边摸免费视频| 免费黄频网站在线观看国产| 久久久久久亚洲精品国产蜜桃av| 国产日韩欧美亚洲二区| 国产麻豆69| 久久精品久久久久久久性| 国产福利在线免费观看视频| av电影中文网址| 国产一区二区三区综合在线观看| 男女边吃奶边做爰视频| av在线app专区| 亚洲熟女精品中文字幕| 亚洲熟女毛片儿| 黄频高清免费视频| 亚洲国产最新在线播放| 国产精品九九99| 又粗又硬又长又爽又黄的视频| 久久性视频一级片| 国产精品99久久99久久久不卡| 两性夫妻黄色片| 麻豆国产av国片精品| 免费在线观看黄色视频的| 精品久久久久久电影网| 日韩大码丰满熟妇| 伊人亚洲综合成人网| av有码第一页| 国产精品久久久久久人妻精品电影 | 亚洲国产av新网站| 精品欧美一区二区三区在线| 亚洲精品久久午夜乱码| 99九九在线精品视频| 国产日韩欧美视频二区| 国产在线一区二区三区精| 91麻豆精品激情在线观看国产 | 91精品三级在线观看| 久久精品人人爽人人爽视色| 亚洲国产欧美一区二区综合| 免费高清在线观看视频在线观看| 亚洲欧洲国产日韩| 亚洲欧洲日产国产| 美国免费a级毛片| 精品国产超薄肉色丝袜足j| 丝袜喷水一区| av有码第一页| 美女主播在线视频| 国产成人啪精品午夜网站| 中文字幕高清在线视频| 麻豆乱淫一区二区| tube8黄色片| 国产在线免费精品| 欧美日韩视频精品一区| 国产亚洲午夜精品一区二区久久| 亚洲av片天天在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 欧美变态另类bdsm刘玥| 伊人亚洲综合成人网| 丰满人妻熟妇乱又伦精品不卡| 99国产精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲精品久久久久5区| 18禁裸乳无遮挡动漫免费视频| 欧美乱码精品一区二区三区| 午夜免费成人在线视频| 人人妻人人澡人人爽人人夜夜| 午夜老司机福利片| 欧美日韩国产mv在线观看视频| 午夜视频精品福利| 国产欧美亚洲国产| bbb黄色大片| 一区二区三区乱码不卡18| 午夜影院在线不卡| 最近最新中文字幕大全免费视频 | 丝瓜视频免费看黄片| 美女高潮到喷水免费观看| 麻豆国产av国片精品| 亚洲成人手机| 亚洲专区中文字幕在线| 亚洲av日韩精品久久久久久密 | 久久99热这里只频精品6学生| 欧美日韩黄片免| 免费高清在线观看日韩| 菩萨蛮人人尽说江南好唐韦庄| 王馨瑶露胸无遮挡在线观看| 国产精品久久久人人做人人爽| 日韩,欧美,国产一区二区三区| 欧美黑人欧美精品刺激| 51午夜福利影视在线观看| 国产一区二区三区av在线| 91老司机精品| 亚洲国产成人一精品久久久| 亚洲精品美女久久av网站| 亚洲色图 男人天堂 中文字幕| 久久99热这里只频精品6学生| 老司机靠b影院| 啦啦啦 在线观看视频| 99热国产这里只有精品6| 国产1区2区3区精品| 乱人伦中国视频| 人人妻人人澡人人爽人人夜夜| 国产97色在线日韩免费| 高清欧美精品videossex| 天天影视国产精品| 1024视频免费在线观看| 脱女人内裤的视频| 国产爽快片一区二区三区| 亚洲九九香蕉| www.精华液| 热re99久久精品国产66热6| 看十八女毛片水多多多| 成人国语在线视频| 一边摸一边抽搐一进一出视频| 黄片小视频在线播放| 黄色片一级片一级黄色片| 超碰成人久久| 宅男免费午夜| 亚洲欧洲国产日韩| 老司机亚洲免费影院| 亚洲第一青青草原| 国产国语露脸激情在线看| 久久精品亚洲熟妇少妇任你| a级片在线免费高清观看视频| 99精品久久久久人妻精品| 91九色精品人成在线观看| 少妇被粗大的猛进出69影院| 国产精品成人在线| 99久久综合免费| 亚洲人成网站在线观看播放| 午夜av观看不卡| 午夜福利影视在线免费观看| 久久天堂一区二区三区四区| 国产极品粉嫩免费观看在线| 热99国产精品久久久久久7| 欧美老熟妇乱子伦牲交| 日本一区二区免费在线视频| 热99国产精品久久久久久7| 精品免费久久久久久久清纯 | 国产黄色免费在线视频| 亚洲欧洲精品一区二区精品久久久| 超碰97精品在线观看| 国产在线视频一区二区| 亚洲国产av影院在线观看| 亚洲精品日本国产第一区| 午夜免费男女啪啪视频观看| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产国语对白av| 国产淫语在线视频| 免费高清在线观看视频在线观看| 国产av国产精品国产| videos熟女内射| 男的添女的下面高潮视频| 国产精品一区二区在线不卡| 欧美 日韩 精品 国产| 天天添夜夜摸| 一级毛片电影观看| 妹子高潮喷水视频| 99精国产麻豆久久婷婷| 欧美另类一区| 飞空精品影院首页| 免费高清在线观看日韩| 国产亚洲午夜精品一区二区久久| 欧美黑人精品巨大| 99香蕉大伊视频| 不卡av一区二区三区| 人妻一区二区av| 2018国产大陆天天弄谢| 日本五十路高清| 国产黄频视频在线观看| 国产成人a∨麻豆精品| 欧美精品高潮呻吟av久久| 精品国产乱码久久久久久男人| 又大又黄又爽视频免费| 国产精品一区二区在线不卡| 国产99久久九九免费精品| 亚洲av片天天在线观看| 久久国产精品影院| 久久久久久免费高清国产稀缺| 伊人久久大香线蕉亚洲五| av有码第一页| 国产成人啪精品午夜网站| 久久久精品区二区三区| 美女脱内裤让男人舔精品视频| 国产精品一区二区在线不卡| 久久久国产欧美日韩av| 日韩大码丰满熟妇| 99久久99久久久精品蜜桃| 久久亚洲国产成人精品v| 最黄视频免费看| 国产片特级美女逼逼视频| 亚洲国产精品999| 99久久精品国产亚洲精品| 亚洲三区欧美一区| 丝袜美腿诱惑在线| 久久99热这里只频精品6学生| 国产又爽黄色视频| 操出白浆在线播放| 亚洲欧美精品综合一区二区三区| 久久久国产精品麻豆| 日韩av免费高清视频| 丝瓜视频免费看黄片| 久久亚洲精品不卡| 在线观看免费午夜福利视频| 久久久久国产精品人妻一区二区| 亚洲av电影在线进入| 国产高清国产精品国产三级| 国产深夜福利视频在线观看| 亚洲五月婷婷丁香| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久av网站| 精品欧美一区二区三区在线| 欧美日韩亚洲高清精品| 免费看十八禁软件| av不卡在线播放| 啦啦啦 在线观看视频| 制服人妻中文乱码| 欧美97在线视频| 性少妇av在线| 天堂俺去俺来也www色官网| 美女大奶头黄色视频| 人妻人人澡人人爽人人| 男女床上黄色一级片免费看| 啦啦啦在线观看免费高清www| 亚洲av综合色区一区| 日韩人妻精品一区2区三区| 高清欧美精品videossex| 亚洲av日韩精品久久久久久密 | 最近手机中文字幕大全| 国产精品.久久久| 黄片播放在线免费| 精品一品国产午夜福利视频| 国产一区亚洲一区在线观看| 精品亚洲成a人片在线观看| 制服诱惑二区| 天堂俺去俺来也www色官网| 欧美精品一区二区大全| 亚洲中文日韩欧美视频| a级片在线免费高清观看视频| 夫妻午夜视频| 国产真人三级小视频在线观看| 日韩大片免费观看网站| 久久精品久久久久久久性| 国产精品 国内视频| 亚洲一区中文字幕在线| 免费女性裸体啪啪无遮挡网站| 亚洲国产欧美网| 亚洲国产日韩一区二区| 久久人人爽av亚洲精品天堂| 赤兔流量卡办理| 欧美另类一区| 成人午夜精彩视频在线观看| 色综合欧美亚洲国产小说| 涩涩av久久男人的天堂| 久久精品久久久久久噜噜老黄| 热re99久久精品国产66热6| 又大又爽又粗| 欧美精品亚洲一区二区| 永久免费av网站大全| 女性被躁到高潮视频| 你懂的网址亚洲精品在线观看| 少妇猛男粗大的猛烈进出视频| 欧美激情极品国产一区二区三区| 高潮久久久久久久久久久不卡| 欧美老熟妇乱子伦牲交| av一本久久久久| 午夜福利一区二区在线看| 午夜久久久在线观看| 日韩精品免费视频一区二区三区| 婷婷丁香在线五月| 男人添女人高潮全过程视频| 又粗又硬又长又爽又黄的视频| videos熟女内射| 国产精品亚洲av一区麻豆| 精品人妻1区二区| 欧美日韩一级在线毛片| 亚洲精品中文字幕在线视频| 这个男人来自地球电影免费观看| 亚洲五月色婷婷综合| 最新的欧美精品一区二区| 亚洲欧美日韩另类电影网站| 久久免费观看电影| 亚洲av电影在线进入| 美女福利国产在线| 久久久国产精品麻豆| 啦啦啦视频在线资源免费观看| 国产欧美亚洲国产| 黄色毛片三级朝国网站| 亚洲,欧美,日韩| 黑人欧美特级aaaaaa片| 高清视频免费观看一区二区| 两个人看的免费小视频| 亚洲男人天堂网一区| 亚洲欧美日韩高清在线视频 | 国产视频首页在线观看| 赤兔流量卡办理| 亚洲av电影在线观看一区二区三区| 欧美老熟妇乱子伦牲交| 黄片播放在线免费| 国产一区二区三区综合在线观看| 黄色怎么调成土黄色| 久久久久精品人妻al黑| 久久国产亚洲av麻豆专区| 黄色一级大片看看| 国产成人a∨麻豆精品| 精品久久久久久久毛片微露脸 | 99香蕉大伊视频| 三上悠亚av全集在线观看| 极品人妻少妇av视频| 久久久国产欧美日韩av| 午夜激情久久久久久久| 精品视频人人做人人爽| 久久久精品94久久精品| 久久亚洲国产成人精品v| 亚洲欧美日韩高清在线视频 | 99九九在线精品视频| 精品国产乱码久久久久久男人| 一本大道久久a久久精品| 欧美在线一区亚洲|