CN104573241A - RP-3 aviation kerosene cavitation numerical simulation method - Google Patents

Abstract

The invention relates to a RP-3 aviation kerosene cavitation numerical simulation method, relates to the method of reducing the problem for cavitation of RP-3 aviation kerosene during the using process in the industrial area and belongs to the fluid machinery engineering, multi-phase flow and computational fluid mechanics technology. The method uses the cavitation numerical simulation result for guiding the RP-3 aviation kerosene actual application working condition design for reducing the cavitation phenomena in industrial application, and further reducing the vibration, noise and the material surface damage problems caused by cavitation during the storing and transporting process of the RP-3 aviation kerosene. The error between the generated RP-3 aviation kerosene physical replacing model and true RP-3 aviation kerosene known attribute is small, the unknown substance attribute of the true RP-3 aviation kerosene can be compensated, and the method has important practical significance for knowing and analyzing the cavitation character of the RP-3 aviation kerosene.

Description

A kind of method for numerical simulation of RP-3 aviation kerosene cavitation

Technical field

The present invention relates to a kind of method for numerical simulation of RP-3 aviation kerosene cavitation, and relate to the problem brought at industrial circle use procedure cavitation for reducing RP-3 aviation kerosene, belong to fluid machinery engineering, polyphasic flow and Fluid Mechanics Computation technical field.

Background technology

When the liquid internal local pressure of uniform temperature is reduced to hold-up vapour pressure, vaporization phenomenon can be produced, the gas be simultaneously dissolved in liquid also can be separated out, form steam bubble (also known as cavity, hole), when behind the place that steam bubble is higher to pressure with liquid flow movement, steam regelation in bubble, steam bubble is crumbled and fall.Cavity Emergence and Development in this liquid stream, the process crumbled and fall, and consequent series of physical and chemical change process are called cavitation.The generation of cavitation often causes machine efficiency to decline and causes the problems such as vibration, noise and material surface destruction, and the instability causing hydraulic to run and the fatigure failure of blade, can make machine cisco unity malfunction time serious.

Cavitation Problems is the key core problem in fluid machinery engineering and polyphasic flow field always, and due to the complicacy of cavitation phenomenon, method for numerical simulation is one of Main Means of cavitation research always.At present, many Fluid Mechanics Computations (CFD) software can realize the cavitation numerical simulation of Common fluids, applies comparatively extensive.But kerosene is the very complicated potpourri of a kind of composition, the constituent of kerosene is many and mostly be hydrocarbon compound, and its material property is different because of the production place of production, production time and production technology.Cavitation numerical simulation needs many material properties of fluid media (medium) liquid phase and vapour phase, and for RP-3 aviation kerosene, the material property of its liquid phase and vapour phase only has part to be known, then directly cannot be carried out the numerical simulation of kerosene cavitation by cfdrc.Therefore, set up a kind of method for numerical simulation that can realize RP-3 aviation kerosene cavitation, to realize the cavitation numerical simulation of RP-3 aviation kerosene in whole temperature range, have important practical significance.

Existing related art method comprises: the Supertrapp Physical Property Analysis software that National Institute of Standards and Technology (NIST) develops can obtain all kinds of material properties (the Ely J F of conventional hydrocarbon, HuberM L.NIST standard reference datebase 4-NIST thermophysical properties ofhydrocarbon mixtures [M] .National Inst.Of standards, Gaithersburg, MD, 1990.); Under " Groupe Europ é en de Recherches Gazieres " supports, Kunz etc. proposed GERG-2004 equation in 2007, this equation represents with dimensionless Helmholtz free energy equation form based on multi-fluid approximation theory, structure (the Kunz O of potpourri can be realized, Klimeck R, Wagner W, et al.TheGERG-2004Wide-Range Equation of State for Natural Gases and Other Mixtures.Fortschr.-Ber.VDI, Reihe 6, Nr.557, VDI Verlag:D ü sseldorf, 2007.).

Summary of the invention

The technical problem to be solved in the present invention is in the numerical simulation to realize RP-3 aviation kerosene cavitation under the condition of data supporting to the material property of the dirty body medium liquid phase of fixed temperature and vapour phase, and utilize the numerical simulation result of cavitation to reduce the appearance of commercial Application cavitation phenomenon, and then reduce RP-3 aviation kerosene in problems such as the vibration stored and transport process cavitation causes, noise and material surface destructions.The method for numerical simulation of a kind of RP-3 aviation kerosene cavitation disclosed by the invention can improve cavitation simulation precision, realizes the simulation that becomes more meticulous of practical application in industry RP-3 aviation kerosene cavitation processes.Described fluid media (medium) liquid phase and the material property of vapour phase comprise: molal weight, density, specific heat capacity, saturated vapor pressure, specific enthalpy, viscosity, thermal conductivity, thermal expansivity.

The object of the invention is to be achieved through the following technical solutions.

The method for numerical simulation of a kind of RP-3 aviation kerosene cavitation disclosed by the invention, concrete steps comprise:

Step one: the formula that RP-3 aviation kerosene physical replacement model is provided, the molar percentage of each composition is:

N-dodecane: 52 ~ 58%

Decane: 28 ~ 32%

Octane: 1 ~ 3%

Methylcyclohexane: 1 ~ 5%

Toluene: 8 ~ 12%

Step 2: the Supertrapp Physical Property Analysis software adopting National Institute of Standards and Technology (NIST) to develop obtains all kinds of material properties of each composition (pure material) of RP-3 aviation kerosene physical replacement model formula.

Step 3: adopt GERG-2004 equation that each composition (pure material) is mixed into a kind of new virtual substance, realize the attribute of described virtual object qualitative attribution approaching to reality RP-3 aviation kerosene as far as possible, be used for representing RP-3 aviation kerosene potpourri, described virtual substance has complete material property data, and described virtual substance is RP-3 aviation kerosene physical replacement model.The basic structure form of GERG-2004 equation is:

α(δ,τ,x)＝α ⁰(ρ,T,x)+α ^r(δ,τ,x) (1)

${\mathrm{\α}}^0(\mathrm{\ρ},T,x)=\underset{a=1}{\overset{N}{\mathrm{\Σ}}}{x}_a({\mathrm{\α}}_{0a}^0(\mathrm{\ρ},T)+\ln x_a)---\left(2\right)$

${\mathrm{\α}}^r(\mathrm{\δ},\mathrm{\τ},x)=\underset{a=1}{\overset{N}{\mathrm{\Σ}}}{x}_a{\mathrm{\α}}_{0a}^r(\mathrm{\δ},\mathrm{\τ})+\underset{a=1}{\overset{N-1}{\mathrm{\Σ}}}\underset{b=a+1}{\overset{N}{\mathrm{\Σ}}}{x}_a{x}_b{F}_{\mathrm{ab}}{\mathrm{\α}}_{\mathrm{ab}}^r(\mathrm{\δ},\mathrm{\τ})---\left(3\right)$

$\mathrm{\δ}=\frac{\mathrm{\ρ}}{{\mathrm{\ρ}}_r\left(x\right)}---\left(4\right)$

$\mathrm{\τ}=\frac{T_r\left(x\right)}{T}---\left(5\right)$

Wherein subscript a, b represent two components respectively, and, subscript 0 and r original bulk and relative quantity respectively, subscript 0 and r represent desirable item and remainder respectively, α during polycomponent then ⁰for the desirable item of Helmholtz free energy equation, α ^rfor remainder, δ is reduced density, and τ is reduced temperature, F _abfor regulatory factor, ρ _r(x) and T _rx () is respectively density and the temperature funtion of potpourri, the mole fraction of x potpourri shared by each component.Utilize above-mentioned system of equations can by solving by dimensionless Helmholtz free energy equation partial differential the material property obtaining potpourri.As solving of pressure p can be obtained by formula (6) and formula (7):

$\frac{p(\mathrm{\δ},\mathrm{\τ},x)}{\mathrm{\ρRT}}=1+\mathrm{\δ}{\mathrm{\α}}_{\mathrm{\δ}}^r---\left(6\right)$

${\mathrm{\α}}_{\mathrm{\δ}}^r={\left(\frac{\∂{\mathrm{\α}}^r}{\∂\mathrm{\δ}}\right)}_{\mathrm{\τ},x}---\left(7\right)$

Step 4: the material property of RP-3 aviation kerosene physical replacement model is imported commercial CFD code material depot or self-editing CFD program, is defined as a kind of new fluid media (medium), liquid phase and vapour phase data are set respectively.

Step 5: adopt 3D sculpting software to carry out the modeling of flow field regions according to practical situations, adopts stress and strain model software to carry out stress and strain model, grid file is imported commercial CFD code or self-editing CFD program.

Step 6: the boundary condition and initialization condition that calculate basin are arranged according to practical situations.

Step 7: adopt the finite volume method based on finite element to carry out discrete to system of equations, wherein convective term adopts High Accuracy Differencing Scheme, and other adopt central difference schemes, adopts fully implicit solution coupling technique to solving of system of equations.Solving equation group comprises:

Continuity equation:

$\frac{\∂{\mathrm{\ρ}}_m}{\∂t}+\frac{\∂\left({\mathrm{\ρ}}_m{u}_j\right)}{{\∂x}_j}=0---\left(8\right)$

The equation of momentum:

$\frac{\∂\left({\mathrm{\ρ}}_m{u}_i\right)}{\∂t}+\frac{\∂\left({\mathrm{\ρ}}_m{u}_i{u}_j\right)}{{\∂x}_j}=-\frac{\∂p}{\∂x_i}+\frac{\∂}{\∂x_j}\left[(\mathrm{\μ}+{\mathrm{\μ}}_t)(\frac{{\∂u}_i}{{\∂x}_j}+\frac{{\∂u}_j}{{\∂x}_i}-\frac{2}{3}\frac{{\∂u}_i}{{\∂x}_j}{\mathrm{\δ}}_{\mathrm{ij}})\right]---\left(9\right)$

Energy equation:

$\frac{\∂}{\∂t}\left({\mathrm{\ρ}}_m\mathrm{CpT}\right)+\frac{\∂}{\∂x_j}\left({\mathrm{\ρ}}_m{u}_j\mathrm{CpT}\right)=\frac{\∂}{\∂x_j}\left[(\frac{\mathrm{\μ}}{{\mathrm{Pr}}_L}+\frac{{\mathrm{\μ}}_t}{{\mathrm{Pr}}_t})\frac{\∂h}{\∂x_j}\right]-\left\{\frac{\∂}{\∂t}\right[{\mathrm{\ρ}}_m\left(f_{v}L\right)]+\frac{\∂}{{\∂x}_j}[{\mathrm{\ρ}}_m{u}_j\left(f_{v}L\right)\left]\right\}---\left(10\right)$

Mass-conservation equation:

$\frac{\∂{\mathrm{\ρ}}_1{\mathrm{\α}}_1}{\∂t}+\frac{\∂\left({\mathrm{\ρ}}_1{\mathrm{\α}}_1{u}_j\right)}{\∂x_j}={\stackrel{\·}{m}}^{+}+{\stackrel{\·}{m}}^{-}---\left(11\right)$

Evaporation source item:

${\stackrel{\·}{m}}^{+}=F_e\frac{3{\mathrm{\α}}_{\mathrm{nuc}}(1-{\mathrm{\α}}_v){\mathrm{\ρ}}_v}{R_B}\sqrt{\frac{2}{3}\frac{|P_v-P|}{{\mathrm{\ρ}}_l}}---\left(12\right)$

Condensation source item:

${\stackrel{\·}{m}}^{-}=F_c\frac{3{\mathrm{\α}}_v{\mathrm{\ρ}}_v}{R_B}\sqrt{\frac{2}{3}\frac{|P_v-P|}{{\mathrm{\ρ}}_l}}---\left(13\right)$

Turbulent Kinetic k equation:

$\frac{d\left(\mathrm{\ρk}\right)}{\mathrm{dt}}=P_t-\mathrm{\ρ\ϵ}+\frac{\∂}{{\∂x}_j}\left[(\mathrm{\μ}+\frac{{\mathrm{\μ}}_t}{{\mathrm{\σ}}_k})\frac{\∂k}{{\∂x}_j}\right]---\left(14\right)$

Dissipation turbulent kinetic energy ε equation:

$\frac{\mathrm{d\ϵ}}{\mathrm{dt}}=C_{\mathrm{\ϵ}1}\frac{\mathrm{\ϵ}}{k}{P}_t-C_{\mathrm{\ϵ}2}\frac{{\mathrm{\ϵ}}^2}{k}+\frac{\∂}{{\∂x}_j}\left[(\mathrm{\μ}+\frac{{\mathrm{\μ}}_t}{{\mathrm{\σ}}_{\mathrm{\ϵ}}})\frac{\∂\mathrm{\ϵ}}{\∂x_j}\right]---\left(15\right)$

Turbulent viscosity:

${\mathrm{\μ}}_t=\frac{C_{\mathrm{\μ}}{\mathrm{\ρ}}_m{k}^2}{\mathrm{\ϵ}}---\left(16\right)$

Wherein, ρ _m=ρ _lα _l+ ρ _v(1-α _l), u and p is respectively the density of mixed phase, speed and pressure, μ and μ _tbe respectively laminar flow and the turbulent kinetic viscosity of blending agent, f _vfor the massfraction of nitrogen vapor, L is the latent heat of vaporization, Pr _land Pr _tbe respectively the Prandtl number of laminar flow and turbulent flow, h is enthalpy, α _lliquid phase volume mark, in energy equation, last is energy source item. be respectively condensation and evaporation source item, subscript i and j represents coordinate direction respectively, subscript m, l and v represents mixed phase respectively, liquid phase and vapour phase.Employing standard k-ε two equation turbulence model achieves closing of system of equations.

Solver solves stopping after calculating the accuracy requirement reaching given, can obtain the cavitating flows situation of RP-3 aviation kerosene under given geometry and boundary condition by corresponding aftertreatment.Using the numerical simulation of realizing RP-3 aviation kerosene cavitation to the material property of the dirty body medium liquid phase of fixed temperature and vapour phase under the condition of data supporting.

Step 8: utilize the analog result of the method for numerical simulation of a kind of RP-3 aviation kerosene cavitation described in step one to seven to instruct RP-3 aviation kerosene practical application operating condition design, change RP-3 aviation kerosene actual operating mode, reduce the appearance of commercial Application cavitation phenomenon, and then reduce RP-3 aviation kerosene in problem appearance such as the vibration stored and transport process cavitation causes, noise and material surface destructions.Described change RP-3 aviation kerosene actual operating mode mainly refers to the structure changing transport and store, and the temperature of actual motion, speed and pressure.

Beneficial effect:

1, the method for numerical simulation of a kind of RP-3 aviation kerosene cavitation of the present invention, between the known attribute of the RP-3 aviation kerosene physical replacement model generated and true RP-3 aviation kerosene, error is little, can make up the unknown materials attribute of true RP-3 aviation kerosene.

2, the method for numerical simulation of a kind of RP-3 aviation kerosene cavitation of the present invention, can realize the cavitation numerical simulation of RP-3 aviation kerosene in whole temperature range, has important practical significance to the Cavitation Characteristics of understanding and analysis RP-3 aviation kerosene.

3, the method for numerical simulation of a kind of RP-3 aviation kerosene cavitation of the present invention, can be used in the numerical simulation of RP-3 aviation kerosene non-cavitating flow field and cavitating flow in kerosene storage, kerosene transport, kerosene pump fuel feeding, kerosene course of injection, there is important engineer applied and be worth.

4, the method for numerical simulation of a kind of RP-3 aviation kerosene cavitation disclosed by the invention, the numerical simulation result of cavitation is utilized to instruct RP-3 aviation kerosene practical application operating condition design, reduce the appearance of commercial Application cavitation phenomenon, and then reduce RP-3 aviation kerosene in problem appearance such as the vibration stored and transport process cavitation causes, noise and material surface destructions.

Accompanying drawing explanation

Fig. 1 is that physical replacement model and RP-3 aviation kerosene saturated vapor pressure contrast;

Fig. 2 is that physical replacement model and RP-3 aviation kerosene liquid phase kinematic viscosity contrast;

Fig. 3 is that physical replacement model and RP-3 aviation kerosene density of liquid phase contrast;

Fig. 4 is the vapour phase density of physical replacement model formula B;

Fig. 5 is the vapor phase movement viscosity of physical replacement model formula B;

Fig. 6 is the physical replacement model formula liquid phase of B and the specific enthalpy of vapour phase;

Fig. 7 is the physical replacement model formula liquid phase of B and the specific heat capacity of vapour phase;

Fig. 8 is the physical replacement model formula liquid phase of B and the thermal conductivity of vapour phase;

Fig. 9 is the physical replacement model formula liquid phase of B and the thermal expansivity of vapour phase;

Figure 10 is geometry and the boundary condition of numerical simulation in embodiment;

Figure 11 is the numerical simulation result of RP-3 aviation kerosene in embodiment;

Figure 12 is the method for numerical simulation process flow diagram of a kind of RP-3 aviation kerosene of the present invention cavitation.

Embodiment

Be described in further detail the present invention below in conjunction with embodiment, concrete steps comprise:

Step one: the formula providing RP-3 aviation kerosene physical replacement model, table 1 is 5 embodiments of the present invention, and component formula is by mole% meter.

The component formula (molar percentage) of table 1 embodiment

Kerosene is filled a prescription

N-dodecane

Decane

Octane

Methylcyclohexane

Toluene

(mol ratio)	C 12H 26	C 10H 22	C 8H 18	C 7H 14	C 7H 8
						A	58％	32％	1％	1％	8％
B	58％	28％	1％	1％	12％
						C	55％	30％	2％	3％	10％
D	52％	32％	3％	5％	8％
						E	52％	28％	3％	5％	12％

Step 2: adopt the Supertrapp Physical Property Analysis software that National Institute of Standards and Technology (NIST) develops ^[1]obtain the material property of n-dodecane, decane, octane, methylcyclohexane and toluene.

Step 3: adopt GERG-2004 equation that each composition (pure material) is mixed into RP-3 aviation kerosene physical replacement model by formula rate, gained RP-3 aviation kerosene physical replacement model has complete material property data.As shown in Figure 1, Figure 2 and Figure 3, for the contrast of physical replacement model and the known material property of true RP-3 aviation kerosene, can find out that physical replacement model can matching RP-3 aviation kerosene is known well material property, wherein density of liquid phase error is larger, be about 4%, but belong to tolerance interval in engineering.Physical replacement model in described formula range all can substitute RP-3 aviation kerosene preferably and be used for cavitation numerical simulation.In described 5 formulation Example, the fitting precision of formula B is the highest, next carries out step explanation for the B that fills a prescription.Determine that formula B is RP-3 aviation kerosene physical replacement model, extract the material property required for cavitation numerical simulation, comprise: molal weight, density, specific heat capacity, saturated vapor pressure, specific enthalpy, viscosity, thermal conductivity, thermal expansivity, wherein moieties attribute is as shown in Fig. 4 ~ 9.

Step 4: described cavitation numerical simulation adopts large commercial cfdrc ANSYS-CFX14.0 to complete.The material property of RP-3 aviation kerosene physical replacement model is imported the material storehouse of ANSYS-CFX-Pre, a kind of new fluid of definition is RP-3 aviation kerosene liquid, and the another kind of new fluid of definition is RP-3 aviation kerosene gas.

Step 5: adopt ANSYS-ICEM software to carry out modeling to flow field regions as shown in Figure 10, complete stress and strain model simultaneously in ANSYS-ICEM software, grid file is saved as suffix and be called cfx5 file, then grid file is imported ANSYS-CFX-Pre.

Step 6: in ANSYS-CFX-Pre, the boundary condition and initialization condition that calculate basin are arranged, as shown in Figure 10, hydrofoil surface and basin up-and-down boundary are set to without slippage wall, and before and after basin, border is set to the plane of symmetry; Inlet boundary is set to speed entrance, speed of incoming flow U _∞be set as 13.89m/s, reynolds number Re is 7 × 10 ⁵, inlet boundary liquid phase component is set to 1, and vapor-phase composition is set to 0, represents that incoming flow is liquid phase; Outlet border is set to pressure export, and absolute pressure p is set as 59267Pa, and outlet border liquid phase component is set to 1, and vapor-phase composition is set to 0, indicates stream and is liquid phase; Initialize installation arranges in boundary condition and is consistent; Arranging compute type is permanent calculating, and arranging turbulence model is k-ε turbulence model, and arranging cavitation model is Zwart cavitation model, and arranging and calculating the condition of convergence is 1e-4; After being provided with, derive the numerical evaluation file that suffix is def.

Step 7: adopt ANSYS-CFX-Solver Manager software to solve the file that suffix is def, solve after reaching the condition of convergence, solve the file that rear automatic output suffix is called res, described suffix is called res file can carry out aftertreatment with ANSYS-CFX-Post software, flow field data needed for extraction, the RP-3 aviation kerosene obtained for numerical simulation is as shown in figure 11 around the vapour phase volume fraction cloud charts of hydrofoil.

Step 8: utilize the analog result of the method for numerical simulation of a kind of RP-3 aviation kerosene cavitation described in step one to seven to instruct RP-3 aviation kerosene practical application operating condition design, change RP-3 aviation kerosene actual operating mode, reduce the appearance of commercial Application cavitation phenomenon, and then reduce RP-3 aviation kerosene in problem appearance such as the vibration stored and transport process cavitation causes, noise and material surface destructions.Described change RP-3 aviation kerosene actual operating mode mainly refers to the structure changing transport and store, and the temperature of actual motion, speed and pressure.

Claims (2)

1. a method for numerical simulation for RP-3 aviation kerosene cavitation, is characterized in that: concrete steps comprise:

Step one: the formula that RP-3 aviation kerosene physical replacement model is provided, the molar percentage of each composition is:

N-dodecane: 52 ~ 58%

Decane: 28 ~ 32%

Octane: 1 ~ 3%

Methylcyclohexane: 1 ~ 5%

Toluene: 8 ~ 12%

Step 2: the Supertrapp Physical Property Analysis software adopting National Institute of Standards and Technology (NIST) to develop obtains all kinds of material properties of each composition (pure material) of RP-3 aviation kerosene physical replacement model formula;

Step 3: adopt GERG-2004 equation that each composition (pure material) is mixed into a kind of new virtual substance, realize the attribute of described virtual object qualitative attribution approaching to reality RP-3 aviation kerosene as far as possible, be used for representing RP-3 aviation kerosene potpourri, described virtual substance has complete material property data, described virtual substance is RP-3 aviation kerosene physical replacement model, and its basic structure form can be expressed as:

α(δ,τ,x)＝α ⁰(ρ,T,x)+α ^r(δ,τ,x) (1)

${\mathrm{\α}}^0(\mathrm{\ρ},T,x)=\underset{a=1}{\overset{N}{\mathrm{\Σ}}}{x}_a({\mathrm{\α}}_{0a}^0(\mathrm{\ρ},T)+\ln x_a)---\left(2\right)$

${\mathrm{\α}}^r(\mathrm{\δ},\mathrm{\τ},x)=\underset{a=1}{\overset{N}{\mathrm{\Σ}}}{x}_a{\mathrm{\α}}_{0a}^r(\mathrm{\δ},\mathrm{\τ})+\underset{a=1}{\overset{N-1}{\mathrm{\Σ}}}\underset{b=a+1}{\overset{N}{\mathrm{\Σ}}}{x}_a{x}_b{F}_{\mathrm{ab}}{\mathrm{\α}}_{\mathrm{ab}}^r(\mathrm{\δ},\mathrm{\τ})$

$\mathrm{\δ}=\frac{\mathrm{\ρ}}{{\mathrm{\ρ}}_r\left(x\right)}---\left(4\right)$

$\mathrm{\τ}=\frac{T_r\left(x\right)}{T}---\left(5\right)$

Wherein subscript a, b represent two components respectively, and, subscript 0 and r original bulk and relative quantity respectively, subscript 0 and r represent desirable item and remainder respectively, α during polycomponent then ⁰for the desirable item of Helmholtz free energy equation, α ^rfor remainder, δ is reduced density, and τ is reduced temperature, F _abfor regulatory factor, ρ _r(x) and T _rx () is respectively density and the temperature funtion of potpourri, the mole fraction of x potpourri shared by each component; Utilize above-mentioned system of equations can by solving by dimensionless Helmholtz free energy equation partial differential the material property obtaining potpourri; As solving of pressure p can be obtained by formula (6) and formula (7):

$\frac{p(\mathrm{\δ},\mathrm{\τ},x)}{\mathrm{\ρRT}}=1+\mathrm{\δ}{\mathrm{\α}}_{\mathrm{\δ}}^r---\left(6\right)$

${\mathrm{\α}}_{\mathrm{\δ}}^r={\left(\frac{\∂{\mathrm{\α}}^r}{\∂\mathrm{\δ}}\right)}_{\mathrm{\τ},x}---\left(7\right)$

Step 4: the material property of RP-3 aviation kerosene physical replacement model is imported commercial CFD code material depot or self-editing CFD program, is defined as a kind of new fluid media (medium), liquid phase and vapour phase data are set respectively;

Step 5: adopt 3D sculpting software to carry out the modeling of flow field regions according to practical situations, adopts stress and strain model software to carry out stress and strain model, grid file is imported commercial CFD code or self-editing CFD program;

Step 6: the boundary condition and initialization condition that calculate basin are arranged according to practical situations;

Step 7: adopt the finite volume method based on finite element to carry out discrete to system of equations, wherein convective term adopts High Accuracy Differencing Scheme, and other adopt central difference schemes, adopts fully implicit solution coupling technique to solving of system of equations; Solving equation group comprises:

Continuity equation:

$\frac{{\∂\mathrm{\ρ}}_m}{\∂t}+\frac{\∂\left({\mathrm{\ρ}}_m{u}_j\right)}{{\∂x}_j}=0---\left(8\right)$

The equation of momentum:

$\frac{\∂\left({\mathrm{\ρ}}_m{u}_i\right)}{\∂t}+\frac{\∂\left({\mathrm{\ρ}}_m{u}_i{u}_j\right)}{{\∂x}_j}=-\frac{\∂p}{{\∂x}_i}+\frac{\∂}{{\∂x}_j}[(\mathrm{\μ}+{\mathrm{\μ}}_t)(\frac{{\∂u}_i}{{\∂x}_j}+\frac{{\∂u}_j}{{\∂x}_i}-\frac{2}{3}\frac{{\∂u}_i}{{\∂x}_j}{\mathrm{\δ}}_{\mathrm{ij}})---\left(9\right)$

Energy equation:

$\frac{\∂}{\∂t}\left({\mathrm{\ρ}}_m\mathrm{CpT}\right)+\frac{\∂}{{\∂x}_j}\left({\mathrm{\ρ}}_m{u}_j\mathrm{CpT}\right)=\frac{\∂}{{\∂x}_j}\left[(\frac{\mathrm{\μ}}{{\mathrm{Pr}}_L}+\frac{{\mathrm{\μ}}_t}{{\mathrm{Pr}}_t})\frac{\∂h}{{\∂x}_j}\right]-\left\{\frac{\∂}{\∂t}\right[{\mathrm{\ρ}}_m\left(f_{v}L\right)]+\frac{\∂}{{\∂x}_j}[{\mathrm{\ρ}}_m{u}_j\left(f_{v}L\right)\left]\right\}---\left(10\right)$

Mass-conservation equation:

$\frac{\∂{\mathrm{\ρ}}_1{\mathrm{\α}}_1}{\∂t}+\frac{\∂\left({\mathrm{\ρ}}_1{\mathrm{\α}}_1{u}_j\right)}{{\∂x}_j}={\stackrel{.}{m}}^{+}+{\stackrel{.}{m}}^{-}---\left(11\right)$

Evaporation source item:

${\stackrel{.}{m}}^{+}=F_e\frac{3{\mathrm{\α}}_{\mathrm{nuc}}(1-{\mathrm{\α}}_v){\mathrm{\ρ}}_v}{R_B}\sqrt{\frac{2}{3}\frac{|P_v-P}{{\mathrm{\ρ}}_l}}---\left(12\right)$

Condensation source item:

${\stackrel{.}{m}}^{-}=F_c\frac{3{\mathrm{\α}}_v{\mathrm{\ρ}}_v}{R_B}\sqrt{\frac{2}{3}\frac{|P_v-P|}{{\mathrm{\ρ}}_l}}---\left(13\right)$

Turbulent Kinetic k equation:

$\frac{d\left(\mathrm{\ρk}\right)}{\mathrm{dt}}=P_t-\mathrm{\ρ\ϵ}+\frac{\∂}{{\∂x}_j}\left[(\mathrm{\μ}+\frac{{\mathrm{\μ}}_t}{{\mathrm{\σ}}_k})\frac{\∂k}{{\∂x}_j}\right]---\left(14\right)$

Dissipation turbulent kinetic energy ε equation:

$\frac{\mathrm{d\ϵ}}{\mathrm{dt}}=C_{\mathrm{\ϵ}1}\frac{\mathrm{\ϵ}}{k}{P}_t-C_{\mathrm{\ϵ}2}\frac{{\mathrm{\ϵ}}^2}{k}+\frac{\∂}{{\∂x}_j}\left[(\mathrm{\μ}+\frac{{\mathrm{\μ}}_t}{{\mathrm{\σ}}_{\mathrm{\ϵ}}})\frac{\∂\mathrm{\ϵ}}{{\∂x}_j}\right]---\left(15\right)$

Turbulent viscosity:

${\mathrm{\μ}}_t=\frac{C_{\mathrm{\μ}}{\mathrm{\ρ}}_m{k}^2}{\mathrm{\ϵ}}---\left(16\right)$

Wherein, ρ _m=ρ _lα _l+ ρ _v(1-α _l), u and p is respectively the density of mixed phase, speed and pressure, μ and μ _tbe respectively laminar flow and the turbulent kinetic viscosity of blending agent, f _vfor the massfraction of nitrogen vapor, L is the latent heat of vaporization, Pr _land Pr _tbe respectively the Prandtl number of laminar flow and turbulent flow, h is enthalpy, α _lliquid phase volume mark, in energy equation, last is energy source item; be respectively condensation and evaporation source item, subscript i and j represents coordinate direction respectively, subscript m, l and v represents mixed phase respectively, liquid phase and vapour phase; Employing standard k-ε two equation turbulence model achieves closing of system of equations;

Solver solves stopping after calculating the accuracy requirement reaching given, can obtain the cavitating flows situation of RP-3 aviation kerosene under given geometry and boundary condition by corresponding aftertreatment; Using the numerical simulation of realizing RP-3 aviation kerosene cavitation to the material property of the dirty body medium liquid phase of fixed temperature and vapour phase under the condition of data supporting.

2. the method for numerical simulation of a kind of RP-3 aviation kerosene cavitation as claimed in claim 1, it is characterized in that: also comprise step 8: utilize the analog result of the method for numerical simulation of a kind of RP-3 aviation kerosene cavitation described in step one to seven to instruct RP-3 aviation kerosene practical application operating condition design, change RP-3 aviation kerosene actual operating mode, reduce the appearance of commercial Application cavitation phenomenon, and then reduce RP-3 aviation kerosene in problem appearance such as the vibration stored and transport process cavitation causes, noise and material surface destructions; Described change RP-3 aviation kerosene actual operating mode mainly refers to the structure changing transport and store, and the temperature of actual motion, speed and pressure.