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

    一種用于間斷裝配的流場結(jié)構(gòu)辨識算法

    2019-08-22 01:07:22陳澤棟鄒東陽
    空氣動力學(xué)學(xué)報 2019年3期
    關(guān)鍵詞:方法

    劉 君,陳澤棟,陳 潔,鄒東陽

    (1.大連理工大學(xué) 航空航天學(xué)院,大連 116024;2.大連理工大學(xué) 工程力學(xué)系,大連 116024)

    0 引 言

    流體微團(tuán)加速到超聲速以后無法感知到達(dá)之前遇到的空間物體,只能在物體前通過急劇壓縮的激波才能把速度降下來。因此,激波是可壓縮流動的基本特征。早期計算流體力學(xué)(CFD)的算法無法模擬這一現(xiàn)象,直到1954年Lax提出處理包含有間斷流動問題的弱解理論后[1-2],國內(nèi)外眾多科學(xué)家開始研究在計算過程中能自動分辨出激波的算法,提出TVD、NND、ENO、WENO等具有里程碑意義的激波捕捉格式[3],解決了航天航空等領(lǐng)域很多工程模擬問題,極大提高了CFD地位,使其發(fā)展為一個重要的學(xué)科分支。盡管建立這些格式時考慮了激波間斷,但是本質(zhì)上還是1950年von Neumann提出的把間斷描述成流動參數(shù)連續(xù)變化過渡區(qū)的思路[4],構(gòu)造出來的格式存在以下難以克服的理論困境。

    數(shù)值計算的離散格式是對偏微分方程的近似,按照CFD收斂性理論,數(shù)值解隨著網(wǎng)格尺度減小趨向理論解,CFD追求的極限是方程的理論解,原則上說對于方程不能描述的流動現(xiàn)象,數(shù)值解和實驗值不具有可比性。海平面環(huán)境下馬赫數(shù)小于10的激波厚度大約1×10-8m,流體在這么小的空間內(nèi)受到劇烈壓縮,暫且不說準(zhǔn)確模擬所需要的網(wǎng)格尺度及其帶來的計算量,從連續(xù)介質(zhì)假設(shè)出發(fā)建立的N-S方程是否適用于激波模擬也存在爭議,因為Stoke假設(shè)(3λ+2μ=0)不再成立[5],計算時需要對體積膨脹黏性系數(shù)λ進(jìn)行修正才能得到與實驗相符的激波內(nèi)部溫度[6]。從事稀薄氣體動力學(xué)研究的學(xué)者認(rèn)為激波內(nèi)部結(jié)構(gòu)涉及到熱力學(xué)非平衡,連續(xù)介質(zhì)假設(shè)本身就不成立,“從理論上說,Navier-Stokes方程不能用來描述激波結(jié)構(gòu)的[7]”。無黏流體力學(xué)把激波表示成間斷,兩側(cè)參數(shù)和運動速度滿足蘭金-許貢紐(Rankine-Hugoniot,R-H)關(guān)系式,這是迄今為止幾乎很少有爭議的理論。所謂間斷在數(shù)學(xué)上就是在同一空間位置對應(yīng)至少兩組(激波前后)流體變量,在多維情況下可能更多,例如二維馬赫反射三岔點有四組流體變量。很多基于弱解理論構(gòu)造的激波捕捉格式考慮了間斷,但是用參數(shù)連續(xù)變化的思路描述激波注定不可能計算出真正的數(shù)學(xué)間斷,捕捉到的激波主要是現(xiàn)象模擬。

    按照上述觀點,常規(guī)N-S方程及其理論解不能用來描述激波內(nèi)部物理現(xiàn)象,那么基于離散格式的數(shù)值解至少在過渡區(qū)內(nèi)沒有物理意義,采用Euler方程產(chǎn)生的后果是無法評價過渡區(qū)內(nèi)數(shù)值解的精度(激波前后參數(shù)相差很大),激波捕捉格式對于激波相交等現(xiàn)象更是無能為力。

    20世紀(jì)60年代Moretti等人將激波裝配方法和Lax-Wendroff格式結(jié)合,在CDC 600型計算機上進(jìn)行,總耗時為360 s,給出了超聲速鈍頭體問題的數(shù)值解,2002年Moretti回顧從業(yè)36年來CFD領(lǐng)域模擬激波的進(jìn)展,針對超聲速鈍頭體繞流比較了捕捉法和裝配法計算結(jié)果,認(rèn)為準(zhǔn)確模擬激波還需要發(fā)展裝配法[8-9]。對于50年前提出的裝配法,文獻(xiàn)[10]的評價具有代表性——“激波裝配法具有計算精度高、物理概念清晰等優(yōu)點,但是計算過程復(fù)雜、計算量大,并且必須事先知道激波的大概位置。一般來說,激波裝配法僅適用于那些激波比較簡單而事先可以估計激波形狀和位置的定常流動問題。”

    造成計算復(fù)雜的主要原因是求解過程中激波運動和變形。裝配法減少了網(wǎng)格量,造成計算量大的原因是激波及其相交點等空間網(wǎng)格點需要判別、標(biāo)記或特殊處理,增加代碼編寫復(fù)雜性和并行計算難度。隨著CFD網(wǎng)格技術(shù)發(fā)展,這些難點也得到緩解。國外文獻(xiàn)把最早Moretti提出的稱為邊界激波裝配法 ,僅能處理形狀較為規(guī)則和位移不大的脫體激波。后來Richtmyer等人[11]提出了浮動激波裝配方法,在靜止網(wǎng)格基礎(chǔ)上嵌入R-H關(guān)系式,允許激波在背景網(wǎng)格上滑動,解決了內(nèi)部激波裝配難題。近年來Zhong等人將這種方法和高精度有限差分格式相結(jié)合有效提高內(nèi)部光滑流場的計算精度,成功應(yīng)用于許多高超聲速流動的求解中[12-13]。基于結(jié)構(gòu)網(wǎng)格的邊界激波裝配法和浮動激波裝配法在模擬激波相交和馬赫反射等復(fù)雜波系結(jié)構(gòu)時面臨一定困難。Paciorri和Bonfigolioli等人[14-16]采用非結(jié)構(gòu)網(wǎng)格技術(shù)和局部網(wǎng)格重構(gòu)提出了激波陣面追蹤法,記錄激波在背景網(wǎng)格上的位置,對激波附近的背景網(wǎng)格進(jìn)行挖洞,然后再利用Delanay方法在洞內(nèi)重構(gòu)以激波陣面分界的新網(wǎng)格,時間推進(jìn)過程中根據(jù)R-H關(guān)系式得到新的激波位置,重新判別、挖洞和填補背景網(wǎng)格。這種算法不需要考慮網(wǎng)格節(jié)點之間的連接關(guān)系,也就減少了對于特殊點進(jìn)行處理的麻煩,增加了代碼的普適性,但是在激波邊界節(jié)點對應(yīng)的網(wǎng)格單元只能采用一階插值,像常用的多塊重疊網(wǎng)格技術(shù)一樣存在插值引起的精度降低問題。其次,不允許激波在一個時間步內(nèi)跨越多個網(wǎng)格節(jié)點,通常選擇足夠小的時間步長,因而對計算效率有所影響。上述這3類裝配法時間推進(jìn)過程中需要判別激波位置,增加了計算量,而且激波前后網(wǎng)格數(shù)目是變化的,不利于并行計算。2015年劉君等人[17]提出了基于非結(jié)構(gòu)動網(wǎng)格的激波(間斷)裝配方法,從能夠描述網(wǎng)格變形的任意拉格朗日-歐拉坐標(biāo)系(Arbitrary Lagrangian-Eulerian,ALE)下的控制方程出發(fā),得到激波位置及其兩側(cè)的參數(shù)后,裝配的激波就變成了一個移動邊界,對于存在內(nèi)部激波的計算域則變成類似于拼接網(wǎng)格的多個分區(qū),激波運動由上下游參數(shù)驅(qū)動,采用成熟的動網(wǎng)格技術(shù)[18]進(jìn)行移動和變形。與Paciorri及Bonfigolioli的裝配法相比,在保留非結(jié)構(gòu)網(wǎng)格的優(yōu)點外,新的裝配法不需要在每一步對計算網(wǎng)格進(jìn)行挖洞重構(gòu),避免了插值誤差,易于并行計算,并且時間步長不受限制。采用這種算法模擬了內(nèi)部存在正規(guī)反射和馬赫反射等復(fù)雜激波結(jié)構(gòu)的算例[19-21],從效果看,較好地解決了傳統(tǒng)裝配法遇到的“計算過程復(fù)雜、計算量大”問題。

    為了突破“激波裝配法僅適用于那些激波比較簡單而事先可以估計激波形狀和位置的定常流動問題”的限制,劉君等人提出了在采用捕捉法計算的流場基礎(chǔ)上辨識并裝配主要激波結(jié)構(gòu)、計算過程中自動裝配內(nèi)部細(xì)致結(jié)構(gòu)和新生激波的解決思路。在技術(shù)途徑上提出了新的數(shù)據(jù)結(jié)構(gòu),在激波邊界通量計算時,如果采用Roe、van Leer、AUSM等格式就是捕捉法,如果采用R-H關(guān)系式就是裝配法,實現(xiàn)兩種方法之間任意轉(zhuǎn)換。按照上述思路,為了解決裝配法的應(yīng)用難題,亟需發(fā)展激波、接觸間斷等流場特征結(jié)構(gòu)的辨識算法。

    目前,常用的激波辨識算法主要有三種[22]:第一種根據(jù)最大變量梯度確定激波面,該方法實現(xiàn)簡單,但需要設(shè)置合適的過濾器才能消除錯誤的結(jié)果;第二種為法向馬赫數(shù)法[24],該方法認(rèn)為在激波面上流動的法向馬赫數(shù)等于1,該方法與第一種方法相比適用性更強,但也需要設(shè)置適當(dāng)?shù)倪^濾器才能獲得良好結(jié)果,且難以用于探測接觸間斷;第三種為基于特征線理論的方法[25-26],該方法與前兩種方法相比更精確、更有效,但實施過程復(fù)雜且計算量大。

    綜上所述,這些激波辨識算法難以在復(fù)雜三維問題中廣泛應(yīng)用,且上述算法辨識出的激波通常為具有一定寬度的激波帶,難以直接用于裝配法。因此,需要發(fā)展新的辨識算法來推動激波裝配法實用化進(jìn)程。

    本文提出了一種光源射線法,該方法可將捕捉法中探測到的具有一定寬度的激波帶擬合成一系列離散線段(二維)或三角形面片(三維),利用曲線/曲面擬合算法進(jìn)一步將其擬合成連續(xù)的曲線/曲面。該曲線/曲面可作為激波裝配法中的激波面,并且可根據(jù)計算需要對其進(jìn)行任意形式的網(wǎng)格劃分。為便于表述,先用二維模型的實際操作過程進(jìn)行原理說明,然后簡單介紹三維推廣過程,并展示實際工程問題的辨識結(jié)果。

    1 流場結(jié)構(gòu)辨識算法

    1.1 二維辨識算法流程

    (a)捕捉法流場 (b)激波單元 (c)激波曲線

    (d)分片擬合

    (e)確定特征點圖1 光源射線法示意圖Fig.1 Sketch of source-rays method

    1.2 脫體激波辨識

    針對二維超聲速鈍頭體繞流問題,采用1.1節(jié)所述方法進(jìn)行激波面辨識。算例中,空氣來流馬赫數(shù)Ma=5.0、無量綱密度ρ=1、無量綱溫度T=1,來流攻角為0°。鈍頭體中心位置坐標(biāo)為(0,0),鈍頭體半徑為0.0254。網(wǎng)格為160×100的四邊形網(wǎng)格,其中沿鈍頭體法向的節(jié)點數(shù)為100。計算網(wǎng)格和捕捉法流場壓力云圖分別如圖2(a)和圖2(b)所示。

    (a)網(wǎng)格 (b)壓力云圖 (c)激波點集 (d)激波面圖2 基于壓力梯度的激波辨識結(jié)果Fig.2 Shock wave surface extraction based on pressure gradient

    考慮到很多高精度格式的限制器選擇壓力梯度來進(jìn)行調(diào)節(jié),在此也將其作為激波特征參數(shù)來探測脫體激波。對于捕捉法得到的壓力場,用格林公式計算出每個網(wǎng)格單元的壓力梯度矢量。作為間斷面,激波厚度理論上處于分子自由程大小這一數(shù)量級,其兩側(cè)壓力的梯度值為無窮大。但是,在捕捉法數(shù)值計算中,激波表現(xiàn)為帶狀分布的過渡區(qū)且具有較大壓力梯度。對于不同的問題以及不同的計算網(wǎng)格,激波帶上的壓力梯度值都有所不同。因此,難以通過一個固定的壓力梯度值將激波帶提取出來。本文為了準(zhǔn)確獲得激波點集的分布,首先計算出流場中的最大壓力梯度值|p|max,然后針對不同問題設(shè)置相應(yīng)的閾值系數(shù)k,壓力梯度閾值|p|thre=k|p|max,所有滿足條件|p|≥|p|thre的單元格心點構(gòu)成激波點集T。Case1取閾值系數(shù)k=0.1,獲得離散點集,如圖2(c)所示。針對上述激波離散點集,采用光源射線法進(jìn)行曲線擬合。光源S位置固定在(x,y)=(0.1,0)處,相鄰兩條射線間的夾角取常值φi=4.5°。辨識出的結(jié)果如圖2(d)所示,從圖中可以看出辨識出的激波曲線與捕捉法流場中的激波點集T吻合度很高,能準(zhǔn)確描述激波形狀,對激波曲線重新進(jìn)行網(wǎng)格離散即可作為激波邊界直接用于激波裝配法。

    離散點集的分布取決于所選擇參數(shù)和閾值,為了比較不同壓力梯度閾值下所辨識出的激波曲線,在捕捉法流場的基礎(chǔ)上另外補充兩個辨識狀態(tài):Case2取k=0.2,Case3取k=0.3。得到的點集及其最終辨識出來的脫體激波曲線如圖3所示。從3個狀態(tài)的辨識結(jié)果可以看出,不同的壓力梯度條件下所獲取的離散點集有所不同,導(dǎo)致最終辨識出的激波曲線長度有所區(qū)別,但三條激波曲線都能與捕捉法流場中激波的位置和形狀保持一致,將其作為激波裝配法中的激波邊界可快速得到收斂解。得到初始激波位置以后的裝配法計算方法和過程參見文獻(xiàn)[17,19-21],這里不再贅述。

    (a)Case1 (b)Case2 (c)Case3

    圖3 不同壓力梯度閾值下的脫體激波辨識結(jié)果比較
    Fig.3 Comparison of shock surface extraction results under different pressure gradient thresholds

    1.3 內(nèi)部激波及接觸間斷辨識

    對波系更為復(fù)雜的馬赫反射問題進(jìn)行數(shù)值模擬,計算區(qū)域及邊界條件如圖4所示。邊長為LAF=0.4、LAB=0.1、LDE=0.312、LEF=0.85,斜楔BC傾角為14°;邊界AF為超聲速入口、DE為超聲速出口,其余為滑移壁。計算網(wǎng)格為300×218的四邊形網(wǎng)格,其中入口和出口邊界上的網(wǎng)格節(jié)點數(shù)為300。來流馬赫數(shù)Ma=1.9,攻角為0°。

    圖4 計算域Fig.4 Calculation domain

    由于斜楔BC的偏轉(zhuǎn)作用,自由來流在點B處產(chǎn)生一道入射激波(Shock1);激波達(dá)到壁面FE處時發(fā)生馬赫反射,產(chǎn)生一道反射激波(Shock2)和一條馬赫桿(Shock3);反射激波在壁面CD處再次發(fā)生反射,產(chǎn)生第二道反射激波(Shock4)。流場穩(wěn)定后的密度云圖以及波系結(jié)構(gòu)如圖5所示。采用激波及接觸間斷探測方法可從流場中獲得主要的四條激波點集和一條接觸間斷點集分布。為了檢驗不同的探測方法對最終結(jié)果的影響,激波點集采用密度梯度作為參考量進(jìn)行探測,接觸間斷點集由密度梯度和壓力梯度共同決定。同樣,先計算出全場最大密度梯度值|ρ|max和最大壓力梯度|p|max。激波探測的密度梯度閾值所有滿足的單元節(jié)點構(gòu)成激波點集,這里取接觸間斷單元由密度梯度和壓力梯度共同決定,滿足密度梯度較大且壓力梯度較小,即其中所有滿足該條件的單元節(jié)點構(gòu)成曲線擬合所需的接觸間斷點集。這里取其中在不同密度梯度閾值下辨識出的接觸間斷長度有一定區(qū)別,但位置和形狀基本一致。激波和接觸間斷離散點集分布如圖5所示。

    獲得離散點集后,通過在區(qū)域內(nèi)合理布置光源可擬合出每條激波和接觸間斷曲線。首先在Ⅱ區(qū)設(shè)置一個光源,光源發(fā)出等間距射線,采用上文光源射線法可辨識出表征激波1和激波2的曲線,如圖6(a)所示。同理,在Ⅲ區(qū)布置光源可辨識出激波4和接觸間斷曲線,如圖6(b)所示。在I區(qū)設(shè)置光源可辨識出馬赫桿曲線,如圖6(c)所示。所有曲線的最終辨識結(jié)果與流場密度云圖的比較如圖7所示。可以看出,辨識出的曲線與采用捕捉法獲得的波系結(jié)構(gòu)高度吻合,直接用于裝配法有助于快速得到收斂解。得到初始激波位置以后的裝配法計算方法和過程參見文獻(xiàn)[17,19-21],后文也會有簡單介紹,這里不再贅述。

    圖5 捕捉法密度云圖及間斷點集分布Fig.5 Density contours and discontinuity point set distribution under shock-capturing schemes

    (a)激波1-2 (b)激波4和接觸間斷 (c)激波3

    圖6 激波和接觸間斷辨識過程
    Fig.6 Shock wave and contact discontinuity surface extraction process

    圖7 激波及接觸間斷辨識結(jié)果與流場比較Fig.7 Comparison between feature surface extraction result and flow field structure

    值得注意的是,這里光源的位置分布并不唯一,可采用不同的分布方式達(dá)到相似的效果,比如可先在Ⅰ區(qū)設(shè)置光源同時擬合出激波1和激波3曲線,再在Ⅲ區(qū)設(shè)置光源得到激波2、激波4和接觸間斷曲線。

    2 三維實現(xiàn)

    2.1 三維辨識算法流程

    本文通過引入單位球可確定各射線之間的連接關(guān)系,如圖9所示。對半徑為1的球面進(jìn)行三角形離散,獲得一系列三角形面片,各頂點連接關(guān)系確定,且三角形面片滿足連續(xù)且不相交條件。

    圖8 射線方向示意圖

    Fig.8 Ray direction

    圖9 單位球

    Fig.9 Unit sphere

    將單位球的球心作為光源S,球心與單位球上三角形面片的三個頂點ai、aj、ak的連線li、lj、lk構(gòu)成一個子區(qū)域Ap,如圖10(a)所示。若Ap內(nèi)離散點的個數(shù)不小于3,則對Ap內(nèi)所有離散點進(jìn)行平面擬合,擬合出的平面與三條射線的交點分別為bi,p、bj,p和bk,p;若該子區(qū)域內(nèi)離散點的個數(shù)小于3,則跳過該區(qū)域。所有區(qū)域擬合完成后,每條射線li與不同子區(qū)域內(nèi)的擬合平面形成若干交點。將各射線li上的交點取平均值獲得特征點bi,則各子區(qū)域Ap的三條射線邊界上的特征點bi、bj和bk構(gòu)成目標(biāo)三角形面片,如圖10(b)所示。最終,整個離散點集T被擬合成一系列連續(xù)的三角形面片,三角形面片頂點的連接關(guān)系由單位球確定。可對上述三角形面片構(gòu)成的空間曲面重新進(jìn)行網(wǎng)格劃分,獲得可用于裝配法的激波面網(wǎng)格。

    (a)分區(qū)擬合 (b)三角面片圖10 子區(qū)域內(nèi)三角形面片擬合過程Fig.10 Triangular fitting process in sub-region

    2.2 三維測試算例

    通過以下兩個算例驗證本文辨識算法在三維問題中的有效性。

    在橢球面附近隨機生成100 000個點作為離散點集,隨機點坐標(biāo)B(x′,y′,z′)滿足:

    其中,

    式中,i、j、k、l、m為區(qū)間[0,1]內(nèi)的隨機數(shù)。

    獲得的離散點集如圖11所示。單位球共有736個節(jié)點、1468個三角形單元(如圖9),球心位置坐標(biāo)為(0,0,0)。單位球位置及橢球面辨識結(jié)果如圖12所示。從結(jié)果可以看出,辨識出的曲面光滑完整且形狀與離散點分布形狀一致。

    圖11 橢球隨機點集圖12 單位球及擬合結(jié)果Fig.11 Random point set of ellipsoidFig.12 Unit sphere and the extracted surface of ellipsoid

    圖13 雙橢球隨機點集圖14 雙橢球面擬合結(jié)果Fig.13 Random point set of double-ellipsoidFig.14 Extracted surface of double-ellipsoid

    2.3 三維繞流激波辨識

    對三維超聲速球頭繞流進(jìn)行模擬,來流馬赫數(shù)Ma=5.0,無量綱密度ρ=1.0,無量綱溫度T=1.0,來流攻角和側(cè)滑角均為0°。球頭的球心位置坐標(biāo)為(0,0,0),球頭半徑為0.0254。計算網(wǎng)格共有508 365個的六面體單元,球面網(wǎng)格如圖15(a)所示。圖15(b)為流場壓力云圖,圖15(c)為根據(jù)流場壓力梯度獲得的激波帶離散點集。算例中,壓力梯度閾值系數(shù)k=0.05離散點為滿足|p|≥|p|thre的所有單元節(jié)點,共有23 640個。單位球共包含1734個節(jié)點,3464個三角形面片,球心坐標(biāo)為(2,0,0)。最終辨識出的激波面如圖15(d)所示,除了邊緣附近外的其他區(qū)域辨識結(jié)果與離散點集吻合較好,將邊緣位置稍作處理即可作為初始激波面用于激波裝配法。對于不同的初始位置的激波,最終收斂解相同,故本例不再展示裝配法的計算過程,計算方法參見文獻(xiàn)[17,19-21]。

    (a)球頭表面網(wǎng)格 (b)捕捉法壓力云圖

    (c)激波點集 (d)激波面圖15 三維超聲速球頭繞流激波辨識過程Fig.15 Shock wave surface extraction process of 3D blunt body flow

    在球頭繞流算例的基礎(chǔ)上,對波系較為復(fù)雜的球柱錐組合體外形進(jìn)行模擬,檢驗特征面提取方法處理三維復(fù)雜問題的能力。計算模型如圖16(a)所示,來流馬赫數(shù)Ma=4.04,來流攻角為20°。波系結(jié)構(gòu)如圖16(b)所示,在最外層形成一道弓形激波,在錐裙與柱體的結(jié)合處形成一個內(nèi)嵌激波。探測出的激波點集和提取出的激波面分別如圖16(c)和圖16(d)所示。

    (a)計算模型 (b)捕捉法壓力云圖

    (c)激波點集 (d)激波面圖16 球柱錐組合體繞流激波辨識過程Fig.16 Shock wave surface extraction process of cylinder with a hemispherical nose and a conical flow

    3 裝配法實例

    利用上述方法獲得初始激波面后,將其作為激波邊界用于激波裝配法,進(jìn)行流動數(shù)值模擬,考查所述激波辨識方法在實際使用中的效果。

    對于定常激波裝配法,初始激波位置及形狀會影響計算的收斂過程,但最終的收斂結(jié)果都相同。因此,本節(jié)直接采用文獻(xiàn)[19]算例的部分結(jié)果說明激波裝配法與捕捉法的區(qū)別。該算例為具有復(fù)雜波系結(jié)構(gòu)的馬赫反射問題,來流馬赫數(shù)2.0,偏轉(zhuǎn)角δ=14°。計算區(qū)域為L×L的正方形,邊界上網(wǎng)格節(jié)點均勻分布,節(jié)點間距為0.01L,內(nèi)部區(qū)域為三角形網(wǎng)格,單元數(shù)為22 429。流場波系結(jié)構(gòu)及邊界條件如圖17所示,入射激波IS和反射激波RS、馬赫桿MS交于三波點TP,SS為接觸間斷。

    首先采用激波捕捉法獲得初始流場,在此基礎(chǔ)上采用本文辨識算法辨識出激波及接觸間斷的位置和形狀,將其作為裝配法的初始間斷面,結(jié)合R-H關(guān)系式進(jìn)行裝配法流場計算。得到初始激波位置以后的計算過程參見文獻(xiàn)[19],計算結(jié)果完全一致。下面給出部分結(jié)果,如圖18所示,其中圖18(a)為采用捕捉法的計算結(jié)果,圖18(b)為采用裝配法的計算結(jié)果。從圖中可以看出,捕捉法流場中激波附近出現(xiàn)等值線的扭曲和震蕩,且激波呈帶狀分布,具有很寬的過渡區(qū)。與此相比,采用裝配法所獲得的流場中各區(qū)域之間界線分明,激波是一個沒有寬度的界面,激波兩側(cè)變量值嚴(yán)格滿足R-H關(guān)系式。

    圖17 邊界條件及流場結(jié)構(gòu)示意圖[19]Fig.17 Boundary conditions and sketch of flow field structure[19]

    (a)捕捉法計算結(jié)果

    (b)裝配法計算結(jié)果圖18 裝配法與捕捉法計算結(jié)果對比(馬赫數(shù)云圖)[19]Fig.18 Comparison between shock-capturing solution and shock-fitting solution[19]

    4 結(jié) 論

    通常對于具有一定“厚度”的空間點集,采用傳統(tǒng)的擬合方法難以將其擬合成形狀一致的光滑曲面,本文提出的光源射線法很好地解決了此類問題。

    對于從捕捉法流場中獲取的具有一定“厚度”的激波帶和接觸間斷帶,采用該方法進(jìn)行空間曲面擬合后可作為激波裝配法中的初始間斷面。尤其對于具有復(fù)雜波系結(jié)構(gòu)的流場,該方法提供了一種裝配初始流場的途徑。本文詳細(xì)介紹了該算法在二維問題中的實現(xiàn)過程,并將其成功推廣應(yīng)用于三維問題。算例結(jié)果表明,辨識出的激波曲面整體上與捕捉法流場中激波帶分布相吻合,但是三維問題中激波曲面邊緣附近的辨識效果有待進(jìn)一步提高。對于復(fù)雜工程問題,往往存在較為復(fù)雜的三維波系結(jié)構(gòu),間斷面的提取過程更為困難和繁瑣,該方法還需進(jìn)一步完善。

    對于傳統(tǒng)激波裝配法,在初始激波面給定后需要經(jīng)過若干步迭代才能收斂到準(zhǔn)確的形狀和位置。尤其對于三維問題,收斂過程更加漫長,甚至出現(xiàn)難以收斂的情況。采用本文方法可提供較為準(zhǔn)確的初始激波,有助于縮短收斂過程。

    本文所提方法除了用于激波裝配法外,未來還可應(yīng)用于自適應(yīng)網(wǎng)格算法。根據(jù)不同的流動特征可擬合出相應(yīng)的空間曲面,進(jìn)而可以任意調(diào)整該曲面附近區(qū)域的網(wǎng)格分布,達(dá)到網(wǎng)格根據(jù)流動特征自動調(diào)整的目的。

    猜你喜歡
    方法
    中醫(yī)特有的急救方法
    中老年保健(2021年9期)2021-08-24 03:52:04
    高中數(shù)學(xué)教學(xué)改革的方法
    河北畫報(2021年2期)2021-05-25 02:07:46
    化學(xué)反應(yīng)多變幻 “虛擬”方法幫大忙
    變快的方法
    兒童繪本(2020年5期)2020-04-07 17:46:30
    學(xué)習(xí)方法
    可能是方法不對
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    最有效的簡單方法
    山東青年(2016年1期)2016-02-28 14:25:23
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    久久精品国产亚洲网站| 亚洲欧美日韩另类电影网站 | 综合色丁香网| 男人爽女人下面视频在线观看| 久久精品国产a三级三级三级| 免费看不卡的av| 成年女人看的毛片在线观看| 一级毛片 在线播放| 男男h啪啪无遮挡| 一级av片app| 国产探花在线观看一区二区| 国产免费一级a男人的天堂| 免费av毛片视频| 色播亚洲综合网| 婷婷色av中文字幕| 中文资源天堂在线| 日本色播在线视频| 久久97久久精品| 色视频www国产| 国产成人午夜福利电影在线观看| 男的添女的下面高潮视频| www.色视频.com| 久久6这里有精品| 国产精品伦人一区二区| 天天一区二区日本电影三级| 日韩精品有码人妻一区| 欧美成人一区二区免费高清观看| 99久久人妻综合| 岛国毛片在线播放| 日韩精品有码人妻一区| 日本一二三区视频观看| 男女那种视频在线观看| 最近中文字幕2019免费版| 又大又黄又爽视频免费| 能在线免费看毛片的网站| 久久精品国产鲁丝片午夜精品| 白带黄色成豆腐渣| 高清日韩中文字幕在线| 亚洲四区av| 老女人水多毛片| 三级男女做爰猛烈吃奶摸视频| 久久久久精品久久久久真实原创| 高清av免费在线| 精品午夜福利在线看| 中文字幕亚洲精品专区| 又大又黄又爽视频免费| 欧美精品国产亚洲| 欧美一区二区亚洲| 一级毛片我不卡| 男人爽女人下面视频在线观看| 婷婷色av中文字幕| tube8黄色片| 亚洲精品一二三| 亚洲成色77777| 午夜福利网站1000一区二区三区| 国产午夜精品久久久久久一区二区三区| www.色视频.com| 国产综合精华液| 少妇熟女欧美另类| 亚洲无线观看免费| 国内揄拍国产精品人妻在线| 好男人视频免费观看在线| 亚洲av不卡在线观看| 免费看光身美女| 日本午夜av视频| 69人妻影院| 亚洲国产成人一精品久久久| 女人久久www免费人成看片| 少妇被粗大猛烈的视频| 一本久久精品| 日韩欧美一区视频在线观看 | av专区在线播放| 久久热精品热| 亚洲怡红院男人天堂| 欧美精品国产亚洲| 97在线视频观看| 在线观看美女被高潮喷水网站| 插阴视频在线观看视频| 成人鲁丝片一二三区免费| 夜夜爽夜夜爽视频| 色综合色国产| 亚洲激情五月婷婷啪啪| 深夜a级毛片| 国产男女内射视频| 97热精品久久久久久| 国产亚洲一区二区精品| 国产色婷婷99| 免费av观看视频| 免费av不卡在线播放| 天堂中文最新版在线下载 | 亚洲美女视频黄频| 看非洲黑人一级黄片| 五月伊人婷婷丁香| 一边亲一边摸免费视频| 黄色怎么调成土黄色| 国产精品熟女久久久久浪| 亚洲国产成人一精品久久久| 国产在线男女| 99视频精品全部免费 在线| 亚洲欧美精品专区久久| 国产男人的电影天堂91| 免费观看性生交大片5| 精品人妻熟女av久视频| 99久久人妻综合| 不卡视频在线观看欧美| 最新中文字幕久久久久| 国语对白做爰xxxⅹ性视频网站| 观看免费一级毛片| 国产精品国产三级专区第一集| 美女xxoo啪啪120秒动态图| 男人狂女人下面高潮的视频| 最近中文字幕2019免费版| 大码成人一级视频| 干丝袜人妻中文字幕| 国产毛片a区久久久久| 日本wwww免费看| 毛片女人毛片| 亚洲无线观看免费| 高清日韩中文字幕在线| 五月天丁香电影| 亚洲av中文字字幕乱码综合| 深夜a级毛片| 美女内射精品一级片tv| 精品视频人人做人人爽| 深夜a级毛片| 六月丁香七月| 国产在线男女| 菩萨蛮人人尽说江南好唐韦庄| 男插女下体视频免费在线播放| 日韩电影二区| 男女下面进入的视频免费午夜| 欧美xxxx黑人xx丫x性爽| videos熟女内射| 永久免费av网站大全| 亚洲最大成人av| 特大巨黑吊av在线直播| 久久久久久久久久人人人人人人| 国国产精品蜜臀av免费| 免费黄网站久久成人精品| 亚洲电影在线观看av| 一级a做视频免费观看| 啦啦啦啦在线视频资源| 亚洲成色77777| 国产成人精品婷婷| 欧美成人精品欧美一级黄| 国产男女超爽视频在线观看| 777米奇影视久久| 秋霞在线观看毛片| 久久鲁丝午夜福利片| 国产精品久久久久久精品电影| 国产亚洲精品久久久com| 国产午夜精品一二区理论片| 99热网站在线观看| 国产人妻一区二区三区在| 尾随美女入室| 亚洲精品色激情综合| 国产精品99久久99久久久不卡 | 成人无遮挡网站| 亚洲丝袜综合中文字幕| 最新中文字幕久久久久| 丝袜脚勾引网站| 人人妻人人爽人人添夜夜欢视频 | 三级经典国产精品| 97超视频在线观看视频| 青春草国产在线视频| 成年版毛片免费区| 成人特级av手机在线观看| 欧美日韩亚洲高清精品| 91久久精品电影网| 欧美xxxx黑人xx丫x性爽| 日韩欧美一区视频在线观看 | 国产黄色免费在线视频| 丰满少妇做爰视频| 久久亚洲国产成人精品v| 日韩不卡一区二区三区视频在线| 国产中年淑女户外野战色| 亚洲最大成人中文| 国产成人福利小说| av线在线观看网站| 自拍偷自拍亚洲精品老妇| 亚洲婷婷狠狠爱综合网| av专区在线播放| 国产一区二区在线观看日韩| 国产老妇伦熟女老妇高清| 国产精品人妻久久久久久| 国产精品秋霞免费鲁丝片| 久久久久网色| 成人亚洲欧美一区二区av| 九色成人免费人妻av| 久久99热6这里只有精品| 亚洲成色77777| 欧美国产精品一级二级三级 | 欧美日韩精品成人综合77777| 午夜福利视频1000在线观看| 日本熟妇午夜| 丰满乱子伦码专区| av线在线观看网站| 久久国产乱子免费精品| 国产 精品1| 中文资源天堂在线| 国产成人福利小说| 久久久久久久久久久丰满| 最近中文字幕2019免费版| 亚洲精品国产成人久久av| 亚洲久久久久久中文字幕| 欧美精品一区二区大全| 欧美日韩国产mv在线观看视频 | 边亲边吃奶的免费视频| 国产免费视频播放在线视频| 欧美一级a爱片免费观看看| 国产精品成人在线| 亚洲综合精品二区| 亚洲av成人精品一区久久| 天堂中文最新版在线下载 | 成人亚洲精品av一区二区| 成人一区二区视频在线观看| 男人舔奶头视频| 超碰av人人做人人爽久久| 熟女人妻精品中文字幕| 国产精品精品国产色婷婷| 观看美女的网站| 99久久人妻综合| 22中文网久久字幕| 超碰av人人做人人爽久久| 高清毛片免费看| 制服丝袜香蕉在线| 国产综合精华液| 亚洲精品一区蜜桃| 在线看a的网站| 视频中文字幕在线观看| 丝袜脚勾引网站| 黄色欧美视频在线观看| 男女啪啪激烈高潮av片| 噜噜噜噜噜久久久久久91| 欧美3d第一页| 国产毛片a区久久久久| 一区二区av电影网| 国产精品久久久久久av不卡| 国产久久久一区二区三区| 免费av不卡在线播放| www.色视频.com| 欧美xxxx黑人xx丫x性爽| 老女人水多毛片| 永久免费av网站大全| 精品一区在线观看国产| 国产视频内射| 久久精品国产亚洲av涩爱| av播播在线观看一区| 亚洲精品久久午夜乱码| 日本三级黄在线观看| 夜夜看夜夜爽夜夜摸| 国产成人福利小说| a级毛色黄片| 人妻制服诱惑在线中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 男插女下体视频免费在线播放| 精品视频人人做人人爽| 国产亚洲5aaaaa淫片| 校园人妻丝袜中文字幕| 日本三级黄在线观看| 亚洲欧美精品自产自拍| 亚洲精品国产色婷婷电影| 久久精品久久久久久久性| 伦精品一区二区三区| 亚洲熟女精品中文字幕| 一二三四中文在线观看免费高清| 国产大屁股一区二区在线视频| 两个人的视频大全免费| 麻豆国产97在线/欧美| 亚洲欧美日韩另类电影网站 | 男女国产视频网站| 国内精品宾馆在线| 一级毛片aaaaaa免费看小| 国产极品天堂在线| 日韩成人av中文字幕在线观看| 欧美高清性xxxxhd video| av.在线天堂| 久久久久久国产a免费观看| 麻豆精品久久久久久蜜桃| 岛国毛片在线播放| 日日撸夜夜添| 精品人妻熟女av久视频| 伦精品一区二区三区| 亚洲av.av天堂| 亚洲欧美日韩另类电影网站 | 少妇人妻精品综合一区二区| a级一级毛片免费在线观看| 黄片无遮挡物在线观看| 校园人妻丝袜中文字幕| 国产日韩欧美在线精品| 黄片无遮挡物在线观看| 午夜福利视频1000在线观看| 欧美潮喷喷水| 午夜福利在线在线| 久久综合国产亚洲精品| 蜜桃亚洲精品一区二区三区| 亚洲成色77777| 高清av免费在线| 久久久久性生活片| 亚洲久久久久久中文字幕| 欧美97在线视频| 亚洲av在线观看美女高潮| 亚洲经典国产精华液单| 最近最新中文字幕免费大全7| 免费黄色在线免费观看| www.av在线官网国产| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最近手机中文字幕大全| 国产精品不卡视频一区二区| 精品久久久久久电影网| 国产片特级美女逼逼视频| 国产成人a区在线观看| 黄色一级大片看看| 97超视频在线观看视频| 国产成人a∨麻豆精品| 精品人妻一区二区三区麻豆| 亚洲欧美日韩东京热| 亚洲国产欧美在线一区| 亚洲不卡免费看| 久久久久久九九精品二区国产| 91在线精品国自产拍蜜月| 国产高清三级在线| 免费高清在线观看视频在线观看| 国产高清国产精品国产三级 | 亚洲高清免费不卡视频| 成人综合一区亚洲| 亚洲一区二区三区欧美精品 | 成人高潮视频无遮挡免费网站| 99热网站在线观看| 国产乱来视频区| 欧美日韩综合久久久久久| 精品久久国产蜜桃| 亚洲国产精品国产精品| 久久久午夜欧美精品| 一本久久精品| 另类亚洲欧美激情| 99热网站在线观看| 建设人人有责人人尽责人人享有的 | 亚洲精品日韩在线中文字幕| 一本久久精品| 99热这里只有是精品50| 少妇被粗大猛烈的视频| 男人爽女人下面视频在线观看| 亚洲色图av天堂| 日本猛色少妇xxxxx猛交久久| 欧美极品一区二区三区四区| 日韩一本色道免费dvd| 青春草亚洲视频在线观看| 极品少妇高潮喷水抽搐| 五月伊人婷婷丁香| 日韩欧美精品免费久久| 亚洲在久久综合| 国产精品三级大全| 亚洲不卡免费看| av播播在线观看一区| 水蜜桃什么品种好| 伊人久久精品亚洲午夜| 噜噜噜噜噜久久久久久91| 免费黄频网站在线观看国产| 热99国产精品久久久久久7| 亚洲欧美成人精品一区二区| 美女国产视频在线观看| 成人亚洲精品av一区二区| 久久99精品国语久久久| 国产精品99久久99久久久不卡 | 久久午夜福利片| 亚洲欧洲日产国产| 免费黄色在线免费观看| 国产老妇女一区| 日韩强制内射视频| 国产高清有码在线观看视频| 亚洲婷婷狠狠爱综合网| 亚洲精品乱码久久久v下载方式| 在线观看av片永久免费下载| 一二三四中文在线观看免费高清| 亚洲av中文字字幕乱码综合| 国产精品一及| 国产乱来视频区| 黄色日韩在线| 少妇人妻久久综合中文| 国产美女午夜福利| 五月开心婷婷网| 成人特级av手机在线观看| 别揉我奶头 嗯啊视频| 国产极品天堂在线| 男插女下体视频免费在线播放| 国产成人freesex在线| 男女边吃奶边做爰视频| 亚洲精品一区蜜桃| 午夜福利在线在线| 韩国高清视频一区二区三区| 国产 一区精品| 嘟嘟电影网在线观看| 日本午夜av视频| 亚洲精品国产av蜜桃| 一区二区av电影网| 日本一本二区三区精品| 亚洲国产欧美人成| 亚洲国产精品999| 国产视频首页在线观看| 99re6热这里在线精品视频| 嫩草影院新地址| 国产毛片a区久久久久| 美女内射精品一级片tv| 制服丝袜香蕉在线| 少妇人妻精品综合一区二区| 日本色播在线视频| 国产91av在线免费观看| 97超视频在线观看视频| 80岁老熟妇乱子伦牲交| 欧美日韩视频高清一区二区三区二| 亚洲国产日韩一区二区| 18禁裸乳无遮挡免费网站照片| 精品熟女少妇av免费看| 97精品久久久久久久久久精品| 只有这里有精品99| 成年av动漫网址| 搡女人真爽免费视频火全软件| 一级黄片播放器| 日韩欧美 国产精品| 美女xxoo啪啪120秒动态图| 精品午夜福利在线看| 99热全是精品| av在线蜜桃| 大片免费播放器 马上看| 国产精品99久久久久久久久| 国产精品爽爽va在线观看网站| 久久精品熟女亚洲av麻豆精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久亚洲中文字幕| 亚洲国产av新网站| 欧美极品一区二区三区四区| 国产一区二区三区av在线| 久久精品国产亚洲网站| 亚洲av欧美aⅴ国产| 亚洲最大成人手机在线| 免费av不卡在线播放| 亚洲,一卡二卡三卡| 伦精品一区二区三区| 听说在线观看完整版免费高清| 看黄色毛片网站| 日韩av免费高清视频| 亚洲av成人精品一区久久| 99久久中文字幕三级久久日本| 亚洲久久久久久中文字幕| 精华霜和精华液先用哪个| 国产69精品久久久久777片| 亚洲精品视频女| 亚洲人与动物交配视频| 久久热精品热| 国产精品一区www在线观看| 欧美xxⅹ黑人| 国产亚洲av嫩草精品影院| 在线天堂最新版资源| 国产av国产精品国产| 亚洲精品一二三| 精品久久久久久久久av| 免费少妇av软件| 日日撸夜夜添| 久久人人爽人人片av| 国产免费一级a男人的天堂| 毛片女人毛片| 欧美日韩一区二区视频在线观看视频在线 | 成年女人在线观看亚洲视频 | 亚洲av不卡在线观看| 国产亚洲5aaaaa淫片| 亚洲国产精品成人久久小说| 制服丝袜香蕉在线| 婷婷色综合www| 女人久久www免费人成看片| 日本av手机在线免费观看| 国产精品不卡视频一区二区| 国产一区二区亚洲精品在线观看| 精品少妇黑人巨大在线播放| 91久久精品国产一区二区三区| 日本熟妇午夜| 一区二区三区免费毛片| 中文天堂在线官网| 一级av片app| 七月丁香在线播放| 午夜免费男女啪啪视频观看| 丝袜美腿在线中文| 麻豆久久精品国产亚洲av| 天堂俺去俺来也www色官网| 日本wwww免费看| 久久精品熟女亚洲av麻豆精品| 国产综合精华液| 欧美区成人在线视频| 亚洲成人久久爱视频| 日韩欧美精品v在线| 精品人妻一区二区三区麻豆| 男女国产视频网站| 亚洲av二区三区四区| 简卡轻食公司| 精品国产乱码久久久久久小说| 国产高清三级在线| 偷拍熟女少妇极品色| 亚洲av二区三区四区| 午夜爱爱视频在线播放| 少妇被粗大猛烈的视频| 联通29元200g的流量卡| 日本爱情动作片www.在线观看| 国产精品久久久久久久电影| 99热这里只有是精品在线观看| 草草在线视频免费看| 国产探花在线观看一区二区| 久久精品久久精品一区二区三区| 日本-黄色视频高清免费观看| 99热这里只有是精品在线观看| 菩萨蛮人人尽说江南好唐韦庄| 久热这里只有精品99| 久久ye,这里只有精品| 亚洲精华国产精华液的使用体验| 中国国产av一级| 麻豆乱淫一区二区| 亚洲在线观看片| 亚洲精品色激情综合| 日本一本二区三区精品| 亚洲国产精品专区欧美| 成人亚洲精品av一区二区| eeuss影院久久| 国产爽快片一区二区三区| 一区二区av电影网| 别揉我奶头 嗯啊视频| 精品一区二区免费观看| 国产在线一区二区三区精| 精品一区二区免费观看| 大又大粗又爽又黄少妇毛片口| 国产精品女同一区二区软件| 精品熟女少妇av免费看| 婷婷色av中文字幕| 啦啦啦在线观看免费高清www| 麻豆成人av视频| 国产亚洲av嫩草精品影院| 日韩欧美精品v在线| 国产成人a区在线观看| 日本av手机在线免费观看| 日本黄大片高清| 如何舔出高潮| 高清av免费在线| 五月开心婷婷网| 亚洲人与动物交配视频| 国产精品久久久久久精品电影| 在线播放无遮挡| 最近手机中文字幕大全| 亚洲熟女精品中文字幕| 精品视频人人做人人爽| 日日摸夜夜添夜夜爱| av网站免费在线观看视频| 丝瓜视频免费看黄片| 亚洲精品国产成人久久av| 国产老妇女一区| 欧美丝袜亚洲另类| 美女国产视频在线观看| 国产伦在线观看视频一区| 亚洲欧美清纯卡通| 欧美日韩在线观看h| 黄色欧美视频在线观看| 亚洲国产精品成人综合色| 男女边摸边吃奶| 国产日韩欧美亚洲二区| 国产综合精华液| 成人毛片60女人毛片免费| 国产欧美日韩一区二区三区在线 | 午夜精品一区二区三区免费看| av在线播放精品| 日韩人妻高清精品专区| 久久精品国产鲁丝片午夜精品| 嫩草影院精品99| 日本免费在线观看一区| 欧美 日韩 精品 国产| 精品久久久久久久久av| 国产男女超爽视频在线观看| 热re99久久精品国产66热6| 自拍欧美九色日韩亚洲蝌蚪91 | 麻豆成人av视频| 好男人在线观看高清免费视频| 欧美bdsm另类| 午夜免费男女啪啪视频观看| 亚洲欧美清纯卡通| eeuss影院久久| 午夜免费男女啪啪视频观看| 禁无遮挡网站| 黄片wwwwww| 国产又色又爽无遮挡免| 国产亚洲精品久久久com| 黄片wwwwww| 午夜免费男女啪啪视频观看| 高清欧美精品videossex| 婷婷色综合大香蕉| 免费看不卡的av| 国内精品宾馆在线| 亚洲欧洲国产日韩| 色5月婷婷丁香| 国产人妻一区二区三区在| 久久久久久伊人网av| 不卡视频在线观看欧美| 51国产日韩欧美| 爱豆传媒免费全集在线观看| 日韩,欧美,国产一区二区三区| 久久精品夜色国产| 亚洲一级一片aⅴ在线观看| 九草在线视频观看| 麻豆乱淫一区二区| 国产成人精品婷婷| 在线a可以看的网站| 啦啦啦啦在线视频资源| 亚洲欧洲日产国产| 夜夜看夜夜爽夜夜摸| 亚洲人成网站在线播| 嫩草影院精品99| 午夜福利网站1000一区二区三区| 尤物成人国产欧美一区二区三区| 青春草亚洲视频在线观看| 日韩欧美 国产精品|