CN106503379A - A kind of gas burst emulation mode that is reacted based on adaptive simplifying with grid subdivision - Google Patents
A kind of gas burst emulation mode that is reacted based on adaptive simplifying with grid subdivision Download PDFInfo
- Publication number
- CN106503379A CN106503379A CN201610971312.6A CN201610971312A CN106503379A CN 106503379 A CN106503379 A CN 106503379A CN 201610971312 A CN201610971312 A CN 201610971312A CN 106503379 A CN106503379 A CN 106503379A
- Authority
- CN
- China
- Prior art keywords
- gas
- explosive
- emulation mode
- gas burst
- adaptive
- 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.)
- Granted
Links
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 28
- 238000006243 chemical reaction Methods 0.000 claims abstract description 43
- 239000002360 explosive Substances 0.000 claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000007246 mechanism Effects 0.000 claims abstract description 23
- 230000035939 shock Effects 0.000 claims abstract description 10
- 238000004880 explosion Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims abstract description 5
- 239000000203 mixture Substances 0.000 claims abstract description 4
- 239000000126 substance Substances 0.000 claims abstract description 4
- 230000008569 process Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000009792 diffusion process Methods 0.000 claims description 4
- 230000006978 adaptation Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000005516 engineering process 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
- 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
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 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
- 238000007716 flux method Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000012545 processing Methods 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 mode that is reacted based on adaptive simplifying with grid subdivision, it is characterised in that the gas burst emulation mode includes:Determine explosive gas composition attribute;Chemical Kinetics analysis is carried out to explosive gas, explosive gas detailed reaction mechanism is given;Using the detailed reaction mechanism of adaptive simplifying reaction method gained, phase parameter threshold value set in the explosive field produced by blast is determined;The region big to pressure in explosive field and density gradient carries out mesh adaption subdivision, catches the temporal-spatial evolution of shock wave and shock temperature, pressure and concentration of component in explosive field.The inventive method can solve the problems, such as that the computational efficiency that the safe large scale emulation simulation of gas burst is run in true explosion accident scene calculates is low.
Description
Technical field
The present invention relates to explosion safety technical field, and in particular to a kind of based on adaptive simplifying reaction and grid subdivision
Gas burst emulation mode.
Background technology
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 neither one 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.
Content of the invention
It is an object of the invention to provide a kind of gas burst emulation side that is reacted based on adaptive simplifying with grid subdivision
Method is low in order to solve the computational efficiency that the simulation of gas burst safe large scale emulation run in true explosion accident scene calculates
Problem.By increase and reduction component and the reaction step of adaptive simplifying elasticity of response formula, computational efficiency and precision is improve,
Compare that simple reaction model is more accurate, and huge amount of calculation will not be introduced.
For achieving 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, it is characterised in that the gas burst emulation mode includes:Determine explosive gas composition 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 for obtaining, determines phase parameter threshold value set in the explosive field produced by blast;Wherein, the phase parameter bag
Include temperature, pressure and gas component concentrations;The region big to pressure in explosive field and density gradient carries out mesh adaption subdivision,
Catch shock wave and shock temperature, pressure and the reactive component distribution that explodes in front structure in explosive field.
Preferably, the phase parameter threshold value set is determined using space time correlation kinetic model.
Preferably, the step of mesh adaption is segmented includes:Determine the big area of pressure and density gradient in explosive field
Domain;Dynamic grid subdivision is carried out in this region;The border that the region is tried to achieve by coarse grid value interpolation, refined net is mapped to
On coarse grid;Carry out cell span.
Preferably, using adaptive simplifying method gained detailed reaction mechanism the step of, includes:The threshold of designated phase parameter
Value, temporally and spatially enters Mobile state letter using the path analysis models that generate more to explosive gas detailed reaction mechanism
Change.
Preferably, the gas burst emulation mode also includes:Using the ASIRK-2B models of the additional semi-implicit scheme of second order,
Hydrodynamics conservation equation, chemical reaction Similarity Model and adaptive simplifying reaction model process Multiple Time Scales and introduce
Stiff problem.
Preferably, spread using five rank WENO form discrete space convective terms and six rank central difference schemes discrete spaces
, the stiff problem that interruption and Multiple Time Scales cause is processed in conjunction with the additional semi-implicit scheme ASIRK-2B of second order.
Preferably, the gas burst emulation mode is used with CHEMKIN softwares overlap joint.
Preferably, using phase parameter threshold value set in explosive field, combining adaptive Reduced mechanisms and grid, realize
The extensive Flow Field Numerical Calculation on a large scale of gas burst.
The inventive method has the advantage that:
Beneficial effect:
1st, the present invention is widely used in gas burst safe simulation research and development of software, and then causes gas burst emulation
Can be achieved in real explosion accident scene is followed the trail of;
2nd, this technology reacts (CODAC) model, the increase of flexible type and reduction component and reaction step by adaptive simplifying,
Computational efficiency and precision is improve, it is more accurate to compare simple reaction model, and will not introduce as detailed reaction model
Huge amount of calculation;
3rd, the present invention is adopted CODAC models and the adaptive mesh divided method based on BLOCK, with flexibly and just
In programming, it is that gas burst analog simulation is realized, and quickly and accurately acquisition explosion wave and thermo parameters method data is carried
For effective approach.
Description of the drawings
Fig. 1 is that the present invention reacts a kind of concrete reality with the gas burst emulation mode of grid subdivision based on adaptive simplifying
Apply the schematic flow sheet of mode.
Fig. 2 is that the present invention reacts a kind of concrete reality with the gas burst emulation mode of grid subdivision based on adaptive simplifying
Apply the schematic flow sheet of mode.
Specific embodiment
Following examples are used for the present invention to be described, but are not limited to the scope of the present invention.
As shown in Figure 1-2, the present invention provides a kind of gas burst emulation that reacts with grid subdivision based on adaptive simplifying
Method, it is characterised in that the gas burst emulation mode includes:Determine that explosive gas is constituted;Explosive gas is carried out
Chemical Kinetics is analyzed, and obtains detailed reaction mechanism;Gained detailed reaction machine is simplified using adaptive simplifying reaction model
Reason, determines phase parameter threshold value set in the explosive field produced by explosive gas blast;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, catches
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
It is well known to those skilled in the art, self adaptation reaction mechanism simplification is carried out using the model.Due in explosive field, blast wave
Front state hardly changes, and also almost identical apart from the special dynamic characteristic of the remote downstream stream of shock surface after exploding, therefore,
CODAC model simplification reaction mechanisms have advantage richly endowed by nature.By specifying the threshold value of a relevant parameter, generate path using more
Analysis (PFA) method, carries out dynamic simplification to detailed reaction mechanism over time and space, not only has the simplification of flexible type anti-
Should, and can guarantee that precision, it is also possible to avoid excessive simplified reaction template.Additional semi-implicit scheme is adopted in time scale simultaneously
ASIRK-2B methods, in conjunction with the stiff problem that the CODAC places of coming together reason Multiple Time Scales cause.
The method of mesh adaption subdivision is well-known to those skilled in the art, in conjunction with reaction is simplified, sets up consideration viscous
Property diffusion and heat transfer chemical reaction hydrodynamics governing equation group, including hydrodynamic equations, state equation with
And reactive material composition equation.Construction conservation form high-order computational scheme, and the card of correlation is carried out to its convergence, stability
Bright, develop based on adaptive simplifying reaction mechanism gas burst calculation code, further develop INERTFACE models, will grind
The model that sends out is overlapped with CHEMKIN, so that software is can be widely used in the emulation of various gas bursts.
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).
Using the gas burst simulation software can true simulated explosion initiation process, and BP feature is real
When provide pressure in explosive field, Temperature Distribution provides guidance for accident rescue.The step of mesh adaption is segmented is wrapped
Include:Determine the big region of pressure and density gradient in explosive field;Dynamic grid subdivision is carried out in this region;Inserted by coarse grid value
Value tries to achieve the border 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
Fried calculating:Second order additional semi-implicit scheme (ASIRK-2B), hydrodynamics conservation equation, chemical reaction Similarity Model are adopted
And adaptive simplifying reaction model processes the stiff problem that Multiple Time Scales are introduced, and adopt the discrete sky of five rank WENO forms
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
Self adaptation reaction mechanism simplifies.
The present invention will be further described with embodiment below in conjunction with the accompanying drawings.
1) model based on CODAC methods is set up, for detonation, the physical characteristic in each stage of BP, is passed through
Simplifying response path flux method carries out reaction mechanism simplification.Space time correlation kinetics mechanism is by specifying a threshold value for phase parameter
Gather to determine, when phase parameter is fallen in the range of changes of threshold, set up correlation dynamics mechanism, and anti-using same simplification
Mechanism is answered, therefore, reactions of the CODAC in processing gas explosive field simplifies advantageous;
2) hydrodynamics conservation equation is set up, and is coupled with reactive material component governing equation group after adaptive simplifying,
By five rank WENO form discrete space convective terms, six rank central difference schemes discrete space diffusion terms, the additional half hidden lattice of second order
The ASIRK-2B methods of formula process the stiff problem that Multiple Time Scales are introduced, and see Fig. 1.
3) region big to pressure in explosive field and density gradient determines BLOCK subdomains, carries out dynamic grid in subdomain
Subdivision, is tried to achieve by coarse grid value interpolation on the border of each BLOCK, and then further interpolation obtains the refined net inside BLOCK
Value, proceeds to, and being mapped to refined net coarse grid on synchronous with coarse grid, carries out cell span to the value of refined net,
Grid subdivision is completed, the final accurate seizure that realizes to shock wave and reaction coupled structure feature.
4) combine 1), 2) with method 3), design program processes of research & development figure, rational structure data structure develops gas
Explosion simulation software.By the reliability of practical proof software, Fig. 2 is seen.
Although, the present invention is described in detail above to have used general explanation and specific embodiment, 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, belong to the scope of protection of present invention without departing from theon the basis of the spirit of the present invention.
Claims (8)
1. a kind of based on adaptive simplifying reaction and grid subdivision gas burst emulation mode, it is characterised in that the gas
Explosion simulation method includes:
Determine explosive gas composition attribute;
Chemical Kinetics analysis is carried out to explosive gas, explosive gas detailed reaction mechanism is obtained;
Using the detailed reaction mechanism of adaptive simplifying reaction method gained, phase parameter in the explosive field produced by blast 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, catches in explosive field shock wave and swashs
Reactive component distribution in ripple temperature, pressure and blast front structure.
2. gas burst emulation mode according to claim 1, it is characterised 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, it is characterised in that the step of mesh adaption is segmented
Including:Determine the big region of pressure and density gradient in explosive field;Dynamic grid subdivision is carried out in this region;By coarse grid value
Interpolation tries to achieve the border in the region, and refined net is mapped on coarse grid;Carry out cell span.
4. gas burst emulation mode according to claim 1, it is characterised in that detailed using adaptive simplifying method gained
The step of thin reaction mechanism, includes:The threshold value of designated phase parameter, generates path analysis model in time and space using more
On dynamic simplification is carried out to explosive gas detailed reaction mechanism.
5. gas burst emulation mode according to claim 1, it is characterised in that the gas burst emulation mode is also wrapped
Include:Adopt the ASIRK-2B models of the additional semi-implicit scheme of second order, hydrodynamics conservation equation, chemical reaction Similarity Model with
And adaptive simplifying reaction model processes the stiff problem that Multiple Time Scales are introduced.
6. gas burst emulation mode according to claim 5, it is characterised in that using five rank WENO form discrete spaces
Convective term and six rank central difference schemes discrete space diffusion terms, process interruption in conjunction with the additional semi-implicit scheme ASIRK-2B of second order
The stiff problem caused with Multiple Time Scales.
7. gas burst emulation mode according to claim 1, it is characterised in that by the gas burst emulation mode with
CHEMKIN softwares overlap joint is used.
8. gas burst emulation mode according to claim 1, it is characterised in that using phase parameter threshold value in explosive field
Set, combining adaptive Reduced mechanisms and grid, realize the extensive Flow Field Numerical Calculation on a large scale of gas burst.
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 true CN106503379A (en) | 2017-03-15 |
CN106503379B 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) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113553739A (en) * | 2021-07-06 | 2021-10-26 | 西安近代化学研究所 | Method for calculating explosion output characteristics of mixed explosive |
CN113591345A (en) * | 2021-07-08 | 2021-11-02 | 北京理工大学 | Explosion reaction flow high-precision prediction method based on generalized Riemann solution method device |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102608284A (en) * | 2011-12-23 | 2012-07-25 | 南京工业大学 | Method for determining explosion limit of multicomponent 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 multicomponent mixed gas |
Non-Patent Citations (4)
Title |
---|
乔瑜等: "基于自适应化学理论甲烷预混燃烧的非主量成分的数值计算", 《中国电机工程学报》 * |
孙晓晖等: "激波绕射碰撞加速诱导爆轰的数值模拟", 《爆炸与冲击》 * |
胡湘渝等: "氢氧混合气体爆轰波的真实化学反应模型数值模拟", 《爆炸与冲击》 * |
蔡晓东等: "基于自适应网格加密的超声速可燃气热射流起爆详细反应数值模拟", 《航空动力学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113553739A (en) * | 2021-07-06 | 2021-10-26 | 西安近代化学研究所 | Method for calculating explosion output characteristics of mixed explosive |
CN113553739B (en) * | 2021-07-06 | 2023-05-12 | 西安近代化学研究所 | Mixed explosive explosion output characteristic calculation method |
CN113591345A (en) * | 2021-07-08 | 2021-11-02 | 北京理工大学 | Explosion reaction flow high-precision prediction method based on generalized Riemann solution method device |
CN113591345B (en) * | 2021-07-08 | 2024-01-23 | 北京理工大学 | Explosion reaction flow high-precision prediction method based on generalized Riemann solver |
Also Published As
Publication number | Publication date |
---|---|
CN106503379B (en) | 2019-11-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 | |
CN109488261B (en) | A kind of virtual construction method of ocean gas hydrate three phase dynamic | |
Park et al. | Evaluation of global atmospheric solvers using extensions of the Jablonowski and Williamson baroclinic wave test case | |
CN110750933B (en) | Accurate interface tracking processing method for coupling Lagrange particles and Euler method | |
De Andrés et al. | Time domain model for a two-body heave converter: Model and applications | |
CN106503379A (en) | A kind of gas burst emulation mode that is reacted based on adaptive simplifying with grid subdivision | |
CN108625828A (en) | Predict the method and device of perforation blast load output size | |
CN104453879B (en) | The Forecasting Methodology of pressure before boring | |
CN106053016A (en) | Method for pushing plate to create random wave | |
CN105468844A (en) | Analogy method of water-gas coupling transient flow in pipeline | |
CN103198356B (en) | A kind of solid propellant formulation optimization design based on genetic algorithm | |
CN108090268A (en) | The integrated adding method of Seismic Time-history ripple under a kind of Visco-spring Boundary | |
CN104699984A (en) | Method for determining grid sizes for underwater explosion numerical simulation | |
CN102566446B (en) | Method for establishing full-envelope mathematical model of unmanned helicopter based on linear model group | |
CN106354016A (en) | Control method for tamping strength of metal tube | |
Yuan-bin | Genetic algorithm’s application for optimization of PID parameters in dynamic positioning vessel | |
CN107122533A (en) | A kind of water surface heat exchange method for numerical simulation based on EFDC program updates | |
CN106503456A (en) | Based on the Ensemble Kalman Filter Reservoir behavior history-matching method that suprasphere is converted | |
CN104318599A (en) | High-precision fluid animation modeling method based on geometrical features | |
Mukhedkar et al. | Effects of Different Meteorological Standards on Projectile Path. | |
CN103136430A (en) | Solid propellant formula optimization method based on genetic algorithm and energy feature graph | |
Yue et al. | Modelling of a CFD microscale model and its application in wind energy resource assessment | |
CN112560326A (en) | Method and device for determining pressure field | |
CN105223833A (en) | A kind of automatic flight control system parameter regulates and integrated approach |
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 |