Method, device and medium for calculating frictional heat flow of turbulent boundary layer of ablation particle wall

Abstract

The application relates to the technical field of aircrafts, and discloses an ablative particle wall turbulent boundary layer frictional resistance and heat flow calculation method, equipment and medium, the method comprising the following steps: acquiring particle two-phase wall turbulent direct numerical simulation data; calculating each sub-effect corresponding to frictional resistance decomposition and each sub-effect corresponding to heat flow decomposition respectively according to the particle two-phase wall turbulent direct numerical simulation data; performing MPI reduction operation summation on each sub-effect corresponding to the frictional resistance decomposition to obtain the frictional resistance of the ablative particle wall turbulent boundary layer; and performing MPI reduction operation summation on each sub-effect corresponding to the heat flow decomposition to obtain the heat flow of the ablative particle wall turbulent boundary layer. The above method has strong compatibility and can be coupled with any compressible particle two-phase flow calculation mode; the method is efficient in parallel and can quickly process particle two-phase wall turbulent direct numerical simulation data to obtain the frictional resistance and heat flow of the ablative particle two-phase flow boundary layer, thereby providing a technical means for diagnosing and optimizing the design of high-speed ablative aircrafts.

Description

Technical Field

[0001] The present invention relates to the field of aircraft technology, and in particular to a method, equipment and medium for calculating the friction heat flow of a turbulent boundary layer of an ablative particle wall. Background Art

[0002] High-speed aircraft are powerful weapons for a wide range of applications. During reentry, the stagnation point at the vehicle's nose is heated to high temperatures, generating ablated particles. These particles, driven by aerodynamic forces, enter the downstream boundary layer. These ablated particles interact with the compressible boundary layer, altering its dynamic and thermodynamic processes.

[0003] In existing technology, full detailed information about the particle-fluid two-phase system is typically obtained through numerical wind tunnel measurements. However, engineering organizations are interested in the frictional drag and heat flux generation mechanisms of aircraft, but there is currently a lack of efficient and rapid analysis of massive amounts of direct numerical simulation wind tunnel measurement data to understand these mechanisms. Summary of the Invention

[0004] In view of this, the present invention aims to provide a method, device, and medium for calculating the friction and heat flux of the turbulent boundary layer of ablative particle walls. This method can efficiently and concurrently obtain the friction and heat flux of the boundary layer of ablative particle two-phase flow, with strong compatibility. The specific scheme is as follows:

[0005] A method for calculating frictional heat flow in a turbulent boundary layer of an ablative particle wall includes:

[0006] Obtain direct numerical simulation data of granular two-phase wall turbulence;

[0007] According to the direct numerical simulation data of the particle two-phase wall turbulence, each sub-effect corresponding to the friction decomposition and each sub-effect corresponding to the heat flow decomposition are calculated respectively;

[0008] Performing an MPI reduction operation on each sub-effect corresponding to the friction decomposition and summing them up to obtain the friction of the turbulent boundary layer of the ablated particle wall;

[0009] The MPI reduction operation is performed on each sub-effect corresponding to the heat flow decomposition and the sum is obtained to obtain the heat flow of the turbulent boundary layer of the ablated particle wall.

[0010] Preferably, in the above-mentioned method for calculating the frictional heat flux of the ablative particle wall turbulent boundary layer provided by the embodiment of the present invention, the calculation of each sub-effect corresponding to the frictional decomposition is performed based on the direct numerical simulation data of the particle two-phase wall turbulence, including:

[0011] According to the direct numerical simulation data of the particle two-phase wall turbulence, the particle two-phase viscosity effect, compressibility effect, turbulence pulse effect, wake effect, pressure gradient effect, particle average size effect and particle pulse size effect are calculated.

[0012] Preferably, in the above-mentioned method for calculating the heat flux of the ablative particle wall turbulent boundary layer provided in an embodiment of the present invention, the calculation of each sub-effect corresponding to the heat flux decomposition is performed based on the direct numerical simulation data of the particle two-phase wall turbulence, including:

[0013] Based on the direct numerical simulation data of the particle two-phase wall turbulence, the normal enthalpy flow and normal pressure energy flow effects, the friction effect, the turbulence pulsation effect, and the particle heat transfer effect are calculated.

[0014] Preferably, in the above-mentioned ablation particle wall turbulent boundary layer friction heat flow calculation method provided in an embodiment of the present invention, according to the direct numerical simulation data of the particle two-phase wall turbulence, the particle two-phase viscosity effect, compressibility effect, turbulence pulse effect, wake effect, pressure gradient effect, particle average size effect and particle pulse size effect are calculated, including:

[0015] Calculate the viscous force effect based on the boundary layer thickness and fluid dynamic viscosity;

[0016] Calculating a compressibility effect based on the boundary layer thickness, the fluid dynamic viscosity, the average velocity of the fluid flow direction, and the fluid velocity divergence;

[0017] Calculating turbulent pulse effects based on the boundary layer thickness, the fluid dynamic viscosity, the fluid density, and the Reynolds stress;

[0018] Calculating a wake effect based on the boundary layer thickness, the fluid density, and the fluid kinetic energy;

[0019] Calculate the pressure gradient effect based on fluid flow and fluid pressure;

[0020] Calculating the average particle size effect based on the boundary layer thickness, the average flow velocity of the fluid, the apparent density of the particle group, the particle inertia relaxation coefficient, and the average flow velocity of the particles;

[0021] The particle pulse scale effect is calculated according to the boundary layer thickness, the particle flow pulsating velocity, the fluid flow pulsating velocity, the particle group apparent density, and the particle inertial relaxation coefficient.

[0022] Preferably, in the above-mentioned method for calculating the friction heat flow of the turbulent boundary layer of the ablation particle wall provided in an embodiment of the present invention, the calculation of the normal enthalpy flow and normal pressure energy flow effects, the friction work effect, the turbulent pulsation effect, and the particle heat transfer effect is performed based on the direct numerical simulation data of the particle two-phase wall turbulence, including:

[0023] Calculating normal enthalpy flow and normal pressure energy flow effects based on the fluid enthalpy, the fluid normal average velocity, and the fluid pressure;

[0024] Calculating the effect of friction based on the boundary layer thickness, the average velocity of the fluid flow, and the fluid shear stress;

[0025] Calculating turbulence pulsation effects based on the boundary layer thickness, the fluid specific heat capacity, the fluid pressure, and normal turbulence temperature pulsation;

[0026] The particle heat transfer effect is calculated based on the average temperature of the particles, the average temperature of the fluid, the apparent density of the particle group, the specific heat capacity of the particles, and the thermal response coefficient of the particles.

[0027] Preferably, in the above-mentioned method for calculating the frictional heat flux of the turbulent boundary layer of the ablation particle wall provided in the embodiment of the present invention, the frictional heat flux of the turbulent boundary layer of the ablation particle wall is obtained by the following formula:

[0028]

[0029] Among them, the overline indicates the statistical average, and the superscript ' indicates the pulsating amount; represents the viscous force effect, represents the compressibility effect, represents the turbulent pulsation effect, represents the wake effect, represents the pressure gradient effect, represents the average particle size effect, represents the particle pulsation scale effect; δ represents the boundary layer thickness, 0≤y≤δ, represents the fluid dynamic viscosity, represents the fluid density, represents the average velocity of the fluid flow, u' represents the pulsating velocity of the fluid flow, θ represents the fluid velocity divergence, represents the Reynolds stress, D represents the body derivative, t represents the time variable, represents the kinetic energy of the fluid, m represents the fluid flow rate, represents the fluid pressure, x represents the displacement variable, ρ s represents the apparent density of the particle group, τ p represents the particle inertia relaxation coefficient, represents the average particle flow velocity, and V' represents the particle flow pulsating velocity.

[0030] Preferably, in the above-mentioned method for calculating the friction heat flux of the turbulent boundary layer of the ablation particle wall provided in the embodiment of the present invention, the heat flux of the turbulent boundary layer of the ablation particle wall is obtained by using the following formula:

[0031]

[0032] in, represents the effects of normal enthalpy flow and normal pressure flow, Indicates that friction is the effect of the effect. represents the turbulent pulsation effect, Represents particles represents the fluid enthalpy, represents the average normal velocity of the fluid, represents the fluid shear stress, C p represents the specific heat capacity of the fluid, represents the normal turbulent temperature fluctuation, represents the average particle temperature, represents the average temperature of the fluid, C p,p represents the particle specific heat capacity, τ p,convec Represents the particle thermal response coefficient.

[0033] An embodiment of the present invention further provides an electronic device comprising a processor and a memory, wherein when the processor executes the computer program stored in the memory, the method for calculating the frictional heat flow of the turbulent boundary layer of the ablated particle wall provided in the embodiment of the present invention is implemented.

[0034] An embodiment of the present invention further provides a computer-readable storage medium for storing a computer program, wherein when the computer program is executed by a processor, the method for calculating the frictional heat flow of the turbulent boundary layer of the ablated particle wall provided in the embodiment of the present invention is implemented.

[0035] It can be seen from the above technical solution that the present invention provides a method for calculating the friction and heat flux of the turbulent boundary layer of the ablative particle wall, including: obtaining direct numerical simulation data of granular two-phase wall turbulence; calculating the sub-effects corresponding to the friction decomposition and the sub-effects corresponding to the heat flux decomposition according to the direct numerical simulation data of the granular two-phase wall turbulence; performing an MPI reduction operation on the sub-effects corresponding to the friction decomposition to obtain the friction of the turbulent boundary layer of the ablative particle wall; performing an MPI reduction operation on the sub-effects corresponding to the heat flux decomposition to obtain the heat flux of the turbulent boundary layer of the ablative particle wall.

[0036] The above-mentioned calculation method for the friction and heat flow of the ablative particle wall turbulent boundary layer provided by the present invention has strong compatibility and can be coupled with any compressible particle two-phase flow calculation method; it is efficient and parallel, and can quickly process the direct numerical simulation data of particle two-phase wall turbulence to obtain the friction and heat flow of the ablative particle two-phase flow boundary layer, providing a technical means for the diagnosis and optimization design of high-speed ablative aircraft, filling the technical gap in this field.

[0037] In addition, the present invention also provides corresponding equipment and computer-readable storage media for the calculation method of the friction heat flow of the turbulent boundary layer of the ablation particle wall, further making the above method more practical. The equipment and computer-readable storage medium have corresponding advantages. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or related technologies, the following briefly introduces the drawings required for use in the embodiments or related technical descriptions. Obviously, the drawings described below are merely embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0039] Figure 1 A flow chart of a method for calculating frictional heat flux in a turbulent boundary layer of an ablated particle wall provided in an embodiment of the present invention;

[0040] Figure 2 A schematic diagram of the process corresponding to the method for calculating the frictional heat flow of the turbulent boundary layer of the ablation particle wall provided by an embodiment of the present invention;

[0041] Figure 3 A schematic diagram of the calculation process of the friction solver provided in an embodiment of the present invention;

[0042] Figure 4 A schematic diagram of the calculation process of the heat flow solver provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0043] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0044] The present invention provides a method for calculating the friction heat flow of the turbulent boundary layer of the ablation particle wall. Figure 1 As shown, the following steps are included:

[0045] S101, obtaining direct numerical simulation data of particle two-phase wall turbulence;

[0046] Specifically, if Figure 2 As shown, any compressible ablative particle two-phase wall turbulence solver can be used to generate information such as the density, velocity, and temperature of the particle-fluid two-phase in parallel calculation to obtain direct numerical simulation data of the particle two-phase wall turbulence.

[0047] S102. Calculate the sub-effects corresponding to the friction decomposition and the sub-effects corresponding to the heat flow decomposition based on the direct numerical simulation data of the particle two-phase wall turbulence;

[0048] S103, performing an MPI reduction operation on each sub-effect corresponding to the friction decomposition and summing them up to obtain the friction of the turbulent boundary layer of the ablated particle wall;

[0049] Specifically, if Figure 2 As shown in the figure, the compressible particle wall turbulence friction solver can be used to read the direct numerical simulation data of particle two-phase wall turbulence in the form of parallel computing, calculate the various sub-effects corresponding to the friction decomposition, and perform MPI reduction operations on the various sub-effects corresponding to the friction decomposition to sum up and obtain the friction of the ablative particle wall turbulence boundary layer.

[0050] S104 , performing an MPI reduction operation on each sub-effect corresponding to the heat flow decomposition and summing the sum to obtain the heat flow of the turbulent boundary layer of the ablated particle wall.

[0051] Specifically, if Figure 2 As shown in the figure, the compressible particle wall turbulence heat flow solver can be used to read the direct numerical simulation data of particle two-phase wall turbulence in the form of parallel computing, calculate the various sub-effects corresponding to the heat flow decomposition, and perform MPI reduction operations on the various sub-effects corresponding to the heat flow decomposition to obtain the heat flow of the ablative particle wall turbulence boundary layer.

[0052] The above-mentioned ablative particle wall turbulence boundary layer friction and heat flow calculation method provided in the embodiment of the present invention has strong compatibility and can be coupled with any compressible particle two-phase flow calculation method; it is efficient and parallel, and can quickly process the direct numerical simulation data of particle two-phase wall turbulence to obtain the friction and heat flow of the ablative particle two-phase flow boundary layer, providing a technical means for the diagnosis and optimization design of high-speed ablative aircraft, filling the technical gap in this field.

[0053] Furthermore, in the specific implementation, in the above-mentioned ablation particle wall turbulence boundary layer friction heat flow calculation method provided in the embodiment of the present invention, step S102 calculates each sub-effect corresponding to the friction decomposition based on the direct numerical simulation data of the particle two-phase wall turbulence, such as Figure 3 As shown, it can specifically include: calculating the particle two-phase viscosity effect, compressibility effect, turbulence pulse effect, wake effect, pressure gradient effect, particle average scale effect and particle pulse scale effect based on the direct numerical simulation data of particle two-phase wall turbulence.

[0054] In specific implementation, the above steps may include calculating the particle two-phase viscosity effect, compressibility effect, turbulence pulse effect, wake effect, pressure gradient effect, particle average size effect, and particle pulse size effect based on the direct numerical simulation data of particle two-phase wall turbulence, which may specifically include:

[0055] Calculate the viscous force effect based on the boundary layer thickness and fluid dynamic viscosity;

[0056] Calculate compressibility effects based on boundary layer thickness, fluid dynamic viscosity, fluid streamwise average velocity, and fluid velocity divergence;

[0057] Calculate turbulent pulse effects based on boundary layer thickness, fluid dynamic viscosity, fluid density, and Reynolds stress;

[0058] Calculate the wake effect based on the boundary layer thickness, fluid density, and fluid kinetic energy;

[0059] Calculate the pressure gradient effect based on fluid flow and fluid pressure;

[0060] The average particle size effect is calculated based on the boundary layer thickness, the average fluid flow velocity, the apparent density of the particle group, the particle inertia relaxation coefficient, and the average particle flow velocity;

[0061] The particle pulse scale effect is calculated based on the boundary layer thickness, particle flow pulsation velocity, fluid flow pulsation velocity, particle group apparent density, and particle inertia relaxation coefficient.

[0062] Specifically, the compressible particle wall turbulent friction solver reads massive data through parallel computing, such as boundary layer thickness, fluid dynamic viscosity, average fluid flow velocity, fluid velocity divergence, fluid density, Reynolds stress, fluid kinetic energy, fluid flow rate, fluid pressure, particle group apparent density, particle inertial relaxation coefficient, particle average flow velocity, particle flow pulsating velocity, fluid flow pulsating velocity and other data, and calculates each sub-effect corresponding to the friction decomposition one by one: viscosity, compressibility, turbulent pulsation, wake, pressure gradient, particle effect, and then performs standard MPI reduction operations to sum up and obtain the total friction and the proportion of each component.

[0063] Furthermore, in the specific implementation, in the above-mentioned ablation particle wall turbulence boundary layer friction heat flux calculation method provided in the embodiment of the present invention, step S102 calculates each sub-effect corresponding to the heat flux decomposition based on the direct numerical simulation data of the particle two-phase wall turbulence, such as Figure 4 As shown, it can specifically include: calculating the normal enthalpy flow and normal pressure energy flow effects, friction work effect, turbulence pulsation effect, and particle heat transfer effect based on direct numerical simulation data of particle two-phase wall turbulence.

[0064] In specific implementation, in the above-mentioned ablation particle wall turbulent boundary layer friction heat flow calculation method provided in the embodiment of the present invention, based on the direct numerical simulation data of the particle two-phase wall turbulence, the normal enthalpy flow and normal pressure energy flow effects, the friction work effect, the turbulent pulsation effect, and the particle heat transfer effect are calculated, which may specifically include:

[0065] Calculate the normal enthalpy flow and normal pressure energy flow effects based on the fluid enthalpy, fluid normal average velocity, and fluid pressure;

[0066] Calculate the friction effect based on the boundary layer thickness, average velocity of the fluid flow, and fluid shear stress;

[0067] Calculate turbulent fluctuation effects based on boundary layer thickness, fluid specific heat capacity, fluid pressure, and normal turbulent temperature fluctuations;

[0068] The particle heat transfer effect is calculated based on the average particle temperature, the average fluid temperature, the apparent density of the particle group, the particle specific heat capacity, and the particle thermal response coefficient.

[0069] Specifically, the compressible particle wall turbulent heat flow solver reads massive data through parallel computing, such as fluid enthalpy, fluid normal average velocity, fluid pressure, boundary layer thickness, fluid flow average velocity, fluid shear stress, fluid specific heat capacity, normal turbulent temperature pulsation, particle average temperature, fluid average temperature, particle group apparent density, particle specific heat capacity, particle thermal response coefficient and other data, and calculates each sub-effect corresponding to the heat flow decomposition one by one: friction work, normal enthalpy flow, normal pressure energy flow, turbulent temperature and velocity pulsation, particle effect, and then performs standard MPI reduction operations to sum up and obtain the total friction and the proportion of each component.

[0070] In the specific implementation, in the above-mentioned calculation method of the friction heat flux of the turbulent boundary layer of the ablation particle wall provided in the embodiment of the present invention, the friction heat flux of the turbulent boundary layer of the ablation particle wall is obtained by the following formula: w :

[0071]

[0072] For the sake of simplicity, all the superscripts - in the formulas of this application represent statistical averages, and the superscript ' represents pulsation quantities; the first term on the right side of the equal sign in formula (1) is represents the viscous force effect, the second term represents the compressibility effect, and the third term Indicates the turbulent pulsation effect, the fourth represents the wake effect, the fifth term represents the pressure gradient effect, the sixth term represents the average particle size effect, the seventh term represents the particle pulsation scale effect; δ represents the boundary layer thickness, 0≤y≤δ, represents the fluid dynamic viscosity, represents the fluid density, represents the average velocity of the fluid flow, u' represents the pulsating velocity of the fluid flow, θ represents the fluid velocity divergence, represents the Reynolds stress, D represents the body derivative, t represents the time variable, represents the kinetic energy of the fluid, m represents the fluid flow rate, represents the fluid pressure, x represents the displacement variable, ρ s represents the apparent density of the particle group, τ p represents the particle inertia relaxation coefficient, represents the average particle flow velocity, and V' represents the particle flow pulsating velocity.

[0073] In the specific implementation, in the above-mentioned calculation method of friction heat flux of turbulent boundary layer of ablation particle wall provided in the embodiment of the present invention, the heat flux q of turbulent boundary layer of ablation particle wall is obtained by the following formula: w :

[0074]

[0075] Among them, the first term on the right side of the equation (2) is represents the effects of normal enthalpy flow and normal pressure flow, and the second term The third term represents the effect of friction. represents the turbulent pulsation effect, the fourth term represents the particle heat transfer effect; represents the fluid enthalpy, represents the average normal velocity of the fluid, represents the fluid shear stress, C p represents the specific heat capacity of the fluid, represents the normal turbulent temperature fluctuation, represents the average particle temperature, represents the average temperature of the fluid, C p,p represents the particle specific heat capacity, τ p,convec Represents the particle thermal response coefficient.

[0076] It should be noted that the error between the calculation method for the frictional heat flux of the turbulent boundary layer of the ablation particle wall provided by the embodiment of the present invention and the theoretical solution is on the order of 0.1%.

[0077] Correspondingly, an embodiment of the present invention further discloses an electronic device comprising a processor and a memory; wherein, when the processor executes the computer program stored in the memory, the method for calculating the frictional heat flow of the turbulent boundary layer of the ablated particle wall disclosed in the aforementioned embodiment is implemented.

[0078] For more specific details about the above method, please refer to the corresponding contents disclosed in the aforementioned embodiments, which will not be described again here.

[0079] Furthermore, the present invention also discloses a computer-readable storage medium for storing a computer program; when the computer program is executed by a processor, the aforementioned method for calculating the frictional heat flow of the turbulent boundary layer of the ablated particle wall is implemented.

[0080] For more specific details about the above method, please refer to the corresponding contents disclosed in the aforementioned embodiments, which will not be described again here.

[0081] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from the other embodiments. References to the same or similar portions of the various embodiments will be sufficient. The devices and storage media disclosed in the embodiments are described briefly because they correspond to the methods disclosed in the embodiments. For relevant details, refer to the method descriptions.

[0082] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the components and steps of each example according to their functions. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians may use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0083] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in a random access memory (RAM), internal memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0084] In summary, an embodiment of the present invention provides a method for calculating the friction and heat flux of an ablative particle wall turbulent boundary layer, comprising: obtaining direct numerical simulation data of granular two-phase wall turbulence; calculating each sub-effect corresponding to the friction decomposition and each sub-effect corresponding to the heat flux decomposition based on the direct numerical simulation data of the granular two-phase wall turbulence; performing an MPI reduction operation on each sub-effect corresponding to the friction decomposition and summing them to obtain the friction of the ablative particle wall turbulent boundary layer; performing an MPI reduction operation on each sub-effect corresponding to the heat flux decomposition and summing them to obtain the heat flux of the ablative particle wall turbulent boundary layer. The above-mentioned method for calculating the friction and heat flux of an ablative particle wall turbulent boundary layer has strong compatibility and can be coupled with any compressible particle two-phase flow calculation method; it is efficient and parallel, and can quickly process the direct numerical simulation data of granular two-phase wall turbulence to obtain the friction and heat flux of the ablative particle two-phase flow boundary layer, providing a technical means for diagnosing and optimizing the design of high-speed ablative aircraft, filling the technical gap in this field. In addition, the present invention also provides corresponding equipment and computer-readable storage media for the calculation method of the friction heat flow of the turbulent boundary layer of the ablation particle wall, further making the above method more practical. The equipment and computer-readable storage medium have corresponding advantages.

[0085] Finally, it should be noted that, in this document, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or device. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article, or device comprising the element.

[0086] The above is a detailed introduction to the method, equipment and medium for calculating the friction heat flux of the turbulent boundary layer of the ablation particle wall provided by the present invention. Specific examples are used in this article to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea; at the same time, for general technical personnel in this field, according to the ideas of the present invention, there will be changes in the specific implementation methods and application scopes. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for calculating the heat flux of the turbulent boundary layer of ablated particle wall, characterized in that: include: Obtain direct numerical simulation data of granular two-phase wall turbulence; According to the direct numerical simulation data of the particle two-phase wall turbulence, each sub-effect corresponding to the friction decomposition and each sub-effect corresponding to the heat flow decomposition are calculated respectively; The MPI reduction operation is performed on each sub-effect corresponding to the friction decomposition to obtain the friction of the turbulent boundary layer of the ablation particle wall using the following formula: ； Among them, the overline Indicates statistical average, superscript Indicates the pulsation amount; represents the viscous force effect, represents the compressibility effect, represents the turbulent pulsation effect, represents the wake effect, represents the pressure gradient effect, represents the average particle size effect, represents the particle pulsation scale effect; represents the boundary layer thickness, , represents the fluid dynamic viscosity, represents the fluid density, represents the average velocity of the fluid flow, Indicates the fluid flow pulsation velocity, represents the fluid velocity divergence, represents the Reynolds stress, D represents the body derivative, t represents the time variable, represents the kinetic energy of the fluid, represents the fluid flow rate, represents the fluid pressure, x represents the displacement variable, represents the apparent density of the particle group, represents the particle inertia relaxation coefficient, represents the average streamwise velocity of the particles, Indicates the particle flow pulsation velocity; The MPI reduction operation is performed on each sub-effect corresponding to the heat flow decomposition, and the heat flow of the turbulent boundary layer of the ablation particle wall is obtained using the following formula: ； in, represents the effects of normal enthalpy flow and normal pressure flow, Indicates that friction is the effect of force. represents the turbulent pulsation effect, represents the particle heat transfer effect; represents the fluid enthalpy, represents the average normal velocity of the fluid, represents the fluid shear stress, represents the specific heat capacity of the fluid, represents the normal turbulent temperature fluctuation, represents the average particle temperature, represents the average temperature of the fluid, represents the specific heat capacity of the particles, Represents the particle thermal response coefficient.

2. The method for calculating the heat flux of the turbulent boundary layer of the ablated particle wall according to claim 1, characterized in that: Based on the direct numerical simulation data of the particle two-phase wall turbulence, the sub-effects corresponding to the friction decomposition are calculated, including: According to the direct numerical simulation data of the particle two-phase wall turbulence, the particle two-phase viscosity effect, compressibility effect, turbulence pulse effect, wake effect, pressure gradient effect, particle average size effect and particle pulse size effect are calculated.

3. The method for calculating the heat flux of the turbulent boundary layer of the ablated particle wall according to claim 2, characterized in that: Based on the direct numerical simulation data of the particle two-phase wall turbulence, the various sub-effects corresponding to the heat flux decomposition are calculated, including: Based on the direct numerical simulation data of the particle two-phase wall turbulence, the normal enthalpy flow and normal pressure energy flow effects, the friction effect, the turbulence pulsation effect, and the particle heat transfer effect are calculated.

4. The method for calculating the heat flux of the turbulent boundary layer of the ablated particle wall according to claim 3, characterized in that: Based on the direct numerical simulation data of the particle two-phase wall turbulence, the particle two-phase viscosity effect, compressibility effect, turbulence pulse effect, wake effect, pressure gradient effect, particle average size effect and particle pulse size effect are calculated, including: Calculate the viscous force effect based on the boundary layer thickness and fluid dynamic viscosity; Calculating a compressibility effect based on the boundary layer thickness, the fluid dynamic viscosity, the average velocity of the fluid flow direction, and the fluid velocity divergence; Calculating turbulent pulse effects based on the boundary layer thickness, the fluid dynamic viscosity, the fluid density, and the Reynolds stress; Calculating a wake effect based on the boundary layer thickness, the fluid density, and the fluid kinetic energy; Calculate the pressure gradient effect based on fluid flow and fluid pressure; Calculating the average particle size effect based on the boundary layer thickness, the average flow velocity of the fluid, the apparent density of the particle group, the particle inertia relaxation coefficient, and the average flow velocity of the particles; The particle pulse scale effect is calculated according to the boundary layer thickness, the particle flow pulsating velocity, the fluid flow pulsating velocity, the particle group apparent density, and the particle inertial relaxation coefficient.

5. The method for calculating the heat flux of the turbulent boundary layer of the ablated particle wall according to claim 4, characterized in that: Based on the direct numerical simulation data of the particle two-phase wall turbulence, the normal enthalpy flow and normal pressure energy flow effects, the friction effect, the turbulent pulsation effect, and the particle heat transfer effect are calculated, including: Calculating normal enthalpy flow and normal pressure energy flow effects based on the fluid enthalpy, the fluid normal average velocity, and the fluid pressure; Calculating the effect of friction based on the boundary layer thickness, the average velocity of the fluid flow, and the fluid shear stress; Calculating turbulence pulsation effects based on the boundary layer thickness, the fluid specific heat capacity, the fluid pressure, and normal turbulence temperature pulsation; The particle heat transfer effect is calculated based on the average temperature of the particles, the average temperature of the fluid, the apparent density of the particle group, the specific heat capacity of the particles, and the thermal response coefficient of the particles.

6. An electronic device, characterized in that: The method comprises a processor and a memory, wherein when the processor executes the computer program stored in the memory, the method for calculating the friction heat flow of the turbulent boundary layer of the ablated particle wall according to any one of claims 1 to 5 is implemented.

7. A computer-readable storage medium, characterized in that Used to store a computer program, wherein when the computer program is executed by a processor, the method for calculating the frictional heat flow of the turbulent boundary layer of the ablated particle wall according to any one of claims 1 to 5 is implemented.

Images