CN114840999B - Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model - Google Patents

Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model Download PDF

Info

Publication number
CN114840999B
CN114840999B CN202210493825.6A CN202210493825A CN114840999B CN 114840999 B CN114840999 B CN 114840999B CN 202210493825 A CN202210493825 A CN 202210493825A CN 114840999 B CN114840999 B CN 114840999B
Authority
CN
China
Prior art keywords
cavitation
revolving body
model
water
slenderness ratio
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.)
Active
Application number
CN202210493825.6A
Other languages
Chinese (zh)
Other versions
CN114840999A (en
Inventor
白瑜亮
王小刚
荣思远
张梓晨
单永志
周宏宇
王帅
韩秀杰
徐天富
张龙
郭建明
王晓光
马晓冬
张欢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202210493825.6A priority Critical patent/CN114840999B/en
Publication of CN114840999A publication Critical patent/CN114840999A/en
Application granted granted Critical
Publication of CN114840999B publication Critical patent/CN114840999B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a method for constructing a large slenderness ratio revolving body water-entering cavity evolution model, wherein a cavitator is arranged at the head part of the large slenderness ratio revolving body, and a steady cavity model of the large slenderness ratio revolving body is established; selecting different revolving body water inlet parameters as initial conditions to carry out numerical simulation to obtain cavitation evolution processes with different cavitation numbers, obtaining cavitation length and maximum cavitation diameter by measuring simulation results, and fitting the cavitation length, diameter and cavitation number to obtain the relation between the cavitation form and the cavitation number; and for the unsteady process of the large slenderness ratio revolving body entering water at high speed, correcting the unsteady cavitation model by utilizing the relation between the cavitation bubble form and the cavitation number obtained by fitting to obtain the unsteady cavitation evolution model. The invention can accurately simulate the unsteady process of the high-speed water entering of the revolving body.

Description

Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model
Technical Field
The invention relates to the technical field of underwater navigation bodies, in particular to a method for constructing a large slenderness ratio revolving body water-entering cavity evolution model.
Background
The projectile body water entering process can be regarded as a high-speed water entering process of a revolving body with a large slenderness ratio, in the process, the high-speed relative motion of the projectile body and water can reduce the pressure between the projectile body and the water to generate a low-pressure area, when the pressure is smaller than the saturated vapor pressure in the current state, the water can be vaporized to form a cavity, the revolving body is wrapped, the revolving body is separated from the water, the resistance of the revolving body is reduced, and therefore the purpose of improving the underwater navigation speed is achieved. In order to research the high-speed water entering rule of the large slenderness ratio revolving body and model the water entering trajectory, the evolution of the water entering vacuole needs to be researched.
For the unsteady process of the high-speed water entry of the revolving body with the large slenderness ratio, namely the unsteady process of the high-speed water entry of the moving body with cavitation bubbles, the prediction effect of the existing cavitation bubble model suitable for the steady condition is poor, and the unsteady process of the high-speed water entry of the revolving body cannot be accurately simulated.
Therefore, the problem to be solved by the technical personnel in the field is how to provide a method for constructing a cavitation evolution model of a large slenderness ratio revolving body entering water, which can accurately predict the cavitation form in an unsteady process.
Disclosure of Invention
In view of the above, the method corrects the air bubble shape prediction model by using a numerical simulation technology, and verifies the effectiveness of the model.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention discloses a method for constructing a large slenderness ratio revolving body water-entering cavity evolution model, wherein a cavitator is arranged at the head part of the large slenderness ratio revolving body, and the method comprises the following steps:
establishing a steady cavitation model of a large slenderness ratio revolving body;
selecting different revolving body water entry parameters as cavitation initial conditions to carry out numerical simulation to obtain cavitation evolution processes with different cavitation numbers, obtaining cavitation length and maximum cavitation diameter by measuring simulation results, and fitting the cavitation length, diameter and cavitation number to obtain the relation between cavitation form and cavitation number:
Figure BDA0003621920770000021
wherein R is n Is the radius of the cavitator, sigma is the cavitation number,
Figure BDA0003621920770000022
p 0 indicating the ambient pressure outside the vacuole, p c Which represents a saturated vapor pressure, typically 3540Pa,
Figure BDA0003621920770000023
is dynamic pressure, ρ is density of water, v is revolving body velocity, R k Is the maximum radius of the vacuole, L k Is the cavitation length;
and for the unsteady process of the large slenderness ratio revolving body entering water at high speed, correcting the steady cavitation model by utilizing the relation between the cavitation bubble form obtained by fitting and the cavitation number to obtain the unsteady cavitation evolution model.
Preferably, the expression of the steady cavitation model is as follows:
Figure BDA0003621920770000024
wherein x is the distance from the section of the cavitation bubble to the head cavitator, R c (x) Is the radius of the section of the cavitation at x, R n Radius of the cavitator, x 1 Indicates the position of the uniform cross section, x 1 =2R n ,R 1 =1.92R n
Preferably, the initial conditions include the water inlet speed and the water inlet angle of the revolving body.
Through the technical scheme, compared with the prior art, the invention has the beneficial effects that:
the method can effectively improve the prediction accuracy of the vacuole form in the unsteady process of the large slenderness ratio revolving body entering water at high speed, and the model fitting result and the numerical simulation result have high goodness of fit, thereby verifying the effectiveness of the model.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts;
FIG. 1 is a flow chart of a method for constructing a large slenderness ratio revolving body water-entering cavity evolution model provided by the embodiment of the invention;
FIG. 2 is a schematic diagram of a model of a large slenderness ratio vehicle according to an embodiment of the present invention;
FIG. 3 is a flow chart of numerical simulation of a cavitation evolution model according to an embodiment of the present invention;
FIG. 4 is a first diagram comparing the result of numerical simulation of the cavitation evolution model and the result of model fitting provided by the embodiment of the invention;
FIG. 5 is a second comparison graph of the result of numerical simulation of the cavitation evolution model and the result of model fitting provided by the embodiment of the invention;
fig. 6 is a third diagram comparing a result of numerical simulation of the cavitation evolution model and a result of model fitting provided by the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The embodiment of the invention discloses a method for constructing an evolution model of large slenderness ratio revolving body water-entering vacuoles, wherein a cavitator is arranged at the head part of the large slenderness ratio revolving body, and firstly, the evolution law of the revolving body water-entering vacuole form is explained:
the method comprises the following steps: deducing the principle of independent expansion
The independent expansion principle is an important theory of supercavity flow, and can reflect the essence and main characteristics of supercavity flow. When the revolving body enters water at a high speed to generate supercavity flow, the shape of the vacuole can be obtained by solving according to the condition that the internal pressure value of the vacuole is constant. The cavitation bubbles are formed on the cavitator and gradually develop along with the movement of the cavitator, and the independent expansion principle is that the cavitation bubbles are decomposed into a plurality of sections, the area of each section is irrelevant to the movement process of the cavitator and only relevant to the pressure difference between the cavitation bubbles and the environment and the movement property of the cavitator in the current state.
According to the energy conservation principle, the work of the revolving body on the fluid is converted into the kinetic energy and the potential energy of the fluid in the state, and therefore, a cavitation bubble section expansion equation is deduced:
Figure BDA0003621920770000031
in the formula: s represents the sectional area of the cavity at a certain moment; Δ p represents the difference between the bubble at infinity and the internal pressure at that time;
Figure BDA0003621920770000032
C x a is a semiempirical constant related to the length of cavitation bubbles and the cavitation number, and is usually between 1.5 and 2; ρ is the density of the void boundary fluid. The cavitation section expansion speed can be obtained by integrating the above formula
Figure BDA0003621920770000033
Expression (c):
Figure BDA0003621920770000041
wherein:
Figure BDA0003621920770000042
v (0) is the speed of movement of the rotor at time 0, R n Is the tail radius of the revolving body. Expansion speed of section of hollow bubble
Figure BDA0003621920770000043
And then integrating to obtain an expression of the sectional area S of the vacuole at a certain moment:
Figure BDA0003621920770000044
in the formula: s 0 Area of cavitator, S 0 =πR 0 2 ,R 0 Is the radius of the cavitator. The area of each section at the same time is drawn along the axis of the revolving bodyThe complete vacuole morphology at a time is obtained.
Step two: establishing a steady cavitation model
According to the independent expansion principle, a Soviet Union scientist gives an approximate solution and provides a Logvinovich vacuole model under a constant condition, and the expression of a vacuole form described by the vacuole model is as follows:
Figure BDA0003621920770000045
wherein: x is the distance from the section of the cavity to the head cavitator, R c (x) Is the radius of the section of the vacuole at x, R n Radius of the cavitator, x 1 Denotes the position of the uniform cross section, x =2R n ,R k Is the maximum radius of the cavitation, L k Is the cavitation length.
For the unsteady process of high-speed water entering of the large slenderness ratio revolving body, the prediction effect of the cavitation model is poor, and therefore the numerical simulation technology is utilized to fit the cavitation form to correct the cavitation model. The following explains the cavitation model correction procedure proposed in the embodiment of the present invention:
the method comprises the following steps:
firstly, a rotator model with a large slenderness ratio is established, wherein the slenderness ratio of the rotator is
Figure BDA0003621920770000046
The profile is shown in figure 2.
Selecting different water entry speeds and water entry angles of the revolving body as initial conditions, and performing numerical simulation by using CFD software to obtain cavitation evolution processes with different cavitation numbers, wherein the flow of the numerical simulation is shown in figure 3, firstly, performing mesh division on a geometrical model of the revolving body by using preprocessing software based on an overlapping mesh technology; then, a solver is set for CFD software, the CFD software comprises a basic control equation, a cavitation model and a turbulence model, and a VOF multi-phase flow model is introduced into the basic control equation; and solving and calculating recorded data, and representing the cavitation evolution process with different cavitation numbers in the form of a simulation picture after post-processing data analysis.
The proportion of cavitation bubbles to the length and the diameter of the revolving body in the simulation picture is calculated by utilizing the length and the diameter of the existing revolving body model to obtain the cavitation bubble length and the maximum cavitation bubble diameter, and the existing software fitting kit is utilized to fit the cavitation bubble length diameter and the cavitation number to obtain the relation between the cavitation bubble form and the cavitation number:
Figure BDA0003621920770000051
wherein R is n Is the radius of the cavitator, sigma is the cavitation number,
Figure BDA0003621920770000052
p 0 indicating the ambient pressure outside the vacuole, p c Which represents a saturated vapor pressure, typically 3540Pa,
Figure BDA0003621920770000053
is dynamic pressure, ρ is density of water, v is revolving body speed, R k Is the maximum radius of the cavitation, L k Is the cavitation length;
and for the unsteady process of the large slenderness ratio revolving body entering water at high speed, correcting the unsteady cavitation model by utilizing the relation between the cavitation bubble form and the cavitation number obtained by fitting to obtain the unsteady cavitation evolution model.
The simulation result is verified for the unsteady cavitation evolution model as follows:
in order to verify the applicability and the accuracy of the fitting model, the conditions of the revolving body under cavitation numbers sigma =0.027, sigma =0.045 and sigma =0.051 are selected, wherein sigma =0.027 represents a characteristic state of the revolving body under the condition of supercavitation, sigma =0.045 represents a characteristic state of the revolving body at the moment when cavitation bubbles are to be contracted to the tail, sigma =0.051 represents a characteristic state of the revolving body when the cavitation bubbles collapse to the column section of the revolving body, and the three characteristic states can describe the whole water entering process of the revolving body. The results of the numerical simulation were compared with the model fitting results, which are shown in fig. 4-6:
fig. 4a is a vacuole morphological model fitting case for σ =0.027, and fig. 4b is a numerical simulation result for σ = 0.027; fig. 5a is a fitting condition of the vacuole morphological model under the condition of σ =0.045, and fig. 5b is a result of numerical simulation under the condition of σ = 0.045; fig. 6a shows the fitting of the vacuolar morphology model for σ =0.051, and fig. 6b shows the results of numerical simulation for σ = 0.051.
The black line in fig. 4b-6b is a predicted form of cavitation, and the red line is a model form of a large aspect ratio rotor, i.e., a navigation body.
From fig. 4, it can be seen that the length of the model cavitation bubble is predicted to be 3.97m and the maximum radius of the cavitation bubble is 0.0913m by fitting under the condition that the cavitation number σ =0.027, and the degree of coincidence is large by comparing with the numerical simulation.
From fig. 5, it can be obtained that the cavitation length of the fitted prediction model is 2.65m and the maximum radius of the cavitation is 0.0827m under the condition that the cavitation number σ =0.045, which is basically consistent with the actual condition.
From fig. 6, it can be seen that when the σ =0.051 condition is met, the cavitation bubbles are collapsed under the revolving body model, the position of the closed cavitation bubbles on the surface of the revolving body is 1.77m away from the head, the maximum cavitation radius is 0.0801m, the prediction of the cavitation bubble model on the wet part of the revolving body can be found to be accurate after comparison with a numerical simulation diagram, and the accuracy of the fitting cavitation bubble prediction model is verified through result analysis of the three conditions.
The method for constructing the evolution model of the large slenderness ratio revolving body cavity in water is described in detail, a specific example is applied in the embodiment to explain the principle and the implementation mode of the invention, and the description of the embodiment is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined in this embodiment may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (3)

1. A method for constructing a large slenderness ratio revolving body water-entering cavity evolution model is characterized by comprising the following steps of:
establishing a steady cavitation model of a large slenderness ratio revolving body;
selecting different revolving body water entry parameters as initial conditions to carry out numerical simulation to obtain cavitation evolution processes with different cavitation numbers, obtaining cavitation length and maximum cavitation diameter by measuring simulation results, and fitting the cavitation length, diameter and cavitation number to obtain the relation between the cavitation form and the cavitation number:
Figure FDA0004066122920000011
wherein R is n Is the radius of the cavitator, sigma is the cavitation number,
Figure FDA0004066122920000012
p 0 indicating the ambient pressure outside the vacuole, p c Represents a saturated vapor pressure of 3540Pa,
Figure FDA0004066122920000013
is dynamic pressure, ρ is density of water, v is revolving body speed, R k Is the maximum radius of the cavitation, L k Is the cavitation length;
and for the unsteady process of the large slenderness ratio revolving body entering water at high speed, correcting the steady cavitation model by utilizing the relation between the cavitation bubble form obtained by fitting and the cavitation number to obtain the unsteady cavitation evolution model.
2. The method for constructing the cavitation evolution model of the large slenderness ratio revolving body entering water according to claim 1, wherein the expression of the steady cavitation model is as follows:
Figure FDA0004066122920000014
wherein x is the distance from the section of the cavitation bubble to the head cavitator, R c (x) Is the radius of the section of the cavitation at x, R n Radius of the cavitator, x 1 Indicates the position of the uniform cross section, x 1 =2R n ,R 1 =1.92R n
3. The method for constructing the evolution model of the cavitation bubbles of the large slenderness ratio revolving body in the water as claimed in claim 1, wherein the initial conditions comprise the revolving body in-water speed and the in-water angle.
CN202210493825.6A 2022-04-28 2022-04-28 Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model Active CN114840999B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210493825.6A CN114840999B (en) 2022-04-28 2022-04-28 Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210493825.6A CN114840999B (en) 2022-04-28 2022-04-28 Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model

Publications (2)

Publication Number Publication Date
CN114840999A CN114840999A (en) 2022-08-02
CN114840999B true CN114840999B (en) 2023-03-14

Family

ID=82567355

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210493825.6A Active CN114840999B (en) 2022-04-28 2022-04-28 Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model

Country Status (1)

Country Link
CN (1) CN114840999B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104280207A (en) * 2014-10-20 2015-01-14 中国运载火箭技术研究院 Novel vacuole form recognition and hydrodynamic measuring method
CN108763800A (en) * 2018-06-04 2018-11-06 北京理工大学 A kind of cavitation compressible flows shock-wave dynamics method for numerical simulation
CN111814366A (en) * 2020-05-29 2020-10-23 大连理工大学 Forecasting method for internal wave and ventilation cavitation effect
CN113139306A (en) * 2020-01-17 2021-07-20 哈尔滨工业大学 Numerical prediction method for cavitation noise of composite propeller
CN113885539A (en) * 2021-09-29 2022-01-04 哈尔滨工业大学 Design method of super-cavity navigation body LPV controller

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107784161A (en) * 2017-09-27 2018-03-09 北京理工大学 A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed
CN113821996B (en) * 2021-07-14 2024-01-30 南京理工大学 Novel rapid calculation method for high-speed entry trajectory of projectile
CN114332035B (en) * 2021-12-30 2024-06-25 中国科学院力学研究所 Method for measuring medium parameters in cavitation bubbles

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104280207A (en) * 2014-10-20 2015-01-14 中国运载火箭技术研究院 Novel vacuole form recognition and hydrodynamic measuring method
CN108763800A (en) * 2018-06-04 2018-11-06 北京理工大学 A kind of cavitation compressible flows shock-wave dynamics method for numerical simulation
CN113139306A (en) * 2020-01-17 2021-07-20 哈尔滨工业大学 Numerical prediction method for cavitation noise of composite propeller
CN111814366A (en) * 2020-05-29 2020-10-23 大连理工大学 Forecasting method for internal wave and ventilation cavitation effect
CN113885539A (en) * 2021-09-29 2022-01-04 哈尔滨工业大学 Design method of super-cavity navigation body LPV controller

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
不同密度回转体高速倾斜入水空化数值模拟;穆青等;《兵工学报》;20200415;118-123 *
曹伟 等.自然超空泡形态特性的射弹试验研究.《工程力学》.2006,第23卷(第12期),175-179,187. *
水下高速射弹空泡形态与阻力特性研究;弋辉等;《鱼雷技术》;20110815(第04期);3-7 *
水下高速射弹超空泡形态特性的数值模拟研究;易文俊等;《舰船科学技术》;20080815(第04期);118-122 *
水下高速航行体自然超空泡形态特性仿真研究;易文俊等;《南京理工大学学报(自然科学版)》;20090630(第03期);58-62 *
自然超空泡形态特性的射弹试验研究;曹伟 等;《工程力学》;20061231;第23卷(第12期);175-179,187 *
超声速细长锥型射弹超空泡流动数值计算方法;张志宏等;《华中科技大学学报(自然科学版)》;20140121(第01期);44-48+114 *
轴对称弹体非定常空泡流数值仿真;胡晓等;《计算机仿真》;20150415(第04期);40-44 *

Also Published As

Publication number Publication date
CN114840999A (en) 2022-08-02

Similar Documents

Publication Publication Date Title
CN109948301B (en) Near-water surface sliding jump fluid-solid coupling numerical value prediction method based on grid control
CN107784161A (en) A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed
CN110222467B (en) Analysis method for swing high-speed jet characteristic of underwater engine spray pipe
CN110321588B (en) Rail vehicle air resistance calculation method based on numerical simulation
CN108319793A (en) A kind of low speed based on CFD tilts into the water the analysis method of free movement
CN114168796B (en) Method for establishing high-altitude aerodynamic database of aircraft
CN116341421B (en) Hypersonic flow field numerical simulation method, hypersonic flow field numerical simulation system, electronic equipment and storage medium
CN105760602A (en) Total flow field numerical simulation method for finite volume weighted essentially non-oscillatory scheme
CN113221478A (en) Numerical analysis method and system for actively controlling rotor wing
CN115310388B (en) Method for calculating three-dimensional asymmetric double-Gaussian wake flow wind speed of wind turbine with space variation
CN111709196A (en) Cavitation erosion resistance assessment method for underwater high-speed navigation body
CN106886666B (en) Boundary layer grid solving method for lift-drag ratio of underwater glider under different attack angles
CN113128096A (en) Method for acquiring direct navigation additional mass of underwater vehicle
CN114840999B (en) Method for constructing large slenderness ratio revolving body water-entering vacuole evolution model
Ma et al. Effects of uncertainties in the launch parameters on the pressure-equalizing air film around a vertically launched underwater vehicle
CN109325288A (en) A kind of Solid Launch Vehicle population parameter based on uncertainty optimization determines method and system
CN110210093B (en) Verification method and device for reliability of impact load engineering model
Cummings et al. Supersonic, turbulent flow computation and drag optimization for axisymmetric afterbodies
CN112613247A (en) Flow field analysis method and system for filling high-porosity porous medium open channel
Soltani et al. Numerical investigation of the unstart suppression in a supersonic air intake
CN116415518A (en) WENO-based high-resolution high-precision shock wave intermittent capturing method
CN110750836B (en) Aircraft pitching damping moment calculation method based on steady-state cone movement
CN111144036B (en) Multi-scale simulation analysis method for propeller noise
CN112182771B (en) Data processing method based on numerical simulation, storage medium and electronic device
CN112446126B (en) Simulation method for tail beat motion state of supercavitation navigation body

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
GR01 Patent grant
GR01 Patent grant