Publications

Abstract not provided.

This study examined a concept for improving the nation's ability to detect unauthorized transport of radiological and nuclear material that involves detecting not the radiological or nuclear material itself, but rather the anomalous presence of materials, such as lead and tungsten, typically used to shield such illicit materials from detection. Specifically, the project studied a potentially inexpensive, non-intrusive, and fast method of detecting anomalous shielding: use of magnetic induction sensors to measure currents induced in conductors by a driver coil. The first phase of this work, documented in this report, identified the signatures of materials of interest, assessed the usefulness of measured signatures for discrimination between materials, and then benchmarked models and calculations. The next task involved assessing the strength of these signatures when the materials of interest were enclosed within other materials, such as aluminum and steel, of various thicknesses. A final task mimicked a real-world scenario by testing the ability to detect material signatures of interest at a scale of a vehicle or cargo vessel against a cluttered backdrop, with and without the shielding of another material, such as aluminum or steel. This final task raised doubts regarding whether this detection scheme will be useful under real-world conditions.

In this report we discuss a new maritime surveillance and detection concept based on Raman scattering of water molecules. Using a scanning lidar that detects Raman scattered photons from water, the absence or change of signal indicates the presence of a non-water object. With sufficient spatial resolution a negative two dimensional image of the object can be generated by the scanning lidar. Because Raman scattering is an inelastic process with a relatively large wavelength shift for water, this concept completely avoids the problematic elastic scattering for objects at or very close to the water surface. Elastic scattering makes it difficult to discriminate between water and dark objects at or near the water surface especially when automated detection is required. It is also difficult to deal with elastic scattering from the bottom surface for shallow waters. The maximum detection depth for this concept is limited by the attenuation of the excitation and return Raman light in water. If excitation in the UV is used, fluorescence can be used for discrimination between organic and non-organic objects. Range gating can be used for this concept for detection of objects below a specified depth. In this report we develop a lidar model for this concept to estimate the number of detected Raman photons for variable lidar parameters and depths in the presence of the solar background. We also report on the results of proof-of-concept measurements using the Sandia Ares lidar with excitation at 355 nm. The measurements show good agreement with the lidar model predictions. The detected number of photons for typical lidar parameter shows the concept is viable and applicable to a variety of day and nighttime detection scenarios. This concept has many potential applications including near-surface mine detection, swimmer detection for security purposes, wide area search, as well as other civilian applications.

Abstract not provided.

Suppression of the ambient gamma background radiation by traffic structure and cargo is a well-understood and studied effect for deployed radiation portal monitors (RPM). For effective analysis of measured RPM profiles with the objective of inferring the spatial characteristics of radiation sources, it is important to account for the effects of background suppression. In this report we analyze background suppression for a test dataset from vehicle RPMs at a sample port and estimate the distributions of suppression amplitudes and shapes. Cluster analysis of standardized and normalized profiles is used to obtain the dominant suppression shapes in the sample field data. We determine that a large fraction of non-alarm RPM occupancies are represented by a small number of suppression shapes. This fraction increases when the signal-to-noise ratio of an occupancy profile is improved by the addition of signals for multiple RPM detectors located at the same height. The calculated suppression shapes from RPM data can be used along with source models in the process of spatial profile analysis both in the field or offline. This background suppression analysis is an important step in improving the effectiveness of the RPM profile analysis methodology which is currently being investigated and may lead to methods that reduce the number of secondary inspections as well as to decision support tools that aid operators in evaluating RPM data that do not contain spectral information.

Abstract not provided.

Abstract not provided.

Conventional full spectrum gamma spectroscopic analysis has the objective of quantitative identification of all the radionuclides present in a measurement. For low-energy resolution detectors such as NaI, when photopeaks alone are not sufficient for complete isotopic identification, such analysis requires template spectra for all the radionuclides present in the measurement. When many radionuclides are present it is difficult to make the correct identification and this process often requires many attempts to obtain a statistically valid solution by highly skilled spectroscopists. A previous report investigated using the targeted principal component analysis method (TPCA) for detection of embedded sources for RPM applications. This method uses spatial/temporal information from multiple spectral measurements to test the hypothesis of the presence of a target spectrum of interest in these measurements without the need to identify all the other radionuclides present. The previous analysis showed that the TPCA method has significant potential for automated detection of target radionuclides of interest, but did not include the effects of shielding. This report complements the previous analysis by including the effects of spectral distortion due to shielding effects for the same problem of detection of embedded sources. Two examples, one with one target radionuclide and the other with two, show that the TPCA method can successfully detect shielded targets in the presence of many other radionuclides. The shielding parameters are determined as part of the optimization process using interpolation of library spectra that are defined on a 2D grid of atomic numbers and areal densities.

Conventional full spectrum gamma spectroscopic analysis has the objective of quantitative identification of all the isotopes present in a measurement. For low energy resolution detectors, when photopeaks alone are not sufficient for complete isotopic identification, such analysis requires template spectra for all the isotopes present in the measurement. When many isotopes are present it is difficult to make the correct identification and this process often requires many trial solutions by highly skilled spectroscopists. This report investigates the potential of a new analysis method which uses spatial/temporal information from multiple low energy resolution measurements to test the hypothesis of the presence of a target spectrum of interest in these measurements without the need to identify all the other isotopes present. This method is referred to as targeted principal component analysis (TPCA). For radiation portal monitor applications, multiple measurements of gamma spectra are taken at equally spaced time increments as a vehicle passes through the portal and the TPCA method is directly applicable to this type of measurement. In this report we describe the method and investigate its application to the problem of detection of a radioactive localized source that is embedded in a distributed source in the presence of an ambient background. Examples using simulated spectral measurements indicate that this method works very well and has the potential for automated analysis for RPM applications. This method is also expected to work well for isotopic detection in the presence of spectrally and spatially varying backgrounds as a result of vehicle-induced background suppression. Further work is needed to include effects of shielding, to understand detection limits, setting of thresholds, and to estimate false positive probability.

Abstract not provided.

In a multiyear research agreement with Tenix Investments Pty. Ltd., Sandia has been developing field deployable technologies for detection of biotoxins in water supply systems. The unattended water sensor or UWS employs microfluidic chip based gel electrophoresis for monitoring biological analytes in a small integrated sensor platform. This instrument collects, prepares, and analyzes water samples in an automated manner. Sample analysis is done using the {mu}ChemLab{trademark} analysis module. This report uses analysis results of two datasets collected using the UWS to estimate performance of the device. The first dataset is made up of samples containing ricin at varying concentrations and is used for assessing instrument response and detection probability. The second dataset is comprised of analyses of water samples collected at a water utility which are used to assess the false positive probability. The analyses of the two sets are used to estimate the Receiver Operating Characteristic or ROC curves for the device at one set of operational and detection algorithm parameters. For these parameters and based on a statistical estimate, the ricin probability of detection is about 0.9 at a concentration of 5 nM for a false positive probability of 1 x 10{sup -6}.

We present a rapid method for the identification of viruses using microfluidic chip gel electrophoresis (CGE) of high-copy number proteins to generate unique protein profiles. Viral proteins are solubilized by heating at 95°C in borate buffer containing detergent (5 min), then labeled with fluorescamine dye (10 s), and analyzed using the μChemLab CGE system (5 min). Analyses of closely related T2 and T4 bacteriophage demonstrate sufficient assay sensitivity and peak resolution to distinguish the two phage. CGE analyses of four additional viruses - MS2 bacteriophage, Epstein - Barr, respiratory syncytial, and vaccinia viruses - demonstrate reproducible and visually distinct protein profiles. To evaluate the suitability of the method for unique identification of viruses, we employed a Bayesian classification approach. Using a subset of 126 replicate electropherograms of the six viruses and phage for training purposes, successful classification with non-training data was 66/69 or 95% with no false positives. The classification method is based on a single attribute (elution time), although other attributes such as peak width, peak amplitude, or peak shape could be incorporated and may improve performance further. The encouraging results suggest a rapid and simple way to identify viruses without requiring specialty reagents such as PCR probes and antibodies. © 2008 American Chemical Society.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report addresses the problem of selection of lidar parameters, namely wavelengths for absorption lidar and excitation fluorescence pairs for fluorescence lidar, for optimal detection of species. Orthogonal spectra and cross sections are used as mathematical representations which provide a quantitative measure of species distinguishability in mixtures. Using these quantities, a simple expression for the absolute error in calculated species concentration is derived and optimization is accomplished by variation of lidar parameters to minimize this error. It is shown that the optimum number of wavelengths for detection of a species using absorption lidar (excitation fluorescence pairs for fluorescence lidar) is the same as the number of species in the mixture. Each species present in the mixture has its own set of optimum wavelengths. There is usually some overlap in these sets. The optimization method is applied to two examples, one using absorption and the other using fluorescence lidar, for analyzing mixtures of four organic compounds. The effect of atmospheric attenuation is included in the optimization process. Although the number of optimum wavelengths might be small, it is essential to do large numbers of measurements at these wavelengths in order to maximize canceling of statistical errors.

Performance of an electromagnetic induction launcher is considered for three types of armatures. These are: solid, 1-element wound and 16-element wound aluminum armatures. The one element wound armature has uniform current density throughout and thus can withstand field reversal (working against embedded armature flux) and still maintain low temperature. Slingshot simulations were performed for several configurations. Best performance was obtained for a single element wound armature with two field reversals. For a 60 kg projectile, 10.5 cm coil inner radius and 5.5 cm coil build, the velocity after 50 meters of launcher length (670 stages) exceeded 3.5 km/sec with an overall efficiency of about 45%. For the same parameters the solid and 16-element wound armatures reach a velocity of about 3.3 km/sec after 800 stages (60 meters of launcher length) but without field reversal. A velocity of 3.5 km/sec is possible after 60 meters of launcher length with the 16-element wound armature with one field reversal, but the temperature is close to the melting temperature of aluminum. In all simulations with a solid armature, melting of some of the surface material occurs. However, it is shown that most of the melting occurs after contribution has been made to the forward going pressure, that is, melting does not affect the electrical performance of the launcher. The effect of coil firing time jitter on launcher performance is also considered and is found to be very small for realistic perturbations. For ± 2 μ–secs random jitter, the reduction in the final velocity for a 60 meter launcher with a solid armature is less than 0.1% and the increase in temperature is only 2%. This holds for all types of armatures. © 1995 IEEE

The transverse motion of a projectile in an electromagnetic induction launcher is considered. The equations of motion for translation and rotation are derived assuming a rigid projectile and a flyway restoring force per unit length that is proportional to the local displacement. Transverse forces and torques due to energized coils are derived for displaced or tilted projectile elements based on a first order perturbation method. The resulting equations of motion for a rigid projectile composed of multiple elements in a multi-coil launcher are analyzed as a coupled oscillator system of equations and a simple stability condition is derived. The equations of motion are incorporated into the 2-D Slingshot code and numerical solutions for the transverse motion are obtained. For the 20 meter navy launcher parameters we find that stability is achieved with a flyway spring constant of k {approx} 1{times} 10{sup 8} N/m{sup 2}. For k {approx} 1.5 {times} 10{sup 8} N/m{sup 2} and sample coil misalignment modeled as a sine wave of I mm amplitude at wavelengths of one or two meters, the projectile displacement grows to a maximum of 4 mm. This growth is due to resonance between the natural frequency of the Projectile transverse motion and the coil displacement wavelength. This resonance does not persist because of the changing axial velocity. Random coil displacement is also found to cause roughly the same projectile displacement. For the maximum displacement a rough estimate of the transverse pressure is 50 bars.

The Recirculating Linear Accelerator (RLA) utilizes the Ion Focused Regime (IFR) of beam transport plus a ramped bending field to guide the beam around the curved sections. Several issues of beam transport are considered. Beam transverse perturbations that could result in growth of the ion hose instability are analyzed. It is found that transverse kicks due to bending field errors, energy mismatches and fringe fields are the most important. The scaling of these perturbations with beam and channel parameters is derived. The effect of ramping of the bending field on the preformed plasma channel is then considered. For RLA experimental parameters the effect is found to be very small. For high energies however, in addition to axial heating, it is found that ramping the field causes compression of the plasma channel along the radius of curvature. This compression results in a quasi-equilibrium plasma electron temperature along the field lines which leads to collisionless transport towards the walls. The analysis of compression is done in an approximate way using a single particle picture and the channel expansion is analyzed using an envelope solution which gives a simple expression for the expansion time. This solution is then verified by Buckshot simulations. For a bending field of 2 kG ramped in 2 {mu}-secs and an argon channel (RLA parameters) we estimate that the channel radius doubling time (along field lines) is of the order of 0.5 {mu}-secs. Finally the effect of electron impact ionization due to axially heated electrons by the action of the inductive field is estimated. It is found that in Argon gas the electron avalanche time could be as low as 0.5 {mu}-sec which is smaller than the field ramp time.

