基于廣義Burgers方程的超聲速客機遠場聲爆高精度預測方法
發(fā)布時間:2021-07-03 06:17
聲爆高精度預測技術(shù)是新一代環(huán)保型超聲速客機設(shè)計的核心關(guān)鍵技術(shù)之一。基于廣義Burgers方程發(fā)展了可考慮"大氣風"效應(yīng)的遠場高精度預測方法,開發(fā)了聲爆預測程序"bBoom",并研究了近場聲爆信號提取位置和"大氣風"對遠場聲爆計算結(jié)果的影響。首先,給出了廣義Burgers方程和聲爆傳播射線的計算方法,重點討論了方程中計算分子弛豫效應(yīng)和熱黏吸收效應(yīng)的關(guān)鍵參數(shù)。其次,通過簡單軸對稱構(gòu)型標模算例、NASA的C25D構(gòu)型和洛馬的LM1021構(gòu)型等復雜超聲速客機算例對所發(fā)展的方法進行了驗證,表明本文發(fā)展的方法在預測遠場聲爆時具有較高可信度。最后,基于所發(fā)展的方法對比了由不同近場提取位置傳播到地面的波形,研究了飛機向不同方向飛行時,"大氣風"對地面聲爆強度和地面影響域的影響。結(jié)果表明:對于類C25D標模構(gòu)型,為了確保遠場聲爆預測結(jié)果具有較高精度,應(yīng)取機身下方約3倍機身長度位置處的近場壓強信號作為傳播方程的輸入;另外,"大氣風"會影響地面聲爆強度及地面影響域,在預測時有必要加以考慮。
【文章來源】:空氣動力學學報. 2019,37(04)北大核心CSCD
【文章頁數(shù)】:12 頁
【部分圖文】:
圖1波音公司擬研發(fā)的下一代超聲速客機[4]Fig.1Thenextgenerationofsupersonictransport[4]
?。對(10)式進行離散,可得:R(i+1)=R(i)+ΔR(i)N(i+1)=N(i)+ΔN(i){(11)其中,ΔR(i)=(c0N+W)|iΔtΔN(i)=[(I-N?NT)!(c0+W·N)]|iΔt烅烄烆(12)式中:Δt為聲爆在空間傳播過程中第i步到第i+1步的時間間隔。在計算聲線管面積時,通常的做法是計算由四條聲線圍成的“管”[41],如圖2所示。圖中Δφ為聲線在飛機周向上的角度增量。則聲線管面積的計算式為:A=12{(R4-R1)×(R3-R2)}·N1(13)1.3數(shù)值求解方法用數(shù)值方法求解方程(1)時,直接對分子弛豫過程進行離散比較困難,通常的求解策略是采用“算子分裂法”[22]。在很小的空間推進步Δσ內(nèi),分別單獨計算式(1)中各個效應(yīng)對聲學壓強的影響,即依次單獨求解式(14)-(18),并將前一方程的解作為后一方程的輸入,從而達到各個效應(yīng)解耦的目的。已經(jīng)證明,當Δσ足夠小時,采用算子分裂法的計算結(jié)果收斂于式(1)的解[42-43]。圖2聲線管面積求解示意圖Fig.2Sketchofraytubeareacalculation?P?σ=∑jCj(1+θj??τ)?2P?τ2(14)?P?σ=1Γ?2P?τ2(15)?P?σ=-12A?A?σ
圖3所開發(fā)的bBoom程序求解廣義Burgers方程的框架Fig.3Frameworkof“bBoom”codesolvingaugmentedBurgersequationusinganoperatorsplittingmethod分裂法將式(1)右邊各項效應(yīng)依次進行單獨求解,每一效應(yīng)求解過程中均采用有限差分方法;最后,沿傳播路徑重復上一過程,直至傳播到地面,即獲得地面聲爆波形。圖4為所開發(fā)的“bBoom”程序中,坐標系定義及聲爆傳播周向角?的定義。坐標原點O定義為聲爆由機體開始向外傳播時的飛機位置,x軸為沿飛機軸線指向飛行方向,y軸為垂直于飛機軸線指向飛行員的右側(cè),z軸根據(jù)右手坐標系定義,其垂直于xOy平面向下。聲爆傳播周向角定義為由z軸向y偏轉(zhuǎn)時為正,?=0°時代表聲爆向飛機正下方傳播(Under-track)。另外,定義x軸正向與正北向的夾角β為飛機飛行方位角。針對有“大氣風”的大氣剖面,該角度對地面聲爆計算具有重要作用。圖4“bBoom”程序所采用的坐標系統(tǒng)及傳播周向角定義Fig.4Definitionofcoordinatesystemandrollangleofsonic-boompropagationinthecode“bBoom”2算例驗證與討論2.1標模算例驗證以AIAA第二屆聲爆預測研討會(SBPW-2)[44-46]提供的軸對稱標模、低聲爆概念機C25D構(gòu)型和LM1021構(gòu)型為例,驗證所
【參考文獻】:
期刊論文
[1]網(wǎng)格對聲爆近場預測影響的數(shù)值研究[J]. 馬博平,王剛,雷知錦,葉正寅. 西北工業(yè)大學學報. 2018(05)
[2]基于增廣Burgers方程的音爆遠場計算及應(yīng)用[J]. 張繹典,黃江濤,高正紅. 航空學報. 2018(07)
[3]典型標模音爆的數(shù)值預測與分析[J]. 王剛,馬博平,雷知錦,任炯,葉正寅. 航空學報. 2018(01)
[4]典型低音爆構(gòu)型的近場音爆計算研究[J]. 徐悅,宋萬強. 航空科學技術(shù). 2016(07)
[5]低聲爆靜音錐設(shè)計方法研究[J]. 馮曉強,宋筆鋒,李占科. 航空學報. 2013(05)
[6]超音速客機音爆問題初步研究[J]. 馮曉強,李占科,宋筆鋒. 飛行力學. 2010(06)
本文編號:3262042
【文章來源】:空氣動力學學報. 2019,37(04)北大核心CSCD
【文章頁數(shù)】:12 頁
【部分圖文】:
圖1波音公司擬研發(fā)的下一代超聲速客機[4]Fig.1Thenextgenerationofsupersonictransport[4]
?。對(10)式進行離散,可得:R(i+1)=R(i)+ΔR(i)N(i+1)=N(i)+ΔN(i){(11)其中,ΔR(i)=(c0N+W)|iΔtΔN(i)=[(I-N?NT)!(c0+W·N)]|iΔt烅烄烆(12)式中:Δt為聲爆在空間傳播過程中第i步到第i+1步的時間間隔。在計算聲線管面積時,通常的做法是計算由四條聲線圍成的“管”[41],如圖2所示。圖中Δφ為聲線在飛機周向上的角度增量。則聲線管面積的計算式為:A=12{(R4-R1)×(R3-R2)}·N1(13)1.3數(shù)值求解方法用數(shù)值方法求解方程(1)時,直接對分子弛豫過程進行離散比較困難,通常的求解策略是采用“算子分裂法”[22]。在很小的空間推進步Δσ內(nèi),分別單獨計算式(1)中各個效應(yīng)對聲學壓強的影響,即依次單獨求解式(14)-(18),并將前一方程的解作為后一方程的輸入,從而達到各個效應(yīng)解耦的目的。已經(jīng)證明,當Δσ足夠小時,采用算子分裂法的計算結(jié)果收斂于式(1)的解[42-43]。圖2聲線管面積求解示意圖Fig.2Sketchofraytubeareacalculation?P?σ=∑jCj(1+θj??τ)?2P?τ2(14)?P?σ=1Γ?2P?τ2(15)?P?σ=-12A?A?σ
圖3所開發(fā)的bBoom程序求解廣義Burgers方程的框架Fig.3Frameworkof“bBoom”codesolvingaugmentedBurgersequationusinganoperatorsplittingmethod分裂法將式(1)右邊各項效應(yīng)依次進行單獨求解,每一效應(yīng)求解過程中均采用有限差分方法;最后,沿傳播路徑重復上一過程,直至傳播到地面,即獲得地面聲爆波形。圖4為所開發(fā)的“bBoom”程序中,坐標系定義及聲爆傳播周向角?的定義。坐標原點O定義為聲爆由機體開始向外傳播時的飛機位置,x軸為沿飛機軸線指向飛行方向,y軸為垂直于飛機軸線指向飛行員的右側(cè),z軸根據(jù)右手坐標系定義,其垂直于xOy平面向下。聲爆傳播周向角定義為由z軸向y偏轉(zhuǎn)時為正,?=0°時代表聲爆向飛機正下方傳播(Under-track)。另外,定義x軸正向與正北向的夾角β為飛機飛行方位角。針對有“大氣風”的大氣剖面,該角度對地面聲爆計算具有重要作用。圖4“bBoom”程序所采用的坐標系統(tǒng)及傳播周向角定義Fig.4Definitionofcoordinatesystemandrollangleofsonic-boompropagationinthecode“bBoom”2算例驗證與討論2.1標模算例驗證以AIAA第二屆聲爆預測研討會(SBPW-2)[44-46]提供的軸對稱標模、低聲爆概念機C25D構(gòu)型和LM1021構(gòu)型為例,驗證所
【參考文獻】:
期刊論文
[1]網(wǎng)格對聲爆近場預測影響的數(shù)值研究[J]. 馬博平,王剛,雷知錦,葉正寅. 西北工業(yè)大學學報. 2018(05)
[2]基于增廣Burgers方程的音爆遠場計算及應(yīng)用[J]. 張繹典,黃江濤,高正紅. 航空學報. 2018(07)
[3]典型標模音爆的數(shù)值預測與分析[J]. 王剛,馬博平,雷知錦,任炯,葉正寅. 航空學報. 2018(01)
[4]典型低音爆構(gòu)型的近場音爆計算研究[J]. 徐悅,宋萬強. 航空科學技術(shù). 2016(07)
[5]低聲爆靜音錐設(shè)計方法研究[J]. 馮曉強,宋筆鋒,李占科. 航空學報. 2013(05)
[6]超音速客機音爆問題初步研究[J]. 馮曉強,李占科,宋筆鋒. 飛行力學. 2010(06)
本文編號:3262042
本文鏈接:http://sikaile.net/kejilunwen/lxlw/3262042.html
最近更新
教材專著