Method and system for predicting bubble dynamic characteristics considering phase change effect near water area boundary

Abstract

The application relates to a method and system for predicting the dynamic characteristics of a bubble near a water boundary considering phase transition effects, and relates to the field of cavitation bubble dynamics. The method solves the problems that existing analytical models are difficult to accurately reflect the deformation and migration track of a bubble when dealing with the bubble near a boundary, and the prediction results deviate from the actual situation, because the wall effect and the non-uniformity of a flow field are ignored. The method comprises the following steps: inputting initial conditions and environmental parameters of the bubble; calculating a boundary-induced flow field pressure; calculating a flow field environmental pressure based on the boundary-induced pressure and the hydrostatic pressure; calculating the internal pressure of the bubble based on the thermodynamic phase transition effect; calculating the enthalpy difference of the bubble surface based on the flow field environmental pressure and the internal pressure of the bubble; substituting the obtained result into a unified form of a bubble equation considering phase transition and migration effects and solving the bubble equation to obtain the equivalent radius and migration distance of the bubble changing with time. The method is also suitable for the field of predicting the equivalent radius calculation of a single cavitation bubble during boundary movement.

Description

Technical Field

[0001] This invention relates to the field of cavitation bubble dynamics technology, specifically to a method and system for predicting bubble dynamic characteristics near the boundary of a water body, taking into account phase change effects. Background Technology

[0002] Cavitation is the process by which liquids generate and develop bubbles under certain special conditions (e.g., local pressure drops below saturated vapor pressure), and it is widely observed in hydraulic machinery, ship propellers, and microfluidic devices. The nucleation, growth, collapse, and interaction of bubbles with surrounding boundaries in fluids not only lead to decreased equipment efficiency but can also corrode material surfaces or cause noise and vibration. Simultaneously, researchers are continuously exploring and utilizing the high temperature, high pressure, and strong turbulence characteristics of cavitation for applications such as ultrasonic cleaning, wastewater treatment, material modification, and microfluidic chips. For example, acoustically induced cavitation can generate intense local high temperatures and pressures at the moment of microbubble collapse, promoting increased chemical reaction rates; while in water treatment or emulsification processes, the turbulence generated by cavitation can enhance mixing and promote dissolution, thereby further improving treatment efficiency. Therefore, how to avoid or mitigate the harms of cavitation while effectively utilizing its advantages has become an important research topic in fluid mechanics and engineering applications.

[0003] Currently, research methods for predicting bubble motion characteristics under cavitation conditions mainly fall into two categories: theoretical analytical models and numerical simulations. Among these, analytical methods based on the Rayleigh–Plesset equation and its improved models can provide a relatively simplified description of the growth and collapse process of a single bubble, offering high computational efficiency. However, in the dynamic process of bubbles near the boundary, due to the non-uniform distribution of the flow field and wall effects, traditional analytical models often struggle to accurately characterize bubble deformation and trajectory. Meanwhile, while commercial CFD software or self-developed numerical simulation methods can capture the dynamic evolution of the gas-liquid interface well, they place higher demands on mesh generation accuracy, boundary condition settings, and the rationality of the phase transition model when dealing with phase change effects (phase change between the liquid and the vapor inside the bubble) and the coupling of wall boundary conditions. These requirements necessitate greater computational resources.

[0004] Current research on simulating bubble motion near boundaries suffers from the following shortcomings: First, traditional analytical models, such as the Rayleigh–Plesset equation, fail to accurately reflect bubble deformation and migration trajectories when dealing with bubbles near boundaries due to neglecting wall effects and flow field inhomogeneities, leading to predictions that deviate from reality. Second, while existing CFD methods have advantages in simulating complex flow fields and gas-liquid interface dynamics, they still face challenges in accurately coupling phase transition effects and require significant computational resources. Furthermore, although existing high-order theoretical models consider some complex factors, they are still insufficient in providing a comprehensive coupled description of phase transition processes and bubble dynamics, and the complexity of their analytical solutions limits their practicality in variable engineering environments. Summary of the Invention

[0005] This invention addresses several issues in existing technologies. Firstly, analytical models like the Rayleigh–Plesset equation, when dealing with bubbles near the boundary, fail to accurately reflect bubble deformation and migration trajectories due to neglecting wall effects and flow field inhomogeneities, leading to predictions that deviate from reality. Secondly, while existing CFD methods excel in simulating complex flow fields and gas-liquid interface dynamics, they still face challenges in accurately coupling phase transition effects and require significant computational resources. Furthermore, while existing high-order theoretical models consider some complex factors, they are insufficient in providing a comprehensive coupled description of phase transition processes and bubble dynamics, and the complexity of their analytical solutions limits their practicality in variable engineering environments.

[0006] To solve the above-mentioned technical problems, the present invention is achieved through the following technical solution:

[0007] Option 1: This invention proposes a method for predicting bubble dynamics characteristics near the boundary of a water body, considering phase transition effects. The method includes the following steps:

[0008] S1. Input the initial conditions and environmental parameters of the bubble;

[0009] S2. Use the mirror method to set the parameters of the mirror bubble and calculate the boundary-induced flow field pressure. Calculate the flow field environment pressure based on the boundary-induced pressure and hydrostatic pressure.

[0010] S3. Calculate the internal pressure of the bubble based on the thermodynamic phase change effect;

[0011] S4. Calculate the enthalpy difference on the bubble surface based on the flow field environment pressure described in S2 and the internal pressure of the bubble described in S3;

[0012] S5. Substitute the results obtained from S2 to S4 into the unified form of the bubble equation that considers phase change and migration effects and solve it to obtain the bubble equivalent radius and migration distance that change with time, thus completing the prediction of bubble dynamics characteristics considering phase change effects near the water boundary.

[0013] Furthermore, a preferred embodiment is provided, wherein the initial conditions of the bubble input in S1 include the initial radius of the bubble, the initial internal pressure of the bubble, the water depth at the center of the bubble, and the distance between the bubble and the boundary; the environmental parameters include fluid density, flow field reference sound velocity, fluid saturated vapor pressure, fluid surface tension, fluid dynamic viscosity, flow field temperature, and gravitational acceleration.

[0014] Furthermore, a preferred embodiment is provided, wherein the method for setting the parameters of the mirror bubble using the mirror method in S2 is as follows:

[0015]

[0016] Where R and v are the radius and migration velocity of the actual bubble, respectively. I v I For mirrored bubble v I The radius and migration speed.

[0017] Furthermore, a preferred embodiment is provided, wherein the method for calculating the flow field environment pressure based on boundary induced pressure in S2 is as follows:

[0018]

[0019] Among them, o R Let ρ be the position vector of the actual bubble, and ρ be the fluid density. For the velocity potential induced by the mirror bubble, u I Let u be the flow velocity induced by the mirror bubble; where u is the velocity of the flow field induced by the mirror bubble. I and The expression is:

[0020]

[0021] Where H is the enthalpy difference and d is the distance between the bubble center and the boundary at time t.

[0022] Furthermore, a preferred embodiment is provided, wherein the method for calculating the internal pressure of the bubble based on the thermodynamic phase transition effect in step 3 is as follows:

[0023]

[0024] In the formula, P g The pressure of the gas inside the bubble (i.e., vapor and non-condensable gas), σ is the surface tension, μ is the viscosity of the fluid, and ρ is the viscosity of the fluid. g This represents the average density of the gas inside the bubble. This represents the evaporation or condensation rate per unit area of ​​the bubble surface.

[0025] Furthermore, a preferred embodiment is provided, wherein the evaporation or condensation rate per unit area of ​​the bubble surface is... The calculation method is as follows:

[0026]

[0027] Where, α m For the adaptation factor; P v (T l ) is the liquid at T l Saturated vapor pressure P at temperature vb β represents the actual vapor pressure inside the bubble; β is the correction coefficient.

[0028] Furthermore, a preferred embodiment is provided, wherein the unified bubble equation considering phase transition and migration effects in step four is:

[0029]

[0030] Where υ is the migration velocity vector and γ is the drag force coefficient.

[0031] Option 2: A bubble dynamics prediction system considering phase change effects near the water boundary, the system comprising:

[0032] The input module is used to input the initial conditions and environmental parameters of the bubble;

[0033] The flow field environment pressure calculation module is used to set the parameters of the mirror bubble using the mirror method and calculate the boundary-induced flow field pressure. It calculates the flow field environment pressure based on the boundary-induced pressure and the hydrostatic pressure.

[0034] The bubble internal pressure calculation module is used to calculate the internal pressure of a bubble based on the thermodynamic phase change effect.

[0035] The bubble surface enthalpy difference calculation module is used to calculate the bubble surface enthalpy difference based on the flow field environment pressure calculated by the flow field environment pressure calculation module and the bubble internal pressure calculated by the bubble internal pressure calculation module.

[0036] The prediction module is used to substitute the results obtained from the flow field environment pressure calculation module and the bubble surface enthalpy difference calculation module into the unified form of the bubble equation that considers phase change and migration effects, and solve it to obtain the bubble equivalent radius and migration distance that change with time, thus completing the prediction of bubble dynamic characteristics near the water boundary considering phase change effects.

[0037] Option 3: A computer device, including a memory and a processor, wherein the memory stores a computer program, and when the processor runs the computer program stored in the memory, the processor executes the method described in any one of Options 1.

[0038] Option 4: A computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the method described in any one of Options 1.

[0039] The advantages of this invention are:

[0040] The present invention presents a method and system for predicting bubble dynamics characteristics near water boundaries considering phase change effects. This method, which incorporates phase change effects, offers a valuable research and application approach for calculating bubble motion near boundaries. The method not only accurately predicts the dynamic changes in bubble radius but also simultaneously simulates the spatiotemporal migration behavior of bubbles, comprehensively considering the impact of phase change processes on bubble dynamics. This significantly improves the accuracy and reliability of cavitation bubble motion prediction while saving substantial computational resources.

[0041] This invention is based on a bubble equation that takes into account phase change and migration effects and has a unified form. It uses the mirror method to design the flow field pressure induced by the boundary on the bubble, while also considering the thermodynamic effects of migration and compressibility. This achieves the goal of accurately predicting the bubble motion behavior near the boundary, providing a theoretical basis and fundamental technical support for engineering construction.

[0042] This invention is also applicable to the field of theoretical calculation methods for accurately predicting the equivalent radius and migration of a single cavitation bubble during its movement near the boundary. Attached Figure Description

[0043] Figure 1 This is a flowchart of the bubble dynamics prediction method considering phase change effects near the water boundary as described in Implementation Method 1.

[0044] Figure 2 This is a schematic diagram of the mirror method described in Implementation Method 1.

[0045] Figure 3 This is a schematic diagram of the phase transition effect on the surface of the air bubbles as described in Embodiment 1.

[0046] Figure 4 This is a time-history curve of the predicted equivalent radius of the electric spark bubble near the free surface as described in Implementation Method Eleven.

[0047] Figure 5 This is a time-history curve of the predicted migration distance of spark bubbles near the free surface as described in Implementation Method 1.

[0048] Figure 6This is a time-history curve of the predicted equivalent radius of a laser pulse bubble near a rigid boundary, as described in Implementation Method 1.

[0049] Figure 7 This is a time-history curve of the predicted laser pulse bubble migration distance near the rigid boundary as described in Implementation Method 1. Detailed Implementation

[0050] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this application, and not all of them.

