Subcooled boiling and boiling critical value numerical simulation method, device, medium and equipment

Abstract

The application discloses a subcooled boiling and boiling critical numerical simulation method, device, medium and equipment, and relates to the technical field of heat transfer. The method comprises the following steps: a mixture model is established, and an inertial acceleration volume force generated in a non-inertial system is added to a momentum conservation equation of the mixture model to obtain a non-inertial system mixture model; a bubble dynamics model is established according to a plurality of acting forces acting on a wall surface bubble and the inertial acceleration volume force, and a plurality of bubble characteristic parameters are obtained according to the bubble dynamics model; a wall surface heat flow distribution model is established, and the bubble characteristic parameters are substituted into the wall surface heat flow distribution model to obtain a wall surface heat flow density model; and the wall surface heat flow density model and the non-inertial system mixture model are coupled to obtain a two-phase boiling heat exchange and boiling critical calculation model, which is used for simulating a boiling heat exchange working condition under a motion condition. The numerical simulation method has the advantages of less calculation resources, wide application range, good convergence, and high prediction accuracy.

Description

Technical Field

[0001] The present application relates to the field of heat transfer technology, and in particular to a method, device, medium and equipment for numerical simulation of subcooled boiling and boiling criticality. Background Art

[0002] In high-end industrial fields such as spacecraft thermal management systems, superconducting magnet cooling devices, and liquid cooling of high-power electronic devices, accurate simulation of boiling heat transfer processes is crucial to equipment safety and efficiency. For example, superconducting magnets need to maintain stable operation at extremely low temperatures, and the boiling behavior of liquid helium in their cooling channels directly affects magnet performance. Spacecraft thermal management systems need to predict phase change heat transfer characteristics in microgravity and variable acceleration environments to ensure thermal protection reliability. Currently, traditional numerical simulation methods have become a key bottleneck restricting technological breakthroughs in related fields because they cannot balance computational efficiency and multi-physics field coupling accuracy.

[0003] Current commercial CFD (Computational Fluid Dynamics) platforms often use an Euler-Euler two-fluid model coupled with a wall heat flux distribution model for boiling simulation. While this model can describe complex gas-liquid two-phase flows by solving six sets of coupled conservation equations and interphase forces, its high degree of equation coupling and computational complexity lead to high model complexity and poor convergence. Furthermore, the bubble detachment characteristic submodel in traditional models relies on empirical relationships, limiting its versatility to the scope of experimental data. This makes it difficult to adapt to new working fluids or changes in physical parameters under dynamic conditions, resulting in a limited scope of applicability. Summary of the Invention

[0004] In view of this, the embodiments of the present application provide a method, device, medium and equipment for numerical simulation of subcooled boiling and boiling criticality, the main purpose of which is to solve the technical problems of high computational overhead, poor model convergence and small scope of application in simulating boiling heat transfer conditions.

[0005] According to one aspect of the present application, a method for numerical simulation of subcooled boiling and boiling criticality is provided, the method comprising:

[0006] Establishing a mixture model, and adding the inertial acceleration body force generated in the non-inertial system to the momentum conservation equation of the mixture model to obtain the non-inertial system mixture model;

[0007] Establishing a bubble dynamics model based on multiple forces acting on the wall bubbles and the inertial acceleration volume force, and obtaining multiple bubble characteristic parameters based on the bubble dynamics model;

[0008] The wall heat flux density model is obtained by establishing a wall heat flux distribution model and substituting the bubble characteristic parameters into the wall heat flux distribution model.

[0009] The two-phase boiling heat exchange and boiling critical calculation model is obtained by coupling the wall heat flux density model and the non-inertial system mixture model, and is used for simulating the boiling heat exchange working condition under the motion condition.

[0010] According to another aspect of the present application, a subcooled boiling and boiling critical numerical simulation device is provided, which comprises:

[0011] The mixture model establishing module is configured to establish a Mixture mixture model, add the inertial acceleration volume force generated in the non-inertial system to a momentum conservation equation of the Mixture mixture model, and obtain a non-inertial system mixture model.

[0012] The bubble dynamics model establishing module is configured to establish a bubble dynamics model according to a plurality of acting forces acting on the wall bubble and the inertial acceleration volume force, and obtain a plurality of bubble characteristic parameters according to the bubble dynamics model.

[0013] The wall heat flux density model establishing module is configured to establish a wall heat flux distribution model, and substitute the bubble characteristic parameters into the wall heat flux distribution model to obtain a wall heat flux density model.

[0014] The boiling heat exchange working condition simulation module is configured to couple the wall heat flux density model and the non-inertial system mixture model to obtain a two-phase boiling heat exchange and boiling critical calculation model, and simulate the boiling heat exchange working condition under the motion condition.

[0015] According to still another aspect of the present application, a storage medium having a computer program stored thereon is provided, and the program is executed by a processor to implement the above-mentioned subcooled boiling and boiling critical numerical simulation method.

[0016] According to still another aspect of the present application, a computer device is provided, which comprises a storage medium, a processor, and a computer program stored in the storage medium and executable on the processor, and the processor executes the program to implement the above-mentioned subcooled boiling and boiling critical numerical simulation method.

[0017] By means of the above technical solution, the embodiment of the present application provides a method and apparatus, medium and equipment for numerical simulation of supercooled boiling and boiling criticality. By adopting the Mixture mixture model as the model framework, it is only necessary to independently solve three conservation equations and one cavitation fraction equation. The model is simple and has good convergence, which can effectively reduce the computational overhead of the model. At the same time, by adopting the bubble dynamics model instead of the bubble detachment model of the empirical relationship, and adding the inertial acceleration volume force generated by the motion to the energy conservation equation and the bubble dynamics model, the adaptability of the model to new working fluids and complex motion conditions can be significantly improved. In addition, by coupling the bubble characteristic parameters with the thermodynamic equations, and coupling the wall heat flux density model with the non-inertial system mixture model, a multi-scale coupling framework is obtained, which can realize two-way feedback between bubble dynamics and macroscopic flow field, thereby significantly improving the prediction accuracy of boiling heat transfer and boiling critical conditions.

[0018] The above description is only an overview of the technical solution of the present application. In order to more clearly understand the technical means of the present application, it can be implemented in accordance with the contents of the specification. In order to make the above and other purposes, features and advantages of the present application more obvious and easy to understand, the specific implementation methods of the present application are listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] The drawings described herein are used to provide a further understanding of the present application and constitute a part of the present application. The illustrative embodiments of the present application and their descriptions are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0020] Figure 1 A schematic flow chart of a method for numerically simulating subcooled boiling and boiling criticality provided in an embodiment of the present application is shown;

[0021] Figure 2 A schematic diagram of the force acting on a static bubble provided in an embodiment of the present application is shown;

[0022] Figure 3 A schematic diagram of bubble sliding provided in an embodiment of the present application is shown;

[0023] Figure 4 A schematic diagram of a process for calculating bubble characteristic parameters provided by an embodiment of the present application is shown;

[0024] Figure 5 A schematic diagram of a process for solving a two-phase boiling heat transfer and boiling criticality calculation model provided in an embodiment of the present application is shown;

[0025] Figure 6 A schematic structural diagram of a subcooled boiling and boiling critical numerical simulation device provided in an embodiment of the present application is shown. DETAILED DESCRIPTION

[0026] The present application will be described in detail below with reference to the accompanying drawings and in combination with embodiments. It should be noted that, unless there is a conflict, the embodiments and features in the embodiments of the present application can be combined with each other.

[0027] Currently, most boiling models available on commercial CFD platforms primarily utilize a two-fluid Euler-Euler multiphase flow model coupled with a wall heat flux distribution model (Rensselaer Polytechnic Institute (RPI)). However, the Euler-Euler two-fluid model requires solving six coupled conservation equations, coupled with multiple interphase forces, making the model complex and prone to poor convergence. Furthermore, the traditional RPI bubble detachment submodel relies heavily on empirical relationships, but these limitations restrict its applicability. For example, a submodel applicable to water is no longer applicable to liquid nitrogen. Furthermore, the lack of corresponding experimental relationships for novel working fluids limits the model's applicability and its mechanics are weak. It cannot comprehensively account for the effects of wall inclination, pressure, and other wall conditions. In particular, under motion conditions, motion generates inertial acceleration, which alters the forces acting on the bubbles, affecting their dynamic behavior and, consequently, boiling heat transfer. These empirical relationships cannot be accounted for in the bubble submodel.

[0028] In order to solve the above problems, in one embodiment, Figure 1 As shown, a numerical simulation method for subcooled boiling and boiling criticality is provided, which is described by taking the application of the method to a computer device as an example, and includes the following steps:

[0029] Step 101: Establish a mixture model, and add the inertial acceleration body force generated in the non-inertial system to the momentum conservation equation of the mixture model to obtain the non-inertial system mixture model.

[0030] The Mixture model is a multiphase flow model that equates gas-liquid two-phase flow to a single mixed working fluid. This model introduces the cavity fraction equation to describe the volume fractions of the two phases. Through three sets of conservation equations (mass conservation equation, momentum conservation equation, and energy conservation equation) and a cavity fraction equation, it effectively simplifies the model's computational complexity. Furthermore, inertial acceleration body forces refer to additional forces generated by motion in a non-inertial frame, such as centrifugal force and Coriolis force. These inertial acceleration body forces can act on the momentum conservation equation of the Mixture model through the body force source term.

[0031] Specifically, the conservation equations of the Mixture mixture model are first established, including the mass conservation equation, momentum conservation equation, energy conservation equation and cavitation fraction equation. Then, according to the motion characteristics of the fluid pipeline, the additional inertial acceleration force equation generated in the non-inertial system is established to calculate the inertial acceleration volume force generated by the motion of the object in the non-inertial system. The inertial acceleration volume force is added to the momentum conservation equation in the Mixture mixture model as a volume force source term to obtain the non-inertial system mixture model.

[0032] This embodiment uses the Mixture model as a framework to equate the gas-liquid two-phase working fluid to a mixture working fluid, so that only three conservation equations and one cavitation fraction equation need to be solved independently. The model is simple and has good convergence, which can effectively reduce computational overhead.

[0033] Step 102 : establishing a bubble dynamics model based on multiple forces acting on the wall bubbles and the inertial acceleration body force, and obtaining multiple bubble characteristic parameters based on the bubble dynamics model.

[0034] Among them, the bubble dynamics model refers to a microscopic mathematical model established based on the bubble force analysis mechanism. This model can be used to predict the generation, growth and detachment behavior of bubbles.

[0035] Specifically, a model for wall bubble growth and the various forces acting on it is first established, including the bubble growth equation, surface tension equation, bubble shear lift equation, bubble growth force equation, bubble drag equation, bubble buoyancy equation, bubble contact pressure equation, bubble dynamic pressure equation, and virtual mass force equation. Then, the inertial acceleration volume force calculated in step 101 is applied to the bubble, thereby accounting for the influence of motion conditions on bubble dynamics. Furthermore, the forces acting on the bubble are decomposed along the flow direction (i.e., streamwise) and the direction perpendicular to the flow (i.e., normal), and a bubble dynamics model is established based on the various forces. Finally, key parameters such as bubble detachment frequency, bubble detachment diameter, slip distance, bubble rise diameter, and bubble nucleation density can be obtained based on the bubble dynamics model.

[0036] This embodiment adopts a bubble dynamics model instead of the bubble detachment model of the empirical relationship, and adds the inertial acceleration volume force generated by the movement to the bubble dynamics model, which can significantly improve the adaptability to new working fluids (such as liquid nitrogen) and complex motion conditions. At the same time, the dynamic force analysis of the bubbles can accurately capture the dynamic behavior of the bubbles, thereby improving the accuracy of the boiling heat transfer simulation.

[0037] Step 103: Establish a wall heat flux distribution model, and substitute the bubble characteristic parameters into the wall heat flux distribution model to obtain a wall heat flux density model.

[0038] The wall heat flux distribution model is a mathematical model that describes the distribution of wall heat flux between the bubble generation zone and the single-phase zone. Its core is to quantify the local heat flux density using bubble characteristic parameters. The wall heat flux density model is a macroscopic heat flux distribution model established by combining bubble behavior with the heat conduction equation.

[0039] Specifically, first, the wall supercooled boiling phenomenon and boiling critical phenomenon are modeled, that is, a wall heat flux distribution model is established, which can specifically include models such as quenching heat flux, liquid relative heat flux, evaporation heat flux, slip heat flux and vapor relative heat flux. Then, the bubble characteristic parameters such as bubble detachment frequency, bubble detachment diameter, slip distance, bubble rising diameter and bubble nucleation density calculated in step 102 can be introduced into the quenching heat flux, liquid relative heat flux, evaporation heat flux, slip heat flux and vapor relative heat flux models, so as to accurately model the wall boiling phenomenon and boiling criticality and obtain the wall heat flux density model under supercooled boiling and critical boiling environments.

[0040] This embodiment couples the bubble characteristic parameters with the thermodynamic equations to achieve dynamic prediction of the wall heat flux density. At the same time, by adding the inertial acceleration volume force to the bubble dynamics model and the energy conservation equation, the adaptability of the model to motion conditions can be enhanced.

[0041] Step 104 : Couple the wall heat flux model and the non-inertial system mixture model to obtain a two-phase boiling heat transfer and boiling criticality calculation model for simulating boiling heat transfer conditions under motion conditions.

[0042] Among them, the two-phase boiling heat transfer and boiling critical calculation model refers to a multi-scale framework constructed by bidirectionally coupling the wall heat flux model and the non-inertial system mixture model. Its core is to obtain the interaction between bubble behavior and the macroscopic flow field through iterative solution.

[0043] Specifically, by coupling the wall heat flux model established in step 103 with the non-inertial mixture model established in step 101, a cross-scale, widely applicable two-phase boiling heat transfer and boiling criticality calculation model based on bubble dynamics and the mixture model under motion conditions can be obtained. Once established, this model can be applied to boiling heat transfer conditions with different inclination angles, pressures, flow rates, and working fluids. It can also be used to simulate boiling heat transfer conditions under motion conditions, thereby evaluating the model's applicability based on the simulation results.

[0044] This embodiment uses the above-mentioned multi-scale coupling framework to achieve two-way feedback between bubble dynamics and macroscopic flow field, thereby significantly improving the prediction accuracy of boiling heat transfer and boiling critical conditions.

[0045] The above embodiment adopts the Mixture mixture model as the model framework, and only needs to independently solve three conservation equations and one cavitation fraction equation. The model is simple and has good convergence, which can effectively reduce the computational overhead of the model. At the same time, by adopting the bubble dynamics model instead of the bubble detachment model of the empirical relationship, and adding the inertial acceleration volume force generated by the movement to the energy conservation equation and the bubble dynamics model, the adaptability of the model to new working fluids and complex motion conditions can be significantly improved. In addition, by coupling the bubble characteristic parameters with the thermodynamic equations, and coupling the wall heat flux density model with the non-inertial system mixture model, a multi-scale coupling framework is obtained, which can realize two-way feedback between bubble dynamics and macroscopic flow field, thereby significantly improving the prediction accuracy of boiling heat transfer and boiling critical conditions.

[0046] In one embodiment, step 101 can be implemented in the following manner: first, based on the motion characteristics of the fluid in the pipeline, an inertial acceleration body force equation generated in a non-inertial system is established to calculate the inertial acceleration body force generated by the motion of the object; then, the mass conservation equation, momentum conservation equation, energy conservation equation, and cavitation fraction equation of the mixture model are established, and the inertial acceleration body force is added to the momentum conservation equation as a body force source term to obtain a non-inertial system mixture model.

[0047] For example, the inertial acceleration body force equation is shown in Formula 1:

[0048] (Formula 1)

[0049] in, is the inertial force density in three orthogonal directions (N / m 3 ), is the density of the mixture (kg / m 3 ), is the angular velocity (rad / s), is the position vector (m), is the angular acceleration (rad / s 2 ), is the relative speed (m / s), are orthogonal unit vectors in three directions.

[0050] For example, the mass conservation equation is shown in Formula 2:

[0051] (Formula 2)

[0052] in, is the density of the mixture (kg / m 3 ), is the vapor phase density (kg / m 3), is the liquid density (kg / m 3 ), is the mixture velocity (m / s), is the vapor phase velocity (m / s), is the liquid phase velocity (m / s), is the vacuolar fraction, is the volume fraction of the liquid phase.

[0053] For example, the momentum conservation equation is shown in Formula 3:

[0054] (Formula 3)

[0055] in, is the density of the mixture (kg / m 3 ), is the k-phase density (kg / m 3 ), is the viscosity of the mixture (kg / (ms)), is the mixture pressure (Pa), is the acceleration due to gravity (m / s 2 ), is the external body force, which is the body force of inertial acceleration due to motion (N / m 3 ), is the mixture velocity (m / s), is the k-phase relative velocity (m / s), is the k-phase velocity (m / s), is the k-phase cavitation fraction, is the k-phase viscosity (kg / (ms)), is the stress tensor (Pa).

[0056] For example, the energy conservation equation is shown in Formula 4:

[0057] (Formula 4)

[0058] in, is the k-phase density (kg / m 3 ), is the k-phase cavitation fraction, is the energy of the mixture (J / kg), is the mixture pressure (Pa), is the effective thermal conductivity (W / (m K)), is the mixture temperature (K), is the volume heat source (J / m 3 ), where is the latent heat generated by wall boiling and the latent heat released by condensation of bubbles in the mainstream, is the k-phase enthalpy (J / (kgK)), is the k-phase velocity (m / s).

[0059] For example, the cavitation fraction equation is shown in Formula 5:

[0060] (Formula 5)

[0061] in, is the k-phase density (kg / m 3 ), is the k-phase cavitation fraction, is the mixture velocity (m / s), is the k-phase relative velocity (m / s), is the mass source term transferred from the q-phase to the k-phase (kg / m3), is the mass source term transferred from k-phase to q-phase (kg / m 3 ), which is mainly due to the mass exchange caused by wall boiling and the mass exchange caused by mainstream bubble condensation.

[0062] In one embodiment, step 102 can be implemented as follows: First, equations for bubble growth on the wall and the various forces acting on it are established, wherein the force equations include the bubble growth equation, the surface tension equation, the bubble shear lift equation, the bubble growth force equation, the bubble drag equation, the bubble buoyancy equation, the bubble contact pressure equation, the bubble dynamic pressure equation, and the virtual mass force equation. Then, a bubble dynamics model is established based on the bubble growth equations and the various forces acting on it, as well as the inertial acceleration volume force equation corresponding to the inertial acceleration volume force. Finally, based on the bubble dynamics model, multiple bubble characteristic parameters are calculated, wherein the bubble characteristic parameters include bubble detachment frequency, bubble detachment diameter, bubble sliding distance, bubble rise diameter, and bubble nucleation density.

[0063] For example, the bubble growth equation is shown in Formula 6:

[0064] (Formula 6)

[0065] in, is the change of bubble radius with time (m / s), is an empirical constant, is the Jakob number, is the thermal diffusivity (m 2 / s), is time (s), is the liquid density (kg / m 3 ), is the liquid specific heat capacity (J / (kg K)), Wall superheat (K), is the vapor phase density (kg / m 3 ), is the latent heat of vaporization (kJ / kg).

[0066] For example, the surface tension equation is shown in Formula 7:

[0067] (Formula 7)

[0068]

[0069] in, is the surface tension in two orthogonal directions (N), represents the surface tension coefficient (m -1 ), is the front contact angle (degrees), is the rear contact angle (degrees), is the bubble contact diameter (m).

[0070] For example, the bubble shear lift equation is shown in Formula 8:

[0071] (Formula 8)

[0072]

[0073] in, is the bubble shear lift (N), is the liquid density (kg / m 3 ), is the relative velocity between the bubble and the liquid phase (m / s), is the bubble radius, is the bubble Reynolds number, is the liquid phase dynamic viscosity (kg / (ms)), is the bubble diameter (m), is the dimensionless liquid shear gradient.

[0074] Exemplarily, the bubble growth force equation is shown in Formula 9:

[0075] (Formula 9)

[0076] in, is the bubble growth force (N), is the shape factor, is the liquid density (kg / m 3 ), is the bubble radius (m), is the bubble growth rate (m / s), is the distance from the bubble to the wall (m).

[0077] For example, the bubble drag equation is shown in Formula 10:

[0078] (Formula 10)

[0079] in, is the bubble drag force (N), is the liquid density (kg / m 3 ), is the liquid phase velocity (m / s), is the bubble velocity (m / s), is the bubble radius, is the bubble Reynolds number, is the drag coefficient.

[0080] Exemplarily, the bubble buoyancy equation is shown in Formula 11:

[0081] (Formula 11)

[0082] in, is the bubble buoyancy (N), is the bubble radius, is the liquid density (kg / m 3 ), is the vapor phase density (kg / m 3 ), is the acceleration due to gravity (m / s 2 ).

[0083] For example, the bubble contact pressure equation is shown in Formula 12:

[0084] (Formula 12)

[0085] in, is the bubble contact pressure (N), is the bubble contact diameter (m), represents the surface tension coefficient (m -1 ), is the bubble radius.

[0086] For example, the bubble dynamic pressure equation is shown in Formula 13:

[0087] (Formula 13)

[0088] in, is the bubble dynamic pressure (N), is the liquid density (kg / m 3 ), is the relative velocity between the bubble and the liquid phase (m / s), is the bubble contact diameter (m).

[0089] For example, the virtual mass force equation is shown in Formula 14:

[0090] (Formula 14)

[0091] in, is the virtual mass force (N), is the liquid density (kg / m 3 ), is the bubble radius, is the liquid phase velocity (m / s), is the bubble velocity (m / s).

[0092] Furthermore, through the above-mentioned force equations (refer to formula 6 to formula 14) and the inertial acceleration volume force equation in step 101 (refer to formula 1), the force acting on the bubble can be obtained, thereby establishing a bubble dynamics model. Figure 2 , the force on the bubble can be decomposed along the flow direction (i.e., streamwise) and the vertical direction (i.e., normal direction). For example, the bubble dynamics model is shown in Formula 15:

[0093] (Formula 15)

[0094] in, is the bubble growth in two orthogonal directions and the various forces it is subjected to, is the surface tension in two orthogonal directions (N), is the bubble growth force in two orthogonal directions (N), is the inertial acceleration body force in two orthogonal directions, is the bubble shear lift (N), is the bubble dynamic pressure (N), is the bubble contact pressure (N), is the bubble buoyancy (N), is the bubble drag force (N).

[0095] Furthermore, after the bubble dynamics model is established, the bubble growth and the various forces it is subjected to can be calculated based on the bubble dynamics model. Then, based on the bubble growth and the various forces it is subjected to and the flow characteristics of the bubble, the forces acting on the bubble can be calculated and the forces acting on the bubble can be decomposed along the flow direction and the normal direction. If the force balance of the bubble in the flow direction is destroyed ( ) and the force balance in the normal direction is not destroyed ( ), then multiple bubble characteristic parameters are obtained through the force equation of the bubble along the flow direction and the preset bubble detachment parameter calculation equation; if the force balance of the bubble in the flow direction is not destroyed ( ) and the force balance in the normal direction is destroyed ( ), multiple bubble characteristic parameters are obtained through the bubble force equation along the normal direction and the preset bubble detachment parameter calculation equation.

[0096] For example, if the force balance of the bubble in the flow direction is destroyed ( ), while the force balance in the normal direction is not destroyed ( ), the bubble will slip after detaching from the nucleation point. At this time, the force equation of the bubble along the flow direction is shown in Formula 16:

[0097] (Formula 16)

[0098] in, is the force on the bubble along the flow direction, is the bubble buoyancy (N), is the bubble drag force (N), is the virtual mass force (N), is the inertial acceleration body force in the slip direction.

[0099] Furthermore, if the force balance of the bubble in the flow direction is not destroyed ( ), while the force balance in the normal direction is destroyed ( ), the bubble will not slip, but will directly detach from the wall and enter the mainstream. When detaching, the bubble contact diameter is 0. For example, the force equation of the bubble along the normal direction is shown in Formula 17:

[0100] (Formula 17)

[0101] in, is the force on the bubble in the normal direction, is the bubble buoyancy (N), is the bubble shear lift (N), is the bubble growth force (N), is the inertial acceleration body force in the normal direction.

[0102] Furthermore, the bubble detachment parameter calculation equation includes a bubble detachment frequency calculation equation and a bubble nucleation density calculation equation. For example, the bubble detachment frequency calculation equation is shown in Formula 18:

[0103] , ,

[0104] , (Formula 18)

[0105] ,

[0106] in, is the bubble detachment frequency (1 / s), is the bubble waiting time (s), is the bubble growth time (s), which indicates the time it takes for the bubble to separate from the nucleation point. is the thermal diffusivity, is the saturation temperature (K), is the mixture temperature (K), is the wall temperature (K), is the surface tension coefficient (m -1 ), is the vapor phase density (kg / m 3 ), is the latent heat of vaporization (kJ / kg), is the vaporization core radius (m), is the liquid phase thermal conductivity (W / (m K)), is the wall heat flux density (W / (m 2 K))、 is a custom constant, The factor F represents the ability of the vaporization core cavity to be immersed in water and the degree of wetting of the heating surface. , then the entire cavity will be flooded. If , then all cavities are not flooded, i.e. they contain gas or vapor.

