Aircraft wall surface ablation response data processing method based on interpolation algorithm

Abstract

The invention discloses an aircraft wall surface ablation response data processing method based on an interpolation algorithm, and relates to the field of ablation response data processing, and the method comprises the following steps: laying ablation sites, and outputting the initial coordinate information of the ablation sites; acquiring coordinate information of a grid center point of an ablated wall surface grid, and calculating coordinates of an ablated site after ablation retreating in combination with the initial coordinate information of the ablated site; ablative information at the ablative site is obtained, interpolation is conducted on the ablative information based on the ablative site coordinates after ablative retreating, and ablative information after interpolation is obtained; and performing wall mapping on the interpolated ablation information to obtain continuous distribution of the material ablation response on the whole ablation wall. According to the method, interpolation processing of various kinds of ablation information in the ablation sites distributed on the complex ablation profile in a limited mode can be achieved, and compared with a typical two-dimensional or three-dimensional interpolation algorithm, the method has the advantages of being simple in algorithm, high in applicability, high in calculation efficiency and the like.

Description

Technical Field

[0001] This invention relates to the field of ablation response data processing, and more specifically to a method for processing ablation response data of aircraft walls based on interpolation algorithms. Background Technology

[0002] Aerodynamic / thermal environment-material ablation response coupled calculations involve complex multi-domain, multi-field, and multi-scale coupled computations. During spacecraft reentry, the surface is subjected to prolonged aerodynamic / thermal loads due to the intense compression of shock waves and the viscous stagnation effect of walls, resulting in intense thermochemical ablation of the surface materials. On one hand, the material ablation response, including ablation gas mass ejection, ablation surface retreat, surface ablation products, and ablation temperature distribution, alters the near-wall flow state and the aerodynamic shape of the spacecraft, thereby changing the aerodynamic / thermal environment. On the other hand, the altered aerodynamic / thermal environment also affects the material ablation response. Therefore, to simplify the computational process and improve the efficiency of coupled calculations, a loosely coupled calculation method based on partitioned calculations of aerodynamic / thermal environment and material ablation response, and bidirectional coupling of aerodynamic / thermal and material ablation response, is often adopted. In the loosely coupled computational process, a crucial step is constructing ablation response boundary conditions suitable for CFD (Computational Fluid Dynamics) to incorporate the material's ablation response into the aerodynamic / thermal environment calculations of the aircraft. Currently, common methods involve constructing wall mass ejection boundary conditions based on the momentum / mass conservation relationship between the ablated wall and the first-layer computational grid, characterizing the distribution of ablation components on the wall using fixed boundary conditions, and employing dynamic meshing techniques to capture the ablation retreat profile and update the flow domain computational grid. The establishment of these ablation response boundary conditions requires the continuous distribution of material ablation response information across the entire ablated wall. Since most theoretically mature and widely applied ablation calculation methods are one-dimensional, meaning they can only calculate the material ablation response at a single site at a time, and considering the efficiency of two-way coupled computation, the number of ablation sites must be limited. Therefore, interpolation algorithms are needed to interpolate the ablation information from the finite number of scattered sites to obtain a continuous distribution of the wall material's ablation response.

[0003] The ablation calculation sites on the ablation walls of typical aircraft are usually a finite number of points scattered on a three-dimensional curved surface in Cartesian coordinate space. Because the ablation wall varies greatly with the aerodynamic shape of the aircraft and recedes in real time during the ablation process, its numerical expression is difficult to determine. Therefore, traditional two-dimensional or three-dimensional linear interpolation algorithms cannot obtain continuously distributed wall ablation information. Furthermore, the types of ablation information that need to be interpolated at each ablation site are numerous, including ablation site coordinates, wall mass flow rate, wall ablation temperature, wall ejection velocity, and the mass distribution of wall ablation components (including various thermochemical ablation products of carbon and various phenolic resin pyrolysis products). Traditional high-dimensional linear interpolation algorithms have significant shortcomings in computational efficiency; the time and resources required for interpolation calculations increase dramatically with the number of data points processed and the amount of interpolated information, greatly limiting their application potential in the refined interpolation calculation of ablation information on large aircraft walls. Summary of the Invention

[0004] To address the aforementioned shortcomings in existing technologies, the interpolation algorithm-based aircraft wall ablation response data processing method provided by this invention solves the problem that existing two-dimensional or three-dimensional linear interpolation algorithms cannot obtain continuously distributed wall ablation information.

[0005] To achieve the above-mentioned objectives, the technical solution adopted by this invention is as follows: A method for processing ablation response data of aircraft walls based on interpolation algorithms is provided, which includes the following steps: Determine the ablation sites and output the initial coordinate information of the ablation sites; Obtain the coordinate information of the grid center point of the ablation wall surface, and calculate the coordinates of the ablation site after the ablation retreat by combining the initial coordinate information of the ablation site. Obtain ablation information at the ablation site, and interpolate the ablation information based on the coordinates of the ablation site after ablation retreat to obtain the interpolated ablation information. The interpolated ablation information is mapped onto the wall surface to obtain the continuous distribution of the material ablation response across the entire ablated wall surface.

[0006] Furthermore, the specific methods for deploying ablation sites include: For the initial blunt-nosed body configuration and typical transition ablation morphologies at different times, a mesh drawing tool was used to arrange six groups of ablation sites along the x-axis for the initial blunt-nosed body configuration and the ablation morphologies at different times during typical transition ablation. The first group had one ablation site, located at the apex of the blunt-nosed body; the remaining groups each had six ablation sites, and the x-axis coordinates of the ablation sites within each group were kept consistent. After the points were arranged, the coordinate information of the ablation sites was output in groups, generating the ablation_position.dat file.

[0007] Furthermore, specific methods for obtaining the coordinate information of the grid center points on the ablated wall surface include: Read the grid center point position coordinate information ablation_position.dat file generated by the mesh drawing software, and read the grid center point coordinates and normal vector information cellcenter_message.dat file given by the ablation wall mesh CFD program for calculating the ablation environment, to obtain the grid center point coordinates and normal vector of the ablation wall mesh.

[0008] Furthermore, the specific method for calculating the coordinates of the ablation site after retreat, based on the initial coordinate information of the ablation site, includes the following steps: Based on the principle of minimizing the Euclidean distance between the grid center and the ablation site on the ablation wall, and using the formula:

[0009] Calculate the Euclidean distance between the current ablation site and the center points of all ablation wall grids. ;in The coordinates of the center point of any ablation wall grid; The coordinates of the current ablation site; Based on the Euclidean distance between the current ablation site and the center points of all ablation wall grids, determine the center point of the ablation wall grid that is closest to the current ablation site, and use the normal vector information at the center point of the ablation wall grid as the initial normal vector of the current ablation site. Based on the initial normal vector of the current ablation site, the coordinates of the ablation site after retreat are calculated, and the expression is as follows:

[0010] in The coordinates of the ablation site after the ablation retreat; The initial coordinates of the current ablation site; This is the initial normal vector for the current ablation site; This represents the ablation retreat at the current ablation site. Output the coordinates of the ablation site after the ablation retreat as a file named new_ablation_position.dat.

[0011] Furthermore, the specific method for obtaining ablation information at the ablation site is as follows: Ablation information at the ablation site is obtained by reading the ablation information file ablation_calculation.dat output by the material ablation analysis software.

[0012] Furthermore, specific methods for interpolating ablation information include: Interpolation is performed only in the x-direction for the ablation information of the two-dimensional curve; interpolation is then performed in the z-direction and x-direction for the ablation information of the three-dimensional surface.

[0013] Furthermore, the specific method for interpolating the ablation information of the two-dimensional curve only in the x-direction includes the following steps: Determine the projection plane and project all the coordinates of the ablation sites after ablation retreat onto the xz plane; assign the ablation information at the ablation site to the corresponding ablation site after ablation retreat. Determine the number of interpolation points between two adjacent ablation retreat sites within the same group and perform linear interpolation of the ablation information; that is, according to the formula or A linear interpolation of ablation information is performed sequentially along the x-direction between two adjacent ablation retreat sites within the group; where... This is the result of a linear interpolation. and These are the ablation information of two adjacent ablation retreat sites within the group; and It is a linear interpolation factor; Output the interpolation results of the ablation sites in the projection plane after ablation and generate the interpolation_ablation_message.dat file to obtain the ablation information of each interpolation point, that is, the ablation information after two-dimensional curve interpolation.

[0014] Furthermore, the specific method for interpolating the three-dimensional surface ablation information sequentially in the z and x directions includes the following steps: Determine the projection plane and project all the coordinates of the ablation sites after ablation retreat onto the xz plane; assign the ablation information at the ablation site to the corresponding ablation site after ablation retreat. Determine the number of groups, the number of ablation retreat sites in the same group, the number of interpolations between two adjacent ablation retreat sites in the x-direction, and the number of interpolations between two adjacent ablation retreat sites in the z-direction. The ablation information of the ablation retreat and ablation site group at the end of the aircraft is copied and expanded to be consistent with the number of other ablation retreat and ablation site groups. For each group of ablation sites after retreat, a linear interpolation of the ablation information within the group is performed, i.e., according to the formula... or A linear interpolation of ablation information is performed sequentially along the z-direction between two adjacent ablation sites after retreat within each group; where... This is the result of a linear interpolation. and These are the ablation information of two adjacent ablation retreat sites within the group; and It is a linear interpolation factor; Along the x-direction, after a linear interpolation, a second interpolation of the ablation information is performed between the corresponding ablation sites in the adjacent two groups of ablation sites. The interpolation results of the ablation sites in the projection plane are output and the interpolation_ablation_message.dat file is generated to obtain the ablation information of each interpolation point, that is, the ablation information after three-dimensional surface interpolation.

[0015] Furthermore, , Pi , The coordinates of the current point to be interpolated along the interpolation direction; The coordinates of the first ablation retreat ablation site among two adjacent ablation retreat ablation sites; The coordinates are the coordinates of the second ablation retreat ablation site among two adjacent ablation retreat ablation sites.

[0016] Furthermore, specific methods for wall mapping of the interpolated ablation information include: The coordinates of the grid center points of the ablation surface after retreat are obtained through the RBF dynamic meshing program for ablation surface retreat. The Euclidean distance between the interpolation point and each grid center point of the ablation surface after retreat is calculated. The grid center point of the ablation surface after retreat with the minimum Euclidean distance to the current interpolation point is found. The ablation information of the current interpolation point is assigned to the grid center point of the ablation surface after retreat, thus completing the wall mapping of the interpolated ablation information and outputting it as the ablation_BC.dat file.

[0017] The beneficial effects of this invention are as follows: (1) This method can realize the interpolation processing of the ablation site information scattered on the ablation surface during the ablation retreat of complex aircraft through four processes: ablation site layout, ablation site coordinate calculation after ablation retreat, ablation information interpolation, and interpolation information wall mapping. This method can obtain the continuous distribution of material ablation response on the entire ablation surface, further expand the application of one-dimensional material ablation response program in ablation calculation of complex-shaped aircraft, and also provide a new coupling scheme for the aerodynamic / thermal environment-material ablation coupled calculation method of aircraft.

[0018] (2) This method is suitable for interpolation processing of ablation site data in full-model or half-model aerodynamic / thermal-material ablation coupling calculation of various two-dimensional or three-dimensional complex aircraft configurations; it is suitable for interpolation of ablation site data when there is no ablation retreat or small ablation retreat on the aircraft surface during short-term ablation, and it is also suitable for interpolation of ablation site data when there is large ablation retreat on the aircraft surface during long-term severe ablation.

[0019] (3) This method is suitable for interpolation processing of ablation site data for complex-shaped aircraft, including laminar ablation, transition ablation and turbulent ablation. At the same time, this method has no restrictions on the ablation thermal protection materials used in the aircraft and can be adapted to the interpolation processing of ablation response information of various thermal protection materials such as carbon-based materials, silicon-based materials, resin materials and ceramic-based materials.

[0020] (4) This method can handle interpolation processing of various ablation information, including ablation site coordinates (x,y or x,y,z), wall mass flow rate, wall ablation temperature, wall ejection velocity, mass fraction of typical fiber thermochemical ablation components (taking carbon fiber as an example: O, O2, N, N2, CO, CO2, CN, C2N, C, C2, C3, C4, C5, NO) and mass fraction of typical matrix pyrolysis components (taking phenolic resin as an example: H2O, C6H6, C6H5OH, CH4, H2).

[0021] (5) In the process of interpolating the ablation information of the ablation sites after ablation, this method performs one-dimensional linear interpolation between adjacent ablation sites within a group and between corresponding ablation sites in two adjacent groups, both along the arrangement direction of the ablation sites within the group and perpendicular to the arrangement direction of the ablation sites within the group. By combining the grouping and layout scheme of the ablation sites and the two successive one-dimensional linear interpolations of the ablation site information in the two directions and the wall mapping of the interpolated information, this method can realize the interpolation processing of multiple ablation information in a finitely scattered ablation sites on a complex ablation surface. Compared with typical two-dimensional or three-dimensional interpolation algorithms, it has the advantages of simple algorithm, strong applicability and high computational efficiency. Attached Figure Description

[0022] Figure 1 This is a flowchart illustrating the method. Figure 2 A schematic diagram of the linear interpolation process for ablation information; Figure 3 The images show the initial blunt-nosed body configuration and typical transitional ablation morphology at different times; where (a) is the initial blunt-nosed body configuration, (b) is the typical transitional ablation morphology at time t1, and (c) is the typical transitional ablation morphology at time t2. Figure 4 Generatrixes representing the initial blunt-nosed body configuration and typical transitional ablation morphologies at different times; Figure 5 The results of heat flux calculation and interpolation for the initial blunt-nosed body configuration are compared; (a) is the calculated wall heat flux cloud map, (b) is the interpolated wall heat flux cloud map, (c) is the calculated wall heat flux data along the centerline, and (d) is the interpolated wall heat flux data along the centerline. Figure 6 The results show the convergence of the aerodynamic / thermal environment-material ablation coupling calculation for typical transitional ablation morphology at different times, where (a) is the data at time t1 and (b) is the data at time t2. Detailed Implementation

[0023] 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.

[0024] like Figure 1 As shown, the method for processing aircraft wall ablation response data based on interpolation algorithms includes the following steps: S1. Deploy ablation sites and output the initial coordinate information of the ablation sites; S2. Obtain the coordinate information of the grid center point of the ablation wall surface, and calculate the coordinates of the ablation site after the ablation retreat by combining the initial coordinate information of the ablation site. S3. Obtain ablation information at the ablation site, and interpolate the ablation information based on the coordinates of the ablation site after ablation retreat to obtain the interpolated ablation information. S4. Perform wall mapping on the interpolated ablation information to obtain the continuous distribution of the material ablation response across the entire ablation wall.

[0025] During the placement of ablation sites, groups of ablation sites are set up along the maximum length direction of the ablation wall of the aircraft (e.g., the x-axis direction). The coordinates of each group of ablation sites remain constant along the maximum length direction of the ablation wall and are evenly distributed circumferentially. The number of ablation sites in each group remains consistent. If the ablation wall covers the apex of the aircraft's nose, only one ablation site needs to be placed at the apex. In areas of intense ablation, such as near the apex of a blunt-nosed body, the distance between adjacent groups of ablation sites is appropriately reduced to better capture the rapidly changing aerodynamic / thermal environment and thus better characterize the changes in the corresponding material ablation response. The spacing between ablation sites within different groups remains relatively consistent. However, in areas of intense ablation, such as the windward side of the aircraft, the number of ablation sites can be appropriately increased, and the distance between adjacent ablation sites can be reduced.

[0026] In this embodiment, the specific method for setting up ablation sites in step S1 includes: like Figure 3 and Figure 4As shown, for the initial blunt-nosed body configuration and typical transitional ablation morphologies at different times, six groups of ablation sites were arranged along the x-axis direction for the initial blunt-nosed body configuration and the ablation morphologies at different times during typical transitional ablation using mesh drawing tools such as pointwise. The first group had one ablation site, located at the vertex of the blunt-nosed body; the remaining groups each had six ablation sites, and the x-axis coordinates of the ablation sites within each group were kept consistent. After the points were arranged, the coordinate information of the ablation sites was output in groups, generating the ablation_position.dat file.

[0027] In this embodiment, the specific method for obtaining the coordinate information of the grid center point of the ablated wall mesh in step S2 includes: reading the grid center point position coordinate information ablation_position.dat file generated by mesh drawing software such as pointwise, and reading the grid center point coordinate and normal vector information cellcenter_message.dat file given by the aircraft aerodynamic / thermal environment calculation CFD program to obtain the grid center point coordinate and normal vector of the ablated wall mesh.

[0028] In this embodiment, the specific method for calculating the coordinates of the ablation site after retreat in step S2, based on the initial coordinate information of the ablation site, includes the following steps: S2-1. Based on the principle of minimizing the Euclidean distance between the grid center point and the ablation site on the ablation wall, and using the formula:

[0029] Calculate the Euclidean distance between the current ablation site and the center points of all ablation wall grids. ;in The coordinates of the center point of any ablation wall grid; The coordinates of the current ablation site; S2-2. Based on the Euclidean distance between the current ablation site and the center points of all ablation wall grids, determine the center point of the ablation wall grid that is closest to the current ablation site, and use the normal vector information at the center point of the ablation wall grid as the normal vector at the current ablation site. S2-3. Based on the normal vector at the current ablation site, calculate the coordinates of the ablation site after the ablation retreat. The expression is as follows:

[0030] in The coordinates of the ablation site after the ablation retreat; The initial coordinates of the current ablation site; This is the initial normal vector for the current ablation site; This represents the ablation retreat at the current ablation site. S2-4. Output the coordinates of the ablation site after the ablation retreat as a file named new_ablation_position.dat.

[0031] In step S3, the specific method for obtaining ablation information at the ablation site is as follows: the ablation information at the ablation site is obtained by reading the ablation information file ablation_calculation.dat output by the material ablation analysis software.

[0032] In this embodiment, as Figure 2 As shown, the specific method for interpolating the ablation information in step S3 includes: interpolating the two-dimensional curve ablation information only in the x direction; and interpolating the three-dimensional surface ablation information in the z and x directions in turn.

[0033] In the specific implementation process, the method for interpolating the two-dimensional curve ablation information only in the x-direction includes the following steps: A1. Determine the projection plane and project all the coordinates of the ablation sites after the ablation retreats onto the xz plane to ensure that the ablation sites are not overlapping in the projection plane; assign the ablation information at the ablation site to the corresponding ablation site after the ablation retreats. A2. Determine the number of interpolation points between two adjacent ablation retreat sites within the same group and perform linear interpolation of the ablation information; that is, according to the formula... or A linear interpolation of ablation information is performed sequentially along the x-direction between two adjacent ablation retreat sites within the group; where... This is the result of a linear interpolation. and These are the ablation information of two adjacent ablation retreat sites within the group; and It is a linear interpolation factor; A3. Output the interpolation results of the ablation sites in the projection plane after ablation and generate the interpolation_ablation_message.dat file to obtain the ablation information of each interpolation point, that is, the ablation information after two-dimensional curve interpolation.

[0034] In the specific implementation process, the method for interpolating the three-dimensional surface ablation information sequentially in the z and x directions includes the following steps: B1. Determine the projection plane and project all the coordinates of the ablation sites after the ablation retreat into the xz plane to ensure that the ablation sites are not overlapping in the projection plane; assign the ablation information at the ablation site to the corresponding ablation site after the ablation retreat. B2. Determine the number of groups, the number of ablation retreat sites in the same group, the number of interpolations between two adjacent ablation retreat sites in the x direction, and the number of interpolations between two adjacent ablation retreat sites in the z direction. B3. Copy and expand the ablation information of the ablation retreat and ablation site group at the end point of the aircraft (i.e. the apex of the blunt body) to be consistent with the number of other ablation retreat and ablation site groups (6). B4. Perform linear interpolation of the ablation information within each group for the ablation sites after retreat, i.e., according to the formula... or A linear interpolation of ablation information is performed sequentially along the z-direction between two adjacent ablation sites after retreat within each group; where... This is the result of a linear interpolation. and These are the ablation information of two adjacent ablation retreat sites within the group; and It is a linear interpolation factor; B5. Along the x-direction, sequentially perform secondary interpolation of the ablation information between the corresponding ablation sites in the two adjacent groups after linear interpolation, output the interpolation results of the ablation sites in the projection plane and generate the interpolation_ablation_message.dat file to obtain the ablation information of each interpolation point, that is, the ablation information after three-dimensional surface interpolation.

[0035] In this embodiment, , Pi , The coordinates of the current point to be interpolated along the interpolation direction; The coordinates of the first ablation retreat ablation site among two adjacent ablation retreat ablation sites; The coordinates are the coordinates of the second ablation retreat ablation site among two adjacent ablation retreat ablation sites.

[0036] As can be seen from the above process, in the ablation information interpolation process, it is first necessary to project the coordinates of the ablation sites after ablation onto the plane containing the maximum length direction of the ablation wall of the aircraft, ensuring that the ablation sites after ablation are not overlapping in the projection plane. When interpolating the ablation site information, firstly, intra-group interpolation is performed, that is, one-dimensional linear interpolation of the ablation information is performed sequentially along a certain direction between two adjacent ablation sites within a group. After the interpolation in this direction is completed, a second interpolation of the ablation information is performed sequentially along another direction perpendicular to it between corresponding sites in two adjacent groups after the first interpolation. Finally, the interpolation results of the ablation sites in the projection plane are output. It should be noted that if the aircraft calculation model is a full model, it is necessary to divide the aircraft into upper and lower parts along the projection plane, and sequentially project the ablation sites of each part onto the projection plane for interpolation calculation. When interpolating the ablation information between two adjacent ablation sites, different linear interpolation algorithms can be selected to better represent the distribution of ablation information.

[0037] In this embodiment, the wall mapping of the interpolated ablation information is performed last. The interpolated information of the ablation points in the projection plane is projected onto the ablation wall of the aircraft using the criterion of minimum Euclidean distance between the interpolation point and the center point of the ablation surface grid after ablation retreat. In this process, for each interpolation point, the center points of the ablation surface grid after ablation retreat are traversed sequentially, the Euclidean distance between the current interpolation point and each center point of the ablation surface grid after ablation retreat is calculated, and the center point of the ablation surface grid after ablation retreat with the current interpolation point is found. The ablation information of the current interpolation point is assigned to this grid point and output as ablation_BC.dat.

[0038] Figure 5 The figure shows a comparison of the calculation and interpolation results of the initial blunt-nosed body configuration heat flux data. As can be seen from the figure, this method can effectively achieve wall interpolation of aerodynamic / ablation response physical information such as heat flux of the ablation wall. It can obtain the continuous distribution of the ablation response of the wall material by interpolation processing based on the ablation information of a limited number of sites. Figure 6 The iterative convergence of the coupled aerodynamic / thermal environment-material ablation calculation for typical transition ablation morphologies at different time points reveals that this method can achieve coupled calculation of the aerodynamic / thermal environment-material ablation response of complex-shaped aircraft. Furthermore, after three rounds of coupled iterations, the surface heat flux of typical transition ablation morphologies at each time point achieves good convergence. This indicates that the method meets the requirements for engineering applications.

[0039] In practical implementation, this method can be applied to the coupled calculation of aerodynamic / thermal environment and material ablation response of complex-shaped aircraft. Complex-shaped aircraft include two-dimensional and three-dimensional configurations. The wall ablation response data includes, but is not limited to, typical carbon-carbon material surface thermochemical ablation data: ablation site coordinates (x,y or x,y,z), wall mass flow rate, wall ablation temperature, wall ejection velocity, and mass fraction of carbothermal ablation components (O, O2, N, N2, CO, CO2, CN, ...). Volumetric ablation data for C2N, C, C2, C3, C4, C5, NO) and typical carbon phenolic materials: ablation site coordinates (x, y or x, y, z), wall mass flow rate, wall ablation temperature, wall ejection velocity, mass fraction of carbothermal ablation components (O, O2, N, N2, CO, CO2, CN, C2N, C, C2, C3, C4, C5, NO) and mass fraction of phenolic pyrolysis components (H2O, C6H6, C6H5OH, CH4, H2).

Claims

1. A method for processing aircraft wall ablation response data based on interpolation algorithm, characterized in that, Includes the following steps: Determine the ablation sites and output the initial coordinate information of the ablation sites; Obtain the coordinate information of the grid center point of the ablation wall surface, and calculate the coordinates of the ablation site after the ablation retreat by combining the initial coordinate information of the ablation site. Obtain ablation information at the ablation site, and interpolate the ablation information based on the coordinates of the ablation site after ablation retreat to obtain the interpolated ablation information. The interpolated ablation information is mapped onto the wall surface to obtain the continuous distribution of the material ablation response across the entire ablated wall surface.

2. The method for processing aircraft wall ablation response data based on interpolation algorithm according to claim 1, characterized in that, Specific methods for setting up ablation sites include: For the initial blunt-nosed body configuration and typical transition ablation morphologies at different times, a mesh drawing tool was used to arrange six groups of ablation sites along the x-axis for the initial blunt-nosed body configuration and the ablation morphologies at different times during typical transition ablation. The first group had one ablation site, located at the apex of the blunt-nosed body; the remaining groups each had six ablation sites, and the x-axis coordinates of the ablation sites within each group were kept consistent. After the points were arranged, the coordinate information of the ablation sites was output in groups, generating the ablation_position.dat file.

3. The method for processing aircraft wall ablation response data based on interpolation algorithm according to claim 1, characterized in that, Specific methods for obtaining the coordinate information of the grid center points of the ablated wall surface include: Read the grid center point position coordinate information ablation_position.dat file generated by the mesh drawing software, and read the grid center point coordinates and normal vector information cellcenter_message.dat file given by the ablation wall mesh CFD program for calculating the ablation environment, to obtain the grid center point coordinates and normal vector of the ablation wall mesh.

4. The method for processing aircraft wall ablation response data based on interpolation algorithm according to claim 1, characterized in that, The specific method for calculating the coordinates of the ablation site after retreat, based on the initial coordinate information of the ablation site, includes the following steps: Based on the principle of minimizing the Euclidean distance between the grid center and the ablation site on the ablation wall, and using the formula: Calculate the Euclidean distance between the current ablation site and the center points of all ablation wall grids. ;in The coordinates of the center point of any ablation wall grid; The coordinates of the current ablation site; Based on the Euclidean distance between the current ablation site and the center points of all ablation wall grids, determine the center point of the ablation wall grid that is closest to the current ablation site, and use the normal vector information at the center point of the ablation wall grid as the initial normal vector of the current ablation site. Based on the initial normal vector of the current ablation site, the coordinates of the ablation site after retreat are calculated, and the expression is as follows: in The coordinates of the ablation site after the ablation retreat; The initial coordinates of the current ablation site; This is the initial normal vector for the current ablation site; This represents the ablation retreat at the current ablation site. Output the coordinates of the ablation site after the ablation retreat as a file named new_ablation_position.dat.

5. The method for processing aircraft wall ablation response data based on interpolation algorithm according to claim 1, characterized in that, The specific method for obtaining ablation information at the ablation site is as follows: Ablation information at the ablation site is obtained by reading the ablation information file ablation_calculation.dat output by the material ablation analysis software.

6. The method for processing aircraft wall ablation response data based on interpolation algorithm according to claim 1, characterized in that, Specific methods for interpolating ablation information include: Interpolation is performed only in the x-direction for the ablation information of the two-dimensional curve; interpolation is then performed in the z-direction and x-direction for the ablation information of the three-dimensional surface.

7. The method for processing aircraft wall ablation response data based on interpolation algorithm according to claim 6, characterized in that, The specific method for interpolating the ablation information of a two-dimensional curve only in the x-direction includes the following steps: Determine the projection plane and project all the coordinates of the ablation sites after ablation retreat onto the xz plane; assign the ablation information at the ablation site to the corresponding ablation site after ablation retreat. Determine the number of interpolation points between two adjacent ablation retreat sites within the same group and perform linear interpolation of the ablation information; that is, according to the formula or A linear interpolation of ablation information is performed sequentially along the x-direction between two adjacent ablation retreat sites within the group; where... This is the result of a linear interpolation. and These are the ablation information of two adjacent ablation retreat sites within the group; and It is a linear interpolation factor; Output the interpolation results of the ablation sites in the projection plane after ablation and generate the interpolation_ablation_message.dat file to obtain the ablation information of each interpolation point, that is, the ablation information after two-dimensional curve interpolation.

8. The method for processing aircraft wall ablation response data based on interpolation algorithm according to claim 6, characterized in that, The specific method for interpolating the ablation information of a three-dimensional surface in the z and x directions sequentially includes the following steps: Determine the projection plane and project all the coordinates of the ablation sites after ablation retreat onto the xz plane; assign the ablation information at the ablation site to the corresponding ablation site after ablation retreat. Determine the number of groups, the number of ablation retreat sites in the same group, the number of interpolations between two adjacent ablation retreat sites in the x-direction, and the number of interpolations between two adjacent ablation retreat sites in the z-direction. The ablation information of the ablation retreat and ablation site group at the end of the aircraft is copied and expanded to be consistent with the number of other ablation retreat and ablation site groups. For each group of ablation sites after retreat, a linear interpolation of the ablation information within the group is performed, i.e., according to the formula... or A linear interpolation of ablation information is performed sequentially along the z-direction between two adjacent ablation sites after retreat within each group; where... This is the result of a linear interpolation. and These are the ablation information of two adjacent ablation retreat sites within the group; and It is a linear interpolation factor; Along the x-direction, after a linear interpolation, a second interpolation of the ablation information is performed between the corresponding ablation sites in the adjacent two groups of ablation sites. The interpolation results of the ablation sites in the projection plane are output and the interpolation_ablation_message.dat file is generated to obtain the ablation information of each interpolation point, that is, the ablation information after three-dimensional surface interpolation.

9. The method for processing aircraft wall ablation response data based on interpolation algorithm according to claim 7 or 8, characterized in that, , Pi , The coordinates of the current point to be interpolated along the interpolation direction; The coordinates of the first ablation retreat ablation site among two adjacent ablation retreat ablation sites; The coordinates are the coordinates of the second ablation retreat ablation site among two adjacent ablation retreat ablation sites.

10. The method for processing aircraft wall ablation response data based on interpolation algorithm according to claim 1, characterized in that, Specific methods for wall mapping of interpolated ablation information include: The coordinates of the grid center points of the ablation surface after retreat are obtained through the RBF dynamic meshing program for ablation surface retreat. The Euclidean distance between the interpolation point and each grid center point of the ablation surface after retreat is calculated. The grid center point of the ablation surface after retreat with the minimum Euclidean distance to the current interpolation point is found. The ablation information of the current interpolation point is assigned to the grid center point of the ablation surface after retreat, thus completing the wall mapping of the interpolated ablation information and outputting it as the ablation_BC.dat file.

Images