CN113340779A - Explosion simulation experiment substance diffusion calculation method and device and storage medium - Google Patents
Explosion simulation experiment substance diffusion calculation method and device and storage medium Download PDFInfo
- Publication number
- CN113340779A CN113340779A CN202110619434.XA CN202110619434A CN113340779A CN 113340779 A CN113340779 A CN 113340779A CN 202110619434 A CN202110619434 A CN 202110619434A CN 113340779 A CN113340779 A CN 113340779A
- Authority
- CN
- China
- Prior art keywords
- light intensity
- formula
- explosion
- equation
- medium
- 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.)
- Pending
Links
- 238000004880 explosion Methods 0.000 title claims abstract description 54
- 239000000126 substance Substances 0.000 title claims abstract description 47
- 238000009792 diffusion process Methods 0.000 title claims abstract description 37
- 238000004088 simulation Methods 0.000 title claims abstract description 36
- 238000004364 calculation method Methods 0.000 title claims description 30
- 239000002245 particle Substances 0.000 claims abstract description 54
- 238000000034 method Methods 0.000 claims abstract description 52
- 238000005474 detonation Methods 0.000 claims description 39
- 239000002360 explosive Substances 0.000 claims description 32
- 230000008569 process Effects 0.000 claims description 28
- 230000001133 acceleration Effects 0.000 claims description 16
- 238000012360 testing method Methods 0.000 claims description 15
- 230000008033 biological extinction Effects 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 12
- 239000006185 dispersion Substances 0.000 claims description 8
- 230000009471 action Effects 0.000 claims description 5
- 239000000443 aerosol Substances 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 5
- 230000035939 shock Effects 0.000 claims description 5
- 230000006870 function Effects 0.000 claims description 4
- 238000000149 argon plasma sintering Methods 0.000 claims description 3
- 239000002085 irritant Substances 0.000 claims description 3
- 231100000021 irritant Toxicity 0.000 claims description 3
- 239000012798 spherical particle Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 238000010521 absorption reaction Methods 0.000 abstract description 5
- 230000002238 attenuated effect Effects 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 7
- 239000000463 material Substances 0.000 description 5
- 239000011257 shell material Substances 0.000 description 5
- VONWDASPFIQPDY-UHFFFAOYSA-N dimethyl methylphosphonate Chemical compound COP(C)(=O)OC VONWDASPFIQPDY-UHFFFAOYSA-N 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- XTFIVUDBNACUBN-UHFFFAOYSA-N 1,3,5-trinitro-1,3,5-triazinane Chemical compound [O-][N+](=O)N1CN([N+]([O-])=O)CN([N+]([O-])=O)C1 XTFIVUDBNACUBN-UHFFFAOYSA-N 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000008187 granular material Substances 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 229920000747 poly(lactic acid) Polymers 0.000 description 2
- 239000004626 polylactic acid Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 231100000167 toxic agent Toxicity 0.000 description 2
- 239000003440 toxic substance Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 231100000419 toxicity Toxicity 0.000 description 1
- 230000001988 toxicity Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
- G01N15/075—Investigating concentration of particle suspensions by optical means
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application provides a method and a device for calculating diffusion of an explosion simulation experiment substance and a storage medium, wherein the method for calculating diffusion of the explosion simulation experiment substance comprises the following steps: determining a scattering area of a medium to be detected, and arranging a laser emitting unit and a light intensity sensor unit in the scattering area; scattering a medium to the scattering area, emitting light by the laser emitting unit, and collecting an incident light intensity signal and an emergent light intensity signal by the light intensity sensor unit; obtaining incident light intensity data and emergent light intensity data according to the incident light intensity signal and the emergent light intensity signal, and obtaining the incident light intensity data and the emergent light intensity data according to a formula 1: i ═ I0e‑τlAnd calculating turbidity data of the throwing area. When the intensity of one beam is I0When the monochromatic parallel light is incident to the medium containing suspended particles, due to the scattering and absorption of the particles to light, the emergent light intensity I can be attenuated to a certain degree, the turbidity of the region to be measured is different, and the distribution condition of the substances in the region can be obtained by utilizing a plurality of light rays at different positions.
Description
Technical Field
The application relates to the technical field of explosion simulation, in particular to a method and a device for calculating diffusion of an explosion simulation experiment substance and a storage medium.
Background
In recent years, the continuous dangerous chemical explosion accidents are integrated with 'burning, explosion and toxicity', and are various, large in quantity and serious in harm. For example, the dangerous chemical super fire and explosion accidents in the new Tianjin coastal region of "8.12" in 2015 and the dangerous chemical warehouse fire and explosion accidents in the lake region of the Shenzhen city of Rou lake region of 8 Rizhen in 1993. The field treatment work of several accidents is powerful, orderly and effective. However, many problems are also exposed in the accident rescue process, such as unclear toxic substance diffusion process under the action of explosion, no grasp on distribution condition and rule of toxic substance concentration, only rough judgment of chemical hazard type and degree, and difficulty in accurate protection in the implementation of chemical emergency rescue task, so that rescue work progresses slowly.
In the prior art, generally, explosion is analyzed in an explosion simulation experiment mode, so that the rescue experience is increased, but in the prior art, a high-speed camera is mostly adopted for shooting, so that the shooting is visual but not specific, but the distribution condition of substances in an environment after explosion is not clear, so that the rapid and efficient development of rescue actions is not facilitated.
Therefore, there is a need in the art for a method, an apparatus and a storage medium for calculating the diffusion of an explosion simulation test substance.
Accordingly, the present application is directed to such a situation.
Disclosure of Invention
The application aims to provide a method for calculating diffusion of an explosion simulation experiment substance, which is used for calculating the distribution condition of the substance in an environment after explosion.
The first aspect of the application provides an explosion simulation experiment substance diffusion calculation method, which comprises the following steps:
determining a scattering area of a medium to be detected, and arranging a laser emitting unit and a light intensity sensor unit in the scattering area;
scattering a medium to the scattering area, emitting light by the laser emitting unit, and collecting an incident light intensity signal and an emergent light intensity signal by the light intensity sensor unit;
the incident light intensity data and the emergent light intensity data are obtained according to the incident light intensity signal and the emergent light intensity signal,
according to equation 1: i ═ I0e-τlCalculating turbidity data of the throwing area; tau-turbidity, l-diameter of the pipe containing particles, I-outgoing light intensity data, I0-incident light intensity data.
By adopting the scheme, when the intensity of one beam is I0When the monochromatic parallel light incides the medium that contains the suspended particles, because the scattering and the absorption of granule light, emergent light intensity I has the decay of certain degree, and the regional turbidity of awaiting measuring is different, then the light attenuation is different, utilizes the formula can calculate turbidity data fast, utilizes many light at different positions, can obtain the regional interior material distribution condition.
The diameter of the particle pipeline is the distance between the laser emission unit and the light intensity sensor unit.
Further, a functional relationship is established, and if there are N spherical particles with a diameter D in a unit volume of the medium to be measured, the turbidity data τ is calculated according to equation 2:
σ — the light-facing area of the particle;
k is extinction coefficient;
d-particle diameter (expressed as D32-Sauter mean diameter).
By adopting the scheme, the turbidity is calculated according to the extinction coefficient, and the calculation accuracy is improved.
Further, the extinction coefficient is a function related to the wavelength λ, the refractive index m, and the particle diameter d, and is characterized by k, λ, m, d, which is calculated according to Mie light scattering theory, and the expression formula 3 is as follows:
in equation 3:
-dimensionless size parameter;
an、bn-Mie coefficients.
By adopting the scheme, the extinction coefficient is calculated.
Further, the calculation of the extinction coefficient can solve the calculation problem by approximating expression formula 4:
ρ — normalized scale factor, defined as equation 5:
by adopting the scheme, the calculation difficulty is reduced, and the calculation efficiency is improved through a formula 4.
Preferably, when 1 < m.ltoreq.1.5, the refractive index correction coefficient k is adoptedmEquation 6:
equation 4 is modified to equation 7: k'e(λ,m,d)＝km×ke(λ,m,d)；
Substituting formula 4, formula 5, and formula 6 into formula 7, formula 8:
by adopting the scheme, the formula is corrected, and the calculation precision is improved.
Preferably, formula 2 and formula 7 are substituted into formula 1, and logarithm is taken at two sides, formula 9:
the particle number N in equation 9 is expressed by the weight concentration C of the particles, i.e., equation 10:
substituting equation 10 into equation 9, equation 11:
the weight concentration C of the particles was calculated using the formula 11.
By adopting the scheme, the medium concentration of each position in the scattering area can be calculated by utilizing the plurality of light beams, the distribution condition of the substances in the environment after the simulated explosion is calculated, the understanding of rescuers on each position in the scattering area is improved, and the rescue efficiency is improved.
Further, the particle size changes with time, and the larger the viscosity of the medium to be measured is, the smaller the size of the dispersed particle size is, the negative correlation is formed; and the more stable and fast the particle size in the whole dispersion process is, the positive correlation is formed, and the rule of the change of the particle size of the medium at different positions along with the time is fitted through a formula 12:
in formula 12, y is the particle size;
x-time, ms;
y0, A1, A2, t1, t 2-constant.
By adopting the scheme, the diameter distribution rule of the mediator in the scattering area can be analyzed, and analysis data can be provided for actual rescue of rescuers.
Further, the distance has no influence on the variation in particle size.
Further, the concentration of the medium changes along with the time, the simulated explosive is detonated in the throwing area, the simulated explosive is filled with the medium to be tested, the concentration of the medium is different at different distances from the center of explosion at the same moment, the change situation of the concentration of the medium along with the time at the same position and different moments is researched, and the medium is fitted through a fitting formula 13:
in formula 13, y is the concentration of the substance, g/cm;
x-time, ms;
y0-initial concentration, g/cm 3;
xc-time to peak, ms;
a, W is a constant.
By adopting the scheme, the change conditions of the medium concentration at the same position and different moments along with the time are fitted through the formula 13, and analysis data are provided for actual rescue of rescuers.
Furthermore, the detonation product does work, the simulated explosive consists of a shell, a stimulant charge and a central explosive charge, and the explosion dispersion process is divided into an acceleration stage, a deceleration stage, a turbulence stage and a diffusion stage; the detonation of the central explosive column and the disintegration time of the shell are very short, the detonation process and the propagation and reflection of shock waves are not considered, the flying of detonation products along the axial direction of the explosive charge are ignored, the expansion of high-temperature and high-pressure detonation products is considered, the overall radial accelerated motion of the stimulant explosive column is promoted, the turbulence of the cloud cluster in the acceleration stage is not significant, so the air action force is mainly the pressure of air on the cloud cluster, in the initial stage, the driving pressure pg of the detonation products is greater than the surrounding pressure pe on the cloud cluster, the cloud cluster does the outward accelerated motion, the sum of the positive work done on the cloud cluster by the driving force and the negative work done on the cloud cluster by the air resistance is the kinetic energy obtained by the cloud cluster, and a formula 14 is obtained according to the work process of the detonation products:
in equation 14: subscript 0 denotes the initial time;
vog is the volume of the central grain;
vg is the volume of detonation product of the central explosive column;
vop is the volume of the stimulant charge;
vp is the sum of the volume Vg of the detonation product and the volume Vl of the aerosol, and the unit is m 3;
ml is the total mass of the irritant and the unit is kg;
ua is the maximum velocity of the stimulant particles, in m/s.
By adopting the scheme, the speed change rule of the radial motion of the explosive cloud cluster in the acceleration stage is analyzed and researched.
Further, from the multi-party equation, the pressure of the detonation product is given by equation 15:
in equation 15: pog is the initial detonation pressure of the center charge in Pa.
By adopting the scheme, the pressure of the detonation product is calculated.
Further, the expression of the maximum velocity ua in the stimulant acceleration phase calculated by combining equation 14 and equation 15 is, equation 16:
by adopting the scheme, the maximum speed of the stimulant acceleration stage is calculated.
Further, since the detonation product volume Vg is much larger than the central charge volume Vog during the acceleration phase, Vg1-n<<Vog1 -nEquation 16 is further expressed as equation 17:
by adopting the scheme, the calculation efficiency is improved, and the explosion result is convenient to analyze.
Further, aerosol refers to a gaseous dispersion system composed of solid or liquid particles suspended in a gaseous medium in an explosion, and the detonation product refers to a substance generated by the explosion.
A second aspect of the present application provides an explosion simulation experimental substance diffusion calculation apparatus, which includes a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor executes the computer program to implement the above-mentioned explosion simulation experimental substance diffusion calculation method.
A third aspect of the present application provides a storage medium comprising one or more programs executable by a processor to perform the above-described explosion simulation experimental material diffusion calculation method.
In summary, the present application has the following beneficial effects:
1. according to the method for calculating the diffusion of the explosion simulation experiment substance, when the intensity of one beam is I0When the monochromatic parallel light is incident to a medium containing suspended particles, due to the scattering and absorption effects of the particles on the light, the emergent light intensity I is attenuated to a certain degree, the turbidities of the areas to be detected are different, the light intensity attenuation is different, the turbidity data can be quickly calculated by using a formula, and the distribution condition of the substances in the areas can be obtained by using a plurality of light rays at different positions;
2. according to the method for calculating the diffusion of the explosion simulation experiment substance, the medium concentration of each position in the scattering area can be calculated by utilizing a plurality of light beams, the substance distribution condition in the environment after the simulation explosion is calculated, the understanding of rescuers on each position in the diffusion area is improved, and the rescue efficiency is improved;
3. according to the method for calculating the diffusion of the explosion simulation experiment substances, the diameter distribution rule of the mediator in the scattering area can be obtained through analysis, and analysis data are provided for actual rescue of rescuers;
4. according to the method for calculating the diffusion of the explosion simulation experiment substances, the change conditions of the medium concentration at the same position and different moments along with the time are fitted through a formula 13, and analysis data are provided for actual rescue of rescuers.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of an embodiment of a method for calculating the diffusion of an explosion simulation test substance according to the present application;
FIG. 2 is a schematic diagram of the experimental apparatus connection according to an embodiment of the method for calculating the diffusion of the experimental substance for explosion simulation of the present application;
FIG. 3 is a partial connection view of FIG. 2;
fig. 4 is a schematic diagram of the structure of the simulated explosive.
Detailed Description
Reference will now be made in detail to the exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, like numbers in different drawings represent the same or similar elements unless otherwise indicated. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with certain aspects of the present application, as detailed in the appended claims.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used in this application and the appended claims, the singular forms "a", "an", and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any and all possible combinations of one or more of the associated listed items.
As shown in fig. 1, a first aspect of the present application provides a method for calculating a diffusion of an explosion simulation test substance, including the following steps:
determining a scattering area of a medium to be detected, and arranging a laser emitting unit and a light intensity sensor unit in the scattering area;
as shown in fig. 2 and 3, in the specific implementation process, the light intensity sensor unit may be an RS485 light intensity sensor, the laser emitting unit may be a laser emitter, and the apparatus in the experiment further includes a controller, a synchronous control device, a high-speed camera, a shock tube for scattering a medium to be measured, an electromagnetic valve, and a high-pressure air supply device.
In the specific implementation process, the medium to be detected can be scattered through the shock tube, the electromagnetic valve and the high-pressure air supply device, and the medium can also be scattered through simulating the explosion of explosives.
The dispersed cloud in fig. 3 is the scattered medium, and when the light path passes through the cloud, the light intensity value is attenuated based on the absorption and scattering of the cloud to the light, and the sensor obtains the intensity values before and after the light intensity attenuation.
Scattering a medium to the scattering area, emitting light by the laser emitting unit, and collecting an incident light intensity signal and an emergent light intensity signal by the light intensity sensor unit;
the incident light intensity data and the emergent light intensity data are obtained according to the incident light intensity signal and the emergent light intensity signal,
in the specific implementation process, the incident light intensity data and the emergent light intensity data are obtained according to the incident light intensity signal and the emergent light intensity signal and are realized through a sensor and a computer.
According to equation 1: i ═ I0e-τlCalculating turbidity data of the throwing area; tau-turbidity, l-diameter of the pipe containing particles, I-outgoing light intensity data, I0-incident light intensity data.
Said formula 1 is based on the Lambert-Beer Light transmission law (Light transmission law).
The particle conduit diameter may be a pre-measured value.
By adopting the scheme, when the intensity of one beam is I0When the monochromatic parallel light incides the medium that contains the suspended particles, because the scattering and the absorption of granule light, emergent light intensity I has the decay of certain degree, and the regional turbidity of awaiting measuring is different, then the light attenuation is different, utilizes the formula can calculate turbidity data fast, utilizes many light at different positions, can obtain the regional interior material distribution condition.
The turbidity is the degree of obstruction of the medium to the passage of light.
In the specific implementation process, a functional relationship is established, and if N spherical particles with the diameter of D exist in the medium to be measured in unit volume, the turbidity data τ is obtained according to the formula 2:
σ — the light-facing area of the particle;
k is extinction coefficient;
d-particle diameter (expressed as D32-Sauter mean diameter).
By adopting the scheme, the turbidity is calculated according to the extinction coefficient, and the calculation accuracy is improved.
In a specific implementation process, the extinction coefficient is a function related to a wavelength λ, a refractive index m and a particle size d, and is characterized by k, λ, m and d, and is calculated according to Mie light scattering theory, and an expression formula 3 is as follows:
in equation 3:
-dimensionless size parameter;
an、bn-Mie coefficients.
By adopting the scheme, the extinction coefficient is calculated.
In a specific implementation, the calculation of the extinction coefficient can solve the calculation problem through an approximate expression formula 4:
ρ — normalized scale factor, defined as equation 5:
by adopting the scheme, the calculation difficulty is reduced, and the calculation efficiency is improved through a formula 4.
In a preferred embodiment of the invention, the refractive index modification factor k is used when 1 < m.ltoreq.1.5mEquation 6:
in specific implementations, m can be 1.2 or 1.5, and the like.
Equation 4 is modified to equation 7: k'e(λ,m,d)＝km×ke(λ,m,d)；
Substituting formula 4, formula 5, and formula 6 into formula 7, formula 8:
by adopting the scheme, the formula is corrected, and the calculation precision is improved.
In a preferred embodiment of the present invention, formula 2 and formula 7 are substituted into formula 1, and logarithm is taken on both sides, formula 9:
the particle number N in equation 9 is expressed by the weight concentration C of the particles, i.e., equation 10:
substituting equation 10 into equation 9, equation 11:
the weight concentration C of the particles was calculated using the formula 11.
By adopting the scheme, the medium concentration of each position in the scattering area can be calculated by utilizing the plurality of light beams, the distribution condition of the substances in the environment after the simulated explosion is calculated, the understanding of rescuers on each position in the scattering area is improved, and the rescue efficiency is improved.
In the specific implementation process, the particle size changes along with time, and the larger the viscosity of the medium to be measured is, the smaller the size of the dispersed particle size is, the negative correlation is formed; and the more stable and fast the particle size in the whole dispersion process is, the positive correlation is formed, and the rule of the change of the particle size of the medium at different positions along with the time is fitted through a formula 12:
in formula 12, y is the particle size;
x-time, ms;
y0, A1, A2, t1, t 2-constant.
By adopting the scheme, the diameter distribution rule of the mediator in the scattering area can be analyzed, and analysis data can be provided for actual rescue of rescuers.
In the specific implementation process, the time-varying exponential fitting parameters of the particle diameters of water and dimethyl methylphosphonate DMMP at different positions under 2.6MPa are shown in table 1:
TABLE 1
In the specific implementation process, the medium concentration changes along with time, the simulated explosive is detonated in the throwing area, the simulated explosive is filled with the medium to be tested, and because the medium concentrations at different distances from the center of explosion at the same moment are different, the change conditions of the medium concentrations along with time at the same position and different moments are researched, and the medium concentration is fitted through a fitting formula 13:
in formula 13, y is the concentration of the substance, g/cm;
x-time, ms;
y0-initial concentration, g/cm 3;
xc-time to peak, ms;
a, W is a constant.
By adopting the scheme, the change conditions of the medium concentration at the same position and different moments along with the time are fitted through the formula 13, and analysis data are provided for actual rescue of rescuers.
Table 2-2.6MPa concentrations of Water and DMMP at different positions over time, fitting parameters for equation 13
Table 3-200ms DMMP at different distances and different driving loads in concentration distribution
In the specific implementation process, the detonation product does work, the simulated explosive consists of a shell, a stimulant charge and a central explosive charge, and the explosion dispersion process is divided into an acceleration stage, a deceleration stage, a turbulence stage and a diffusion stage; the detonation of the central explosive column and the disintegration time of the shell are very short, the detonation process and the propagation and reflection of shock waves are not considered, the flying of detonation products along the axial direction of the explosive charge are ignored, the expansion of high-temperature and high-pressure detonation products is considered, the overall radial accelerated motion of the stimulant explosive column is promoted, the turbulence of the cloud cluster in the acceleration stage is not significant, so the air action force is mainly the pressure of air on the cloud cluster, in the initial stage, the driving pressure pg of the detonation products is greater than the surrounding pressure pe on the cloud cluster, the cloud cluster does the outward accelerated motion, the sum of the positive work done on the cloud cluster by the driving force and the negative work done on the cloud cluster by the air resistance is the kinetic energy obtained by the cloud cluster, and a formula 14 is obtained according to the work process of the detonation products:
in equation 14: subscript 0 denotes the initial time;
vog is the volume of the central grain;
vg is the volume of detonation product of the central explosive column;
vop is the volume of the stimulant charge;
vp is the sum of the volume Vg of the detonation product and the volume Vl of the aerosol, and the unit is m 3;
ml is the total mass of the irritant and the unit is kg;
ua is the maximum velocity of the stimulant particles, in m/s.
By adopting the scheme, the speed change rule of the radial motion of the explosive cloud cluster in the acceleration stage is analyzed and researched.
In the specific implementation process, when an explosive simulation experiment is carried out, a movable shelter, a detonation control device, a wind direction anemoscope, a curtain, a support and the like can be installed for assistance.
As shown in fig. 4, in a specific implementation, the shell material of the simulated explosive may be PLA (polylactic acid) material, the thickness is 2.3mm, the stimulant charge may be CS agent, the central explosive charge may be passivated hexogen, and specific parameters may be as shown in table 4:
TABLE 4
In the specific implementation process, the pressure of the detonation product is obtained by the multi-party equation, and the formula 15 is as follows:
in equation 15: pog is the initial detonation pressure of the center charge in Pa.
By adopting the scheme, the pressure of the detonation product is calculated.
In the specific implementation process, the expression of the maximum velocity ua in the stimulant acceleration phase is obtained by the combined calculation of formula 14 and formula 15, and is formula 16:
by adopting the scheme, the maximum speed of the stimulant acceleration stage is calculated.
In the specific implementation process, as the volume Vg of detonation products in the acceleration stage is far larger than the volume Vog of the central explosive column, the volume Vg1-n<<Vog1-nEquation 16 is further expressed as equation 17:
by adopting the scheme, the calculation efficiency is improved, and the explosion result is convenient to analyze.
In the specific implementation process, the explosion is recorded by shooting through a high-speed camera, and the recording can be used for recording explosion pictures at various time points in the explosion.
In the specific implementation process, the maximum speed value of the acceleration stage is 1495m/s obtained by analyzing the picture, the volume of the cloud cluster is approximately processed into the volume of a cylinder by reading the diameter and the height of the cloud cluster photo when the ground explosion is 0.7ms, and V is obtainedpIs 0.35m3From the basic parameters of the test projectile, the ammunition volume VopIs 9.9 multiplied by 10-5m3Volume V of central columnogIs 6.8 multiplied by 10-6m3The detonation pressure p is selected according to the content and packing density of the passivated hexogen usedogIs 28.7X 109Pa, pe may be 1X 105Pa, ml can be 0.06kg, the polytropic index n is combined with the practical approximation 3 to solve the maximum speed u of the acceleration stageaThe theoretical calculation value is 1444m/s, and the theoretical calculation value is better matched with the test value.
In the specific implementation, aerosol refers to a gaseous dispersion system composed of solid or liquid particles suspended in a gaseous medium in an explosion, and the detonation product refers to a substance generated by the explosion.
A second aspect of the present application provides an explosion simulation experimental substance diffusion calculation apparatus, which includes a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor executes the computer program to implement the above-mentioned explosion simulation experimental substance diffusion calculation method.
A third aspect of the present application provides a storage medium comprising one or more programs executable by a processor to perform the above-described explosion simulation experimental material diffusion calculation method.
It should be noted that, for those skilled in the art, without departing from the principle of the present application, several improvements and modifications can be made to the present application, and these improvements and modifications also fall into the protection scope of the claims of the present application.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
It should be understood that the technical problems can be solved by combining and combining the features of the embodiments from the claims.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application or portions thereof that substantially contribute to the prior art may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the application. Thus, the present application is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (10)
1. A method for calculating the diffusion of an explosion simulation experiment substance is characterized by comprising the following steps:
determining a scattering area of a medium to be detected, and arranging a laser emitting unit and a light intensity sensor unit in the scattering area;
scattering a medium to the scattering area, emitting light by the laser emitting unit, and collecting an incident light intensity signal and an emergent light intensity signal by the light intensity sensor unit;
the incident light intensity data and the emergent light intensity data are obtained according to the incident light intensity signal and the emergent light intensity signal,
according to equation 1: i ═ I0e-τlCalculating turbidity data of the throwing area; tau-turbidity, l-diameter of the pipe containing particles, I-outgoing light intensity data, I0-incident light intensity data.
2. The method for calculating diffusion of an explosion-simulation test substance according to claim 1, wherein: establishing a functional relation, wherein if N spherical particles with the diameter D exist in the medium to be detected in unit volume, the turbidity data tau is obtained according to a formula 2:
σ — the light-facing area of the particle;
k is extinction coefficient;
d-particle diameter (expressed as D32-Sauter mean diameter).
3. The method for calculating diffusion of an explosion-simulation test substance according to claim 2, wherein: the extinction coefficient is a function related to wavelength lambda, refractive index m and particle size d, is characterized by k, lambda, m and d, and is calculated according to Mie light scattering theory, and the expression formula 3 is as follows:
in equation 3:
-dimensionlessA dimensional parameter;
an、bn-Mie coefficients.
4. The method for calculating diffusion of an explosion-simulation test substance according to any one of claims 1 to 3, wherein: the calculation of the extinction coefficient can solve the calculation problem by approximating expression formula 4:
ρ — normalized scale factor, defined as equation 5:
5. the method for calculating diffusion of an explosion-simulation test substance according to claim 4, wherein: when m is more than 1 and less than or equal to 1.5, a refractive index correction coefficient k is adoptedmEquation 6:
equation 4 is modified to equation 7: k'e(λ,m,d)＝km×ke(λ,m,d)；
Substituting formula 4, formula 5, and formula 6 into formula 7, formula 8:
6. the method for calculating diffusion of an explosion-simulation test substance according to claim 5, wherein: calculating the weight concentration, substituting formula 2 and formula 7 into formula 1, and taking logarithm at two sides to obtain formula 9:
the particle number N in equation 9 is expressed by the weight concentration C of the particles, i.e., equation 10:
substituting equation 10 into equation 9, equation 11:
the weight concentration C of the particles was calculated using the formula 11.
7. The method for calculating diffusion of an explosion-simulation test substance according to claim 2 or 6, wherein: the particle size changes along with time, and is in negative correlation as the viscosity of the medium to be measured is higher and the size of the dispersed particle size is smaller; and the more stable and fast the particle size in the whole dispersion process is, the positive correlation is formed, the rule of the change of the particle size of the medium at different positions along with the time is calculated by fitting according to a formula 12:
y-particle size;
x-time, ms;
y0, A1, A2, t1, t 2-constant.
8. The method for calculating diffusion of an explosion-simulation test substance according to claim 7, wherein: the concentration of the medium changes with time, the simulated explosive is detonated in the throwing area, the simulated explosive is filled with the medium to be tested, because the concentrations of the medium are different at different distances from the center of explosion at the same moment, the change situation of the concentration of the medium with time at the same position and different moments is researched, fitting is carried out through a fitting formula 13, and the change situation of the concentration of the medium with time at the same position and different moments is calculated:
in formula 13, y is the concentration of the substance, g/cm;
x-time, ms;
y0-initial concentration, g/cm 3;
xc-time to peak, ms;
a, W is a constant.
9. The method for calculating diffusion of an explosion-simulation test substance according to claim 8, wherein: the detonation product does work, the simulated explosive consists of a shell, a stimulant charge column and a central explosive column, and the explosion dispersion process is divided into an acceleration stage, a deceleration stage, a turbulence stage and a diffusion stage; the detonation of the central explosive column and the disintegration time of the shell are very short, the detonation process and the propagation and reflection of shock waves are not considered, the flying of detonation products along the axial direction of the explosive charge are ignored, the expansion of high-temperature and high-pressure detonation products is considered, the overall radial accelerated motion of the stimulant explosive column is promoted, the turbulence of the cloud cluster in the acceleration stage is not significant, so the air action force is mainly the pressure of air on the cloud cluster, in the initial stage, the driving pressure pg of the detonation products is greater than the surrounding pressure pe on the cloud cluster, the cloud cluster does the outward accelerated motion, the sum of the positive work done on the cloud cluster by the driving force and the negative work done on the cloud cluster by the air resistance is the kinetic energy obtained by the cloud cluster, and a formula 14 is obtained according to the work process of the detonation products:
in equation 14: subscript 0 denotes the initial time;
vog is the volume of the central grain;
vg is the volume of detonation product of the central explosive column;
vop is the volume of the stimulant charge;
vp is the sum of the volume Vg of the detonation product and the volume Vl of the aerosol, and the unit is m 3;
ml is the total mass of the irritant and the unit is kg;
ua is the maximum velocity of the stimulant particles, in m/s.
10. An explosion simulation test substance diffusion calculation apparatus, comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor executes the program to implement the explosion simulation test substance diffusion calculation method according to any one of claims 1 to 9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110619434.XA CN113340779A (en) | 2021-06-03 | 2021-06-03 | Explosion simulation experiment substance diffusion calculation method and device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110619434.XA CN113340779A (en) | 2021-06-03 | 2021-06-03 | Explosion simulation experiment substance diffusion calculation method and device and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113340779A true CN113340779A (en) | 2021-09-03 |
Family
ID=77474732
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110619434.XA Pending CN113340779A (en) | 2021-06-03 | 2021-06-03 | Explosion simulation experiment substance diffusion calculation method and device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113340779A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116818887A (en) * | 2023-07-07 | 2023-09-29 | 大唐环境产业集团股份有限公司 | Particle two-phase flow concentration measuring method and device and computer equipment |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101929940A (en) * | 2009-06-26 | 2010-12-29 | 上海市上海中学 | Embedded system-based method and device for online measurement of micron particles |
| CN103592103A (en) * | 2013-11-11 | 2014-02-19 | 浙江大学 | Mini-channel liquid-solid two-phase flow parameter measurement device and method based on laser extinction method |
| CN203616232U (en) * | 2013-12-09 | 2014-05-28 | 太原科技大学 | Device for measuring average particle size and concentration of particles in atmosphere |
| CN105067489A (en) * | 2015-08-07 | 2015-11-18 | 中国科学院计算技术研究所 | Dynamic light scattering technology-based suspended particle particle size measuring device, and method |
| CN107132545A (en) * | 2017-04-28 | 2017-09-05 | 北京理工大学 | The analysis method that particulate matter in aerosol influences on the quality of remotely sensed image |
| CN210953732U (en) * | 2019-10-19 | 2020-07-07 | 山西奥斯帕克信息科技有限公司 | Mining dust concentration sensor |
| CN112129439A (en) * | 2020-10-16 | 2020-12-25 | 上海欣诺通信技术股份有限公司 | Torsion sensor based on optical fiber light scattering induction |
-
2021
- 2021-06-03 CN CN202110619434.XA patent/CN113340779A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101929940A (en) * | 2009-06-26 | 2010-12-29 | 上海市上海中学 | Embedded system-based method and device for online measurement of micron particles |
| CN103592103A (en) * | 2013-11-11 | 2014-02-19 | 浙江大学 | Mini-channel liquid-solid two-phase flow parameter measurement device and method based on laser extinction method |
| CN203616232U (en) * | 2013-12-09 | 2014-05-28 | 太原科技大学 | Device for measuring average particle size and concentration of particles in atmosphere |
| CN105067489A (en) * | 2015-08-07 | 2015-11-18 | 中国科学院计算技术研究所 | Dynamic light scattering technology-based suspended particle particle size measuring device, and method |
| CN107132545A (en) * | 2017-04-28 | 2017-09-05 | 北京理工大学 | The analysis method that particulate matter in aerosol influences on the quality of remotely sensed image |
| CN210953732U (en) * | 2019-10-19 | 2020-07-07 | 山西奥斯帕克信息科技有限公司 | Mining dust concentration sensor |
| CN112129439A (en) * | 2020-10-16 | 2020-12-25 | 上海欣诺通信技术股份有限公司 | Torsion sensor based on optical fiber light scattering induction |
Non-Patent Citations (6)
| Title |
|---|
| 孔明等: "改进的LMS算法在光全散射法粒度反演中的应用", 《光学技术》 * |
| 李廷等: "基于分装结构的CS刺激剂爆炸分散过程研究", 《含能材料》 * |
| 赵建华等: "多波长激光全散射法探测火灾烟雾的研究", 《应用激光》 * |
| 赵建华等: "烟雾浓度及粒径的激光探测实验研究", 《激光与红外》 * |
| 郑刚等: "颗粒浓度及粒度的光散射在线测量", 《中国激光》 * |
| 郑刚等: "颗粒浓度在线监测的双波长消光法", 《仪器仪表学报》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116818887A (en) * | 2023-07-07 | 2023-09-29 | 大唐环境产业集团股份有限公司 | Particle two-phase flow concentration measuring method and device and computer equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Settles | High-speed imaging of shock waves, explosions and gunshots: new digital video technology, combined with some classic imaging techniques, reveals shock waves as never before | |
| Knock et al. | Predicting blast waves from the axial direction of a cylindrical charge | |
| CN113340779A (en) | Explosion simulation experiment substance diffusion calculation method and device and storage medium | |
| Wang et al. | An experimental study on thermal radiation of fire whirl | |
| Tachajapong et al. | Experimental modelling of crown fire initiation in open and closed shrubland systems | |
| Schunck et al. | Blast mitigation by water mist: the effect of the detonation configuration | |
| Liang et al. | Improvement of aerosol spray scavenging efficiency with water mist | |
| Gouhier et al. | The geometry of Strombolian explosions: Insights from Doppler radar measurements | |
| Goroshin et al. | Measurement of particle density during explosive particle dispersal | |
| Vigne et al. | Influence of variability of soot yield parameter in assessing the safe evacuation conditions in advanced modeling analysis. Results of physical and numerical modeling comparison | |
| CN205449753U (en) | Two dimension smog concentration field measuring device based on sheet laser | |
| Phylaktou et al. | Flame speed measurements in dust explosions | |
| Fu et al. | Dynamic concentration measurement of micro/nano aluminum powder | |
| Berg et al. | Investigation of dust dispersion in a modified Hartmann tube using positron emission particle tracking and simulations | |
| Avelar et al. | Wind tunnel simulation of the atmospheric boundary layer for studying the Wind pattern at Centro de Lançamento de Alcântara | |
| Bagheri | Numerical and experimental investigation of particle terminal velocity and aggregation in volcanic plumes | |
| Schweitzer et al. | Partial Confinement of Detonation Products by Shock Reflection from a Convergent Nozzle Opening | |
| Wang et al. | Aerodynamic characteristics calculation and diffusion law analysis of rectangular-chaff clouds under airflow | |
| Cheng et al. | Physicochemical characterization of Capstone depleted uranium aerosols II: particle size distributions as a function of time | |
| Patel et al. | Theoretical and experimental evaluation of a compact aerosol wind-tunnel and its application for performance investigation of particulate matter instruments | |
| Chen et al. | Predictive model for concentration distribution of explosive dispersal | |
| MIKART et al. | Capture of particles dispersed by detonation using an aqueous foam confinement | |
| Wang et al. | A Computing Model of Foil-surface-type IR Decoy Radiation Based on Pyrophoric Multi-hole Activated Metal Combustion Model | |
| Igra | My Walk in the Field of Shock Waves | |
| Dygdała et al. | Aerosol produced by explosive detonation |
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 |