CN103136450B - Aircraft surface erosion amount balancing method under supersonic speed - Google Patents
Aircraft surface erosion amount balancing method under supersonic speed Download PDFInfo
- Publication number
- CN103136450B CN103136450B CN201310049187.XA CN201310049187A CN103136450B CN 103136450 B CN103136450 B CN 103136450B CN 201310049187 A CN201310049187 A CN 201310049187A CN 103136450 B CN103136450 B CN 103136450B
- Authority
- CN
- China
- Prior art keywords
- rho
- velocity
- alpha
- air
- particle
- 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
Links
- 230000003628 erosive effect Effects 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000002245 particle Substances 0.000 claims abstract description 51
- 238000004364 calculation method Methods 0.000 claims description 21
- 239000000463 material Substances 0.000 claims description 7
- 239000007787 solid Substances 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000004907 flux Effects 0.000 description 7
- 239000006185 dispersion Substances 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000010349 pulsation Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241000533950 Leucojum Species 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses aircraft surface erosion amount balancing method under a kind of supersonic speed, can based on the flow field situation on Euler method calculating aircraft surface and mass particle distribution situation, calculate aircraft surface erosion amount situation, the three-dimensional aircraft erosion amount that under the supersonic speed designed by the present invention, aircraft surface erosion amount balancing method calculates has good accuracy.
Description
Technical Field
The invention relates to a method for calculating the surface erosion amount of an aircraft at supersonic speed.
Background
When an aircraft flies in the atmosphere at supersonic speed, particles such as dust, raindrops, ice crystals, snowflakes and the like can collide with each other, on one hand, the particles influence the boundary layer characteristics of the collision surface to increase convective heat transfer, on the other hand, the particles directly collide at high speed to cause the quality loss of the collision surface, the phenomenon is called as erosion phenomenon, the safety of the aircraft flying at the hypersonic speed is seriously damaged by the erosion of the surface, and therefore the erosion phenomenon is an important subject of the hypersonic speed protection design of the aircraft.
According to classical erosion theory, the erosion receding rate is believed to be a function of particle velocity, particle density, angle of incidence, and the properties of the material. The erosion calculation model that can be seen at present can be written in the form:
Se=Ep/(ρtCN)
wherein,which represents the density of the particles of the incoming flow,representing the incoming particle velocity, VPThe speed of the particles at the object plane is shown,denotes the initial mass of the particle, mPRepresenting the object plane particle mass, θ is the object plane angle.
The above formula can simulate the solid particle erosion condition more effectively, but also has limitations: (1) the formula is suitable for calculating the particle orbit by the Lagrange method, the Lagrange method has advantages when the trace is in two dimensions, but in the simulation of a three-dimensional entity, the Eulerian method is a frequently adopted method, and the formula is difficult to realize in a three-dimensional flow field; (2) sin theta in the formula has better adaptability to specific two-dimensional conditions, and needs to be changed under three-dimensional conditions; (3) the formula does not take into account the impact of the impingement surface on the gas flow and there is an error in calculating the erosion.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a method for measuring the surface erosion amount of the aircraft at supersonic speed, which can calculate the surface erosion amount of the aircraft based on the flow field condition and the particle mass distribution condition of the aircraft surface calculated by the Euler method.
The invention adopts the following technical scheme to solve the technical problems: the invention designs a method for measuring the surface erosion amount of an aircraft at supersonic speed, which comprises the following specific steps:
step (1): acquiring the speed V, the density rho and the ambient pressure p of the air flowing on the surface of the aircraft in a supersonic speed state;
step (2): according to the mass and momentum conservation principle, the velocity V, the density rho and the ambient pressure p of the incoming air collected in the step (1) are used for listing the following three-dimensional Reynolds time-average equation, and the air velocity u at any position in the space is obtained by solving the equation:
wherein, subscript i, j takes on 1, 2, 3, respectively representing x, y, z direction, i.e. u1Denotes the component velocity of the air in the x-direction, u2Representing the component velocity of the air in the y-direction, u3Represents the component velocity of the air in the z direction; x is the number of1Representing a parameter on the x-axis, x2Representing a parameter, x, on the y-axis3Representing a parameter in the z-axis; μ represents the kinetic viscosity of air;
the derivative with respect to time is represented as,the derivative of x is represented as a function of,the derivative to y is represented as a function of,represents the derivative to z;
and (3) based on the air velocity calculation result in the step (2), listing mass and momentum equations for the particles and solving to obtain the volume fraction α of the particles and the velocity u of the particless:
Where ρ issDefined as the density of the solid particles, usDefined as the particle velocity, F as the gravitational acceleration, K as the inertia factor, α as the particle volume fraction at any position in space, uaDefined as the air velocity;
and (4): calculating the erosion receding rate by using the following formula according to the particle volume fraction alpha and the particle speed result obtained by calculation in the step (3):
Se=kEknV/(ρtCN)
wherein S iseDefined as the erosion receding Rate, EknDefined as the normal momentum, p, of the particle striking the surfacetDefined as the density of the impinging material, CNDefined as the erosion coefficient of the material, and k is defined as a parameter related to the profile of the eroded surface, αnDefined as the volume fraction of particles striking the surface, V, in the calculation of the trajectory of the solid particlesnDefined as the normal velocity of the particle when it strikes the surface.
Compared with the prior art, the invention has the following advantages:
the method for measuring the erosion amount of the surface of the aircraft at the supersonic speed can calculate the mass distribution condition of the particles on the surface of the aircraft by using the Euler method, is convenient to popularize to the three-dimensional condition, and has better accuracy.
Drawings
FIG. 1 is a comparison graph of erosion amount simulation calculation results.
Detailed Description
The invention is further described in detail below with reference to the accompanying drawings:
when calculating the three-dimensional erosion situation, the lagrangian method is inconvenient to calculate, and therefore, the main calculation method is the eulerian method, but a formula for calculating the erosion amount of the three-dimensional eulerian method is not seen yet, so that a new erosion amount calculation method is needed.
The invention designs a method for measuring the surface erosion amount of an aircraft at supersonic speed, which comprises the following specific steps:
step (1): acquiring the speed V, the density rho and the ambient pressure p of the air flowing on the surface of the aircraft in a supersonic speed state;
step (2): according to the mass and momentum conservation principle, the velocity V, the density rho and the ambient pressure p of the incoming air collected in the step (1) are used for listing the following three-dimensional Reynolds time-average equation, and the air velocity u at any position in the space is obtained by solving the equation:
wherein, subscript i, j takes on 1, 2, 3, respectively representing x, y, z direction, i.e. u1Denotes the component velocity of the air in the x-direction, u2Representing the component velocity of the air in the y-direction, u3Represents the component velocity of the air in the z direction; x is the number of1Representing a parameter on the x-axis, x2Representing a parameter, x, on the y-axis3Representing a parameter in the z-axis; μ represents the kinetic viscosity of air;
to representThe derivative with respect to time is that of,the derivative of x is represented as a function of,the derivative to y is represented as a function of,represents the derivative to z;
establishing a three-dimensional unsteady Navier-Stokes equation according to three conservation relations of mass, momentum and energy, and writing the equation into an integral form aiming at a certain control body omega as follows:
wherein Q is a conservative variable, FC is a non-viscous flux, Fv is a viscous flux, n is a control body in vitro vector, and S is a control body boundary:
in viscous flux, the expressions for stress and heat flux are:
in order to make the N-S equation system closed, some physical relations are also needed.
For a complete gas, there is the equation of state:
P=ρRT
let the spatial derivatives qx, qy, qz be solved by the following system of equations:
in the above equation set, the flow parameters to the respective grid nodes are used, which solve the formula as follows:
wherein: r isi=(x0,ixn)2+(y0,i-yn)2+(z0,i-zn)2]1/2
q0,,iIs the flow parameter at the grid center where node n is located, ri is the distance from point n to the grid center point.
The viscosity coefficient is composed of a laminar viscosity coefficient and a turbulent viscosity coefficient, and the determination of the vortex viscosity coefficient is illustrated by taking an RNG two-equation turbulence model as an example. Assuming that the Boussinesq vortex viscosity relationship is satisfied between the turbulent shear stress and the average strain rate tensor, the Reynolds stress can be written as:
according to dimensional analysis, μtCan be represented by k and:
wherein C isμAn empirical constant of 0.09; k is turbulent pulsation kinetic energy; is the dissipation ratio of turbulent pulsating kinetic energy. Transport method of turbulent flow pulsation kinetic energy kThe process is as follows:
after introducing a transport equation, an N-S equation is closed, a central difference format is adopted in spatial dispersion, and a multi-step Runge-Kutta format is adopted in time dispersion. The solution adopts a famous Jameson center difference method, and has two characteristics: on spatial dispersion, the flux at the inner field edge is obtained by averaging the flux values at the center of the left and right cells; in order to better capture discontinuity points and improve the stability of a format, an artificial dissipation term is added in flux calculation, the artificial dissipation term consists of second-order and fourth-order difference terms of conservation variables, the calculation of a mass, momentum and energy three-large fluid control equation is carried out on the basis of a structural grid, a constant solution is obtained by multi-step Runge-Kutta iteration in time, in addition, an accelerated convergence and residual fairing technology is adopted, the CFL number can be increased to 6-8, the convergence speed is improved by 50%, the convergence precision is controlled to be 10-8 orders of magnitude, and the far-field undisturbed boundary condition is set to be 15-fold chord length, namely a non-reflection boundary condition.
In order to obtain reliable flow field calculation parameters, the developed flow field solver is mainly based on an Euler or Navier-Stokes equation, an S-A model, A k-model and the like can be selected when A viscous flow field is solved, an explicit format-Runge-KuttA (Runge-KuttA) multi-step format and A multi-step fourth-order Runge-KuttA method can be adopted for solving the constant differential equation of time, and the format is as follows:
QI,1=QI,0+α1ΔtRI,0
QI,2=QI,0+α2ΔtRI,1
QI,3=QI,0+α3ΔtRI,2
QI,4=QI,0+α4ΔtRI,3
wherein:
for the explicit propulsion format of the Runge-Kutta method, the time step is limited for stability requirements. However, since the unsteady equation is used to solve the unsteady flow problem, and the time advance and the space dispersion are independent in the solution process, the unsteady solution is independent of the time step, so that the local time step can be adopted on each grid unit, and the flow field is advanced with the time step close to the stability limit everywhere, thereby accelerating the convergence process of the whole iterative computation.
The time step per grid cell can be calculated using the following formula:
in the above formula, Ω is the volume of the grid, u, v, w are velocity components, SX, SY, SZ are projection components of the area along the coordinate axis, c is the local sound velocity, and the value of the CFL number is determined by the specific flow condition. The explicit format has the advantages of small calculation amount and storage amount and simple program every time step is advanced. The disadvantages are that the time step length is limited by stability conditions, the number of calculation is too small, the efficiency is low, the number of steps required for convergence is large, and the total calculation time is very long, so that the characteristics of the explicit format can be described by 'small step fast running'.
And (3) based on the air velocity calculation result in the step (2), listing mass and momentum equations for the particles and solving to obtain the volume fraction α of the particles and the velocity u of the particless:
Where ρ issDefined as the density of the solid particles, usDefined as the particle velocity, F as the gravitational acceleration, K as the inertia factorα is defined as the volume fraction of particles, u, at any position in spaceaDefined as the air velocity;
and (4): calculating the erosion receding rate by using the following formula according to the particle volume fraction alpha and the particle speed result obtained by calculation in the step (3):
Se=kEknV/(ρtCN)
wherein S iseDefined as the erosion receding Rate, EknDefined as the normal momentum, p, of the particle striking the surfacetDefined as the density of the impinging material, CNDefined as the erosion coefficient of the material, and k is defined as a parameter related to the profile of the eroded surface, αnDefined as the volume fraction of particles striking the surface, V, in the calculation of the trajectory of the solid particlesnDefined as the normal velocity of the particle when it strikes the surface.
In order to verify the feasibility of the method, the erosion amount of a specific shape is calculated by the method and is compared with the test result under the same working condition and the equation simulation result in budget, as shown in fig. 1, the comparison graph of the erosion amount result calculated by the formula and the erosion amount result calculated by the original formula and the test result under the same working condition is shown in the graph, and the calculation result of the method is closer to the test value than the calculation result of the original method in the overall trend.
At the geometric leading edge point (Y is 0.00 m), the experimental value is 0.02544m, the calculated value of the original formula is 0.0299m, and the calculated value of the patent is 0.02722m, so that the method adopted by the patent is slightly smaller than the calculated value of the original formula, is closer to the experimental result, and is in good agreement.
Claims (1)
1. A method for measuring the surface erosion amount of an aircraft at supersonic speed is characterized by comprising the following specific steps:
step (1): acquiring the speed V, the density rho and the ambient pressure p of the air flowing on the surface of the aircraft in a supersonic speed state;
step (2): according to the mass and momentum conservation principle, the velocity V, the density rho and the ambient pressure p of the incoming air collected in the step (1) are used for listing the following three-dimensional Reynolds time-average equation, and the air velocity u at any position in the space is obtained by solving the equation:
wherein, subscript i, j takes on 1, 2, 3, respectively representing x, y, z direction, i.e. u1Denotes the component velocity of the air in the x-direction, u2Representing the component velocity of the air in the y-direction, u3Represents the component velocity of the air in the z direction; x is the number of1Representing a parameter on the x-axis, x2Representing a parameter, x, on the y-axis3Representing a parameter in the z-axis; μ represents the kinetic viscosity of air;
the derivative with respect to time is represented as,the derivative of x is represented as a function of,the derivative to y is represented as a function of,represents the derivative to z;
and (3) based on the air velocity calculation result in the step (2), listing mass and momentum equations for the particles and solving to obtain the volume fraction α of the particles and the velocity u of the particless:
Where ρ issDefined as the density of the solid particles, usDefining as particle velocity, F as gravitational acceleration, K as an inertia factor, α as particle volume fraction at any position in space, u as air velocity;
and (4): calculating the erosion receding rate by using the following formula according to the particle volume fraction alpha and the particle speed result obtained by calculation in the step (3):
Se=kEknV/(ρtCN)
wherein S iseDefined as the erosion receding Rate, EknDefined as the normal momentum, p, of the particle striking the surfacetDefined as the density of the impinging material, CNDefined as the erosion coefficient of the material, and k is defined as a parameter related to the profile of the eroded surface, αnDefined as the volume fraction of particles striking the surface, V, in the calculation of the trajectory of the solid particlesnDefined as the normal velocity of the particle when it strikes the surface.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310049187.XA CN103136450B (en) | 2013-02-07 | 2013-02-07 | Aircraft surface erosion amount balancing method under supersonic speed |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310049187.XA CN103136450B (en) | 2013-02-07 | 2013-02-07 | Aircraft surface erosion amount balancing method under supersonic speed |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103136450A CN103136450A (en) | 2013-06-05 |
CN103136450B true CN103136450B (en) | 2016-02-17 |
Family
ID=48496268
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310049187.XA Active CN103136450B (en) | 2013-02-07 | 2013-02-07 | Aircraft surface erosion amount balancing method under supersonic speed |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103136450B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113505551B (en) * | 2021-09-09 | 2021-12-07 | 中国空气动力研究与发展中心计算空气动力研究所 | Simulation method, system, storage medium and terminal for inducing unusual changes in incoming flow |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009008907B4 (en) * | 2009-02-13 | 2014-07-24 | Airbus Operations Gmbh | Process for plasma treatment and painting of a surface |
CN102073279A (en) * | 2010-12-30 | 2011-05-25 | 清华大学 | Composite recognition method for flight dynamics model of unmanned helicopter |
-
2013
- 2013-02-07 CN CN201310049187.XA patent/CN103136450B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN103136450A (en) | 2013-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104699947B (en) | A kind of method of use RANS/LES hybrid technologies simulated flight device rock and roll motion | |
López et al. | Verification and Validation of HiFiLES: a High-Order LES unstructured solver on multi-GPU platforms | |
Li et al. | Hydrodynamic analysis of the energy dissipation of droplets on vibrating superhydrophobic surfaces | |
Huang et al. | Study of control effects of vortex generators on a supercritical wing | |
CN103136450B (en) | Aircraft surface erosion amount balancing method under supersonic speed | |
US9390205B2 (en) | Vorticity-refinement based numerical method for simulating aircraft wing-tip vortex flows | |
Xu et al. | Delayed detached eddy simulations of fighter aircraft at high angle of attack | |
Lee et al. | Predicting aerodynamic rotor-fuselage interactions by using unstructured meshes | |
Ou et al. | Computational sports aerodynamics of a moving sphere: simulating a ping pong ball in free flight | |
Ghoreyshi et al. | CFD modeling for trajectory predictions of a generic fighter configuration | |
Denison et al. | Evaluation of CFD Predictions of Cobra-MRV Control Surface Effectiveness at the NASA Langley Unitary Plan Wind Tunnel | |
Yamashita et al. | Full-field sonic boom simulation in real atmosphere | |
Osman et al. | Numerical analysis of an external store separation from an airplane | |
Misaka et al. | Numerical simulation of jet-wake vortex interaction | |
Thomas et al. | A hybrid cfd methodology to model the two-phase flowfield beneath a hovering laboratory scale rotor | |
Aogaki et al. | Computational Study of Aerodynamic Characteristics of Reusable Rocket at High-Angle-of-Attack | |
Sathyanarayana et al. | Trajectory prediction using coupled CFD-RBD with dynamic meshing | |
Voevodenko et al. | Numerical and experimental studies of the flow on HEXAFLY-INT experimental flight test vehicle air intake | |
Sugaya et al. | Grid metrics modification approach for flow simulation around 3D geometries on Cartesian CFD method | |
Pulok et al. | Aerodynamic and vibration analysis of the morphing wings of a hypersonic vehicle | |
Palmer et al. | Comparing particle flow regimes in the L2K Arcjet with Martian entry conditions | |
Demir | Computational Fluid Dynamics Modelling of Store Separation for Transonic Generic Store | |
Zhou et al. | Simulation of water droplet trajectory in icing research of aircraft | |
Wenzhao et al. | Mesh Impact Analysis of Eulerian Method for Droplet Impingement Characteristics Under Aircraft Icing Conditions. | |
Oki et al. | Numerical simulation of transonic flow past an F-16A aircraft configuration using CASPER |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |