Cavitation numerical simulation method based on consideration of multiphase flow shear state
By constructing a numerical simulation method for cavitation under multiphase flow shear state, the problems of insufficient convergence and stability in the calculation of cavitation two-phase flow are solved, and efficient and refined simulation of the cavitation process of a vehicle is realized, which is applicable to axisymmetric vehicles.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-12-26
- Publication Date
- 2026-05-19
AI Technical Summary
Existing computational methods for cavitation two-phase flow have shortcomings in terms of computational convergence and stability, especially in cases where the water vapor interface has a large density ratio in cavitation flow, making it difficult to effectively improve the convergence speed and stability of the computation.
A numerical simulation method for cavitation based on multiphase flow shear states is adopted, including constructing a mathematical model of the preprocessing matrix, turbulence model equations, and flux calculation model of the upwind mode. Combined with the wall mesh construction and flow field initialization of the vehicle body, the existing preprocessing method is improved by constructing an axisymmetric cavitation flow calculation model.
It improves the computational efficiency and flow field refinement of cavitation initiation and evolution processes, ensuring computational accuracy while enhancing computational stability and convergence speed, and is suitable for cavitation simulation of axisymmetric-shaped vehicles.
Smart Images
Figure CN122065707A_ABST
Abstract
Description
Technical Field
[0001] This disclosure relates to the field of multiphase flow numerical computation, and more specifically, to a numerical simulation method based on considering cavitation in multiphase flow shear states. Background Technology
[0002] When a vehicle moves at high speed in water, the local pressure on its surface and even in the wake region drops to the saturated vapor pressure. Water then vaporizes, forming cavitation. Under sufficiently low cavitation numbers, hypercavitation can occur, covering the entire vehicle. Currently, the mainstream method for calculating cavitation two-phase flow is the SIMPLE-type method. The SIMPLE-type method uses the continuity equation and momentum equation to form an approximate pressure correction equation to solve the velocity field, resulting in weak convergence and poor computational stability for high density ratios at the water-vapor interface in cavitation flow.
[0003] Currently, preprocessing methods are also used in the calculation of cavitation two-phase flow. These methods were originally used to calculate low-speed aerodynamic problems. By adding a preprocessing matrix to the imaginary time term, the inviscid flux Jacobian matrix is modified, making the eigenvalues of this matrix more consistent in magnitude under low-speed flow conditions, thus accelerating the calculation convergence speed. Therefore, preprocessing methods are now also used to calculate cavitation two-phase flow problems. However, decoupling the pressure term in cavitation two-phase flow is a challenge. How to construct a reasonable preprocessing matrix to solve the pressure term in two-phase flow while improving the calculation convergence speed is also a difficult problem. Summary of the Invention
[0004] To address the problems existing in the prior art, this disclosure proposes a numerical simulation method based on cavitation considering multiphase flow shear states, to solve at least one of the technical problems listed in the background art. The technical solution adopted in this disclosure is as follows: This disclosure provides a numerical simulation method based on cavitation considering multiphase flow shear state, the method comprising: S100. Initialize the basic parameters of the vehicle body, including its shape, length, underwater speed, and global cavitation number. S200. Construct a mathematical model based on a preprocessing matrix. The mathematical model based on the preprocessing matrix includes a continuity equation, a momentum equation, a component transport equation based on the volume fraction of the cavitation vapor phase, a turbulence model equation, and a cavitation model considering the shear state of multiphase flow. S300, Construct a flux calculation model using an upwind pattern; S400, Construct a mesh for the wall of the vehicle; S500, Initialize the flow field settings for the vehicle body; S600. Based on the basic parameters of the vehicle, the mathematical model based on the preprocessed matrix, the flux calculation model, the mesh constructed for the wall of the vehicle, and the flow field settings initialized for the vehicle, cavitation numerical simulation calculations are performed on the vehicle.
[0005] Preferably, the formula for calculating the global cavitation number is as follows: ; in, The global emptying number; For incoming flow pressure; To disregard the phase change pressure under multiphase flow shear conditions, it can be determined based on fluid properties, which is common knowledge. The density of the liquid; The incoming flow velocity.
[0006] Preferably, the preprocessing matrix takes the following form: ; in, This is a free parameter, with a value range of 5-10; This refers to the volume fraction of the liquid phase. The density of the mixed phase, and ; The density of the liquid phase; This is the density of the vapor phase; This represents the density difference between the liquid and vapor phases.
[0007] Preferably, the continuity equation takes the following form:
[0008] in, For pressure, It is a velocity vector. =50, =0.1; ; ; ; in, To disregard the phase change pressure under multiphase flow shear conditions, it can be determined based on fluid properties. For the incoming flow velocity, For the time scale of the incoming flow, m l It serves as the source phase for the transformation from liquid to vapor phase. m v It serves as the source phase for the transformation of the vapor phase into the liquid phase. s 1 is an axisymmetric source phase.
[0009] Preferably, the momentum equation takes the following form: ; ; ; in, The model feature scale; The value is the liquid phase viscosity, which is 0.001.
[0010] Preferably, the component transport equation based on the volume fraction of the cavitation vapor phase takes the following form: ; ; ; ; in, The viscosity coefficient of the liquid mixture. ; This is the vapor viscosity coefficient; is the turbulent viscosity coefficient.
[0011] Preferably, the turbulence model equations take the following form: ; ; , ; , ; , ; , ; ; ; ; ; .
[0012] Preferably, the cavitation model considering the multiphase flow shear state takes the following form: ; , All are normal stresses; This is the tangential stress.
[0013] During implementation, it is necessary to calculate the normal stress and shear stress at each grid point and determine whether cavitation has formed at each grid point.
[0014] Preferably, the flux calculation model using the upwind pattern is as follows: ; ; in, and Let the variables on the left and right sides of the mesh interface obtained by the selected limiter be the variables, and the eigenvalues of the Jacobian coefficient matrix of the two-phase flow system be: , , , ; ; ; ; .
[0015] Preferably, step S400, which involves constructing a mesh for the wall of the vehicle, specifically includes: constructing a first layer of mesh and a cavitation interface mesh for the wall of the vehicle. Wherein, the normal distance of the first layer of mesh is ≤0.00005 meters, the height growth rate is ≤1.2, and the number of mesh layers in the first layer of mesh that show an increasing relationship is ≤30; the normal distance of the cavitation interface mesh is ≤0.0005 meters.
[0016] During the construction of the cavitation interface mesh, the computational domain can be tested first without considering the influence of the cavitation interface. After the test is completed, the range of the fine mesh is calculated or determined based on the cavitation interface. The mesh within the range is then redrawn to ensure that the mesh normal distance is within 0.0005 meters.
[0017] Preferably, S500 initializes the flow field settings for the vehicle body, specifically including: calculating the physical parameters and incoming flow parameters of each grid point at the initial moment, wherein the incoming flow parameters include incoming flow velocity, incoming flow pressure, liquid phase composition, turbulent kinetic energy, and dissipation rate.
[0018] Preferably, the cavitation numerical simulation calculation of the aircraft specifically includes: S601, will , Set the value to zero and calculate the cavitation-free flow field; S602, Through , The cavitation flow field is calculated using the formula.
[0019] Preferably, in the cavitation numerical simulation calculation, the calculation convergence criterion is that the maximum residual of the liquid phase volume fraction within the calculation domain is within 0.0001.
[0020] The formula for calculating the residual liquid volume fraction at each grid point is as follows: ,in For the first The current liquid phase volume fraction at each grid point For the first The liquid phase volume fraction for each grid point in the next step, with the maximum residual being: .
[0021] The beneficial effects of this disclosure are as follows: This disclosure proposes a numerical simulation method for cavitation based on multiphase flow shear states. By comprehensively considering computational resources and numerical accuracy while taking into account the shape characteristics of the aircraft body, it constructs an axisymmetric cavitation flow calculation model and solution method, thereby improving the computational efficiency and flow field refinement of the initiation and evolution processes of multiphase flow cavitation. This invention, by constructing a cavitation model considering multiphase flow shear states, fully characterizes the cavitation development process while ensuring computational accuracy.
[0022] Compared with the prior art, the present invention has the following beneficial effects: (1). This invention fully considers the numerical simulation method of cavitation in multiphase flow shear state and fully considers the underwater multiphase flow shear effect of the vehicle body, so as to simulate the cavitation initiation phenomenon more accurately.
[0023] (2). Based on the shape characteristics of the aircraft, this invention establishes an axisymmetric mathematical model and performs two-dimensional mesh calculations, thereby saving computing resources and obtaining calculation results more quickly.
[0024] (3). This invention improves existing preprocessing methods by constructing a mathematical model based on a preprocessing matrix, thereby enhancing computational stability and convergence speed compared to traditional preprocessing methods. Attached Figure Description
[0025] The accompanying drawings, which form part of this application, are used to provide a further understanding of this disclosure. The illustrative embodiments of this disclosure and their descriptions are used to explain this disclosure and do not constitute an undue limitation of this disclosure.
[0026] Figure 1 This is a flowchart of a numerical simulation method based on cavitation considering multiphase flow shear state as described in this disclosure.
[0027] Figure 2 This is a schematic diagram of an axisymmetric cavitation multiphase flow resulting from the implementation of this disclosure.
[0028] Figure 3 This is an example diagram of the first layer mesh and the cavitation interface mesh described in this disclosure.
[0029] Figure 4 This is an example diagram of the cavitation-free flow field described in this disclosure.
[0030] Figure 5 for Figure 4 A schematic diagram showing the agreement between the pressure distribution of the corresponding busbar and the test data.
[0031] Figure 6 This is the residual convergence curve for the liquid phase volume fractions described in this disclosure.
[0032] Figure 7 This is an example diagram of the cavitation flow field described in this disclosure. Detailed Implementation
[0033] The present disclosure will now be described in detail with reference to the accompanying drawings and embodiments. It should be noted that, unless otherwise specified, the embodiments and features described in the present application can be combined with each other.
[0034] The following detailed descriptions are exemplary and intended to provide further detailed explanation of this disclosure. Unless otherwise specified, all technical terms used in this disclosure have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains. The terminology used in this disclosure is for the purpose of describing particular embodiments only and is not intended to limit the exemplary embodiments according to this disclosure.
[0035] Accurately simulating the initiation and evolution of underwater cavitation is crucial for evaluating the dynamic characteristics and overall hydrodynamic design of underwater vehicles. This is because cavitation bubbles surrounding the vehicle significantly impact its dynamic characteristics, and significantly reducing drag can dramatically increase its underwater speed. Therefore, this invention proposes a numerical simulation method for cavitation based on multiphase flow shear states. By constructing a cavitation model that considers multiphase flow shear states, it fully characterizes the cavitation development process while ensuring computational accuracy. This invention provides computational means for accurately predicting the underwater dynamic characteristics of vehicles, including a mathematical model based on a preprocessed matrix, encompassing continuity equations, momentum equations, cavitation vapor volume fraction transport equations, turbulence model equations, and a cavitation model considering multiphase flow shear states; analogous to flux calculation methods for upwind scenarios; the first layer of mesh height on the vehicle wall is ≤0.00005 meters, and the mesh thickness at the cavitation interface is ≤0.0005 meters.
[0036] This invention is applicable to the conventional axisymmetric shape of torpedoes. Through schematic diagrams of cavitation multiphase flow in axisymmetric shapes resulting from the implementation of this invention, it can be seen that this invention achieves more thorough capture of cavitation generated by high-shear flow at the bottom of axisymmetric shapes, particularly effective. Figure 2 As shown. This invention constructs a first-layer mesh and a cavitation interface mesh on the wall of the vehicle, forming a two-dimensional mesh mechanism, thereby saving computational resources and quickly obtaining cavitation simulation results.
[0037] like Figure 1 As shown, this disclosure provides a numerical simulation method based on cavitation considering multiphase flow shear state, the method including steps S100 to S600.
[0038] Step S100: Initialize the basic parameters of the vehicle body, including its shape, length, underwater speed, and global cavitation number.
[0039] Furthermore, the formula for calculating the global cavitation number can be as follows: ; in, The global emptying number; For incoming flow pressure; To disregard the phase change pressure under multiphase flow shear conditions, it can be determined based on fluid properties, which is the existing technology. The density of the liquid; The incoming flow velocity.
[0040] Step S200: Construct a mathematical model based on a preprocessing matrix. The mathematical model based on the preprocessing matrix includes a continuity equation, a momentum equation, a component transport equation based on the volume fraction of the cavitation vapor phase, a turbulence model equation, and a cavitation model considering the shear state of multiphase flow.
[0041] Furthermore, the preprocessing matrix takes the following form: ; in, This is a free parameter, with a value range of 5-10; This refers to the volume fraction of the liquid phase. The density of the mixed phase, and ; The density of the liquid phase; This is the density of the vapor phase; This represents the density difference between the liquid and vapor phases.
[0042] Furthermore, the continuity equation takes the following form: ; in, For pressure, It is a velocity vector. =50, =0.1; ; ; ; in, To disregard the phase change pressure under multiphase flow shear conditions, it can be determined based on fluid properties. For the incoming flow velocity, For the time scale of the incoming flow, m l It serves as the source phase for the transformation from liquid to vapor phase. m v It serves as the source phase for the transformation of the vapor phase into the liquid phase. s 1 is an axisymmetric source phase.
[0043] Furthermore, the momentum equation takes the following form: ; ; ; in, The model feature scale; The value is the liquid phase viscosity, which is 0.001.
[0044] Furthermore, the component transport equation based on the volume fraction of the cavitation vapor phase takes the following form: ; ; ; ; in, The viscosity coefficient of the liquid mixture. ; This is the vapor viscosity coefficient; is the turbulent viscosity coefficient.
[0045] Furthermore, the turbulence model equations are in the following form: ; ; , ; , ; , ; , ; ; ; ; ; .
[0046] Furthermore, the phase transformation pressure considering the multiphase flow shear state can be determined using the generalized Newtonian stress formula, specifically as follows: Step (1): Construct the generalized Newtonian stress formula as follows: ; ; in, The viscosity coefficient of the liquid; Step (2), and Substituting 0.001 into the generalized Newtonian stress formula, the pressure at each point is divided into normal stress and tangential stress, and the criteria for liquid cavitation are set as follows: ; Step (3): The phase change pressure considering the multiphase flow shear state is calculated as follows: .
[0047] Furthermore, the cavitation model considering the multiphase flow shear state can take the following form: ; , All are normal stresses; This is the tangential stress.
[0048] During implementation, it is necessary to calculate the normal stress and shear stress at each grid point and determine whether cavitation has formed at each grid point.
[0049] Step S300: Construct a flux calculation model using an upwind approach.
[0050] Furthermore, the flux calculation model using the upwind pattern is as follows: ; ; in, and Let the variables on the left and right sides of the mesh interface obtained by the selected limiter be the variables, and the eigenvalues of the Jacobian coefficient matrix of the two-phase flow system be: , , , ; ; ; ; .
[0051] Step S400: Construct a mesh for the wall of the vehicle.
[0052] Furthermore, the step of constructing a mesh for the wall of the vehicle specifically includes: constructing a first layer mesh and a cavitation interface mesh for the wall of the vehicle; Wherein, the normal distance of the first layer of mesh is ≤0.00005 meters, the height growth rate is ≤1.2, and the number of mesh layers in the first layer of mesh that show an increasing relationship is ≤30; the normal distance of the cavitation interface mesh is ≤0.0005 meters.
[0053] Furthermore, during the construction of the cavitation interface mesh, the computational domain can be tested first without considering the influence of the cavitation interface. After the test is completed, the range of the finer mesh is calculated or determined based on the cavitation interface, and the mesh within the range is redrawn to ensure that the mesh normal distance is within 0.0005 meters. Step S500: Initialize the flow field settings for the vehicle.
[0054] Furthermore, the initial flow field settings for the vehicle body specifically include: calculating the physical parameters and incoming flow parameters of each grid point at the initial moment, wherein the incoming flow parameters include incoming flow velocity, incoming flow pressure, liquid phase composition, turbulent kinetic energy, and dissipation rate.
[0055] Step S600: Based on the basic parameters of the vehicle, the mathematical model based on the preprocessing matrix, the flux calculation model, the mesh constructed for the wall of the vehicle, and the flow field settings initialized for the vehicle, perform cavitation numerical simulation calculations on the vehicle.
[0056] Furthermore, the cavitation numerical simulation calculation of the aircraft specifically includes: S601, will , Set the value to zero and calculate the cavitation-free flow field; S602, Through , The cavitation flow field is calculated using the formula.
[0057] Furthermore, in the cavitation numerical simulation calculation, the convergence criterion is that the maximum residual of the liquid phase volume fraction within the computational domain is within 0.0001. The determination step is as follows: the formula for calculating the residual of the liquid phase volume fraction at each grid point is... ,in For the first Current liquid phase volume fraction at each grid point For the first The next liquid phase volume fraction for each grid point, with the maximum residual being: .
[0058] To better understand the technical principles of this disclosure, the numerical simulation method based on cavitation considering multiphase flow shear state provided by this invention can be understood through the following example steps.
[0059] Step (100): Initialize the basic parameters of the vehicle body, specifically: set the shape of the vehicle body to a hemispherical cylinder, the length to 1m, the underwater speed to 10m / s, and the global cavitation number to 0.3.
[0060] Step (200): Construct a mathematical model based on a preprocessing matrix. The mathematical model based on the preprocessing matrix includes a continuity equation, a momentum equation, a component transport equation based on the volume fraction of the cavitation vapor phase, a turbulence model equation, and a cavitation model considering the shear state of multiphase flow.
[0061] Furthermore, the preprocessing matrix takes the following form: ; In the formula This is a free parameter, with a value of 10. This refers to the volume fraction of the liquid phase. For the mixed phase density, , =1000kg / m 3 , =0.59kg / m 3 , 999.41 kg / m 3 .
[0062] Furthermore, the continuity equation takes the following form: ; in, For pressure, It is a velocity vector. =50, =0.1.
[0063] ; ; ; in, =86325Pa, =10m / s, =0.1s, m l It serves as the source phase for the transformation from liquid to vapor phase. m v It serves as the source phase for the transformation of the vapor phase into the liquid phase. s 1 is an axisymmetric source phase.
[0064] Furthermore, the momentum equation takes the following form: ; ; ; In the formula For model feature scale, The viscosity of the liquid phase is taken as 0.001.
[0065] Furthermore, the component transport equation based on the volume fraction of the cavitation vapor phase takes the following form: ; ; ; ; in, The viscosity coefficient of the liquid mixture. ; This is the vapor viscosity coefficient; is the turbulent viscosity coefficient.
[0066] Furthermore, the turbulence model equations are taken in the following form: ; ; , ; , ; , ; , ; ; ; ; ; ; .
[0067] Furthermore, the component transport equation based on the volume fraction of the cavitation vapor phase is adopted in the following form: ; ; ; ; in, The viscosity coefficient of the liquid mixture. ; This is the vapor viscosity coefficient; is the turbulent viscosity coefficient.
[0068] Furthermore, the phase transformation pressure considering the multiphase flow shear state can be determined using the generalized Newtonian stress formula, specifically: Step (1): Construct the generalized Newtonian stress formula as follows: ; ; in, The viscosity coefficient of the liquid; Step (2), and Substituting 0.001 into the generalized Newtonian stress formula, the pressure at each point is divided into normal stress and tangential stress, and the criteria for liquid cavitation are set as follows: ; Step (3): The phase change pressure considering the multiphase flow shear state is obtained by calculation: .
[0069] Furthermore, the cavitation model considering the multiphase flow shear state can take the following form: ; , All are normal stresses; This is the tangential stress.
[0070] During implementation, the normal stress and shear stress of each grid point are calculated, and it is determined whether cavitation has formed at each grid point.
[0071] Step (300): Construct a flux calculation model using an upwind pattern, wherein the upwind flux calculation model takes the following form: ; ; in, and Let the variables on the left and right sides of the mesh interface obtained by the selected limiter be the variables, and the eigenvalues of the Jacobian coefficient matrix of the two-phase flow system be: , , , ; ; ; ; .
[0072] Step (400): Construct a mesh for the wall of the vehicle, specifically including: constructing a first layer mesh and a cavitation interface mesh for the wall of the vehicle.
[0073] like Figure 3 As shown, the normal distance of the first layer mesh is set to 0.00005 meters, the height growth rate is ≤1.2, and the number of mesh layers in the first layer mesh that show an increasing relationship is 20; the normal distance of the cavitation interface mesh is set to 0.0005 meters.
[0074] Step (500): Initialize the flow field settings for the vehicle body, specifically as follows: The physical parameters of each grid point at the initial moment are equal to the incoming flow parameters, and the incoming flow velocity is... =10m / s, number of liquid phase components turbulent kinetic energy Dissipation rate Incoming pressure =101325.
[0075] Step (600): Based on the basic parameters of the vehicle, the mathematical model based on the preprocessing matrix, the flux calculation model, and the cavitation numerical simulation calculation of the vehicle, perform cavitation numerical simulation calculation.
[0076] Furthermore, the cavitation numerical simulation calculation of the aircraft specifically includes: Step (601) will , Set the value to zero and calculate the cavitation-free flow field. After 3000 calculation steps, the cavitation-free flow field is obtained, as shown below. Figure 4 As shown; the pressure distribution of the busbar agrees well with the test data, such as Figure 5 As shown.
[0077] Step (602) is passed , The cavitation flow field was calculated using the formula. The maximum residual for the liquid phase volume fraction was within 0.0001, and the calculation converged after 4000 steps (e.g., ...). Figure 6 As shown), thus obtaining the cavitation flow field (as shown). Figure 7 (As shown). The global cavitation number is 0.3, and the cavitation flow field exhibits a cavitation bubble morphology.
[0078] In summary, this disclosure proposes a numerical simulation method for cavitation based on multiphase flow shear states. By comprehensively considering computational resources and numerical accuracy while taking into account the shape characteristics of the aircraft body, it constructs an axisymmetric cavitation flow calculation model and solution method, thereby improving the computational efficiency and flow field refinement of the initiation and evolution processes of multiphase flow cavitation. This invention, by constructing a cavitation model considering multiphase flow shear states, fully characterizes the cavitation development process while ensuring computational accuracy.
[0079] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this disclosure and not to limit them. Although this disclosure has been described in detail with reference to the above embodiments, those skilled in the art should understand that modifications or equivalent substitutions can still be made to the specific implementation of this disclosure. Any modifications or equivalent substitutions that do not depart from the spirit and scope of this disclosure should be covered within the protection scope of the claims of this disclosure.
Claims
1. A numerical simulation method based on cavitation considering multiphase flow shear state, characterized in that, The method includes: S100. Initialize the basic parameters of the vehicle body, including its shape, length, underwater speed and global cavitation number; S200. Construct a mathematical model based on a preprocessing matrix. The mathematical model based on the preprocessing matrix includes a continuity equation, a momentum equation, a component transport equation based on the volume fraction of the cavitation vapor phase, a turbulence model equation, and a cavitation model considering the shear state of multiphase flow. S300, Construct a flux calculation model using an upwind pattern; S400, Construct a mesh for the wall of the vehicle; S500. Initialize the flow field settings for the vehicle body, specifically including: calculating the physical parameters and incoming flow parameters of each grid point at the initial moment, wherein the incoming flow parameters include incoming flow velocity, incoming flow pressure, liquid phase composition, turbulent kinetic energy and dissipation rate; S600. Based on the basic parameters of the vehicle, the mathematical model based on the preprocessed matrix, the flux calculation model, the mesh constructed for the wall of the vehicle, and the flow field settings initialized for the vehicle, cavitation numerical simulation calculations are performed on the vehicle. The formula for calculating the global cavitation number is as follows: ； in, The global emptying number; For incoming flow pressure; Phase change pressure without considering the shear state of multiphase flow; The density of the liquid; The incoming flow velocity.
2. The numerical simulation method based on cavitation considering multiphase flow shear state as described in claim 1, characterized in that, The preprocessing matrix takes the following form: ； in, This is a free parameter, with a value range of 5-10; This refers to the volume fraction of the liquid phase. The density of the mixed phase, and ; The density of the liquid phase; This is the density of the vapor phase; This represents the density difference between the liquid and vapor phases.
3. The numerical simulation method based on cavitation considering multiphase flow shear state as described in claim 1, characterized in that, The continuity equation takes the following form: ； in, For pressure, It is a velocity vector. =50, =0.1; ； ； ； in, Phase change pressure without considering the shear state of multiphase flow; For the incoming flow velocity, For the time scale of the incoming flow, m l It serves as the source phase for the transformation from liquid to vapor phase. m v It serves as the source phase for the transformation of the vapor phase into the liquid phase. s 1 is an axisymmetric source phase.
4. The numerical simulation method based on cavitation considering multiphase flow shear state as described in claim 1, characterized in that, The momentum equation takes the following form: ； ； ； in, The model feature scale; The value is the liquid phase viscosity, which is 0.
001.
5. The numerical simulation method based on cavitation considering multiphase flow shear state as described in claim 1, characterized in that, The component transport equation based on the volume fraction of cavitation vapor phase takes the following form: ； ； ； ； in, The viscosity coefficient of the liquid mixture. ; This is the vapor viscosity coefficient; is the turbulent viscosity coefficient.
6. The numerical simulation method based on cavitation considering multiphase flow shear state as described in claim 1, characterized in that, The turbulence model equations take the following form: ； ； ， ； ， ； ， ； ， ； ； ； ； ； 。 7. The numerical simulation method based on cavitation considering multiphase flow shear state as described in claim 1, characterized in that, The cavitation model considering the multiphase flow shear state takes the following form: ； , All are normal stresses; This is the tangential stress.
8. The numerical simulation method based on cavitation considering multiphase flow shear state as described in claim 1, characterized in that, The flux calculation model using the upwind pattern is as follows: ； ； in, and Let the variables on the left and right sides of the mesh interface obtained by the selected limiter be the variables, and the eigenvalues of the Jacobian coefficient matrix of the two-phase flow system be: ， ， ， ； ； ； ； 。 9. The numerical simulation method based on cavitation considering multiphase flow shear state as described in claim 1, characterized in that, S400, constructing a mesh for the wall of the vehicle, specifically includes: constructing a first layer mesh and a cavitation interface mesh for the wall of the vehicle; Wherein, the normal distance of the first layer of mesh is ≤0.00005 meters, the height growth rate is ≤1.2, and the number of mesh layers in the first layer of mesh that show an increasing relationship is ≤30; the normal distance of the cavitation interface mesh is ≤0.0005 meters.
10. The numerical simulation method based on cavitation considering multiphase flow shear state as described in claim 3, characterized in that, The cavitation numerical simulation calculation of the aircraft specifically includes: S601, will , Set the value to zero and calculate the cavitation-free flow field; S602, Through , The cavitation flow field is calculated using the formula. In the cavitation numerical simulation calculation, the convergence criterion is that the maximum residual of the liquid phase volume fraction within the computational domain is within 0.0001. The determination steps are as follows: the formula for calculating the residual of the liquid phase volume fraction at each grid point is... ,in For the first Current liquid phase volume fraction at each grid point For the first The next liquid phase volume fraction for each grid point, with the maximum residual being: 。