Publications

This report examines the localization of high frequency electromagnetic fields in general three-dimensional convex walled cavities along periodic paths between opposing sides of the cavity. The report examines the three-dimensional case where the mirrors at the end of the orbit have two different radii of curvature. The cases where these orbits lead to unstable localized modes are known as scars.

Use of insensitive high explosives (IHEs) has significantly improved ammunition safety because of their remarkable insensitivity to violent cook-off, shock and impact. Triamino-trinitrobenzene (TATB) is the IHE used in many modern munitions. Previously, lightning simulations in different test configurations have shown that the required detonation threshold for standard density TATB at ambient and elevated temperatures (250 C) has a sufficient margin over the shock caused by an arc from the most severe lightning. In this paper, the Braginskii model with Lee-More channel conductivity prescription is used to demonstrate how electrical arcs from lightning could cause detonation in TATB. The steep rise and slow decay in typical lightning pulse are used in demonstrating that the shock pressure from an electrical arc, after reaching the peak, falls off faster than the inverse of the arc radius. For detonation to occur, two necessary detonation conditions must be met: the Pop-Plot criterion and minimum spot size requirement. The relevant Pop-Plot for TATB at 250 C was converted into an empirical detonation criterion, which is applicable to explosives subject to shocks of variable pressure. The arc cross-section was required to meet the minimum detonation spot size reported in the literature. One caveat is that when the shock pressure exceeds the detonation pressure the Pop-Plot may not be applicable, and the minimum spot size requirement may be smaller.

This paper provides an overview of the electromagnetic frequency domain simulation capabilities of the Electromagnetic Theory department at Sandia National Laboratories via a description of two of its codes. EIGER is a Method of Moments code for electromagnetic simulations, but it only runs on traditional CPUs, not on new architectures. Gemma is in development to replace EIGER and will run on many architectures, including CPUs, GPUs, and MICs, by leveraging the Kokkos library.

Abstract not provided.

Abstract not provided.

In this paper we investigate full-wave simulations of realistic implementations of multifunctional nanoantenna enabled detectors (NEDs). We focus on a 2x2 pixelated array structure that supports two wavelengths of operation. We design each resonating structure independently using full-wave simulations with periodic boundary conditions mimicking the whole infinite array. We then construct a supercell made of a 2x2 pixelated array with periodic boundary conditions mimicking the full NED; in this case, however, each pixel comprises 10-20 antennas per side. In this way, the cross-talk between contiguous pixels is accounted for in our simulations. We observe that, even though there are finite extent effects, the pixels work as designed, each responding at the respective wavelength of operation. This allows us to stress that realistic simulations of multifunctional NEDs need to be performed to verify the design functionality by taking into account finite extent and cross-talk effects.

Abstract not provided.

A block base sparse approximate inverse preconditioner for the electric field integral equations is documented and tested. It utilized the Kokkos library for performance portability and shows superior performance when compared to a direct method, 36x faster for a 112K DOF problem. Furthermore, due to the abstractions available in the Kokkos library it allows one to migrate from CPU to GPU in a trivial way.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The goal of this paper is to investigate full-wave simulations of realistic implementations of multifunctional nanoantenna enabled detectors (NEDs). We realize a 2×2 pixelated array structure that supports two wavelengths of operation. After designing each resonating structure independently using full-wave simulations with periodic boundary conditions mimicking the whole infinite array, we construct a supercell made of a 2×2 pixelated array with periodic boundary conditions mimicking the full NED. In the NED, each pixel comprises 10-20 nanoantennas. Our simulations account for the cross-talk between contiguous pixels. We observe that, even though there are finite extent effects, the pixels work as designed, each responding at the respective wavelength of operation. We want to stress that realistic simulations of multifunctional NEDs need to be performed to verify the design functionality by taking into account finite extent and cross-talk effects.

The growth of a cylindrical s park discharge channel in water and Lexan is studied using a series of one - dimensional simulations with the finite - element radiation - magnetohydrodynamics code ALEGRA. Computed solutions are analyzed in order to characterize the rate of growth and dynamics of the spark c hannels during the rising - current phase of the drive pulse. The current ramp rate is varied between 0.2 and 3.0 kA/ns, and values of the mechanical coupling coefficient K p are extracted for each case. The simulations predict spark channel expansion veloc ities primarily in the range of 2000 to 3500 m/s, channel pressures primarily in the range 10 - 40 GPa, and K p values primarily between 1.1 and 1.4. When Lexan is preheated, slightly larger expansion velocities and smaller K p values are predicted , but the o verall behavior is unchanged.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report examines bounds on the penetrant power through ports of entry into a conductive cavity. We first replace the cavity by a load and consider the maximum power transfer properties of an antenna or an aperture. We consider how limitations on the load quality factor place limits on received power. For general frequency ranges we model the backing region by means of a uniformly distributed matched load along a slot aperture and adjust its value for maximum power transfer. This result is derived in closed form using a transmission line model for the aperture. This result illustrates the reduction in received power for low frequencies with finitely conducting wall materials. At high frequencies it approaches the receiving cross section of a linear array having the slot length dimension. Next we examine a slot aperture in a conducting rectangular enclosure and determine how the cavity wall losses and resulting quality factor limit the penetrant power. Detailed simulations and experimental measurements are compared with each other and with the bounding results to assess the accuracy of the bounds. These comparisons also indicate limitations on the accuracy of the models due to perturbing influences in construction, such as bolted joints.

Streamers are a type of ionization wave occurring during the early time phase of a gas discharge. They are typically launched when the evolving space charge of an electron avalanche reaches a certain critical level, beyond which the fi eld of the space charge itself is su ffi cient to drive further evolution of the ionization process. One of the most common ways to model streamers is known as a 1.5D model where the fi eld of a uniformly charged set of discs of chosen radius is evaluated along the cylinder axis. This fi eld drives a one-dimensional kinetic ionization process, which results in the nonlinear evolution of the streamer. This model is e ffi cient, but has the drawback of fi xing the radius and requiring it as an input parameter. Previously, we tried to extend t he 1.5D model to include evolution of its radius by developing a two-step process of axial and radial exp ansion but we encountered stability issues with the model that we thought could have been due to decoupling the two steps. In this report we introduce a new formulation of a streamer model that includes radial expansion. The goal is to take radial moments of the starting axisymmetric fl uid equations and thereby include the radial evolution of the streamer naturally and self-consistently from the beginning. We fi rst develop the fl uid model moments without electron attachment. We review the calculation of the electric fi elds required for the model and investigate approximations to improve computational e ffi ciency. We discuss the code implementation of the model and fi nally, we add attachment to allow the treatment of electronegative gases. Intentionally Left Blank

The di ff usion through shells consisting of either a single conducting or double conducting layers are examined. Exterior drives resulting from Electromagn etic Radiation (EMR), Electromagnetic Pulse (EMP), nearby (indirect) lightning, and DC (low frequency) magnetic fi eldsareused. Boththeinterior fi eld and the induced voltage from a maximally oriented and sized single turn loop are estimated. It is shown that the loop voltage with the empty cavity bounds the case where the center region is excluded by a conducting object. The cases of interior magnetic and electric fi elds from an exterior magnetic drive and the interior electric fi eld from an exterior electric drive are both solved; the magnetic interior fi eldfromanexterior magnetic drive is the only case that results in a nonzero low frequency penetration. Intentionally Left Blank

Abstract not provided.

Abstract not provided.

This report examines the localization of high frequency electromagnetic fi elds in three-dimensional axisymmetric cavities along periodic paths between opposing sides of the cavity. The cases where these orbits lead to unstable localized modes are known as scars. This report treats both the case where the opposing sides, or mirrors, are convex, where there are no interior foci, and the case where they are concave, leading to interior foci. The scalar problem is treated fi rst but the approximations required to treat the vector fi eld components are also examined. Particular att ention is focused on the normalization through the electromagnetic energy theorem. Both projections of the fi eld along the scarred orbit as well as point statistics are examined. Statistical comparisons are m ade with a numerical calculation of the scars run with an axisymmetric simulation. This axisymmetric cas eformstheoppositeextreme(wherethetwomirror radii at each end of the ray orbit are equal) from the two -dimensional solution examined previously (where one mirror radius is vastly di ff erent from the other). The enhancement of the fi eldontheorbitaxiscanbe larger here than in the two-dimensional case. Intentionally Left Blank

When emitters of electromagnetic energy are operated in the vicinity of sensitive components, the electric field at the component location must be kept below a certain level in order to prevent the component from being damaged, or in the case of electro-explosive devices, initiating. The V-Curve is a convenient way to set the electric field limit because it requires minimal information about the problem configuration. In this report we will discuss the basis for the V-Curve. We also consider deviations from the original V-Curve resulting from inductive versus capacitive antennas, increases in directivity gain for long antennas, decreases in input impedance when operating in a bounded region, and mismatches dictated by transmission line losses. In addition, we consider mitigating effects resulting from limited antenna sizes.