CN107808065A - 三维复杂外形高速飞行器流‑固‑热快速计算方法 - Google Patents
三维复杂外形高速飞行器流‑固‑热快速计算方法 Download PDFInfo
- Publication number
- CN107808065A CN107808065A CN201711177571.2A CN201711177571A CN107808065A CN 107808065 A CN107808065 A CN 107808065A CN 201711177571 A CN201711177571 A CN 201711177571A CN 107808065 A CN107808065 A CN 107808065A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- rsqb
- lsqb
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种三维复杂外形高速飞行器流‑固‑热快速计算方法,属飞行器气动计算领域。本发明依据边界层理论,采用CFD数值计算与工程算法相结合的方法,静气动热的计算简化为绕飞行器的无粘外流数值解和边界层内热流求解两份部分,同时耦合防热系统结构传热计算模型。本发明结合了全长数值模拟与工程算法各自的优势,可以快速高效对三维复杂外形高速飞行器的多种飞行状态下气动加热与结构耦合传热特性的计算方法,给出飞行器全机表面热流密度与防热结构温度场时变特性,弥补了直接数值模拟方法计算代价高、效率低、周期长等缺陷,同时拓展了工程算法的应用范围。
Description
技术领域
本发明属飞行器气动计算技术领域,具体指代一种三维复杂外形高速飞行器流-固- 热快速计算方法。
背景技术
高超声速飞行器在大气层中持续飞行时,飞行器表面将承受巨大的气动热(王江峰, 伍贻兆.高超声速复杂气动问题数值方法研究进展[J].航空学报.2015,33(1):159-175)。气 动热会导致飞行器结构温度升高,改变结构内温度场分布,从而改变结构材料的物理属 性,而结构材料特性的改变会引起结构刚度和结构模态的变化,对飞行器的飞行安全造 成极大隐患。因此热防护问题对于高超声速飞行器设计显得尤为重要(程克明,吕英伟. 飞行器持续气动加热的耦合性分析[J].南京航空航天大学学报,2000,32(2):150-155),越 来越受到世界各国的高度重视。
气动加热问题(卞荫贵,徐立功.气动热力学[M].中国科技大学出版社.1997:15-20)的 主要研究方法包括:数值模拟方法(SinhaK,K.Reddy D S.Effect of ChemicalReaction Rates on Aeroheating Predictions of Reentry Flows[J].Journal ofThermophysics and Heat Transfer,2011,25(1):21-33)、工程近似计算方法(Hamilton HH,WeilmuensterKJ, DeJarnette F R.Ap-proximate Method for Computing ConvectiveHeating on Hypersonic Vehicles Using Unstructured Grids[J].Journal ofspacecraft and rockets.2014,1288:1305)、 地面风洞实验(田旭昂,王成鹏,程克明.Ma5斜激波串动态特性实验研究[J].推进技 术.2014,35(8):1030-1039)及自由飞实验等,其中后两种方法因代价高昂等因素不适于工 程设计的初期与选型、改型阶段。
数值模拟方法主要是对Euler方程、N-S(Navier-Stokes)方程及相关简化形式的控制 方程进行求解,具有计算精度高、可处理复杂流动和全机外形等优点,但在气动热与结 构传热耦合问题的求解方面对计算资源需求巨大且非常耗时。国内外在求解高超声速气 动热相关问题的数值计算方法方面做了大量研究,主要工作集中在计算效率与精度等方 面。
气动加热工程计算方法(Zhao J S,Gu L X,Ma H Z.A rapid approach toconvective aeroheating prediction of hypersonic vehicles[J].Science ChinaTechnological Sciences,2013: 1-15)对简单外形的气动热求解具有高效、准确的特点,因此在实际工程应用中率先得 到发展。但在复杂气动外形的气动热分析方面适应性较差,并且需要基于大量实验数据 对计算方法与结果进行人工修正。在这方面比较著名的是NASA兰利研究中心研发的一 套气动加热预测程序(MINIVER)(李会萍.高超声速飞行器气动加热特性及其计算方法 研究[D].上海:上海交通大学,2010),该程序中驻点区域使用经典的Fay-Riddle公式,采 用参考焓方法来计算高速流动压缩性效应,另外可对过渡流区气动热进行计算。 Hamilton(Hamilton H H,Weilmuenster KJ,DeJarnette FR.Approximate method for computing laminar and turbulent convective heatingon hypersonic vehicles using unstructured grids[C].The 41st Annual AIAAThermophysics Conference,AIAA Paper.2009, 4310:2009.)针对三维钝头体发展了一种适用于空气平衡气体的高超声速流动热流密度 计算方法,该方法可计算不同高速流动状态(层流、转捩、湍流)下的热流密度。Zoby(Zoby E V,Simmonds A L.Engineeringflowfield method with angle-of-attack applications[J]. Journal of spacecraftand rockets,1985,22(4):398-404)等人开发了LATCH方法,该方法 基于参考焓和修正雷诺比拟计算热流密度,可用于有化学反应参与的钝头体气动热计 算。李建林(李建林,唐乾刚,霍霖,程兴华.复杂外形高超声速飞行器气动热快速工程估算 [J].国防科技大学学报.2012,34(6):89-93)等人对气动热工程计算方法进行拓展应用,对乘 波体构型高速飞行器进行气动热快速估算,得到较好结果。总体而言,工程计算方法存 在需要大量人工干涉及以庞大实验数据为基础等缺陷,在处理复杂外形飞行器方面的通 用性仍需完善。
发明内容
针对于上述现有技术的不足,本发明的目的在于提供一种三维复杂外形高速飞行器 流-固-热快速计算方法,以解决直接数值模拟方法计算代价高、效率低、周期长等缺陷,同时拓展了工程算法的应用范围。
为达到上述目的,本发明采用的技术方案如下:
一种三维复杂外形高速飞行器流-固-热快速计算方法,其特征在于,包括步骤如下:
步骤1、无粘外流场数值解;
步骤2、工程表面热流率计算;
步骤3、工程结构传热耦合计算;
步骤4、弹道状态动态插值方法;
步骤5、高温化学非平衡效应。
步骤1具体包括:
控制方程采用三维欧拉方程对无粘外流场求解,采用块结构网格机翼基于分布式并 行计算技术,取无粘流场解结果中的物面参数作为边界层外缘参数提供给工程方法中的 边界层简化算法。输出的物面参数包括:坐标(m)、物面速度分量(m/s)、马赫数、压 强(Pa)、密度(kg/m3)、静温(K)、总能(J/m3)。而边界层的工程算法可以得到用于 防热系统中结构传热计算所需的飞行器表面热流及传热系数等参数。
步骤2具体包括:
本发明针对的是全机外形,根据热流密度工程计算方法将飞行器表面热流的计算划 分为驻点区和非驻点区两个区域。驻点区热流计算采用目前广泛使用的Fay-Riddell公式 (中国人民解放军总装备部军事训练教材编辑工作委员会.高超声速气动热和热防护[M].北京:国防工业出版社,2003:46-136.):
式中:ρw、μw、hw分别表示物面密度、物面粘性系数和物面焓,ρs、μs分别为驻点密 度及驻点粘性系数,hD为平均空气电离焓,计算中取普朗特数Pr=0.71、路易斯数 Le=1.0。
非驻点区的热流密度计算公式(卞荫贵,钟家康.高温边界层传热[M].科学出版社.1986.),本发明采用的是工程中常用的平板传热模型。
步骤3具体包括:
热防护结构内表面采用绝热壁边界条件,依据防热结构材料的毕奥数Bi的大小,采 不同的传热计算模型(航天工业部部标准,“战术导弹气动加热工程计算方 法”QJ1734-89,1989.):
式中α为传热系数,δ为材料结构厚度,λδ为当前材料热传导系数。
当Bi<0.1时,采用如下式所示的热薄壁传热模型:
初始条件为:
Tw|t=0=T0 (4)
式中ρδ、cδ、ε分别为表示分别表示防热材料的密度、比热容和表面辐射系数。采用差 分方法对式(3)进行加热时间t推进求解,即可得到热薄壁条件下热防护层温度随时间的 变化。
当Bi>0.1时,采用热厚壁传热模型。将热厚壁由内而外分为j层进行逐层推进求解:
表层:
最内层:
中间层:
初始条件为:
Tj|t=0=T1|t=0=Tn|t=0=T0 (8)
式中下标m为为材料类型,n表示层数,λm为材料m的热传导系数。利用差分离散可计 算求得热厚壁条件下热防护结构各层材料的温度随时间变化的结果。
气动加热与结构传热过程的耦合特性(吴洁,阎超.气动热与热响应的耦合研究[J]. 导弹与航天运载技术,2009(004):35-39.)以物面温度Tw作为物面边界输入条件,然后通 过简化边界层工程方法计算得到表面热流密度qn;而在计算结构传热时,又以热流密度 qn作为热边界条件,计算得到物面温度。在时间步tn的求解中,以上一时间步tn-1的物 面温度作为本时间步的物面温度边界条件,采用工程算法进行边界层的气动加热计 算,得出tn时间步的热流密度然后以作为热边界输入条件提供给结构传热计算 模块,计算得到tn时间步的物面温度如此循环以实现结构传热和气动热环境的耦合 计算。
步骤4具体包括:
计算非定常弹道状态的气动加热与结构传热耦合计算需要足够多的无粘外流的数值 解作为其输入参数。
本发明利用无粘外流解动态插值方法提高耦合算法计算效率:通过已经预先完成的 有限个无粘外流解的流场结果,插值得到当前弹道时间点上飞行条件下的流场解,以供气动加热与结构传热计算所需,具体步骤如下:
1)根据飞行器物面上的当地流场参数与来流参数的比值在不同飞行高度下几乎保持不 变(Bova S W,Howard M A.Coupling strategies for high-speed aeroheatingproblems[R]. Sandia National Laboratories,2011.)这一规律,对同一飞行器,已经某高度H0下的 边界层外缘参数P0以及该高度下的流动参数Pinf0,同时已知任意高度下的大气参数Pinfx,可以通过Px=Pinfx×P0/Pinf0的简单换算获得该飞行器在任意高度上的边界层 外缘参数。
2)由1),对于整个弹道状态,无粘外流解只需计算设定参考飞行高度下的两个马赫数 和两个迎角,即4个飞行状态,其他飞行状态则可根据插值方法得到。这四个计算 状态分别为设定飞行高度下的两个马赫数和两个迎角的组合,即: (Ma1,α1),(Ma1,α2),(Ma2,α1),(Ma2,α2)。与这四个状态所对应的飞行器表面任一点 处的流动参数的值分别记作P11、P12、P21、P22,则某个弹道时间点上飞行条件 (Max,αx)状态下的飞行器表面同一位置处的流动参数Pxx可通过如下插值算法得到:
步骤5具体包括:
在实际高超声速飞行条件下,空气的高温非平衡效应对气动热影响显著(欧阳水吾, 谢中强.高温非平衡空气扰流[M].国防工业出版社,2003)。高温空气非平衡边界层驻点 传热与平衡边界层驻点传热表面热通量的关系式表示为:
式中,下标O、N分别表示氧原子和氮原子,φ为壁面催化因子;h0为生成焓;CO.s和CN.s分别表示氧原子和氮原子的质量浓度;Le为路易斯数,取值范围1~2。
对于组元参数特性,使用组元参数简化计算模型,在流场温度9000K以下时,只考虑O2、O、O+、N2、N、N+、NO、NO+、e-等9种组元及以下6个化学反应方程式:
其中Ki为摩尔密度平衡常数。本发明暂不考虑NO、NO+、e-、N+、O+。
氧原子和氮原子的壁面催化因子有如下关系成立:
其中,[O2]0,[N2]0分别表示为空气中氧分子和氮分子的摩尔密度。联立式(11)、(12) 与反应方程①、②求解,得到氧原子和氮原子的摩尔密度:
取一级近似φO=φN=0。联立反应方程①②③,可得到[O2]、[N2]和[NO]的摩尔密度计算 式为:
联立反应方程④⑤⑥和电荷守恒方程(唐国庆.高超声速多体分离动态气动加热计算技 术[D].南京:南京航空航天大学,2013.),可以得到电子的摩尔密度:
考虑高温化学非平衡效应下的气动加热计算方法可修正为:首先由Fay-Riddell驻点 热流密度计算公式(1)得到平衡条件下热流密度qeq,然后根据化学反应计算模型得到各组元摩尔密度和质量分数,最后由驻点化学非平衡边界层与平衡边界层热流密度关系式(10)求得高温非平衡条件下的热流密度qne
本发明的有益效果:
本发明依据边界层理论,采用CFD(Computational Fluid Dynamics)数值计算与工 程算法相结合的方法,将气动热的计算简化为绕飞行器的无粘外流数值解和边界层内热 流求解两个部分,同时耦合防热系统结构传热计算模型,发展了可用于快速计算三维复杂外形高速飞行器气动加热与结构传热特性的计算方法,实现了复杂飞行条件下飞行器全机表面热流密度与防热结构温度场时变特性的计算。
该方法的优点在于:综合了全流场数值模拟与工程算法各自的优势,可以快速高效 对三维复杂外形高速飞行器的多种飞行状态下气动加热与结构耦合传热特性进行计算与分析,给出热流密度与温度等参数分布,弥补了直接数值模拟方法计算代价高、效率 低、周期长等缺陷,同时拓展了工程算法的应用范围。
附图说明
图1a为沿子午线的组元摩尔浓度分布(t=1000s)。
图1b沿子午线表面温度分布(t=1000s)。
图1c为表面温度:不考虑化学非平衡效应(t=1000s)。
图1d为表面温度:考虑化学非平衡效应(t=1000s)。
图2为计算模型及物面网格和TPS方案。
图3a为第1000s时得飞行器防热层外表面温度分布。
图3b为第1000s时得飞行器防热层内表面温度分布。
图3c为飞行器在计算初始时刻(t=0s)的防护层外表面热流密度分布。
图3d为飞行器在计算结束时刻(t=1000s)的防护层外表面热流密度分布。
图3e为驻点内/外层温度随时间的变化曲线。
图3f为驻点热流密度随时间的变化曲线。
图4为弹道状态动态插值方法计算流程图。
图5a为弹道状态t=0s时的表面温度分布。
图5b为弹道状态t=100s时的表面温度分布。
图5c为弹道状态t=350s时的表面温度分布。
图5d为弹道状态t=600s时的表面温度分布。
图5e为弹道状态t=1000s时的表面温度分布。
图6为本发明计算技术总体流场示意图
具体实施方式
为了便于本领域技术人员的理解,下面结合实施例与附图对本发明作进一步的说明,实施方式提及的内容并非对本发明的限定。
1、RAMC-II巡航状态长时加热(高温非平衡效应)
1.1技术参数
计算外形采用典型高速飞行器RAMC-II球锥体试验模型,其几何参数为:头部曲率半径0.152m,半顶角9°,模型长度1.3m。无粘外流场解的计算网格为总单元数约 40万、物面单元约1.5万的块结构网格。设定的巡航状态为:Ma∞=12,迎角α=0°,飞 行高度H=33km,飞行时间t=1000s。
长时加热中考虑结构传热,飞行器热防护结构材料为金属Ti,厚度2mm,表面发 射率为0.8。结构传热计算中,时间步长取为0.05s(即总时间推进步数为2万步),采用 热薄壁传热模型.该算例在普通微机上计算耗时约17minCPU机时。
2.2数值模拟结果
图1给出了该算例数值计算结果。图1a为九组元模型O、N、O2、N2、NO、e-、 N+、O+、NO+在长时加热第1000秒时刻沿子午线分布的组元摩尔浓度分布。由于在本 算例计算状态下空气电离度很低,因此计算中未考虑电离反应对气动热计算的影响。因 此,图中可以看出第1000秒时,带电粒子摩尔浓度都基本为零,沿子午线组元成分主 要以O、N、N2为主,O2基本分解完毕,N2部分分解,计算中O2与N2的分解特性与 空气高温非平衡流理论相符。图1b为子午线上考虑空气高温非平衡效应对比完全气体 的表面温度分布,由图可见两者的温度差异明显,在驻点处化学非平衡效应导致了约 300K的温降,这与理论分析趋势一致。图1c、1d分别给出了第1000秒时不考虑化学 反应效应和考虑化学反应效应时的表面温度分布云图。在本算例条件下多组元的高温化 学反应效应较明显,化学非平衡效应使得壁面热流密度降低,飞行器表面温度下降,这 与理论分析相符。该算例表明本文所发展的计算方法可以用于高温化学非平衡流的气动 加热问题。
2、X-37B巡航状态及弹道状态长时传热
2.1、技术参数
本发明的计算模型选取类X-37B外形,如图2所示,图2同时给出了计算所用的表面网格。计算状态如表1所示。
表1
用于无粘外流求解的块结构网格总单元数约512万,物面网格单元约25.1万。在热防 护结构方面,该算例进行了更接近于工程设计的复杂方案,即飞行器头部区域设定为热防护区1,以TPS1表示,全机其他部位设为热防护区2,以TPS2表示(见图1),不同的热 防护区域使用表2给出的不同的热防护方案,其中热防护结构材料的最外层(飞行器表面) 和最内层温度初值设为飞行高度上的大气环境温度245K。根据该算例的飞行条件,不考 虑高温非平衡效应影响特性。该算例耗时约28minCPU机时。
表2
2.2、数值模拟结果
图3a和3b分别为第1000秒时飞行器防热层外表面(TWsurf)和内表面(TWin)温 度分布,飞行器在计算初始时刻(第0s)和计算结束时(第1000s)的防护层外表面热流 密度分布如图3c、3d所示。由计算结果可见,巡航1000秒时,类X-37B热防护层内 表面温度在不同的热防护区域(TPS)出现明显差异(图3b),而热防护层整个外表面 温度分布均匀。这是由于TPS1区使用的是热薄壁与隔热层的组合防热结构,并且隔热 层采用热沉性好的二氧化硅材料(SiO2),故在该热防护区域内表面温度出现显著降低, 大部分区域温度在410K左右,说明防热层隔热效果良好,与工程实际情况相符。而TPS2 区仅采用了热薄壁结构,防热材料为0.002m的金属Ti,其热传导性很好,故防护层内 外表面温度很快达到平衡,与外表面温度差别很小,且防护区域热流密度很小,该区域 温度大部分在600K以上。图3e、3f分别为驻点内外层温度与热流密度随时间变化曲线, 可以发现TPS1内外表面温度变化差别较大,计算结果显示驻点区域防热层外表面温度 为最高温约为979K,内表面最高温度约为521K,造成此温度值差异的原因仍然是由于 在这个区域使用了SiO2隔热层。
2.3、弹道状态长时加热
在以上算例的基础上,考虑弹道状态,即变高度、迎角、马赫数,弹道状态更加接近飞行器实际飞行包线状态。计算流程如图4所示。
本文根据相关文献设定的弹道参数如表3所示,弹道参数考虑了马赫数、飞行高速及迎角的变化,弹道飞行时间1000s,在此期间飞行高度由60km降至30km,飞行马赫 数由Ma=6降至Ma=4,迎角变化范围为±5°。该算例全机表面采用同一种组合热防 护方案模型,具体参数上述算例中TPS1区相同。该算例耗时约31minCPU机时(根据 弹道飞行参数,不考虑高温非平衡效应)。
表3
图5a~图5e给出了该算例计算结果,从图中可以清晰的看到飞行器表面沿弹道飞行 时逐渐加热到平衡状态的过程,并且随迎角从-5度至5度变化过程中飞行器头部与机翼前缘处的高温区往下表面移动,第1000秒时,最高温出现在飞行器头部下表面,约为1050K,这与高超声速飞行器防热结构材料特性理论分析相符。由于没有在公开发表文 献中找到与此类似的算例进行对比,因此这里仅给出本文算例计算结果的分析。该算例 表明本文所发展的计算方法可对复杂高超声速飞行器在弹道飞行状态下的气动加热与 结构传热耦合进行快速高效计算与分析,可为高速飞行器的热防护设计与选型提高技术 参考。
本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。
Claims (8)
1.一种三维复杂外形高速飞行器流-固-热快速计算方法,其特征在于,包括步骤如下:
步骤1、无粘外流场数值解;
步骤2、工程表面热流率计算;
步骤3、工程结构传热耦合计算;
步骤4、弹道状态动态插值方法;
步骤5、高温化学非平衡效应。
2.如权利要求1所述的三维复杂外形高速飞行器流-固-热快速计算方法,其特征在于:
步骤1具体包括:
控制方程采用三维欧拉方程对无粘外流场求解,采用块结构网格机翼基于分布式并行计算技术,取无粘流场解结果中的物面参数作为边界层外缘参数提供给工程方法中的边界层简化算法,输出的物面参数包括:坐标(m)、物面速度分量(m/s)、马赫数、压强(Pa)、密度(kg/m3)、静温(K)、总能(J/m3),而边界层的工程算法可以得到用于防热系统中结构传热计算所需的飞行器表面热流及传热系数等参数。
3.如权利要求1所述的三维复杂外形高速飞行器流-固-热快速计算方法,其特征在于:
步骤2具体包括:
本发明针对的是全机外形,根据热流密度工程计算方法将飞行器表面热流的计算划分为驻点区和非驻点区两个区域。驻点区热流计算采用Fay-Riddell公式:
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0.763</mn>
<mo>*</mo>
<msup>
<mi>Pr</mi>
<mrow>
<mo>-</mo>
<mn>0.6</mn>
</mrow>
</msup>
<mo>*</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>w</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>0.1</mn>
</msup>
<mo>*</mo>
<msqrt>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>s</mi>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>du</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>s</mi>
</msub>
</mrow>
</msqrt>
<mo>*</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>Le</mi>
<mn>0.52</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<msub>
<mi>h</mi>
<mi>D</mi>
</msub>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>w</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:ρw、μw、hw分别表示物面密度、物面粘性系数和物面焓,ρs、μs分别为驻点密度及驻点粘性系数,hD为平均空气电离焓,计算中取普朗特数Pr=0.71、路易斯数Le=1.0;
非驻点区的热流密度计算公式采用的是工程中常用的平板传热模型。
4.如权利要求1所述的三维复杂外形高速飞行器流-固-热快速计算方法,其特征在于:
步骤3具体包括:
热防护结构内表面采用绝热壁边界条件,依据防热结构材料的毕奥数Bi的大小,采不同的传热计算模型:
<mrow>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&alpha;</mi>
<mi>&delta;</mi>
</mrow>
<msub>
<mi>&lambda;</mi>
<mi>&delta;</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中α为传热系数,δ为材料结构厚度,λδ为当前材料热传导系数。
当Bi<0.1时,采用如下式所示的热薄壁传热模型:
<mrow>
<msub>
<mi>&rho;</mi>
<mi>&delta;</mi>
</msub>
<msub>
<mi>c</mi>
<mi>&delta;</mi>
</msub>
<mi>&delta;</mi>
<mfrac>
<mrow>
<mi>d</mi>
<mi>T</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>&alpha;</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>&epsiv;&sigma;T</mi>
<mi>w</mi>
<mn>4</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
初始条件为:
Tw|t=0=T0 (4)
式中ρδ、cδ、ε分别为表示分别表示防热材料的密度、比热容和表面辐射系数。采用差分方法对式(3)进行加热时间t推进求解,即可得到热薄壁条件下热防护层温度随时间的变化。
当Bi>0.1时,采用热厚壁传热模型。将热厚壁由内而外分为j层进行逐层推进求解:
表层:
最内层:
中间层:
初始条件为:
Tj|t=0=T1|t=0=Tn|t=0=T0 (8)
式中下标m为为材料类型,n表示层数,λm为材料m的热传导系数。利用差分离散可计算求得热厚壁条件下热防护结构各层材料的温度随时间变化的结果。
5.如权利要求4所述的三维复杂外形高速飞行器流-固-热快速计算方法,其特征在于:
气动加热与结构传热过程的耦合特性以物面温度Tw作为物面边界输入条件,然后通过简化边界层工程方法计算得到表面热流密度qn;而在计算结构传热时,又以热流密度qn作为热边界条件,计算得到物面温度。在时间步tn的求解中,以上一时间步tn-1的物面温度作为本时间步的物面温度边界条件,采用工程算法进行边界层的气动加热计算,得出tn时间步的热流密度然后以作为热边界输入条件提供给结构传热计算模块,计算得到tn时间步的物面温度如此循环以实现结构传热和气动热环境的耦合计算。
6.如权利要求1所述的三维复杂外形高速飞行器流-固-热快速计算方法,其特征在于:
步骤4具体包括:
计算非定常弹道状态的气动加热与结构传热耦合计算需要足够多的无粘外流的数值解作为其输入参数;
利用无粘外流解动态插值方法提高耦合算法计算效率:通过已经预先完成的有限个无粘外流解的流场结果,插值得到当前弹道时间点上飞行条件下的流场解,以供气动加热与结构传热计算所需。
7.如权利要求6所述的三维复杂外形高速飞行器流-固-热快速计算方法,其特征在于:
1)根据飞行器物面上的当地流场参数与来流参数的比值在不同飞行高度下几乎保持不变这一规律,对同一飞行器,已经某高度H0下的边界层外缘参数P0以及该高度下的流动参数Pinf0,同时已知任意高度下的大气参数Pinfx,可以通过Px=Pinfx×P0/Pinf0的简单换算获得该飞行器在任意高度上的边界层外缘参数;
2)由1),对于整个弹道状态,无粘外流解只需计算设定参考飞行高度下的两个马赫数和两个迎角,即4个飞行状态,其他飞行状态则可根据插值方法得到。这四个计算状态分别为设定飞行高度下的两个马赫数和两个迎角的组合,即:(Ma1,α1),(Ma1,α2),(Ma2,α1),(Ma2,α2)。与这四个状态所对应的飞行器表面任一点处的流动参数的值分别记作P11、P12、P21、P22,则某个弹道时间点上飞行条件(Max,αx)状态下的飞行器表面同一位置处的流动参数Pxx可通过如下插值算法得到:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>x</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mn>11</mn>
</msub>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>Ma</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Ma</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>P</mi>
<mn>21</mn>
</msub>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>M</mi>
<mi>x</mi>
</msub>
<mo>-</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>Ma</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Ma</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>x</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mn>12</mn>
</msub>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>Ma</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Ma</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>P</mi>
<mn>22</mn>
</msub>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>M</mi>
<mi>x</mi>
</msub>
<mo>-</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>Ma</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>Ma</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>x</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>x</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>x</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
8.如权利要求1所述的三维复杂外形高速飞行器流-固-热快速计算方法,其特征在于:
步骤5具体包括:
在实际高超声速飞行条件下,空气的高温非平衡效应对气动热影响显著,高温空气非平衡边界层驻点传热与平衡边界层驻点传热表面热通量的关系式表示为:
<mrow>
<msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>q</mi>
<mrow>
<mi>N</mi>
<mi>e</mi>
</mrow>
</msub>
<msub>
<mi>q</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Le&phi;</mi>
<mi>O</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msubsup>
<mi>h</mi>
<mi>O</mi>
<mn>0</mn>
</msubsup>
<msub>
<mi>C</mi>
<mrow>
<mi>O</mi>
<mo>,</mo>
<mi>s</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Le&phi;</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msubsup>
<mi>h</mi>
<mi>N</mi>
<mn>0</mn>
</msubsup>
<msub>
<mi>C</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>s</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mi>e</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mfrac>
<msub>
<mi>h</mi>
<mi>D</mi>
</msub>
<mrow>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mn>0.5</mn>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,下标O、N分别表示氧原子和氮原子,φ为壁面催化因子;h0为生成焓;CO.s和CN.s分别表示氧原子和氮原子的质量浓度;Le为路易斯数,取值范围1~2。对于组元参数特性,使用组元参数简化计算模型,在流场温度9000K以下时,只考虑O2、O、O+、N2、N、N+、NO、NO+、e-等9种组元及以下6个化学反应方程式:
①②
③④
⑤⑥
其中Ki为摩尔密度平衡常数。本发明暂不考虑NO、NO+、e-、N+、O+。
氧原子和氮原子的壁面催化因子有如下关系成立:
<mrow>
<msub>
<mi>&phi;</mi>
<mi>O</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>N</mi>
<mi>O</mi>
<mo>&rsqb;</mo>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<msup>
<mi>NO</mi>
<mo>+</mo>
</msup>
<mo>&rsqb;</mo>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<msup>
<mi>O</mi>
<mo>+</mo>
</msup>
<mo>&rsqb;</mo>
</mrow>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>O</mi>
<mn>2</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>0</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mo><</mo>
<mo><</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&phi;</mi>
<mi>N</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>N</mi>
<mi>O</mi>
<mo>&rsqb;</mo>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<msup>
<mi>NO</mi>
<mo>+</mo>
</msup>
<mo>&rsqb;</mo>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<msup>
<mi>N</mi>
<mo>+</mo>
</msup>
<mo>&rsqb;</mo>
</mrow>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>0</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mo><</mo>
<mo><</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,[O2]0,[N2]0分别表示为空气中氧分子和氮分子的摩尔密度。联立式(11)、(12)与反应方程①、②求解,得到氧原子和氮原子的摩尔密度:
<mrow>
<mo>&lsqb;</mo>
<mi>O</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mn>4</mn>
</mfrac>
<mo>&lsqb;</mo>
<mo>-</mo>
<mn>1</mn>
<mo>+</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mn>16</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mi>O</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>O</mi>
<mn>2</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
</mfrac>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&lsqb;</mo>
<mi>N</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mn>4</mn>
</mfrac>
<mo>&lsqb;</mo>
<mo>-</mo>
<mn>1</mn>
<mo>+</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mn>16</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
取一级近似φO=φN=0。联立反应方程①②③,可得到[O2]、[N2]和[NO]的摩尔密度计算式为:
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>O</mi>
<mn>2</mn>
</msub>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>O</mi>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>,</mo>
<mo>&lsqb;</mo>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>N</mi>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>,</mo>
<mo>&lsqb;</mo>
<mi>N</mi>
<mi>O</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mi>N</mi>
<mo>&rsqb;</mo>
<mo>&lsqb;</mo>
<mi>O</mi>
<mo>&rsqb;</mo>
</mrow>
<msub>
<mi>K</mi>
<mn>3</mn>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
联立反应方程④⑤⑥和电荷守恒方程,可以得到电子的摩尔密度:
<mrow>
<mo>&lsqb;</mo>
<msup>
<mi>e</mi>
<mo>-</mo>
</msup>
<mo>&rsqb;</mo>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>K</mi>
<mn>4</mn>
</msub>
<mo>&lsqb;</mo>
<mi>N</mi>
<mi>O</mi>
<mo>&rsqb;</mo>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>5</mn>
</msub>
<mo>&lsqb;</mo>
<mi>O</mi>
<mo>&rsqb;</mo>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>6</mn>
</msub>
<mo>&lsqb;</mo>
<mi>N</mi>
<mo>&rsqb;</mo>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
考虑高温化学非平衡效应下的气动加热计算方法可修正为:首先由Fay-Riddell驻点热流密度计算公式(1)得到平衡条件下热流密度qeq,然后根据化学反应计算模型得到各组元摩尔密度和质量分数,最后由驻点化学非平衡边界层与平衡边界层热流密度关系式(10)求得高温非平衡条件下的热流密度qne。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711177571.2A CN107808065B (zh) | 2017-11-23 | 2017-11-23 | 三维复杂外形高速飞行器流-固-热快速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711177571.2A CN107808065B (zh) | 2017-11-23 | 2017-11-23 | 三维复杂外形高速飞行器流-固-热快速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107808065A true CN107808065A (zh) | 2018-03-16 |
CN107808065B CN107808065B (zh) | 2019-12-31 |
Family
ID=61581312
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711177571.2A Active CN107808065B (zh) | 2017-11-23 | 2017-11-23 | 三维复杂外形高速飞行器流-固-热快速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107808065B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108629089A (zh) * | 2018-04-17 | 2018-10-09 | 西南交通大学 | 一种基于三维连续-离散单元耦合的强夯加固地基模拟方法 |
CN109029907A (zh) * | 2018-07-18 | 2018-12-18 | 大连理工大学 | 一种气动热环境试验模拟条件的参数相似方法 |
CN109583067A (zh) * | 2018-11-22 | 2019-04-05 | 北京空天技术研究所 | 基于温度平衡的高速飞行器转捩位置测量传感器设计方法 |
CN109726432A (zh) * | 2018-11-26 | 2019-05-07 | 北京空天技术研究所 | 飞行器底部结构温度计算方法 |
CN109918808A (zh) * | 2019-03-13 | 2019-06-21 | 北京宇航系统工程研究所 | 一种气热弹三场耦合仿真分析方法 |
CN110119535A (zh) * | 2019-04-11 | 2019-08-13 | 上海交通大学 | 一种轴对称气固耦合传热模型、分析方法及应用系统 |
CN111780949A (zh) * | 2020-07-10 | 2020-10-16 | 南京航空航天大学 | 基于cfd分析的高速进气道前体风洞实验总压修正方法 |
CN111859532A (zh) * | 2020-06-16 | 2020-10-30 | 空气动力学国家重点实验室 | 考虑高超声速化学非平衡效应的改进热壁修正方法 |
CN112329148A (zh) * | 2020-11-18 | 2021-02-05 | 西北工业大学 | 一种高超声速飞行器热防护结构优化方法及系统 |
CN113065275A (zh) * | 2021-03-05 | 2021-07-02 | 浙江大学 | 一种飞行器飞行过程中驻点热流的在线预示方法 |
CN113867381A (zh) * | 2021-12-02 | 2021-12-31 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种飞行器姿态控制方法 |
CN114021499A (zh) * | 2021-11-05 | 2022-02-08 | 沈阳飞机设计研究所扬州协同创新研究院有限公司 | 基于fvm-tlbfs方法的飞行器热防护结构热传导计算方法 |
CN114186330A (zh) * | 2021-11-04 | 2022-03-15 | 北京机电工程研究所 | 高速流场环境下轴类结构温度估算方法 |
CN114626313A (zh) * | 2022-03-04 | 2022-06-14 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种可解析时变热响应的高速气动热cfd求解方法 |
CN116341397A (zh) * | 2023-05-30 | 2023-06-27 | 湖南光华防务科技集团有限公司 | 一种基于深度神经网络的灭火弹温度场模拟方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
WO2013078628A1 (zh) * | 2011-11-30 | 2013-06-06 | 天津空中代码工程应用软件开发有限公司 | 直升机旋翼飞行结冰的数值模拟方法 |
CN103366052A (zh) * | 2013-06-27 | 2013-10-23 | 中国航天空气动力技术研究院 | 一种高超声速飞行器热静气弹分析方法 |
CN106528990A (zh) * | 2016-10-27 | 2017-03-22 | 中国运载火箭技术研究院 | 一种基于泛函优化的高超声速尖锥外形热流密度建模方法 |
CN106682392A (zh) * | 2016-11-24 | 2017-05-17 | 南京航空航天大学 | 复杂高超声速飞行器烧蚀效应快速计算技术 |
-
2017
- 2017-11-23 CN CN201711177571.2A patent/CN107808065B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
WO2013078628A1 (zh) * | 2011-11-30 | 2013-06-06 | 天津空中代码工程应用软件开发有限公司 | 直升机旋翼飞行结冰的数值模拟方法 |
CN103366052A (zh) * | 2013-06-27 | 2013-10-23 | 中国航天空气动力技术研究院 | 一种高超声速飞行器热静气弹分析方法 |
CN106528990A (zh) * | 2016-10-27 | 2017-03-22 | 中国运载火箭技术研究院 | 一种基于泛函优化的高超声速尖锥外形热流密度建模方法 |
CN106682392A (zh) * | 2016-11-24 | 2017-05-17 | 南京航空航天大学 | 复杂高超声速飞行器烧蚀效应快速计算技术 |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108629089A (zh) * | 2018-04-17 | 2018-10-09 | 西南交通大学 | 一种基于三维连续-离散单元耦合的强夯加固地基模拟方法 |
CN108629089B (zh) * | 2018-04-17 | 2021-08-03 | 西南交通大学 | 一种基于连续-离散单元耦合的强夯加固地基模拟方法 |
WO2020015375A1 (zh) * | 2018-07-18 | 2020-01-23 | 大连理工大学 | 一种气动热环境试验模拟条件的参数相似方法 |
CN109029907A (zh) * | 2018-07-18 | 2018-12-18 | 大连理工大学 | 一种气动热环境试验模拟条件的参数相似方法 |
US11454566B2 (en) | 2018-07-18 | 2022-09-27 | Dalian University Of Technology | Parameter similarity method for test simulation conditions of aerodynamic heating environment |
CN109583067A (zh) * | 2018-11-22 | 2019-04-05 | 北京空天技术研究所 | 基于温度平衡的高速飞行器转捩位置测量传感器设计方法 |
CN109726432B (zh) * | 2018-11-26 | 2023-05-16 | 北京空天技术研究所 | 飞行器底部结构温度计算方法 |
CN109726432A (zh) * | 2018-11-26 | 2019-05-07 | 北京空天技术研究所 | 飞行器底部结构温度计算方法 |
CN109918808A (zh) * | 2019-03-13 | 2019-06-21 | 北京宇航系统工程研究所 | 一种气热弹三场耦合仿真分析方法 |
CN109918808B (zh) * | 2019-03-13 | 2023-10-20 | 北京宇航系统工程研究所 | 一种气热弹三场耦合仿真分析方法 |
CN110119535B (zh) * | 2019-04-11 | 2023-05-26 | 上海交通大学 | 一种轴对称气固耦合传热模型、分析方法及应用系统 |
CN110119535A (zh) * | 2019-04-11 | 2019-08-13 | 上海交通大学 | 一种轴对称气固耦合传热模型、分析方法及应用系统 |
CN111859532B (zh) * | 2020-06-16 | 2023-11-28 | 空气动力学国家重点实验室 | 考虑高超声速化学非平衡效应的改进热壁修正方法 |
CN111859532A (zh) * | 2020-06-16 | 2020-10-30 | 空气动力学国家重点实验室 | 考虑高超声速化学非平衡效应的改进热壁修正方法 |
CN111780949A (zh) * | 2020-07-10 | 2020-10-16 | 南京航空航天大学 | 基于cfd分析的高速进气道前体风洞实验总压修正方法 |
CN112329148A (zh) * | 2020-11-18 | 2021-02-05 | 西北工业大学 | 一种高超声速飞行器热防护结构优化方法及系统 |
CN113065275B (zh) * | 2021-03-05 | 2022-05-24 | 浙江大学 | 一种飞行器飞行过程中驻点热流的在线预示方法 |
CN113065275A (zh) * | 2021-03-05 | 2021-07-02 | 浙江大学 | 一种飞行器飞行过程中驻点热流的在线预示方法 |
CN114186330A (zh) * | 2021-11-04 | 2022-03-15 | 北京机电工程研究所 | 高速流场环境下轴类结构温度估算方法 |
CN114021499A (zh) * | 2021-11-05 | 2022-02-08 | 沈阳飞机设计研究所扬州协同创新研究院有限公司 | 基于fvm-tlbfs方法的飞行器热防护结构热传导计算方法 |
CN114021499B (zh) * | 2021-11-05 | 2024-06-11 | 沈阳飞机设计研究所扬州协同创新研究院有限公司 | 基于fvm-tlbfs方法的飞行器热防护结构热传导计算方法 |
CN113867381A (zh) * | 2021-12-02 | 2021-12-31 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种飞行器姿态控制方法 |
CN114626313A (zh) * | 2022-03-04 | 2022-06-14 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种可解析时变热响应的高速气动热cfd求解方法 |
CN116341397A (zh) * | 2023-05-30 | 2023-06-27 | 湖南光华防务科技集团有限公司 | 一种基于深度神经网络的灭火弹温度场模拟方法 |
CN116341397B (zh) * | 2023-05-30 | 2023-09-05 | 湖南光华防务科技集团有限公司 | 一种基于深度神经网络的灭火弹温度场模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107808065B (zh) | 2019-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107808065A (zh) | 三维复杂外形高速飞行器流‑固‑热快速计算方法 | |
Tsien | Superaerodynamics, mechanics of rarefied gases | |
Xiao et al. | Numerical investigation of transpiration cooling for porous nose cone with liquid coolant | |
Yen et al. | Synthetic jets as a boundary vorticity flux control tool | |
Xiong et al. | Numerical simulation of transpiration cooling for sintered metal porous strut of the scramjet combustion chamber | |
Sebastian et al. | Numerical analysis of hypersonic flow past blunt bodies with aerospikes | |
CN106682392A (zh) | 复杂高超声速飞行器烧蚀效应快速计算技术 | |
CN107944137A (zh) | 高超声速飞行器弹道状态多场耦合的热气动弹性计算技术 | |
CN109800488A (zh) | 关于液体火箭高空环境下底部热环境的数值计算方法 | |
Wei | The development and application of CFD technology in mechanical engineering | |
Jiang et al. | Optimization on conventional and electric air-cycle refrigeration systems of aircraft: A short-cut method and analysis | |
Fan et al. | Prediction of hypersonic boundary layer transition with variable specific heat on plane flow | |
Huang et al. | Coupled Fluid–Structural Thermal Numerical Methods for Thermal Protection System | |
Ferrero et al. | A numerical method for conjugate heat transfer problems in hypersonic flows | |
Mei et al. | Coupled simulation for reentry ablative behavior of hypersonic vehicles | |
Aogaki et al. | High angle-of-attack pitching moment characteristics of slender-bodied reusable rocket | |
Wang et al. | Research on the cooling performance of the discontinuous transpiration surface structure for the leading edge of a hypersonic vehicle | |
Li et al. | An engineering method of aerothermodynamic environments prediction for complex reentry configurations | |
Su et al. | Fluid–Thermal–Structure Interaction of Hypersonic Inlets Under Different Aspect Ratios | |
Meng et al. | A hypersonic aeroheating calculation method based on inviscid outer edge of boundary layer parameters | |
Joseph et al. | Coolant gas injection on a blunt-nosed re-entry vehicle | |
Saha et al. | Numerical prediction of surface heat flux during multiple jets firing for missile control | |
Wang et al. | Variation character of stagnation point heat flux for hypersonic pointed bodies from continuum to rarefied flow states and its bridge function study | |
Clarey et al. | Three-Temperature Thermochemical Nonequilibrium Model with Application to Slender-Body Wakes | |
Viviani et al. | Analysis of flowfield past a sample return capsule with heatshield ablation |
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 |