CN109960878A - 一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法 - Google Patents

一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法 Download PDF

Info

Publication number
CN109960878A
CN109960878A CN201910229264.7A CN201910229264A CN109960878A CN 109960878 A CN109960878 A CN 109960878A CN 201910229264 A CN201910229264 A CN 201910229264A CN 109960878 A CN109960878 A CN 109960878A
Authority
CN
China
Prior art keywords
thermal protection
active
protection system
passive thermal
passive
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
Application number
CN201910229264.7A
Other languages
English (en)
Other versions
CN109960878B (zh
Inventor
苟建军
龚春林
陈兵
王健磊
谷良贤
胡嘉欣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201910229264.7A priority Critical patent/CN109960878B/zh
Publication of CN109960878A publication Critical patent/CN109960878A/zh
Application granted granted Critical
Publication of CN109960878B publication Critical patent/CN109960878B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Laminated Bodies (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,首先以飞行器构型和网格、任务弹道以及主动冷却参数为输入,基于热平衡方程,利用工程算法求解获得气动热;其次选择热防护概念得到整个高超声速飞行器的被动热防护概念分布;之后确定被动热防护系统的规模包括厚度和质量,以及主动热防护系统的容量包括冷却工质的质量流量和允许的温升;最后进行评估以及迭代设计直至满足飞行器总体设计需求。本发明将主动冷却系统置于被动热防护系统壁面之下,从而在进行飞行器表面气动热计算、被动热防护概念分布以及被动热防护系统规模计算时,考虑了主动冷却的影响,实现了飞行器被动和主动热防护系统的耦合设计。

Description

一种面向高超声速飞行器总体的主/被动热防护系统耦合设 计方法
技术领域
本发明涉及飞行器热防护系统设计领域,特别是涉及一种面向高超声速飞行器的主动冷却与被动热防护系统耦合设计方法。
背景技术
目前,随着高超声速工程的发展,特别是航天飞行器可重复使用、宽包线飞行和长时间高速巡航需求的出现,飞行器热环境非常恶劣,单一的被动热防护技术无法满足防热需求,因此需要开发被动与主动相结合的热防护系统。被动热防护系统通常由不同的热防护概念组成,如NASA基于航天飞机工程研发了一系列金属、陶瓷防热瓦和陶瓷隔热毡(D EMyers,C J Martin,and M L Blosser,Parametric Weight Comparison of AdvancedMetallic,Ceramic Tile,and Ceramic Blanket Thermal Protection Systems.2000,NASA Langley Technical Report Server.),依靠防热概念的防、隔热性能阻隔气动热进入机体;而主动热防护系统可由冷却工质的对流网络等组成,将热量从高温区域带至低温区域,以保护飞行器安全。主动冷却会影响飞行器壁面的气动加热,进而影响被动热防护概念的分布、尺寸和质量,但目前的研究中,主动和被动热防护系统间的设计通常耦合程度很低,很难获得飞行器热防护系统的最优设计区间。
发明内容
为解决现有技术存在的问题,本发明提出一种面向高超声速飞行器总体的主动与被动热防护系统的耦合设计方法,通过建立等效热平衡模型,即将主动冷却系统置于被动热防护系统壁面之下,并利用等效换热系数表征主动冷却的强度,通过迭代设计,实现了面向高超声速飞行器总体的被动与主动热防护系统的耦合设计。
本发明的技术方案为:
所述一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:包括以下步骤:
步骤1:以高超声速飞行器的三维构型及划分网格、任务弹道和主动冷却参数为输入参数,基于热平衡方程利用工程算法求解获得高超声速飞行器气动热,输出气动热载、壁面温度和主动换热量;所述主动冷却参数包括等效换热系数、冷却工质温度;
步骤2:根据步骤1中得出的壁面温度,确定分析区域在全弹道下的最大壁面温度;以分析区域在全弹道下的最大壁面温度为输入条件,从事先建立的热防护概念数据库中为高超声速飞行器壁面的分析区域选择热防护概念,要求该分析区域在全弹道下的最大壁面温度小于所选择热防护概念的许用温度,进而得到整个高超声速飞行器的被动热防护概念分布;
步骤3:针对步骤2确定的被动热防护概念建立一维非稳态传热数值模型,并以步骤1得到的气动热载和主动换热量分别为外、内边界条件进行传热分析,计算被动热防护概念的温度分布;而后以内边界的设定温度和被动热防护概念各层材料的许用温度为约束,,以被动热防护系统质量最轻为目标,以被动热防护概念中隔热层的厚度为设计参数,进行结构优化,计算获得最优的被动热防护系统的厚度和质量;以冷却工质的比热以及步骤1得到的主动换热量为输入,计算获得冷却工质的质量流量和允许的温升,
步骤4:判断步骤3得到的被动热防护系统的质量和厚度、主动热防护系统所需的冷却工质质量流量是否满足高超声速飞行器总体设计约束,如果满足,则当前系统可作为飞行器的主/被动热防护系统;若不满足,则修改主动冷却参数,重复步骤1至3,进行新一轮的设计,如此迭代直至最终满足高超声速飞行器总体设计约束。
进一步的优选方案,所述一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:步骤1中任务弹道包括飞行时间、高度、速度、攻角和侧滑角。
进一步的优选方案,所述一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:步骤1中所述热平衡方程为气动加热量等于辐射耗散量与主动换热量之和;所述主动换热量等于壁面温度与冷却工质温度之差乘以等效换热系数。
进一步的优选方案,所述一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:步骤1中采用的工程算法基于普朗特边界层理论,边界层外视作理想气体,边界层内为有粘气体,通过边界层外的计算获得边界层外气动参数,之后利用参考温度法计算获得边界层内的气动热。
进一步的优选方案,所述一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:步骤2中热防护概念数据库包括热防护概念各层的密度、最大许用温度、辐射发射率、导热系数、比热和初始厚度。
有益效果
本发明在进行飞行器气动热计算、被动热防护概念分布确定、被动热防护系统质量和厚度计算时,耦合考虑了主动冷却的影响,并可同时输出主动热防护系统所需的工质质量流量和工质温升,实现了被动与主动热防护系统间的耦合关系。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为主/被动热防护系统耦合设计流程。
图2为热防护系统热平衡示意图。
图3为等效的热防护系统热平衡示意图。
图4为一维平板有限差分模型。
具体实施方式
目前国内外飞行器被动热防护系统设计时,无法耦合考虑主动冷却的影响,本发明针对这一问题,提出了一种被动与主动热防护系统的耦合设计方法。实际工程中,主动冷却位于被动热防护系统的内边界底部,本发明建立等效热平衡模型,即将主动冷却系统置于被动热防护系统壁面之下,从而在进行飞行器表面气动热计算、被动热防护概念分布以及被动热防护系统规模计算时,考虑了主动冷却的影响,实现了飞行器被动和主动热防护系统的耦合设计。
本发明的技术方案包括四个步骤,如图1所示:
第一步为计算气动热。以飞行器构型和网格,任务弹道包括飞行时间、高度、速度、攻角和侧滑角,以及主动冷却参数包括等效换热系数、冷却工质温度,为输入参数,基于热平衡方程,即气动加热量等于辐射耗散量与主动换热量之和,利用工程算法求解获得气动热。本发明中,主动换热量等于壁面温度与冷却工质温度之差乘以等效换热系数,实际热平衡模型与等效热平衡模型分别见图2和图3。本发明中工程算法基于普朗特边界层理论,边界层外视作理想气体,边界层内为有粘气体,通过边界层外的计算获得边界层外缘气动参数,进而利用参考温度法计算获得边界层内的气动热。本步将输出气动热载,即热流密度对时间的积分,壁面温度,主动换热量即主动冷却带走的热流密度。
这一步基于工程算法求解飞行器表面热平衡方程获得气动热。飞行器表面平衡方程为:
qaero=α(Tr-Tw)=εσTw 4+qin
式中qaero为气动热热流密度,α为飞行器表面对流换热系数,Tr为气体恢复温度,Tw为飞行器表面温度,ε为飞行器表面黑度,σ为斯特藩-玻尔兹曼常数。本发明中的方法在计算气动加热时,添加了向飞行器舱内传递的热量项qin,即为主动冷却带走的热流密度。实际工程中,主动冷却系统位于被动热防护系统的内边界即冷端底部,系统的热量平衡示意图如图2所示,主动换热项qin为:
qin=h(Tb-Tc)
式中h、Tb以及Tc分别为主动冷却系统与被动热防护结构间的换热系数、被动热防护系统冷端许用温度以及冷却工质的温度。主/被动耦合的热防护系统,主动换热的强度应该与飞行器表面气动加热的强度相对应,即在壁面温度Tw高的区域,加大主动换热的强度,反之则减弱换热强度。但上式中,qin与Tw无关,若将上式带入表面平衡方程,需要以h为设计参数,在气动加热强度不同的区域,需要不同的h,因此,使得设计非常复杂,无法顺利完成。针对该问题,本发明假定主动冷却位于被动热防护系统的壁面之下,如图3所示,此时,主动换热项qin等于:
qin=he(Tw-Tc)
式中,qin与Tw有关,此时等效换热系数he可作为设计变量,对于整个飞行器表面,相同的he即可以模拟气动加热强度不同区域的主动换热项。
等效换热系数he与实际的换热系数h之间的关系由下式确定:
基于上述模型和假设,以弹道、飞行器构型和等效换热系数为输入条件,利用普朗特边界层理论、修正牛顿法及参考温度法构建工程算法程序,通过计算即可获得气动加热热载(热流密度对时间的积分)、壁面温度以及主动冷却换热量。
第二步为确定被动热防护概念的分布。本实施例中以NASA研制的金属防热瓦、陶瓷防热瓦和陶瓷隔热毡构建热防护概念数据库,包含热防护概念各层的密度、最大许用温度、辐射发射率、导热系数、比热和初始厚度。对于飞行器表面某一区域而言,以第一步中得出的全弹道下该区域的最大壁面温度为输入条件,从热防护概念数据库中为飞行器壁面区域选择热防护概念,基本原则是,该区域的最大壁面温度应小于热防护概念的许用温度。针对飞行器表面所有区域确定合适的热防护概念,最终获得飞行器被动热防护概念的分布。
第三步为确定被动热防护系统的规模包括厚度和质量,以及主动热防护系统的容量包括冷却工质的质量流量和允许的温升。针对第二步确定的被动热防护概念建立一维非稳态传热数值模型,如图4所示,以第一步输出的气动热载和主动换热量分别为外、内边界条件,进行传热分析,计算获得被动热防护概念的温度分布即各层材料及边界的温度分布。然后,以内边界即冷端的设定温度Tb和各层材料的许用温度为约束,以热防护系统质量最轻为目标,以被动热防护概念中隔热层的厚度为设计参数,进行结构优化,最终计算获得最优的被动热防护系统的厚度和质量。以第一步输出的主动换热量和冷却工质的比热为输入,计算获得工质的质量流量和允许的温升,即为主动热防护系统的容量。
这一步是计算被动热防护系统规模和主动热防护系统的容量。针对被动热防护系统,建立一维非稳态热传导的有限差分数值模型,以热防护概念的冷端设计温度及各层材料最大许用温度为约束,通过数值分析和优化计算获得飞行器热防护概念的规模,即厚度和质量。热防护概念由多层材料组成,基于其结构特点,建立一维多层平板传热的有限差分模型如图4所示。分析节点沿厚度方向选取,相邻节点间的中垂线为界面,相邻界面间的区域即为当前节点所代表的控制体积。整个计算区域的离散,需要按照先节点位置、后界面位置的方法,最终得到多个子区域组成的空间区域。此外,由于各层材料的厚度和物性互不相同,各层内部需要进行细致离散,而单元的尺寸由结构、温度以及计算规模权衡选取,内外边界以及空隙边界上的单元厚度为其它单元的一半。另外,区域离散还需满足两个条件:热防护概念的内、外表面上必须布置节点,两种材料的交界面必须在单元的边界上。
在区域离散的基础上,则可对控制方程进行离散。对于控制体i(0<i<N),其热平衡方程为:
式中,t为时间,xi为节点的坐标位置,Δi为控制体i的厚度,ρ为该层材料的密度,c为该层材料的比热,k为该层材料的热传导系数。将该式写作差分形式,则有:
考虑到热防护结构的层状特征,可以假设在控制体内,T、ρ、c是关于x阶梯式变化,则可认为控制单元内部的T、ρ、c值即为单元节点处的值,因此:
以隐式格式进行离散,取对时间作隐式阶跃变化,即在t到t+Δt时刻的变化用t+Δt时的数值表征,得到下式:
不同材料层之间界面上的导热系数不连续,则界面上的当量导热系数以调和平均法计算获得,从而有:
其中Wi为热阻系数,取值为:
当节点i位于边界上时δ=1,其余情况取0,最终得到一维平板传热分析模型的数值离散方程为:
重新整理可得:
对于外边界,本文采用第三类边界条件,热平衡方程为:
则第一个控制体的能量平衡方程可写为:
从而得到外边界上的离散方程:
整理后得到:
对于内边界,以主动换热量qin为边界条件,为统一格式,将其写成qin=h(TN-Tc),即第三类边界条件,从而得到内边界控制体上能量平衡方程可写为:
从而内边界离散方程可写为:
整理后得到:
综上所述,得到考虑主动冷却的一维平板瞬态传热有限差分方程:
通过迭代求解,可获得热防护概念结构的温度分布,然后以各层材料的许用温度以及系统冷端的设计温度为约束,以结构质量最低为目标,以隔热层厚度为设计参数,基于序列二次规划方法进行结构优化,最终获得最优的热防护系统。该系统优化问题的数学描述为:
s.t.Ti(x)<Timax,i=1,2,…,N
TN(x)<Tb
式中,n为热防护结构的层数,x为各层的厚度,ρi、Ti(x)、Tb分别为各层材料的密度、最大温度和系统冷端温度,其中Tb为整个系统设计的输入约束。
针对主动热防护系统,可基于主动换热量与工质比热,求解下式即可获得主动热防护系统的容量,即工质流量和温升:
qin=cmfΔT
式中,c为冷却工质的比热,mf为质量流量,△T为工质温升。
第四步为输出飞行器主/被动热防护系统。基于第三步获得的被动和主动热防护系统后,评估是否满足飞行器总体设计需求,即被动热防护系统的质量和尺寸、主动热防护系统所需的工质流量是否满足飞行器总体设计约束,如果满足,则当前系统可作为飞行器的主/被动热防护系统;若不满足,则修改主动冷却参数包括等效换热系数及工质温度,重复步骤一至三,进行新一轮的设计,如此迭代直至最终满足飞行器总体设计需求。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (5)

1.一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:包括以下步骤:
步骤1:以高超声速飞行器的三维构型及划分网格、任务弹道和主动冷却参数为输入参数,基于热平衡方程利用工程算法求解获得高超声速飞行器气动热,输出气动热载、壁面温度和主动换热量;所述主动冷却参数包括等效换热系数、冷却工质温度;
步骤2:根据步骤1中得出的壁面温度,确定分析区域在全弹道下的最大壁面温度;以分析区域在全弹道下的最大壁面温度为输入条件,从事先建立的热防护概念数据库中为高超声速飞行器壁面的分析区域选择热防护概念,要求该分析区域在全弹道下的最大壁面温度小于所选择热防护概念的许用温度,进而得到整个高超声速飞行器的被动热防护概念分布;
步骤3:针对步骤2确定的被动热防护概念建立一维非稳态传热数值模型,并以步骤1得到的气动热载和主动换热量分别为外、内边界条件进行传热分析,计算被动热防护概念的温度分布;而后以内边界的设定温度和被动热防护概念各层材料的许用温度为约束,,以被动热防护系统质量最轻为目标,以被动热防护概念中隔热层的厚度为设计参数,进行结构优化,计算获得最优的被动热防护系统的厚度和质量;以冷却工质的比热以及步骤1得到的主动换热量为输入,计算获得冷却工质的质量流量和允许的温升,
步骤4:判断步骤3得到的被动热防护系统的质量和厚度、主动热防护系统所需的冷却工质质量流量是否满足高超声速飞行器总体设计约束,如果满足,则当前系统可作为飞行器的主/被动热防护系统;若不满足,则修改主动冷却参数,重复步骤1至3,进行新一轮的设计,如此迭代直至最终满足高超声速飞行器总体设计约束。
2.根据权利要求1所述一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:步骤1中任务弹道包括飞行时间、高度、速度、攻角和侧滑角。
3.根据权利要求1所述一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:步骤1中所述热平衡方程为气动加热量等于辐射耗散量与主动换热量之和;所述主动换热量等于壁面温度与冷却工质温度之差乘以等效换热系数。
4.根据权利要求1所述一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:步骤1中采用的工程算法基于普朗特边界层理论,边界层外视作理想气体,边界层内为有粘气体,通过边界层外的计算获得边界层外气动参数,之后利用参考温度法计算获得边界层内的气动热。
5.根据权利要求1所述一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法,其特征在于:步骤2中热防护概念数据库包括热防护概念各层的密度、最大许用温度、辐射发射率、导热系数、比热和初始厚度。
CN201910229264.7A 2019-03-25 2019-03-25 一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法 Expired - Fee Related CN109960878B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910229264.7A CN109960878B (zh) 2019-03-25 2019-03-25 一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910229264.7A CN109960878B (zh) 2019-03-25 2019-03-25 一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法

Publications (2)

Publication Number Publication Date
CN109960878A true CN109960878A (zh) 2019-07-02
CN109960878B CN109960878B (zh) 2022-07-05

Family

ID=67024882

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910229264.7A Expired - Fee Related CN109960878B (zh) 2019-03-25 2019-03-25 一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法

Country Status (1)

Country Link
CN (1) CN109960878B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114491791A (zh) * 2021-12-28 2022-05-13 中国航天空气动力技术研究院 一种用于飞行器的可控烧蚀防热结构设计方法
CN114647260A (zh) * 2022-03-19 2022-06-21 西北工业大学 基于变密度结构微小通道热沉性能调控的飞行器伺服电机主动热控制方法
CN114839864A (zh) * 2022-07-04 2022-08-02 中国飞机强度研究所 空天飞机热强度试验控制系统的辐射干涉解耦控制方法
CN116242878A (zh) * 2023-02-23 2023-06-09 北京天兵科技有限公司 热防护材料性能一体化测试系统及方法
CN118228521A (zh) * 2024-05-24 2024-06-21 西安现代控制技术研究所 一种适应气动热耦合内热的固体发动机复合壳体设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5322725A (en) * 1988-08-31 1994-06-21 The Boeing Company Thermal protection system
CN105956286A (zh) * 2016-05-06 2016-09-21 北京航空航天大学 一种高超声速飞行器前舱热防护系统全弹道温度边界预测方法
CN107792392A (zh) * 2017-09-26 2018-03-13 北京航天长征飞行器研究所 一种飞行器控制舵尖化前缘主动式热疏导试验系统及方法
CN107914862A (zh) * 2017-11-20 2018-04-17 北京临近空间飞行器系统工程研究所 一种全主动冷却高超声速飞行器

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5322725A (en) * 1988-08-31 1994-06-21 The Boeing Company Thermal protection system
CN105956286A (zh) * 2016-05-06 2016-09-21 北京航空航天大学 一种高超声速飞行器前舱热防护系统全弹道温度边界预测方法
CN107792392A (zh) * 2017-09-26 2018-03-13 北京航天长征飞行器研究所 一种飞行器控制舵尖化前缘主动式热疏导试验系统及方法
CN107914862A (zh) * 2017-11-20 2018-04-17 北京临近空间飞行器系统工程研究所 一种全主动冷却高超声速飞行器

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
P. ZHANG等: "Compensation for supply and motor unbalances and line resistance for a DC signal injection-based thermal protection method", 《2009 IEEE 6TH INTERNATIONAL POWER ELECTRONICS AND MOTION CONTROL CONFERENCE》 *
高效伟等: "集成单元边界元法及其在主动冷却热防护系统分析中的应用", 《力学学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114491791A (zh) * 2021-12-28 2022-05-13 中国航天空气动力技术研究院 一种用于飞行器的可控烧蚀防热结构设计方法
CN114647260A (zh) * 2022-03-19 2022-06-21 西北工业大学 基于变密度结构微小通道热沉性能调控的飞行器伺服电机主动热控制方法
CN114647260B (zh) * 2022-03-19 2023-08-29 西北工业大学 基于变密度结构微小通道热沉性能调控的飞行器伺服电机主动热控制方法
CN114839864A (zh) * 2022-07-04 2022-08-02 中国飞机强度研究所 空天飞机热强度试验控制系统的辐射干涉解耦控制方法
CN114839864B (zh) * 2022-07-04 2022-09-13 中国飞机强度研究所 空天飞机热强度试验控制系统的辐射干涉解耦控制方法
CN116242878A (zh) * 2023-02-23 2023-06-09 北京天兵科技有限公司 热防护材料性能一体化测试系统及方法
CN118228521A (zh) * 2024-05-24 2024-06-21 西安现代控制技术研究所 一种适应气动热耦合内热的固体发动机复合壳体设计方法

Also Published As

Publication number Publication date
CN109960878B (zh) 2022-07-05

Similar Documents

Publication Publication Date Title
CN109960878A (zh) 一种面向高超声速飞行器总体的主/被动热防护系统耦合设计方法
Catchpole-Smith et al. Thermal conductivity of TPMS lattice structures manufactured via laser powder bed fusion
CN104267062B (zh) 一种气动热模拟试验中冷壁热流转换为热壁热流的方法
Tang et al. Optimal cooling system design for multi-cavity injection molding
Scozia et al. Natural convection in slender cavities with multiple fins attached to an active wall
Thakur et al. Some implementational issues of convection schemes for finite-volume formulations
CN106570252B (zh) 一种基于概率技术的热防护系统设计方法
CN105843996A (zh) 一种防热结构优化方法
Kong et al. Convective heat transfer enhancement by novel honeycomb-core in sandwich panel exchanger fabricated by additive manufacturing
CN106897520A (zh) 一种含有模糊参数的传热系统可靠性分析方法
Alvarado et al. Interaction between natural convection and surface thermal radiation in tilted slender cavities
Zhao et al. Inverse design of indoor environment using an adjoint RNG k‐ε turbulence model
He et al. An improved unit circuit model for transient heat conduction performance analysis and optimization in multi-layer materials
CN114444413A (zh) 一种板状燃料堆芯亚通道级三维热工水力分析方法
CN105087882A (zh) 一种立式淬火炉热处理阶段划分方法
CN111680441B (zh) 一种适用于热力工况的梯度点阵夹芯板结构
Mouhtadi et al. Natural convection in a horizontal channel provided with heat generating blocks: Discussion of the isothermal blocks validity
CN109598059A (zh) 一种基于代理模型的热防护系统优化设计方法及设计系统
Gawas et al. Axisymmetric thermal-lattice Boltzmann method for Rayleigh-Bénard convection with anisotropic thermal diffusion
Ferrero et al. A numerical method for conjugate heat transfer problems in hypersonic flows
Rong et al. Numerical investigation of heat transfer and parameter coupling characteristics for Spider web microchannel topological structure
Maughan et al. Fully developed mixed convection in a horizontal channel heated uniformly from above and below
Chordiya et al. Numerical analysis of the longitudinal size of the partition on natural convection heat transfer and fluid flow within a differentially heated porous enclosure
Chen et al. Aerothermoelastic analysis of a hypersonic vehicle based on thermal modal reconstruction
Cheng et al. Numerical investigation of channel for helium cooled China HCCB-TBM first wall by application enhance heat transfer method

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220705