CN106682392B - 复杂高超声速飞行器烧蚀效应快速计算方法 - Google Patents
复杂高超声速飞行器烧蚀效应快速计算方法 Download PDFInfo
- Publication number
- CN106682392B CN106682392B CN201611056358.1A CN201611056358A CN106682392B CN 106682392 B CN106682392 B CN 106682392B CN 201611056358 A CN201611056358 A CN 201611056358A CN 106682392 B CN106682392 B CN 106682392B
- Authority
- CN
- China
- Prior art keywords
- ablation
- heat
- carbon
- formula
- layer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000002679 ablation Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 66
- 230000008569 process Effects 0.000 title claims abstract description 47
- 238000004364 calculation method Methods 0.000 title claims abstract description 39
- 239000003575 carbonaceous material Substances 0.000 claims abstract description 11
- 238000004088 simulation Methods 0.000 claims abstract description 11
- 239000000463 material Substances 0.000 claims description 50
- 239000010410 layer Substances 0.000 claims description 44
- 238000006243 chemical reaction Methods 0.000 claims description 25
- 238000010438 heat treatment Methods 0.000 claims description 23
- 238000007254 oxidation reaction Methods 0.000 claims description 22
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 19
- 230000003647 oxidation Effects 0.000 claims description 18
- 229910052799 carbon Inorganic materials 0.000 claims description 17
- 238000004821 distillation Methods 0.000 claims description 12
- 238000012546 transfer Methods 0.000 claims description 10
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 9
- 229910052760 oxygen Inorganic materials 0.000 claims description 9
- 239000001301 oxygen Substances 0.000 claims description 9
- 239000002344 surface layer Substances 0.000 claims description 9
- 239000012530 fluid Substances 0.000 claims description 8
- 238000004134 energy conservation Methods 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 7
- CKUAXEQHGKSLHN-UHFFFAOYSA-N [C].[N] Chemical compound [C].[N] CKUAXEQHGKSLHN-UHFFFAOYSA-N 0.000 claims description 6
- 239000000470 constituent Substances 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- IJGRMHOSHXDMSA-UHFFFAOYSA-N nitrogen Substances N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 4
- 229910002091 carbon monoxide Inorganic materials 0.000 claims description 4
- 229910052757 nitrogen Inorganic materials 0.000 claims description 4
- 238000006213 oxygenation reaction Methods 0.000 claims description 4
- 238000005094 computer simulation Methods 0.000 claims description 3
- 239000004615 ingredient Substances 0.000 claims description 3
- 238000000859 sublimation Methods 0.000 claims description 3
- 230000008022 sublimation Effects 0.000 claims description 3
- 238000005092 sublimation method Methods 0.000 claims description 3
- 125000004430 oxygen atom Chemical group O* 0.000 claims description 2
- 230000005855 radiation Effects 0.000 claims 1
- 238000009826 distribution Methods 0.000 abstract description 32
- 230000004907 flux Effects 0.000 abstract description 22
- 238000013461 design Methods 0.000 abstract description 4
- 239000007770 graphite material Substances 0.000 description 20
- 238000004458 analytical method Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012407 engineering method Methods 0.000 description 3
- 230000001681 protective effect Effects 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 229910002804 graphite Inorganic materials 0.000 description 2
- 239000010439 graphite Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- XCWPUUGSGHNIDZ-UHFFFAOYSA-N Oxypertine Chemical compound C1=2C=C(OC)C(OC)=CC=2NC(C)=C1CCN(CC1)CCN1C1=CC=CC=C1 XCWPUUGSGHNIDZ-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- CREMABGTGYGIQB-UHFFFAOYSA-N carbon carbon Chemical compound C.C CREMABGTGYGIQB-UHFFFAOYSA-N 0.000 description 1
- 239000011203 carbon fibre reinforced carbon Substances 0.000 description 1
- 229910052681 coesite Inorganic materials 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 229910052906 cristobalite Inorganic materials 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010494 dissociation reaction Methods 0.000 description 1
- 208000018459 dissociative disease Diseases 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 239000012774 insulation material Substances 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000000197 pyrolysis Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 229910052682 stishovite Inorganic materials 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000006557 surface reaction Methods 0.000 description 1
- 239000005439 thermosphere Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 229910052905 tridymite Inorganic materials 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种复杂高超声速飞行器烧蚀效应快速计算方法,属飞行器气动设计领域。该方法采用了基于碳基材料的烧蚀模型和基于温度分区的烧蚀效应计算方法,结合普朗特理论的无粘外流解与边界层理论,发展了一种碳基材料烧蚀效应快速数值计算技术,给出了复杂外形高超声速飞行器热防护系统表面烧蚀的热流密度分布,烧蚀率以及热防护结构温度分布的时变特性。本发明同时增加了高超声速飞行器在高速飞行弹道过程中,热防护系统表面发生烧蚀情况的模拟能力。该方法计算效率高,鲁棒性好,结果准确,填补了国内在高超声速飞行器发生烧蚀情况下,高效、快速、准确的给出数值模拟结果方面的空白,为高超声速飞行器初期设计阶段提供有效的技术支持。
Description
技术领域:
本发明涉及一种复杂高超声速飞行器烧蚀效应快速计算方法,属飞行器气动设计领域。
背景技术:
随着高超声速飞行器的迅速发展,飞行速度越来越快,飞行条件越来越严苛,碳-碳复合材料逐渐取代了早期防护材料。热防护系统方案设计的关键问题之一是防热材料发生烧蚀的情况下,对气动加热的快速准确预测。
目前,对于高超声速飞行器的烧蚀防热问题,工程方法、数值方法和工程与数值相结合的方法的发展已相对成熟。实验是研究空气动力学问题最早的手段,各种研究方法所得的结果,一般以实验结果为准,其地位无可取代。目前,国内外仍是以实验方法为主要研究手段,因此,发展数值模拟计算预测变得十分重要。
本发明采用了基于碳基材料的烧蚀模型和基于温度分区的烧蚀效应计算方法,结合普朗特理论的无粘外流解与边界层理论,发展了一种碳基材料烧蚀效应快速数值计算技术,给出了复杂外形高超声速飞行器热防护系统表面烧蚀的热流密度分布,烧蚀率以及热防护结构温度分布的时变特性。
该方法增加了高超声速飞行器在高速飞行过程中,热防护系统表面发生烧蚀情况的模拟能力。该方法计算效率高,鲁棒性好,结果准确,填补了国内在高超声速飞行器发生烧蚀情况下,高效、快速、准确的给出数值模拟结果方面的空白,在高超声速飞行器初期设计阶段,提供有效的技术支持。
发明内容:
本发明主要采用数值计算手段,围绕高超声速气动加热和烧蚀效应,发展了一种复杂高超声速飞行器烧蚀效应快速计算技术。
本发明采用如下技术方案:一种复杂高超声速飞行器烧蚀效应快速计算方法,包括如下步骤:
步骤1、高超声速飞行器气动加热;
步骤2、材料的烧蚀模型;
步骤3、烧蚀效应快速计算;
步骤4、利用动态插值技术,对弹道状态进行准动态模拟。
进一步地,步骤1具体包括:
步骤1-1、无粘外流,采用欧拉方程作为控制方程求解无粘流动,取物面参数近似作为边界层的外缘参数,其中,所需要的边界层外缘参数有:网格信息、速度分量、马赫数、静压、密度、静温以及总能;
步骤1-2、高速边界层传热,计算表面热流密度需分为两个部分:驻点区和非驻点区,在驻点区,在任何情况下,驻点的流动都存在精确相似解,采用Fay-Riddell方法,驻点热流密度计算公式表示如下:
式中,Pr=0.71,Le=1.0;式中ρw为物面密度,μw为物面粘性系数,hw为物面焓值,
ρs为驻点处密度,μs为驻点处粘性系数,hD为平均空气电离焓。
步骤1-3、气动加热流固耦合计算,对于结构传热,根据热防护结构的毕奥数Bi大小,防热层可以分为热薄壁和热厚壁两种类型;
对于热薄壁,采用一维热传导方程,计算公式如下:
初始条件:Tw|t=0=T0
式中,ρδ为材料密度;cδ为材料比热容;δ为材料厚度;α为传热系数;ε为材料辐射系数;
对上式构建差分方程进行求解,并使用代替如下:
对于热厚壁,在计算中,由内向外热厚壁分为j层,简化后的热传导方程为:
(1)表层:
(2)最内层:
(3)中间层:
初始条件:Tj|t=0=T1|t=0=Tn|t=0=T0
其中下标m为材料类型,n为第n层,λm为热传导系数;
同样构建差分方程,并做近似处理的近似代替,的近似代替,如下:
(1)表层:
(2)最内层:
(3)中间层:
在tn时间步的求解中,先以上一时间步tn-1的壁面温度计算结构作为本时间步的壁面热边界条件,进行边界层的气动加热计算,得到tn时间步的热流密度接着,以该热流密度作为热边界条件,进行结构传热的计算,得到tn时间步的壁面温度然后又作为下一时间步tn+1流场计算的物面热边界条件。
进一步地,步骤2具体包括:
碳材料在较低温度下首先发生氧化,随着温度升高,氧化急剧增加,进入氧化扩散控制区,氧化率受边界层内氧气的输运快慢程度的控制,在更高温度下,碳-氮反应和碳的升华反应出现,引入质量烧蚀率和B无因此质量烧蚀率系数,下列计算公式中,C代表组元浓度,其下标表示组元的成分,同时下标中e与w分别表示边界层外缘条件和壁面条件,MO表示氧原子的分子量,Tw为表面温度,对于层流,对于湍流,qor为无烧蚀情况下的热流,hr为总焓,氧化受速率控制的过程Tw<1700K,忽略该温度范围内的碳氮反应以及碳的升华反应,按照下式计算:
其中,为氧化速率系数,
为当量扩散系数,αc(0)=ψqor/hr,
氧化受氧气运输快慢控制1700K≤Tw<3300K,忽略该温度范围被的碳氮反应以及碳的升华反应,按照下式计算:
在升华过程中3300K≤Tw,可按照下式计算:
根据能量守恒原则可以得到:
式中ΔHCO表示一氧化碳的生成热,h表示空气焓,下标r表示总焓,
在升华烧蚀阶段,空气温度较高,碳的氧化反应、氮化反应以及升华反应都已经出现,同样根据能量守恒原则得到:
其中,
ΔH表示反应热,表示平均升华热,下标V表示碳蒸汽。
进一步地,步骤3具体包括:在步骤2中引入qab表示考虑烧蚀效应的热流密度,在获得了qab之后,将其作为第二类边界条件带入结构传热计算当中
对于热薄壁,可以改写为:
对于热厚壁,可以改写为:
表层:
式中,对于氧化烧蚀:
对于升华烧蚀:
进一步地,步骤4具体包括:在不同时刻,沿整个弹道飞行仅需要四次无粘外流解的计算状态,即4个插值文件,假设在该时刻各个状态对应的马赫数和攻角分别为:(Ma1,α1),(Ma1,α2),(Ma2,α1),(Ma2,α2),对应的飞行器表面某点处某参数的值分别为A11,A12,A21,A22,则在所需要求解的马赫数攻角状态(Max,αx)下在参考高度上的该点处的该参数可以用以下一组公式进行计算:
本发明具有如下有益效果:
1.发明了碳基防热材料烧蚀模型及数值计算方法
基于温度分区的方法,建立了快速高效的碳基材料烧蚀工程计算方法。针对不同的材料,只需要给出相关物理参数,采用本发明的计算方法,就可对高速流动下其烧蚀特性进行研究。
2.发明了弹道状态动态插值技术
在不同时刻,沿整个弹道飞行仅需要四次无粘外流解的计算状态,即4个插值文件,采用该弹道状态动态插值技术,可以求解所需要求解的任意马赫数攻角状态下相同位置所对应的参数值。该技术增加了高超声速飞行器在高速飞行弹道过程中,热防护系统表面发生烧蚀情况的模拟能力。
3.发明了复杂高超声速飞行器烧蚀效应快速计算技术
发明了一种无粘欧拉解与工程方法相结合的数值计算方法用于快速计算与分析,复杂高超声速飞行器复杂高速流动(弹道状态)状态下,烧蚀效应对气动加热与防热结构材料传热特性的影响,给出了复杂外形高超声速飞行器热防护系统表面烧蚀的热流密度分布,烧蚀率以及热防护结构温度分布的时变特性。所发明的计算方法具有快速高效及工程设计可接受的计算精度,鲁棒性好等特点,填补了国内在高超声速飞行器发生烧蚀情况下,高效、快速和准确的给出数值模拟结果方面的空白,可为高超声速飞行器的热防护系统设计及热环境特性分析等提供技术参考与支持。
附图说明:
图1(a)为第0.05s时石墨材料表面热流密度分布。
图1(b)为第0.05s时509材料表面热流密度分布。
图2(a)为第100s时石墨材料表面热流密度分布。
图2(b)为第100s时509材料表面热流密度分布。
图3(a)为第200s时石墨材料表面热流密度分布。
图3(b)为第200s时509材料表面热流密度分布。
图4(a)为第300s时石墨材料表面热流密度分布。
图4(b)为第300s时509材料表面热流密度分布。
图5(a)为第400s时石墨材料表面热流密度分布。
图5(b)为第400s时509材料表面热流密度分布。
图6(a)为第500s时石墨材料表面热流密度分布。
图6(b)为第500s时509材料表面热流密度分布。
图7(a)为第0.05s时石墨材料表面温度分布。
图7(b)为第0.05s时509材料表面温度分布。
图8(a)为第100s时石墨材料表面温度分布。
图8(b)为第100s时509材料表面温度分布。
图9(a)为第200s时石墨材料表面温度分布。
图9(b)为第200s时509材料表面温度分布。
图10(a)为第300s时石墨材料表面温度分布。
图10(b)为第300s时509材料表面温度分布。
图11(a)为第400s时石墨材料表面温度分布。
图11(b)为第400s时509材料表面温度分布。
图12(a)为第500s时石墨材料表面温度分布。
图12(b)为第500s时509材料表面温度分布。
图13为复杂高超声速飞行器烧蚀效应快速计算技术方法计算流程图。
图14为第500s时对称面内外层温度分布。
图15为RAMC-II对称面热流密度分布。
图16为RAMC-II对称面温度分布。
图17为RAMC-II驻点内外层温度随时间的变化。
图18为烧蚀气动加热和结构传热耦合计算步骤。
具体实施方式:
本发明复杂高超声速飞行器烧蚀效应快速计算技术,其具体包括如下步骤:
步骤1、高超声速飞行器气动加热问题。本发明采用工程方法与数值模拟计算相结合的算法:无粘外流解、高速边界层传热以及结构传热三者耦合的方法。其大致分为两个部分:其一是分析飞行器外部高超声速流场的传热,得到飞行器表面热流密度的分布;其二是表面热流进入飞行器被捕之后,分析在热防护结构内的传热问题。
步骤2、材料的烧蚀模型。本发明针对热防护材料的烧蚀问题,采用了基于碳基材料的烧蚀模型和基于温度分区的烧蚀效应计算方法,该方法能正确地表征在烧蚀过程的传质和传热机理,给出烧蚀效应较好的模拟结果。
步骤3、烧蚀效应快速计算方法。将烧蚀效应的计算方法耦合到气动加热计算方法之中,提出一种能够考虑碳基材料烧蚀效应的气动加热计算技术,从而可以得到高超声速飞行器热防护系统表面烧蚀的热流密度分布,烧蚀率以及防热结构的温度分布。
步骤4、利用动态插值技术,对弹道状态进行准动态模拟。
进一步地,步骤1具体包括:
步骤1-1、无粘外流。由于高超声速边界层厚度很薄,在求解边界层外的无粘流动时忽略边界层的厚度,采用Euler方程作为控制方程求解无粘流动,取物面参数近似作为边界层的外缘参数。其中,所需要的边界层外缘参数有:网格信息(物面网格点坐标和网格单元信息)、速度分量(物面网格点上三个方向的速度分量)、马赫数、静压、密度、静温以及总能。
步骤1-2、高速边界层传热。计算表面热流密度需分为两个部分:驻点区和非驻点区。
在驻点区,在任何情况下,驻点的流动都存在精确相似解,本发明采用Fay-Riddell方法,驻点热流密度计算公式表示如下:
式中,Pr=0.71,Le=1.0。
在非驻点区,热力学参数以及输运参数沿物面方向不变化的假设不成立,从而得不到精确的相似解。但是在高冷壁条件下,仍可假设,沿法向各热力学参数的变化远大于沿物面方程的变化,作局部相似。本发明主要是应用平板传热模型计算传热系数和热流密度,对于较复杂外形的计算进行适当的简化处理。
步骤1-3、气动加热流固耦合计算。飞行器表面以及热防护层的温度处在不断变化之中,而热流密度的计算需要物面温度等参数作为条件,两个重要物理量的求解各自需要彼此作为求解条件,这样就使得表面温度与热流密度直接存在耦合关系。
对于结构传热,根据热防护结构的毕奥数Bi大小,防热层可以分为热薄壁和热厚壁两种类型。
对于热薄壁,采用一维热传导方程,计算公式如下:
初始条件:Tw|t=0=T0
式中,ρδ为材料密度;cδ为材料比热容;δ为材料厚度;α为传热系数;ε为材料辐射系数。对上式构建差分方程进行求解,并使用代替如下:
对于热厚壁,在计算中,由内向外热厚壁分为j层,简化后的热传导方程为:
(1)表层:
(2)最内层:
(3)中间层:
初始条件:Tj|t=0=T1|t=0=Tn|t=0=T0
其中下标m为材料类型,n为第n层,λm为热传导系数。
同样构建差分方程,并做近似处理(为的近似代替,为的近似代替),如下:
(1)表层:
(2)最内层:
(3)中间层:
而气动加热和结构传热要进行耦合计算——气动加热与结构传热的交叉迭代,即在tn时间步的求解中,先以上一时间步tn-1的壁面温度计算结构作为本时间步的壁面热边界条件,进行边界层的气动加热计算,得到tn时间步的热流密度接着,以该热流密度作为热边界条件,对结构传热进行计算,得到tn时间步的壁面温度然后又作为下一时间步tn+1流场计算的物面热边界条件。如此循环往复,实现结构传热和气动热环境计算的交替进行和相互影响,达到综合考虑气动加热与结构传热的目的。
进一步地,步骤2具体包括:
碳材料在较低温度下首先发生的化学反应是氧化,起初氧化是速率控制的过程,表面反应动力学条件决定其氧化率。随着温度逐渐升高,氧化反应急剧增加,氧气逐渐供应不足,进入氧化扩散控制区,这时,氧化率受边界层内氧气的输运快慢程度的控制。随着温度进一步升高,碳-氮反应和碳的升华反应逐渐出现。引入质量烧蚀率和B无因此质量烧蚀率系数。下列计算公式中,C代表组元浓度,其下标表示组元的成分,同时下标中e与w分别表示边界层外缘条件和壁面条件,MO表示氧原子的分子量,Tw为表面温度。对于层流,对于湍流,qor为无烧蚀情况下的热流,hr为总焓。
氧化受速率控制的过程Tw<1700K,忽略该温度范围内的碳氮反应以及碳的升华反应。可按照下式计算:
其中,为氧化速率系数。
为当量扩散系数。αc(0)=ψqor/hr。
氧化受氧气运输快慢控制1700K≤Tw<3300K,忽略该温度范围被的碳氮反应以及碳的升华反应。可按照下式计算:
在升华过程中3300K≤Tw,可按照下式计算:
作为简化假设,认为在碳基材料的氧化烧蚀阶段,氧化产物只有,而此时空气温度较低,没有出现氧气和氮气等的解离反应。因此,根据能量守恒原则可以得到:
式中ΔHCO表示一氧化碳的生成热,h表示空气焓,下标r表示总焓。
在升华烧蚀阶段,空气温度较高,碳的氧化反应、氮化反应以及升华反应都已经出现,同样根据能量守恒原则得到:
其中,
ΔH表示反应热,表示平均升华热,下标V表示碳蒸汽。
进一步地,步骤3具体包括:
在步骤2中提出了基于温度分区基于温度分区的烧蚀效应计算方法,给出了基于能量守恒原则表面热流密度计算的表达式。引入qab表示考虑烧蚀效应的热流密度,在获得了qab之后,可以将其作为第二类边界条件带入结构传热计算当中。
对于热薄壁,可以改写为:
对于热厚壁,可以改写为:
表层:
式中,对于氧化烧蚀:
对于升华烧蚀:
在给出初始时刻温度和材料烧蚀物理属性后,按照时间步长推进,就可以得到考虑烧蚀效应的气动加热下的结构温度分布。
进一步地,步骤4具体包括:
弹道状态是指飞行器在沿着弹道飞行的过程中,飞行器随着时间的变化改变其高度、迎角、马赫数三个参数(可将侧滑角忽略)。对整个弹道状态,选取不同的时间,通过对烧蚀效应的计算,得到各物理量在不同材料的飞行器表面的变化。下面,简要描述对无粘外流解弹道状态插值技术:
在不同时刻,均对高度效应作简化处理,沿整个弹道飞行仅需要四次无粘外流解的计算状态(根据弹道的时间、马赫数和迎角参数等的变化),即4个插值文件。假设在该时刻各个状态对应的马赫数和攻角分别为:(Ma1,α1),(Ma1,α2),(Ma2,α1),(Ma2,α2),对应的飞行器表面某点处某参数的值分别为A11,A12,A21,A22,则在所需要求解的马赫数攻角状态(Max,αx)下在参考高度上的该点处的该参数可以用以下一组公式进行计算:
下面结合算例对本发明进行说明。
1.技术参数
采用的算例的工程背景是巡航状态下的高超声速钝头体气动加热问题。计算模型是RAMC-II,是一个总长1.295m,头部半径为0.1525m,半锥角9°的顿锥。设定巡航状态:马赫数15、巡航高度70km、巡航时间500;热防护材料:热薄壁(石墨、厚度2mm)+隔热层(SiO2、厚度20mm),表面发射率0.8。该算例所得带烧蚀效应的气动热结果,将对比不同烧蚀材料对于烧蚀效应计算的影响。热防护设置如下表所示:
2.数值模拟结果
整个算例的数值模拟在普通个人PC上运行,算例模拟计算时间为CPU耗时29分钟,国内暂无出现同等高效的高超声速飞行器烧蚀效应计算技术,属国内首次,下面为计算结果分析:
从图1(a)~图6(b)中可以发现,不论石墨材料还是509材料,由于碳基材料以及隔热层的导热性较差,表面热流密度在很短的时间内就达到了平衡;对比相同时间的表面热流密度,石墨材料的表面热流密度要高于509材料的。从图7(a)~图12(b)可以发现,在烧蚀过程中,石墨材料的烧蚀过程是一个放热过程,固态燃烧产生一氧化碳,释放热量,导致表面热流密度高于没有烧蚀的表面热流密度。相反,509材料的烧蚀过程是一个吸热过程,材料热解吸收的热量不仅超过了燃烧产生的热量,还包括了没有烧蚀的气动加热量,导致了表面热流密度的减小。
从图14中可以看出,在对称面上,表面温度明显高于内层温度,头部区域的最高温差达到1700K,在头部与身部的过渡区,温度梯度最大。
从图15和图16中可以看出,热流密度分布的不同导致了表面温度分布的不同。在初始时刻,两种烧蚀材料的表面温度都很低,随着时间的增加,加热效果明显,高温区逐步向下游地区延伸;在任意相同的时刻、相同位置,由于石墨烧蚀是放热过程,而509烧蚀是吸热过程,石墨材料的表面温度要高于509材料。
从图17可以看出,驻点表层温度在巡航的初始阶段就达到了平衡,石墨材料表面温度升高的趋势较509材料更剧烈一些,石墨材料的平衡温度在2500K左右,而509材料的平衡温度较低,在2000K左右。两种材料的内层温度升高趋势较表层温度要缓和得多,且都还没有到达平衡值。因为509材料的表层温度要低于石墨材料,所以509材料的内层温度也低于石墨材料。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下还可以作出若干改进,这些改进也应视为本发明的保护范围。
Claims (2)
1.一种复杂高超声速飞行器烧蚀效应快速计算方法,其特征在于:包括如下步骤
步骤1、高超声速飞行器气动加热;
步骤2、材料的烧蚀模型;
步骤3、烧蚀效应快速计算;
步骤4、利用动态插值技术,对弹道状态进行准动态模拟;步骤1具体包括:
步骤1-1、无粘外流,采用欧拉方程作为控制方程求解无粘流动,取物面参数近似作为边界层的外缘参数,其中,所需要的边界层外缘参数有:网格信息、速度分量、马赫数、静压、密度、静温以及总能;
步骤1-2、高速边界层传热,计算表面热流密度需分为两个部分:驻点区和非驻点区,在驻点区,在任何情况下,驻点的流动都存在精确相似解,采用Fay-Riddell方法,驻点热流密度计算公式表示如下:
式中,Pr=0.71,Le=1.0;式中ρw为物面密度,μw为物面粘性系数,hw为物面焓值,ρs为驻点处密度,μs为驻点处粘性系数,hD为平均空气电离焓;
步骤1-3、气动加热流固耦合计算,对于结构传热,根据热防护结构的毕奥数Bi大小,
防热层可以分为热薄壁和热厚壁两种类型;
对于热薄壁,采用一维热传导方程,计算公式如下:
初始条件:Tw|t=0=T0
式中,ρδ为材料密度;cδ为材料比热容;δ为材料厚度;α为传热系数;ε为材料辐射系数;
对上式构建差分方程进行求解,并使用代替如下:
对于热厚壁,在计算中,由内向外热厚壁分为j层,简化后的热传导方程为:
(1)表层:
(2)最内层:
(3)中间层:
初始条件:Tj|t=0=T1|t=0=Tn|t=0=T0
其中下标m为材料类型,n为第n层,λm为热传导系数;
同样构建差分方程,并做近似处理的近似代替,为的近似代替,如下:
(1)表层:
(2)最内层:
(3)中间层:
在tn时间步的求解中,先以上一时间步tn-1的壁面温度计算结构作为本时间步的壁面热边界条件,进行边界层的气动加热计算,得到tn时间步的热流密度接着,以该热流密度作为热边界条件,进行结构传热的计算,得到tn时间步的壁面温度然后又作为下一时间步tn+1流场计算的物面热边界条件;
步骤2具体包括:
碳材料在较低温度下首先发生氧化,随着温度升高,氧化急剧增加,进入氧化扩散控制区,氧化率受边界层内氧气的输运快慢程度的控制,在更高温度下,碳-氮反应和碳的升华反应出现,引入质量烧蚀率和B无因此质量烧蚀率系数,下列计算公式中,C代表组元浓度,其下标表示组元的成分,同时下标中e与w分别表示边界层外缘条件和壁面条件,MO表示氧原子的分子量,Tw为表面温度,对于层流,对于湍流,qor为无烧蚀情况下的热流,hr为总焓,氧化受速率控制的过程Tw<1700K,忽略该温度范围内的碳氮反应以及碳的升华反应,按照下式计算:
其中,为氧化速率系数,为当量扩散系数,αc(0)=ψqor/hr,
氧化受氧气运输快慢控制1700K≤Tw<3300K,忽略该温度范围被的碳氮反应以及碳的升华反应,按照下式计算:
在升华过程中3300K≤Tw,可按照下式计算:
根据能量守恒原则可以得到:
式中ΔHCO表示一氧化碳的生成热,h表示空气焓,下标r表示总焓,
在升华烧蚀阶段,空气温度较高,碳的氧化反应、氮化反应以及升华反应都已经出现,同样根据能量守恒原则得到:
其中,
ΔH表示反应热,表示平均升华热,下标V表示碳蒸汽;
步骤3具体包括:在步骤2中重新引入qab表示考虑烧蚀效应的热流密度,在获得了qab之后,将其作为第二类边界条件带入结构传热计算当中
对于热薄壁,可以改写为:
对于热厚壁,可以改写为:
表层:
式中,对于氧化烧蚀:
对于升华烧蚀:
2.如权利要求1所述的复杂高超声速飞行器烧蚀效应快速计算方法,其特征在于:在不同时刻,沿整个弹道飞行仅需要四次无粘外流解的计算状态,即4个插值文件,假设在该时刻各个状态对应的马赫数和攻角分别为:(Ma1,α1),(Ma1,α2),(Ma2,α1),(Ma2,α2),对应的飞行器表面某点处某参数的值分别为A11,A12,A21,A22,则在所需要求解的马赫数攻角状态(Max,αx)下在参考高度上的该点处的该参数可以用以下一组公式进行计算:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611056358.1A CN106682392B (zh) | 2016-11-24 | 2016-11-24 | 复杂高超声速飞行器烧蚀效应快速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611056358.1A CN106682392B (zh) | 2016-11-24 | 2016-11-24 | 复杂高超声速飞行器烧蚀效应快速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106682392A CN106682392A (zh) | 2017-05-17 |
CN106682392B true CN106682392B (zh) | 2019-07-19 |
Family
ID=58865900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611056358.1A Active CN106682392B (zh) | 2016-11-24 | 2016-11-24 | 复杂高超声速飞行器烧蚀效应快速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106682392B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107832494B (zh) * | 2017-10-13 | 2021-02-19 | 南京航空航天大学 | 高超声速飞行器前缘流-热-固一体化计算方法 |
CN108333035B (zh) * | 2017-10-20 | 2019-02-15 | 北京空天技术研究所 | 一种用于热防护结构的低温特性分析方法和系统 |
CN107808065B (zh) * | 2017-11-23 | 2019-12-31 | 南京航空航天大学 | 三维复杂外形高速飞行器流-固-热快速计算方法 |
CN107944137B (zh) * | 2017-11-23 | 2020-05-26 | 南京航空航天大学 | 高超声速飞行器弹道状态多场耦合的热气动弹性计算技术 |
CN108268702A (zh) * | 2017-12-27 | 2018-07-10 | 中国航天空气动力技术研究院 | 超声速湍流导管内试验件表面冷壁热流快速计算方法 |
CN110309552B (zh) * | 2019-06-10 | 2023-04-14 | 中国航天空气动力技术研究院 | 一种考虑质量引射效应的飞行器湍流预测方法及系统 |
CN110567413B (zh) * | 2019-08-16 | 2022-04-08 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种获取复合材料氧化膜层厚度的方法、装置及电子设备 |
CN111458366B (zh) * | 2020-04-17 | 2024-06-11 | 北京空天技术研究所 | 一种烧蚀热防护系统结构气动热/传热耦合分析方法 |
CN115544675B (zh) * | 2022-12-01 | 2023-04-07 | 北京航空航天大学 | 高超声速飞行器防热材料表面催化特性多尺度预测方法 |
CN115587551B (zh) * | 2022-12-12 | 2023-02-17 | 北京航空航天大学 | 一种高超声速飞行器防热结构烧蚀行为多尺度预测方法 |
CN117972279B (zh) * | 2024-03-28 | 2024-05-31 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种基于边界层理论计算球头驻点热流的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
CN106508022B (zh) * | 2010-12-31 | 2014-09-10 | 上海机电工程研究所 | 一种烧蚀式防热结构三维温度场分析计算方法 |
-
2016
- 2016-11-24 CN CN201611056358.1A patent/CN106682392B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106508022B (zh) * | 2010-12-31 | 2014-09-10 | 上海机电工程研究所 | 一种烧蚀式防热结构三维温度场分析计算方法 |
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
Non-Patent Citations (4)
Title |
---|
Approximate method for calculating heating rates on three-dimensional vehicles;Hamilton H H等;《Journal of Spacecraft and Rockets》;19940331;第31卷(第3期);345-354 |
带烧蚀效应的高超声速气动加热计算技术;施昆等;《2015年第二届中国航空科学技术大会论文集》;20150915;35-37 |
真实气体效应影响下的高超声速气动加热计算方法;孙伟;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20150215;C039-16 |
高超声速飞行器气动加热计算技术;王杰;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20130215;13-42 |
Also Published As
Publication number | Publication date |
---|---|
CN106682392A (zh) | 2017-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106682392B (zh) | 复杂高超声速飞行器烧蚀效应快速计算方法 | |
CN107808065B (zh) | 三维复杂外形高速飞行器流-固-热快速计算方法 | |
WO2020015375A1 (zh) | 一种气动热环境试验模拟条件的参数相似方法 | |
Hedlund et al. | Local swirl chamber heat transfer and flow structure at different Reynolds numbers | |
Coulthard et al. | Effect of jet pulsing on film cooling—part I: effectiveness and flow-field temperature results | |
Mayle et al. | Effect of streamline curvature on film cooling | |
CN109145388A (zh) | 航空发动机部件的热分析方法 | |
Rutledge et al. | Computational Fluid Dynamics Evaluations of Unconventional Film Cooling Scaling Parameters on a Simulated Turbine Blade Leading Edge | |
CN111458366B (zh) | 一种烧蚀热防护系统结构气动热/传热耦合分析方法 | |
CN110309552B (zh) | 一种考虑质量引射效应的飞行器湍流预测方法及系统 | |
CN106845072A (zh) | 多组分防热材料多反应机制控制下的烧蚀速率确定方法 | |
Bruce-Black et al. | Practical slot configurations for turbine film cooling applications | |
CN106960089A (zh) | 含内部复杂边界结构体的温度场和热流同时重构方法 | |
CN106897537A (zh) | 含三维或曲面外形结构体的温度场与热流同时重构方法 | |
Meng et al. | The performance evaluation for thermal protection of turbine vane with film cooling and thermal barrier coating | |
White | Numerical investigation of condensing steam flow in boundary layers | |
Yao et al. | Effect of density ratio on film-cooling effectiveness distribution and its uniformity for several hole geometries on a flat plate | |
Kovalnogov et al. | Modeling, research and development the technology of cooling of turbine engine blades | |
Womack et al. | Combined effects of wakes and jet pulsing on film cooling | |
Chen et al. | Roughness implementation and convective heat transfer coefficient computation toward ice accretion simulation | |
Reller | Heat transfer to blunt nose shapes with laminar boundary layers at high supersonic speeds | |
Ingram et al. | a superposition technique for multiple-row film-cooling for calculation of 2-d effectiveness and heat transfer coefficients | |
Ma et al. | Arrangement of Sensors for Measuring Temperature in the Test of Autoclave | |
Ames et al. | Effects of Catalytic and Dry Low NO x Combustor Turbulence on Endwall Heat Transfer Distributions | |
Sutton | The temperature history in a thick skin subjected to laminar heating during entry into the atmosphere |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |