Multi-scale solving method for cavitation flow field numerical simulation
By combining the bidirectional coupling method of Euler frame and Lagrangian frame, the steam structures at different scales in the cavitation flow field are solved, and the problem of inconvenient steam structure conversion in the prior art is achieved, and a higher accuracy and flexibility of numerical simulation of the cavitation flow field is achieved.
Patent Information
- Application Number
- CN202411526117.3
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-30
- Publication Date
- 2025-02-25
- Estimated Expiration
- Not applicable · inactive patent
AI Technical Summary
The existing two-scale method has limitations when dealing with the transformation of steam structures at different scales, making it difficult to ensure the reasonable conversion of bubbles between different frames, resulting in the impact of the continuity and accuracy of the simulation.
Using a two-way coupling method combining Euler frame and Lagrangian frame, the vacuole scale is used to determine the vacuole scale, and the large-scale vacuole is processed in the Euler frame and small-scale vacuoles are processed in the Lagrangian frame respectively, and the reasonable transformation of the steam structure between different scales is achieved through the scale conversion criteria.
It improves the accuracy and flexibility of numerical simulation of cavitation flow field, can accurately track the evolution of small-scale bubbles, and enhances the ability to capture cavitation noise and cavitation decapitation phenomena.
Smart Images
Figure CN119514407A_ABST
Abstract
Description
Technical Field
[0001] The invention relates to the technical field of numerical simulation of cavitation flow field, and in particular to multi-scale solution of numerical simulation of cavitation flow field. Background Art
[0002] Cavitation flow phenomenon is widely present in high-speed hydraulic machinery and equipment, such as marine propellers, turbines, water pumps, fuel nozzles, etc. When the local pressure in the fluid flow drops below the saturated vapor pressure, the liquid will vaporize, forming cavitation bubbles. This cavitation phenomenon will cause noise, vibration, efficiency loss and cavitation erosion of the equipment, so the study of cavitation phenomenon has important engineering significance. In recent years, with the development of computer simulation technology, numerical simulation methods have gradually become an important means to study cavitation phenomena, helping engineers to deeply understand the formation and evolution of cavitation and its impact on the flow field.
[0003] Existing numerical simulation methods for cavitation flows are mainly based on the Euler-Euler model, which discretizes the flow field through a grid, captures the vapor-liquid interface based on the volume fraction (VOF) or Level-Set method, and uses the cavitation model to describe the phase change mass transfer between the vapor and the liquid. The Euler-Euler model has high efficiency and accuracy when dealing with large-scale cavitation bubbles and is suitable for simulating macroscopic cavitation phenomena. However, the calculation accuracy of the Euler-Euler model is heavily dependent on the grid resolution, especially when simulating small-scale bubbles and complex flow structures. Due to the limitation of the grid scale, it is difficult to accurately capture these small-scale phenomena. The existence of small-scale bubbles is the main cause of cavitation noise and cavitation erosion. Although they have little effect on macroscopic flow field characteristics such as lift and drag coefficients, they have a direct impact on the noise characteristics and life of the equipment. Therefore, how to accurately simulate these small-scale cavitation phenomena is an important challenge in the existing technology.
[0004] Compared with the Euler framework, the Lagrangian method is an explicit method that is not limited by the grid resolution and is suitable for tracking small-scale bubble groups in the flow field. The Lagrangian model regards each bubble as a discrete particle and describes its motion and deformation through Newton's second law and the bubble dynamics equation. However, when using the Lagrangian method alone to deal with a large number of bubbles or large-scale cavitation, its computational complexity and consumption of computing resources increase significantly. Therefore, the existing technology usually faces the dilemma of balancing computational efficiency and simulation accuracy in application.
[0005] In order to overcome these limitations, researchers have begun to try to combine the Euler method and the Lagrangian method in recent years, using their respective advantages to deal with cavitation phenomena at different scales. For example, some studies deal with large-scale cavitation in the Euler framework, while tracking some small-scale bubbles in the Lagrangian framework. This dual-scale processing method can theoretically improve simulation accuracy and reduce computational burden. However, due to the lack of an effective steam structure conversion mechanism, the existing dual-scale method still has limitations when dealing with the conversion of steam structures at different scales. It is difficult to ensure the reasonable conversion of bubbles between different frameworks, which affects the continuity and accuracy of the simulation. Summary of the invention
[0006] In order to solve the technical problem that the existing dual-scale method in the prior art still has limitations when processing the conversion of steam structures of different scales, it is difficult to ensure the reasonable conversion of bubbles between different frames, resulting in the continuity and accuracy of the simulation being affected, the technical solution provided by the present invention is:
[0007] A multi-scale solution method for numerical simulation of cavitation flow field, the method comprising:
[0008] Steps to determine cavitation size;
[0009] The steps to deal with large-scale cavitation in the Eulerian framework include identifying the vapor structure, capturing the vapor-liquid interface, and solving the macroscopic flow field;
[0010] The steps for treating small-scale cavitation in the Lagrangian framework include solving the trajectory, velocity, and pulsation.
[0011] Furthermore, a preferred embodiment is provided, in which a flood filling algorithm identifies the steam structure, a volume fraction method captures the vapor-liquid interface, and a turbulence model solves the macroscopic flow field.
[0012] Furthermore, a preferred embodiment is provided to solve the macroscopic flow field by using RANS and LES turbulence models.
[0013] Furthermore, a preferred embodiment is provided to solve the motion trajectory, velocity and pulsation based on Newton's second law and bubble dynamics equations.
[0014] Furthermore, a preferred embodiment is provided, which also includes the step of collecting a scale threshold. When the cavitation scale is larger than the scale threshold, it is regarded as a large-scale cavitation, and when the cavitation scale is lower than the scale threshold, it is regarded as a small-scale cavitation.
[0015] Furthermore, a preferred embodiment is provided, wherein the large-scale cavitation is recorded as a Lagrangian bubble, and the small-scale cavitation is recorded as an Euler cavitation; when the Lagrangian bubble and the Euler cavitation come into contact, the Lagrangian bubble is converted into a cavity under the Euler framework.
[0016] Based on the same inventive concept, the present invention also provides a multi-scale solution device for numerical simulation of cavitation flow field, the device comprising:
[0017] A module for determining the size of cavitation bubbles;
[0018] A module for handling large-scale cavitation in an Eulerian framework, including identifying vapor structures, capturing vapor-liquid interfaces, and solving macroscopic flow fields;
[0019] A module for processing small-scale cavitation in a Lagrangian framework, including solving for trajectory, velocity, and pulsation.
[0020] Based on the same inventive concept, the present invention also provides a computer storage medium for storing a computer program, and when the computer reads the computer program, the computer executes the described method.
[0021] Based on the same inventive concept, the present invention also provides a computer, including a processor and a storage medium. When the processor reads a computer program stored in the storage medium, the computer executes the method described.
[0022] Based on the same inventive concept, the present invention also provides a computer program product, which is a computer program. When the computer program is executed, the method described above is implemented.
[0023] Compared with the prior art, the technical solution provided by the present invention is beneficial in that:
[0024] The bidirectional coupling method combining the Euler framework and the Lagrangian framework effectively solves the problem that small-scale bubbles in the cavitation flow field cannot be accurately captured by the traditional Euler-Euler method. By processing large-scale cavitation bubbles in the Euler framework and small-scale bubbles in the Lagrangian framework, the evolution of small-scale bubbles can be accurately tracked without increasing the grid resolution while maintaining the macroscopic characteristics of the overall flow field. This multi-scale processing method makes the simulation of cavitation flow more accurate, especially for capturing cavitation noise and cavitation erosion.
[0025] Compared with the traditional Euler-Euler method, the present invention realizes the reasonable conversion of steam structure between different scales through the scale conversion criterion. By evaluating the scale of cavitation and bubbles at each time step, small-scale cavitation can be converted into Lagrangian bubbles, and Lagrangian bubbles with too large scale can be converted into the Euler framework. Such a conversion mechanism ensures the adaptability of steam structures of different scales during the simulation process and improves the flexibility and accuracy of numerical simulation. In particular, compared with the limitations of single framework processing, this method shows significant advantages in capturing the multi-scale cavitation phenomenon of dynamic changes in the flow field.
[0026] The application of the Lagrangian method makes the tracking of small-scale bubbles unrestricted by the grid resolution. Compared with the traditional Euler framework, it avoids the loss of small-scale bubble information due to insufficient grid scale. By solving the motion and pulsation characteristics of bubbles in the external environment through Newton's second law and bubble dynamics equations, the present invention simulates the behavior of bubble groups more accurately. This method is particularly suitable for capturing the details of noise and bubble vibration in cavitation flow, which are difficult to achieve with traditional methods, and therefore has important significance for improving the authenticity of simulation results.
[0027] The present invention uses a flood filling algorithm to identify discrete cavity structures in the Euler framework, so that large-scale and small-scale steam structures can be effectively identified and distinguished. Compared with the traditional identification method based on volume fraction method or other discrete methods, the flood filling algorithm can more accurately distinguish the vapor-liquid interface, especially in the complex cavitation flow field, further enhancing the simulation accuracy and the ability to characterize phenomena of different scales. The high efficiency of this identification method makes the overall simulation computing resource utilization more reasonable and avoids the computing burden caused by over-refining the grid.
[0028] It is suitable for multi-scale solution of numerical simulation of cavitation flow field. BRIEF DESCRIPTION OF THE DRAWINGS
[0029] Figure 1 A schematic diagram of a multi-scale solution method for numerical simulation of cavitation flow fields;
[0030] Figure 2 Schematic diagram of discrete cavity identification algorithm;
[0031] Figure 3 Schematic diagram of the conversion criteria of steam structure between Eulerian and Lagrangian frameworks;
[0032] Figure 4 Schematic diagram of the solution idea of multi-scale method. DETAILED DESCRIPTION
[0033] In order to make the advantages and benefits of the technical solution provided by the present invention more clearly reflected, the technical solution provided by the present invention is now further described in detail in conjunction with the accompanying drawings, specifically:
[0034] Embodiment 1: This embodiment provides a multi-scale solution method for numerical simulation of cavitation flow field, the method comprising:
[0035] Steps to determine cavitation size;
[0036] The steps to deal with large-scale cavitation in the Eulerian framework include identifying the vapor structure, capturing the vapor-liquid interface, and solving the macroscopic flow field;
[0037] The steps for treating small-scale cavitation in the Lagrangian framework include solving the trajectory, velocity, and pulsation.
[0038] Specifically, the specific implementation of this embodiment is based on the numerical simulation of the cavitation flow field, adopts a bidirectionally coupled multi-scale solver, and combines the Euler framework and the Lagrangian framework to process the large-scale cavitation and small-scale bubbles in the cavitation flow field respectively. The specific implementation process first identifies the cavitation under the Euler framework and performs scale evaluation, then calculates and updates the state of the small-scale bubbles in the Lagrangian framework, and combines the scale evaluation criteria to perform structural conversion between different frameworks, and finally completes the numerical simulation of the macroscopic flow field.
[0039] In detail:
[0040] Step 1: Identification and scale assessment of large-scale cavitation in the Euler framework
[0041] First, the large-scale cavitation in the flow field is identified using the Euler framework. Based on the flood filling algorithm, the vapor phase structure in the Euler grid is gradually identified. Through the grid cell index, the information of adjacent cells is obtained to determine the boundary conditions of the vapor-liquid interface, thereby identifying the discrete cavity structure in the Euler framework. After the cavity is identified, its scale is evaluated to determine whether it is suitable for further solution in the Euler framework.
[0042] In the specific implementation, the volume fraction method is used to determine the vapor-liquid distribution in each grid cell, and the number of grid cells occupied by the cavity is calculated. If the number of grid cells occupied by the cavity is lower than the threshold, it is judged as a small-scale cavity and marked as needing to be converted to the Lagrangian framework for processing. The output of this step is the scale information of the cavity and a mark indicating whether the framework conversion is required.
[0043] Step 2: Treatment of small-scale bubbles in the Lagrangian framework
[0044] For the cavities marked as small in step 1, they are converted from the Euler framework to the Lagrangian framework for processing. In the Lagrangian framework, each bubble is regarded as a discrete particle, and the pulsation, velocity and motion trajectory of the bubble under the action of external pressure are calculated based on Newton's second law and the bubble dynamics equation.
[0045] First, the initial state of the bubble is calculated to determine its initial position, volume, and velocity. These initial states are obtained by averaging the small-scale cavities in the Euler grid, including volume, center of mass position, and average velocity. Then, the velocity and trajectory of the bubble are updated based on Newton's second law, taking into account the effects of external pressure, buoyancy, drag, etc. At the same time, the expansion and contraction process of the bubble under pressure changes is described based on the Rayleigh-Plesset equation. The output of this step is the state of the Lagrangian bubble, including its position, velocity, radius, and other parameters.
[0046] Step 3: Lagrangian bubble scale assessment and Eulerian frame conversion
[0047] After processing the bubble in the Lagrangian frame, the size of the bubble needs to be evaluated to determine whether it needs to be converted back to the Eulerian frame. Specifically, if the size of the Lagrangian bubble exceeds a threshold, it is considered too large to continue to be solved as a Lagrangian particle, and it needs to be converted to the Eulerian frame for processing.
[0048] The conversion process first removes the Lagrangian bubbles and then maps their volumes into Eulerian grid cells to maintain mass conservation. At the same time, when the Lagrangian bubbles and Eulerian cavities come into contact with each other, the two will merge and transform the Lagrangian bubbles into cavities in the Eulerian framework to better describe the evolution of the steam structure. The output of this step is the steam structure information in the updated Eulerian framework.
[0049] Step 4: Solving the macroscopic flow field under the Euler framework
[0050] After completing the bidirectional conversion between bubbles and cavities, the macroscopic flow field is solved in the Euler framework. The velocity field, pressure field, and volume fraction field of the flow field are updated using the RANS or LES turbulence model. The VOF method is used to update the volume fraction field to capture the changes in the vapor-liquid interface, and the cavitation model is combined to describe the mass transfer between the phases to ensure that the generation and dissipation of bubbles in the flow field conform to physical laws.
[0051] In the specific implementation, the volume fraction field in the macroscopic flow field is first updated, and the deformation of the vapor-liquid interface is captured based on the VOF method. Then, the velocity field and pressure field are solved, and the turbulence model is used to calculate the motion state of the fluid, considering the influence of large-scale cavitation on the flow field. The output of this step is the updated flow field state, including the volume fraction field, velocity field, and pressure field.
[0052] Step 5: Iterate and solve until convergence
[0053] The above steps are iterated repeatedly until the solution of the flow field reaches the convergence condition. The identification, conversion and solution process of cavitation and bubbles are repeated at each time step to capture the multi-scale dynamic evolution in the cavitation flow field. Through continuous iteration, the steady-state or quasi-steady-state numerical solution of cavitation flow is finally obtained.
[0054] In the implementation, this implementation method relies on the open source CFD platform OpenFOAM, and through secondary development of its two-phase flow solver InterFoam and Lagrangian library, realizes the multi-scale numerical simulation of cavitation flow. Each time step includes the identification and conversion of the Euler cavity, the state update of the Lagrangian bubble, and the solution of the macroscopic flow field until the convergence of the solution meets the preset conditions.
[0055] Implementation method 2: This implementation method further limits the multi-scale solution method for numerical simulation of cavitation flow field provided in implementation method 1. The flood filling algorithm identifies the steam structure, the volume fraction method captures the vapor-liquid interface, and the turbulence model solves the macroscopic flow field.
[0056] Implementation method three: This implementation method further limits the multi-scale solution method for numerical simulation of cavitation flow field provided in implementation method two, and solves the macroscopic flow field through RANS and LES turbulence models.
[0057] Implementation method 4: This implementation method further limits the multi-scale solution method for numerical simulation of cavitation flow field provided in implementation method 1, and solves the motion trajectory, velocity and pulsation based on Newton's second law and bubble dynamics equations.
[0058] Implementation method five. This implementation method is a further limitation of the multi-scale solution method for numerical simulation of cavitation flow field provided in implementation method one, and also includes the step of collecting scale thresholds. When the cavitation scale is larger than the scale threshold, it is regarded as a large-scale cavitation, and when the cavitation scale is lower than the scale threshold, it is regarded as a small-scale cavitation.
[0059] Implementation method six. This implementation method is a further limitation of the multi-scale solution method for numerical simulation of cavitation flow field provided in implementation method five. The large-scale cavitation is recorded as Lagrangian bubble, and the small-scale cavitation is recorded as Euler cavitation. When the Lagrangian bubble and the Euler cavitation come into contact, the Lagrangian bubble is converted into a cavity under the Euler framework.
[0060] Embodiment 7: This embodiment provides a multi-scale solution device for numerical simulation of cavitation flow field, the device comprising:
[0061] A module for determining the size of cavitation bubbles;
[0062] A module for handling large-scale cavitation in an Eulerian framework, including identifying vapor structures, capturing vapor-liquid interfaces, and solving macroscopic flow fields;
[0063] A module for processing small-scale cavitation in a Lagrangian framework, including solving for trajectory, velocity, and pulsation.
[0064] Embodiment 8: This embodiment provides a computer storage medium for storing a computer program. When the computer reads the computer program, the computer executes the method provided in embodiment 1.
[0065] Embodiment 9: This embodiment provides a computer, including a processor and a storage medium. When the processor reads a computer program stored in the storage medium, the computer executes the method provided in embodiment 1.
[0066] Embodiment 10: This embodiment provides a computer program product, which is a computer program. When the computer program is executed, the method provided in embodiment 1 is implemented.
[0067] Implementation Method XI: Combination Figure 1-4 This embodiment is described in detail by using specific examples to further describe the technical solution provided above. Specifically:
[0068] The main problems in the existing technology are: the existing Euler-Euler method is not accurate enough when capturing small-scale bubbles, the Lagrangian method has high computational complexity when dealing with a large number of bubbles or large-scale cavitation, and the existing dual-scale method has deficiencies in the conversion mechanism of bubbles and cavitation, making it difficult to ensure a comprehensive and accurate simulation of cavitation phenomena. These technical problems limit the in-depth understanding of cavitation flow, especially in capturing cavitation noise and cavitation erosion.
[0069] The purpose of this embodiment is to address the shortcomings of the Euler-Euler method commonly used in cavitation flow simulation and to propose a multi-scale cavitation numerical simulation method that can capture both large-scale cavitation and small-scale bubbles in the flow field.
[0070] In order to achieve the above-mentioned object of the invention, the concept of the present embodiment is:
[0071] See also Figure 4 , based on the Euler method to deal with large-scale cavitation in the flow field, the Lagrangian method is used to track small-scale bubbles in the flow field. In addition, the Euler frame and the Lagrangian frame are bidirectionally coupled, and the steam structure is reasonably transformed between the Euler frame and the Lagrangian frame according to the criteria.
[0072] According to the above invention concept, this embodiment adopts the following technical solution:
[0073] For large-scale cavitation, the VOF method is used to capture the vapor-liquid interface, the cavitation model is used to describe the interphase mass transfer, and the RANS or LES series turbulence model is used to solve the macroscopic flow field; for small-scale bubbles, the bubble pulsation under the external environmental pressure is solved based on the bubble dynamics equation, and the bubble velocity and motion trajectory are updated using Newton's second law. In addition, the environmental pressure and velocity required for the calculation process are obtained by interpolation on the macroscopic flow field. Considering that the scale of the steam structure will change during the evolution of the flow field, the scale of the cavitation treated in the Euler framework and the bubble treated in the Lagrangian framework are evaluated at each time step during the calculation. If the scale of a cavitation treated in the Euler framework is too small, it will be converted to the Lagrangian framework; on the contrary, if the scale of a bubble treated in the Lagrangian framework is too large, it will be converted to the Euler framework; in addition, if a Lagrangian bubble and an Euler cavitation contact each other, the two will merge.
[0074] Specifically, the transformation of steam structure between Eulerian framework and Lagrangian framework
[0075] The core of the multi-scale method is to allow cavitation to be reasonably converted between the Eulerian frame and the Lagrangian frame, which mainly includes: the identification of discrete cavities, the transformation criteria of steam structures between frames, and the transformation methods of steam structures between frames. This first requires the identification of discrete cavities in the flow field. In the current work, the flood filling algorithm is used to gradually identify the vapor phase coherent structure, and the grid cell adjacency information is obtained based on the grid cell index. Figure 2 A two-dimensional example is given to illustrate the discrete cavity identification process.
[0076] In multiscale methods, the scale of a steam structure is usually measured based on the number of Eulerian grids it occupies.
[0077] like Figure 3 (a) If a cavity occupies a grid number N bubble Less than a threshold N LE , it is considered that the cavity cannot be effectively resolved by the Euler framework and it will be transformed into the Lagrangian framework; on the contrary, Figure 3 (b) If the number of grids occupied by a Lagrangian bubble is N bubble At a threshold N LE , it is considered that the Lagrangian bubble is too large to be solved based on the seven rounds of spherical bubbles for its motion and pulsation, so it will be converted to the Euler framework; in addition, if Figure 3 As shown in (c), when the Lagrangian bubble comes into contact with the Eulerian cavity, the two will merge with each other, that is, the Lagrangian bubble is transformed into the Eulerian frame.
[0078] Once an Eulerian cavity is too small, it will be transformed into a Lagrangian bubble, whose initial velocity, position and volume are equal to the average velocity, center of mass position and volume of the small-scale cavity, that is:
[0079]
[0080] Among them, V b0 , R b0 、u b0 、x b0 are the initial volume, radius, velocity and spatial coordinates of the Lagrangian bubble, α v,i 、V cell,i 、u cell,i 、x cell,i are the vapor volume fraction, volume, velocity and centroid coordinates of the i-th grid cell occupied by the Euler cavity. In addition, the Lagrangian bubble is in equilibrium at the initial moment of injection into the flow field, so the partial pressure of the non-condensable gas in the Lagrangian bubble p g0 for:
[0081]
[0082] where p c is the ambient pressure at the bubble, p v is the saturated vapor pressure, and σ is the surface tension coefficient. On the contrary, when the Lagrangian bubble with too large scale is transformed into the Eulerian frame, the Lagrangian bubble is first removed from the flow field, and its volume is simultaneously mapped to the Eulerian grid.
[0083] In a specific embodiment, the specific implementation is carried out on the open source CFD platform: OpenFOAM, and the multi-scale calculation of cavitation flow is realized by secondary development of the incompressible two-phase flow solver InterFoam and the Lagrangian library. Figure 1 In the calculation process, this implementation method first identifies the cavitation bubbles in the Euler framework and evaluates their scale. The Euler cavitation bubbles with too small scale will be converted to the Lagrangian framework; then the velocity, position coordinates and radius of the Lagrangian bubbles in the flow field are updated based on Newton's second law and the bubble dynamics equation respectively; then the scale of the Lagrangian bubbles is evaluated, and the ones with too large scale will be converted to the Euler framework; then the macroscopic flow field in the Euler framework is calculated, and the volume fraction field, velocity field and pressure field are updated respectively.
[0084] The technical solution provided by the present invention is further described in detail above through several specific implementation modes in order to highlight the advantages and benefits of the technical solution provided by the present invention. However, the several specific implementation modes described above are not intended to be used as limitations on the present invention. Any reasonable modification and improvement of the present invention, combination of implementation modes and equivalent substitution within the spirit and principle of the present invention shall be included in the protection scope of the present invention.
Claims
1. A multi-scale solution method for numerical simulation of cavitation flow field, characterized in that: The method comprises: Steps to determine cavitation size; The steps to deal with large-scale cavitation in the Eulerian framework include identifying the vapor structure, capturing the vapor-liquid interface, and solving the macroscopic flow field; The steps for treating small-scale cavitation in the Lagrangian framework include solving the trajectory, velocity, and pulsation.
2. A multi-scale solution method for numerical simulation of cavitation flow field according to claim 1, characterized in that: The flood filling algorithm identifies the steam structure, the volume fraction method captures the vapor-liquid interface, and the turbulence model solves the macroscopic flow field.
3. A multi-scale solution method for numerical simulation of cavitation flow field according to claim 2, characterized in that: The macroscopic flow field is solved by RANS and LES turbulence models.
4. The multi-scale solution method for numerical simulation of cavitation flow field according to claim 1, characterized in that: The motion trajectory, velocity and pulsation are solved based on Newton's second law and bubble dynamics equations.
5. The multi-scale solution method for numerical simulation of cavitation flow field according to claim 1, characterized in that: It also includes the step of collecting a scale threshold. When the cavitation scale is greater than the scale threshold, it is regarded as a large-scale cavitation. When the cavitation scale is lower than the scale threshold, it is regarded as a small-scale cavitation.
6. A multi-scale solution method for numerical simulation of cavitation flow field according to claim 5, characterized in that: The large-scale cavitation is recorded as a Lagrangian bubble, and the small-scale cavitation is recorded as an Euler cavitation; when the Lagrangian bubble and the Euler cavitation come into contact, the Lagrangian bubble is converted into a cavity under the Euler framework.
7. A multi-scale solution device for numerical simulation of cavitation flow field, characterized in that: The device comprises: A module for determining the size of cavitation bubbles; A module for handling large-scale cavitation in an Eulerian framework, including identifying vapor structures, capturing vapor-liquid interfaces, and solving macroscopic flow fields; A module for processing small-scale cavitation in a Lagrangian framework, including solving for trajectory, velocity, and pulsation.
8. A computer storage medium for storing a computer program, characterized in that: When the computer reads the computer program, the computer executes the method of claim 1.
9. A computer, comprising a processor and a storage medium, characterized in that: When the processor reads the computer program stored in the storage medium, the computer executes the method of claim 1 .
10. A computer program product, being a computer program, characterized in that When the computer program is executed, the method of claim 1 is implemented.