Next Article in Journal
Aircraft Trajectory Clustering in Terminal Airspace Based on Deep Autoencoder and Gaussian Mixture Model
Next Article in Special Issue
Aerodynamic Characteristics of Re-Entry Capsules with Hyperbolic Contours
Previous Article in Journal
Prediction of Aircraft Noise Impact with Application to Hong Kong International Airport
Previous Article in Special Issue
Numerical Study on Low-Temperature Region as Heat Sink and Its Heat Dissipation Capacity for Hypersonic Vehicle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Rapid Data-Driven Prediction Method of Coupled Fluid–Thermal–Structure for Hypersonic Vehicles

1
School of Aeronautic Science and Engineering, Beihang University (BUAA), Beijing 100191, China
2
Science and Technology on Space Physics Laboratory, Beijing 100076, China
*
Author to whom correspondence should be addressed.
Aerospace 2021, 8(9), 265; https://doi.org/10.3390/aerospace8090265
Submission received: 22 August 2021 / Revised: 13 September 2021 / Accepted: 14 September 2021 / Published: 16 September 2021
(This article belongs to the Special Issue Hypersonics: Emerging Research)

Abstract

:
This work demonstrates the use of Latin Hypercube Sampling and Proper Orthogonal Decomposition in combination with a Radial Basis Function model to perform on vehicle prediction coupled fluid–thermal–structure. We regarded the Mach number, flight altitude and angle of attack as input parameters and established a rapid prediction model. The basic process of numerical simulation of the hypersonic vehicle coupled fluid–thermal–structure was studied to obtain the database of pressure coefficient, heat flux, structural temperature and structural stress as the sample data to train this prediction method. The prediction error was analyzed. The prediction results showed that the data-driven method proposed in this paper based on proper orthogonal decomposition and radial basis function could be used for predicting vehicle coupled fluid–thermal–structure with good efficiency.

1. Introduction

A hypersonic vehicle is a challenging research area. This is due to the fact that researchers need to deal with a large flight envelope, rapidly changing aerodynamic coefficients, as well as coupling between the structure and the aerodynamics of the vehicle [1]. Design of structures [2] and thermal protection systems [3] for hypersonic vehicles depend on accurate predictions of the aerothermal loads, structural temperatures, and structural stresses [4,5,6]. Traditionally, a rigid isothermal body is assumed to analyze the surface pressures and heating rates [7]. These aerodynamic heating rates are used to analyze the structural temperature. The temperature and aerodynamic pressures are used to analyze the structural deformations and stresses [8,9]. The traditional independent approaches are inefficient because several iterations are required between the different analyses [10]. Aerodynamic heat is generated during the process of flight, which is transferred to the structure [11]. Thermal deformation of the structure is generated due to aerodynamic heat, which in turn affects the flow field. The problem should be simulated by the method of coupling fluid–thermal–structure. Pramote Dechaumphai et al. [12] carried out a study of a numerical simulation on leading edges coupled fluid–thermal–structure, and compared this with the experiment result. Many researchers such as Hirschel [2], Anderson [3], and so on had gone into studies on the aerodynamic heat of hypersonic vehicles. Apart from coupling, there is a large flight envelope. It is impossible to simulate the results of all flight conditions in the stage of engineering design [13,14], which will cost a lot of time. Therefore, it is necessary to develop a prediction method that is equivalent to CFD accuracy on vehicle coupled fluid–thermal–structure.
The Proper Orthogonal Decomposition (POD), also known as Principle Components Analysis (PCA), has been widely used for a broad range of applications. POD analysis yields a set of empirical modes, which describes the dominant behavior or dynamics of a given problem. This technique can be used for a variety of applications, including derivation of reduced-order dynamical models, steady analysis and design of inviscid airfoils, image processing, and pattern recognition. Sirovich [15] introduced the method of snapshots as a way to efficiently determine the POD modes for large problems. In particular, the method of snapshots has been widely applied to computational fluid-dynamic (CFD) formulations to obtain reduced-order models for unsteady aerodynamic applications. A set of instantaneous flow solutions, or “snapshots” is obtained from a simulation of the CFD method. The POD process then computes a set of modes from these snapshots, which is optimal in the sense that, for any given basis size, the error between the original and reconstructed data is minimized. Reduced-order models can be derived by projecting the CFD model onto the reduced space spanned by the POD modes. Duan et al. [16] solved the inverse design problem of two-dimensional airfoil based on this method. Adding disturbance to the reference airfoil and performing CFD calculations to obtain sampling solutions. The airfoil inverse design can be performed based on the method of fitting the POD coefficients of the data, the airfoil shape corresponding to target pressure distribution can be achieved. Bui Thanh Tan [17] proposed to interpolate the POD coefficient under the condition of the known incoming flow state (such as Mach number, angle of attack) corresponding to the pressure coefficient distribution. This method can quickly obtain the approximate pressure coefficient distribution of the unknown flow field and greatly improves the calculation efficiency. The research object of these articles is very nearly flow field, and there are very few studies on the prediction of physical quantities coupled fluid–thermal–structure.
In this paper, we demonstrate the use of Latin Hypercube Sampling (LHS) and Proper Orthogonal Decomposition (POD) in combination with a Radial Basis Function (RBF) model to perform a prediction of the physical quantities coupled fluid–thermal–structure.

2. Rapid Data-Driven Prediction Method

In this section, we briefly review the theory of the POD and RBF Interpolation, followed by the rapid data-driven prediction method of a vehicle coupled fluid–thermal–structure, which are employed throughout this paper.

2.1. Proper Orthogonal Decomposition

The basic POD procedure is summarized briefly here. POD is based on the factorization of a spatial–temporal data matrix obtained from simulations or experiments [18]. The method is based on singular value decomposition (SVD) and attempts to find a linear lower-dimensional subspace than the original data. Through SVD, it is possible to obtain computed POD modes which are ranked in an energetic sense according to their residual variance [19].
Let Y = [ y 1 , y 2 , , y n ] p × n be the response matrix, which is obtained through experiments or simulations. In general, for numerical simulations, p represents the number of model nodes, n represents the number of snapshots. The purpose of POD is to find the optimal basis vector. The optimal POD basis vectors u are chosen to maximize the cost:
{ max   1 n i = 1 n | ( Y ,   u ) | 2 ( u , u ) = 1
where ( Y , u ) is the inner product of the basis vector u with the field Y . Equation (1) can be solved by the Lagrange multiplier method:
J ( u ) = i = 1 n ( ( Y i ,   u ) 2 λ ( u 2 1 ) )
where λ is the Lagrange multiplier, then take the derivative of u :
d d u J ( u ) = 2 Y Y T u 2 λ u  
If we let C = Y Y T , C is a symmetric positive definite matrix. In general, for numerical simulations that generate big data problems, the snapshot POD method is more appropriate since it builds the covariance matrix with the Y transposed data matrix. This leads to a covariance matrix that has a reduced size, but with the same eigenvalues and eigenvectors. Hence, the basis vectors are found by the eigen-decomposition problem:
C u = λ u
where the eigenvalues λ are positive and the eigenvectors u form an orthonormal basis.
Consider energy, the larger the eigenvalue, the more energy it represents, and in terms of containing information, the larger the eigenvalue, the more information it has. The relative “energy” captured by the i th basis vector is given by:
I ( r ) = i = 1 r λ i i = 1 n λ i  
where I ( r ) is energy ratio, the approximate prediction of the field Y is then given by a linear combination of the eigenfunctions.
The basic POD procedure (just outlined) considers time-varying flows by taking a series of flow solutions at different instants in time [20]. In this paper, the procedure is applied in parameter space, that is, obtaining snapshots while allowing a parameter to vary. The parameter of interest could, for example, be the Mach number, angle of attack, and altitude.

2.2. Radial Basis Function Interpolation

Let f ^ ( x ) be the original function to be modeled, and f ( x i ) be the values at N discrete points x i ,   i = 1 ,   2 , ,   n , where x i is the vector of inputs at the i th sample point. The set of data points X = { x 1 ,   x 2 ,   ,   x n } is confined to a domain Ω in n -dimensional space, which is normalized to the unit hypercube [ 0 ,   1 ] n for interpolation [21].
The RBF model is a linear combination of basis functions, whose argument is the Euclidean distance between the point x at which the interpolation is made and all the other points in the known dataset [22]. In other words, the interpolation at an untried site is a sum of contributions from all the known function values, the influence of which is controlled by a basis function that depends on the distance they are from the new site [23]. It provides a comprehensive and flexible interpolation method. Due to its good approximation, robustness and ease of implementation, the RBF model is used in many fields such as statistics, engineering optimization, and so on [24].
If φ ( x ) is the chosen basis function and · is used to denote the Euclidean norm, the interpolation model has the form:
f ^ ( x ) = i = 1 n ω i · φ ( x x i )
where ω i ,   i = 1 ,   2 , , n is the weight coefficient. The coefficients are found by requiring exact recovery of the original data f ^ ( x ) = f ( x ) for all points in the training dataset X .
When there are no repeated sample points, the radial function φ ( x ) is positive The interpolation coefficients λ 1 ,   λ 1 , , λ n can be found by solving the linear system problem:
[ φ ( x 1 x 1 ) φ ( x 1 x 2 ) φ ( x 1 x n ) φ ( x 2 x 1 ) φ ( x 2 x 2 ) φ ( x 2 x n ) φ ( x n x 1 ) φ ( x n x 2 ) φ ( x n x n ) ] A [ ω 1 ω 2 ω n ]   = [ f ( x 1 ) f ( x 2 ) f ( x n ) ]
where A is the interpolation matrix. The obtained RBF interpolant f ^ ( x ) can be used to approximate the given function f ( x ) .

2.3. Rapid Prediction Method Based on Data-Driven

In our work, we made predictions on the hypersonic vehicle coupled fluid–thermal–structure, which concerned heat flux of flow field, pressure coefficient of the flow field, structural temperature and structural stress. POD provides a method of processing the dataset [25,26,27,28]. Thereby the specific process of the prediction method [29,30,31,32] is introduced in detail next, and the corresponding flow chart is shown in Figure 1.
Firstly, the state variables and their range is selected according to the flight status. Latin Hypercube Sampling is used to obtain sample points X in the design space as shown in Figure 2. The state variables and their range is shown in Table 1.
Secondly, numerical simulation is adopted to obtain the response matrix of all sample points Y . The response matrix Y can be expressed as the superposition of the average value Y ¯ and the pulsation value Y , that is Y = Y ¯ + Y . The pulsation value Y is used to construct matrix M :
M = [ Y 1 ( 1 ) Y 2 ( 1 ) Y 1 ( 2 ) Y 2 ( 2 ) Y 1 ( n ) Y 2 ( n )   Y p ( 1 ) Y p ( 2 ) Y p ( n ) ]      
where n is the number of sampling points; p is the number of model nodes; Y is pulsation response value of model nodes, which could be heat flux, pressure coefficient, structural temperature and structural stress in this paper.
Thirdly, pulsation response matrix M is used to construct the correlation matrix C . Then, we decompose the eigenvalue of C , C = M T M , C A = λ A , u 0 = 1 / λ M A , where λ and A respectively correspond to eigenvalue and eigenvector, u 0 is the orthogonal basis vector.
Fourthly, eigenvalues are sorted from large to small, and the larger the eigenvalue, the more energy it represents. POD basis vector u , which is m dimensional ( m n ) , is selected according to λ .
Fifthly, the POD basis vector is obtained in the fourth step. From the approximate relationship Y ( i ) j = 1 m a j ( i ) u ( j ) , the coefficients a of the POD basis vector under all sample points can be calculated by the least square method. Due to the relationship between sample points X and coefficients a , RBF is used to establish their approximate fitting relationship [33,34,35,36,37,38].
Sixthly, the argument of RBF is the Euclidean distance between the point x at which the interpolation is made and all the other points in the known dataset. Given a test sample point X ( e ) in the design space, the corresponding coefficient a ( e ) can be obtained, followed by response value of physical quantities can be obtained according to Y ( e ) = Y ¯ ( i ) + j = 1 m a j ( e ) u ( j ) .

3. Numerical Simulation

The numerical simulation coupled fluid–thermal–structure was studied to obtain the response value of sampling points to establish a database set. The simulation should be verified to ensure that the database we build can accurately contain the information of the flow field and structure.
The equations for unsteady compressible flow are described by the conservation of mass, momentum, and energy equations [12]. These equations can be written in the conservation form as:
t { Q F } + x { E F I + E F V } + y { G F I + G F V } = 0  
where { Q F } is the vector of the conservation variables; { E F I } and { G F I } are vectors of the inviscid flux components in the x and y directions; { E F V } and { G F V } are vectors of the viscous flux components in the x and y directions.
For transient heat conduction without an internal heat source, the thermal response of the structure is described by the energy equation that can be written in conservation form as:
t ( Q T ) + x ( E T ) + y ( F T ) = 0
where Q T is the conservation variable, E T and F T are the component of flux.
The structural response is described by the quasistatic equilibrium equations that can be written in conservation form as:
t ( Q s ) + x ( E s ) + y ( F s ) = 0
where Q s is the displacement vector; E s and F s are vectors of the stress components.
An experiment was performed in the NASA Langley 8-ft High-Temperature Tunnel in 1987. The NASA Langley Research Center 8-ft High-Temperature Tunnel (8-ft HTT) is a hypersonic blowdown tunnel in which a high energy level for simulating hypersonic flight is obtained by burning methane and air in a high-pressure combustor. A 3-inch diameter and 0.5-inch thick, stainless-steel cylinder was mounted on the panel holder and subjected to a uniform high-enthalpy Mach 6.47 flow. The solution to Mach 6.47 flow over a cylinder is used to demonstrate the integrated fluid–thermal–structure analysis approach. The schematic diagram of this 8-ft HTT can be found in Figures 8 and 9 of [10]. And more details of the experimental configurations, the tunnel flow conditions, and the experimental results can be found in [10].
Because of symmetry, the model is reduced to two dimensions. Only one-half of the incoming flow domain and cylinder is modeled. In this paper, the finite-element model representing the flow domain and the cylinder is shown in Figure 3. The number of grids of fluid is 150 × 400. The number of grids of structure is 100 × 13. The grid near the wall of fluid is encrypted, and the minimum height is 10 6 m. In order to better capture the characteristics of the flow field, the grid at the location of the shock wave is encrypted.
The hypersonic aerodynamic heat is calculated by FLUENT and the structural heat transfer is calculated by ABAQUS. The pressure far field boundary condition is used in the flow domain and the S S T   k ω turbulence model is used in the whole flow computational domain. A no-slip condition is specified at the surface of the wall. The analysis step of coupling temperature-displacement is used in the cylinder and the initial temperature of the cylinder is 294.4 K. The code is used to realize data exchange in the process of simulation through the surface of the wall between ABAQUS and FLUENT. The incoming condition is shown in Table A1. The material property is shown in Table A2.
The result of the numerical simulation is shown in Figure 4 and Table 2. Density contours are compared in Figure 4 with an experiment result [10]. The comparison of the shock shape and position indicates the global flow field is reasonably well obtained. We also compared the result of the maximum temperature of structure, maximum heat flux and the maximum pressure of flow field between the numerical simulation in this paper and experiment, which is shown in Table 2. We could conclude that the numerical simulation used in this paper could obtain the information of flow field and structure for establishing a database. For the experimental results and how they were obtained please refer to [10].

4. Result of Prediction

Five test conditions were selected according to the Latin Hypercube Sampling method. We take the arithmetic average of these predicted results. The error analysis is performed on an arithmetic average. The test conditions are shown in Table 3.
Heat flux and pressure coefficient of the wall of the fluid domain are usually of interest to researchers. By comparing the actual value with predicted value, the accuracy of the prediction method could be judged. The wall consists of 355 nodes, the actual pressure coefficient and predicted pressure coefficient of these nodes are shown in Figure 5. The actual heat flux and predicted heat flux of these nodes are shown in Figure 6. It could be seen that the predicted value is almost consistent with the actual value regardless of the pressure coefficient or the heat flux.
The contours of actual and predicted results, including temperature and stress, are shown in Figure 7 and Figure 8. Through comparing actual and predicted results, it could be seen that the actual result is similar to the predicted result in terms of value and distribution.
In order to make the results more intuitive, the relative error of the structural temperature and stress are shown in Figure 9.
It could be seen that the maximum relative error of structural temperature is only 2%, which is practically located at the front of the structure. The maximum relative error of structural stress is only 3.8%, which appears in very few nodes. However, the test sets of data were simply to verify the prediction effect of the method. There might have been a few errors between the different sets of data that could be reduced by obtaining more sets of verification data. The prediction error can be reduced by increasing the sample data.
The comparison of efficiency between the predicted method and numerical simulation is shown in Table 4. The computer is configured as AMD Ryzen 7 4800 H, 2.9 GHz, 16.0 GB RAM. The calculation with the software took nearly 12 h. It is difficult to obtain the information of flow field and structure quickly. However, it only took 0.14 s to complete the calculation by using the prediction method proposed in this paper.

5. Conclusions

In this paper, the prediction accuracy of the data-driven method on vehicle coupled fluid–thermal–structure was investigated for pressure coefficient, heat flux, structural temperature and structural stress by use of the POD and RBF models. According to the analysis of Section 4, it could be seen that the data-driven method proposed in this paper could be used for predicting the physical quantity of a hypersonic vehicle coupled fluid–thermal–structure with good efficiency. Compared with the numerical simulation method, the data-driven prediction method had the advantages of low cost and high speed. If sufficient data can be provided, the calculation result will be consistent with the actual values.
In future work, in order to improve the accuracy of the prediction method, we will study the relationship between the accuracy of the prediction result and the size of the database. The combination of proper orthogonal decomposition with other machine learning models will be studied, and we hope to improve the generalization ability of the prediction methods.

Author Contributions

Conceptualization, J.L., M.W. and S.L.; methodology, J.L. and M.W.; formal analysis, J.L.; writing—original draft preparation, J.L.; writing—review and editing, J.L. and M.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The incoming flow conditions and material properties of the stainless-steel cylinder are given in the following table:
Table A1. The incoming flow conditions.
Table A1. The incoming flow conditions.
P ( P a ) T ( K ) M a R e
648.1 241.5 6.47 1.31 × 10 6
Table A2. The material property of the stainless-steel cylinder.
Table A2. The material property of the stainless-steel cylinder.
Elastic Modulus
( G P a )
Poisson’s RatioDensity
( K g / m 3 )
Coefficient of Linear Expansion
× 10 6 ( K 1 )
Thermal Conductivity
( W / ( m · K ) )
Specific Heat Capacity
( J / ( k g · K ) )
2060.3803017.516.27502.48

References

  1. Zel’dovich, Y.B.; Raizer, Y.P. Elements of Gasdynamics and the Classical Theory of Shock Waves. In Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena; Hayes, W.D., Probstein, R.F., Eds.; Dover Publications: Mineola, NY, USA, 1967; pp. 1–104. [Google Scholar]
  2. Hirschel, E.H. The Flight Environment. In Basics of Aerothermodynamics; Springer: Berlin, Germany, 2005; pp. 15–26. [Google Scholar]
  3. Anderson, J.D. Hypersonic Shock and Expansion-Wave Relations. In Hypersonic and High-Temperature Gas Dynamic; AIAA Press: Palo Alto, CA, USA, 2019; pp. 39–53. [Google Scholar]
  4. Shang, J.S.; Surzhikov, S.T. Radiative Energy Transfer. In Plasma Dynamics for Aerospace Engineering; Cambridge University Press: New York, NY, USA, 2018; pp. 304–336. [Google Scholar]
  5. Gui, Y.; Yuan, X. Numerical simulation on the coupling phenomena of aerodynamic heating with thermal response in the region of the leading edge. J. Eng. Thermophys. 2002, 23, 733–735. [Google Scholar]
  6. Chen, F.; Liu, H.; Zhang, S. Time-adaptive loosely coupled analysis on fluid–thermal–structural behaviors of hypersonic wing structures under sustained aeroheating. Aerosp. Sci. Technol. 2018, 78, 620–636. [Google Scholar] [CrossRef]
  7. Geng, X.; Liu, L.; Gui, Y.; Tang, W.; Wang, A. Studying the Temperature Field of Hypersonic Vehicle Structure with Aero-Thermo-Elasticity Deformation. Int. Sch. Sci. Res. Innov. 2015, 9, 1434–1437. [Google Scholar]
  8. Zhao, X.; Sun, Z.; Tang, L.; Zheng, G. Coupled flow-thermal-structural analysis of hypersonic aerodynamically heated cylindrical Leading Edge. Eng. Appl. Comput. Fluid Mech. 2011, 5, 170–179. [Google Scholar] [CrossRef] [Green Version]
  9. Korkegi, R.H. Survey of Viscous Interactions Associated with High Mach Number Flight. AIAA J. 1971, 9, 771–784. [Google Scholar] [CrossRef]
  10. Wieting, A.R. Experimental Study of Shock Wave Interference Heating on a Cylindrical Leading Edge. Ph.D. Thesis, Old Dominion University, Norfolk, VA, USA, 1987. [Google Scholar]
  11. Thornton, E.A.; Dechaumphai, P. Coupled Flow, Thermal and Structural Analysis of Aerodynamically heated Panels. J. Aircr. 1988, 25, 1052–1059. [Google Scholar] [CrossRef]
  12. Dechaumphai, P.; Thornton, E.A.; Wieting, A.R. Flow-Thermal-Structural Study of Aerodynamic Heated Leading Edges. J. Aircr. 1989, 26, 201–209. [Google Scholar]
  13. Li, J.; Wang, J.; Yang, L.; Shu, C. A hybrid lattice Boltzmann flux solver for integrated hypersonic fluid-thermal-structural analysis. Chin. J. Aeronaut. 2020, 33, 2295–2312. [Google Scholar] [CrossRef]
  14. Liu, H.; Wang, Z. Fluid–thermal–structural coupling investigations of opposing jet in hypersonic flows. Int. Commun. Heat Mass Transf. 2021, 120, 105017. [Google Scholar]
  15. Sirovich, L. Turbulence and the dynamics of coherent structures. Part 1: Coherent structures. Quart. Appl. Math. 1987, 45, 561–571. [Google Scholar] [CrossRef] [Green Version]
  16. Duan, Y.; Cai, J.; Li, Y. Gappy proper orthogonal decomposition-based two-step optimization for airfoil design. AIAA J. 2012, 50, 968–971. [Google Scholar] [CrossRef]
  17. Tan, B.; Damodaran, M.; Willcox, K. Aerodynamic data reconstruction and inverse design using proper orthogonal decomposition. AIAA J. 2004, 42, 1505–1516. [Google Scholar]
  18. Marc, O.; Markus, B.; Kilian, O.; Sonja, S. Statistical characterization of horizontal slug flow using snapshot proper orthogonal decomposition. Int. J. Multiph. Flow 2021, 134, 103453. [Google Scholar]
  19. Huang, D.; Jan, N.; Christian, W.; Peter, W. A machine learning based plasticity model using proper orthogonal decomposition. Comput. Methods Appl. Mech. Eng. 2020, 365, 113008. [Google Scholar] [CrossRef] [Green Version]
  20. Legresley, P.; Alonso, J. Airfoil design optimization using reduced order models based on proper orthogonal decomposition. Proc. AIAA Space 2000, 40, 1954–1960. [Google Scholar]
  21. Esmaeilbeigi, M.; Garmanjani, G. Gaussian radial basis function interpolant for the different data sites and basis centers. Calcolo 2017, 54, 155–166. [Google Scholar] [CrossRef]
  22. Wang, X.D.; Ding, Y.H.; Shao, H.H. The improved radial basis function neural network and its application. Artif. Life Robot. 1998, 2, 8–11. [Google Scholar] [CrossRef]
  23. Tyshchenko, O. A reservoir radial-basis function neural network in prediction tasks. Autom. Control Comput. Sci. 2016, 50, 65–71. [Google Scholar] [CrossRef]
  24. Cheng, Q.S.; Wu, L.W.; Wang, S.Z. Fusion prediction based on the attribute clustering network and the radial basis function. Chin. Sci. Bull. 2001, 46, 789–792. [Google Scholar] [CrossRef]
  25. Liu, Z.; Ning, F.; Zhai, Q.; Ding, H.; Wei, J. Study on the flow characteristics in the supersonic morphing cavities using direct numerical simulation and proper orthogonal decomposition. Wave Motion 2021, 104, 102751. [Google Scholar] [CrossRef]
  26. Di, G.; Dakun, S.; Ruize, X.; Daniel, B.; Sun, X.F.; Ni, S.L.; Du, J.; Zhao, D. Experimental investigation on axial compressor stall phenomena using aeroacoustics measurements via empirical mode and proper orthogonal decomposition methods. Aerosp. Sci. Technol. 2021, 112, 106655. [Google Scholar]
  27. Lukas, K.; Cyrille, V.; Christian, B. Using a Proper Orthogonal Decomposition representation of the aerodynamic forces for stochastic buffeting prediction. J. Fluids Struct. 2020, 99, 103178. [Google Scholar]
  28. Chen, X.; Liu, L.; Yue, Z. Reduced order aerothermodynamic modeling research for hypersonic vehicles based on proper orthogonal decomposition and surrogate method. Acta Aeronaut. Astronaut. Sin. 2015, 36, 462–472. [Google Scholar]
  29. Kou, J.; Zhang, W. Dynamic mode decomposition and its applications in fluid dynamics. Acta Aerodyn. Sin. 2018, 36, 163–179. [Google Scholar]
  30. Wang, Y.; Han, R.; Liu, Z. Progress of deep learning modeling technology for fluid mechanics. Acta Aeronaut. Astronaut. Sin. 2021, 42, 524779. [Google Scholar]
  31. Cui, E. MEMS and Intelligent Fluid Mechanics. Acta Aerodyn. Sin. 2000, 18, 52–59. [Google Scholar]
  32. Karami, S.; Soria, J. Analysis of Coherent Structures in an Under-Expanded Supersonic Impinging Jet Using Spectral Proper Orthogonal Decomposition (SPOD). Aerospace 2018, 5, 73. [Google Scholar] [CrossRef] [Green Version]
  33. Hu, Z. Research on Intelligent Prediction and Machine Learning Based on Big Data. Digit. Technol. Appl. 2021, 39, 84–86. [Google Scholar]
  34. Ren, F.; Gao, C.; Tang, H. Machine learning for flow control: Applications and development trends. Acta Aerodyn. Sin. 2021, 42, 524686. [Google Scholar]
  35. Lin, Y.; Guan, Z. The Use of Machine Learning for the Prediction of the Uniformity of the Degree of Cure of a Composite in an Autoclave. Aerospace 2021, 8, 130. [Google Scholar] [CrossRef]
  36. Hashemi, S.; Botez, R.; Grigorie, T. New Reliability Studies of Data-Driven Aircraft Trajectory Prediction. Aerospace 2020, 7, 145. [Google Scholar] [CrossRef]
  37. Uzun, M.; Demirezen, M.; Inalhan, G. Physics Guided Deep Learning for Data-Driven Aircraft Fuel Consumption Modeling. Aerospace 2021, 8, 44. [Google Scholar] [CrossRef]
  38. Lerro, A.; Brandl, A.; Battipede, M.; Gili, P. A Data-Driven Approach to Identify Flight Test Data Suitable to Design Angle of Attack Synthetic Sensor for Flight Control Systems. Aerospace 2020, 7, 63. [Google Scholar] [CrossRef]
Figure 1. Flowchart of prediction of the data-driven method.
Figure 1. Flowchart of prediction of the data-driven method.
Aerospace 08 00265 g001
Figure 2. Sampling points.
Figure 2. Sampling points.
Aerospace 08 00265 g002
Figure 3. Finite-element model: (a) flow domain; (b) cylinder.
Figure 3. Finite-element model: (a) flow domain; (b) cylinder.
Aerospace 08 00265 g003
Figure 4. Comparison of density contours.
Figure 4. Comparison of density contours.
Aerospace 08 00265 g004
Figure 5. Predicted and actual pressure coefficients.
Figure 5. Predicted and actual pressure coefficients.
Aerospace 08 00265 g005
Figure 6. Predicted and actual heat fluxes.
Figure 6. Predicted and actual heat fluxes.
Aerospace 08 00265 g006
Figure 7. The contour of structural temperature: (a) actual temperature; (b) predicted temperature.
Figure 7. The contour of structural temperature: (a) actual temperature; (b) predicted temperature.
Aerospace 08 00265 g007
Figure 8. The contour of structural stress: (a) actual stress; (b) predicted stress.
Figure 8. The contour of structural stress: (a) actual stress; (b) predicted stress.
Aerospace 08 00265 g008
Figure 9. Relative error contour of the structural temperature and stress: (a) temperature; (b) stress.
Figure 9. Relative error contour of the structural temperature and stress: (a) temperature; (b) stress.
Aerospace 08 00265 g009
Table 1. State variables and bounds of variables.
Table 1. State variables and bounds of variables.
Mach   Number   ( M a ) Altitude   ( H / k m ) Angle   of   Attack   ( α / ° )
3 M a 8 20 H 40 8 α 8
Table 2. Comparison of numerical simulation in this paper with the experiment.
Table 2. Comparison of numerical simulation in this paper with the experiment.
MethodMaximum Temperature of Structure/K Maximum   Heat   Flux / ( k W / m 2 ) Maximum Pressure of Flow Field/Pa
Numerical Simulation in this paper43666335,386
Experiment [10]46567037,815
Table 3. Test sampling point.
Table 3. Test sampling point.
Test Condition Altitude   ( H / k m ) Mach   Number   ( M a ) Angle   of   Attack   ( α / ° )
1384.2−3.2
2265−6.4
3303.43.2
4346.60
5225.86.4
Table 4. Comparision of efficiency with prediction method and numerical simulation.
Table 4. Comparision of efficiency with prediction method and numerical simulation.
MethodNumber of SnapshotsCPU Time for One Snapshot/hCPU Time for Predicting one Test Condition/s
Prediction method5 12 0.14
Numerical simulation60 12
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, J.; Wang, M.; Li, S. The Rapid Data-Driven Prediction Method of Coupled Fluid–Thermal–Structure for Hypersonic Vehicles. Aerospace 2021, 8, 265. https://doi.org/10.3390/aerospace8090265

AMA Style

Liu J, Wang M, Li S. The Rapid Data-Driven Prediction Method of Coupled Fluid–Thermal–Structure for Hypersonic Vehicles. Aerospace. 2021; 8(9):265. https://doi.org/10.3390/aerospace8090265

Chicago/Turabian Style

Liu, Jing, Meng Wang, and Shu Li. 2021. "The Rapid Data-Driven Prediction Method of Coupled Fluid–Thermal–Structure for Hypersonic Vehicles" Aerospace 8, no. 9: 265. https://doi.org/10.3390/aerospace8090265

APA Style

Liu, J., Wang, M., & Li, S. (2021). The Rapid Data-Driven Prediction Method of Coupled Fluid–Thermal–Structure for Hypersonic Vehicles. Aerospace, 8(9), 265. https://doi.org/10.3390/aerospace8090265

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop