A method, device, equipment and storage medium for numerical simulation of internal and external flow coupling

Abstract

The present invention discloses a method, device, equipment and storage medium for numerical simulation of internal and external flow coupling, which relates to the field of numerical calculation technology, determines the flow state at the initial moment to obtain the external field wall temperature, solves the external field multi-component chemical non-equilibrium N-S equation according to the external field fluid boundary condition, obtains the external fluid boundary heat flow at the next moment and the external field pressure at the next moment, which are used as the coupling interface boundary condition for solving the internal field, and solves the internal flow of the porous sweating structure through the enthalpy method equation group of TPMM and LTNE to obtain the parameters at the coupling boundary at the next moment, and finally obtains the high enthalpy non-equilibrium external field parameters at the termination moment. By using the N-S equation under the assumption of multi-component chemical non-equilibrium to solve the supersonic outflow, the chemical reaction between the active cooling medium and the supersonic incoming flow under supersonic conditions is simulated, and the wall heat flow and structural thermal response of the porous sweating cooling system are obtained.

Description

Technical Field

[0001] The present invention relates to the field of numerical computing technology, and in particular to a method, device, equipment and storage medium for numerical simulation of internal and external flow coupling. Background Art

[0002] When using numerical simulation methods to study the mechanism and design the scheme of the sweating cooling system, it is necessary to couple the heat and mass transport process of the high Mach number supersonic (Ma>10) external field and the porous structure internal field. At present, the control equations of different regions are usually used to solve the fluid and solid regions respectively, and the internal and external coupling is realized through the interface parameter transfer.

[0003] In the existing technology, the calculation method does not consider the chemical reaction between the active cooling medium and the supersonic flow under supersonic conditions when calculating the external field, which leads to large prediction deviations in the calculation results of the heat flux and thermal response of the solid porous structure wall, which seriously affects the fine design of high Mach number aircraft. Therefore, it is necessary to develop a numerical simulation of the internal and external flow coupling that considers the chemical non-equilibrium effect for the heat and mass transport process of the active sweating cooling system under supersonic conditions. Summary of the invention

[0004] The purpose of the present invention is to provide a method, device, equipment and storage medium for numerical simulation of internal and external flow coupling in order to overcome the existing technical defects. By adopting the NS equation under the assumption of multi-component chemical non-equilibrium to solve the supersonic outflow, the chemical reaction between the active cooling medium and the supersonic incoming flow under supersonic conditions is simulated to obtain accurate wall heat flux and structural thermal response of the porous sweating cooling system.

[0005] The object of the present invention is achieved through the following technical solutions:

[0006] In a first aspect, the present invention provides a method for numerical simulation of internal and external flow coupling, the method comprising:

[0007] Step S1, determining the flow state at the initial moment, and obtaining the corresponding physical parameters at the gas-solid interface, wherein the physical parameters include the internal field fluid temperature, the internal field solid temperature, the coupling interface pressure, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction, and the cooling medium gas phase mass fraction;

[0008] Step S2, using the internal field fluid temperature and the internal field solid temperature to obtain the external field wall temperature by weighted average of the porosity of the porous structure, and taking the external field wall temperature, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction and the cooling medium gas phase mass fraction at the initial moment as the external field fluid boundary conditions;

[0009] Step S3, solving the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary condition, obtaining the external fluid boundary heat flow at the next moment and the external field pressure at the next moment, and using them as the coupling interface boundary condition for solving the internal field;

[0010] Step S4, based on the coupling interface boundary conditions solved by the internal field, the internal flow of the porous sweating structure is solved by the enthalpy method equation group of TPMM and LTNE to obtain the internal field fluid temperature, internal field solid temperature, cooling medium mass flow rate, cooling medium liquid phase mass fraction and cooling medium gas phase mass fraction at the coupling boundary at the next moment;

[0011] Step S5, iterating steps S2 to S4 until a predetermined termination time is reached or a convergence condition is satisfied, and obtaining high enthalpy non-equilibrium external field parameters at the termination time.

[0012] In a possible implementation, the external wall temperature for ,in is the porosity of the porous structure, is the internal solid temperature, is the internal fluid temperature.

[0013] In one possible implementation, the external field multi-component chemical non-equilibrium NS equation includes: , , , , ,in, For components The density of For components The diffusion coefficient, For components The quality score, For components The mass generation rate source term, For components The molar mass of is the density of the mixture, , are the velocity components in different directions, For pressure, is the temperature, is the universal gas constant, is the total energy of the mixture, is the viscous stress tensor, is the heat flux, is the Kronecker symbol, represents the total number of body components, For time, and They represent different direction coordinates respectively.

[0014] In the case of chemical reactions, components The mass generation rate source term for: ,in and Represent the components in the i-th reaction The stoichiometric number in the reaction and the product, and are the forward and reverse reaction rate coefficients of the elementary reaction, is the total number of chemical reactions;

[0015] In the case of gas-liquid phase transition, the components The mass generation rate source term for: ,in is the model coefficient, is the liquid density, is the saturation temperature.

[0016] In a possible implementation, the steps of solving the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary condition, obtaining the external fluid boundary heat flux at the next moment and the external field pressure at the next moment, and using them as the coupling interface boundary condition for the internal field solution include:

[0017] Solve the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary conditions to obtain the external fluid boundary heat flow and the external field pressure at the next moment;

[0018] Decompose the external fluid boundary heat flux into convection heat flux at the next moment , diffusion heat flow and radiant heat flux ;

[0019] Based on convection heat flow , diffusion heat flow and radiant heat flux , using the porosity of the porous structure and the weighted gas-solid thermal conductivity to calculate the fluid boundary heat flow in the internal field and solid boundary heat flux , and denote the fluid and solid temperatures at the boundary, respectively. represents the porosity;

[0020] The fluid boundary heat flux and solid boundary heat flux of the internal field are used as the coupled interface boundary conditions for solving the internal field.

[0021] In one possible implementation, the enthalpy method equations of TPMM and LTNE are: , , and ,in represents enthalpy, represents the permeability of porous materials, is the porosity of the porous structure, is the convection correction coefficient of the enthalpy equation, is the effective thermal diffusion coefficient of the enthalpy equation, represents the heat transfer source term between the internal solid and fluid, is the dynamic viscosity of the mixture, is the specific heat capacity of porous materials, is the porous material temperature, is the equivalent thermal conductivity of the porous material.

[0022] In a possible implementation, the high enthalpy non-equilibrium external field parameters at the termination time include external field velocity, density, temperature, pressure, component distribution, and fluid temperature distribution and solid temperature distribution in the porous medium.

[0023] In a second aspect, the present invention provides a device for numerical simulation of internal and external flow coupling, the device comprising:

[0024] A determination module, used to determine the flow state at the initial moment, and obtain the corresponding physical parameters at the gas-solid interface, wherein the physical parameters include the internal field fluid temperature, the internal field solid temperature, the coupling interface pressure, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction, and the cooling medium gas phase mass fraction;

[0025] The weighted average module is used to obtain the external field wall temperature by weighted average of the internal field fluid temperature and the internal field solid temperature through the porosity of the porous structure, and the external field wall temperature, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction and the cooling medium gas phase mass fraction at the initial moment are used as the external field fluid boundary conditions;

[0026] The first solution module is used to solve the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary conditions, obtain the external fluid boundary heat flow at the next moment and the external field pressure at the next moment, and use them as the coupling interface boundary conditions for the internal field solution;

[0027] The second solution module is used to solve the coupling interface boundary conditions based on the internal field solution, and to obtain the internal field fluid temperature, internal field solid temperature, cooling medium mass flow rate, cooling medium liquid phase mass fraction and cooling medium gas phase mass fraction at the coupling boundary at the next moment through the enthalpy method equation group of TPMM and LTNE;

[0028] The iteration module is used to iterate until a predetermined termination time is reached or a convergence condition is met, and obtain the high enthalpy non-equilibrium external field parameters at the termination time.

[0029] In a third aspect, the present invention further proposes a computer device, comprising a processor and a memory, wherein the memory stores a computer program, and the computer program is loaded and executed by the processor to implement the internal and external flow coupling numerical simulation method as described in any one of the first aspects.

[0030] In a fourth aspect, the present invention further proposes a computer-readable storage medium, in which a computer program is stored. The computer program is loaded and executed by a processor to implement the internal and external flow coupling numerical simulation method as described in any one of the first aspects.

[0031] The above main scheme of the present invention and its further options can be freely combined to form multiple schemes, all of which are schemes that can be adopted and claimed for protection by the present invention; and in the present invention, (non-conflicting options) options and other options can also be freely combined. After understanding the scheme of the present invention, those skilled in the art can understand that there are many combinations based on the existing technology and common knowledge, all of which are technical schemes to be protected by the present invention, and they are not exhaustively listed here.

[0032] The present invention discloses a method, device, equipment and storage medium for numerical simulation of internal and external flow coupling, firstly determining the flow state at the initial moment, obtaining the external field wall temperature, using the external field wall temperature, cooling medium mass flow rate, cooling medium liquid phase mass fraction and cooling medium gas phase mass fraction at the initial moment as the external field fluid boundary conditions to solve the external field multi-component chemical non-equilibrium NS equation, obtaining the external fluid boundary heat flow at the next moment and the external field pressure at the next moment, using them as the coupling interface boundary conditions for solving the internal field, solving the internal flow of the porous sweating structure through the enthalpy method equation group of TPMM and LTNE to obtain the parameters at the coupling boundary at the next moment, and finally obtaining the high enthalpy non-equilibrium external field parameters at the termination moment. By using the NS equation under the assumption of multi-component chemical non-equilibrium to solve the supersonic outflow, simulating the chemical reaction between the active cooling medium and the supersonic incoming flow under supersonic conditions, accurate wall heat flow and structural thermal response of the porous sweating cooling system are obtained. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments are briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without creative work.

[0034] Figure 1 A schematic flow chart of a method for numerical simulation of internal and external flow coupling proposed in an embodiment of the present invention is shown.

[0035] Figure 2A schematic diagram of the interface coupling strategy proposed in an embodiment of the present invention is shown. DETAILED DESCRIPTION

[0036] The following describes the embodiments of the present invention by specific examples, and those skilled in the art can easily understand other advantages and effects of the present invention from the contents disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and the details in this specification can also be modified or changed in various ways based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that the following embodiments and features in the embodiments can be combined with each other without conflict.

[0037] Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in the field without making any creative work shall fall within the scope of protection of the present invention.

[0038] In the prior art, when using numerical simulation methods to study the sweating cooling system, it is necessary to couple the heat and mass transport processes of the supersonic external field and the porous structure internal field. At present, different control equations are usually used to solve the fluid and solid regions respectively, and the internal and external coupling is achieved through interface parameter transfer. The existing calculation method does not consider the chemical reaction between the active cooling fluid and the supersonic incoming flow under supersonic conditions when calculating the external field, which leads to deviations in the calculation results of the heat flux and thermal response of the solid porous structure wall.

[0039] Therefore, in order to solve the above-mentioned technical problems, an embodiment of the present invention proposes a method, device, equipment and storage medium for numerical simulation of internal and external flow coupling, and adopts different control equations to solve the supersonic outflow and the internal flow field and temperature field of the porous structure respectively, and realizes internal and external coupling through parameter transfer at the interface, wherein the supersonic outflow is solved by the Navier-Stokes (NS) equation under the assumption of multi-component chemical non-equilibrium, and the internal flow field and temperature field of the porous structure are solved by the enthalpy method based on the two-phase mixture model (Two-Phase Mixture Model, TPMM) and the porous medium thermal non-equilibrium model (Local Thermal Non-Equilibrium, LTNE), simulating the chemical reaction between the active cooling working fluid and the supersonic incoming flow under supersonic conditions, and obtaining accurate wall heat flux and structural thermal response of the porous sweating cooling system, which will be described in detail below.

[0040] The main application scenarios of the internal and external flow coupling numerical simulation method proposed in the embodiment of the present invention may include the following aspects:

[0041] First, during high-speed flight, the surface of the aircraft will be subjected to extremely high heat flux density due to aerodynamic heating. The active sweating cooling system sprays cooling medium to reduce the wall temperature to protect the aircraft structure. This method can simulate the working state of the thermal protection system in an extreme thermal environment, evaluate the cooling efficiency and optimize the spraying method and parameters of the cooling medium.

[0042] Second, this method can accurately capture the complex coupling effects between internal and external fields (i.e., the interaction between the cooling medium and the high-temperature airflow), and predict the surface heat flux distribution and temperature field of the aircraft under different flight conditions, thus providing a basis for thermal load assessment and material selection.

[0043] Third, in addition to the thermal load of the external flow field, this method can also simulate the impact of the cooling system on the internal structure of the aircraft, help analyze the temperature distribution and thermal stress inside the structure, and provide support for structural strength and life assessment.

[0044] Please refer to Figure 1 and Figure 2 , Figure 1 A schematic flow chart of a method for numerical simulation of internal and external flow coupling proposed in an embodiment of the present invention is shown, and the method includes the following steps:

[0045] Step S1, determine the flow state at the initial moment and obtain the corresponding physical parameters at the gas-solid interface.

[0046] The physical parameters include internal field fluid temperature, internal field solid temperature, coupled interface pressure, cooling medium mass flow rate, cooling medium liquid phase mass fraction and cooling medium gas phase mass fraction.

[0047] Figure 2 The schematic diagram of the interface coupling strategy proposed in the embodiment of the present invention is shown. Figure 2 As a basis for further explanation of the internal and external flow coupling numerical simulation method, the internal field fluid temperature is the temperature of the fluid inside the porous structure, and the solid temperature is the temperature of the solid part inside the porous structure, and the external field pressure is the pressure in the supersonic outflow region, Represents the pressure of the internal fluid boundary and the mass flow rate of the cooling medium Refers to the mass flow rate of the cooling medium passing through the porous structure, the mass fraction of the cooling medium liquid phase Refers to the mass fraction of the liquid phase in the cooling medium and the mass fraction of the gas phase in the cooling medium It is the mass fraction of the gas phase in the cooling fluid.

[0048] First, it is necessary to set the flow state at the initial moment, including parameters such as the initial temperature, pressure and flow velocity inside and outside the porous structure. At the interface between the porous structure and the external flow field, the above six physical parameters are measured or calculated. They are used in subsequent numerical simulations to ensure the coupled solution of the internal and external flow fields.

[0049] Step S2, using the internal field fluid temperature and the internal field solid temperature to obtain the external field wall temperature by weighted average of the porosity of the porous structure, and taking the external field wall temperature, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction and the cooling medium gas phase mass fraction at the initial moment as the external field fluid boundary conditions.

[0050] External wall temperature for ,in is the porosity of the porous structure, is the internal solid temperature, is the internal fluid temperature.

[0051] Then, some key physical parameters at the initial moment are set as boundary conditions of the external field fluid, including the external field wall temperature, the mass flow rate of the cooling medium, the mass fraction of the cooling medium liquid phase, and the mass fraction of the cooling medium gas phase. The external field wall temperature is obtained by weighted average of the internal field fluid temperature and solid temperature and the porosity of the porous structure. By extrapolation, the normal velocity of the external wall is obtained. .

[0052] Step S3, solving the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary conditions, obtaining the external fluid boundary heat flux at the next moment and the external field pressure at the next moment, and using them as the coupling interface boundary conditions for solving the internal field.

[0053] Using the external fluid boundary conditions (external wall temperature, cooling fluid mass flow rate, cooling fluid liquid phase mass fraction, and cooling fluid gas phase mass fraction) determined above as initial conditions, the NS equations are discretized in time and space, and numerical methods (Runge-Kutta method, implicit method, etc.) are used to solve the discretized NS equations to obtain the flow field variables (such as velocity, pressure, temperature, component concentration, etc.) at the next moment. Finally, the boundary heat flux of the external fluid at the next moment is extracted from the solution results. and pressure , which can accurately simulate multi-component chemical non-equilibrium flow and heat transfer processes.

[0054] The external field multi-component chemical non-equilibrium NS equations include: , , , , ,in, For components The density of For components The diffusion coefficient, For components The quality score, For components The mass generation rate source term, For components The molar mass of is the density of the mixture, , are the velocity components in different directions, For pressure, is the temperature, is the universal gas constant, is the total energy of the mixture, is the viscous stress tensor, is the heat flux, is the Kronecker symbol, represents the total number of body components, For time, and They represent different direction coordinates respectively.

[0055] In the case of chemical reactions, components The mass generation rate source term for: ,in and Represent the components in the i-th reaction The stoichiometric number in the reaction and the product, and are the forward and reverse reaction rate coefficients of the elementary reaction, is the total number of chemical reactions;

[0056] In the case of gas-liquid phase transition, the components The mass generation rate source term for: ,in is the model coefficient, is the liquid density, is the saturation temperature.

[0057] The liquid and gas phases of the cooling medium are considered as two different components for calculation. Specifically, the gas phase medium participates in the chemical non-equilibrium reaction system in the ultrasonic flow field, while the liquid phase medium only participates in the phase change reaction. For the gas-liquid phase change, the mass source term is obtained through the gas-liquid phase change model Lee model, and the saturation temperature It can be obtained by fitting and interpolating the REFPROP database according to the pressure p. , which can be approximately considered .

[0058] The steps of solving the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary condition, obtaining the external fluid boundary heat flux at the next moment and the external field pressure at the next moment, and using them as the coupling interface boundary condition for the internal field solution include:

[0059] Solve the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary conditions to obtain the external fluid boundary heat flow and the external field pressure at the next moment;

[0060] Decompose the external fluid boundary heat flux into convection heat flux at the next moment , diffusion heat flow and radiant heat flux ;

[0061] Based on convection heat flow , diffusion heat flow and radiant heat flux , using the porosity of the porous structure and the weighted gas-solid thermal conductivity to calculate the fluid boundary heat flow in the internal field and solid boundary heat flux , and denote the fluid and solid temperatures at the boundary, respectively. represents the porosity;

[0062] The fluid boundary heat flux and solid boundary heat flux of the internal field are used as the coupled interface boundary conditions for solving the internal field.

[0063] The external field heat flux and wall pressure are obtained, and the wall heat flux is decomposed into convection heat flux, diffusion heat flux and radiation heat flux. The fluid boundary heat flux and solid boundary heat flux of the internal field are calculated by weighting the porosity of the porous structure and the gas-solid thermal conductivity. The calculated fluid boundary heat flux and solid boundary heat flux of the internal field are used as the coupled interface boundary conditions for solving the internal field.

[0064] Step S4, based on the coupling interface boundary conditions solved by the internal field, the internal flow of the porous sweating structure is solved by the enthalpy method equation group of TPMM and LTNE to obtain the internal field fluid temperature, internal field solid temperature, cooling fluid mass flow rate, cooling fluid liquid phase mass fraction and cooling fluid gas phase mass fraction at the coupling boundary at the next moment.

[0065] Using the internal fluid and solid boundary conditions determined above, the internal flow of the porous sweating structure is solved by the enthalpy method equations of TPMM and LTNE to obtain the internal fluid temperature at the coupled boundary at the next moment. , solid temperature , cooling medium mass flow , Cooling fluid liquid mass fraction and cooling medium gas phase mass fraction .

[0066] The enthalpy method equations for TPMM and LTNE are: , , and ,in represents enthalpy, represents the permeability of porous materials, is the porosity of the porous structure, is the convection correction coefficient of the enthalpy equation, is the effective thermal diffusion coefficient of the enthalpy equation, represents the heat transfer source term between the internal solid and fluid, is the dynamic viscosity of the mixture, is the specific heat capacity of porous materials, is the porous material temperature, is the equivalent thermal conductivity of the porous material.

[0067] Step S5, iterating steps S2 to S4 until a predetermined termination time is reached or a convergence condition is satisfied, and obtaining high enthalpy non-equilibrium external field parameters at the termination time.

[0068] The high enthalpy non-equilibrium external field parameters at the termination time include external field velocity, density, temperature, pressure, component distribution, and the fluid temperature distribution and solid temperature distribution in the porous medium.

[0069] The coupled iterative method is used, that is, the enthalpy method equations of TPMM and LTNE are combined with the coupled boundary conditions of the internal and external fields to perform iterative solutions. The iterative calculation will continue. In each iteration, the boundary conditions are updated according to the current solution, and the equations are re-solved until the predetermined termination time is reached or the convergence conditions are met. Until the predetermined termination time is reached, the velocity, density, temperature, pressure and component distribution of the high enthalpy non-equilibrium external field are obtained at the termination time. These parameters describe the state of the external field at the termination time. At the same time, the fluid temperature distribution and solid temperature distribution of the internal field of the porous medium can also be obtained, which describes the temperature state inside the porous medium at the termination time.

[0070] Compared with the prior art, the embodiments of the present invention have the following beneficial effects:

[0071] First, by solving the supersonic outflow using the Navier-Stokes (NS) equations under the assumption of multi-component chemical non-equilibrium, the chemical reaction between the cooling medium and the incoming flow under supersonic conditions can be simulated more accurately. This assumption takes into account the dynamic process of the chemical reaction rather than a simple equilibrium state, thereby improving the accuracy of the simulation.

[0072] Second, inside the porous structure, the enthalpy method based on TPMM and LTNE is used to solve the flow field and temperature field. These models can more realistically reflect the flow and heat transfer process inside the porous medium, especially considering the interaction between the gas and liquid phases and the thermal non-equilibrium effect.

[0073] Third, by transferring parameters at the interface to achieve the coupled solution of internal and external fields, the simulation of the entire system is more complete and consistent. It is not only applicable to single-phase fluids, but can also handle complex multiphase fluid systems, enhancing the applicability of the model.

[0074] Fourth, the liquid and gas phases of the coolant are considered as two different components for calculation. The gas phase works in the chemical non-equilibrium reaction system in the ultrasonic flow field, and the liquid phase works only in the phase change reaction. This treatment method can more accurately simulate the behavior of the coolant in the porous sweating cooling system, especially under high temperature and high speed conditions.

[0075] Fifth, by iterating the calculation until the termination time, the convergence and stability of the simulation results are ensured. In each iteration, the boundary conditions of the internal and external fields are updated to gradually approach the final solution, which improves the efficiency of the calculation.

[0076] Sixth, this method can provide detailed internal and external flow field and temperature field distribution, which helps to optimize the structural design of the porous sweating cooling system and improve the cooling efficiency and thermal protection performance.

[0077] A possible implementation of an internal and external flow coupling numerical simulation device is provided below, which is used to execute each execution step and corresponding technical effect of the internal and external flow coupling numerical simulation method shown in the above embodiment and possible implementation, and the device includes:

[0078] A determination module, used to determine the flow state at the initial moment, and obtain the corresponding physical parameters at the gas-solid interface, wherein the physical parameters include the internal field fluid temperature, the internal field solid temperature, the coupling interface pressure, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction, and the cooling medium gas phase mass fraction;

[0079] The weighted average module is used to obtain the external field wall temperature by weighted average of the internal field fluid temperature and the internal field solid temperature through the porosity of the porous structure, and the external field wall temperature, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction and the cooling medium gas phase mass fraction at the initial moment are used as the external field fluid boundary conditions;

[0080] The first solution module is used to solve the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary conditions, obtain the external fluid boundary heat flow at the next moment and the external field pressure at the next moment, and use them as the coupling interface boundary conditions for the internal field solution;

[0081] The second solution module is used to solve the coupling interface boundary conditions based on the internal field solution, and to obtain the internal field fluid temperature, internal field solid temperature, cooling medium mass flow rate, cooling medium liquid phase mass fraction and cooling medium gas phase mass fraction at the coupling boundary at the next moment through the enthalpy method equation group of TPMM and LTNE;

[0082] The iteration module is used to iterate until a predetermined termination time is reached or a convergence condition is met, and obtain the high enthalpy non-equilibrium external field parameters at the termination time.

[0083] This preferred embodiment provides a computer device, which can implement the steps in any embodiment of the internal and external flow coupling numerical simulation method provided in the embodiment of the present invention. Therefore, the beneficial effects of the internal and external flow coupling numerical simulation method provided in the embodiment of the present invention can be achieved. Please refer to the previous embodiments for details and will not be repeated here.

[0084] Those skilled in the art can understand that all or part of the steps in the various methods of the above embodiments can be completed by instructions, or by controlling related hardware through instructions, and the instructions can be stored in a computer-readable storage medium and loaded and executed by a processor. To this end, an embodiment of the present invention provides a storage medium, which stores multiple instructions, and the instructions can be loaded by a processor to execute the steps of any embodiment of the internal and external flow coupling numerical simulation method provided by the embodiment of the present invention.

[0085] The storage medium may include: a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, etc.

[0086] Since the instructions stored in the storage medium can execute the steps in any embodiment of the internal and external flow coupling numerical simulation method provided in the embodiments of the present invention, the beneficial effects that can be achieved by any internal and external flow coupling numerical simulation method provided in the embodiments of the present invention can be achieved. Please refer to the previous embodiments for details and will not be repeated here.

[0087] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for numerical simulation of internal and external flow coupling, characterized in that: The method comprises: Step S1, determining the flow state at the initial moment, and obtaining the corresponding physical parameters at the gas-solid interface, wherein the physical parameters include the internal field fluid temperature, the internal field solid temperature, the coupling interface pressure, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction, and the cooling medium gas phase mass fraction; Step S2, using the internal field fluid temperature and the internal field solid temperature to obtain the external field wall temperature by weighted average of the porosity of the porous structure, and taking the external field wall temperature, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction and the cooling medium gas phase mass fraction at the initial moment as the external field fluid boundary conditions; Step S3, solving the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary condition to obtain the external fluid boundary heat flow at the next moment and the external field pressure at the next moment; Decompose the external fluid boundary heat flux into convection heat flux at the next moment , diffusion heat flow and radiant heat flux ; Based on convection heat flow , diffusion heat flow and radiant heat flux , using the porosity of the porous structure and the weighted gas-solid thermal conductivity to calculate the fluid boundary heat flow in the internal field and solid boundary heat flux , and denote the fluid and solid temperatures at the boundary, respectively. is the porosity of the porous structure; The fluid boundary heat flux and solid boundary heat flux of the internal field are used as the coupled interface boundary conditions for solving the internal field. The external field multi-component chemical non-equilibrium NS equations include: , , , , ,in, For components The density of For components The diffusion coefficient, For components The quality score, For components The mass generation rate source term, For components The molar mass of is the density of the mixture, , are the velocity components in different directions, For pressure, is the temperature, is the universal gas constant, is the total energy of the mixture, is the viscous stress tensor, is the heat flux, is the Kronecker symbol, represents the total number of body components, For time, and Represent different direction coordinates respectively; In the case of chemical reactions, components The mass generation rate source term for: ,in and Represent the components in the i-th reaction The stoichiometric number in the reaction and the product, and are the forward and reverse reaction rate coefficients of the elementary reaction, is the total number of chemical reactions; In the case of gas-liquid phase transition, the components The mass generation rate source term for: ,in is the model coefficient, is the liquid density, is the saturation temperature; Step S4, based on the coupling interface boundary conditions solved by the internal field, the internal flow of the porous sweating structure is solved by the enthalpy method equation group of TPMM and LTNE to obtain the internal field fluid temperature, internal field solid temperature, cooling medium mass flow rate, cooling medium liquid phase mass fraction and cooling medium gas phase mass fraction at the coupling boundary at the next moment; Step S5, iterating steps S2 to S4 until a predetermined termination time is reached or a convergence condition is satisfied, and obtaining high enthalpy non-equilibrium external field parameters at the termination time.

