Publications

Abstract not provided.

Abstract not provided.

This brief paper explores the development of scalable, nonlinear, fully-implicit solution methods for a stabilized unstructured finite element (FE) discretization of the 2D incompressible (reduced) resistive MHD system. The discussion considers the stabilized FE formulation in context of a fully-implicit time integration and direct-to-steady-state solution capability. The nonlinear solver strategy employs Newton-Krylov methods, which are preconditioned using fully-coupled algebraic multilevel (AMG) techniques and a new approximate block factorization (ABF) preconditioner. The intent of these preconditioners is to enable robust, scalable and efficient solution approaches for the large-scale sparse linear systems generated by the Newton linearization. We present results for the fully-coupled AMG preconditioner for two prototype problems, a low Lundquist number MHD Faraday conduction pump and moderately-high Lundquist number incompressible magnetic island coalescence problem. For the MHD pump results we explore the scaling of the fully-coupled AMG preconditioner for up to 4096 processors for problems with up to 64M unknowns on a CrayXT3/4. Using the island coalescence problem we explore the weak scaling of the AMG preconditioner and the influence of the Lundquist number on the iteration count. Finally we present some very recent results for the algorithmic scaling of the ABF preconditioner.

Abstract not provided.

Phase transformation between the ferroelectric (FE) and the antiferroelectric (AFE) phases in tin modified lead zirconate titanate (PSZT) ceramics can be influenced by pressure and electric field. Increasing the pressure has the tendency to favor the AFE phase while electric field favors the FE phase. In this study, these phase transformations are studied as functions of external pressure, temperature, and dc bias. The shifting of transformation temperature and the relative phase stability between FE and AFE with respect to these external parameters will be presented. Results will be compared to a pressure-induced depoling behavior (or FE-to-AFE phase transformation) for the PSZT ceramic. Fundamental issues relates to the relative phase stability will be discussed from the perspective of lattice dynamics theory.

Grid-tied PV energy smoothing was implemented by using a valve regulated lead-acid (VRLA) battery as a temporary energy storage device to both charge and discharge as required to smooth the inverter energy output from the PV array. Inverter output was controlled by the average solar irradiance over the previous 1h time interval. On a clear day the solar irradiance power curve is offset by about 1h, while on a variable cloudy day the inverter output power curve will be smoothed based on the average solar irradiance. Test results demonstrate that this smoothing algorithm works very well. Battery state of charge was more difficult to manage because of the variable system inefficiencies. Testing continued for 30-days and established consistent operational performance for extended periods of time under a wide variety of resource conditions. Both battery technologies from Exide (Absolyte) and East Penn (Advanced Valve Regulated Lead-Acid) proved to cycle well at a partial state of charge over the time interval tested.

Dielectric resonators are an effective means to realize isotropic, low-loss optical metamaterials. As proof of this concept, a cubic resonator is analytically designed and then tested in the long-wave infrared.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Dynamical systems theory provides a powerful framework for understanding the behavior of complex evolving systems. However applying these ideas to large-scale dynamical systems such as discretizations of multi-dimensional PDEs is challenging. Such systems can easily give rise to problems with billions of dynamical variables, requiring specialized numerical algorithms implemented on high performance computing architectures with thousands of processors. This talk will describe LOCA, the Library of Continuation Algorithms, a suite of scalable continuation and bifurcation tools optimized for these types of systems that is part of the Trilinos software collection. In particular, we will describe continuation and bifurcation analysis techniques designed for large-scale dynamical systems that are based on specialized parallel linear algebra methods for solving augmented linear systems. We will also discuss several other Trilinos tools providing nonlinear solvers (NOX), eigensolvers (Anasazi), iterative linear solvers (AztecOO and Belos), preconditioners (Ifpack, ML, Amesos) and parallel linear algebra data structures (Epetra and Tpetra) that LOCA can leverage for efficient and scalable analysis of large-scale dynamical systems.

Near Real Time Accountability (NRTA) of actinides at high precision in reprocessing plants has been a long sought-after goal in the safeguards community. Achieving this goal is hampered by the difficulty of making precision measurements in the reprocessing environment, equipment cost, and impact to plant operations. Thus the design of future reprocessing plants requires an optimization of different approaches. The Separations and Safeguards Performance Model, developed at Sandia National Laboratories, was used to evaluate a number of NRTA strategies in a UREX+ reprocessing plant. Strategies examined include the incorporation of additional actinide measurements of internal plant vessels, more use of process monitoring data, and the option of periodic draining of inventory to key tanks. Preliminary results show that the addition of measurement technologies can increase the overall measurement uncertainty due to additional error propagation, so care must be taken when designing an advanced system. Initial results also show that relying on a combination of different NRTA techniques will likely be the best option. The model provides a platform for integrating all the data. The modeling results for the different NRTA options under various material loss conditions will be presented.

Abstract not provided.

Simulating the detailed evolution of microstructure at the mesoscale is increasingly being addressed by a number of methods. Discrete element modeling and Potts kinetic Monte Carlo have achieved success in capturing different aspects of sintering well. Discrete element cannot treat the details of neck formation and other shape evolution, especially when considering particles of arbitrary shapes. Potts kMC treats the micorstructural evolution very well, but cannot incorporate complex stress states that form especially during differential sintering. A model that is capable of simulating microstructural evolution during sintering at the mesoscale and can incorporate differential stresses is being developed. This multi-physics model that can treat both interfacial energies and the inter-particle stresses will be introduced. It will be applied to simulate microstructural evolution while resolving individual particles and the stresses that develop between them due to local shrinkage. Results will be presented and the future development of this model will be discussed.

Abstract not provided.

The National Nuclear Security Administration (NNSA) has submitted an application to the Nuclear Regulatory Commission (NRC) for the air shipment of plutonium metal within the Plutonium Air Transportable (PAT-1) packaging. The PAT-1 packaging is currently authorized for the air transport of plutonium oxide in solid form only. The INMM presentation will provide a limited overview of the scope of the plutonium metal initiative and provide a status of the NNSA application to the NRC.

The Standing Committee on International Security of Radioactive and Nuclear Materials in the Nonproliferation and Arms Control Division conducted its fourth annual workshop in February 2010 on Reducing the Risk from Radioactive and Nuclear Materials. This workshop examined new technologies in real-time tracking of radioactive materials, new risks and policy issues in transportation security, the best practices and challenges found in addressing illicit radioactive materials trafficking, industry leadership in reducing proliferation risk, and verification of the Nuclear Nonproliferation Treaty, Article VI. Technology gaps, policy gaps, and prioritization for addressing the identified gaps were discussed. Participants included academia, policy makers, radioactive materials users, physical security and safeguards specialists, and vendors of radioactive sources and transportation services. This paper summarizes the results of this workshop with the recommendations and calls to action for the Institute of Nuclear Materials Management (INMM) membership community.

Abstract not provided.

The data in many disciplines such as social networks, web analysis, etc. is link-based, and the link structure can be exploited for many different data mining tasks. In this paper, we consider the problem of temporal link prediction: Given link data for time periods 1 through T, can we predict the links in time period T + 1? Specifically, we look at bipartite graphs changing over time and consider matrix- and tensor-based methods for predicting links. We present a weight-based method for collapsing multi-year data into a single matrix. We show how the well-known Katz method for link prediction can be extended to bipartite graphs and, moreover, approximated in a scalable way using a truncated singular value decomposition. Using a CANDECOMP/PARAFAC tensor decomposition of the data, we illustrate the usefulness of exploiting the natural three-dimensional structure of temporal link data. Through several numerical experiments, we demonstrate that both matrix- and tensor-based techniques are effective for temporal link prediction despite the inherent difficulty of the problem.

We present the bandwidth enhancement of an EAM monolithically integrated with two mutually injection-locked lasers. An improvement in the modulation efficiency and bandwidth are shown with mutual injection locking.

Abstract not provided.

Abstract not provided.

Doppler radars can distinguish targets from clutter if the target's velocity along the radar line of sight is beyond that of the clutter. Some targets of interest may have a Doppler shift similar to that of clutter. The nature of sea clutter is different in the clutter and exo-clutter regions. This behavior requires special consideration regarding where a radar can expect to find sea-clutter returns in Doppler space and what detection algorithms are most appropriate to help mitigate false alarms and increase probability of detection of a target. This paper studies the existing state-of-the-art in the understanding of Doppler characteristics of sea clutter and scattering from the ocean to better understand the design and performance choices of a radar in differentiating targets from clutter under prevailing sea conditions.

Abstract not provided.