[0107] Further, exemplary, the bubble nucleation density calculation equation is shown in Formula 19:

[0108] ,

[0109] (Formula 19)

[0110] in, is the nucleation density (1 / m 2 ), is the bubble rising diameter (m), represents the surface tension coefficient (m -1 ), is the saturation temperature (K), is the effective superheat (K), is the vapor phase density (kg / m3), is the latent heat of vaporization (kJ / kg), is the bubble detachment diameter (m), is the relative density, is the liquid density (kg / m3), As an inhibitory factor, is the wall superheat (K).

[0111] For example, the calculation process of bubble characteristic parameters is as follows Figure 4 As shown, first, the bubble growth model is selected according to the wall fluid parameters and wall temperature, that is, the bubble radius changes with time Then, the various forces acting on the static bubble at different times are calculated based on the bubble diameter and fluid parameters. Among them, the static bubble is not affected by the virtual mass force, and its force diagram is as follows Figure 2 As shown. Then, the forces can be decomposed along the flow direction (x) and the vertical direction (y). If the normal force balance is broken ( ) before the force balance in the flow direction breaks down ( ), the bubble will detach directly without slipping. At this time, the bubble diameter is , bubble rising diameter Diameter of bubble separation Equal, bubble sliding distance is 0, the bubble detachment frequency It is related to the bubble detachment time and waiting time, and the time is the bubble growth time Finally, the bubble detachment frequency can be calculated based on the bubble detachment frequency calculation equation and the bubble nucleation density calculation equation in the bubble detachment parameter calculation equation. and nucleation density .

[0112] Further, refer to Figure 3 , if the stream force balance is broken ( ) before the force balance in the normal direction breaks down ( ), the bubble slips and the flow direction force balance is broken ( ) is the bubble detachment diameter At this time, according to the force equation of the bubble along the normal direction, when the normal force balance is broken ( ) is the diameter of the bubble rising The bubble slip velocity can be obtained by solving the ordinary differential equation based on the force equation of the bubble along the flow direction, and the bubble slip distance can be obtained by integrating the velocity with time. , is the convective heat flux, is the share of other heat flux densities in the wall heat flux density except the convection heat flux density. Finally, the bubble detachment frequency can be calculated according to the bubble detachment frequency calculation equation and the bubble nucleation density calculation equation in the bubble detachment parameter calculation equation. and nucleation density .

[0113] In one embodiment, the wall heat flux model includes a subcooled boiling heat flux model and a boiling critical heat flux model. Based on this, step 103 can be implemented in the following manner: establish a wall heat flux distribution model, and substitute the bubble characteristic parameters into the wall heat flux distribution model to obtain multiple wall heat flux components, wherein the wall heat flux components include a quenching heat flux component, a convection heat flux component, an evaporation heat flux component, a slip heat flux component, and a single-phase vapor heat flux component. In a subcooled boiling environment, the quenching heat flux component, the convection heat flux component, the evaporation heat flux component, and the slip heat flux component can be combined to obtain a subcooled boiling heat flux model; in a boiling critical environment, the total wall dry area can be calculated based on the bubble statistical parameters, and the boiling critical heat flux model can be obtained based on the total wall dry area, the quenching heat flux component, the convection heat flux component, the evaporation heat flux component, the slip heat flux component, and the single-phase vapor heat flux component.

[0114] For example, the quenching heat flux equation in the wall heat flux distribution model is shown in Formula 20:

[0115] (Formula 20)

[0116] in, is the quenching heat flux (kW / m 2 ), is the liquid phase thermal conductivity (W / (m 2 K))、 Liquid density (kg / m 3 ), Liquid phase heat capacity (J / (kg K)), is the bubble waiting time (s), is the wall temperature (K), is the mixture temperature (K), To reduce the factor, is the nucleation density (1 / m 2 ), is the bubble influence factor, is the bubble detachment diameter (m), is the bubble detachment frequency (1 / s).

[0117] For example, the convection heat flux equation in the wall heat flux distribution model is shown in Formula 21:

[0118] (Formula 21)

[0119]

[0120] in, is the convective heat flux (kW / m 2 ), is the convective heat transfer coefficient, which can be obtained through the wall function, is the wall temperature (K), is the mixture temperature (K), is the share of other heat fluxes in the wall heat flux density except the convection heat flux density, To reduce the factor, is the nucleation density (1 / m 2 ), is the bubble influence factor, is the bubble detachment diameter (m), is the bubble waiting time (s), is the bubble detachment frequency (1 / s), is the slip distance (m), is the average diameter of bubbles detached from the wall (m), is the slip time (s).

[0121] For example, the evaporation heat flux equation in the wall heat flux distribution model is shown in Formula 22:

[0122] (Formula 22)

[0123] in, is the evaporation heat flux (kW / m 2 ), R f To reduce the factor, is the bubble volume (m 3 ), Vapor phase density (kg / m 3 ), is the nucleation density (1 / m 2 ), is the latent heat of vaporization (kJ / kg), is the bubble detachment frequency (1 / s), is the bubble rising diameter (m).

[0124] For example, the slip heat flow equation in the wall heat flux distribution model is shown in Formula 23:

[0125] (Formula 23)

[0126] in, is the slip heat flow (kW / m 2 ), is the liquid phase thermal conductivity (W / (m 2 K))、 Liquid density (kg / m 3 ), Liquid phase heat capacity (J / (kg K)), is the bubble waiting time (s), is the wall temperature (K), is the mixture temperature (K), To reduce the factor, is the nucleation density (1 / m 2 ), is the slip distance (m), is the bubble influence factor, is the average diameter of bubbles detached from the wall (m), is the bubble detachment frequency (1 / s), is the slip time (s).

[0127] In the subcooled boiling scenario, the total wall heat flux can be calculated from other heat fluxes. For example, the subcooled boiling heat flux model is shown in Equation 24:

[0128] (Formula 24)

[0129] in, is the total wall heat flux (kW / m 2 ), is the evaporation heat flux (kW / m 2 ), is the quenching heat flux (kW / m 2 ), is the slip heat flow (kW / m 2 ), is the convective heat flux (kW / m 2 ).

[0130] In one embodiment, under a boiling critical environment, a boiling critical heat flux density model can be obtained by the following method: under a boiling critical environment, first, based on a preset data statistical rule, the probability of at least a preset number of bubbles existing in a set influence area is calculated, then based on the probability of at least a preset number of bubbles existing in a set influence area, the average density of hot spots is calculated, and based on the average density of hot spots, the occupancy share of all hot spots is calculated, then the total dry area of ​​the wall is calculated based on the occupancy share of all hot spots, and based on the total dry area of ​​the wall, a first weight coefficient and a second weight coefficient are set, and finally, the product of the first weight coefficient and the single-phase steam heat flux component and the product of the second weight coefficient and the sum of the quenching heat flux component, the convection heat flux component, the evaporation heat flux component and the slip heat flux component are combined to obtain a boiling critical heat flux density model.

[0131] For example, in the boiling critical scenario, the total wall heat flux can be calculated using the bubble statistical parameters. Assuming that the bubbles on the wall follow a Poisson distribution, then The probability that at least one activated bubble exists within the impact area The probability that at least one activated bubble exists within the impact area

[0132] (Formula Twenty-Five)

[0133] wherein, The probability that at least one activated bubble exists within the impact area is the number of activated bubbles, is the impact area (m 2 ), is the average nucleation density (1 / m 2 ), is Euler's number. Based on this, the probability that at least bubbles exist within the impact area can be calculated by Formula Twenty-Six:

[0134] (Formula Twenty-Six)

[0135] wherein, is the probability that at least bubbles exist within the impact area , is the number of activated bubbles, is the number of critical bubbles when dryout occurs, is the impact area, is the average nucleation density (1 / m 2 ), is Euler's number. Assuming that all hotspots are mutually independent, the total hotspot share can be calculated by Formula Twenty-Seven:

[0136] , (Formula Twenty-Seven)

[0137] wherein, is the total hotspot share, is the average hotspot density (1 / m 2 ), is the area of one hotspot (m 2 ), is the average nucleation density (1 / m 2 ), is the probability that at least bubbles exist within the impact area .

[0138] In this scenario, the total dryout area can be calculated by Formula Twenty-Eight:

[0139] ​​ (Formula 28)

[0140] in, is the total dry area, Occupy share for all hot spots, is the Euler number.

[0141] Furthermore, when the heat flux is high, that is, the wall heat flux increases and reaches a certain critical point, the steam convection term will account for the vast majority of the wall heat flux, causing the wall temperature to rise sharply and boiling criticality to occur. Based on this, considering the influence of steam convection, the boiling critical heat flux model is shown in Equation 29:

[0142] , (Formula 29)

[0143] in, is the total wall heat flux (kW / m 2 ), is the evaporation heat flux (kW / m 2 ), is the quenching heat flux (kW / m 2 ), is the slip heat flow (kW / m 2 ), is the convective heat flux (kW / m 2 ), is the single-phase steam heat flow (kW / m2), is the total dry area, is the steam convection heat transfer coefficient, which can be obtained based on the wall function. is the wall temperature (K), is the mixture temperature (K).

[0144] In one embodiment, referring to Figure 5 Step 104 can be implemented in the following manner: first, based on the multiple bubble characteristic parameters calculated by the bubble dynamics model, the total wall heat flux is obtained through the wall heat flux density model, and the total wall heat flux is converted into the volume heat source of the wall; then, the energy exchange term is calculated according to the volume heat source, and the mass exchange term and the momentum exchange term are calculated; then, based on the energy exchange term, the mass exchange term and the momentum exchange term, the non-inertial system mixture model is iteratively solved, and it is determined whether the solution value of the non-inertial system mixture model converges; if the solution value converges, the simulation results of subcooled boiling and boiling critical value are obtained through the solution value of the non-inertial system mixture model.

[0145] Further, as Figures 1 to 5 The specific implementation of the method, the embodiment of the present application provides a subcooled boiling and boiling critical numerical simulation device, such as Figure 6As shown, the device includes:

[0146] The mixture model establishment module 21 can be used to establish a mixture model and add the inertial acceleration volume force generated in the non-inertial system to the momentum conservation equation of the mixture model to obtain the non-inertial system mixture model;

[0147] The bubble dynamics model establishment module 22 is configured to establish a bubble dynamics model based on multiple forces acting on the wall bubbles and the inertial acceleration volume force, and obtain multiple bubble characteristic parameters based on the bubble dynamics model;

[0148] The wall heat flux density model establishment module 23 can be used to establish a wall heat flux distribution model and substitute the bubble characteristic parameters into the wall heat flux distribution model to obtain a wall heat flux density model;

[0149] The boiling heat transfer condition simulation module 24 can be used to couple the wall heat flux model and the non-inertial system mixture model to obtain a two-phase boiling heat transfer and boiling criticality calculation model for simulating the boiling heat transfer condition under motion conditions.

[0150] In a specific application scenario, the mixture model establishment module 21 can be used to establish the inertial acceleration volume force equation generated in the non-inertial system according to the motion characteristics of the fluid in the pipeline, which is used to calculate the inertial acceleration volume force generated by the movement of the object; establish the mass conservation equation, momentum conservation equation, energy conservation equation and cavitation fraction equation of the mixture model, and add the inertial acceleration volume force to the momentum conservation equation as a volume force source term to obtain the non-inertial system mixture model.

[0151] In a specific application scenario, the bubble dynamics model establishment module 22 can be specifically used to establish wall bubble growth and various force equations, wherein the force equations include bubble growth equation, surface tension equation, bubble shear lift equation, bubble growth force equation, bubble drag equation, bubble buoyancy equation, bubble contact pressure equation, bubble dynamic pressure equation and virtual mass force equation; based on the bubble growth and various force equations and the inertial acceleration volume force equation corresponding to the inertial acceleration volume force, a bubble dynamics model is established; based on the bubble dynamics model, a plurality of bubble characteristic parameters are calculated, wherein the bubble characteristic parameters include bubble detachment frequency, bubble detachment diameter, bubble sliding distance, bubble rise diameter and bubble nucleation density.

