CN114626313B - 一种可解析时变热响应的高速气动热cfd求解方法 - Google Patents
一种可解析时变热响应的高速气动热cfd求解方法 Download PDFInfo
- Publication number
- CN114626313B CN114626313B CN202210246933.3A CN202210246933A CN114626313B CN 114626313 B CN114626313 B CN 114626313B CN 202210246933 A CN202210246933 A CN 202210246933A CN 114626313 B CN114626313 B CN 114626313B
- Authority
- CN
- China
- Prior art keywords
- heat flow
- time
- temperature
- heat
- surface temperature
- 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
Images
Classifications
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提供一种可解析时变热响应的高速气动热CFD求解方法,所述方法为:根据传热学基础理论,将防热结构表面受热问题假定为半无限大平板非稳态导热问题,从而构建加热表面附近时变温度与热流之间的积分关系,将该积分关系代入基于给定表面温度和给定热流的CFD计算获知的表面温度‑热流线性关联式中,积分获得可解析时变热响应的表面热流结果。本发明不需要空气流动与结构导热耦合的大规模非稳态计算,只需要单独采用CFD计算,即可获得可解析时变热响应条件下的气动加热结果,采用本发明的成本大大降低,有利于气动热环境的快速评估,支撑飞行器热防护***设计和工程应用。
Description
技术领域
本发明涉及数值模拟技术领域,具体而言,涉及一种可解析时变热响应的高速气动热CFD求解方法。
背景技术
CFD被广泛应用于飞行器设计与评估领域,在高速可压缩流动气动热CFD数值模拟中,为确保Navier-Stokes方程的适定性,需要确定热学边界条件。目前通常采用单独给定温度(第一类边界条件)或热流(第二类边界条件),或者给定辐射平衡(广义第三类)边界条件。基于上述边界条件的单独CFD计算难以获得解析材料时变热响应条件下的气动加热结果,需要通过高速流动与结构导热的非稳态耦合计算来获取。但是,相比于非耦合气动加热计算,耦合计算成本大大增加,不利于热防护评估与飞行器工程应用。因此,需要根据材料热响应的物理特征和数学内涵,建立一种用于CFD热学边界的时变材料热响应的解析方案,实现可避免大规模耦合计算的高速气动热计算流体力学(CFD)求解方法。
发明内容
本发明旨在提供一种可解析时变热响应的高速气动热CFD求解方法,以避免通过高速流动与结构导热的非稳态耦合计算来获取气动加热结果的大规模耦合计算导致成本大大增加的问题。
本发明提供的一种可解析时变热响应的高速气动热CFD求解方法,包括如下步骤:
步骤A,根据自由来流条件设定给定温度Tw(x,0)和给定热流qw,g(x),并通过两次CFD计算得到表面温度Tw,g(x)和初始热流qw(x,0);
步骤B,基于步骤A的两次CFD计算结果,构建t时刻的表面温度-热流线性关联式;
步骤C,给定温度Tw(x,0)和初始热流qw(x,0)已知,或基于时刻t表面温度-热流线性关联式,经过i次循环步后ti时刻表面温度Tw(x,ti)和热流qw(x,ti)已知,则将热流qw(x,ti)减去表面辐射散热热流qr得到导入防热结构内部的净热流qc;
步骤D,忽略防热结构内部横向导热和舱内对流传热的影响,将防热结构表面受热问题假定为半无限大平板非稳态导热问题,由此利用净热流qc建立时变温度-热流积分关系;然后在时间步Δti=ti+1-ti内对时变温度-热流积分关系进行有限差分离散,求解ti+1时刻的表面温度;至此,本时间步Δti内的积分计算结束;
步骤E,令i=i+1,若i<n,返回步骤B进行下一时刻积分计算,直到求解至tn时刻的表面温度Tw(x,tn)和热流qw(x,tn)。
进一步的,步骤A中给定温度Tw(x,0)为热响应初始时刻温度,给定热流qw,g(x)为绝热条件;然后分别对设定温度Tw(x,0)和给定热流qw,g(x)进行CFD计算:
根据给定温度Tw(x,0)进行CFD计算得到初始热流qw(x,0);
根据给定热流qw,g(x)进行CFD计算得到表面温度Tw,g(x)。
进一步的,步骤B中所述t时刻的表面温度-热流线性关联式如下:
其中,t时刻表面温度Tw(x,t)和热流qw(x,t)为未知变量,但二者互相关联,且表面温度Tw(x,t)所对应的t时刻未知待求。
进一步的,步骤C中,表面辐射散热热流qr表示为:
其中,∈表示表面发射率;σ表示斯忒藩-玻尔兹曼常数,单位为W/(m2·K4);Tw表示表面温度,单位为K、T∞表示环境温度,单位为K。
进一步的,步骤D中利用净热流qc建立的时变温度-热流积分关系如下:
其中,α表示热扩散系数、κ表示绝热指数、s表示中间变量。
进一步的,步骤D中在时间步Δti=ti+1-ti内对时变温度-热流积分关系进行有限差分离散,求解ti+1时刻的表面温度的公式如下:
其中,Tw(x,ti+1)表示t+1时刻的表面温度。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
本发明不需要空气流动与结构导热耦合的大规模非稳态计算,只需要单独采用CFD计算,即可获得可解析时变热响应条件下的气动加热结果,采用本发明的成本大大降低,有利于气动热环境的快速评估,支撑飞行器热防护***设计和工程应用。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明实施例的可解析时变热响应的高速气动热CFD求解方法的流程图。
图2为本发明实施例的表面温度和热流的线性关系示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
为了避免通过高速流动与结构导热的非稳态耦合计算来获取气动加热结果的大规模耦合计算导致成本大大增加的问题,本发明的构思为:根据传热学基础理论,将防热结构表面受热问题假定为半无限大平板非稳态导热问题,从而构建加热表面附近时变温度与热流之间的积分关系,将该积分关系代入基于给定表面温度和给定热流的CFD计算获知的表面温度-热流线性关联式中,积分获得可解析时变热响应的表面热流结果。由此以通过高速可压缩流动气动热CFD定常解来求解tn时刻的表面温度Tw(x,tn)和热流qw(x,tn)为例,如图1所示,本实施例提出一种可解析时变热响应的高速气动热CFD求解方法,包括如下步骤:
步骤A,如图1所示的第01~04步,根据自由来流条件设定给定温度Tw(x,0)和给定热流qw,g(x),并通过两次CFD计算得到表面温度Tw,g(x)和初始热流qw(x,0);
其中,给定温度Tw(x,0)为热响应初始时刻温度,给定热流qw,g(x)为绝热条件;然后分别对设定温度Tw(x,0)和给定热流qw,g(x)进行CFD计算:
根据给定温度Tw(x,0)进行CFD计算得到初始热流qw(x,0);
根据给定热流qw,g(x)进行CFD计算得到表面温度Tw,g(x)。
步骤B,如图1所示的第05~06步,基于步骤A的两次CFD计算结果,构建t时刻的表面温度-热流线性关联式:
其中,t时刻表面温度Tw(x,t)和热流qw(x,t)为未知变量,但二者互相关联,满足如图2所示的线性关系,且表面温度Tw(x,t)所对应的t时刻未知待求。
步骤C,如图1所示的第07~09步,给定温度Tw(x,0)和初始热流qw(x,0)已知,或基于时刻t表面温度-热流线性关联式,经过i次循环步后ti时刻表面温度Tw(x,ti)和热流qw(x,ti)已知,则将热流qw(x,ti)减去表面辐射散热热流qr得到导入防热结构内部的净热流qc。
所述表面辐射散热热流表示为:
其中,∈表示表面发射率;σ表示斯忒藩-玻尔兹曼常数,单位为W/(m2·K4);Tw表示表面温度,单位为K;T∞表示环境温度,单位为K。
则在表面上满足如下能量平衡关系:
qc(x,ti)=qw(x,ti)-qr(x,ti)
其中,qc(x,ti)即i次循环步后ti时刻的净热流;qr(x,ti)即i次循环步后ti时刻的表面辐射散热热流。
步骤D,如图1所示的第10步,忽略防热结构内部横向导热和舱内对流传热的影响,将防热结构表面受热问题假定为半无限大平板非稳态导热问题,由此利用净热流qc建立时变温度-热流积分关系:
其中,α表示热扩散系数、κ表示绝热指数、s表示中间变量。
然后在时间步Δti=ti+1-ti内对时变温度-热流积分关系进行有限差分离散,求解ti+1时刻的表面温度:
其中,Tw(x,ti+1)表示t+1时刻的表面温度;
至此,本时间步Δti内的积分计算结束。
步骤E,如图1所示的第10~11步,令i=i+1,若i<n,返回步骤B进行下一时刻积分计算,直到求解至tn时刻的表面温度Tw(x,tn)和热流qw(x,tn),至此求解结束。
通过上述求解过程可见,本发明不需要空气流动与结构导热耦合的大规模非稳态计算,只需要单独采用CFD计算,即可获得可解析时变热响应条件下的气动加热结果,采用本发明的成本大大降低,有利于气动热环境的快速评估,支撑飞行器热防护***设计和工程应用。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种可解析时变热响应的高速气动热CFD求解方法,其特征在于,包括如下步骤:
步骤A,根据自由来流条件设定给定温度Tw(x,0)和给定热流qw,g(x),并通过两次CFD计算得到表面温度Tw,g(x)和初始热流qw(x,0);
步骤B,基于步骤A的两次CFD计算结果,构建t时刻的表面温度-热流线性关联式;
步骤C,给定温度Tw(x,0)和初始热流qw(x,0)已知,或基于时刻t表面温度-热流线性关联式,经过i次循环步后ti时刻表面温度Tw(x,ti)和热流qw(x,ti)已知,则将热流qw(x,ti)减去表面辐射散热热流qr得到导入防热结构内部的净热流qc;
步骤D,忽略防热结构内部横向导热和舱内对流传热的影响,将防热结构表面受热问题假定为半无限大平板非稳态导热问题,由此利用净热流qc建立时变温度-热流积分关系;然后在时间步Δti=ti+1-ti内对时变温度-热流积分关系进行有限差分离散,求解ti+1时刻的表面温度;至此,本时间步Δti内的积分计算结束;
步骤E,令i=i+1,若i<n,返回步骤B进行下一时刻积分计算,直到求解至tn时刻的表面温度Tw(x,tn)和热流qw(x,tn);
步骤D中利用净热流qc建立的时变温度-热流积分关系如下:
其中,α表示热扩散系数、κ表示绝热指数、s表示中间变量;
步骤D中在时间步Δti=ti+1-ti内对时变温度-热流积分关系进行有限差分离散,求解ti+1时刻的表面温度的公式如下:
其中,Tw(x,ti+1)表示t+1时刻的表面温度。
2.根据权利要求1所述的可解析时变热响应的高速气动热CFD求解方法,其特征在于,步骤A中给定温度Tw(x,0)为热响应初始时刻温度,给定热流qw,g(x)为绝热条件;然后分别对设定温度Tw(x,0)和给定热流qw,g(x)进行CFD计算:
根据给定温度Tw(x,0)进行CFD计算得到初始热流qw(x,0);
根据给定热流qw,g(x)进行CFD计算得到表面温度Tw,g(x)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210246933.3A CN114626313B (zh) | 2022-03-04 | 2022-03-04 | 一种可解析时变热响应的高速气动热cfd求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210246933.3A CN114626313B (zh) | 2022-03-04 | 2022-03-04 | 一种可解析时变热响应的高速气动热cfd求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114626313A CN114626313A (zh) | 2022-06-14 |
CN114626313B true CN114626313B (zh) | 2023-05-16 |
Family
ID=81901645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210246933.3A Active CN114626313B (zh) | 2022-03-04 | 2022-03-04 | 一种可解析时变热响应的高速气动热cfd求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114626313B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095603A (zh) * | 2015-09-09 | 2015-11-25 | 哈尔滨工业大学 | 一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法 |
CN109408926A (zh) * | 2018-10-12 | 2019-03-01 | 大连理工大学 | 求解复杂结构多维瞬态非线性热传导反问题的方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598602B (zh) * | 2009-07-10 | 2011-05-18 | 北京航空航天大学 | 导弹高温气动热试验弹体前表面温度的非接触式红外测量装置 |
CN101916314B (zh) * | 2010-08-16 | 2012-12-05 | 北京理工大学 | 高速飞行器升力面气动热结构多学科优化设计平台 |
GB2591972B8 (en) * | 2014-12-11 | 2022-03-16 | Fulton Group N A Inc | Ribbed tubeless heat exchanger for fluid heating systems including a rib component and methods of manufacture thereof |
CN109804217B (zh) * | 2016-10-14 | 2024-05-28 | 达纳加拿大公司 | 具有空气动力学特征以改善性能的热交换器 |
CN107368661B (zh) * | 2017-07-27 | 2020-06-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种高超声速飞行器热气动弹性特性的耦合计算方法 |
CN107808065B (zh) * | 2017-11-23 | 2019-12-31 | 南京航空航天大学 | 三维复杂外形高速飞行器流-固-热快速计算方法 |
-
2022
- 2022-03-04 CN CN202210246933.3A patent/CN114626313B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095603A (zh) * | 2015-09-09 | 2015-11-25 | 哈尔滨工业大学 | 一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法 |
CN109408926A (zh) * | 2018-10-12 | 2019-03-01 | 大连理工大学 | 求解复杂结构多维瞬态非线性热传导反问题的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114626313A (zh) | 2022-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Marco et al. | A note on limiting laminar Nusselt number in ducts with constant temperature gradient by analogy to thin-plate theory | |
Aziz et al. | Convective‐radiative fins with simultaneous variation of thermal conductivity, heat transfer coefficient, and surface emissivity with temperature | |
Chen et al. | Inverse problem in determining convection heat transfer coefficient of an annular fin | |
CN111859532B (zh) | 考虑高超声速化学非平衡效应的改进热壁修正方法 | |
Ahmad et al. | Stochastic numerical computing with Levenberg-Marquardt backpropagation for performance analysis of heat sink of functionally graded material of the porous fin | |
CN108053330B (zh) | 一种污染物无组织排放量的计算方法及*** | |
Weng et al. | The complex variable reproducing kernel particle method for two-dimensional inverse heat conduction problems | |
Wang et al. | Forced convection gaseous slip flow in a porous circular microtube: An exact solution | |
Mills-Brown et al. | The development of a high temperature tensile testing rig for composite laminates | |
CN114626313B (zh) | 一种可解析时变热响应的高速气动热cfd求解方法 | |
Thunman et al. | Separation of drying and devolatilization during conversion of solid fuels | |
Wu et al. | Insulation performance of heat-resistant material for high-speed aircraft under thermal environments | |
Marek | Basic research of thermal transfer simulations | |
Huang et al. | Fluid flow and heat transfer across a curved hollow fiber membrane tube bank (CHFMTB): Effects of the tube deformations | |
Lee et al. | Analytic solutions for heat conduction in functionally graded circular hollow cylinders with time‐dependent boundary conditions | |
Zhu et al. | High-order compact difference scheme of 1D nonlinear degenerate convection–reaction–diffusion equation with adaptive algorithm | |
Su et al. | Improved lumped models for transient heat conduction in a slab with temperature-dependent thermal conductivity | |
CN111898263B (zh) | 工程含热源结构热传导情况求解方法以及装置 | |
Hasan | Heat transfer analysis in thermal entrance region under turbulent flow conditions | |
Conley et al. | An assessment of the accuracy of numerical solutions to the Graetz problem | |
Juncu | Unsteady heat transfer from an elliptic cylinder | |
Wessies et al. | Cooling of heated solid cylinder supported on bedded and embedded substrates by impinging air jet | |
Yanou et al. | Estimation of thermal conductivity for model with radiative heat transfer by extended Kalman filter | |
KR20060091633A (ko) | 휠름 블로잉 공정의 동적 모사와 과도해 계산 도출 방법 | |
Subramanian | " Unsteady Heat Transfer: Lumped Thermal Capacity Model |
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 |