[0051] Implementation Method 1: This implementation method proposes a method for predicting bubble dynamics characteristics near the boundary of a water body, considering the phase transition effect. The method includes the following steps:

[0052] S1. Input the initial conditions and environmental parameters of the bubble;

[0053] S2. Use the mirror method to set the parameters of the mirror bubble and calculate the boundary-induced flow field pressure. Calculate the flow field environment pressure based on the boundary-induced pressure and hydrostatic pressure.

[0054] S3. Calculate the internal pressure of the bubble based on the thermodynamic phase change effect;

[0055] S4. Calculate the enthalpy difference on the bubble surface based on the flow field environment pressure described in S2 and the internal pressure of the bubble described in S3;

[0056] S5. Substitute the results obtained from S2 to S4 into the unified form of the bubble equation that considers phase change and migration effects and solve it to obtain the bubble equivalent radius and migration distance that change with time, thus completing the prediction of bubble dynamics characteristics considering phase change effects near the water boundary.

[0057] Implementation Method 2: This implementation method further defines the bubble dynamics prediction method considering phase change effects near the water boundary described in Implementation Method 1. The initial conditions of the bubble input in S1 include the initial radius of the bubble, the initial internal pressure of the bubble, the water depth at the center of the bubble, and the distance between the bubble and the boundary. The environmental parameters include fluid density, flow field reference sound velocity, fluid saturated vapor pressure, fluid surface tension, fluid dynamic viscosity, flow field temperature, and gravitational acceleration.

[0058] Implementation Method 3: This implementation method further defines the bubble dynamics prediction method considering phase change effects near the water boundary described in Implementation Method 1. The method for setting the parameters of the mirror bubble using the mirror method in S2 is as follows:

[0059]

[0060] Where R and v are the radius and migration velocity of the actual bubble, respectively. I v I For mirrored bubble v I The radius and migration speed.

[0061] Implementation Method Four: This implementation method further defines the bubble dynamics prediction method considering phase change effects near the water boundary described in Implementation Method One. The method for calculating the flow field environmental pressure based on boundary-induced pressure in S2 is as follows:

[0062]

[0063] Among them, o R Let ρ be the position vector of the actual bubble, and ρ be the fluid density. For the velocity potential induced by the mirror bubble, u I Let u be the flow velocity induced by the mirror bubble; where u is the velocity of the flow field induced by the mirror bubble. I and The expression is:

[0064]

[0065] Where H is the enthalpy difference and d is the distance between the bubble center and the boundary at time t.

[0066] Implementation Method 5: This implementation method further defines the bubble dynamics prediction method considering phase change effects near the water boundary described in Implementation Method 1. The method for calculating the internal pressure of the bubble based on the thermodynamic phase change effect in step 3 is as follows:

[0067]

[0068] In the formula, P g The pressure of the gas inside the bubble (i.e., vapor and non-condensable gas), σ is the surface tension, μ is the viscosity of the fluid, and ρ is the viscosity of the fluid. g This represents the average density of the gas inside the bubble. This represents the evaporation or condensation rate per unit area of ​​the bubble surface.

[0069] Implementation Method Six: This implementation method further defines the bubble dynamics prediction method considering phase change effects near the water boundary described in Implementation Method Five, specifically the evaporation or condensation rate per unit area of ​​the bubble surface. The calculation method is as follows:

[0070]

[0071] Where, α m For the adaptation factor; P v (T l) is the liquid at T l Saturated vapor pressure P at temperature vb β represents the actual vapor pressure inside the bubble; β is the correction coefficient.

[0072] Implementation Method Seven: This implementation method further defines the bubble dynamics prediction method considering phase change effects near the water boundary described in Implementation Method One. The unified bubble equation considering phase change and migration effects in step four is as follows:

[0073]

[0074] Where υ is the migration velocity vector and γ is the drag force coefficient.

[0075] Implementation Method 8: This implementation method proposes a bubble dynamics prediction system considering phase change effects near the boundary of a water body. The system includes:

[0076] The input module is used to input the initial conditions and environmental parameters of the bubble;

[0077] The flow field environment pressure calculation module is used to set the parameters of the mirror bubble using the mirror method and calculate the boundary-induced flow field pressure. It calculates the flow field environment pressure based on the boundary-induced pressure and the hydrostatic pressure.

[0078] The bubble internal pressure calculation module is used to calculate the internal pressure of a bubble based on the thermodynamic phase change effect.

[0079] The bubble surface enthalpy difference calculation module is used to calculate the bubble surface enthalpy difference based on the flow field environment pressure calculated by the flow field environment pressure calculation module and the bubble internal pressure calculated by the bubble internal pressure calculation module.

[0080] The prediction module is used to substitute the results obtained from the flow field environment pressure calculation module and the bubble surface enthalpy difference calculation module into the unified form of the bubble equation that considers phase change and migration effects, and solve it to obtain the bubble equivalent radius and migration distance that change with time, thus completing the prediction of bubble dynamic characteristics near the water boundary considering phase change effects.

[0081] Implementation Method Nine: This implementation method provides a computer device, including a memory and a processor. The memory stores a computer program. When the processor runs the computer program stored in the memory, the processor executes the method described in any one of Implementation Methods One to Seven.

[0082] Implementation Method 10: This implementation method provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the method described in any one of Implementation Methods 1 to 7.

[0083] Implementation Method Eleven: This implementation method provides an example, which is used to explain the above-described implementation methods one through eight. See also... Figures 1 to 7 This embodiment is described below, and the specific examples are as follows:

[0084] Step 1: Input the initial conditions and environmental parameters for the bubble.

[0085] Step 2: Based on the boundary type and the initial distance between the bubble and the boundary, design an expression for the environmental pressure at the bubble center under the influence of the boundary.

[0086] Step 3: Based on the phase transition effect, calculate the internal pressure of the bubble, and combine this with the calculation of the enthalpy difference on the bubble surface in Step 2.

[0087] Step four: Substitute the results obtained in steps two and three into the bubble equation with a unified form that considers phase transition and migration effects, and solve it to obtain the bubble equivalent radius and migration distance that change with time.

[0088] Specifically, in step one, the initial conditions and environmental parameters of the bubble are input. The initial conditions of the bubble include the initial radius of the bubble, the initial internal pressure of the bubble, the water depth at the center of the bubble, and the distance between the bubble and the boundary; the environmental parameters include the fluid density, the reference sound velocity of the flow field, the saturated vapor pressure of the fluid, the surface tension of the fluid, the dynamic viscosity of the fluid, the temperature of the flow field, and the acceleration due to gravity. By modifying the above initial information parameters, the initial state of the bubble and the water environment conditions can be set.

[0089] Furthermore, in step two, based on the boundary type, the boundary reflection coefficient α is set, a coordinate system is established and the x-axis coincides with the boundary surface, and the parameters of the mirror bubble are set using the mirror method:

[0090]

[0091] Where R and v are the radius and migration velocity of the actual bubble, respectively. I v I For mirrored bubble v I The radius and migration speed. The principle of the mirror method is as follows: Figure 2 As shown.

[0092] Furthermore, based on the Bernoulli equation, the flow field pressure P induced at the bubble center by the boundary at time t is obtained. B expression:

[0093]

[0094] Among them, o R Let ρ be the position vector of the actual bubble, and ρ be the fluid density. For the velocity potential induced by the mirror bubble, u ILet u be the flow velocity induced by the mirror bubble. I and The expression is:

[0095]

[0096] Where H is the enthalpy difference and d is the distance between the bubble center and the boundary at time t. According to equation (2), the ambient pressure P at the bubble center can be obtained. a =P B +P ∞ P ∞ =ρgh is the hydrostatic pressure.

[0097] Furthermore, in step three, a thermodynamic model is introduced to consider the phase change effect. Based on the pressure balance relationship on the bubble surface, the formula for calculating the bubble surface pressure taking into account the phase change is obtained:

[0098]

[0099] In the formula, P g Let σ be the pressure of the gas (vapor and non-condensable gas) inside the bubble, σ be the surface tension, μ be the viscosity of the fluid, and ρ be the viscosity of the fluid. g This represents the average density of the gas inside the bubble. Let be the evaporation or condensation rate per unit area of ​​the bubble surface. Both the vapor and non-condensable gas in the bubble are considered ideal gases, and both are treated as being uniformly distributed within the bubble. The principle of the phase change effect on the bubble surface is as follows: Figure 3 As shown. P g Calculated based on the ideal gas law:

[0100]

[0101] Where M is the total mass of steam and non-condensable gas, and T b R is the internal temperature of the bubble, V is the volume of the bubble, and R is the internal temperature of the bubble. g T is the constant of the gas mixture inside the bubble. b and R g Calculate using the following formula:

[0102]

[0103] Among them, c v and c p Specific heat capacity of gas under constant volume and constant pressure, S b This represents the surface area of ​​the bubble. κ is the heat exchange rate per unit time on the bubble surface, and T is the heat transfer coefficient. l M is the temperature of the surrounding fluid; v and M n R represents the mass of steam and the mass of non-condensable gas, respectively.v and R n These are the gas constants representing the masses of steam and non-condensable gas, respectively.

[0104] Furthermore, the phase transition rate in equation (4) Calculate using the following formula:

[0105]

[0106] Where, α m For the adaptation factor; P v (T l ) is the liquid at T l Saturated vapor pressure P at temperature vb Γ represents the actual steam pressure inside the bubble; Γ is the correction coefficient.

[0107] Furthermore, combining equation (4) and the environmental pressure P obtained in step two... a Calculate the enthalpy difference H using the following formula:

[0108]

[0109] Steps two and three together constitute an iterative solution process. In step two, P is initially solved. a In equation (3), the enthalpy difference is calculated as H = (P b -P ∞ ) / ρ is calculated, and then P is solved. B Then, proceed to step three to obtain the enthalpy difference after the first iteration. Subsequently, the initial enthalpy difference within each time step is calculated based on the enthalpy difference from the previous iteration.

[0110] Furthermore, in step four, the bubble equation with a unified form used to describe bubble pulsation and migration is:

[0111]

[0112] Where υ is the migration velocity vector and γ is the drag force coefficient. Solving equation (9) yields the equivalent radius and migration distance of the bubble as a function of time during its motion. Figure 4 , Figure 5 The figures show the time-history curves of the predicted equivalent radius and migration distance of spark bubbles near the free surface, respectively, according to the present invention. Figure 6 , Figure 7 The figures show the time-history curves of the equivalent radius and migration distance of the laser pulse bubble near the rigid boundary predicted by the present invention.

[0113] Those skilled in the art will understand that the above description is merely a preferred embodiment of the present invention, and the features described in the various embodiments and / or claims of this disclosure can be combined or combined in various ways, even if such combinations or combinations are not explicitly described in this disclosure. This is not intended to limit the present invention. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions described in the foregoing embodiments or make equivalent substitutions for some of the technical features. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

[0114] Although preferred embodiments of the invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including both the preferred embodiments and all changes and modifications falling within the scope of the invention. Clearly, those skilled in the art can make various alterations and modifications to the invention without departing from its spirit and scope. Thus, if these modifications and modifications of the invention fall within the scope of the claims and their equivalents, the invention is also intended to include these modifications and modifications.

Claims

1. A method for predicting the dynamic behavior of a bubble near a water body boundary taking into account phase change effects, characterized in that, The method comprises the following steps: S1, input initial conditions and environmental parameters of the bubble; S2, set mirror bubble parameters by using mirror method and calculate boundary-induced flow field pressure, calculate flow field environmental pressure based on boundary-induced pressure and hydrostatic pressure; S3, calculate bubble internal pressure based on thermodynamic phase change effect; S4, calculate bubble surface enthalpy difference based on the flow field environmental pressure in S2 and the bubble internal pressure in S3; S5, substitute the results obtained from S2 to S4 into the unified form of bubble equation considering phase change and migration effect and solve, to obtain bubble equivalent radius and migration distance changing with time, and complete prediction of bubble dynamic characteristics near the water area boundary considering phase change effect; The method for calculating bubble internal pressure based on thermodynamic phase change effect in S3 is: wherein Pgas is the pressure of the gas inside the bubble, i.e. of the steam and the non-condensable gas, Pgas is the pressure of the gas inside the bubble, i.e. of the steam and the non-condensable gas, Pgas is the pressure of the gas inside the bubble, i.e. of the steam and the non-condensable gas, Pgas is the pressure of the gas inside the bubble, i.e. of the steam and the non-condensable gas; Pgas is the pressure of the gas inside the bubble, i.e. of the steam and the non-condensable gas; Evaporation or condensation rate per unit area of the bubble surface The calculation method is: wherein, is the adaptation factor; is the saturation vapor pressure of the liquid at temperature, is the actual vapor pressure inside the bubble; is the correction factor; The unified form of bubble equation considering phase change and migration effect in S4 is: wherein is the migration velocity vector, is the drag coefficient.

2. The method of claim 1, wherein, The initial conditions of the bubble input in S1 include initial radius of the bubble, initial internal pressure of the bubble, water depth at the center of the bubble, and distance between the bubble and the boundary; the environmental parameters include fluid density, flow field reference sound speed, fluid saturated vapor pressure, fluid surface tension, fluid dynamic viscosity, flow field temperature, and gravity acceleration.

3. The method of claim 1, wherein, The method for setting mirror bubble parameters by using mirror method in S2 is: wherein, , is the radius and migration velocity magnitude of the real bubble, , is the radius and migration velocity magnitude of the mirror bubble .

4. The method of claim 1, wherein, The method for calculating flow field environmental pressure based on boundary-induced pressure in S2 is: wherein, is the position vector of the real bubble, is the fluid density, is the mirror bubble-induced velocity potential, is the mirror bubble-induced flow field velocity; wherein, and the expression for is: wherein, is the enthalpy difference, is the distance of the bubble center from the boundary at the moment.

5. A system for predicting the dynamic behavior of a bubble near a water boundary taking into account phase change effects, characterized in that, The system is realized based on the method in claim 1, and the system comprises: An input module for inputting initial conditions and environmental parameters of the bubble; A flow field environmental pressure calculation module for setting mirror bubble parameters by using mirror method and calculating boundary-induced flow field pressure, and calculating flow field environmental pressure based on boundary-induced pressure and hydrostatic pressure; A bubble internal pressure calculation module for calculating bubble internal pressure based on thermodynamic phase change effect; A bubble surface enthalpy difference calculation module for calculating bubble surface enthalpy difference based on the flow field environmental pressure calculated by the flow field environmental pressure calculation module and the bubble internal pressure calculated by the bubble internal pressure calculation module; A prediction module for substituting the results obtained from the flow field environmental pressure calculation module to the bubble surface enthalpy difference calculation module into the unified form of bubble equation considering phase change and migration effect and solving, to obtain bubble equivalent radius and migration distance changing with time, and complete prediction of bubble dynamic characteristics near the water area boundary considering phase change effect.

6. Computer device comprising a memory and a processor, characterized in that The memory stores a computer program, and when the processor executes the computer program stored in the memory, the processor executes the method in any one of claims 1 to 4.

7. A computer readable storage medium characterized by, The computer readable storage medium stores a computer program, and when the processor executes the computer program, the steps of the method in any one of claims 1 to 4 are realized.

Images