The equations of motion of beam and channel particles are analyzed in the ion focused regime. Using the paraxial approximation and assuming only transverse electro-magnetostatic interactions between beam and channel particles for the equations of motion (the same equations solved in the BUCKSHOT code) are written in non-dimensional form and scaling relations are derived for propagation parameters, magnetic erosion, and evaporation. These relations are very useful in doing parameter studies with a limited number of computer simulations.

Phase mixing of transverse oscillations changes the nature of the ion hose instability from an absolute to a convective instability. The stronger the phase mixing, the faster an electron beam reaches equilibrium with the guiding ion channel. This is important for long distance propagation of relativistic electron beams where it is desired that transverse oscillations phase mix within a few betatron wavelengths of injection and subsequently an equilibrium is reached with no further beam emittance growth. In the linear regime phase mixing is well understood and results in asymptotic decay of transverse oscillations as 1/Z{sup 2} for a Gaussian beam and channel system, Z being the axial distance measured in betatron wavelengths. In the nonlinear regime (which is likely mode of propagation for long pulse beams) results of the spread mass model indicate that phase mixing is considerably weaker than in the regime. In this paper we consider this problem of phase mixing in the nonlinear regime. Results of the spread mass model will be shown along with a simple analysis of phase mixing for multiple oscillator models. Particle simulations also indicate that phase mixing is weaker in nonlinear regime than in the linear regime. These results will also be shown. 3 refs., 4 figs.

Charge exchange results in the transfer of momentum from an ion to a neutral atom, thus it is equivalent to a slowing down mechanism for ions. This might be helpful in IFR propagation of relativistic electron beams. The effect of charge exchange on the ion hose instability and channel expansion after beam passage is investigated in this report. For parameters of interest to current experiments, it is found that the effect on growth of the ion hose instability is insignificant, but there is significant reduction in the energy of the expanding ions after passage of the electron beam. 5 refs., 5 figs.

A new gridless electrostatic field solver which utilities Fourier decomposition in the azimuthal coordinate has been developed and tested. The scaling with the number of simulation particles is N log N. This algorithm has been implemented in the BUCKSHOT code, which originally used a direct summation algorithm with N{sup 2} scaling. The Fourier decomposition in the new algorithm is done about the center of mass of each species, thus nonlinear ion hose physics is included in the m = O mode. Higher order modes describe non-axisymmetric profile changes. The breakeven point between the new solver and the direct summation algorithm is about N = 64 particles per species when up to m = 2 Fourier modes are kept. For a typical ion hose simulation with 256 particles per species the new solver is faster by a factor of about 2.7. 8 refs., 11 figs.

The Buckshot code is a gridless particle simulation code which is used extensively at Sandia to study relativistic electron beam propagation in the ion-focused regime. The present version of the code evaluates the force on a particle by summing over all the other particles, thus the execution time is proportional to N{sup 2} where N is the number of simulation particles. A new gridless field solver has been developed with time scaling of N log N and the breakeven point between the old and new code is around N = 64 on the Cray X-MP computer. For N = 1000, the new solver is about nine times faster than the old version. The new solver is based on a solution of the two-dimensional Poisson equation which if Fourier decomposed in the azimuthal direction and the solution of the radial equation is represented by integrals over the charge density. These integrals are then replaced by sums over the simulation particles which are assumed to be point particles. The near-field singularity is removed by the Fourier decomposition so long as the number of Fourier modes is much less than the number of simulation particles. The algorithm is written in such a way that the field due to a given species is Fourier decomposed with respect to the center of mass of that species, thus it is possible to study linear and nonlinear ion hose physics with a very small number of azimuthal modes. Typically M = 2 is found to be sufficient for most IFR problems. The old and new solvers have been compared and the agreement is excellent. 2 refs., 3 figs.