2. The internal and external flow coupling numerical simulation method according to claim 1, characterized in that: External wall temperature for ,in is the porosity of the porous structure, is the internal solid temperature, is the internal fluid temperature.

3. The internal and external flow coupling numerical simulation method according to claim 1, characterized in that: The enthalpy method equations for TPMM and LTNE are: , , and ,in represents enthalpy, represents the permeability of porous materials, is the porosity of the porous structure, is the convection correction coefficient of the enthalpy equation, is the effective thermal diffusion coefficient of the enthalpy equation, represents the heat transfer source term between the internal solid and fluid, is the dynamic viscosity of the mixture, is the specific heat capacity of porous materials, is the porous material temperature, is the equivalent thermal conductivity of the porous material.

4. The internal and external flow coupling numerical simulation method according to claim 1, characterized in that: The high enthalpy non-equilibrium external field parameters at the termination time include external field velocity, density, temperature, pressure, component distribution, and the fluid temperature distribution and solid temperature distribution in the porous medium.

5. A numerical simulation device for coupling internal and external flow, characterized in that: The device comprises: A determination module, used to determine the flow state at the initial moment, and obtain the corresponding physical parameters at the gas-solid interface, wherein the physical parameters include the internal field fluid temperature, the internal field solid temperature, the coupling interface pressure, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction, and the cooling medium gas phase mass fraction; The weighted average module is used to obtain the external field wall temperature by weighted average of the internal field fluid temperature and the internal field solid temperature through the porosity of the porous structure, and the external field wall temperature, the cooling medium mass flow rate, the cooling medium liquid phase mass fraction and the cooling medium gas phase mass fraction at the initial moment are used as the external field fluid boundary conditions; The first solving module is used to solve the external field multi-component chemical non-equilibrium NS equation based on the external field fluid boundary conditions to obtain the external fluid boundary heat flow at the next moment and the external field pressure at the next moment; Decompose the external fluid boundary heat flux into convection heat flux at the next moment , diffusion heat flow and radiant heat flux ; Based on convection heat flow , diffusion heat flow and radiant heat flux , using the porosity of the porous structure and the weighted gas-solid thermal conductivity to calculate the fluid boundary heat flow in the internal field and solid boundary heat flux , and denote the fluid and solid temperatures at the boundary, respectively. is the porosity of the porous structure; The fluid boundary heat flux and solid boundary heat flux of the internal field are used as the coupled interface boundary conditions for solving the internal field. The external field multi-component chemical non-equilibrium NS equations include: , , , , ,in, For components The density of For components The diffusion coefficient, For components The quality score, For components The mass generation rate source term, For components The molar mass of is the density of the mixture, , are the velocity components in different directions, For pressure, is the temperature, is the universal gas constant, is the total energy of the mixture, is the viscous stress tensor, is the heat flux, is the Kronecker symbol, represents the total number of body components, For time, and Represent different direction coordinates respectively; In the case of chemical reactions, components The mass generation rate source term for: ,in and Represent the components in the i-th reaction The stoichiometric number in the reaction and the product, and are the forward and reverse reaction rate coefficients of the elementary reaction, is the total number of chemical reactions; In the case of gas-liquid phase transition, the components The mass generation rate source term for: ,in is the model coefficient, is the liquid density, is the saturation temperature; The second solution module is used to solve the coupling interface boundary conditions based on the internal field solution, and to obtain the internal field fluid temperature, internal field solid temperature, cooling medium mass flow rate, cooling medium liquid phase mass fraction and cooling medium gas phase mass fraction at the coupling boundary at the next moment through the enthalpy method equation group of TPMM and LTNE; The iteration module is used to iterate until a predetermined termination time is reached or a convergence condition is met, and obtain the high enthalpy non-equilibrium external field parameters at the termination time.

6. A computer device, characterized in that: The computer device includes a processor and a memory, wherein a computer program is stored in the memory, and the computer program is loaded and executed by the processor to implement the internal and external flow coupling numerical simulation method as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: The storage medium stores a computer program, which is loaded and executed by a processor to implement the internal and external flow coupling numerical simulation method as described in any one of claims 1 to 4.

Images