Aircraft aerodynamic force / thermal environment-material ablation response coupling calculation method

Abstract

The invention provides an aircraft aerodynamic force / thermal environment-material ablation response coupling calculation method, and relates to the technical field of aerospace, and the method comprises the steps: dispersing an aircraft reentry trajectory into a plurality of trajectory points, obtaining the trajectory parameter of each trajectory point, determining the temperature and wall surface temperature of an aircraft heat-proof material at the trajectory starting moment, and calculating the temperature of the aircraft heat-proof material at the trajectory starting moment; calculating an initial external flow field of each ballistic point according to the ballistic parameters of each ballistic point, forming initial thermal environment data along the ballistic, calculating the ablation response of the material, obtaining an initial material thermal response calculation result, and calculating the thermal response of the material through multiple coupling iterative calculations. Thermal environment data at each ballistic point moment and a thermal response calculation result at each ballistic point moment are obtained, and when the heat flow relative deviation of the current coupling iteration M and the coupling iteration M-1 does not exceed a threshold value, a convergence standard is reached, and a material ablation response coupling calculation result is obtained. The problem that the material ablation response of the trajectory of the aircraft is difficult to calculate accurately in real time is solved.

Description

Technical Field

[0001] This specification relates to the field of aerospace technology, and in particular to a coupled calculation method for aircraft aerodynamic / thermal environment-material ablation response. Background Technology

[0002] When an aircraft flies at hypersonic speeds, the surrounding gas undergoes intense friction with the aircraft's surface, converting a large amount of the gas's kinetic energy into heat, causing the surface temperature to rise rapidly. For example, under typical hypersonic flight conditions, the temperature at the aircraft's stagnation point can reach 2000K or even higher. Such extreme high-temperature environments place enormous thermal loads on the aircraft's thermal protection system, leading to ablation of the thermal protection materials, including thermochemical ablation of the material surface, high-temperature sublimation, and internal pyrolysis. The ablation process not only alters the material's geometry and physical properties but can also cause structural deformation, seriously threatening the aircraft's flight safety and service life.

[0003] Existing aerodynamic / thermal environment-material ablation response coupling algorithms mostly employ loosely coupled strategies based on flow field and solid heat transfer partitioning calculations. However, these loosely coupled strategies are often unidirectional couplings between the aerodynamic / thermal environment and the material ablation response calculation. That is, during the material ablation response calculation, the aerodynamic / thermal environment at the ablation moment is assumed to be constant, without considering the impact of the material ablation response on the near-wall flow pattern and aerodynamic shape of the aircraft. While this strategy is simple, it violates the coupling principle of mutual correlation and influence between the aerodynamic / thermal environment and the material ablation response, resulting in low algorithm accuracy and often leading to redundant thermal protection for the aircraft. Existing loosely coupled algorithms that consider the bidirectional coupling between aerodynamic / thermal environment and material ablation response often employ a "single-time coupling, point-by-point progression" computational strategy. This involves first calculating the aerodynamic / thermal environment at time T0 and using it as input for material ablation response calculation. The calculation then progresses to the next ablation time T1, using the material ablation response at T1 as input for that time's aerodynamic / thermal environment calculation. This process is repeated until the ablation terminates. While this method considers the bidirectional coupling between aerodynamic / thermal environment and material ablation response, it assumes the aerodynamic / thermal environment remains constant between adjacent time points, resulting in insufficient algorithm accuracy. Furthermore, the algorithm's computational efficiency is low due to the limitation on the time step size, making it unsuitable for the coupled calculation of aerodynamic / thermal environment and material ablation response throughout the entire ballistic trajectory during spacecraft reentry. Summary of the Invention

[0004] To address the aforementioned shortcomings in existing technologies, this invention provides a coupled calculation method for aerodynamic / thermal environment and material ablation response of aircraft, which solves the problem of difficulty in calculating the material ablation response of aircraft trajectories in real time and accurately.

[0005] To achieve the aforementioned objectives, the technical solution adopted by this invention is: a coupled calculation method for aerodynamic / thermal environment and material ablation response of an aircraft, comprising: S1. Discretize the reentry trajectory of the spacecraft into multiple trajectory points, obtain the trajectory parameters of each trajectory point, and determine the temperature of the heat shield material and the wall temperature of the spacecraft at the start of the trajectory based on the trajectory parameters. S2. Calculate the initial external flow field of each ballistic point based on the ballistic parameters of each ballistic point to form the initial thermal environment data along the ballistic trajectory. S3. Using the initial thermal environment data as boundary conditions, and interpolating the wall heat flux and wall pressure between adjacent time points, the initial material thermal response calculation results are obtained. S4. Based on the initial material thermal response calculation results, generate the material ablation response wall boundary conditions, transfer the material ablation retreat amount to the mesh deformation program for mesh update, and calculate the flow field at each ballistic point according to the ballistic parameters to obtain the thermal environment data at each ballistic point. S5. Using thermal environment data as boundary conditions, interpolate the wall heat flux and pressure between adjacent time points to calculate the material thermal response and obtain the thermal response calculation results at each ballistic point. S6. Based on the thermal environment data obtained in S4 and the thermal response calculation results obtained in S5, determine whether the relative deviation of heat flow between the Mth coupling iteration and the (M-1)th coupling iteration exceeds the threshold. If so, return to S4; otherwise, the coupling calculation of the aerodynamic / thermal environment-material ablation response of the aircraft's full ballistic trajectory has reached the convergence condition, the calculation is stopped, and the coupling calculation results of the material ablation response are obtained.

[0006] Further, S2 includes: The wall temperature at the start of the trajectory is T0. w 0=T ini =300K, ablation ejection amount m dot 0=0, the wall temperature at times T1~TN is all set to T w 1~N = Radiation equilibrium wall temperature condition, excitation quantity m dot Set all numbers from 1 to N to 0; Using an aerodynamic / thermal environment CFD calculation program, the initial external flow field of each ballistic point is calculated based on the ballistic parameters of each ballistic point, forming the initial thermal environment data along the ballistic trajectory. The thermal environment data at the ballistic start time T0 remains unchanged during the coupled iteration process.

[0007] Further, S3 includes: Set the temperature of the vehicle's thermal protection material at the ballistic start point T0 as T. s 0=T ini =300K; Using the initial thermal environment data as boundary conditions, and based on the temperature of the spacecraft's heat shield material at the start of the ballistic trajectory, the ablation response calculation program is used to interpolate the wall heat flux and wall pressure between adjacent time points to calculate the ablation response of the material, thus obtaining the initial material thermal response calculation result. The initial material thermal response calculation result remains unchanged during the coupled iteration process.

[0008] Furthermore, the step of using the initial thermal environment data as boundary conditions to interpolate the wall heat flux and wall pressure between adjacent time points, calculating the ablation response of the material, and obtaining the initial material thermal response calculation results includes: The thermal environment data, i.e., the heat flux Q of the aircraft wall, is received from the initial flow field calculation by the material ablation response calculation program. w Wall pressure P w As boundary conditions, and by interpolating the wall heat flux and wall pressure between adjacent time points, the ablation response of the material is calculated, yielding the ablation wall temperature, ablation component mass fraction, ablation mass flow rate, and ablation retreat at each trajectory point. The material wall temperature T at the trajectory start time t0 is also calculated. w 0=T ini =300K, ejection mass m dot 0=0, mass fraction of wall ablation component S0=0, ablation retreat amount a bl 0 = 0, and remains unchanged in subsequent iterations.

[0009] By considering the bidirectional coupling between the aerodynamic / thermal environment of the aircraft and the ablation response of the material, and by using an interpolation algorithm to interpolate the aerodynamic and thermal information between adjacent ballistic points, the calculation of the ablation response of the material is closer to the result of the tightly coupled calculation. Therefore, the coupling algorithm provided by this invention has higher calculation accuracy than the existing unidirectional coupling algorithm and single coupling algorithm.

[0010] Further, S4 includes: Based on the initial material ablation response calculation results, the material ablation response wall boundary conditions are generated, and the material ablation retreat is transferred to the mesh deformation program for mesh update. Based on the mesh update results, the flow field of each ballistic point is recalculated according to the ballistic parameters to obtain the thermal environment data of each ballistic point at any given time.

[0011] Furthermore, the process of updating the computational grid specifically includes: The aerodynamic / thermal environment CFD calculation program receives the initial material thermal response calculation results, namely ablation wall temperature, mass fraction of ablation components on the wall, ablation mass flow rate, and ablation retreat, generates the material ablation response wall boundary conditions, and transmits the material ablation retreat to the mesh deformation program for mesh update calculation.

[0012] Further, S5 includes: The material ablation response calculation program receives thermal environment data, namely wall heat flux and wall pressure, as boundary conditions and performs interpolation processing on the wall heat flux and pressure between adjacent time points to calculate the material thermal response and obtain the thermal response calculation results at each ballistic point.

[0013] Further, S6 includes: The aerodynamic / thermal environment CFD calculation program receives the material thermal response calculation results of the M-1th iteration, namely the ablation wall temperature, the mass fraction of the wall ablation component, the ablation mass flow rate, and the ablation retreat, generates the wall boundary conditions of the material ablation response, and transmits the material ablation retreat to the mesh deformation program for calculation mesh update. It also calculates the flow field of each ballistic point according to the ballistic parameter conditions to obtain the thermal environment data of each ballistic point at the Mth iteration. The material ablation response calculation program receives the thermal environment data of the Mth iteration, namely the wall heat flow and wall pressure, as boundary conditions and interpolates the wall heat flow and pressure between adjacent time points to calculate the material thermal response and obtain the thermal response calculation results of each ballistic point in the Mth iteration. The relative deviation of heat flow between the Mth iteration and the (M-1)th iteration is analyzed. When the relative deviation of heat flow between the Mth iteration and the (M-1)th iteration does not exceed the threshold, the coupling calculation of the aerodynamic / thermal environment-material ablation response of the entire ballistic trajectory of the aircraft reaches the convergence criterion, the iteration is stopped, the coupling calculation result of material ablation response is obtained, and the coupling calculation of aerodynamic / thermal environment-material ablation response of the aircraft is completed.

[0014] The beneficial effects of the present invention are as follows: The present invention provides a coupling calculation method for aerodynamic / thermal environment-material ablation response of an aircraft. (1) It adopts a convergence method of iterative calculation of the entire trajectory, that is, it can realize the coupling calculation of aerodynamic / thermal-material ablation response of all ballistic points of the entire trajectory in one go, which can achieve higher calculation efficiency; (2) By coupling and iterating the thermal environment data and thermal response calculation results of each ballistic point, the coupling calculation of material ablation response of the aircraft trajectory is realized, which makes the construction of material ablation response coupling edge more comprehensive and has higher coupling accuracy. Attached Figure Description

[0015] This specification will be further described by way of exemplary embodiments, which will be described in detail with reference to the accompanying drawings. These embodiments are not limiting; in these embodiments, the same reference numerals denote the same structures, wherein: Figure 1 This is an exemplary flowchart illustrating a coupled calculation method for aerodynamic / thermal environment-material ablation response of an aircraft, according to some embodiments of this specification. Figure 2 This is an exemplary flowchart illustrating an iterative calculation method for the coupled aerodynamic / thermal-material ablation response of a hypersonic reentry vehicle throughout its entire trajectory, based on some embodiments of this specification. Figure 3 This is an exemplary schematic diagram of the inputs and outputs of each calculation stage of a coupled calculation method for the full ballistic aerodynamic / thermal-material ablation response of a hypersonic reentry vehicle, as shown in some embodiments of this specification. Figure 4 This is an exemplary schematic diagram of the reentry trajectory and computational grid division of the Stardust return capsule in the calculation working condition section, according to some embodiments of this specification. Figure 5 This is an exemplary schematic diagram illustrating the calculation of heat flow on the ablation-iterative wall surface of each wheel at the T34 ballistic point of the Stardust return capsule, according to some embodiments of this specification. Figure 6 This is an exemplary schematic diagram illustrating the calculation of the ablation iteration wall temperature of each wheel at the T34 ballistic point of the Stardust return capsule according to some embodiments of this specification; Figure 7 This is an exemplary schematic diagram illustrating the calculation of heat flow on the ablation-iterative wall surface of each wheel at the T42 ballistic point of the Stardust return capsule, according to some embodiments of this specification. Figure 8 This is an exemplary schematic diagram illustrating the calculation of the ablation iteration wall temperature of each wheel at the T42 ballistic point of the Stardust return capsule according to some embodiments of this specification; Figure 9 This is an exemplary schematic diagram illustrating the calculation of heat flow on the ablation-iterative wall surface of each wheel at the T48 trajectory point of the Stardust return capsule, according to some embodiments of this specification. Figure 10 This is an exemplary schematic diagram illustrating the calculation of the ablation iteration wall temperature of each wheel at the T48 ballistic point of the Stardust return capsule according to some embodiments of this specification; Figure 11 This is an exemplary schematic diagram illustrating the calculation of heat flow on the ablation-iterative wall surface of each wheel at the T54 ballistic point of the Stardust return capsule, according to some embodiments of this specification. Figure 12 This is an exemplary schematic diagram illustrating the calculation of the ablation iteration wall temperature of each wheel at the T54 ballistic point of the Stardust return capsule according to some embodiments of this specification; Figure 13 This is an exemplary schematic diagram illustrating the calculation of heat flow on the ablation-iterative wall surface of each wheel at the T60 trajectory point of the Stardust return capsule, according to some embodiments of this specification. Figure 14 This is an exemplary schematic diagram illustrating the calculation of the ablation iteration wall temperature of each wheel at the T60 ballistic point of the Stardust return capsule according to some embodiments of this specification; Figure 15 This is an exemplary schematic diagram illustrating the calculation of heat flow on the ablation-iterative wall surface of each wheel at the T66 ballistic point of the Stardust return capsule, according to some embodiments of this specification. Figure 16 This is an exemplary schematic diagram illustrating the calculation of the ablation iteration wall temperature of each wheel at the T66 ballistic point of the Stardust return capsule according to some embodiments of this specification; Figure 17 This is an exemplary schematic diagram illustrating the calculation of heat flow on the ablation-iterative wall surface of each wheel at the T76 ballistic point of the Stardust return capsule, according to some embodiments of this specification. Figure 18 This is an exemplary schematic diagram illustrating the calculation of the ablation iteration wall temperature of each wheel at the T76 ballistic point of the Stardust return capsule according to some embodiments of this specification; Figure 19 This is an exemplary schematic diagram showing the ablation profile busbar receding over time according to some embodiments of this specification. Detailed Implementation

[0016] The specific embodiments of the present invention are described below to enable those skilled in the art to understand the present invention. However, it should be understood that the present invention is not limited to the scope of the specific embodiments. For those skilled in the art, various changes are obvious as long as they are within the spirit and scope of the present invention as defined and determined by the appended claims. All inventions utilizing the concept of the present invention are protected.

[0017] Example Figure 1 This is an exemplary flowchart illustrating a coupled calculation method for aircraft aerodynamic / thermal environment-material ablation response, based on some embodiments of this specification. Figure 1 As shown, the process includes the following steps.

[0018] In some embodiments, the aerodynamic / thermal environment-material ablation response coupling calculation method for the full trajectory of the Stardust reentry capsule configuration can be used to perform aerodynamic / thermal environment-material ablation response coupling calculations on the Stardust reentry capsule configuration.

[0019] S1: Discretize the reentry trajectory of the spacecraft into multiple trajectory points, obtain the trajectory parameters of each trajectory point, and determine the temperature of the heat shield material and the wall temperature of the spacecraft at the start of the trajectory based on the trajectory parameters.

[0020] The ballistic parameters at each trajectory point are discrete parameters of the hypersonic reentry vehicle's entire trajectory. For example, the trajectory parameters may include altitude, Mach number, and angle of attack, and the trajectory may include a portion of the Stardust reentry capsule's reentry trajectory. Specifically, the calculation conditions (partial reentry trajectory) and mesh generation for the Stardust reentry capsule are as follows: Figure 4 As shown.

[0021] In some embodiments, the reentry trajectory of the aircraft can be discretized into N+1 trajectory points (i.e., T0~TN) to obtain the trajectory parameters of the N+1 trajectory points and the temperature T of the aircraft's heat shield material at the trajectory start time T0. S 0 and wall temperature T w 0.

[0022] In some embodiments, the wall ablation ejection amount m at the start of the trajectory can be considered as... dot 0 is negligible, that is, mdot0 = 0.

[0023] S2: Calculate the initial external flow field of each ballistic point based on the ballistic parameters of each ballistic point, forming the initial thermal environment data along the ballistic trajectory, and the thermal environment data at the start of the ballistic trajectory T0 remains unchanged during the coupling iteration process.

[0024] The initial thermal environment data along the trajectory is the thermal environment data of the initial external flow field at each trajectory point.

[0025] In some embodiments, the wall temperature T at the start of the trajectory can be... w 0=T ini =300K, ablation ejection amount m dot 0=0, the wall temperature at times T1~TN is all set to T w 1~N=300K, ejection quantity m dot All numbers from 1 to N are set to 0. Using the aerodynamic / thermal environment CFD calculation program, the initial external flow field of each ballistic point is calculated based on the ballistic parameters of each ballistic point, forming the initial thermal environment data along the ballistic trajectory. The thermal environment data at the ballistic start time T0 remains unchanged during the coupled iteration process.

[0026] In some embodiments, the wall temperature T at the start of the trajectory can be... w 0=T ini =300K, ablation ejection amount m dot 0=0, the wall temperature at times T1~TN is all set to T w 1~N=300K (or radiation equilibrium temperature), excitation quantity mdot All numbers from 1 to N are set to 0; the initial external flow field at each trajectory point is calculated by the aerodynamic / thermal environment CFD calculation program based on the trajectory parameters at each trajectory point, forming the initial thermal environment data (t0, Q) along the trajectory. w 0, P w 0) ini 、(t1、Q w 1. P w 1) ini ... (tN, Q) w N, P w N) ini , where (tN, Q) w N, P w N) ini This represents the initial wall heat flux and wall pressure of the reentry capsule at the Nth trajectory point, and the thermal environment data (t0, Q) at the start of the trajectory at time T0. w 0, P w 0) ini It remains unchanged in subsequent iterations.

[0027] S3: The initial thermal environment data is used as the boundary condition, and the wall heat flux and wall pressure between adjacent time points are interpolated to obtain the initial material thermal response calculation results. The initial material thermal response calculation results remain unchanged during the coupled iteration process.

[0028] The initial material thermal response calculation results are the calculation results of the initial ablation response of the material.

[0029] The external flow field of a ballistic trajectory is the external flow field data of the aircraft's trajectory.

[0030] In some embodiments, the temperature of the vehicle's thermal protection material at the ballistic start time T0 can be set to T. s 0=T ini =300K; Using the initial thermal environment data along the trajectory as boundary conditions, based on the temperature of the spacecraft's heat shield material at the start of the trajectory, the wall heat flux and wall pressure between adjacent time points are interpolated using a material ablation response calculation program to calculate the ablation response of the material and obtain the initial material thermal response calculation results.

[0031] In some embodiments, the temperature of the vehicle's thermal protection material at the ballistic start time T0 can be set to T. s 0=T ini =300K; The thermal environment data, i.e., the heat flux Q of the aircraft wall, is received from the initial flow field calculation program by the material ablation response calculation program. w Wall pressure P wAs boundary conditions, and by interpolating the wall heat flux and wall pressure between adjacent time points, the ablation response of the material is calculated, yielding the ablation wall temperature, ablation component mass fraction, ablation mass flow rate, and ablation retreat (t1, T) at each trajectory point. w 1, S1, m dot 1. a bl 1) ini ... (tN, T) w N, SN, m dot N, a bl N) ini , where (tN, T) w N, SN, m dot N, a bl N) ini , respectively, represent the ablation wall temperature, ablation component mass fraction, ablation mass flow rate, and ablation retreat amount calculated from the ablation response of the return capsule at time N, and the material wall temperature T at the start of the trajectory t0. w 0=T ini =300K, ejection mass m dot 0=0, mass fraction of wall ablation component S0=0, ablation retreat amount a bl 0 = 0, and remains unchanged in subsequent iterations.

[0032] S4: Based on the initial material thermal response calculation results, generate the material ablation response wall boundary conditions, transfer the material ablation retreat amount to the mesh deformation program for mesh update, and calculate the flow field at each ballistic point according to the ballistic parameters to obtain the thermal environment data at each ballistic point.

[0033] In some embodiments, based on the initial material ablation response calculation results, material ablation response wall boundary conditions are generated, and the material ablation retreat is transmitted to the mesh deformation program for mesh update calculation; based on the mesh update results, the flow field of each ballistic point is recalculated according to the ballistic parameters to obtain the thermal environment data of each ballistic point at any given time.

[0034] In some embodiments, such as Figure 2 As shown, the initial material thermal response calculation results, namely ablation wall temperature, mass fraction of ablation components on the wall surface, ablation mass flow rate, and ablation retreat, can be received by the aerodynamic / thermal environment CFD calculation program. This generates the material ablation response wall boundary conditions, and the material ablation retreat is transmitted to the mesh deformation program for mesh updating. Then, the flow field at each trajectory point is recalculated based on the ballistic parameters, obtaining the thermal environment data (t1, Q) at each trajectory point during the first iteration. w 1. P w 1) iter1 ...(tN, Q) w N, P w N)iter1 .

[0035] In some embodiments, such as Figure 3 As shown, the aerodynamic / thermal environment CFD calculation program can receive the material thermal response calculation results from the first iteration, namely the ablation wall temperature, the mass fraction of ablation components on the wall surface, the ablation mass flow rate, and the ablation retreat, to generate the material ablation response wall boundary conditions. The material ablation retreat is then transferred to the mesh deformation program for mesh updating. Afterward, the flow field at each trajectory point is recalculated based on the ballistic parameters, yielding the thermal environment data (t1, Q) at each trajectory point in the second iteration. w 1. P w 1) iter2 ...(tN, Q) w N, P w N) iter2 .

[0036] S5: Using thermal environment data as boundary conditions, interpolate the wall heat flux and pressure between adjacent time points to calculate the material thermal response and obtain the thermal response calculation results at each trajectory point.

[0037] The thermal environment data at each ballistic point in the first iteration is the thermal environment-related data of the return capsule at the ballistic point after the first iteration. For example, the thermal environment data of the first iteration may include the wall heat flux and the wall pressure of the first iteration.

[0038] In some embodiments, such as Figure 2 As shown, the material ablation response calculation program can receive the thermal environment data from the first iteration, namely the wall heat flux and wall pressure, as boundary conditions, and interpolate the wall heat flux and pressure between adjacent time points to calculate the material thermal response, obtaining the thermal response calculation results for each trajectory point in the first iteration: ablation wall temperature, wall ablation component mass fraction, ablation mass flow rate, and ablation retreat (t1, T). w 1, S1, m dot 1. a bl 1) iter1 ... (tN, T) w N, SN, m dot N, a bl N) iter1 .

[0039] In some embodiments, the thermal environment data at each ballistic point can be used as boundary conditions to interpolate the wall heat flux and pressure between adjacent time points, calculate the material thermal response, and obtain the thermal response calculation results at each ballistic point.

[0040] In some embodiments, the material ablation response calculation program can receive the second iteration thermal environment data, namely wall heat flux and wall pressure, as boundary conditions and perform interpolation processing on the wall heat flux and pressure between adjacent time points to calculate the material thermal response, obtaining the thermal response calculation results at each trajectory point in the second iteration: ablation wall temperature, wall ablation component mass fraction, ablation mass flow rate, and ablation retreat amount (t1, T). w 1, S1, m dot 1. a bl 1) iter2 ... (tN, T) w N, SN, m dot N, a bl N) ite2 .

[0041] S6: Based on the thermal environment data and thermal response calculation results at each ballistic point, determine whether the relative deviation of heat flow between the Mth coupling iteration and the (M-1)th coupling iteration exceeds the threshold. If yes, return to S4; otherwise, stop the calculation when the coupling calculation of the aerodynamic / thermal environment-material ablation response of the entire ballistic trajectory of the aircraft reaches the convergence condition, obtain the coupling calculation result of the material ablation response, and complete the coupling calculation.

[0042] The thermal environment data obtained at each ballistic point moment through iteration is the thermal environment-related data of the return capsule at the ballistic point moment relative to the previous iteration process.

[0043] The thermal response calculation results obtained at each ballistic point are the results of the ablation response of the return capsule at the ballistic point relative to the previous iteration.

[0044] The relative deviation of heat flow is the deviation of heat flow calculated from thermal environment data between adjacent iterations.

[0045] In some embodiments, the processor can calculate the relative heat flux deviation based on the thermal environment data and thermal response calculation results at each ballistic point. For example, the formula for calculating the relative heat flux deviation can be: Relative heat flux deviation = |current heat flux - previous heat flux| / current heat flux; where current heat flux can include the thermal environment data and thermal response calculation results at each ballistic point in the current iteration, and previous heat flux can include the thermal environment data and thermal response calculation results at each ballistic point in the previous iteration.

[0046] The results of the coupled calculation of the material ablation response are the results of the coupled calculation of the ablation response of the return capsule at the final ballistic point.

[0047] In some embodiments, the aerodynamic / thermal environment CFD calculation program can receive the material thermal response calculation results of the (M-1)th iteration, namely the ablation wall temperature, the mass fraction of ablation components on the wall surface, the ablation mass flow rate, and the ablation retreat, generate the material ablation response wall boundary conditions, and transmit the material ablation retreat to the mesh deformation program for computational mesh update. The flow field at each trajectory point is then recalculated based on the ballistic parameter conditions to obtain the thermal environment data (t1, Q) at each trajectory point in the Mth iteration. w 1. P w 1) iterM ...(tN, Q) w N, P w N) iterM .

[0048] In some embodiments, the material ablation response calculation program can receive the thermal environment data of the Mth iteration, namely the wall heat flow and wall pressure, as boundary conditions, and perform interpolation processing on the wall heat flow and pressure between adjacent time points to calculate the material thermal response, obtaining the thermal response calculation results at each trajectory point in the Mth iteration: ablation wall temperature, wall ablation component mass fraction, ablation mass flow rate, and ablation retreat amount (t1, T). w 1, S1, m dot 1. a bl 1) iterM ... (tN, T) w N, SN, m dot N, a bl N) iterM .

[0049] In some embodiments, the relative deviation of heat flow between the Mth iteration and the (M-1)th iteration can be analyzed. When the relative deviation of heat flow between the Mth iteration and the (M-1)th iteration does not exceed 5%, the coupling calculation of the aerodynamic / thermal environment-material ablation response of the entire ballistic trajectory of the aircraft reaches the convergence criterion, the iteration is stopped, the coupling calculation result of the material ablation response is obtained, and the coupling calculation of the aerodynamic / thermal environment-material ablation response of the aircraft is completed.

[0050] In some embodiments, the initial CFD calculation results and the wall heat flux and wall temperature from each CFD-ablation coupled iteration are extracted and analyzed, such as... Figures 5 to 18 As shown. After four rounds of coupled iterations, the calculated results of the total wall heat flux and wall temperature for all working conditions all tended to stabilize. The maximum CFD heat flux error between the fourth and third iterations for each working condition was less than 5%, reaching the convergence criterion. Further, the ablation profile generatrix of the calculation model corresponding to each ballistic point working condition in the final convergence result (fourth iteration result) was extracted, and the results are shown below. Figure 19As shown in the figure, the ablation profile of the model moves backward as the ablation time progresses, with the largest ablation at the stagnation point, approximately 5 mm; the ablation gradually decreases from the stagnation point to the shoulder of the return capsule.

[0051] In some embodiments of this specification, a method for coupled calculation of aerodynamic / thermal environment-material ablation response of an aircraft is provided. (1) A convergence method of iterative calculation of the entire trajectory is adopted, that is, the coupled calculation of aerodynamic / thermal-material ablation response of all ballistic points of the entire trajectory can be realized at one time, which can achieve higher calculation efficiency; (2) The coupled calculation of material ablation response of the aircraft trajectory is realized, which makes the construction of material ablation response coupling strip more comprehensive and has higher coupling accuracy.

Claims

1. A coupled calculation method for aerodynamic / thermal environment and material ablation response of an aircraft, characterized in that, include: S1. Discretize the reentry trajectory of the spacecraft into multiple trajectory points, obtain the trajectory parameters of each trajectory point, and determine the temperature of the heat shield material and the wall temperature of the spacecraft at the start of the trajectory based on the trajectory parameters. S2. Calculate the initial external flow field of each ballistic point based on the ballistic parameters of each ballistic point of the aircraft, and form the initial thermal environment data along the trajectory. S3. Using the initial thermal environment data as boundary conditions, and interpolating the wall heat flux and pressure, the initial material thermal response calculation results are obtained. S4. Based on the initial material thermal response calculation results, generate the material ablation response wall boundary conditions, update the computational mesh, calculate the flow field at each trajectory point according to the ballistic parameters, and obtain the thermal environment data of each trajectory point at any time. S5. Using thermal environment data as boundary conditions, interpolate the wall heat flux and pressure to calculate the material thermal response and obtain the thermal response calculation results at each ballistic point. S6. Based on the thermal environment data obtained in S4 and the thermal response calculation results obtained in S5, determine whether the relative deviation of heat flow between the Mth coupling iteration and the (M-1)th coupling iteration exceeds the threshold. If so, return to S4. Otherwise, the coupled calculation of the aerodynamic / thermal environment and material ablation response of the entire ballistic trajectory of the aircraft reaches the convergence condition, the calculation stops, and the results of the coupled calculation of material ablation response are obtained.

2. The coupled calculation method for aerodynamic / thermal environment-material ablation response of an aircraft according to claim 1, characterized in that, S2 includes: The wall temperature at the start of the trajectory is T0. w 0=T ini =300K, ablation ejection amount m dot 0=0, the wall temperature at times T1~TN is all set to T w 1~N = Radiation equilibrium wall temperature condition, excitation quantity m dot Set all numbers from 1 to N to 0; Using an aerodynamic / thermal environment CFD calculation program, the initial external flow field of each ballistic point is calculated based on the ballistic parameters of each ballistic point, forming the initial thermal environment data along the ballistic trajectory. The thermal environment data at the ballistic start time T0 remains unchanged during the coupled iteration process.

