CN105095603B - 一种高超声速流动‑传热与结构响应的多场耦合瞬态数值的方法 - Google Patents
一种高超声速流动‑传热与结构响应的多场耦合瞬态数值的方法 Download PDFInfo
- Publication number
- CN105095603B CN105095603B CN201510570860.3A CN201510570860A CN105095603B CN 105095603 B CN105095603 B CN 105095603B CN 201510570860 A CN201510570860 A CN 201510570860A CN 105095603 B CN105095603 B CN 105095603B
- Authority
- CN
- China
- Prior art keywords
- fluid
- displacement
- condition
- hypersonic
- data exchange
- 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
Classifications
-
- 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
技术领域
本发明涉及现代高速飞行器设计技术领域,特别涉及一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法。
背景技术
高超声速飞行器的快速发展给热防护设计带来了更为严峻的挑战。准确的预测气动热/力环境、结构温度和应力状态,能够在提高飞行器安全性能的同时减小热防护系统设计冗余,对提高飞行器的性能有着极为重要的意义。
在现有技术中,传统的高超声速飞行器的热/力载荷环境的预测与热防护结构性能的分析基本还处于分离状态。现有技术中的方法一般都是首先在给定的等温壁面条件下进行流体计算,得到壁面热流或传热系数;然后将得到的热载荷作为边界条件加载到结构上进行热分析得到固体热分布;最后再根据该固体热分布计算得到结构的应力和应变。
现有技术中的上述方法实际上是把多物理场耦合的事实人为的分割成多个独立的物理场,并且也没有考虑各个物理场之间的相互作用。因此,在这种情况下,既无法得到精确的气动热/力载荷环境,也无法正确地评价热防护材料及其结构的服役特征。
高超声速飞行热防护设计是一个涉及到真实气体效应、耦合传热和结构热力响应的复杂的多物理场耦合问题,必须采用多场耦合的方法求解。但是,由于多物理场耦合问题的复杂性,还需要进一步开展分析方法研究,深刻把握防热系统多场耦合规律及其效应。
发明内容
有鉴于此,本发明提供一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法,从而可以实现高超声速非平衡流动求解器与结构热/力全耦合求解器相耦合的多场耦合计算,使得对高超声速飞行器的气动热力环境和结构热力响应的预测更符合物理实际,并可保证计算精度。
本发明的技术方案具体是这样实现的:
一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法,该方法包括:
A、预先建立多物理场耦合模型并设置当前边界条件;
B、根据结构确定壁面温度和位移边界条件,在流体-固体耦合界面进行数据交换,得到流体区域的当前温度和位移边界条件;
C、根据所述流体区域的当前温度和位移边界条件,在流体区域同时求解预设的各个守恒方程的耦合解格式,在计算一个时间步Δt之后,得到当前的热流和压力;
D、根据当前的热流和压力,在流体-固体耦合界面进行数据交换,得到固体区域的边界条件;
E、根据所述固体区域的边界条件,在固体区域通过热力全耦合的方法进行求解,计算一个时间步Δt之后,得到壁面温度Tw和结构位移us;
F、判断是否满足预设的停止条件,如果是,则停止整个流程;否则,返回执行步骤B。
较佳的,所述设置当前边界条件包括:
由结构初始条件确定壁面温度和位移边界条件,进行高超声速稳态流动计算,同时求解预设的各个守恒方程的耦合解格式,得到初始热流和初始压力;
将所述初始热流和初始压力作为当前边界条件。
较佳的,所述数据交换包括:关联和插值。
较佳的,在进行数据交换时:
对于位移和温度采用最近邻居搜索方法计算;
对于压力和热流通量载荷则采用守恒插值方法。
较佳的,所述停止条件为:当前的计算周期为最后一个计算周期。
较佳的,所述各个守恒方程包括:
连续守恒方程、动量守恒方程、能量守恒方程和组分守恒方程。
如上可见,通过使用本发明的高超声速流动-传热与结构响应的多场耦合瞬态数值的方法,可以解决传统方法中热/力载荷环境的预测与热防护结构性能分析互相分离的问题,从而实现高超声速非平衡流动求解器与结构热/力全耦合求解器相耦合的多场耦合计算,使得对高超声速飞行器的气动热力环境和结构热力响应的预测更符合物理实际;而且,由于本发明中考虑了真实气体效应、耦合传热和结构热力全耦合等复杂的效应,从而还可以大大提高高超声速飞行器热/力载荷环境和结构响应预测的计算精度。
附图说明
图1为本发明实施例中的高超声速流动-传热与结构响应的多场耦合瞬态数值的方法的流程示意图。
图2为本发明实施例中的多物理场耦合模型的示意图。
图3为本发明实施例中的耦合策略示意图。
图4为本发明实施例中的圆柱绕流计算模型。
图5为本发明实施例中计算得到的驻点温度随时间变化图。
图6为本发明实施例中计算得到的驻点热流随时间变化图。
图7为本发明实施例中计算得到的表面温度分布。
图8为本发明实施例中计算得到的结构Mises应力和位移图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明。
本实施例提供了一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法。
图1为本发明实施例中的高超声速流动-传热与结构响应的多场耦合瞬态数值的方法的流程示意图。如图1所示,本发明实施例中的高超声速流动-传热与结构响应的多场耦合瞬态数值的方法可以包括如下所述的步骤:
步骤11,预先建立多物理场耦合模型并设置当前边界条件。
在本发明的技术方案中,首先需要预先建立相应的多物理场耦合模型。
在本发明的技术方案中,可以根据实际应用的需要,建立所需的多物理场耦合模型。一般来说,可以使用多种方法来建立所述多物理场耦合。以下将以其中的一种具体实施方式为例,对本发明的技术方案进行详细的介绍。
例如,图2为本发明实施例中的多物理场耦合模型的示意图。如图2所示,本发明中的多物理场耦合模型具体可以包括:用于流体区域的流体分析模型、用于固体区域的热-结构分析模型和用于数据交换的数据交换模型。
在所述用于流体区域的流体分析模型中,主要是利用同时求解连续、动量、能量守恒方程的耦合解格式进行流体区域中的气动热与气动力分析,得到壁面热流qw和壁面压力Pw。
在用于固体区域的热-结构分析模型中,主要是利用热力全耦合方法来进行热-结构分析,得到壁面温度Tw和结构位移us。
在用于数据交换的数据交换模型中,主要是利用关联和插值等方法,将流体分析模型的结果和热-结构分析模型的结果在网格间进行参数交换。
图3为本发明实施例中的耦合策略示意图。如图3所示,在本发明的技术方案中,可以采用分区求解方法完成对高超声速流动-传热和结构响应的耦合分析。其中,在流体区域和固体区域的求解器均为瞬态求解,每个求解器所需要的数据在耦合界面上将进行反复交换。以下,将对上述耦合策略进行详细的介绍。
另外,在本发明的技术方案中,在进行上述耦合策略之前,还需要预先设置当前边界条件。
在本发明的技术方案中,可以根据实际应用的需要,使用多种方法来建预先设置当前边界条件。以下将以其中的一种具体实施方式为例,对本发明的技术方案进行详细的介绍。
较佳的,在本发明的具体实施例中,所述设置当前边界条件可以包括:
步骤11a,由结构初始条件确定壁面温度和位移边界条件,进行高超声速稳态流动计算,同时求解预设的各个守恒方程的耦合解格式,得到初始热流和初始压力。
在本步骤中,可以求解瞬态耦合分析的初始条件。即由结构初始条件确定壁面温度和位移边界条件,进行高超声速稳态流动计算,并将计算果作为瞬态耦合分析的初始条件,即初始热流和初始压力。
在本发明的技术方案中,可以使用多种具体实施方式来实现上述的步骤11a,以下将以其中的一种具体实施方式为例,对本发明的技术方案进行详细的介绍。
例如,在本发明的较佳实施例中,所述各个守恒方程包括:连续守恒方程、动量守恒方程和能量守恒方程。
考虑到高超声速飞行器周围的流场为化学非平衡、热力学平衡的粘性可压缩连续流动,所述连续守恒方程、动量守恒方程和能量守恒方程分别可以表示为:
其中,ρ为密度,v为速度,p为压力,E为比总能量,Ji为扩散通量,为应力张量。
较佳的,在本发明的具体实施例中,对于包含组分混合或反应的流动,所述各个守恒方程中还进一步还需要包括:组分守恒方程。该组分守恒方程可以表示为:
其中,Yi为组分质量分数,Ri为组分i的产物。
另外,在本发明的较佳实施例中,组分的分压可以由组分密度和混合气体的温度求得,而混合气体的压力则可由Dalton定律给出:
pi=ρiRiT (5)
其中,pi为组分i的压力。p为混合气体的总压力。
另外,在本发明的较佳实施例中,在高超声速化学非平衡流动计算中,对于上述的连续守恒方程、动量守恒方程、能量守恒方程和组分守恒方程采用有限体积法同时求解;而空间离散格式可采用AUSM+格式,时间上采用交替方向隐式求解法。
此外,在本发明的技术方案中,由于存在动能耗散和激波,高超声速飞行器周围的空气会达到极高的温度,而高温将使得气体离解甚至电离。因此,在本发明的技术方案中,化学非平衡假设即特征化学反应时间与流动的特征时间相当。可以采用Park 5组分(O,N,NO,O2,N2)17反应化学动力学模型,并考虑第三体效应。反应机制包括三个分解反应和两个交换反应。
其中,第r个反应的一般形式可以表示为:
其中,N表示化学组分数量,v′i,r表示反应物化学计量系数,v″i,r表示生成物化学计量系数,Mi表示组分i。
另外,正反应速率常数kf,r可以用Arrhenius形式表示,而逆反应速率常数kb,r则可通过正反应速率常数求出:
其中,Ar表示指前因子,βr表示温度指数,Er表示活化能,T表示温度,R表示通用气体常数,Kr表示平衡常数。
单一组分的粘性系数、热导率和扩散系数可由气体分子动力论给出:
其中,μi为粘性系数,λi为热导率,Dij为扩散系数。
而混合气体的系数可以由半经验Wilke公式给出,即:
式中,Xi表示组分i的摩尔分数,Mi表示组分i的分子质量。
混合气体的多组分扩散系数则可由以下近似表达式计算:
有限速率化学反应流动的计算需要每种单一组分的热力学属性(比热Cp,i和焓hi)。本文中只考虑流动为化学非平衡、热力学平衡状态的反应气体,每种组分的热力学属性由当地静温计算。高温气体的比热Cp,i和焓hi可以通过温度的多项式曲线拟合函数计算:
Cp,i=R(A1+A2T+A3T2+A4T3+A5T4) (16)
则混合气体的属性可以由单一组分的属性确定:
另外,根据傅里叶定律,热流可以表示为
其中,λ是热导率。
因此,根据上述的公式(19),通过对上述各个公式的求解,即可得到初始热流;根据上述的公式(6),并通过对上述各个公式的求解,即可得到初始压力。
步骤11b,将所述初始热流和初始压力作为当前边界条件。
通过上述的步骤11a~11b,即可通过计算得到当前边界条件。
步骤12,根据结构确定壁面温度和位移边界条件,在流体-固体耦合界面进行数据交换,得到流体区域的当前温度和位移边界条件。
在本步骤中,在获得当前边界条件之后,即可根据结构确定壁面温度和位移边界条件,并通过在流体-固体耦合界面进行数据交换,从而可得到流体区域的当前温度和位移边界条件,即经过数据交换之后的当前温度Tw和位移us。其中,Tw表示在边界处的温度,而us则表示在边界处的位移。
在本发明的技术方案中,可以根据实际应用的需要,使用多种方法来进行数据交换。以下将以其中的一种具体实施方式为例,对本发明的技术方案进行详细的介绍。
在本发明的技术方案中,设置了两个区域:流体区域和固体区域。在上述两个区域中分别进行相关计算时,需要将这两个区域分别进行网格划分,各个数据均被定义在某种类型的网格上。然而,由于这两个区域的网格划分是分别进行的,因此在两个区域的交界处,虽然都是使用网格描述同一个几何体,但两个区域中的单元尺寸和网格节点的位置一般都不相同,这被称之为“不匹配网格”。
因此,在将流体区域的数值作为固体区域的边界条件进行计算时,需要先将流体区域中的边界条件(即网格中的数值)进行数据交换,转换成固体区域中的网格中的数值,然后才能在固体区域中进行相应的计算。同理,在流体区域中进行相应计算之前,也需要先将固体区域中的边界条件通过数据交换转换成流体区域中的网格中的数值。
较佳的,在本发明的具体实施例中,在进行数据交换时,对于位移和温度等场变量可采用最近邻居搜索方法计算;而对于压力和热流通量等载荷则需要满足守恒性,即需要采用守恒插值方法。
较佳的,在本发明的具体实施例中,上述的所述数据交换可以包括:关联和插值。对于每一个目标网格中的点,可以在源网格内搜索离其最近的单元,产生一个节点-单元关系;而邻居搜索则可基于Kd-tree进行。
另外,在本发明的较佳实施例中,可以采用形函数算法完成上述的关联和插值。如果使用的是线性单元,那么可以精确的映射线性函数,而精确的映射二次函数则需要使用二次单元。
步骤13,根据流体区域的当前温度和位移边界条件,在流体区域同时求解预设的各个守恒方程的耦合解格式,在计算一个时间步Δt之后,得到当前的热流和压力。
由于在步骤12中可以得到流体区域的当前温度和位移边界条件,因此在本步骤中即可根据流体区域的当前温度和位移边界条件进行高超声速稳定流动计算,通过同时求解预设的各个守恒方程的耦合解格式的方式,在计算一个时间步Δt之后,得到当前的热流和压力。
较佳的,在本发明的较佳具体实施例中,本步骤13中的“同时求解预设的各个守恒方程的耦合解格式”的具体实现方式与上述步骤11a中的具体实现方式可以是相同或相类似的,因此在此不再赘述。
步骤14,根据当前的热流和压力,在流体-固体耦合界面进行数据交换,得到固体区域的边界条件。
在本步骤中,在获得当前的热流和压力之后,即可通过在流体-固体耦合界面对当前边界条件(即当前的热流和压力)进行数据交换,从而可得到固体区域(即结构求解器)的边界条件,即经过数据交换之后的热流qw和压力Pw。其中,qw表示在边界处的热流,而Pw则表示在边界处的压力。
较佳的,在本发明的较佳具体实施例中,本步骤14中的数据交换与步骤12中的数据交换可以使用相同的具体实现方式,因此在此不再赘述。
步骤15,根据所述固体区域的边界条件,在固体区域通过热力全耦合的方法进行求解,计算一个时间步Δt之后,得到壁面温度Tw和结构位移us。
由于在上述步骤14中获得了固体区域的边界条件(即经过数据交换后的热流qw和压力Pw),因此在本步骤中,即可根据所述固体区域的边界条件,在固体区域(即热流、压力联合载荷下的结构场)中采用热力全耦合的方法进行求解,在计算一个时间步Δt后,即可得到壁面温度Tw和结构位移us。
此外,在本发明的技术方案中,可以根据实际应用的需要,使用多种方法来求解得到壁面温度Tw和结构位移us。以下将以其中的一种具体实施方式为例,对本发明的技术方案进行详细的介绍。
例如,在本发明的技术方案中,基于能量守恒定律和Fourier定律,可以得到结构瞬态热传导方程:
当固体结构受到加热时,其体内温度将发生变化,此时在固体结构内部将会有温差存在。由于热膨胀以及结构的约束作用,固体结构将产生热应力,而热应力又导致了热变形的发生。对于二维结构的响应方程为:
[σ]=[D][B]δe=[S]δe (22)
σx、σy和τxy为固体结构的应力分量;[D]是弹性矩阵,[B]和[S]分别为应变、应力矩阵;δe为位移矩阵,其中,位移矩阵中的元素为结构位移us在x、y、z方向上的分量。因此,根据该位移矩阵δe即可获知相应的结构位移us;然后,根据上述的公式(20)~(22),即可求解得到壁面温度Tw。
在固体区域内,可以采用有限元法进行热-结构分析。考虑应力与温度分布之间的双向耦合关系,进行热力全耦合分析。温度用向后差分格式积分,非线性耦合系统采用牛顿法求解。因此,在本发明的技术方案中,热力与结构问题是同时求解的。
步骤16,判断是否满足预设的停止条件,如果是,则停止整个流程;否则,返回执行步骤12。
在本发明的技术方案中,所述停止条件可以根据实际应用的需要而预先设置。例如,较佳的,在本发明的具体实施例中,所述停止条件为:
当前的计算周期为最后一个计算周期。
在本发明的技术方案中,可以根据实际应用需要预先设置N个计算周期,并在每一个计算周期中都执行上述的步骤12~15。因此,在步骤16中,即可判断当前的计算周期是否为最后一个计算周期,如果是,则停止整个流程;如果不是最后一个计算周期(即不是第N个计算周期),则可以返回执行步骤12,执行下一个计算周期。
另外,在本发明的技术方案中,所述时间步Δt的大小可以根据实际应用需要而预先设置。例如,在本发明的较佳实施例中,所述时间步Δt的值可以是:0.01秒(s)、0.001s或0.0001s等。
为验证本发明中所提出的高超声速流动-传热与结构响应的多场耦合瞬态数值的方法的正确性,以下将以经典的圆柱绕流实验为例,对本发明中的高超声速流动-传热与结构响应的多场耦合瞬态数值的方法进行实验。
在本发明中的上述实验中,所使用的来流马赫数、压力和温度分别为6.47、648.1Pa和241.5K(如图4所示);实验中所使用的不锈钢管的长度、直径和厚度分别为0.6096m、0.0762m和0.0127m。流体区域和固体区域的网格均为结构网格,但在耦合界面处不匹配。对边界层网格进行加密以保证计算结果的网格无关性,并使其具有足够的分辨率。
图4为本发明实施例中的圆柱绕流计算模型。本发明中可以使用如图4所示的流场和圆柱的计算模型。在本发明的技术方案中,计算过程中可以采用SSTk-ω湍流模型进行计算。
图5为本发明实施例中计算得到的驻点温度随时间变化图。图5中提供了不同耦合时间步长时驻点处的温度随时间变化曲线。
图6为本发明实施例中计算得到的驻点热流随时间变化图。图6中所示为驻点热流随时间变化的趋势,其中图6的纵坐标qstag即为边界处的热流qw。如图5和图6所示,随着加热时间的增长,驻点温度也将逐渐升高,而热流则趋势相反。另外,温度和热流都在初始的一小段时间里变化剧烈,而后变化趋于缓和。
在本发明的技术方案中,基于不同耦合时间步长得到了不同的温度分布结果。其中,Δt=0.0001s时得到t=2s时的结果(384.9K)与实验结果(即图5中所示的实验点处的数值,约为388.72K)最为接近,而对于Δt=0.001s和Δt=0.01s,得到了略偏高的温度结果。因此可以得出结论,随着时间步长偏离某个合适的特定值逐渐增大,预测的结构温度相应会逐渐偏高。
图7为本发明实施例中计算得到的表面温度分布。图7给出了Δt=0.0001s时不同时间点的温度分布。其中,图7中的横坐标为曲线长度(Curve Length)。如图7所示,在气动加热开始后的前0.5s内,温度变化比较明显。
图8为本发明实施例中计算得到的结构Mises应力和位移图。如图8所示,在热力联合载荷下,驻点处的位移(displacement)和Mises应力都随时间逐渐增大,但是在当前工况下位移量较小,因此不足以明显影响流场。
通过将上述的实验结果与风洞试验的结果进行对比,验证了本发明中的高超声速流动-传热与结构响应的多场耦合瞬态数值的方法的正确性和有效性。
综上可知,通过使用本发明的高超声速流动-传热与结构响应的多场耦合瞬态数值的方法,可以解决传统方法中热/力载荷环境的预测与热防护结构性能分析互相分离的问题,从而实现高超声速非平衡流动求解器与结构热/力全耦合求解器相耦合的多场耦合计算,使得对高超声速飞行器的气动热力环境和结构热力响应的预测更符合物理实际;而且,由于本发明中考虑了真实气体效应、耦合传热和结构热力全耦合等复杂的效应,从而还可以大大提高高超声速飞行器热/力载荷环境和结构响应预测的计算精度。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (6)
1.一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法,其特征在于,该方法包括:
A、预先建立多物理场耦合模型并设置当前边界条件;
B、根据结构确定壁面温度和位移边界条件,在流体-固体耦合界面进行数据交换,得到流体区域的当前温度和位移边界条件;
C、根据流体区域的当前温度和位移边界条件,在流体区域同时求解预设的各个守恒方程的耦合解格式,在计算一个时间步Δt之后,得到当前的热流和压力;
D、根据当前的热流和压力,在流体-固体耦合界面进行数据交换,得到固体区域的边界条件;
E、根据所述固体区域的边界条件,在固体区域通过热力全耦合的方法进行求解,计算一个时间步Δt之后,得到壁面温度Tw和结构位移us;
F、判断是否满足预设的停止条件,如果是,则停止整个流程;否则,返回执行步骤B。
2.根据权利要求1所述的方法,其特征在于,所述设置当前边界条件包括:
由结构初始条件确定壁面温度和位移边界条件,进行高超声速稳态流动计算,同时求解预设的各个守恒方程的耦合解格式,得到初始热流和初始压力;
将所述初始热流和初始压力作为当前边界条件。
3.根据权利要求1所述的方法,其特征在于,
所述数据交换包括:关联和插值。
4.根据权利要求3所述的方法,其特征在于,在进行数据交换时:
对于位移和温度采用最近邻居搜索方法计算;
对于压力和热流通量载荷则采用守恒插值方法。
5.根据权利要求1所述的方法,其特征在于,所述停止条件为:
当前的计算周期为最后一个计算周期。
6.根据权利要求2所述的方法,其特征在于,所述各个守恒方程包括:
连续守恒方程、动量守恒方程、能量守恒方程和组分守恒方程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510570860.3A CN105095603B (zh) | 2015-09-09 | 2015-09-09 | 一种高超声速流动‑传热与结构响应的多场耦合瞬态数值的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510570860.3A CN105095603B (zh) | 2015-09-09 | 2015-09-09 | 一种高超声速流动‑传热与结构响应的多场耦合瞬态数值的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105095603A CN105095603A (zh) | 2015-11-25 |
CN105095603B true CN105095603B (zh) | 2017-11-14 |
Family
ID=54576028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510570860.3A Active CN105095603B (zh) | 2015-09-09 | 2015-09-09 | 一种高超声速流动‑传热与结构响应的多场耦合瞬态数值的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105095603B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105677994B (zh) * | 2016-01-12 | 2019-02-05 | 北京航空航天大学 | 流体-固体耦合传热的松耦合建模方法 |
CN106742060B (zh) * | 2017-01-03 | 2019-06-18 | 北京临近空间飞行器系统工程研究所 | 一种气动热与材料催化特性耦合效应的地面预示方法 |
CN108304595B (zh) * | 2017-05-04 | 2021-04-02 | 北京空天技术研究所 | 一种用于高超声速飞行器半封闭区域的结构温度分析方法 |
CN107368661B (zh) * | 2017-07-27 | 2020-06-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种高超声速飞行器热气动弹性特性的耦合计算方法 |
CN108255781B (zh) * | 2018-01-04 | 2021-05-11 | 北京环境特性研究所 | 一种高超声速目标表面动态温度建模方法 |
CN108416101A (zh) * | 2018-02-05 | 2018-08-17 | 山东理工大学 | 一种轮毂驱动系统多场耦合建模方法 |
CN111460578A (zh) * | 2020-03-23 | 2020-07-28 | 南京航空航天大学 | 高超声速飞行器头锥热环境高精度流-固耦合计算方法 |
CN112651156B (zh) * | 2020-12-21 | 2022-03-04 | 中国科学院力学研究所 | 一种高速来流条件下激光毁伤效应的多场耦合分析方法 |
CN113076670B (zh) * | 2021-03-26 | 2022-10-18 | 贵州航天电子科技有限公司 | 一种多物理场耦合的相控阵制导微系统协同优化方法 |
CN114626313B (zh) * | 2022-03-04 | 2023-05-16 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种可解析时变热响应的高速气动热cfd求解方法 |
CN115544675B (zh) * | 2022-12-01 | 2023-04-07 | 北京航空航天大学 | 高超声速飞行器防热材料表面催化特性多尺度预测方法 |
CN115828436B (zh) * | 2023-01-09 | 2023-04-28 | 中国空气动力研究与发展中心计算空气动力研究所 | 超声速飞行器进气道的总压损失评估方法及相关组件 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
-
2015
- 2015-09-09 CN CN201510570860.3A patent/CN105095603B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
Non-Patent Citations (5)
Title |
---|
modeling and simulation methodology of multifield coupling;Ping Ma etal.;《Asia Simulation Conference 2012》;20121030;第177-184 * |
热防护材料表面催化特性研究进展;孟松鹤 等;《航空学报》;20140225;第35卷(第2期);第287-302页 * |
高超声速复杂气动问题数值方法研究进展;王江峰 等;《航空学报》;20150125;第36卷(第1期);第159-175页 * |
高超声速飞行器气动热结构耦合问题研究;李昱霖 等;《飞航导弹》;20140515(第5期);第84-92页 * |
高超声速飞行器表面测热技术综述;孟松鹤 等;《航空学报》;20140725;第35卷(第7期);第1759-1774页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105095603A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105095603B (zh) | 一种高超声速流动‑传热与结构响应的多场耦合瞬态数值的方法 | |
Kim et al. | Thermal–hydraulic physical models for a Printed Circuit Heat Exchanger covering He, He–CO2 mixture, and water fluids using experimental data and CFD | |
Bao et al. | Numerical analysis of flowing cracked hydrocarbon fuel inside cooling channels in view of thermal management | |
Gao et al. | Application and theoretical analysis of the flamelet model for supersonic turbulent combustion flows in the scramjet engine | |
Muhammad et al. | OpenFOAM for computational combustion dynamics | |
Bansal et al. | Simulation of hypersonic flow and radiation over a mars reentry vehicle using openfoam | |
Huismann et al. | HyperCODA–extension of flow solver CODA towards hypersonic flows | |
Chen et al. | Direct numerical simulation of chemical non-equilibrium turbulent flow | |
CN106815427A (zh) | 一种火箭推进剂的热物性计算方法和装置 | |
Orhan | Investigation of the effect of turbulence on entropy generation in turbomachinery | |
Kahraman et al. | Coupling of volume of fluid and level set methods in condensing heat transfer simulations | |
Marineau et al. | Turbulent Navier-Stokes simulations of heat transfer with complex wall temperature variations | |
Zhu et al. | Advanced modelling and iterative approach for high-accuracy PCHE design in SCO2 Brayton cycle | |
Li et al. | Study of hybrid NS-DSMC simulation method with chemical non-equilibrium for transitional hypersonic flow | |
Spiegel et al. | A cell-centered finite volume method for chemically reacting flows on hybrid grids | |
Blizard et al. | A Graph-Based Technique for the Automated Control-Oriented Modeling of District Heating Networks | |
Heiser et al. | The Brayton cycle using real air and polytropic component efficiencies | |
Wilson | Assessment and analysis of a novel intake for a ramjet engine | |
Gehre et al. | Revised model coefficients for vibrational relaxation in a nitrogen–oxygen gas mixture | |
Verhoff et al. | Extension of a hybrid particle-continuum method for a mixture of chemical species | |
CN115798625B (zh) | 一种耦合单步化学机理的反应流数值求解方法 | |
Zeng et al. | Computational simulation of reentry flows over hypersonic vehicles using nonlinear coupled constitutive relations | |
Sajjadi et al. | Numerical analysis of turbulent natural convection in a square cavity using large-Eddy simulation in lattice Boltzmann method | |
Grossman et al. | Upwind methods for flows with non-equilibrium chemistry and thermodynamics | |
Furusawa et al. | Numerical Simulation of Supercritical Octane Flows in a Heated Circular Tube With Simple Thermal Cracking Model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |