A gas-liquid-ice three-phase flow analysis method, device, equipment and medium

Abstract

The present application discloses a gas-liquid-ice three-phase flow analysis method, device, equipment and medium, which relates to the field of phase state analysis, including: in a calculation domain determined based on industrial equipment, determining the velocity vector of the fluid at the boundary in the calculation domain based on the current liquid-ice phase function to update the current velocity vector of the fluid; inputting the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into a coupling relationship to update the current vortex vector and the current velocity vector of the fluid and predict the new current gas-liquid phase function, and then updating the current temperature of the fluid in combination with the current liquid ice phase function and the current mass flux to update the current mass flux, the current liquid ice phase function and the current gas-liquid phase function based on the updated current temperature, thereby determining the phase state at any position in the calculation domain in real time, realizing autonomous analysis of the dynamic freezing process, filling the gap in the field of gas-liquid-ice three-phase flow analysis, and providing support for analyzing freezing problems.

Description

Technical Field

[0001] The present invention relates to the field of phase analysis, and in particular to a gas-liquid-ice three-phase flow analysis method, device, equipment and medium. Background Art

[0002] The three-phase flow problems of gas, liquid and ice include single-phase flow, gas-liquid two-phase flow, and liquid-ice two-phase flow. The gas phase includes air and water vapor. Water vapor converts between gas and liquid through evaporation and condensation, and the liquid phase converts between liquid and ice through freezing and melting. Since most industrial equipment has a complex shape, the complex vortex system it brings affects the conversion process of gas-liquid and liquid-ice. When the industrial equipment is an aircraft, an important problem faced by the aircraft is ice-tolerant flight. The research and development of ice-tolerant flight requires independent and controllable icing simulation software. The icing problem is actually a three-phase flow problem of gas, liquid and ice. Therefore, how to realize the three-phase flow analysis of gas, liquid and ice is a problem that needs to be solved. Summary of the invention

[0003] In view of this, the purpose of the present invention is to provide a gas-liquid-ice three-phase flow analysis method, device, equipment and medium, which can determine the phase state at any position in the calculation domain in real time, realize the autonomous analysis of the dynamic freezing process, fill the gap in the field of gas-liquid-ice three-phase flow analysis, and provide strong support for the analysis of freezing problems. The specific scheme is as follows:

[0004] In a first aspect, the present application provides a gas-liquid-ice three-phase flow analysis method, comprising:

[0005] In a calculation domain predetermined based on industrial equipment, a velocity vector of a fluid at a boundary in the calculation domain is determined based on a current liquid-ice phase function to update a current velocity vector of the fluid;

[0006] Inputting the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to update the current vortex vector and the current velocity vector of the fluid and predict a new current gas-liquid phase function; wherein the current mass flux includes the current gas-liquid evaporation and condensation mass flux and the current liquid ice freezing and melting mass flux;

[0007] The current temperature of the fluid is updated using the current liquid ice phase function, the current mass flux, the updated current velocity vector and the predicted current gas-liquid phase function, and the current mass flux, the current liquid ice phase function and the current gas-liquid phase function are updated based on the updated current temperature to determine the phase state at any position in the calculation domain and analyze the dynamic freezing process.

[0008] Optionally, determining the velocity vector of the fluid at the boundary in the computational domain based on the current liquid-ice phase function includes:

[0009] When the current liquid-ice phase function characterizes that the current phase state is an ice state, determining the velocity vector of the fluid at the boundary in the calculation domain to be zero;

[0010] When the current liquid-ice phase function indicates that the current phase state is not an ice state, a velocity vector of the fluid at the boundary in the calculation domain is determined based on the current liquid-ice phase function and the current velocity vector of the fluid.

[0011] Optionally, the inputting the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to update the current vortex vector and the current velocity vector of the fluid and predict the new current gas-liquid phase function comprises:

[0012] Inputting the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function, so as to determine the current vortex vector based on the current displacement vector, the current velocity vector and the curl operator of the fluid, and updating the current vortex vector using the apparent kinematic viscosity and the surface tension vector of the fluid;

[0013] A flow function vector is determined based on a current displacement vector and an updated current vortex vector, and a velocity potential function is determined based on the current displacement vector, the current velocity vector and the current mass flux, so as to update the current velocity vector using the flow function vector and the velocity potential function, and predict a new current gas-liquid phase function using the current displacement vector, the current mass flux and the updated current velocity vector.

[0014] Optionally, before updating the current vortex vector using the apparent kinematic viscosity and surface tension vector of the fluid, the method further includes:

[0015] Determining the apparent dynamic viscosity of the fluid according to the dynamic viscosity of the liquid and the saturated water vapor and the current gas-liquid phase function, and determining the apparent density of the fluid according to the density of the liquid and the saturated water vapor and the current gas-liquid phase function, so as to determine the apparent kinematic viscosity of the fluid based on the ratio of the apparent dynamic viscosity to the apparent density;

[0016] The gas-liquid interface normal vector is determined according to the current displacement vector and the current gas-liquid phase function, and the curvature coefficient is determined according to the gas-liquid interface normal vector and the current displacement vector, so as to determine the surface tension vector based on the gas-liquid interface normal vector, the curvature coefficient, the current displacement vector and the current gas-liquid phase function.

[0017] Optionally, the coupling relationship includes:

[0018] ;

[0019] ;

[0020] ;

[0021] ;

[0022] ;

[0023] Wherein, the coupling relation adopts the Einstein summation convention, and Indicates three different directions. Indicates three different directions. Indicates three different directions. are vortex vectors, represents the curl operator, are velocity vectors, All represent displacement vectors, t represents time, represents the apparent kinematic viscosity, represents the surface tension vector, denotes the stream function vector, represents the velocity potential function, represents the density of saturated water vapor, represents the density of the liquid, represents the density of ice, represents the gas-liquid evaporation condensation mass flux, represents the mass flux of liquid ice freezing and melting, represents the gas-liquid phase function.

[0024] Optionally, updating the current temperature of the fluid using the current liquid-ice phase function, the current mass flux, the updated current velocity vector and the predicted current gas-liquid phase function comprises:

[0025] Determine the corresponding apparent heat capacity according to the heat capacity of the liquid and saturated water vapor and the predicted current gas-liquid phase function, and determine the corresponding apparent thermal conductivity according to the thermal conductivity of the liquid, saturated water vapor and ice and the current liquid-ice phase function and the predicted current gas-liquid phase function;

[0026] Determining a corresponding apparent thermal diffusion coefficient based on the apparent density, the apparent heat capacity and the apparent thermal conductivity; the apparent density is a density determined according to the density of the liquid and the saturated water vapor and the predicted current gas-liquid phase function;

[0027] updating a current temperature of the fluid using a current displacement vector of the fluid, the apparent thermal diffusivity, a current mass flux, and an updated current velocity vector;

[0028] Among them, the temperature update formula includes: ; T represents temperature, t represents time, represents the displacement vector, represents the velocity vector, is the apparent thermal diffusion coefficient, represents the gas-liquid evaporation condensation mass flux, represents the mass flux of liquid ice freezing and melting, represents the density of the liquid; and the updated formula for temperature uses the Einstein summation convention, Indicates three different directions.

[0029] Optionally, updating the current mass flux, the current liquid-ice phase function, and the current gas-liquid phase function based on the updated current temperature includes:

[0030] updating the current gas-liquid evaporation-condensation mass flux based on the gas-liquid phase change latent heat and the updated current temperature, and updating the current liquid ice freezing-melting mass flux based on the liquid ice phase change latent heat and the updated current temperature;

[0031] updating the current liquid-ice phase function based on the heat capacity and density of ice and the updated current temperature, and updating the current gas-liquid phase function based on the heat capacity and density of ice, the density of saturated water vapor, and the updated current temperature;

[0032] Among them, the update formulas for mass flux and phase function include:

[0033] ; ;

[0034] ; ;

[0035] represents the gas-liquid evaporation condensation mass flux, represents the latent heat of gas-liquid phase change, represents the apparent thermal conductivity, T represents the temperature, represents the displacement vector, represents the mass flux of liquid ice freezing and melting, is the latent heat of liquid-to-ice phase change, represents the density of ice, is the heat capacity of ice, represents the liquid-ice phase function, t represents time, represents the Laplace operator, represents the gas-liquid phase function, represents the density of saturated water vapor; and the updated formulas for mass flux and phase function use the Einstein summation convention, Indicates three different directions.

[0036] In a second aspect, the present application provides a gas-liquid-ice three-phase flow analysis device, comprising:

[0037] A first updating module is used to determine, in a calculation domain predetermined based on industrial equipment, a velocity vector of a fluid at a boundary in the calculation domain based on a current liquid-ice phase function, so as to update a current velocity vector of the fluid;

[0038] A second updating module is used to input the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to update the current vortex vector and the current velocity vector of the fluid and predict the new current gas-liquid phase function; wherein the current mass flux includes the current gas-liquid evaporation and condensation mass flux and the current liquid ice freezing and melting mass flux;

[0039] A phase state determination module is used to update the current temperature of the fluid using the current liquid ice phase function, the current mass flux, the updated current velocity vector and the predicted current gas-liquid phase function, and to update the current mass flux, the current liquid ice phase function and the current gas-liquid phase function based on the updated current temperature to determine the phase state at any position in the calculation domain and analyze the dynamic freezing process.

[0040] In a third aspect, the present application provides an electronic device, including:

[0041] Memory, used to store computer programs;

[0042] A processor is used to execute the computer program to implement the aforementioned gas-liquid-ice three-phase flow analysis method.

[0043] In a fourth aspect, the present application provides a computer-readable storage medium for storing a computer program, which, when executed by a processor, implements the aforementioned gas-liquid-ice three-phase flow analysis method.

[0044] In the present application, in a calculation domain predetermined based on industrial equipment, the velocity vector of the fluid at the boundary in the calculation domain is determined based on the current liquid ice phase function to update the current velocity vector of the fluid; the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid are input into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to update the current vortex vector and the current velocity vector of the fluid and predict the new current gas-liquid phase function; wherein the current mass flux includes the current gas-liquid evaporation and condensation mass flux and the current liquid ice freezing and melting mass flux; the current temperature of the fluid is updated using the current liquid ice phase function, the current mass flux, the updated current velocity vector and the predicted current gas-liquid phase function, and the current mass flux, the current liquid ice phase function and the current gas-liquid phase function are updated based on the updated current temperature to determine the phase state at any position in the calculation domain and analyze the dynamic freezing process. It can be seen that the present application gradually solves the vortex vector, velocity vector, gas-liquid phase function, temperature and mass flux at each moment to determine the gas-liquid phase function and liquid-ice phase function at each moment, thereby determining the phase state at any position in the calculation domain according to the gas-liquid phase function and liquid-ice phase function, and then realizing autonomous analysis of the dynamic freezing process. It not only fills the gap in the field of gas-liquid-ice three-phase flow analysis, realizes the autonomous control of autonomous source code software, but also helps to better study the mutual transformation between the three phases of gas, liquid and ice, and provides strong support for the analysis of freezing problems. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.

[0046] Figure 1 A flow chart of a gas-liquid-ice three-phase flow analysis method disclosed in this application;

[0047] Figure 2 This is a schematic structural diagram of a gas-liquid-ice three-phase flow analysis device disclosed in this application;

[0048] Figure 3 This is a structural diagram of an electronic device disclosed in this application. DETAILED DESCRIPTION

[0049] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. 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 creative work are within the scope of protection of the present invention.

[0050] Since most industrial equipment has a complex shape, the complex vortex system it brings affects the transformation process of gas-liquid two-phase and liquid-ice two-phase. When the industrial equipment is an aircraft, an important problem faced by the aircraft is ice-tolerant flight. The research and development design of ice-tolerant flight requires autonomous and controllable icing simulation software. The icing problem is actually a three-phase flow problem of gas, liquid and ice. To this end, this application provides a gas-liquid-ice three-phase flow analysis method, which realizes autonomous analysis of the dynamic icing process by real-time determination of the phase state at any position in the calculation domain, filling the gap in the field of gas-liquid-ice three-phase flow analysis and providing strong support for the analysis of icing problems.

[0051] See also Figure 1 As shown, the embodiment of the present invention discloses a gas-liquid-ice three-phase flow analysis method, comprising:

[0052] Step S11: In a calculation domain predetermined based on industrial equipment, a velocity vector of a fluid at a boundary in the calculation domain is determined based on a current liquid-ice phase function to update a current velocity vector of the fluid.

[0053] In this embodiment, in a calculation domain predetermined based on industrial equipment, the velocity vector of the fluid at the boundary in the calculation domain is first determined based on the current liquid-ice phase function, and the complete current velocity vector of the fluid is updated based on the velocity vector of the fluid at the boundary to obtain an updated current velocity vector.

[0054] Specifically, when the current liquid ice phase function represents that the current phase state is an ice state, the velocity vector of the fluid at the boundary in the calculation domain is determined to be zero; when the current liquid ice phase function represents that the current phase state is not an ice state, the velocity vector of the fluid at the boundary in the calculation domain is determined based on the current liquid ice phase function and the current velocity vector of the fluid.

[0055] Among them, the formula for determining the velocity vector of the fluid at the boundary includes:

[0056] ;

[0057] represents the velocity vector of the fluid at the boundary, represents the liquid-ice phase function, represents the complete velocity vector of the fluid, Indicates that the current phase is ice. indicates that the current phase is not ice; and the formula for determining the velocity vector of the fluid at the boundary adopts the Einstein summation convention, while Indicates three different directions corresponding to the computational domain.

[0058] Step S12, inputting the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to update the current vortex vector and the current velocity vector of the fluid and predict the new current gas-liquid phase function; wherein the current mass flux includes the current gas-liquid evaporation and condensation mass flux and the current liquid ice freezing and melting mass flux.

[0059] In this embodiment, the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid are input into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function, so as to determine the current vortex vector based on the current displacement vector, the current velocity vector and the curl operator of the fluid, and update the current vortex vector using the apparent kinematic viscosity and the surface tension vector of the fluid; then the stream function vector is determined based on the current displacement vector and the updated current vortex vector, and the velocity potential function is determined based on the current displacement vector, the current velocity vector and the current mass flux, so as to update the current velocity vector using the stream function vector and the velocity potential function, and predict the new current gas-liquid phase function using the current displacement vector, the current mass flux and the updated current velocity vector. It should be noted that the current mass flux includes the current gas-liquid evaporation and condensation mass flux and the current liquid ice freezing and melting mass flux.

[0060] Among them, the coupling relationship between the velocity vector, vortex vector and gas-liquid phase function includes:

[0061] ;

[0062] ;

[0063] ;

[0064] ;

[0065] ;

[0066] The coupling relation adopts the Einstein summation convention, and Indicates three different directions. Indicates three different directions. Indicates three different directions. are vortex vectors, represents the curl operator, are velocity vectors, All represent displacement vectors, t represents time, represents the apparent kinematic viscosity, represents the surface tension vector, denotes the stream function vector, represents the velocity potential function, represents the density of saturated water vapor, represents the density of the liquid, represents the density of ice, represents the gas-liquid evaporation condensation mass flux, represents the mass flux of liquid ice freezing and melting, It should be noted that all vectors in this application are expressed in the form of Einstein indicators.

[0067] by Take the example of Einstein's summation convention:

[0068] .

[0069] In addition, this embodiment uses the Runge-Kutta third-order format (Runge-Kutta third-order format) coupled with the Crank-Nicolson second-order format (Crank-Nicholson second-order format) to numerically discretize the unsteady time term The nonlinear convection term is numerically discretized using a hybrid format of sixth-order compact differential coupling WENO5 (Weighted Essentially Non-Oscillatory) and convection splitting technology. , the dissipative term is numerically discretized using the second-order central difference scheme and vortex stretching term , using the second-order central difference format to discretize the curl term , the second-order derivative and Both use the second-order central difference format to discretize the curl term. .

[0070] More specifically, this embodiment inputs the current gas-liquid phase function of the fluid, the current gas-liquid evaporation and condensation mass flux, the current liquid ice freezing and melting mass flux, the current velocity vector and the current displacement vector into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to determine the current vortex vector based on the current displacement vector, the current velocity vector and the curl operator of the fluid, and then uses the apparent kinematic viscosity and surface tension vector of the fluid, the current displacement vector, the current velocity vector, the curl operator and the vortex vector at the current moment to determine the vortex vector at the next moment, so as to update the current vortex vector to obtain an updated current vortex vector; wherein the updated current vortex vector is the vortex vector at the next moment. Furthermore, the stream function vector is determined based on the current displacement vector and the updated current vortex vector, and the velocity potential function is determined based on the current displacement vector, the current velocity vector, the current gas-liquid evaporation-condensation mass flux, the current liquid ice freezing and melting mass flux, and the densities of saturated water vapor, liquid and ice. The current velocity vector is then updated using the stream function vector, the velocity potential function, the curl operator and the current displacement vector to obtain the updated current velocity vector. Finally, the new current gas-liquid phase function is predicted using the current displacement vector, the current gas-liquid evaporation-condensation mass flux, the current liquid ice freezing and melting mass flux, the density of the liquid and the updated current velocity vector.

[0071] Before using the apparent kinematic viscosity and surface tension vector of the fluid to update the current vortex vector, it is necessary to first determine the apparent kinematic viscosity and surface tension vector of the fluid. Among them, for the determination of the apparent kinematic viscosity of the fluid, the apparent dynamic viscosity of the fluid can be determined based on the dynamic viscosity of the liquid and saturated water vapor and the current gas-liquid phase function, and the apparent density of the fluid can be determined based on the density of the liquid and saturated water vapor and the current gas-liquid phase function, and then the apparent kinematic viscosity of the fluid can be determined based on the ratio of the apparent dynamic viscosity and the apparent density of the fluid. The formula involved is as follows:

[0072] ;

[0073] represents the apparent density, represents the gas-liquid phase function, represents the density of the liquid, represents the density of saturated water vapor, represents the apparent dynamic viscosity, represents the dynamic viscosity of the liquid, represents the dynamic viscosity of saturated water vapor, Indicates apparent kinematic viscosity.

[0074] To determine the surface tension vector of the fluid, the gas-liquid interface normal vector can be determined based on the current displacement vector and the current gas-liquid phase function, and the curvature coefficient can be determined based on the gas-liquid interface normal vector and the current displacement vector. Then, the surface tension vector of the fluid can be determined based on the gas-liquid interface normal vector, the curvature coefficient, the current displacement vector and the current gas-liquid phase function. The formulas involved are as follows:

[0075] ;

[0076] ;

[0077] The above formulas all use the Einstein summation convention, and Indicates three different directions. represents the surface tension vector, represents the surface tension coefficient, represents the curvature coefficient, represents the displacement vector, represents the gas-liquid phase function, represents the normal vector of the air-liquid interface.

[0078] Step S13, using the current liquid ice phase function, the current mass flux, the updated current velocity vector and the predicted current gas-liquid phase function to update the current temperature of the fluid, and based on the updated current temperature, updating the current mass flux, the current liquid ice phase function and the current gas-liquid phase function to determine the phase state at any position in the calculation domain and analyze the dynamic freezing process.

[0079] In this embodiment, the apparent heat capacity, apparent thermal conductivity and apparent density are first determined based on the current liquid ice phase function and the predicted current gas-liquid phase function, thereby determining the corresponding apparent thermal diffusion coefficient based on the apparent density, apparent heat capacity and apparent thermal conductivity, and then the current temperature of the fluid is updated using the apparent thermal diffusion coefficient, the current mass flux and the updated current velocity vector to obtain the updated current temperature. Further, the current mass flux, the current liquid ice phase function and the current gas-liquid phase function are updated based on the apparent thermal conductivity and the updated current temperature to obtain the updated current mass flux, the updated current liquid ice phase function and the updated current gas-liquid phase function, thereby determining the phase state at any position in the calculation domain based on the updated current liquid ice phase function and the updated current gas-liquid phase function, and realizing the analysis of the dynamic freezing process.

[0080] For updating the current temperature, the corresponding apparent heat capacity can be determined based on the heat capacity of the liquid and saturated water vapor and the predicted current gas-liquid phase function, and the corresponding apparent thermal conductivity can be determined based on the thermal conductivity of the liquid, saturated water vapor and ice, the current liquid-ice phase function and the predicted current gas-liquid phase function, and the apparent density can be determined based on the density of the liquid and saturated water vapor and the predicted current gas-liquid phase function, and then the corresponding apparent thermal diffusion coefficient can be determined based on the apparent density, apparent heat capacity and apparent thermal conductivity, and finally the current displacement vector, apparent thermal diffusion coefficient, current gas-liquid evaporation and condensation mass flux, current liquid-ice freezing and melting mass flux and updated current velocity vector of the fluid can be used to determine the temperature at the next moment based on the current moment temperature, so as to update the current temperature of the fluid and obtain the updated current temperature. It should be noted that the updated current temperature is the temperature at the next moment.

[0081] Among them, the formula for determining the apparent thermal diffusion coefficient includes: ; is the apparent thermal diffusion coefficient, represents the apparent thermal conductivity, represents the apparent heat capacity, represents the apparent density. The formula for determining is: ; is the heat capacity of the liquid, is the heat capacity of saturated water vapor, Represents the gas-liquid phase function. Apparent thermal conductivity The formula for determining is: ; represents the thermal conductivity of the liquid, represents the thermal conductivity of saturated water vapor, is the thermal conductivity of ice, represents the gas-liquid phase function, represents the liquid-ice phase function.

[0082] Among them, the temperature update formula includes: ; T represents temperature, t represents time, represents the displacement vector, represents the velocity vector, is the apparent thermal diffusion coefficient, represents the gas-liquid evaporation condensation mass flux, represents the mass flux of liquid ice freezing and melting, represents the density of the liquid; and the updated formula for temperature uses the Einstein summation convention, Indicates three different directions.

[0083] For the update of mass flux and phase function, the current gas-liquid evaporation-condensation mass flux can be updated based on the apparent thermal conductivity, the gas-liquid phase change latent heat and the updated current temperature, and the current liquid ice freezing-melting mass flux can be updated based on the apparent thermal conductivity, the liquid ice phase change latent heat and the updated current temperature. Then, the current liquid ice phase function is updated based on the heat capacity and density of ice and the updated current temperature, and the current gas-liquid phase function is updated based on the heat capacity and density of ice, the density of saturated water vapor and the updated current temperature.

[0084] Among them, the update formulas for mass flux and phase function include:

[0085] ; ;

[0086] ; ;

[0087] represents the gas-liquid evaporation condensation mass flux, represents the latent heat of gas-liquid phase change, represents the apparent thermal conductivity, T represents the temperature, represents the displacement vector, represents the mass flux of liquid ice freezing and melting, is the latent heat of liquid-to-ice phase change, represents the density of ice, is the heat capacity of ice, represents the liquid-ice phase function, t represents time, represents the Laplace operator, represents the gas-liquid phase function, represents the density of saturated water vapor; and the updated formulas for mass flux and phase function use the Einstein summation convention, Indicates three different directions.

[0088] Specifically, for the update of the liquid ice phase function, the heat capacity and density of ice, the apparent thermal conductivity, and the updated current temperature can be used to determine the liquid ice phase function at the next moment based on the liquid ice phase function at the current moment, so as to update the current liquid ice phase function and obtain the updated current liquid ice phase function; wherein the updated current liquid ice phase function is the liquid ice phase function at the next moment. For the update of the gas-liquid phase function, the heat capacity and density of ice, the density of saturated water vapor, the apparent thermal conductivity, and the updated current temperature can be used to determine the gas-liquid phase function at the next moment based on the gas-liquid phase function at the current moment, so as to update the current gas-liquid phase function and obtain the updated current gas-liquid phase function; wherein the updated current gas-liquid phase function is the gas-liquid phase function at the next moment.

[0089] It can be seen that the present application gradually solves the vortex vector, velocity vector, gas-liquid phase function, temperature and mass flux at each moment to determine the gas-liquid phase function and liquid-ice phase function at each moment, thereby determining the phase state at any position in the calculation domain according to the gas-liquid phase function and liquid-ice phase function, and then realizing autonomous analysis of the dynamic freezing process. It not only fills the gap in the field of gas-liquid-ice three-phase flow analysis, realizes the autonomous control of autonomous source code software, but also helps to better study the mutual transformation between the three phases of gas, liquid and ice, and provides strong support for the analysis of freezing problems.

[0090] See also Figure 2 As shown, the embodiment of the present invention discloses a gas-liquid-ice three-phase flow analysis device, comprising:

[0091] A first updating module 11 is used to determine, in a calculation domain predetermined based on industrial equipment, a velocity vector of a fluid at a boundary in the calculation domain based on a current liquid-ice phase function, so as to update a current velocity vector of the fluid;

[0092] The second updating module 12 is used to input the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to update the current vortex vector and the current velocity vector of the fluid and predict the new current gas-liquid phase function; wherein the current mass flux includes the current gas-liquid evaporation and condensation mass flux and the current liquid ice freezing and melting mass flux;

[0093] The phase state determination module 13 is used to update the current temperature of the fluid using the current liquid ice phase function, the current mass flux, the updated current velocity vector and the predicted current gas-liquid phase function, and to update the current mass flux, the current liquid ice phase function and the current gas-liquid phase function based on the updated current temperature to determine the phase state at any position in the calculation domain and analyze the dynamic freezing process.

[0094] It can be seen that the present application gradually solves the vortex vector, velocity vector, gas-liquid phase function, temperature and mass flux at each moment to determine the gas-liquid phase function and liquid-ice phase function at each moment, thereby determining the phase state at any position in the calculation domain according to the gas-liquid phase function and liquid-ice phase function, and then realizing autonomous analysis of the dynamic freezing process. It not only fills the gap in the field of gas-liquid-ice three-phase flow analysis, realizes the autonomous control of autonomous source code software, but also helps to better study the mutual transformation between the three phases of gas, liquid and ice, and provides strong support for the analysis of freezing problems.

[0095] In some specific embodiments, the first updating module 11 includes:

[0096] A first vector determination unit is used to determine the velocity vector of the fluid at the boundary in the calculation domain to be zero when the current liquid ice phase function characterizes the current phase state as an ice state;

[0097] The second vector determination unit is used to determine the velocity vector of the fluid at the boundary of the calculation domain based on the current liquid ice phase function and the current velocity vector of the fluid when the current liquid ice phase function indicates that the current phase state is not an ice state.

[0098] In some specific embodiments, the second updating module 12 includes:

[0099] a first vector updating unit, for inputting the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into a coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function, so as to determine the current vortex vector based on the current displacement vector, the current velocity vector and the curl operator of the fluid, and to update the current vortex vector using the apparent kinematic viscosity and the surface tension vector of the fluid;

[0100] The second vector updating unit is used to determine the flow function vector based on the current displacement vector and the updated current vortex vector, and determine the velocity potential function based on the current displacement vector, the current velocity vector and the current mass flux, so as to update the current velocity vector using the flow function vector and the velocity potential function, and predict a new current gas-liquid phase function using the current displacement vector, the current mass flux and the updated current velocity vector.

[0101] In some specific embodiments, the gas-liquid-ice three-phase flow analysis device further includes:

[0102] a kinematic viscosity determination unit, for determining the apparent dynamic viscosity of the fluid according to the dynamic viscosity of the liquid and the saturated water vapor and the current gas-liquid phase function, and determining the apparent density of the fluid according to the density of the liquid and the saturated water vapor and the current gas-liquid phase function, so as to determine the apparent kinematic viscosity of the fluid based on the ratio of the apparent dynamic viscosity to the apparent density;

[0103] A tension vector determination unit is used to determine the gas-liquid interface normal vector according to the current displacement vector and the current gas-liquid phase function, and to determine the curvature coefficient according to the gas-liquid interface normal vector and the current displacement vector, so as to determine the surface tension vector based on the gas-liquid interface normal vector, the curvature coefficient, the current displacement vector and the current gas-liquid phase function.

[0104] In some specific embodiments, the phase determination module 13 includes:

[0105] a first coefficient determination unit, for determining corresponding apparent heat capacities according to the heat capacities of the liquid and saturated water vapor and the predicted current gas-liquid phase function, and determining corresponding apparent thermal conductivity coefficients according to the thermal conductivity coefficients of the liquid, saturated water vapor and ice and the current liquid-ice phase function and the predicted current gas-liquid phase function;

[0106] A second coefficient determination unit, configured to determine a corresponding apparent thermal diffusion coefficient based on an apparent density, the apparent heat capacity and the apparent thermal conductivity; the apparent density is a density determined according to the density of the liquid and the saturated water vapor and a predicted current gas-liquid phase function;

[0107] a temperature updating unit, configured to update a current temperature of the fluid using a current displacement vector of the fluid, the apparent thermal diffusion coefficient, a current mass flux and an updated current velocity vector;

[0108] Among them, the temperature update formula includes: ; T represents temperature, t represents time, represents the displacement vector, represents the velocity vector, is the apparent thermal diffusion coefficient, represents the gas-liquid evaporation condensation mass flux, represents the mass flux of liquid ice freezing and melting, represents the density of the liquid; and the updated formula for temperature uses the Einstein summation convention, Indicates three different directions.

[0109] In some specific embodiments, the phase determination module 13 includes:

[0110] A mass flux updating unit, used to update the current gas-liquid evaporation-condensation mass flux based on the gas-liquid phase change latent heat and the updated current temperature, and to update the current liquid ice freezing-melting mass flux based on the liquid ice phase change latent heat and the updated current temperature;

[0111] a phase function updating unit, configured to update a current liquid-ice phase function based on the heat capacity and density of ice and the updated current temperature, and to update a current gas-liquid phase function based on the heat capacity and density of ice, the density of saturated water vapor, and the updated current temperature;

[0112] Among them, the update formulas for mass flux and phase function include:

[0113] ; ;

[0114] ; ;

[0115] represents the gas-liquid evaporation condensation mass flux, represents the latent heat of gas-liquid phase change, represents the apparent thermal conductivity, T represents the temperature, represents the displacement vector, represents the mass flux of liquid ice freezing and melting, is the latent heat of liquid-to-ice phase change, represents the density of ice, is the heat capacity of ice, represents the liquid-ice phase function, t represents time, represents the Laplace operator, represents the gas-liquid phase function, represents the density of saturated water vapor; and the updated formulas for mass flux and phase function use the Einstein summation convention, Indicates three different directions.

[0116] Furthermore, the present application also discloses an electronic device. Figure 3 This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content in the diagram cannot be regarded as any limitation on the scope of use of the present application.

[0117] Figure 3 A schematic diagram of the structure of an electronic device 20 provided in an embodiment of the present application. The electronic device 20 may specifically include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is used to store a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the gas-liquid-ice three-phase flow analysis method disclosed in any of the aforementioned embodiments. In addition, the electronic device 20 in this embodiment may specifically be an electronic computer.

[0118] In this embodiment, the power supply 23 is used to provide working voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and the external device, and the communication protocol it follows is any communication protocol that can be applied to the technical solution of the present application, and is not specifically limited here; the input and output interface 25 is used to obtain external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs and is not specifically limited here.

[0119] In addition, the memory 22, as a carrier for storing resources, can be a read-only memory, a random access memory, a disk or an optical disk, etc. The resources stored thereon can include an operating system 221, a computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0120] The operating system 221 is used to manage and control the hardware devices and computer program 222 on the electronic device 20, and can be Windows Server, Netware, Unix, Linux, etc. In addition to the computer program that can be used to complete the gas-liquid-ice three-phase flow analysis method performed by the electronic device 20 disclosed in any of the aforementioned embodiments, the computer program 222 can further include computer programs that can be used to complete other specific tasks.

[0121] Furthermore, the present application also discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, the gas-liquid-ice three-phase flow analysis method disclosed above is implemented. For the specific steps of the method, reference may be made to the corresponding contents disclosed in the above embodiments, and no further description will be given here.

[0122] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0123] 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 composition and steps of each example have been generally described in the above description according to function. 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 to be beyond the scope of this application.

[0124] The steps of the method or algorithm 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), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

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

[0126] The technical solution provided by the present application is introduced in detail above. Specific examples are used in this article to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method of the present application and its core idea. At the same time, for general technicians in this field, according to the idea of ​​the present application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present application.

Claims

1. A gas-liquid-ice three-phase flow analysis method, characterized in that: include: In a calculation domain predetermined based on industrial equipment, a velocity vector of a fluid at a boundary in the calculation domain is determined based on a current liquid-ice phase function to update a current velocity vector of the fluid; Inputting the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to update the current vortex vector and the current velocity vector of the fluid and predict a new current gas-liquid phase function; wherein the current mass flux includes the current gas-liquid evaporation and condensation mass flux and the current liquid ice freezing and melting mass flux; updating the current temperature of the fluid using the current liquid ice phase function, the current mass flux, the updated current velocity vector and the predicted current gas-liquid phase function, and updating the current mass flux, the current liquid ice phase function and the current gas-liquid phase function based on the updated current temperature to determine the phase state at any position in the computational domain and analyze the dynamic freezing process; The step of inputting the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to update the current vortex vector and the current velocity vector of the fluid and predict a new current gas-liquid phase function comprises: Inputting the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function, so as to determine the current vortex vector based on the current displacement vector, the current velocity vector and the curl operator of the fluid, and updating the current vortex vector using the apparent kinematic viscosity and the surface tension vector of the fluid; Determine a stream function vector based on a current displacement vector and an updated current vortex vector, and determine a velocity potential function based on the current displacement vector, the current velocity vector, and the current mass flux, so as to update the current velocity vector using the stream function vector and the velocity potential function, and predict a new current gas-liquid phase function using the current displacement vector, the current mass flux, and the updated current velocity vector; The coupling relationship includes: ； ； ； ； ； Wherein, the coupling relation adopts the Einstein summation convention, and Indicates three different directions. Indicates three different directions. Indicates three different directions. are vortex vectors, represents the curl operator, are velocity vectors, All represent displacement vectors, t represents time, represents the apparent kinematic viscosity, represents the surface tension vector, denotes the stream function vector, represents the velocity potential function, represents the density of saturated water vapor, represents the density of the liquid, represents the density of ice, represents the gas-liquid evaporation condensation mass flux, represents the mass flux of liquid ice freezing and melting, represents the gas-liquid phase function; Wherein, updating the current temperature of the fluid using the current liquid-ice phase function, the current mass flux, the updated current velocity vector and the predicted current gas-liquid phase function comprises: Determine the corresponding apparent heat capacity according to the heat capacity of the liquid and saturated water vapor and the predicted current gas-liquid phase function, and determine the corresponding apparent thermal conductivity according to the thermal conductivity of the liquid, saturated water vapor and ice and the current liquid-ice phase function and the predicted current gas-liquid phase function; Determining a corresponding apparent thermal diffusion coefficient based on the apparent density, the apparent heat capacity and the apparent thermal conductivity; the apparent density is a density determined according to the density of the liquid and the saturated water vapor and the predicted current gas-liquid phase function; updating a current temperature of the fluid using a current displacement vector of the fluid, the apparent thermal diffusivity, a current mass flux, and an updated current velocity vector; Among them, the temperature update formula includes: ; T represents temperature, t represents time, represents the displacement vector, represents the velocity vector, is the apparent thermal diffusion coefficient, represents the gas-liquid evaporation condensation mass flux, represents the mass flux of liquid ice freezing and melting, represents the density of the liquid; and the updated formula for temperature uses the Einstein summation convention, Indicates three different directions.

2. The gas-liquid-ice three-phase flow analysis method according to claim 1, characterized in that: The determining of the velocity vector of the fluid at the boundary in the computational domain based on the current liquid-ice phase function comprises: When the current liquid-ice phase function characterizes that the current phase state is an ice state, determining the velocity vector of the fluid at the boundary in the calculation domain to be zero; When the current liquid-ice phase function indicates that the current phase state is not an ice state, a velocity vector of the fluid at the boundary in the calculation domain is determined based on the current liquid-ice phase function and the current velocity vector of the fluid.

3. The gas-liquid-ice three-phase flow analysis method according to claim 1, characterized in that: Before the current vortex vector is updated by using the apparent kinematic viscosity and surface tension vector of the fluid, the method further includes: Determining the apparent dynamic viscosity of the fluid according to the dynamic viscosity of the liquid and the saturated water vapor and the current gas-liquid phase function, and determining the apparent density of the fluid according to the density of the liquid and the saturated water vapor and the current gas-liquid phase function, so as to determine the apparent kinematic viscosity of the fluid based on the ratio of the apparent dynamic viscosity to the apparent density; The gas-liquid interface normal vector is determined according to the current displacement vector and the current gas-liquid phase function, and the curvature coefficient is determined according to the gas-liquid interface normal vector and the current displacement vector, so as to determine the surface tension vector based on the gas-liquid interface normal vector, the curvature coefficient, the current displacement vector and the current gas-liquid phase function.

4. The gas-liquid-ice three-phase flow analysis method according to claim 1, characterized in that: The updating of the current mass flux, the current liquid-ice phase function and the current gas-liquid phase function based on the updated current temperature includes: updating the current gas-liquid evaporation-condensation mass flux based on the gas-liquid phase change latent heat and the updated current temperature, and updating the current liquid ice freezing-melting mass flux based on the liquid ice phase change latent heat and the updated current temperature; updating the current liquid-ice phase function based on the heat capacity and density of ice and the updated current temperature, and updating the current gas-liquid phase function based on the heat capacity and density of ice, the density of saturated water vapor, and the updated current temperature; Among them, the update formulas for mass flux and phase function include: ； ； ； ； represents the gas-liquid evaporation condensation mass flux, represents the latent heat of gas-liquid phase change, represents the apparent thermal conductivity, T represents the temperature, represents the displacement vector, represents the mass flux of liquid ice freezing and melting, is the latent heat of liquid-to-ice phase change, represents the density of ice, is the heat capacity of ice, represents the liquid-ice phase function, t represents time, represents the Laplace operator, represents the gas-liquid phase function, represents the density of saturated water vapor; and the updated formulas for mass flux and phase function use the Einstein summation convention, Indicates three different directions.

5. A gas-liquid-ice three-phase flow analysis device, characterized in that: The gas-liquid-ice three-phase flow analysis device is used to implement the gas-liquid-ice three-phase flow analysis method according to any one of claims 1 to 4, and the gas-liquid-ice three-phase flow analysis device comprises: A first updating module is used to determine, in a calculation domain predetermined based on industrial equipment, a velocity vector of a fluid at a boundary in the calculation domain based on a current liquid-ice phase function, so as to update a current velocity vector of the fluid; A second updating module is used to input the current gas-liquid phase function, the current mass flux and the updated current velocity vector of the fluid into the coupling relationship between the velocity vector, the vortex vector and the gas-liquid phase function to update the current vortex vector and the current velocity vector of the fluid and predict the new current gas-liquid phase function; wherein the current mass flux includes the current gas-liquid evaporation and condensation mass flux and the current liquid ice freezing and melting mass flux; A phase state determination module is used to update the current temperature of the fluid using the current liquid ice phase function, the current mass flux, the updated current velocity vector and the predicted current gas-liquid phase function, and to update the current mass flux, the current liquid ice phase function and the current gas-liquid phase function based on the updated current temperature to determine the phase state at any position in the calculation domain and analyze the dynamic freezing process.

6. An electronic device, characterized in that: include: Memory, used to store computer programs; A processor, configured to execute the computer program to implement the gas-liquid-ice three-phase flow analysis method according to any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: Used to store a computer program, which, when executed by a processor, implements the gas-liquid-ice three-phase flow analysis method according to any one of claims 1 to 4.

Images