3. The coupled calculation method for aerodynamic / thermal environment-material ablation response of an aircraft according to claim 1, characterized in that, S3 includes: Set the temperature of the vehicle's thermal protection material at the ballistic start point T0 as T. s 0=T ini =300K; Using the initial thermal environment data as boundary conditions, interpolation is performed on the wall heat flux and wall pressure between adjacent time points to calculate the ablation response of the material, thereby obtaining the initial material thermal response calculation result. The initial material thermal response calculation result remains unchanged during the coupled iteration process.

4. The coupled calculation method for aerodynamic / thermal environment-material ablation response of an aircraft according to claim 3, characterized in that, The step of using the initial thermal environment data as boundary conditions to interpolate the wall heat flux and wall pressure between adjacent time points, calculating the ablation response of the material, and obtaining the initial material thermal response calculation results includes: The thermal environment data, i.e., the heat flux Q of the aircraft wall, is received from the initial flow field calculation by the material ablation response calculation program. w Wall pressure P w As boundary conditions, and by interpolating the wall heat flux and wall pressure between adjacent time points, the ablation response of the material is calculated, yielding the ablation wall temperature, ablation component mass fraction, ablation mass flow rate, and ablation retreat at each trajectory point. The material wall temperature T at the trajectory start time t0 is also calculated. w 0=T ini =300K, ejection mass m dot 0=0, mass fraction of wall ablation component S0=0, ablation retreat amount a bl 0 = 0, and remains unchanged in subsequent iterations.

5. The coupled calculation method for aerodynamic / thermal environment-material ablation response of an aircraft according to claim 1, characterized in that, S4 includes: Based on the initial material ablation response calculation results, the material ablation response wall boundary conditions are generated, and the material ablation retreat is transferred to the mesh deformation program for mesh update. Based on the mesh update results, the flow field of each ballistic point is recalculated according to the ballistic parameters to obtain the thermal environment data of each ballistic point at any given time.

6. The coupled calculation method for aerodynamic / thermal environment-material ablation response of an aircraft according to claim 5, characterized in that, The process of updating the computational grid specifically includes: The aerodynamic / thermal environment CFD calculation program receives the initial material thermal response calculation results, namely ablation wall temperature, mass fraction of ablation components on the wall, ablation mass flow rate, and ablation retreat, generates the material ablation response wall boundary conditions, and transmits the material ablation retreat to the mesh deformation program for mesh update calculation.

7. The coupled calculation method for aerodynamic / thermal environment-material ablation response of an aircraft according to claim 5, characterized in that, S5 includes: The material ablation response calculation program receives thermal environment data, namely wall heat flux and wall pressure, as boundary conditions and performs interpolation processing on the wall heat flux and pressure between adjacent time points to calculate the material thermal response and obtain the thermal response calculation results at each ballistic point.

8. The coupled calculation method for aerodynamic / thermal environment-material ablation response of an aircraft according to claim 1, characterized in that, S6 includes: The aerodynamic / thermal environment CFD calculation program receives the material thermal response calculation results of the M-1th iteration, namely the ablation wall temperature, the mass fraction of the wall ablation component, the ablation mass flow rate, and the ablation retreat, generates the wall boundary conditions of the material ablation response, and transmits the material ablation retreat to the mesh deformation program for calculation mesh update. It also calculates the flow field of each ballistic point according to the ballistic parameter conditions to obtain the thermal environment data of each ballistic point at the Mth iteration. The material ablation response calculation program receives the thermal environment data of the Mth iteration, namely the wall heat flow and wall pressure, as boundary conditions and interpolates the wall heat flow and pressure between adjacent time points to calculate the material thermal response and obtain the thermal response calculation results of each ballistic point in the Mth iteration. The relative deviation of heat flow between the Mth iteration and the (M-1)th iteration is analyzed. When the relative deviation of heat flow between the Mth iteration and the (M-1)th iteration does not exceed the threshold, the coupling calculation of the aerodynamic / thermal environment-material ablation response of the entire ballistic trajectory of the aircraft reaches the convergence criterion, the iteration is stopped, the coupling calculation result of material ablation response is obtained, and the coupling calculation of aerodynamic / thermal environment-material ablation response of the aircraft is completed.

Images