[0152] In a specific application scenario, the bubble dynamics model establishment module 22 can be specifically used to calculate the bubble growth and the various forces it is subjected to according to the bubble dynamics model; based on the bubble growth and the various forces it is subjected to and the flow characteristics of the bubble, calculate the force acting on the bubble, and decompose the force acting on the bubble along the flow direction and the normal direction; if the force balance of the bubble in the flow direction is destroyed and the force balance in the normal direction is not destroyed, then a plurality of bubble characteristic parameters are obtained through the force equation of the bubble along the flow direction and the preset bubble detachment parameter calculation equation; if the force balance of the bubble in the flow direction is not destroyed and the force balance in the normal direction is destroyed, then a plurality of bubble characteristic parameters are obtained through the force equation of the bubble along the normal direction and the preset bubble detachment parameter calculation equation.

[0153] In a specific application scenario, the wall heat flux density model includes a subcooled boiling heat flux density model and a boiling critical heat flux density model; the wall heat flux density model establishment module 23 can be specifically used to establish a wall heat flux distribution model, and substitute the bubble characteristic parameters into the wall heat flux distribution model to obtain multiple wall heat flux components, wherein the wall heat flux components include quenching heat flux component, convection heat flux component, evaporation heat flux component, slip heat flux component and single-phase steam heat flux component; in a subcooled boiling environment, the quenching heat flux component, the convection heat flux component, the evaporation heat flux component and the slip heat flux component are combined to obtain the subcooled boiling heat flux density model; in a boiling critical environment, the total wall dry area is calculated based on the bubble statistical parameters, and based on the total wall dry area, the quenching heat flux component, the convection heat flux component, the evaporation heat flux component, the slip heat flux component and the single-phase steam heat flux component, the boiling critical heat flux density model is obtained.

[0154] In a specific application scenario, the wall heat flux density model establishment module 23 is used to calculate the probability of at least a preset number of bubbles existing in a set influence area based on preset data statistical rules under a boiling critical environment; calculate the average hot spot density based on the probability of at least a preset number of bubbles existing in the set influence area, and calculate the occupancy share of all hot spots based on the hot spot average density; calculate the total wall dry area based on the occupancy share of all hot spots, and set a first weight coefficient and a second weight coefficient based on the total wall dry area; combine the product of the first weight coefficient and the single-phase steam heat flux component and the product of the second weight coefficient and the sum of the quenching heat flux component, the convection heat flux component, the evaporation heat flux component and the slip heat flux component to obtain the boiling critical heat flux density model.

[0155] In a specific application scenario, the boiling heat exchange condition simulation module 24 can be specifically used to obtain the total wall heat flux through the wall heat flux density model based on the multiple bubble characteristic parameters calculated by the bubble dynamics model, and convert the total wall heat flux into the volume heat source of the wall; calculate the energy exchange term according to the volume heat source, and calculate the mass exchange term and the momentum exchange term; iteratively solve the non-inertial system mixture model based on the energy exchange term, the mass exchange term and the momentum exchange term, and judge whether the solution value of the non-inertial system mixture model converges; if the solution value converges, the simulation results of subcooled boiling and boiling critical value are obtained through the solution value of the non-inertial system mixture model.

[0156] It should be noted that for other corresponding descriptions of the functional units involved in the subcooled boiling and boiling critical numerical simulation device provided in the embodiment of the present application, reference can be made to Figures 1 to 5 The corresponding description in the method will not be repeated here.

[0157] The embodiment of the present application also provides a computer device, which can be specifically a personal computer, a server, a network device, etc. The computer device includes a bus, a processor, a memory and a communication interface, and may also include an input and output interface and a display device. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and computer program in the non-volatile storage medium. The database of the computer device is used to store location information. The network interface of the computer device is used to communicate with an external terminal via a network connection. When the computer program is executed by the processor, the steps in each method embodiment are implemented.

[0158] Those skilled in the art will understand that the structure of the above-mentioned computer device is only a partial structure related to the solution of the present application and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components, or combine certain components, or have a different component arrangement.

[0159] In one embodiment, a computer-readable storage medium is provided. The computer-readable storage medium may be non-volatile or volatile, and stores a computer program thereon. When the computer program is executed by a processor, the steps in the above-mentioned method embodiments are implemented.

[0160] In one embodiment, a computer program product is provided, including a computer program, which implements the steps in the above method embodiments when executed by a processor.

[0161] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the above-mentioned embodiments. In particular, any reference to memory, database, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The databases involved in the various embodiments provided herein may include at least one of a relational database and a non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the various embodiments provided herein may be, but are not limited to, general-purpose processors, graphics processors, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, and the like.

[0162] The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0163] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present application. It should be noted that a person of ordinary skill in the art may make various modifications and improvements without departing from the spirit of the present application, and these modifications and improvements fall within the scope of protection of the present application. Therefore, the scope of protection of the present application shall be determined by the appended claims.

Claims

1. A method for numerical simulation of subcooled boiling and boiling criticality, characterized in that: The method comprises: Establishing a Mixture model, and adding the inertial acceleration body force generated in the non-inertial system to the momentum conservation equation of the Mixture model to obtain the non-inertial system mixture model, specifically comprising: establishing an inertial acceleration body force equation generated in the non-inertial system according to the motion characteristics of the fluid in the pipeline, for calculating the inertial acceleration body force generated by the motion of the object; establishing a mass conservation equation, a momentum conservation equation, an energy conservation equation, and a cavitation fraction equation for the Mixture model, and adding the inertial acceleration body force to the momentum conservation equation as a body force source term to obtain the non-inertial system mixture model; Establishing a bubble dynamics model based on multiple forces acting on the wall bubbles and the inertial acceleration volume force, and obtaining multiple bubble characteristic parameters based on the bubble dynamics model; Establishing a wall heat flux distribution model, and substituting the bubble characteristic parameters into the wall heat flux distribution model to obtain a wall heat flux density model; The wall heat flux model and the non-inertial system mixture model are coupled to obtain a two-phase boiling heat transfer and boiling criticality calculation model for simulating boiling heat transfer conditions under motion conditions, specifically including: based on multiple bubble characteristic parameters calculated by the bubble dynamics model, obtaining a total wall heat flux through the wall heat flux model, and converting the total wall heat flux into a volume heat source of the wall; calculating an energy exchange term based on the volume heat source, and calculating a mass exchange term and a momentum exchange term; based on the energy exchange term, the mass exchange term and the momentum exchange term, iteratively solving the non-inertial system mixture model, and judging whether the solution value of the non-inertial system mixture model converges; if the solution value converges, obtaining simulation results of subcooled boiling and boiling critical values ​​through the solution value of the non-inertial system mixture model.

2. The method for numerical simulation of subcooled boiling and boiling criticality according to claim 1, characterized in that: The bubble dynamics model is established based on the multiple forces acting on the wall bubbles and the inertial acceleration volume force, and multiple bubble characteristic parameters are obtained based on the bubble dynamics model, including: Establishing wall bubble growth and multiple force equations, including bubble growth equation, surface tension equation, bubble shear lift equation, bubble growth force equation, bubble drag equation, bubble buoyancy equation, bubble contact pressure equation, bubble dynamic pressure equation, and virtual mass force equation; Establishing a bubble dynamics model based on the bubble growth and multiple force equations and the inertial acceleration body force equation corresponding to the inertial acceleration body force; According to the bubble dynamics model, a plurality of bubble characteristic parameters are calculated, wherein the bubble characteristic parameters include bubble detachment frequency, bubble detachment diameter, bubble sliding distance, bubble rising diameter and bubble nucleation density.

3. The method for numerical simulation of subcooled boiling and boiling criticality according to claim 2, characterized in that: According to the bubble dynamics model, a plurality of bubble characteristic parameters are calculated, including: Calculating bubble growth and various forces acting on the bubble according to the bubble dynamics model; Based on the bubble growth and the multiple forces acting on it and the flow characteristics of the bubble, the force acting on the bubble is calculated, and the force acting on the bubble is decomposed along the flow direction and the normal direction; If the force balance of the bubble in the flow direction is destroyed but the force balance in the normal direction is not destroyed, multiple bubble characteristic parameters are obtained through the force equation of the bubble along the flow direction and the preset bubble detachment parameter calculation equation; If the force balance of the bubble in the flow direction is not destroyed and the force balance in the normal direction is destroyed, multiple bubble characteristic parameters are obtained through the force equation of the bubble along the normal direction and the preset bubble detachment parameter calculation equation.

4. The method for numerical simulation of subcooled boiling and boiling criticality according to claim 1, characterized in that: The wall heat flux density model includes a subcooled boiling heat flux density model and a boiling critical heat flux density model; then the wall heat flux distribution model is established, and the bubble characteristic parameters are substituted into the wall heat flux distribution model to obtain the wall heat flux density model, including: Establishing a wall heat flux distribution model, and substituting the bubble characteristic parameters into the wall heat flux distribution model to obtain multiple wall heat flux components, wherein the wall heat flux components include a quenching heat flux component, a convection heat flux component, an evaporation heat flux component, a slip heat flux component, and a single-phase vapor heat flux component; Under a subcooled boiling environment, the extinction heat flux component, the convection heat flux component, the evaporation heat flux component, and the slip heat flux component are combined to obtain the subcooled boiling heat flux density model; Under a boiling critical environment, the total wall dry-up area is calculated based on bubble statistical parameters, and the boiling critical heat flux density model is obtained based on the total wall dry-up area, the quenching heat flux component, the convection heat flux component, the evaporation heat flux component, the slip heat flux component and the single-phase vapor heat flux component.

5. The method for numerical simulation of subcooled boiling and boiling criticality according to claim 4, characterized in that: The method comprises: calculating the total wall dry-up area based on bubble statistical parameters under the boiling critical environment, and obtaining the boiling critical heat flux density model based on the total wall dry-up area, the quenching heat flux component, the convection heat flux component, the evaporation heat flux component, the slip heat flux component, and the single-phase vapor heat flux component. Under a boiling critical environment, based on preset data statistical rules, the probability of at least a preset number of bubbles existing within a set impact area is calculated; Calculating an average density of hot spots based on a probability that at least a preset number of bubbles exist within the set impact area, and calculating occupancy shares of all hot spots based on the average density of hot spots; Calculating a total dry area of ​​the wall based on all hot spot occupancy shares, and setting a first weight coefficient and a second weight coefficient based on the total dry area of ​​the wall; The boiling critical heat flux density model is obtained by combining the product of the first weight coefficient and the single-phase vapor heat flux component and the product of the second weight coefficient and the sum of the quenching heat flux component, the convection heat flux component, the evaporation heat flux component and the slip heat flux component.

6. A numerical simulation device for supercooled boiling and boiling criticality, characterized in that: The device comprises: A mixture model establishment module is used to establish a mixture model and add the inertial acceleration body force generated in the non-inertial system to the momentum conservation equation of the mixture model to obtain the non-inertial system mixture model. Specifically, the module includes: establishing the inertial acceleration body force equation generated in the non-inertial system according to the motion characteristics of the fluid in the pipeline, which is used to calculate the inertial acceleration body force generated by the motion of the object; establishing the mass conservation equation, momentum conservation equation, energy conservation equation and cavitation fraction equation of the mixture model, and adding the inertial acceleration body force to the momentum conservation equation as a body force source term to obtain the non-inertial system mixture model; a bubble dynamics model establishment module, configured to establish a bubble dynamics model based on multiple forces acting on the wall bubbles and the inertial acceleration volume force, and obtain multiple bubble characteristic parameters based on the bubble dynamics model; a wall heat flux density model establishment module, for establishing a wall heat flux distribution model, and substituting the bubble characteristic parameters into the wall heat flux distribution model to obtain a wall heat flux density model; The boiling heat transfer condition simulation module is used to couple the wall heat flux model and the non-inertial system mixture model to obtain a two-phase boiling heat transfer and boiling criticality calculation model, which is used to simulate the boiling heat transfer condition under motion conditions. Specifically, the module includes: based on multiple bubble characteristic parameters calculated by the bubble dynamics model, the total wall heat flux is obtained through the wall heat flux model, and the total wall heat flux is converted into a volume heat source of the wall; based on the volume heat source, the energy exchange term is calculated, and the mass exchange term and the momentum exchange term are calculated; based on the energy exchange term, the mass exchange term and the momentum exchange term, the non-inertial system mixture model is iteratively solved, and it is determined whether the solution value of the non-inertial system mixture model converges; if the solution value converges, the simulation results of subcooled boiling and boiling critical values ​​are obtained through the solution value of the non-inertial system mixture model.

7. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 5 is implemented.

8. A computer device comprising a storage medium, a processor, and a computer program stored on the storage medium and executable on the processor, wherein: When the processor executes the computer program, the method according to any one of claims 1 to 5 is implemented.

Images