CN116306364B - Wave-absorbing simulation method for explosion water mist in cabin - Google Patents
Wave-absorbing simulation method for explosion water mist in cabin Download PDFInfo
- Publication number
- CN116306364B CN116306364B CN202310248044.5A CN202310248044A CN116306364B CN 116306364 B CN116306364 B CN 116306364B CN 202310248044 A CN202310248044 A CN 202310248044A CN 116306364 B CN116306364 B CN 116306364B
- Authority
- CN
- China
- Prior art keywords
- water mist
- mist particles
- calculating
- cabin
- particles
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 294
- 239000003595 mist Substances 0.000 title claims abstract description 279
- 238000004880 explosion Methods 0.000 title claims abstract description 83
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000004088 simulation Methods 0.000 title claims abstract description 20
- 239000002245 particle Substances 0.000 claims abstract description 204
- 230000035939 shock Effects 0.000 claims abstract description 44
- 238000001704 evaporation Methods 0.000 claims abstract description 28
- 230000008020 evaporation Effects 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 230000008878 coupling Effects 0.000 claims abstract description 21
- 238000010168 coupling process Methods 0.000 claims abstract description 21
- 238000005859 coupling reaction Methods 0.000 claims abstract description 21
- 238000003379 elimination reaction Methods 0.000 claims abstract description 18
- 230000008030 elimination Effects 0.000 claims abstract description 17
- 238000006243 chemical reaction Methods 0.000 claims abstract description 12
- 230000001133 acceleration Effects 0.000 claims description 9
- 238000009834 vaporization Methods 0.000 claims description 6
- 230000008016 vaporization Effects 0.000 claims description 6
- 230000004907 flux Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000012805 post-processing Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 2
- 239000007788 liquid Substances 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 8
- 230000001629 suppression Effects 0.000 abstract description 8
- 238000011161 development Methods 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000013461 design Methods 0.000 abstract description 3
- 238000012827 research and development Methods 0.000 abstract description 3
- 239000000243 solution Substances 0.000 description 9
- 238000002347 injection Methods 0.000 description 6
- 239000007924 injection Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000001808 coupling effect Effects 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011800 void material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/28—Fuselage, exterior or interior
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method for simulating wave elimination of explosion water mist in a cabin, and belongs to the technical field of explosion simulation in a cabin of a ship. The method for simulating the wave elimination of the explosion water mist in the cabin comprises the following steps: s1: creating a geometric model of an explosion flow field calculation domain in the cabin; s2: compiling MATLAB program to generate a water mist information file; s3: setting solving parameters and boundary conditions, and calculating biphase coupling values of explosion shock waves and water mist in the cabin to obtain flow field characteristics; s4: calculating the acting force of explosion shock waves and water mist particles, calculating deformation and crushing of the water mist particles, calculating evaporation of the water mist particles, and updating the state of the water mist particles; s5: according to the state of the water mist particles, calculating a reaction source item of the water mist particles to the flow field; s6: and updating the N-S equation to obtain a new continuous phase flow field of the explosion shock wave. The invention can reflect the attenuation condition of the explosion shock wave in the cabin, provides a basis for the attenuation analysis of the explosion shock wave in the cabin, and provides theoretical guidance for the research and development design of explosion suppression equipment in the cabin and the development of shock wave suppression technology.
Description
Technical Field
The invention relates to the technical field of explosion simulation in a cabin, in particular to a method for simulating wave elimination of explosion water mist in a cabin.
Background
In modern sea warfare, the cabin penetrating and implosion of the guided missile can cause serious damage to ships, and the water mist wave-absorbing technology can effectively attenuate the power of explosion shock waves, so that the novel protection technology is realized. In early studies, water mist wave-absorbing studies were performed by scholars in terms of experiments and numerical simulations, and the experiments were performed earlier than the numerical simulations. The three-dimensional explosion visualization is difficult to realize through the research of experiments, but the cost is high. Computer technology and computational fluid dynamics have evolved rapidly since the twenty-first century, making numerical simulation calculations another important way to study intra-cabin explosions. For the simulation of the explosion and water mist environment in the cabin, the traditional explosion software is adopted to simulate the water mist environment in the cabin, which cannot be simulated.
Disclosure of Invention
The invention aims to provide a method for simulating wave elimination of explosion water mist in a cabin, aiming at the defects of the prior art.
The invention provides a method for simulating wave elimination of explosion water mist in a cabin, which comprises the following steps:
s1: creating a geometric model of an explosion flow field calculation domain in the cabin based on Euler-Lagrangian theory, performing grid division on the geometric model to generate a grid file, and importing the grid file into Fluent software;
s2: compiling a MATLAB program to generate a water mist information file, and importing the water mist information file into Fluent software;
s3: setting solving parameters and boundary conditions in Fluent software based on a water mist information file, and calculating the biphase coupling numerical value of the explosion shock wave and the water mist in the cabin through the Fluent software for the biphase coupling N-S equation to obtain the flow field characteristic of the current time step;
s4: calculating the acting force of explosion shock waves and water mist particles based on the flow field characteristics, and calculating the resultant force and acceleration of discrete phases of the water mist particles; calculating deformation and crushing of the water mist particles according to the resultant force born by the water mist particles; calculating the evaporation of water mist particles, and updating the state of the water mist particles;
s5: according to the state of the water mist particles, calculating a reaction source item of the water mist particles to the flow field;
s6: updating the N-S equation of the biphase coupling according to the reaction source term to obtain a new continuous phase flow field of the explosion shock wave;
further, step S6 further includes:
s7: and stopping calculation when the simulation reaches the termination time, and performing post-processing on a calculation result.
Further, the step S4 specifically includes:
s41: calculating the acting force of the explosion shock wave-water mist particles based on the flow field characteristics;
s42: calculating the resultant force and acceleration of the discrete phase of the water mist particles according to the acting force of the explosion shock wave and the water mist particles, and calculating the deformation and crushing of the water mist particles;
s43: calculating the speed and the position of the water mist particles after passing a time step by adopting an Euler integration method according to the resultant force and the acceleration;
s44: and (5) calculating the evaporation of the water mist particles and updating the state of the water mist particles.
Further, in step S3, setting the solution parameters and the boundary conditions in Fluent software specifically includes: the setting of the solving parameters comprises the following steps: setting a turbulence model as an SST-k-omega turbulence model, setting a solving algorithm as Coupled, setting a discrete format as second-order windward, and setting solving precision and time step; the setting of the boundary conditions includes: the cabin boundary is set as a wall and the cabin DPM boundary is set as a reflection.
Further, the bi-directionally coupled N-S equation in step S3 includes: equation of continuityWhere α is the volume fraction of the gas phase, ρ is the gas phase density, t is the time, u is the gas phase velocity vector, S mass Is a mass source term due to evaporation of droplets; momentum equationWhere p is the pressure, T is the viscous stress tensor,mu is dynamic viscosity, I is unit matrix, F d The resistance of the water mist particles is that alpha is the volume fraction of the gas phase, rho is the gas phase density, t is the time, and u is the gas phase velocity vector; energy equationWhere E is the total energy of the gas phase, e=e+|u| 2 2, e is the gas phase internal energy; j is the diffusion heat flux; s is S energy Alpha is the volume fraction of the gas phase, which is the energy loss source term due to the evaporation of water droplets; p is the gas phase density; t is time; u is the gas phase velocity vector, p is the pressure; t is the viscous stress tensor and,the reaction source item in the step S5 and the step S6 is a quality source item S mass And energy loss source term S energy In step S6, the quality source items S are respectively passed mass Correcting the continuous equation by the energy loss source term S energy The energy equation is modified to implement the N-S equation that updates the biphase coupling.
Further, the pass-through formula in step S41Calculating the force of explosion shock wave-water mist particles, wherein m pn Is the mass of the nth water mist particle, N p The number of water mist particles contained in the gas per unit volume, u pn Is the velocity vector of the nth water mist particle, u is the velocity of the gas phaseDegree vector, C Dn Re is the resistance coefficient of the nth water mist particle pn For the nth water mist particle Reynolds number ρ p The density of water mist particles, mu is dynamic viscosity, d pn The diameter of the nth water mist particle.
Further, in step S43, the formula is passedCalculating the velocity of the water mist particles, wherein u pn Is the velocity vector of the nth water mist particle, u is the velocity vector of the gas phase, C Dn Re is the resistance coefficient of the nth water mist particle pn For the nth water mist particle Reynolds number ρ p The density of water mist particles, mu is dynamic viscosity, d pn The diameter of the nth water mist particle is the diameter, and t is the time.
Further, in step S43, the formula is passedCalculating the position of the water mist particles, wherein r pn The n-th water mist particle position, t is time, u pn Is the velocity vector of the nth water mist particle, N p Is the number of water mist particles contained in a unit volume of gas.
Further, in step S5, a quality source term S is calculated mass The method specifically comprises the following steps: according to the mass and formula of the water mist particlesCalculate the evaporation rate of the nth water mist particle +.>According to the evaporation rate and formula of the water mist particlesAcquiring quality source item S mass Wherein m is pn Is the mass of the nth water mist particle, N p The number of water mist particles contained in a unit volume of gas is given, and t is time.
Further, step S5 calculating energy loss Source term S energy The method specifically comprises the following steps: according to the mass and formula of the water mist particlesCalculate the evaporation rate of the nth water mist particle +.>According to the evaporation rate and formula of the water mist particlesCalculating the energy balance of the nth water mist particle, and according to the energy balance of the water mist particle and the formula +.>Acquiring energy loss source term S energy Wherein m is pn The mass of the nth water mist particle, t is time, C pd Specific heat of nth water mist particle, h pn The heat transfer coefficient of the nth water mist particle is T, the local gas temperature at the water mist particle is T pn Is the nth water mist particle temperature, d pn Is the diameter of the nth water mist particle, h fg Is the vaporization latent heat of water->T is the specific heat of water vapor and the vaporization latent heat of water ref Is the relative temperature (300K), ρ p Is the density of the water mist particles.
The method for simulating wave elimination of explosion water mist in the cabin has the following beneficial effects:
according to the invention, the in-cabin explosion suppression is taken as a research background, a geometric model of an in-cabin explosion flow field calculation domain is created based on Euler-Lagrange theory, a grid file is imported into Fluent software, the Fluent software is subjected to secondary development based on a self-compiling MATLAB program to obtain a water mist information file, solution parameters and boundary conditions are set in the Fluent software based on the water mist information file, numerical calculation is performed on the biphase coupling of explosion shock waves and water mist through the Fluent software, the coupling calculation of explosion shock waves and water mist bad environments is realized through the data interaction and cyclic iteration solution between the continuous phase of the explosion shock waves and the discrete phase of the water mist through the Fluent software, the attenuation condition of the explosion shock waves in the cabin can be accurately reflected by combining with the result processing, a foundation is provided for the attenuation analysis of the explosion shock waves in the cabin, and theoretical guidance is provided for the research and development design of in-cabin explosion suppression equipment and development of shock wave suppression technology.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description, serve to explain the principles of the invention. In the drawings, like reference numerals are used to identify like elements. The drawings, which are included in the description, illustrate some, but not all embodiments of the invention. Other figures can be derived from these figures by one of ordinary skill in the art without undue effort.
FIG. 1 is a schematic flow chart of a method for simulating wave elimination of explosion water mist in a cabin, which is provided by the embodiment of the invention;
FIG. 2 is a schematic diagram of a method in a method for simulating wave elimination of explosion water mist in a cabin according to an embodiment of the invention;
FIG. 3 is a flow chart of adjusting a flow field according to a source term in the method for simulating wave elimination of explosion water mist in a cabin, which is provided by the embodiment of the invention;
fig. 4 is a schematic diagram of a three-dimensional geometric model of the wave-absorbing of the explosion water mist in the cabin in the wave-absorbing simulation method of the explosion water mist in the cabin provided by the embodiment of the invention;
FIG. 5 is a three-dimensional in-cabin water mist environmental effect diagram in the in-cabin explosion water mist wave-breaking simulation method provided by the embodiment of the invention;
FIG. 6 is a diagram of the coupling effect of three-dimensional intra-cabin explosion shock wave water mist in the intra-cabin explosion water mist wave-elimination simulation method provided by the embodiment of the invention;
FIG. 7 is a pressure history chart of three-dimensional in-cabin explosion water mist wave-breaking monitoring points in the in-cabin explosion water mist wave-breaking simulation method provided by the embodiment of the invention;
fig. 8 is a graph of pressure peaks and quasi-static pressure at measuring points under different water mist quality states of three-dimensional cabin explosion water mist in the cabin explosion water mist wave-absorbing simulation method provided by the embodiment of the invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention. It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be arbitrarily combined with each other.
Please refer to fig. 1-8. The invention provides a method for simulating wave elimination of explosion water mist in a cabin, which comprises the following steps:
s1: creating a geometric model of an explosion flow field calculation domain in the cabin based on Euler-Lagrangian theory, performing grid division on the geometric model to generate a grid file, and importing the grid file into Fluent software;
s2: compiling a MATLAB program to generate a water mist information file, and importing the water mist information file into Fluent software;
s3: setting solving parameters and boundary conditions in Fluent software based on a water mist information file, and calculating the biphase coupling numerical value of the explosion shock wave and the water mist in the cabin through the Fluent software for the biphase coupling N-S equation to obtain the flow field characteristic of the current time step;
s4: calculating the acting force of explosion shock waves and water mist particles based on the flow field characteristics, and calculating the resultant force and acceleration of discrete phases of the water mist particles; calculating deformation and crushing of the water mist particles according to the resultant force born by the water mist particles; calculating the evaporation of water mist particles, and updating the state of the water mist particles;
s5: according to the state of the water mist particles, calculating a reaction source item of the water mist particles to the flow field;
s6: and updating the N-S equation of the biphase coupling according to the reaction source term to obtain a new continuous phase flow field of the explosion shock wave.
S7: and stopping calculation when the simulation reaches the termination time, and performing post-processing on a calculation result.
Here, the DPM model in Fluent software can simulate the water mist environment in the cabin. According to the equivalent conversion principle, the high-temperature high-pressure shock wave generated by explosion is converted into the initial condition of high temperature and high pressure, so that the method for absorbing the explosion water mist in the cabin can be realized in Fluent software. Aiming at the defects of the prior art, the invention solves a flow field based on Fluent software, secondarily develops Fluent functions by compiling MATLAB programs, and discloses a three-dimensional intra-cabin explosion water mist wave elimination calculation method: based on the Euler-Lagrange theory of fluid mechanics and a finite element method, numerical simulation is carried out on biphase coupling between a continuous phase of explosion shock waves and a discrete phase of water mist particles in a cabin. The continuous phase calculation of the shock wave discretely solves the Reynolds average compressible N-S equation by a finite volume method; the state of water mist discrete phases is solved through a Lagrange tracking technology, and the numerical simulation of explosion water mist wave elimination in the three-dimensional cabin is realized.
Specifically, in step S1, ANSYS space model software is used to build an intra-cabin explosion flow field geometric model, and the size is 1m×0.75m×0.7m. Adopting ANSYS Fluent Meshing software to carry out grid division on the geometric model; and generating a grid file and importing the grid file into ANSYS Fluent software.
In step S2, specifically, in order to realize the uniform water mist environment in the cabin and the wave-absorbing coupling response of explosion water mist in the cabin, a water mist information program is compiled and developed independently based on MATLAB language, and an information file for uniformly distributing water mist in the cabin is generated through the MATLAB program; parameters such as coordinates, speed, particle size, temperature, quality and the like of the single water mist particles can be defined in the program. The water mist information file contains water mist particle coordinates (x, y, z), water mist particle speeds (u, v, w), water mist particle diameters (diameter), water mist particle temperatures (temperature), and mass flow (mass-flow) of water mist particles (mass flow times time is mass). The water mist information format generated by the injection name (name) (e.g., injection:0; injection:1; injection:2 to distinguish the water mist information to be injected) is: (x y z u v w diameter temperature mass-flow) name. The water mist information file is imported in Fluent software by the following steps: 1. a Models panel; 2. a Discrete Phase sub-panel; 3. clicking Create in the Injections secondary sub-panel creates an Injection-0; 4. selecting a file mode from the Injection Type in the pop-up window; 5. the File … button is then imported directly under the window.
The step S3 includes: the setting of the solving parameters comprises the following steps: the turbulence model is set as an SST-k-omega turbulence model, the solving algorithm is set as Coupled, the discrete format is set as second order windward, and the solving accuracy and time step are set. The setting of the boundary conditions includes: the cabin boundary is set as a wall and the cabin DPM boundary is set as a reflection.
Specifically, a monitoring point P (-0.125, -0.374,0) and a monitoring section are arranged, and component transmission is started; setting a discrete item model (importing generated water mist information file, inj; medium, start/end time of setting water mist); material parameters (density of water, specific heat capacity, etc.) are set.
And starting calculation, and generating a water mist environment in the cabin. And then locally initializing, filling pressure and temperature, which are shown in the following chart, and continuously calculating to obtain the flow field characteristic of the current time step.
The step S4 specifically includes:
s41: calculating the acting force of the explosion shock wave-water mist particles based on the flow field characteristics;
s42: and calculating the resultant force and acceleration of the discrete phase of the water mist particles according to the acting force of the explosion shock wave and the water mist particles, and calculating the deformation and crushing of the water mist particles.
S43: calculating the speed and the position of the water mist particles after passing a time step by adopting an Euler integration method according to the resultant force and the acceleration;
s44: and (5) calculating the evaporation of the water mist particles and updating the state of the water mist particles.
The water mist particles can move, break and evaporate after being subjected to high-temperature and high-pressure impact, and the movement, breaking and evaporation of the water mist particles can be calculated by selecting the modes in the Fluent-DPM model.
The step S6 includes: and counting and updating the void ratio in the grid of the explosion shock wave continuous phase flow field, calculating the average speed of water mist particles in the grid, correcting an N-S equation of the biphase coupling through the calculated void ratio, the average speed of the water mist particles and a reaction source item, and then continuously solving shock wave-water mist coupling in the next calculation step.
The processing in step S7 specifically includes: leading out a case file and a data file through Fluent software; analyzing and drawing a cabin and a coupling effect diagram of water mist particles and shock waves by using post-processing software such as Tecplot and the like; and analyzing and drawing a measuring point pressure history curve chart by using post-processing software such as Origin and the like.
Specifically, the bi-directionally coupled N-S equation in step S3 includes: equation of continuityWhere α is the volume fraction of the gas phase, ρ is the gas phase density, t is the time, u is the gas phase velocity vector, S mass Is a mass source term due to evaporation of droplets; momentum equation->Wherein p is pressure, T is viscous stress tensor, +.>Mu is dynamic viscosity, I is unit matrix, F d The resistance of the water mist particles is that alpha is the volume fraction of the gas phase, rho is the gas phase density, t is the time, and u is the gas phase velocity vector; energy equationWhere E is the total energy of the gas phase, e=e+|u| 2 2, e is the gas phase internal energy; j is the diffusion heat flux; s is S energy Alpha is the volume fraction of the gas phase, which is the energy loss source term due to the evaporation of water droplets; ρ is the gas phase density; t is time; u is the gas phase velocity vector, p is the pressure; t is the viscous stress tensor; reaction in step S5 and step S6The source item is a quality source item S mass And energy loss source term S energy In step S6, the quality source items S are respectively passed mass Correcting the continuous equation by the energy loss source term S energy The energy equation is modified to implement the N-S equation that updates the biphase coupling.
Specifically, the continuous phase control equation further includes: component equation, ideal gas state equation. The component equations:wherein Y is m Is the mass fraction of the gas m, s m For mass flux of gas m, S species，m Is a mass change source term of the gas m, alpha is the volume fraction of the gas phase, rho is the gas phase density, and t is time; u is the gas phase velocity vector. The supersonic flow generated by the gas from the centre to the bulkhead is mainly the effect of the pressure difference, so the gravity effect of the gas is negligible. Ideal gas state equation: p=ρrt, where R is the universal gas constant, T is the temperature, ρ is the gas phase density.
Specifically, in step S41, the formula of the pass-through isCalculating the force of explosion shock wave-water mist particles, wherein m pn Is the mass of the nth water mist particle, N p The number of water mist particles contained in the gas per unit volume, u pn Is the velocity vector of the nth water mist particle, u is the velocity vector of the gas phase, C Dn Re is the resistance coefficient of the nth water mist particle pn For the nth water mist particle Reynolds number ρ p The density of water mist particles, mu is dynamic viscosity, d pn The diameter of the nth water mist particle.
Specifically, in step S42, the water mist particles are deformed and crushed by force, and the TAB model (model in Fluent) is used to simulate the crushing of the water mist particles. In the model, damping, forcing, springs and oscillators are used for simulating deformation r of water mist particles d Resistance, liquid viscosity μ p And surface tension. The resultant force of the water mist particles isIn the middle of Wherein C is F 、C β 、C ζ Is a constant, σ is the surface tension of water, ρ is the gas density, ρ p Density of water, mu p Is the viscosity of the water mist particles. The resistance of the high velocity air stream against the water mist particles causes the water mist particles to deform and vibrate, while the air stream is also inhibited by the various forces in the water mist particles. The surface tension of the mist particles also suppresses the deformation of the mist particles. The following solution describes the variation of the mist particle twist over time: />Wherein the Weber numberIs the ratio of aerodynamic force to surface tension; the critical value of Weber number isThe characteristic time of the distortion change is +.>The oscillation frequency of the water mist particles isDimensionless distortion of y=r d /r c (initial value of distortion of Water mist particles y 0 =0). The initial value of the distortion of the water mist particles is 0. It is assumed that when the distortion value of the water mist particles is larger than the critical value r c At this time, the water mist particles start to disintegrate (break up). Thus, when y > 1, the water mist particles start to disintegrate. Critical distortion value r c Is defined as follows: r is (r) c ＝C bk d pn Wherein C bk =0.25, other constantsObtained from experimental data; c (C) F ＝1/3；C β ＝8；C ξ =5. According to the conservation of energy (including surface potential energy and kinetic energy) of the mother water mist particles and the child water mist particles, the diameter expression of the child water mist after the nth water mist particles are crushed can be obtained:wherein d is pn Is the diameter of the nth water mist particle, y is dimensionless distortion, ρ pn Is the density of the nth water mist particle (water), σ is the surface tension of the water.
Specifically, in step S43, the formula of the pass-through isCalculating the velocity of the water mist particles, wherein u pn Is the velocity vector of the nth water mist particle, u is the velocity vector of the gas phase, C Dn Re is the resistance coefficient of the nth water mist particle pn For the nth water mist particle Reynolds number ρ p The density of water mist particles, mu is dynamic viscosity, d pn The diameter of the nth water mist particle is the diameter, and t is the time.
Specifically, in step S43, the formula of the pass-through isCalculating the position of the water mist particles, wherein r pn The n-th water mist particle position, t is time, u pn Is the velocity vector of the nth water mist particle, N p Is the number of water mist particles contained in a unit volume of gas.
Specifically, in step S44, the water mist particle mass and formula are usedThe evaporation term of the nth water mist particle is calculated.
Specifically, in step S5, a quality source term S is calculated mass The method specifically comprises the following steps: according to the mass and formula of the water mist particlesCalculate the firstEvaporation rate of n water mist particles +.>According to the evaporation rate and formula of the water mist particlesAcquiring quality source item S mass Wherein m is pn Is the mass of the nth water mist particle, N p The number of water mist particles contained in a unit volume of gas is given, and t is time.
Specifically, the energy loss source term S is calculated in step S5 energy The method specifically comprises the following steps: according to the mass and formula of the water mist particlesCalculate the evaporation rate of the nth water mist particle +.>According to the evaporation rate and formula of the water mist particlesCalculating the energy balance of the nth water mist particle, and according to the energy balance of the water mist particle and the formula +.>Acquiring energy loss source term S energy Wherein m is pn The mass of the nth water mist particle, t is time, C pd Specific heat of nth water mist particle, h pn The heat transfer coefficient of the nth water mist particle is T, the local gas temperature at the water mist particle is T pn Is the nth water mist particle temperature, d pn Is the diameter of the nth water mist particle, h fg Is the vaporization latent heat of water->T is the specific heat of water vapor and the vaporization latent heat of water ref Is the relative temperature (300K), ρ p Is the density of the water mist particles.
Mass exchange S between water mist particles and gas phase m At the same time, the gas phase density ρ and the gas phase mass fraction Y are affected, and the specific heat of water vapor is twice that of air, thereby suppressing the gas temperature. In the gas phase composition equationIn addition to the water vapor species, S of the air species，air =0 due to evaporation of water mist particles, and +.>
The invention is based on commercial software Ansys/Fluent, takes in-cabin explosion suppression as a research background, carries out secondary development on Fluent software based on self-compiling MATLAB program, carries out numerical simulation on biphase coupling between explosion shock wave and water mist, and researches the effect of water mist wave elimination. The explosive impulse field calculation is used for discretely solving a Reynolds average compressible N-S equation through a finite volume method, and the impulse continuous phase calculation is used for discretely solving the Reynolds average compressible N-S equation through the finite volume method; the state of water mist particles is solved by the water mist discrete phase through a Lagrange tracking technology. And realizing the coupling calculation of the explosion shock wave and the water mist environment and the numerical simulation of the explosion water mist wave in the three-dimensional cabin by data interaction and cyclic iteration solution between the explosion shock wave continuous phase and the water mist discrete phase in the Fluent. Compared with the prior art, the method for calculating the wave elimination of the explosion water mist in the cabin, provided by the invention, can accurately reflect the attenuation condition of the explosion shock wave in the cabin by combining the result processing, provides a basis for the attenuation analysis of the explosion shock wave in the cabin, and provides theoretical guidance for the research and development design of explosion suppression equipment in the cabin and the development of shock wave suppression technology.
The above description may be implemented alone or in various combinations and these modifications are within the scope of the present invention.
It should be noted that, in the description of the present application, the terms "upper end," "lower end," and "bottom end" of the indicated orientation or positional relationship are based on the orientation or positional relationship shown in the drawings, or the orientation or positional relationship in which the product of the application is conventionally put in use, merely for convenience of describing the present application and simplifying the description, and do not indicate or imply that the device to be referred to must have a specific orientation, be configured and operated in a specific orientation, and therefore should not be construed as limiting the present application. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and are not limiting. Although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (3)
1. The method for simulating wave elimination of explosion water mist in the cabin is characterized by comprising the following steps:
s1: creating a geometric model of an explosion flow field calculation domain in the cabin based on Euler-Lagrangian theory, performing grid division on the geometric model to generate a grid file, and importing the grid file into Fluent software;
s2: compiling a MATLAB program to generate a water mist information file, and importing the water mist information file into Fluent software;
s3: setting solving parameters and boundary conditions in Fluent software based on a water mist information file, and calculating the biphase coupling numerical value of the explosion shock wave and the water mist in the cabin through the Fluent software for the biphase coupling N-S equation to obtain the flow field characteristic of the current time step;
s4: calculating the acting force of explosion shock waves and water mist particles based on the flow field characteristics, and calculating the resultant force and acceleration of discrete phases of the water mist particles; calculating deformation and crushing of the water mist particles according to the resultant force born by the water mist particles; calculating the evaporation of water mist particles, and updating the state of the water mist particles;
s5: according to the state of the water mist particles, calculating a reaction source item of the water mist particles to the flow field;
s6: updating the N-S equation of the biphase coupling according to the reaction source term to obtain a new continuous phase flow field of the explosion shock wave; the bi-directionally coupled N-S equation includes: equation of continuityWhere α is the volume fraction of the gas phase, ρ is the gas phase density, t is the time, u is the gas phase velocity vector, S mass Is a mass source term due to evaporation of droplets; momentum equationWhere p is the pressure, T is the viscous stress tensor,mu is dynamic viscosity, I is unit matrix, F d The resistance of the water mist particles is that alpha is the volume fraction of the gas phase, rho is the gas phase density, t is the time, and u is the gas phase velocity vector; energy equationWhere E is the total energy of the gas phase, e=e+|u| 2 2, e is the gas phase internal energy, j is the diffusion heat flux, S energy Alpha is the volume fraction of the gas phase, which is the energy loss source term due to the evaporation of water droplets; ρ is the gas phase density; t is time, u is the gas phase velocity vector, p is pressure, T is the viscous stress tensor,the reaction source item in the step S5 and the step 6 is a quality source item S mass And energy loss source term S energy In step S6, the quality source items S are respectively passed mass Correcting the continuous equation by the energy loss source term S energy Correcting the energy equation to achieve an N-S equation that updates the biphase coupling; calculate quality source item S mass The method specifically comprises the following steps: according to the quality and formula of the water mist particle->Calculate the evaporation rate of the nth water mist particle +.>According to the evaporation rate and formula of the water mist particles +.>Acquiring quality source item S mass Wherein m is pn Is the mass of the nth water mist particle, N p The quantity of water mist particles contained in a unit volume of gas is represented by t, and the time is represented by t; calculating energy loss source term S energy The method specifically comprises the following steps: according to the mass and formula of the water mist particlesCalculate the evaporation rate of the nth water mist particle +.>According to the evaporation rate and formula of the water mist particlesCalculating the energy balance of the nth water mist particle, and according to the energy balance of the water mist particle and the formula +.>Acquiring energy loss source term S energy Wherein m is pn The mass of the nth water mist particle, t is time, C pd Specific heat of nth water mist particle, h pn Heat transfer system for nth water mist particleThe number T is the local gas temperature at the water mist particles, T pn Is the nth water mist particle temperature, d pn Is the diameter of the nth water mist particle, h fg Is the vaporization latent heat of water->T is the specific heat of water vapor and the vaporization latent heat of water ref Is the relative temperature (300K), ρ p Is the density of the water mist particles; the step S4 specifically includes: s41: calculating the acting force of the explosion shock wave-water mist particles based on the flow field characteristics; s42: calculating the resultant force and acceleration of the discrete phase of the water mist particles according to the acting force of the explosion shock wave and the water mist particles, and calculating the deformation and crushing of the water mist particles; s43: calculating the speed and the position of the water mist particles after passing a time step by adopting an Euler integration method according to the resultant force and the acceleration; s44: calculating the evaporation of water mist particles, and updating the state of the water mist particles; in step S41, the formula ∈ ->Calculating the force of explosion shock wave-water mist particles, wherein m pn Is the mass of the nth water mist particle, N p The number of water mist particles contained in the gas per unit volume, u pn Is the velocity vector of the nth water mist particle, u is the velocity vector of the gas phase, C Dn Re is the resistance coefficient of the nth water mist particle pn For the nth water mist particle Reynolds number ρ p The density of water mist particles, mu is dynamic viscosity, d pn Is the diameter of the nth water mist particle; step S42, stress deformation and crushing of the water mist particles, simulating the crushing of the water mist particles by using a TAB model, wherein the deformation r of the water mist particles is simulated by using damping, forcing, springs and oscillators d Resistance, liquid viscosity μ p And surface tension, the resultant force of the water mist particles isIn-> Wherein C is F 、C β 、C ξ Is a constant, σ is the surface tension of water, ρ is the gas density, ρ p Density of water, mu p Is the viscosity of water mist particles; in step S43, the formula->Calculating the velocity of the water mist particles, wherein u pn Is the velocity vector of the nth water mist particle, u is the velocity vector of the gas phase, C Dn Re is the resistance coefficient of the nth water mist particle pn For the nth water mist particle Reynolds number ρ p The density of water mist particles, mu is dynamic viscosity, d pn The diameter of the nth water mist particle is the diameter, and t is the time; in step S43, the formula->Calculating the position of the water mist particles, wherein r pn The n-th water mist particle position, t is time, u pn Is the velocity vector of the nth water mist particle, N p The number of water mist particles contained in a unit volume of gas; in step S44, according to the mass and formula of the water mist particlesThe evaporation term of the nth water mist particle is calculated.
2. The method for simulating wave elimination of explosion water mist in a cabin according to claim 1, wherein the method comprises the following steps: the step S6 further comprises the following steps:
s7: and stopping calculation when the simulation reaches the termination time, and performing post-processing on a calculation result.
3. A method for simulating wave elimination of explosion water mist in a cabin as set forth in claim 1 or 2, wherein: in step S3, setting solution parameters and boundary conditions in Fluent software specifically includes:
the setting of the solving parameters comprises the following steps: setting a turbulence model as an SST-k-omega turbulence model, setting a solving algorithm as Coupled, setting a discrete format as second-order windward, and setting solving precision and time step; the setting of the boundary conditions includes: the cabin boundary is set as a wall and the cabin DPM boundary is set as a reflection.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310248044.5A CN116306364B (en) | 2023-03-13 | 2023-03-13 | Wave-absorbing simulation method for explosion water mist in cabin |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310248044.5A CN116306364B (en) | 2023-03-13 | 2023-03-13 | Wave-absorbing simulation method for explosion water mist in cabin |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116306364A CN116306364A (en) | 2023-06-23 |
| CN116306364B true CN116306364B (en) | 2024-03-22 |
Family
ID=86788297
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310248044.5A Active CN116306364B (en) | 2023-03-13 | 2023-03-13 | Wave-absorbing simulation method for explosion water mist in cabin |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116306364B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105865283A (en) * | 2016-06-16 | 2016-08-17 | 中铁上海工程局集团有限公司 | Blasting dedusting method |
| CN109214082A (en) * | 2018-09-03 | 2019-01-15 | 北京理工大学 | A kind of high resolution numerical simulation method of near field underwater blast wave load |
| CN111474208A (en) * | 2020-04-17 | 2020-07-31 | 南京工业大学 | Simulation experiment and suppression device for gas cloud explosion of non-uniform gradient concentration formed by gas leakage in urban underground pipe gallery cabin |
| CN113177376A (en) * | 2021-05-18 | 2021-07-27 | 国网陕西省电力公司电力科学研究院 | Method and system for determining parameters of superfine water mist fire extinguishing system based on FLUENT calculation |
| CN113221473A (en) * | 2020-10-12 | 2021-08-06 | 西北工业大学 | Numerical simulation method for gas-liquid drop two-phase flow characteristics in engine combustion chamber |
| CN114861396A (en) * | 2022-03-30 | 2022-08-05 | 西北核技术研究所 | Mathematical model considering fine sand impact compression characteristic and heat absorption phase change and modeling method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7090028B2 (en) * | 2001-09-19 | 2006-08-15 | Nanomist Systems, Llc | Fire suppression using water mist with ultrafine size droplets |
| US20180239848A1 (en) * | 2017-02-21 | 2018-08-23 | Livermore Software Technology Corporation | Numerical Blast Simulation Methods and Systems Thereof |
-
2023
- 2023-03-13 CN CN202310248044.5A patent/CN116306364B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105865283A (en) * | 2016-06-16 | 2016-08-17 | 中铁上海工程局集团有限公司 | Blasting dedusting method |
| CN109214082A (en) * | 2018-09-03 | 2019-01-15 | 北京理工大学 | A kind of high resolution numerical simulation method of near field underwater blast wave load |
| CN111474208A (en) * | 2020-04-17 | 2020-07-31 | 南京工业大学 | Simulation experiment and suppression device for gas cloud explosion of non-uniform gradient concentration formed by gas leakage in urban underground pipe gallery cabin |
| CN113221473A (en) * | 2020-10-12 | 2021-08-06 | 西北工业大学 | Numerical simulation method for gas-liquid drop two-phase flow characteristics in engine combustion chamber |
| CN113177376A (en) * | 2021-05-18 | 2021-07-27 | 国网陕西省电力公司电力科学研究院 | Method and system for determining parameters of superfine water mist fire extinguishing system based on FLUENT calculation |
| CN114861396A (en) * | 2022-03-30 | 2022-08-05 | 西北核技术研究所 | Mathematical model considering fine sand impact compression characteristic and heat absorption phase change and modeling method |
Non-Patent Citations (3)
| Title |
|---|
| Fire suppression performance of water mist under diverse desmokingand ventilation conditions;Wei-Yi Liu 等;Process Safety and Environmental Protection;20201231;第230-241页 * |
| 水介质对舱内爆炸抑制作用的实验研究;李营;任广为;张玮;赵鹏铎;张磊;杜志鹏;;爆炸与冲击;20171125(06);全文 * |
| 腔体内置水、粉对瓦斯爆炸冲击波消减机制数值模拟研究;王雪;中国硕士学位论文全文数据库 工程科技Ⅰ辑;20200115(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116306364A (en) | 2023-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Verbeke et al. | The evolution of coronal mass ejections in the inner heliosphere: Implementing the spheromak model with EUHFORIA | |
| MacCormack | A numerical method for solving the equations of compressible viscous flow | |
| Guilcher et al. | Simulations of breaking wave impacts on a rigid wall at two different scales with a two-phase fluid compressible SPH model | |
| Guan et al. | Numerical study on liquid droplet internal flow under shock impact | |
| Venkatachari et al. | Numerical investigation of the interaction of counterflowing jets and supersonic capsule flows | |
| Parmar et al. | On the unsteady inviscid force on cylinders and spheres in subcritical compressible flow | |
| Haider et al. | Mathematical analysis of flow passing through a rectangular nozzle | |
| Sung et al. | On the assessment of compressibility effects of two-equation turbulence models for supersonic transition flow with flow separation | |
| Jiang et al. | Space–time conservation element and solution element method and its applications | |
| Gao et al. | On the dispersion mechanism of the flutter boundary of the AGARD 445.6 wing | |
| CN116306364B (en) | Wave-absorbing simulation method for explosion water mist in cabin | |
| Islam et al. | Plasma formation in ambient fluid from hypervelocity impacts | |
| Trinh et al. | Numerical simulation of liquid jet atomization including turbulence effects | |
| CN116204985A (en) | Method, device, equipment and storage medium for sampling wall stress tensor of aircraft | |
| West et al. | Development of modeling capabilities for launch pad acoustics and ignition transient environment prediction | |
| Zhang et al. | Hypersonic nonequilibrium flow simulations of a hemispherical nose with a counterflowing jet | |
| Hinkle et al. | Efficient Solution of Surface Erosion in Particle-Laden Hypersonic Flows | |
| Saito et al. | Shock stand-off distance of a solid sphere decelerating in transonic velocity range | |
| Hari et al. | Sloshing behavior in rigid and flexible propellant tanks: computations and experimental validation | |
| Tsutsumi et al. | Numerical analysis of ignition overpressure effect on H-IIB launch vehicle | |
| Zarchi et al. | Computational Fluid Dynamics Validation and Post-test Analysis of Supersonic Retropropulsion in the Ames 9x7 Unitary Tunnel | |
| Housman et al. | Time-derivative preconditioning methods for multicomponent flows—part ii: Two-dimensional applications | |
| Muralidharan et al. | Numerical studies of detonation initiation by supersonic projectiles using a high-order adaptive cut-cell method | |
| Kinzel et al. | An examination of thermal modeling affects to the numerical prediction of large-scale cavitating fluid flow | |
| Waters et al. | Turbulent Reactive Flow Modeling in Engines: Simulating Engine Dynamics Using fearce |
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 |