CN107784161A - A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed - Google Patents
A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed Download PDFInfo
- Publication number
- CN107784161A CN107784161A CN201710887696.8A CN201710887696A CN107784161A CN 107784161 A CN107784161 A CN 107784161A CN 201710887696 A CN201710887696 A CN 201710887696A CN 107784161 A CN107784161 A CN 107784161A
- Authority
- CN
- China
- Prior art keywords
- compressible
- water
- high speed
- speed
- dynamic characteristic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed disclosed by the invention, belongs to fluid machinery engineering field.The present invention comprises the following steps：Step 1：Establish geometrical model；Step 2：The geometrical model established to step 1 and fluid calculation domain carry out mesh generation；Step 3：Boundary condition and computational methods are set；Step 4：The compressible amendment of the material property of water is realized by carrying out secondary development to CFD software；Step 5：Carry out the numerical Analysis of the compressible supercavity flow dynamic characteristic of high speed, that is, realize a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed.Also include step 6：Numerical Simulation Analysis is carried out, numerical simulation result and experimental result are contrasted, verifies the accuracy of analysis method, feasibility.The present invention can disclose the compressible supercavitating flow rule of high speed and mechanism, be provided fundamental basis so as to navigate by water the structure design of body for underwater high-speed, and can solve the problem that underwater high-speed navigates by water body practical application engineering problem.
Description
Technical field
The invention belongs to fluid machinery engineering field, is related to a kind of compressible supercavitating flow characteristic analysis method of high speed,
A kind of more particularly to compressible supercavitating flow characteristic analysis method of high speed based on CFD software secondary development.
Background technology
Cavitation refers to a kind of phase transition phenomena under certain condition, occurred in liquid phase medium.When the local pressure drop in water
During low saturated vapor pressure to water, part water will become vapor, and so as to form the hole of steam in current, i.e., there occurs sky
Change.Flox condition is different, and the form of cavitation bubble also has very big difference.It is broadly divided into：Nascent vacuole, partial cavity and supercavity.
Whole object plane is covered when vacuole area develops into, and is extended to behind object, or even cavity dimensions can also be considerably beyond object chi
Very little, here it is supercavity.
Sail body supercavity flow in underwater exercise has two types, and a kind of is that nature vaporizes the natural supercavitation to be formed
Stream, another kind are the ventilated supercavitation stream that artificial ventilation is formed.No matter any form, be all to be subtracted using supercavity
Resistance, its high speed supercavity underwater high-speed navigation body is exactly a kind of underwater weapons using natural supercavitation drag reduction technology.
Different from conventional sail body underwater navigation problem, supercavity underwater high-speed navigates by water body because headway is very high, up to across sound
Speed even supersonic speed.It is related to the compressible supercavitating flow problem of complexity in motion process, including it is vapour-liquid compressibility, unsteady
Flowing, water shock wave, phase transformation and strong turbulence etc. hydrodynamics problem.At present, for these problems also without more accurately
Analysis means and numerical computation method, it is necessary to study a kind of number for being capable of the compressible supercavity flow dynamic characteristic of accurate simulation high speed
Value calculating method, and be applied under the operating modes such as subsonic speed, transonic speed and supersonic speed, to disclose the compressible supercavity flow of high speed
Dynamic rule and mechanism, provide fundamental basis so as to navigate by water the structure design of body for underwater high-speed.
The content of the invention
Disclosed by the invention to be searched regarding to the issue above and difficult point based on grid, a kind of high speed disclosed by the invention is compressible super
The analysis method technical problems to be solved of vacuole flow behavior be to disclose the compressible supercavitating flow rule of high speed and mechanism,
Provided fundamental basis so as to navigate by water the structure design of body for underwater high-speed, and can solve the problem that underwater high-speed navigates by water body practical application work
Cheng Wenti.
Described underwater high-speed navigation body practical application engineering problem includes underwater high-speed navigation body navigation stability, navigation
Deposit speed.
The purpose of the present invention is achieved through the following technical solutions.
A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed disclosed by the invention, comprises the following steps：
Step 1：Establish geometrical model.
Body is navigated by water as research object using the underwater high-speed of rotary structure, the geometric parameter of revolving body is determined, using modeling
Software carries out Geometric Modeling.
Described modeling software preferred three-dimensional modeling software, described 3D sculpting software include Solidworks or UG.
Step 2：The geometrical model established to step 1 and fluid calculation domain carry out mesh generation.
The underwater high-speed obtained in step 1 navigation body Model is exported into curved surface to mesh generation software, and in mesh generation
Fluid calculation domain is established in software, piecemeal watershed carries out mesh generation, and carries out grid independence checking, obtains Bestgrid
Number.
Step 2 concrete methods of realizing comprises the following steps：
Step 2.1：The underwater high-speed obtained in step 1 navigation body Model is exported into curved surface to mesh generation software, and
Establish sufficiently large fluid calculation domain in mesh generation software, sufficiently large fluid calculation domain refers to big to can ignore that Basin Boundary
Influence to flow field around underwater high-speed navigation body, then navigates by water body to underwater high-speed and each face in basin is defined；
Step 2.2：By basin piecemeal partition structure grid, consider wall effect to grid around underwater high-speed navigation body
Re-encryption, in order to ensure computational accuracy, it is necessary to assure y+Within zone of reasonableness；
Step 2.3：Grid independence checking is carried out, using resistance coefficient and vacuole profile dual indexes as validation criteria,
Take the basin of multigroup different grid numbers to carry out the calculating of identical operating mode, for result of calculation take vacuole profile in steady flow condition and
Resistance coefficient, grid reaches computational accuracy when front and rear prediction result twice is less than preset value with experimental result difference, now for
Optimal value, that is, obtain Bestgrid number.
Y described in step 2.2+=yuτ/νlTo weigh the index of grid precision, wherein y is the thickness of first layer grid,
uτFor wall friction speed, νlFor the dynamic viscosity of water；Described zone of reasonableness is for underwater high-speed navigation body preferably 30≤y+≤
60。
Preset value preferably 5% described in step 2.3.
Mesh generation software described in step 2 includes ICEM or POINTWISE softwares.
Step 3：Boundary condition and computational methods are set.
Grid file is imported in FLUENT, boundary condition, the computation model used, computational methods are configured.
Step 3 concrete methods of realizing comprises the following steps：
Step 3.1：Grid file is imported in FLUENT.
Step 3.2：Boundary condition is set.
Basin condition for import uses speed condition for import；Exit condition uses pressure export condition, and pressure considers ambient water
It is deep；Remaining is using without sliding wall boundary condition.
Step 3.3：The computation model used is set.
Cavitating flows calculating is carried out based on N-S equations, VOF interface capturing methods are used under homogeneous phase model, opens energy
Equation is measured, turbulence model uses SST k- ω models, and cavitation model uses Zwart models.
Step 3.4：Computational methods are set.
Body flow field is navigated by water to underwater high-speed using the non-fixed length solver of double precision and carries out numerical simulation, from PISO algorithms,
Pressure difference value form selects PRESTO！, the discrete method of the equation of momentum, energy equation and Equations of Turbulence inserted from Second-order Up-wind
It is worth form, CICSAM difference forms are selected in the calculating of volume fraction.
Step 4：The compressible amendment of the material property of water is realized by carrying out secondary development to CFD software.
Aqueous water and vapourous water is selected to consider aqueous water in high speed as principal phase as fluid working substance, wherein aqueous water
Compressive characteristics, according to the Tait state equation compiling user custom formula UDF of compressible water, the density attributes of water are carried out
Amendment, and define effective sound velocity in water；Steam-like water is arranged to the second phase, and steam-like water substance attribute is arranged to preferable
Gas.
Step 4 concrete methods of realizing comprises the following steps：
Step 4.1：Aqueous water is selected in Fluent as principal phase working medium；Usual water body is incompressible, but works as water
Lower high speed operation body headway is then considered as compressible fluid close to water body during velocity of sound, from Tait equations as compressible liquid
The state equation of body, the Tait equations for pressing water for being free from temperature influence below are：
Wherein, p∞、ρ∞And a∞For reference pressure, reference density and the velocity of sound of aqueous water at stationary stream field, k and
For attribute coefficients, the value of these variables is depending on actual condition；
Preferred p at normal temperatures and pressures∞=101325Pa, ρ∞=998kg/m3, a∞=1430m/s, k=7.15,
Containing temperature influence the Tait equations for pressing water be：
Wherein, saturation pressure and density are calculated as the function of temperature with Oldenbourg coefficients；
Wherein, θ=1-T/Tc, pc、Tc,ρcIt is critical condition (2.264 × 107Pa, 647.14K, the 332kg/m of water3)；It is public
Coefficient a in formula (3), (4)i(i=1~6), bi(i=1~6) are respectively constant；Wherein, ai=(- 7.85823,
1.83991, -11.7811,22.6705, -15.9393,1.77516)；bi=(1.99206,1.10123, -0.512506, -
1.75263, -45.4485, -6.75615E+05)；
Supercavitation process is related to vehicle repair major mixing, it is necessary to define the effective sound velocity of mixed phase；It is related to vehicle repair major mixing
The expression formula of problem effective sound velocity：
Wherein, αvFor steam volume fraction, ρmTo mix phase density, ρlAnd ρvThe density of respectively pure liquid phase and vapour phase water,
clFor the pure liquid phase velocity of sound, cvFor the velocity of sound of pure gas phase；ρl、ρv、clAnd cvValue depending on actual condition；
Preferred ρ at normal temperatures and pressuresl=998kg/m3, ρv=0.5542kg/m3, cl=1430m/s, cv=429m/s；
The state equation of compressible water and effective sound velocity are compiled into Fluent softwares by UDF, realized to material category
The amendment of property；
Step 4.2：Steam-like water is selected in Fluent as the second phase working medium；Steam-like water substance attribute is arranged to
Perfect gas.
Step 5：The numerical computations of the compressible supercavity flow dynamic characteristic of high speed are carried out in setting based on step 1 to step 4
Analysis, that is, realize a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed.
Also include step 6：Utilize a kind of point of the compressible supercavity flow dynamic characteristic of high speed described in step 1 to step 5
Analysis method carries out Numerical Simulation Analysis to the compressible supercavitating flow process of high speed, by numerical simulation result and experimental result pair
Than verifying the accuracy of the analysis method of the compressible supercavity flow dynamic characteristic of high speed a kind of, feasibility.
Step 6 concrete methods of realizing comprises the following steps：
Step 6.1：The result of calculation in Fluent is post-processed using the poster processing soft, obtains pressure, density, speed
The information of flow such as degree, volume fraction cloud atlas and quantitative distribution situation；
Step 6.2：Step 6.1 Numerical results and experimental result or Classic Experiments data are contrasted, checking is tested
Demonstrate,prove the accuracy of numerical computation method, feasibility.
Also include step 7：Utilize a kind of point of the compressible supercavity flow dynamic characteristic of high speed described in step 1 to step 5
Analysis method carries out Numerical Simulation Analysis to the compressible supercavitating flow process of high speed, is advised with disclosing the compressible supercavitating flow of high speed
Rule and mechanism, provided fundamental basis so as to navigate by water the structure design of body for underwater high-speed, and can solve the problem that underwater high-speed navigates by water body
Practical application engineering problem.
Beneficial effect：
1st, a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed disclosed by the invention, it is contemplated that water during high speed
Compressibility, by the secondary development to Fluent softwares, liquid phase substance attribute is modified, Fluent softwares can be broken through
In can only single-phase compressible limitation, result of calculation is more accurate.
2nd, the analysis method of the compressible supercavity flow dynamic characteristic of a kind of high speed disclosed by the invention, by result and Classic Experiments
Data are contrasted, and demonstrate the accuracy of numerical computation method, feasibility.
3rd, the analysis method of the compressible supercavity flow dynamic characteristic of a kind of high speed disclosed by the invention, super sky compressible to high speed
Steep flow process and carry out Numerical Simulation Analysis, to disclose the compressible supercavitating flow rule of high speed and mechanism, so as to be underwater high
The structure design of fast sail body is provided fundamental basis, and can solve the problem that underwater high-speed navigates by water body practical application engineering problem.
Brief description of the drawings
A kind of flow chart of the analysis method of the compressible supercavity flow dynamic characteristic of high speed of Fig. 1 present invention；
Underwater high-speed navigates by water volume geometric model schematic diagram used in Fig. 2 embodiment of the present invention；
Mesh generation schematic diagram in Fig. 3 embodiment of the present invention；
Boundary condition sets schematic diagram in Fig. 4 embodiment of the present invention；
Numerical method Accuracy Verification figure in Fig. 5 embodiment of the present invention；
Numerical Simulation Results density cloud atlas under different Mach number in Fig. 6 embodiment of the present invention；
Numerical Simulation Results pressure cloud atlas under different Mach number in Fig. 7 embodiment of the present invention；
Numerical Simulation Results vacuole contoured profile figure in Fig. 8 embodiment of the present invention.
Embodiment
In order to better illustrate a kind of numerical computation method of the compressible supercavity flow dynamic characteristic of high speed of the present invention,
Calculating is tested in conjunction with the accompanying drawings and embodiments using the method for the present invention so that technical scheme and beneficial effect are more clear
Chu.
Embodiment 1：
The present embodiment navigates by water body as research object using underwater high-speed disclosed in foreign countries, underwater high-speed navigation body total length
The slenderness ratio of 157.4mm, head cavitator diameter 1.42mm, bullet length and underwater high-speed navigation body maximum gauge is 12.
As shown in figure 1, the flow chart of the analysis method for a kind of compressible supercavity flow dynamic characteristic of high speed of the present embodiment,
A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed, specific implementation step are as follows disclosed in the present embodiment：
Step 1：Establish geometrical model.
Body is navigated by water as research object using the underwater high-speed of rotary structure, it is determined that the geometric parameter of revolving body, using building
Mould software has carried out Geometric Modeling, shows as shown in Fig. 2 navigating by water volume geometric model for underwater high-speed used in the embodiment of the present invention
It is intended to.
Described modeling software selects 3D sculpting software Solidworks.
Step 2：The geometrical model established to step 1 and fluid calculation domain carry out mesh generation.
Step 2 concrete methods of realizing comprises the following steps：
Step 2.1：The underwater high-speed obtained in step 1 navigation body Model is exported into curved surface to mesh generation software, and
Establish sufficiently large fluid calculation domain in mesh generation software, sufficiently large fluid calculation domain refers to big to can ignore that Basin Boundary
Influence to flow field around underwater high-speed navigation body, then navigates by water body to underwater high-speed and each face in basin is defined；
Step 2.2：By basin piecemeal partition structure grid, consider wall effect to grid around underwater high-speed navigation body
Re-encryption, in order to ensure computational accuracy, it is necessary to assure y+Within zone of reasonableness；
Step 2.3：Grid independence checking is carried out, using resistance coefficient and vacuole profile dual indexes as validation criteria,
Take the basin of multigroup different grid numbers to carry out the calculating of identical operating mode, for result of calculation take vacuole profile in steady flow condition and
Resistance coefficient, grid reaches computational accuracy when front and rear prediction result twice is less than preset value 5% with experimental result difference, now
For optimal value, that is, obtain Bestgrid number.
Y described in step 2.2+=yuτ/νlTo weigh the index of grid precision, wherein y is the thickness of first layer grid,
uτFor wall friction speed, νlFor the dynamic viscosity of water；Described zone of reasonableness is for underwater high-speed navigation body preferably 30≤y+≤
60。
Mesh generation software selection ICEM described in step 2.1, the grid finally determined is as shown in figure 3, be the present invention
Mesh generation schematic diagram in embodiment.
Preset value preferably 5% described in step 2.3.
Step 3：Boundary condition and computational methods are set.
Step 3 concrete methods of realizing comprises the following steps：
Step 3.1：Grid file is imported in FLUENT.
Step 3.2：Boundary condition is set.
Basin condition for import uses speed condition for import；Exit condition uses pressure export condition, and pressure considers ambient water
It is deep；Remaining is used without sliding wall boundary condition, specific as shown in figure 4, boundary condition sets schematic diagram in the embodiment of the present invention.
Step 3.3：The computation model used is set.
Cavitating flows calculating is carried out based on N-S equations, VOF interface capturing methods are used under homogeneous phase model, opens energy
Equation is measured, turbulence model uses SST k- ω models, and cavitation model uses Zwart models.
Step 3.4：Computational methods are set.
Body flow field is navigated by water to underwater high-speed using the non-fixed length solver of double precision and carries out numerical simulation, from PISO algorithms,
Pressure difference value form selects PRESTO！, the discrete method of the equation of momentum, energy equation and Equations of Turbulence inserted from Second-order Up-wind
It is worth form, CICSAM difference forms are selected in the calculating of volume fraction.
Step 4：The compressible amendment of the material property of water is realized by carrying out secondary development to CFD software.
Step 4 concrete methods of realizing comprises the following steps：
Step 4.1：Aqueous water is selected in Fluent as principal phase working medium；Usual water body is incompressible, but works as water
Lower high speed operation body headway is then considered as compressible fluid close to water body during velocity of sound, from Tait equations as compressible liquid
The state equation of body, the Tait equations for pressing water for being free from temperature influence below are：
Wherein, p∞、ρ∞And a∞For reference pressure, reference density and the velocity of sound of aqueous water at stationary stream field, k and
For attribute coefficients, the value of these variables is depending on actual condition；
Preferred p at normal temperatures and pressures∞=101325Pa, ρ∞=998kg/m3, a∞=1430m/s, k=7.15,
Containing temperature influence the Tait equations for pressing water be：
Wherein, saturation pressure and density are calculated as the function of temperature with Oldenbourg coefficients；
Wherein, θ=1-T/Tc, pc、Tc,ρcIt is critical condition (2.264 × 107Pa, 647.14K, the 332kg/m of water3)；It is public
Coefficient a in formula (3), (4)i(i=1~6), bi(i=1~6) are respectively constant；Wherein, ai=(- 7.85823,
1.83991, -11.7811,22.6705, -15.9393,1.77516)；bi=(1.99206,1.10123, -0.512506, -
1.75263, -45.4485, -6.75615E+05)；
Supercavitation process is related to vehicle repair major mixing, it is necessary to define the effective sound velocity of mixed phase；It is related to vehicle repair major mixing
The expression formula of problem effective sound velocity：
Wherein, αvFor steam volume fraction, ρmTo mix phase density, ρlAnd ρvThe density of respectively pure liquid phase and vapour phase water,
clFor the pure liquid phase velocity of sound, cvFor the velocity of sound of pure gas phase；ρl、ρv、clAnd cvValue depending on actual condition；
Preferred ρ at normal temperatures and pressuresl=998kg/m3, ρv=0.5542kg/m3, cl=1430m/s, cv=429m/s；
The state equation of compressible water and effective sound velocity are compiled into Fluent softwares by UDF, realized to material category
The amendment of property；
Step 4.2：Steam-like water is selected in Fluent as the second phase working medium；Steam-like water substance attribute is arranged to
Perfect gas.
Step 5：The numerical computations of the compressible supercavity flow dynamic characteristic of high speed are carried out in setting based on step 1 to step 4
Analysis, that is, realize a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed.
Also include step 6, step 6 concrete methods of realizing comprises the following steps：
Step 6.1：The result of calculation in Fluent is post-processed using the poster processing soft, obtains pressure, density, speed
The information of flow such as degree, volume fraction cloud atlas and quantitative distribution situation；
Step 6.2：Step 6.1 Numerical results and experimental result or Classic Experiments data are contrasted, specifically such as
It is numerical method Accuracy Verification figure in inventive embodiments, numerical result and experimental result fitting are preferable, verify number shown in Fig. 5
The accuracy of value calculating method, feasibility.
Also include step 7：Numerical result visualizes as shown in figs 6-8, different horses respectively in the embodiment of the present invention
Density cloud atlas, pressure cloud atlas and the vacuole outline drawing of conspicuous several lower Numerical Simulation Results, for these flow field results, utilize step
Rapid one to the compressible supercavity flow dynamic characteristic of a kind of high speed described in step 5 analysis method to the compressible supercavity flow of high speed
Dynamic process carries out Numerical Simulation Analysis, to disclose the compressible supercavitating flow rule of high speed and mechanism, so as to be underwater high-speed aircraft
The structure design of row body is provided fundamental basis, and can solve the problem that underwater high-speed navigates by water body practical application engineering problem.
Above-described specific descriptions, the purpose, technical scheme and beneficial effect of invention are carried out further specifically
It is bright, it should be understood that the specific embodiment that the foregoing is only the present invention, the protection model being not intended to limit the present invention
Enclose, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc., should be included in the present invention
Protection domain within.
Claims (10)
- A kind of 1. analysis method of the compressible supercavity flow dynamic characteristic of high speed, it is characterised in that：Comprise the following steps,Step 1：Establish geometrical model；Body is navigated by water as research object using the underwater high-speed of rotary structure, the geometric parameter of revolving body is determined, using modeling software Carry out Geometric Modeling；Step 2：The geometrical model established to step 1 and fluid calculation domain carry out mesh generation；The underwater high-speed obtained in step 1 navigation body Model is exported into curved surface to mesh generation software, and in mesh generation software In establish fluid calculation domain, piecemeal watershed carries out mesh generation, and carries out grid independence checking, obtains Bestgrid number；Step 3：Boundary condition and computational methods are set；Grid file is imported in FLUENT, boundary condition, the computation model used, computational methods are configured；Step 4：The compressible amendment of the material property of water is realized by carrying out secondary development to CFD software；Aqueous water and vapourous water is selected to consider that aqueous water can press in high speed as principal phase as fluid working substance, wherein aqueous water Contracting characteristic, according to the Tait state equation compiling user custom formula UDF of compressible water, the density attributes of water are repaiied Just, and effective sound velocity in water is defined；Steam-like water is arranged to the second phase, and steam-like water substance attribute is arranged to preferable gas Body；Step 5：The numerical computations point of the compressible supercavity flow dynamic characteristic of high speed are carried out in setting based on step 1 to step 4 Analysis, that is, realize a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed.
- A kind of 2. analysis method of the compressible supercavity flow dynamic characteristic of high speed as claimed in claim 1, it is characterised in that：Also wrap Step 6 is included, using a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed described in step 1 to step 5 at a high speed Compressible supercavitating flow process carries out Numerical Simulation Analysis, and numerical simulation result and experimental result are contrasted, and verifies a kind of high The accuracy of the analysis method of the compressible supercavity flow dynamic characteristic of speed, feasibility.
- A kind of 3. analysis method of the compressible supercavity flow dynamic characteristic of high speed as claimed in claim 2, it is characterised in that：Step Six concrete methods of realizing comprise the following steps,Step 6.1：The result of calculation in Fluent is post-processed using the poster processing soft, obtain pressure, density, speed, The information of flow such as volume fraction cloud atlas and quantitative distribution situation；Step 6.2：Step 6.1 Numerical results and experimental result or Classic Experiments data are contrasted, authentication number The accuracy of value calculating method, feasibility.
- A kind of 4. analysis method of the compressible supercavity flow dynamic characteristic of high speed as claimed in claim 1, it is characterised in that：Also wrap Include step 7：Using a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed described in step 1 to step 5 at a high speed Compressible supercavitating flow process carries out Numerical Simulation Analysis, to disclose the compressible supercavitating flow rule of high speed and mechanism, from And be underwater high-speed navigation body structure design provide fundamental basis, and can solve the problem that underwater high-speed navigate by water body practical application engineering Problem.
- 5. a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed as described in claim 1,2,3 or 4, its feature exist In：Step 2 concrete methods of realizing comprises the following steps,Step 2.1：The underwater high-speed obtained in step 1 navigation body Model is exported into curved surface to mesh generation software, and in grid Sufficiently large fluid calculation domain is established in division software, and sufficiently large fluid calculation domain refers to big to can ignore that Basin Boundary to water The influence in flow field around lower high speed operation body, body then is navigated by water to underwater high-speed and each face in basin is defined；Step 2.2：By basin piecemeal partition structure grid, consider that wall effect is focused on to grid around underwater high-speed navigation body Encryption, in order to ensure computational accuracy, it is necessary to assure y+ is within zone of reasonableness；Step 2.3：Grid independence checking is carried out, using resistance coefficient and vacuole profile dual indexes as validation criteria, is taken more The basin of the different grid numbers of group carries out the calculating of identical operating mode, and vacuole profile and the resistance in steady flow condition are taken for result of calculation Coefficient, grid reaches computational accuracy when front and rear prediction result twice is less than preset value with experimental result difference, now to be optimal Value, that is, obtain Bestgrid number.
- A kind of 6. analysis method of the compressible supercavity flow dynamic characteristic of high speed as claimed in claim 5, it is characterised in that：Y described in step 2.2+=yuτ/νlFor weigh grid precision index, wherein y be first layer grid thickness, uτFor Wall friction speed, νlFor the dynamic viscosity of water；Described zone of reasonableness is for underwater high-speed navigation body preferably 30≤y+≤60；Preset value described in step 2.3 selects 5%；Mesh generation software described in step 2 includes ICEM or POINTWISE softwares.
- 7. a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed as described in claim 1,2,3 or 4, its feature exist In：Step 3 concrete methods of realizing comprises the following steps,Step 3.1：Grid file is imported in FLUENT；Step 3.2：Boundary condition is set；Basin condition for import uses speed condition for import；Exit condition uses pressure export condition, and pressure considers the environment depth of water；Its It is remaining to use without sliding wall boundary condition；Step 3.3：The computation model used is set；Cavitating flows calculating is carried out based on N-S equations, VOF interface capturing methods are used under homogeneous phase model, opens energy side Journey, turbulence model use SST k- ω models, and cavitation model uses Zwart models；Step 3.4：Computational methods are set；Body flow field is navigated by water to underwater high-speed using the non-fixed length solver of double precision and carries out numerical simulation, from PISO algorithms, pressure Difference form selects PRESTO！, the discrete method of the equation of momentum, energy equation and Equations of Turbulence selects Second-order Up-wind interpolation lattice CICSAM difference forms are selected in formula, the calculating of volume fraction.
- 8. a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed as described in claim 1,2,3 or 4, its feature exist In：Step 4 concrete methods of realizing comprises the following steps,Step 4.1：Aqueous water is selected in Fluent as principal phase working medium；Usual water body is incompressible, but when underwater high Fast sail body headway is then considered as compressible fluid close to water body during velocity of sound, from Tait equations as compressible liquid State equation, the Tait equations for pressing water for being free from temperature influence below are：Wherein, p∞、ρ∞And a∞For reference pressure, reference density and the velocity of sound of aqueous water at stationary stream field, k andFor category Property coefficient, the value of these variables is depending on actual condition；Containing temperature influence the Tait equations for pressing water be：Wherein, saturation pressure and density are calculated as the function of temperature with Oldenbourg coefficients；Wherein, θ=1-T/Tc, pc、Tc,ρcIt is critical condition (2.264 × 107Pa, 647.14K, the 332kg/m of water3)；Formula (3), the coefficient a in (4)i(i=1~6), bi(i=1~6) are respectively constant；Wherein, ai=(- 7.85823,1.83991 ,- 11.7811,22.6705, -15.9393,1.77516)；bi=(1.99206,1.10123, -0.512506, -1.75263, - 45.4485-6.75615E+05)；Supercavitation process is related to vehicle repair major mixing, it is necessary to define the effective sound velocity of mixed phase；It is related to vehicle repair major mixed problem The expression formula of effective sound velocity：Wherein, αvFor steam volume fraction, ρmTo mix phase density, ρlAnd ρvThe density of respectively pure liquid phase and vapour phase water, clFor The pure liquid phase velocity of sound, cvFor the velocity of sound of pure gas phase；ρl、ρv、clAnd cvValue depending on actual condition；The state equation of compressible water and effective sound velocity are compiled into Fluent softwares by UDF, realized to material property Amendment；Step 4.2：Steam-like water is selected in Fluent as the second phase working medium；Steam-like water substance attribute is arranged to preferable Gas.
- A kind of 9. analysis method of the compressible supercavity flow dynamic characteristic of high speed as claimed in claim 8, it is characterised in that：Step P is selected in four at normal temperatures and pressures∞=101325Pa, ρ∞=998kg/m3, a∞=1430m/s, k=7.15,ρ is selected in step 4 at normal temperatures and pressuresl=998kg/m3, ρv=0.5542kg/m3, cl=1430m/s, cv=429m/s.
- 10. a kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed as described in claim 1,2,3 or 4, its feature It is：Modeling software described in step 1 selects 3D sculpting software, and described 3D sculpting software includes Solidworks or UG.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710887696.8A CN107784161A (en) | 2017-09-27 | 2017-09-27 | A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710887696.8A CN107784161A (en) | 2017-09-27 | 2017-09-27 | A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107784161A true CN107784161A (en) | 2018-03-09 |
Family
ID=61434024
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710887696.8A Pending CN107784161A (en) | 2017-09-27 | 2017-09-27 | A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107784161A (en) |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108595751A (en) * | 2018-03-14 | 2018-09-28 | 北京理工大学 | A kind of compressible cavitating flows numerical computation method |
| CN108763800A (en) * | 2018-06-04 | 2018-11-06 | 北京理工大学 | A kind of cavitation compressible flows shock-wave dynamics method for numerical simulation |
| CN109241690A (en) * | 2018-11-12 | 2019-01-18 | 扬州大学 | A kind of sluice overcurrent method of calculating flux based on CFD |
| CN109374254A (en) * | 2018-11-21 | 2019-02-22 | 北京理工大学 | A kind of analysis method of sail body water-entry cavity characteristic |
| CN109376463A (en) * | 2018-11-16 | 2019-02-22 | 重庆科技学院 | Cylinder galloping wind-structure interaction method is flexibly supported under a kind of wake flow |
| CN109740281A (en) * | 2019-01-15 | 2019-05-10 | 哈尔滨工程大学 | A method of liquid drop movement phase transition parameter is calculated based on steam flow field |
| CN109829250A (en) * | 2019-03-06 | 2019-05-31 | 大连理工大学 | The optimum design method of arm is rushed in a kind of submerged jet ditcher spray |
| CN110222467A (en) * | 2019-06-20 | 2019-09-10 | 北京理工大学 | A kind of underwater engine jet pipe swings the analysis method of high-speed jet characteristic |
| CN110472343A (en) * | 2019-08-19 | 2019-11-19 | 沈阳工程学院 | A kind of simulation method of roadbed ice-melt characteristic |
| CN110543705A (en) * | 2019-08-19 | 2019-12-06 | 西安交通大学 | Boiling simulation solving acceleration method in typical channel of nuclear reactor |
| CN110968971A (en) * | 2019-11-13 | 2020-04-07 | 中国舰船研究设计中心 | Real-scale ship sonar air guide sleeve cavitation numerical prediction method |
| CN111046580A (en) * | 2019-12-27 | 2020-04-21 | 上海理工大学 | Weakly compressible model of water for centrifugal pumps and based on the Tait equation |
| CN111159883A (en) * | 2019-12-27 | 2020-05-15 | 中国科学技术大学 | Method for analyzing micro-lens thermal reflow process in complex phase change characteristic fluid environment |
| CN111307408A (en) * | 2020-03-13 | 2020-06-19 | 海仿(上海)科技有限公司 | Cavitation flow pressure simulation method and device |
| CN111709196A (en) * | 2020-06-09 | 2020-09-25 | 哈尔滨工程大学 | Cavitation erosion resistance assessment method for underwater high-speed navigation body |
| CN111776132A (en) * | 2020-06-01 | 2020-10-16 | 中国科学院力学研究所 | Hydrofoil ship with novel supercavitation wing profiles |
| CN111814366A (en) * | 2020-05-29 | 2020-10-23 | 大连理工大学 | Forecasting method for internal wave and ventilation cavitation effect |
| CN111913390A (en) * | 2020-07-08 | 2020-11-10 | 北京电子工程总体研究所 | Steady-state bifurcation control method for supercavitation navigation body |
| CN112162957A (en) * | 2020-10-13 | 2021-01-01 | 中国空气动力研究与发展中心计算空气动力研究所 | Multi-block structure grid data compression storage method, decompression method and device |
| CN112363394A (en) * | 2020-11-02 | 2021-02-12 | 中国科学院力学研究所 | Ventilation supercavitation form control method based on machine learning |
| CN113723030A (en) * | 2021-10-18 | 2021-11-30 | 山东大学 | Actual gas physical property simulation method and system based on computational fluid dynamics |
| CN113947039A (en) * | 2021-09-24 | 2022-01-18 | 哈尔滨工程大学 | Method for forecasting movement of water outlet of navigation body and development of tail vacuole of navigation body |
| CN114840999A (en) * | 2022-04-28 | 2022-08-02 | 哈尔滨工业大学 | Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530482A (en) * | 2013-11-04 | 2014-01-22 | 东南大学 | Method for forecasting noise value of propeller in uneven inflow |
| CN104102783A (en) * | 2014-07-18 | 2014-10-15 | 安徽工业大学 | Method for forecasting numerical value of cavitation of underwater propeller tip vortex |
-
2017
- 2017-09-27 CN CN201710887696.8A patent/CN107784161A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530482A (en) * | 2013-11-04 | 2014-01-22 | 东南大学 | Method for forecasting noise value of propeller in uneven inflow |
| CN104102783A (en) * | 2014-07-18 | 2014-10-15 | 安徽工业大学 | Method for forecasting numerical value of cavitation of underwater propeller tip vortex |
Non-Patent Citations (2)
| Title |
|---|
| 张鹏: "跨超音速射弹的超空泡数值模拟", 《万方数据》 * |
| 李达钦等: "高速可压缩超空泡流动数值与理论研究", 《中国力学大会-2017暨庆祝中国力学学会成立60周年大会》 * |
Cited By (31)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108595751A (en) * | 2018-03-14 | 2018-09-28 | 北京理工大学 | A kind of compressible cavitating flows numerical computation method |
| CN108763800A (en) * | 2018-06-04 | 2018-11-06 | 北京理工大学 | A kind of cavitation compressible flows shock-wave dynamics method for numerical simulation |
| CN108763800B (en) * | 2018-06-04 | 2021-07-27 | 北京理工大学 | Cavitation compressible flow shock wave dynamics numerical simulation method |
| CN109241690A (en) * | 2018-11-12 | 2019-01-18 | 扬州大学 | A kind of sluice overcurrent method of calculating flux based on CFD |
| CN109376463A (en) * | 2018-11-16 | 2019-02-22 | 重庆科技学院 | Cylinder galloping wind-structure interaction method is flexibly supported under a kind of wake flow |
| CN109374254A (en) * | 2018-11-21 | 2019-02-22 | 北京理工大学 | A kind of analysis method of sail body water-entry cavity characteristic |
| CN109740281B (en) * | 2019-01-15 | 2022-06-21 | 哈尔滨工程大学 | Method for calculating liquid drop motion phase change parameters based on steam flow field |
| CN109740281A (en) * | 2019-01-15 | 2019-05-10 | 哈尔滨工程大学 | A method of liquid drop movement phase transition parameter is calculated based on steam flow field |
| CN109829250A (en) * | 2019-03-06 | 2019-05-31 | 大连理工大学 | The optimum design method of arm is rushed in a kind of submerged jet ditcher spray |
| CN110222467A (en) * | 2019-06-20 | 2019-09-10 | 北京理工大学 | A kind of underwater engine jet pipe swings the analysis method of high-speed jet characteristic |
| CN110543705A (en) * | 2019-08-19 | 2019-12-06 | 西安交通大学 | Boiling simulation solving acceleration method in typical channel of nuclear reactor |
| CN110472343A (en) * | 2019-08-19 | 2019-11-19 | 沈阳工程学院 | A kind of simulation method of roadbed ice-melt characteristic |
| CN110472343B (en) * | 2019-08-19 | 2023-09-08 | 辽宁石油化工大学 | Simulation calculation method for ice melting characteristics of roadbed |
| CN110968971A (en) * | 2019-11-13 | 2020-04-07 | 中国舰船研究设计中心 | Real-scale ship sonar air guide sleeve cavitation numerical prediction method |
| CN110968971B (en) * | 2019-11-13 | 2023-09-22 | 中国舰船研究设计中心 | Real-scale ship sonar dome cavitation numerical forecasting method |
| CN111159883A (en) * | 2019-12-27 | 2020-05-15 | 中国科学技术大学 | Method for analyzing micro-lens thermal reflow process in complex phase change characteristic fluid environment |
| CN111159883B (en) * | 2019-12-27 | 2023-03-10 | 中国科学技术大学 | Method for analyzing micro-lens thermal reflow process in complex phase change characteristic fluid environment |
| CN111046580A (en) * | 2019-12-27 | 2020-04-21 | 上海理工大学 | Weakly compressible model of water for centrifugal pumps and based on the Tait equation |
| CN111307408A (en) * | 2020-03-13 | 2020-06-19 | 海仿(上海)科技有限公司 | Cavitation flow pressure simulation method and device |
| CN111814366A (en) * | 2020-05-29 | 2020-10-23 | 大连理工大学 | Forecasting method for internal wave and ventilation cavitation effect |
| CN111776132A (en) * | 2020-06-01 | 2020-10-16 | 中国科学院力学研究所 | Hydrofoil ship with novel supercavitation wing profiles |
| CN111709196A (en) * | 2020-06-09 | 2020-09-25 | 哈尔滨工程大学 | Cavitation erosion resistance assessment method for underwater high-speed navigation body |
| CN111709196B (en) * | 2020-06-09 | 2022-07-15 | 哈尔滨工程大学 | Cavitation erosion resistance assessment method for underwater high-speed navigation body |
| CN111913390A (en) * | 2020-07-08 | 2020-11-10 | 北京电子工程总体研究所 | Steady-state bifurcation control method for supercavitation navigation body |
| CN112162957B (en) * | 2020-10-13 | 2022-05-27 | 中国空气动力研究与发展中心计算空气动力研究所 | Multi-block structure grid data compression storage method, decompression method and device |
| CN112162957A (en) * | 2020-10-13 | 2021-01-01 | 中国空气动力研究与发展中心计算空气动力研究所 | Multi-block structure grid data compression storage method, decompression method and device |
| CN112363394B (en) * | 2020-11-02 | 2022-04-15 | 中国科学院力学研究所 | Ventilation supercavitation form control method based on machine learning |
| CN112363394A (en) * | 2020-11-02 | 2021-02-12 | 中国科学院力学研究所 | Ventilation supercavitation form control method based on machine learning |
| CN113947039A (en) * | 2021-09-24 | 2022-01-18 | 哈尔滨工程大学 | Method for forecasting movement of water outlet of navigation body and development of tail vacuole of navigation body |
| CN113723030A (en) * | 2021-10-18 | 2021-11-30 | 山东大学 | Actual gas physical property simulation method and system based on computational fluid dynamics |
| CN114840999A (en) * | 2022-04-28 | 2022-08-02 | 哈尔滨工业大学 | Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107784161A (en) | A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed | |
| Anderson et al. | Computational fluid mechanics and heat transfer | |
| Cheng et al. | Large eddy simulation of the tip-leakage cavitating flow with an insight on how cavitation influences vorticity and turbulence | |
| Pendar et al. | Investigation of cavitation around 3D hemispherical head-form body and conical cavitators using different turbulence and cavitation models | |
| Pendar et al. | Cavitation characteristics around a sphere: An LES investigation | |
| CN110222467B (en) | Analysis method for swing high-speed jet characteristic of underwater engine spray pipe | |
| Mölder | Curved shock theory | |
| Geng et al. | Assessment of RANS turbulence models and Zwart cavitation model empirical coefficients for the simulation of unsteady cloud cavitation | |
| Long et al. | Numerical simulations of cavitating turbulent flow around a marine propeller behind the hull with analyses of the vorticity distribution and particle tracks | |
| Rashidi et al. | Numerical and experimental study of a ventilated supercavitating vehicle | |
| Katz et al. | Meshless scheme based on alignment constraints | |
| Xu et al. | Analysis of Near-wall effect on cloud cavitating flow that surrounds an axisymmetric projectile using large eddy simulation with Cartesian cut-cell mesh method | |
| CN111191347B (en) | Prediction method for underwater cavitation flow | |
| Larsen | Impact loads on circular cylinders | |
| Chen et al. | CFD-based simulation of the flow around a ship in turning motion at low speed | |
| Cao et al. | Numerical investigation of submarine hydrodynamics and flow field in steady turn | |
| Mu et al. | Numerical simulation on the cavitation flow of high speed oblique water entry of revolution body | |
| Mölder | Reflection of curved shock waves | |
| Wan et al. | Numerical simulations of viscous flows around surface ship by level set method | |
| Pendar et al. | Investigation of cavitation around disk and conical cavitators using different sets of turbulence and cavitation models | |
| Yin et al. | Numerical study of a laser generated cavitation bubble based on FVM and CLSVOF method | |
| Liu et al. | The numerical predicted of SMP11 propeller performance with and without cavitation | |
| Susan-Resiga et al. | Numerical investigation of 3D cavitating flow in Francis turbines | |
| Hejranfar et al. | Assessment of a central difference finite volume scheme for modeling of cavitating flows using preconditioned multiphase Euler equations | |
| Akbarnejad Baie et al. | Simulation of the Internal Wave of a Subsurface Vehicle in a Two-Layer Stratified Fluid |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180309 |