Advances in Condensed Matter Physics

Volume 2014 (2014), Article ID 219547, 8 pages

http://dx.doi.org/10.1155/2014/219547

Research Article

Defects-Induced Hot Spots in TATB

Zhonghua Yan,1,2 Chuanchao Zhang,¹ Hongwei Yan,¹ Zhijie Li,² Li Li,² Ming Huang,³ Bisheng Tan,³ and Xiaodong Yuan¹

¹Research Center of Laser Fusion, China Academy of Engineering Physics, Mianyang 621900, China

²School of Physical Electronics, University of Electronic Science and Technology of China, Chengdu 610054, China

³Institute of Chemical Materials, China Academy of Engineering Physics, Mianyang 621900, China

Received 1 March 2014; Accepted 23 May 2014; Published 17 June 2014

Academic Editor: Kai Sun

Copyright © 2014 Zhonghua Yan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We investigate the interaction between the laser and energetic materials with different defects. The three-dimensional models of triaminotrinitrobenzene (TATB) explosives containing spherical pores, craters, and cracks are established, respectively. The laser ignition process of TATB is simulated with three-dimensional finite difference time domain (3D-FDTD) method to study the electromagnetic field distribution surrounding these defects with 355 nm laser incidence. It indicates that the larger defects in the TATB energetic materials have the stronger electric field modulations to initial incident laser for all the three defects, which is easier to lead to the generation of hot spots. Furthermore, TATB materials with spherical pore defects and crater defects are easier to form hot spots than those with narrow crack defects.

1. Introduction

The ignition of explosives has been studied for so many years. There are several conventional means for the ignition systems, such as electric spark ignition and electric bridge-wire ignition. However, due to the interference of static electricity, radio frequency, electromagnetic field, and other unpredictable factors, traditional ignition technologies of explosives materials tend to explode by mistake or fail to work.

Laser ignition technology, introduced by Brish et al. [1] and Yang and Menichelli [2] in the late sixties and early seventies of last century, has been recognized as a safe and reliable means which has extensive potential applications in the aerospace and military filed owing to the laser’s high intensity and good direction. The first discussion of laser ignition mechanisms is also discussed in [1]. Generally speaking, there are two possible ignition mechanisms used to interpret the explosion phenomenon: thermal ignition mechanism and impulse ignition mechanism. In the thermal mechanism, the intense local heating by the laser pulse causes the rate of chemical reactions within the materials to accelerate sharply. It requires flux densities in the neighborhood of 1 kW/cm² or less and pulse durations of about 1 ms or longer. But in the impulse ignition mechanism, where the chemical reactions are accelerated by a strong laser-generated shock wave, it needs more intense and shorter duration by factors of about 10⁶ with the similar energy content than that in the thermal ignition mechanism [3]. The term “hot spots” represents the regions which can couple efficiently the locally higher temperature and ultimately initiation in materials [4]. Because of their modulations for the laser electromagnetic field, defects (voids, bubbles, etc.) in the explosives materials can result in the generation of hot spots.

Based on the thermal ignition mechanism, we investigated the interaction between the laser and energetic materials with different defects in our work. Triaminotrinitrobenzene (TATB) explosives were firstly characterized with the Scanning Electron Microscopy (SEM). Three-dimensional (3D) models of TATB energetic materials containing spherical pore, crater, and crack defects were established, respectively. The laser ignition process of TATB was simulated with three-dimensional finite difference time domain (3D-FDTD) method [5–8] to study the electromagnetic field distribution (in proportion to the flux density ) surrounding these defects with 355 nm laser incidence. The simulated results of electric field amplitude distributions and the number of high-intensity spots in TATB materials were analyzed to investigate the effect of the three kinds of defects on the generation of hot spots.

2. Models and FDTD Method

The TATB explosives with many microdefects, including spherical pores, craters, and cracks, were characterized with the Scanning Electron Microscopy (SEM) shown in Figure 1. It indicated that all the defect sizes are in the micron or submicron level (shown in Table 1). As it is difficult to know the light modulation of only one of these defects through experiments, in this work we utilized the 3D-FDTD methods to simulate the modulation of each defect to the incident laser. According to these SEM images, we constructed the three-dimensional sphere model, rotated parabolic body model, and straight triangular prism model, respectively, to fit the air area as encircled in Figure 1, as their three-dimensional profiles showed in Figure 2. We supposed that all the defects lie in the top surface of the TATB explosives. Their equations in Cartesian coordinates are as follows.

Table 1: The size range of different defects.

Defects
Size range

SEM
FDTD simulations

Spherical pores
–1.5 μm
–60δ

Craters ()
–3 μm
–100δ

Cracks ()
–3 μm
–100δ,

(a)

(b)

(c)

Figure 1: SEM of TATB explosives with different defects: (a) spherical pores, (b) craters, and (c) cracks.

(a)

(b)

(c)

Figure 2: The three-dimensional profiles of models for the microdefects in Figure 1: (a) sphere, (b) rotated parabolic body, and (c) straight triangular prism.

Sphere

Rotated parabolic body

Straight triangular prism

where is the coordinate of the symmetrical center of the top surface. and represent the radius and the depth of the spherical pores, respectively. and are the width and depth of the craters. , and represent the length, the depth and the width of the cracks respectively. From the SEM, we assumed that , . For easy description, in our study we supposed that the depth/radius ratio of spherical pore defects , the depth/width ratio of crater defects , and the depth/width ratio of crack defects .

The finite difference time domain (FDTD) method, proposed by Yee in 1966 [9], has been proven to be one of the most popular and powerful numerical tools for simulating the propagation and scattering of classical electromagnetic fields in materials. It is a solution method of the partial derivatives in Maxwell’s equations by discretizing these equations in time and space domains using the difference scheme on edges of the “Yee-cell” [9], which need to divide space and time into a regular grid. The field solution at the current time step is deduced from the field values at the previous time steps with the recursive time-marching algorithm.

In our study, the overall lattice in the TATB explosives is a uniformly gridded cuboid simulation domain as a dimension of cells with the same unit cell size for the three directions and a time step to satisfy computational stability and numerical dispersion conditions, where is the 355 nm laser wavelength and is the velocity of light in the air. The time step number is 1600 for all three models. The normal incidence laser was from the bottom of the model with the electric field amplitude () 1.0 V/m shown in Figure 2. The relative dielectric constants () of air and TATB are 1.0 and 4.0, respectively [10]. The boundaries are terminated by a thickness perfectly matched layer (PML) [11].

3. Results and Discussion

Equation (4) expounds the relationship of the flux density and electric field amplitude . As we know, the flux density distribution can be analyzed by the electric field amplitude in the TATB explosives

3.1. Spherical Pore Defects

In order to explore effects of the radius of spherical pore defects and their depth embedded in TATB explosives on the flux density distribution, we investigated the relationship between the electric field distribution and the depth/radius ratio Rs. The size range of FDTD simulations is shown in Table 1. Figure 3 shows the maximum electric field of some planes at different Rs when . As the location of plane changes, fluctuates up and down at different planes. However, the general trend increases gradually. Figure 4(a) shows the maximum electric field of Plane at the radius and the whole simulation domain at , , and changing with the increasing Rs. It is obvious that   of all the simulation domain is not obtained definitely at Plane . It is even much larger than that in the Plane . As Rs increases from 0.125 to 2.5, for the whole space changes in the range of 2.0–2.7 V/m, and the flux density reaches a factor (light intensity enhancement factor, LIEF) of 4–7.29 compared with the initial incident laser. It is clear that the larger the radius of sphere defects is, the stronger its modulation to incident laser is. When the size of spherical pores is small, the modulation effect is weak because of the minor scattering surface. However, when gets to , the reflecting surface increases owing to the increasing surface area. At this time, even total reflection may appear. So the modulation to initial laser is stronger when sphere is larger.

Figure 3: distributions at different planes with various Rs at .

(a)

(b)

Figure 4: (a) and (b) the number of high-intensity spots at various Rs.

Besides, we use the term “high-intensity spot” in this work to express the region where the electric field amplitude  V/m and the number of high-intensity spots was counted and exported by computer program. Figure 4(b) shows the number of high-intensity spots in the whole simulation space. It revealed that the number increases firstly with the increase of Rs and then decreases to a stable range when , in which the defects have been completely embedded into the explosives. As the radius increases, the number of high-intensity spots increases. It shows that the number of high-intensity spots induced by spherical pore defects when  μm () has reached about 3 times of that when  μm .

3.2. Crater Defects

The interaction between laser and TATB explosives containing crater defects is also studied. Figure 5 shows the effects of the width and depth of crater defects (shown in Figure 2(b)) on the electric field distributions. Generally speaking, the maximum electric field increases firstly with the increasing depth/width ratio Rcr and then decreases slowly. The larger the width of crater defects is, the stronger its modulation to incident laser is. When  μm , keeps basically in a stable value around 2.08 V/m and there is a small number of high-intensity spots. However, when gets to 1.8 μm (), is mainly in the range of 2.1–2.3 V/m with the LIEF from 4.41 to 5.29 as Rcr changes. The number of high-intensity spots increases to dozens. Finally, when reaches up to 3 μm , most of the maximum electric field keep in the level above 2.55 V/m with the LIEF more than 6.5 times compared to the initial incident laser. Similar to the spherical pore defects, as the width of crater defects increases, the number of high-intensity spots increases. It shows that the number of high-intensity spots when  μm has reached about 4 times of that when  μm.

(a)  

(b)

Figure 5: (a) and (b) the number of high-intensity spots as the function of Rcr.

Figure 6, respectively, displays the electric field amplitude distributions of Plane and Plane z198 (where the maximum electric field amplitude occurs) when , . As the defects lie in the top surface, the electric field near the top surface is higher. And the maximum electric field appears in Plane z198 which is nearly close to crater defects.

(a)  

(b)

Figure 6: Electric field amplitude distributions when and : (a) Plane and (b) Plane z198 where the maximum electric field amplitude occurs.

3.3. Crack Defects

According to Figure 2(c), we investigated the effect of the depth/width ratio Rck of crack defects with various length , , and on the propagation and scattering of electromagnetic fields within the TATB explosives. From the SEM image in Figure 1(c), we assumed the width . Figure 7(a) displays the maximum electric field amplitude changing with the increasing Rck. As the length of crack defects changes, most lie in the range of 2.1–2.3 V/m with the number of high-intensity spots less than 30. As the length of crack defects increases, there are also some enhanced electric field areas. However, it is not so many as those of crater and crack defects because of the narrow width.

(a)  

(b)

Figure 7: (a) and (b) the number of high-intensity spots as the function of Rck.

The Rck for gaining the maximum of is different for different lengths. When is about 0.6 μm , the maximum, 2.47 V/m, is obtained at . However, for  μm and 3 μm , Rck changes to 4.5 and 2.5 with their values to 2.73 V/m and 2.54 V/m, respectively. Even so, the maximal number of high intensity spots is acquired at when . Figure 8 shows the electric field amplitude distributions in the plane of at the length and width . It indicates that every lies in the plane which is on (Plane z200) or near (Plane z187, z193) the top surface. Especially, the lie in the plane of the bottom of crack defects (Plane z141 for , Plane z121 for ).

(a)

(b)

(c)

(d)

(e)

Figure 8: Electric field amplitude distributions in the plane of of cracks with different Rck at the length and width : (a) , Plane z200; (b) , Plane z187; (c) , Plane z141; (d) , Plane z121; (e) , Plane z193.

From Figures 4, 5, and 7, we can deduce that the larger the defect in the TATB energetic material is, the stronger modulation to incident laser is, which is easier to generate hot spots. Moreover, spherical pore defects when  μm with Rs from 0.75 to 1.75 and crater defects when  μm with Rcr from 0.4 to 1.0 are easier to lead to the generation of hot spots than crack defects with because of the narrow width.

4. Conclusions

The laser interactions with TATB explosives including different defects (spherical pores, craters, and cracks) were investigated, respectively, with the three-dimensional finite difference time domain method. The modulations of these defects to incident laser were analyzed with the maximum electric field amplitude and the number of high-intensity spots.

The results indicate that the maximum electric field of spherical pores changes in the range of 2.0–2.7 V/m as the radius changes from 0.6 to 1.8 μm with the depth/radius ratio Rs increasing from 0.125 to 2.5. And the light intensity enhancement factor (LIEF) reaches 4–7.29 compared with the initial laser. As the radius increases, electric field amplitude strengthens and the number of high-intensity spots increases. The number of high-intensity spots induced by spherical pore defects when  μm reaches about 3 times of that when  μm. Similar to the spherical pore defects, as the width of crater defects increases from 0.6 to 3 μm, the number of high-intensity spots increases to several hundred. When  μm, the number of high-intensity spots has reached about 4 times of that when  μm with most of the maximum electric field in the level above 2.55 V/m (LIEF ≥ 6.5). However, as the length of crack defects changes from 0.6 to 3 μm with the width , most lie in the range of 2.1–2.3 V/m with the number of high-intensity spots less than 30 because of the narrow width.

For all the three defects, the larger the defect in the TATB energetic material is, the stronger the modulation to incident laser is, which is easier to lead to the generation of hot spots. Furthermore, spherical pore defects when  μm with depth/radius ratio from 0.75 to 1.75 and crater defects when  μm with depth/width ratio from 0.5 to 1.0 are easier ways to generate hot spots than crack defects with .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study was financially supported by the National Natural Science Foundation of China (Grant no. 11372289) and the NSAF Joint Foundation of China (Grant no. U1330108).

References

A. A. Brish, I. A. Galeev, B. N. Zaitsev, E. A. Sbitnev, and L. V. Tatarintsev, “Mechanism of initiation of condensed explosives by laser radiation,” Combustion, Explosion, and Shock Waves, vol. 5, no. 4, pp. 326–328, 1972.

L. C. Yang and V. J. Menichelli, “Detonation of insensitive high explosives by a q-switched ruby laser,” Applied Physics Letters, vol. 19, no. 11, pp. 473–475, 1971.

R. J. Harrach, “Estimates on the ignition of high-explosives by laser pulses,” Journal of Applied Physics, vol. 47, no. 6, pp. 2473–2482, 1976.

M. W. Chen, S. You, K. S. Suslick, and D. D. Dlott, “Hot spot generation in energetic materials created by long-wavelength infrared radiation,” Applied Physics Letters, vol. 104, no. 6, Article ID 061907, 2014.

C.-L. Zhang, X.-D. Yuan, X. Xiang et al., “Rear-surface light intensification caused by a Hertzian-conical crack in 355-nm silica optics,” Chinese Physics B, vol. 21, no. 9, Article ID 094213, 2012.

P. Kosmas, Y. Wang, and C. Rappaport, “Three-dimensional FDTD model for GPR detection of objects buried in realistic dispersive soil,” in Detection and Remediation Technologies for Mines and Minelike Targets IV, vol. 4742 of Proceedings of SPIE, pp. 330–338, April 2002.

L. Li, X. Xiang, X. Zu, X. Jiang, X. Yuan, and W. Zheng, “The electric field modulation by three-dimensional crack on fused silica subsurface,” Optik, vol. 122, no. 16, pp. 1423–1425, 2011.

J.-R. Hua, L. Li, X. Xiang, and X.-T. Zu, “Three-dimensional numerical simulation of light field modulation in the vicinity of inclusions in silica subsurface,” Wuli Xuebao/Acta Physica Sinica, vol. 60, no. 4, Article ID 044206, 2011.

K. Yee, “Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media,” IEEE Transactions on Antennas and Propagation, vol. 14, no. 3, pp. 302–307, 1966.

J. W. Pyper, H. M. Buettner, C. J. Cerjan, J. S. Hallam, and R. J. King, “The measurement of bound and free moisture in organic materials by microwave methods,” in Proceedings of the International Symposium on Moisture and Humidity, Washington, DC, USA, 1985.

J.-P. Berenger, “A perfectly matched layer for the absorption of electromagnetic waves,” Journal of Computational Physics, vol. 114, no. 2, pp. 185–200, 1994.