CN106503379B - A kind of gas burst emulation mode based on adaptive simplifying reaction and grid subdivision - Google Patents
A kind of gas burst emulation mode based on adaptive simplifying reaction and grid subdivision Download PDFInfo
- Publication number
- CN106503379B CN106503379B CN201610971312.6A CN201610971312A CN106503379B CN 106503379 B CN106503379 B CN 106503379B CN 201610971312 A CN201610971312 A CN 201610971312A CN 106503379 B CN106503379 B CN 106503379B
- Authority
- CN
- China
- Prior art keywords
- gas
- explosive
- emulation mode
- gas burst
- reaction
- 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.)
- Expired - Fee Related
Links
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 49
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 30
- 239000002360 explosive Substances 0.000 claims abstract description 37
- 230000007246 mechanism Effects 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000004880 explosion Methods 0.000 claims abstract description 12
- 230000035939 shock Effects 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 239000000126 substance Substances 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 5
- 238000009792 diffusion process Methods 0.000 claims description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012827 research and development Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000005474 detonation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007716 flux method Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000000547 structure data Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012546 transfer 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of gas burst emulation modes based on adaptive simplifying reaction and grid subdivision, which is characterized in that the gas burst emulation mode comprises determining that explosive gas forms attribute；Chemical Kinetics analysis is carried out to explosive gas, provides explosive gas detailed reaction mechanism；Using the resulting detailed reaction mechanism of adaptive simplifying reaction method, phase parameter threshold value set in explosive field caused by exploding is determined；The region big to pressure in explosive field and density gradient carries out mesh adaption subdivision, captures the temporal-spatial evolution of shock wave and shock temperature, pressure and concentration of component in explosive field.The method of the present invention can solve the safe large scale emulation of gas burst and simulate the low problem of the computational efficiency encountered in true explosion accident scene calculates.
Description
Technical field
The present invention relates to explosion safety technical fields, and in particular to a kind of based on adaptive simplifying reaction and grid subdivision
Gas burst emulation mode.
Background technique
The country is still within the starting stage in gas burst emulation mode research and development side, and research level is relatively low, does not mark
The business software of standardization, also none can the pervasive universal method for the extensive gas burst of wide scope, exploitation provides
There is the emulation mode of independent intellectual property right significant.
Summary of the invention
The purpose of the present invention is to provide a kind of gas burst emulation sides based on adaptive simplifying reaction and grid subdivision
Method is low to solve the computational efficiency that the safe large scale emulation simulation of gas burst encounters in true explosion accident scene calculates
The problem of.By the increase of adaptive simplifying elasticity of response formula and reduction component and reaction step, computational efficiency and precision are improved,
It is more accurate compared to simple reaction model, without introducing huge calculation amount.
To achieve the above object, the present invention provides a kind of imitative based on adaptive simplifying reaction and the gas burst of grid subdivision
True method, which is characterized in that the gas burst emulation mode comprises determining that explosive gas forms attribute；To fiery gas
Body carries out Chemical Kinetics analysis, obtains explosive gas detailed reaction mechanism；Using adaptive simplifying reaction method institute
The detailed reaction mechanism obtained determines phase parameter threshold value set in explosive field caused by exploding；Wherein, the phase parameter packet
Include temperature, pressure and gas component concentrations；The region big to pressure in explosive field and density gradient carries out mesh adaption subdivision,
Capture the reactive component distribution in shock wave and shock temperature in explosive field, pressure and explosion front structure.
Preferably, the phase parameter threshold value set is determined using space time correlation kinetic model.
Preferably, the step of mesh adaption segments comprises determining that the area that pressure and density gradient are big in explosive field
Domain；Dynamic grid subdivision is carried out in this region；The boundary that the region is acquired by coarse grid value interpolation, refined net is mapped to
On coarse grid；Carry out cell span.
It preferably, include: the threshold of designated phase parameter using the step of detailed reaction mechanism obtained by adaptive simplifying method
Value generates path analysis model temporally and spatially to the progress dynamic letter of explosive gas detailed reaction mechanism using more
Change.
Preferably, the gas burst emulation mode further include: using second order add semi-implicit scheme ASIRK-2B model,
Fluid dynamics conservation equation, chemical reaction Similarity Model and adaptive simplifying reaction model processing Multiple Time Scales introduce
Stiff problem.
Preferably, it is spread using five rank WENO format discrete space convective terms and six rank central difference schemes discrete spaces
, semi-implicit scheme ASIRK-2B is added in conjunction with second order to handle stiff problem caused by interruption and Multiple Time Scales.
Preferably, the gas burst emulation mode and CHEMKIN software overlap joint are used.
Preferably, it using phase parameter threshold value set in explosive field, combining adaptive Reduced mechanisms and grid, realizes
Gas burst a wide range of Flow Field Numerical Calculation on a large scale.
The method of the present invention has the advantages that
The utility model has the advantages that
1, the present invention is widely used in gas burst safe simulation research and development of software, so that gas burst emulates
It can be achieved in true explosion accident scene is tracked；
2, this technology reacts (CODAC) model, the increase of flexible type and reduction component and reaction step by adaptive simplifying,
Computational efficiency and precision are improved, it is more accurate compared to simple reaction model, without being introduced as detailed reaction model
Huge calculation amount；
3, the CODAC model that the present invention uses and the adaptive mesh divided method based on BLOCK, with flexibly and just
It in programming, is realized for gas burst analog simulation, and quickly and accurately acquisition explosion wave and thermo parameters method data mention
For effective approach.
Detailed description of the invention
Fig. 1 is a kind of specific reality of the gas burst emulation mode the present invention is based on adaptive simplifying reaction and grid subdivision
Apply the flow diagram of mode.
Fig. 2 is a kind of specific reality of the gas burst emulation mode the present invention is based on adaptive simplifying reaction and grid subdivision
Apply the flow diagram of mode.
Specific embodiment
The following examples are used to illustrate the present invention, but are not intended to limit the scope of the present invention..
As shown in Figs. 1-2, the present invention provides a kind of based on adaptive simplifying reaction and the gas burst of grid subdivision emulation
Method, which is characterized in that the gas burst emulation mode comprises determining that explosive gas forms；Explosive gas is carried out
Chemical Kinetics analysis, obtains detailed reaction mechanism；Gained detailed reaction machine is simplified using adaptive simplifying reaction model
Reason determines phase parameter threshold value set in explosive field caused by explosive gas explosion；Wherein, the phase parameter includes temperature
Degree, pressure and gas component concentrations；The region big to pressure in explosive field and density gradient carries out mesh adaption subdivision, captures
Shock wave and shock temperature, pressure and concentration of component and its temporal-spatial evolution in explosive field.
Adaptive simplifying of the present invention reacts (CODAC, Correlated Dynamic Adaptive Chemistry) model
Be it is well known to those skilled in the art, adaptive reaction mechanism is carried out using the model and is simplified.Due in explosive field, blast wave
Preceding state hardly happens variation, and the special dynamic characteristic of downstream stream remote apart from shock surface after explosion is also almost the same, therefore,
CODAC model simplification reaction mechanism has advantage richly endowed by nature.By specifying the threshold value of a relevant parameter, generate path using more
(PFA) method of analysis carries out dynamic simplification to detailed reaction mechanism over time and space, and not only the simplification with flexible type is anti-
It answers, and can guarantee precision, it can also be to avoid excessive simplification reaction template.Simultaneously using additional semi-implicit scheme in time scale
ASIRK-2B method, in conjunction with CODAC come together processing the stiff problem as caused by Multiple Time Scales.
The method of mesh adaption subdivision be it is well-known to those skilled in the art, reacted in conjunction with simplifying, establish consider it is viscous
Property diffusion and heat transfer chemical reaction fluid dynamics governing equation group, including hydrodynamic equations, state equation with
And reactive material composition equation.Conservation form high-order computational scheme is constructed, and relevant card is carried out to its convergence, stability
It is bright, it develops based on adaptive simplifying reaction mechanism gas burst calculation code, further develops INERTFACE model, will grind
The model and CHEMKIN of hair overlap, and are can be widely used in the emulation of various gas bursts so as to software.
Using self-adapting subdividing grid method, it is thin that the region big to pressure in explosive field and density gradient carries out dynamic grid
Point, accurately capturing shock and reaction coupled structure feature (temperature, pressure and concentration of component).
It can be real with real simulation explosion initiation process and blast wave propagation feature using the gas burst simulation software
When provide pressure in explosive field, Temperature Distribution provides guidance for accident rescue.The step of mesh adaption segments is wrapped
It includes: determining the region that pressure and density gradient are big in explosive field；Dynamic grid subdivision is carried out in this region；It is inserted by coarse grid value
Value acquires the boundary in the region, and refined net is mapped on coarse grid；Carry out cell span.
As shown in Figure 1, the gas burst emulation mode can also include that the gas based on dynamic simplification reaction mechanism is quick-fried
It is fried to calculate: semi-implicit scheme (ASIRK-2B), fluid dynamics conservation equation, chemical reaction Similarity Model being added using second order
And the stiff problem that adaptive simplifying reaction model processing Multiple Time Scales introduce, and use the five discrete skies of rank WENO format
Between convective term and six rank central difference schemes discrete space diffusion terms.For example, can be by the n-th step by above-mentioned model and program
Component score calculates median MnCalculate the component score of the (n+1)thth step.
As shown in Fig. 2, a kind of specific embodiment, is realized in explosive field using phase parameter threshold value set in explosive field
Adaptive reaction mechanism simplifies.
The present invention will be further described with embodiment with reference to the accompanying drawing.
1) model based on CODAC method is established, for detonation, the physical characteristic in blast wave propagation each stage passes through
Simplify response path flux method progress reaction mechanism to simplify.Space time correlation kinetics mechanism is by specifying a threshold value for phase parameter
Set is to determine, when phase parameter is fallen within the scope of changes of threshold, establishes correlation dynamics mechanism, and using it is same simplify it is anti-
Mechanism is answered, therefore, reaction simplification of the CODAC in processing gas explosive field is advantageous；
2) fluid dynamics conservation equation is established, and is coupled with reactive material component governing equation group after adaptive simplifying,
Half hidden lattice are added by five rank WENO format discrete space convective terms, six rank central difference schemes discrete space diffusion terms, second order
The stiff problem that the ASIRK-2B method processing Multiple Time Scales of formula introduce, is shown in Fig. 1.
3) region big to pressure in explosive field and density gradient determines BLOCK subdomain, and dynamic grid is carried out in subdomain
Subdivision, is acquired on the boundary of each BLOCK by coarse grid value interpolation, then further interpolation finds out the refined net inside BLOCK
Value, proceeds to, and being mapped to refined net coarse grid on synchronous with coarse grid to the value of refined net, progress cell span,
Complete grid subdivision, the final accurate capture realized to shock wave and reaction coupled structure feature.
4) combine 1), 2) and method 3), design program processes of research & development figure, rational structure data structure develop gas
Explosion simulation software.By the reliability of practical proof software, Fig. 2 is seen.
Although above having used general explanation and specific embodiment, the present invention is described in detail, at this
On the basis of invention, it can be made some modifications or improvements, this will be apparent to those skilled in the art.Therefore,
These modifications or improvements without departing from theon the basis of the spirit of the present invention are fallen within the scope of the claimed invention.
Claims (5)
1. a kind of gas burst emulation mode based on adaptive simplifying reaction and grid subdivision, which is characterized in that the gas
Explosion simulation method includes:
Determine that explosive gas forms attribute；
Chemical Kinetics analysis is carried out to explosive gas, obtains explosive gas detailed reaction mechanism；
Using the resulting detailed reaction mechanism of adaptive simplifying reaction method, phase parameter in explosive field caused by exploding is determined
Threshold value set；Wherein, the phase parameter includes temperature, pressure and gas component concentrations；
The region big to pressure in explosive field and density gradient carries out mesh adaption subdivision, captures shock wave in explosive field and swashs
Reactive component distribution in wave temperature, pressure and explosion front structure；
The step of mesh adaption segments comprises determining that the region that pressure and density gradient are big in explosive field；In the region
Middle progress dynamic grid subdivision；The boundary that the region is acquired by coarse grid value interpolation, refined net is mapped on coarse grid；Into
Row cell span；
It include: the threshold value of designated phase parameter using the step of detailed reaction mechanism obtained by adaptive simplifying method, use is mostly raw
Dynamic simplification temporally and spatially is carried out to explosive gas detailed reaction mechanism at path analysis model；
The gas burst emulation mode further include: ASIRK-2B model, the fluid dynamics of semi-implicit scheme are added using second order
The rigidity that conservation equation, chemical reaction Similarity Model and adaptive simplifying reaction model processing Multiple Time Scales introduce is asked
Topic.
2. gas burst emulation mode according to claim 1, which is characterized in that true using space time correlation kinetic model
The fixed phase parameter threshold value set.
3. gas burst emulation mode according to claim 1, which is characterized in that use five rank WENO format discrete spaces
Convective term and six rank central difference schemes discrete space diffusion terms add semi-implicit scheme ASIRK-2B in conjunction with second order to handle interruption
With stiff problem caused by Multiple Time Scales.
4. gas burst emulation mode according to claim 1, which is characterized in that by the gas burst emulation mode with
CHEMKIN software overlap joint uses.
5. gas burst emulation mode according to claim 1, which is characterized in that using phase parameter threshold value in explosive field
Set, combining adaptive Reduced mechanisms and grid realize gas burst a wide range of Flow Field Numerical Calculation on a large scale.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610971312.6A CN106503379B (en) | 2016-10-28 | 2016-10-28 | A kind of gas burst emulation mode based on adaptive simplifying reaction and grid subdivision |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610971312.6A CN106503379B (en) | 2016-10-28 | 2016-10-28 | A kind of gas burst emulation mode based on adaptive simplifying reaction and grid subdivision |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106503379A CN106503379A (en) | 2017-03-15 |
| CN106503379B true CN106503379B (en) | 2019-11-15 |
Family
ID=58323109
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610971312.6A Expired - Fee Related CN106503379B (en) | 2016-10-28 | 2016-10-28 | A kind of gas burst emulation mode based on adaptive simplifying reaction and grid subdivision |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106503379B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113553739B (en) * | 2021-07-06 | 2023-05-12 | 西安近代化学研究所 | Mixed explosive explosion output characteristic calculation method |
| CN113591345B (en) * | 2021-07-08 | 2024-01-23 | 北京理工大学 | Explosion reaction flow high-precision prediction method based on generalized Riemann solver |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102608284A (en) * | 2011-12-23 | 2012-07-25 | 南京工业大学 | Method for determining explosion limit of multi-component mixed gas |
-
2016
- 2016-10-28 CN CN201610971312.6A patent/CN106503379B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102608284A (en) * | 2011-12-23 | 2012-07-25 | 南京工业大学 | Method for determining explosion limit of multi-component mixed gas |
Non-Patent Citations (4)
| Title |
|---|
| 基于自适应化学理论甲烷预混燃烧的非主量成分的数值计算;乔瑜等;《中国电机工程学报》;20050731;第25卷(第13期);正文第1节 * |
| 基于自适应网格加密的超声速可燃气热射流起爆详细反应数值模拟;蔡晓东等;《航空动力学报》;20141031;第29卷(第10期);正文第2节 * |
| 氢氧混合气体爆轰波的真实化学反应模型数值模拟;胡湘渝等;《爆炸与冲击》;20020131;第22卷(第1期);正文第2-3节 * |
| 激波绕射碰撞加速诱导爆轰的数值模拟;孙晓晖等;《爆炸与冲击》;20110731;第31卷(第4期);第407-412页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106503379A (en) | 2017-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Farneti et al. | An assessment of Antarctic Circumpolar Current and Southern Ocean meridional overturning circulation during 1958–2007 in a suite of interannual CORE-II simulations | |
| RU2012102394A (en) | METHOD FOR CALCULATING PHYSICAL VALUES, METHOD FOR NUMERICAL ANALYSIS, PROGRAM FOR CALCULATING PHYSICAL VALUES, PROGRAM FOR NUMERICAL ANALYSIS, DEVICE FOR CALCULATING PHYSICAL VALUES AND DEVICES FOR NUMERICAL ANALYSIS | |
| Riley et al. | Interplanetary signatures of unipolar streamers and the origin of the slow solar wind | |
| CN106503379B (en) | A kind of gas burst emulation mode based on adaptive simplifying reaction and grid subdivision | |
| CN109488261B (en) | A kind of virtual construction method of ocean gas hydrate three phase dynamic | |
| Warnecke et al. | Spoke-like differential rotation in a convective dynamo with a coronal envelope | |
| CN106053016B (en) | A kind of method that random wave is made in push plate | |
| Pomoell et al. | Influence of solar wind heating formulations on the properties of shocks in the corona | |
| CN106382853A (en) | Singular perturbation suboptimal guidance law with terminal ballistic inclined angle and attack angle constraints | |
| Yuan-bin | Genetic algorithm’s application for optimization of PID parameters in dynamic positioning vessel | |
| CN114595946A (en) | Sea area theory minimum tide level calculation method, system, equipment and medium | |
| CN104369875B (en) | Spacecraft guidance control method and the system calculated based on non-linear track | |
| PETRÁŠ | Chaos in fractional-order population model | |
| CN101727514B (en) | Method for simulating passive radar receiver | |
| CN106446396A (en) | Method and device for determining influences of distributions of fractures and karst caves on reservoir permeability | |
| CN106330425B (en) | Production method based on one-dimensional displacement multiplication of chaotic maps sequence | |
| CN104156557A (en) | High-order correcting technology of boundary conditions in movement fixed wall problem | |
| Rhee et al. | Three dimensional meshless point generation technique for complex geometry | |
| Ananth et al. | Self‐consistent local composition model of electrolyte solutions | |
| CN103440406B (en) | A kind of off-line iteration particle trajectory computational methods | |
| Beyene et al. | Thermodynamics of climate change | |
| Liang et al. | Detonation structure under chain branching kinetics | |
| Zheng et al. | AUV buoyancy regulating device design and simulation analysis | |
| Mukhedkar et al. | Effects of Different Meteorological Standards on Projectile Path. | |
| CN110263462A (en) | A kind of propellant powder generation plasma rule model value emulation mode |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191115 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |