Publications

Abstract not provided.

Abstract not provided.

This report outlines the development of load-mitigating feedback control for wave energy converters. A simple, self-tuning multi-objective controller is demonstrated in simulation for a 3-DOF (surge, heave, pitch) point absorber. In previous work, the proposed control architecture has been shown to be effective in experiment for a variety of device archetypes for the single objective of the maximization of electrical power capture: here this architecture is extended to reduce device loading as well. In particular, PTO actuation forces and the minimization of fatigue damage (determined from the sum of wave-exerted and PTO forces) are considered as additional objectives for the self-tuning controller. This controller is demonstrated for two similar, but distinct systems: one described by the identified linear models from physical testing of the WaveBot device, and another based upon a WEC-Sim simulation that expands upon boundary element method data from the WaveBot device. In both cases, because the power surface is consistently fairly flat in the vicinity of control parameters that maximize power capture in contrasting sea-states, it is found to be generally possible to mitigate either fatigue damage or PTO load. However, PTO load is found to conflict with fatigue damage in some sea-states, limiting the efficacy of control objectives that attempt to mitigate both simultaneously. Additionally, coupling between the surge and pitch DOFs also limits the extent to which fatigue damage can be mitigated for both DOFs in some sea-states. Because control objectives can be considered a function of the sea-state (e.g., load mitigation may not be a concern until the sea is sufficiently large) a simple transition strategy is proposed and demonstrated. This transition strategy is found to be effective with some caveats: firstly, it cannot circumvent the aforementioned objective contradictions. Secondly, this objective transition is too slow to act as a system constraint, and objective thresholds must thus be considered quite conservatively. Improvement of the adjustment strategy is demonstrated through the addition of an integral term. Selection of well-performing transition parameters can be a function of sea-state. While a simple selection procedure is proposed, it is non-optimal, and a more robust selection procedure is suggested for future work.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

National Technology & Engineering Solutions of Sandia, LLC (NTESS) has recently amended the NTESS Retirement Income Plan (Pension Plan). The updated Summary Plan Description (SPD) for the Pension Plan effective January 1, 2022 is provided.

Abstract not provided.

Abstract not provided.

This report details a method to estimate the energy content of various types of seismic body waves. The method is based on the strain energy of an elastic wavefield and Hooke’s Law. We present a detailed derivation of a set of equations that explicitly partition the seismic strain energy into two parts: one for compressional (P) waves and one for shear (S) waves. We posit that the ratio of these two quantities can be used to determine the relative contribution of seismic P and S waves, possibly as a method to discriminate between earthquakes and buried explosions. We demonstrate the efficacy of our method by using it to compute the strain energy of synthetic seismograms with differing source characteristics. Specifically, we find that explosion-generated seismograms contain a preponderance of P wave strain energy when compared to earthquake-generated synthetic seismograms. Conversely, earthquake-generated synthetic seismograms contain a much greater degree of S wave strain energy when compared to explosion-generated seismograms.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report details how to successfully use the Fairfield Nodal ZLand seismic instruments to collect data, including preparation steps prior to deploying the instruments, how to record data during a field campaign, and how to retrieve recorded data from the instruments after their deployment. This guide will walk through each step for the novice user, as well as provide a checklist of critical steps for the advanced user to ensure successful, efficient field campaigns and seismic data collection. Currently, use of the seismic nodal instruments is highly limited due to the detailed nature and prior knowledge required to successfully set up, use, and retrieve data from these instruments. With this guide, all interested users will have the knowledge required to perform a seismic deployment and collect data with the Fairfield Nodal instruments.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Savannah River Site plans to reprocess defense spent nuclear fuel currently stored in their L-Basin via the Accelerated Basin Deinventory (ABD) Program. The previous plan for the L-Basin spent nuclear fuel was to dispose of it directly in the federal repository without reprocessing. Implementing the ABD Program will result in final disposal of approximately 900 fewer canisters of defense spent nuclear fuel and the production of approximately 521 more canisters of vitrified high-level waste glass with some specific differences from the planned high-level waste glass. Because the 235U in the L-Basin spent nuclear fuel is not intended to be recovered, the fissile mass loading of the vitrified high-level glass waste form to be produced must be increased above the current value of 897 g/m³ to a maximum of 2,500 g/m³. Therefore, implementing the ABD Program would produce a variant of high-level waste glass—the ABD glass—that needs to be evaluated for future repository licensing, which includes both preclosure safety and postclosure performance. This report describes the approach to and summarizes the results of an evaluation of the potential effects of implementing the ABD Program at the Savannah River Site on the technical basis for future repository licensing for a generic repository that is similar to Yucca Mountain and for one that is fully generic. This evaluation includes the effects on preclosure safety analyses and postclosure performance assessment for both repository settings. The license application for the proposed Yucca Mountain repository (DOE 2008), which is serving as a framework for this evaluation, concluded that the proposed Yucca Mountain repository would meet all applicable regulatory requirements. The evaluation documented in this report found that implementing the ABD Program is not expected to change that conclusion for a generic repository similar to Yucca Mountain or for a generic repository with respect to the preclosure safety analyses. With respect to the postclosure performance of a generic repository, no concerns were identified.

Filtration, pressure drop and quantitative fit of N95 respirators were robust to several decontamination methods including vaporous hydrogen peroxide, wet heat, bleach, and ultraviolet light. Bleach may not have penetrated the hydrophobic outer layers of the N95 respirator. Isopropyl alcohol and detergent both severely degraded the electrostatic charge of the electret filtration layer. First data in N95 respirators that the loss of filtration efficiency was directly correlated with loss of surface potential on the filtration layer. The pressure drop was unchanged, so loss of filtration efficacy would not be apparent during a user seal check. Mechanical straps degrade with repeated mechanical cycling during extended use. Decontamination did not appear to degrade the elastic straps. Significant loss of strap elasticity would be apparent during a user negative pressure seal check.

Here, we explore the dimensionality of the U.S. Department of Agriculture’s household food security survey module among households with children. Using a novel methodological approach to measuring food security, we find that there is multidimensionality in the module for households with children that is associated with the overall household, adult, and child dimensions of food security. Additional analyses suggest official estimates of food security among households with children are robust to this multidimensionality. However, we also find that accounting for the multidimensionality of food security among these households provides new insights into the correlates of food security at the household, adult, and child levels of measurement.

A myriad of phenomena in materials science and chemistry rely on quantum-level simulations of the electronic structure in matter. While moving to larger length and time scales has been a pressing issue for decades, such large-scale electronic structure calculations are still challenging despite modern software approaches and advances in high-performance computing. The silver lining in this regard is the use of machine learning to accelerate electronic structure calculations – this line of research has recently gained growing attention. The grand challenge therein is finding a suitable machine-learning model during a process called hyperparameter optimization. This, however, causes a massive computational overhead in addition to that of data generation. We accelerate the construction of machine-learning surrogate models by roughly two orders of magnitude by circumventing excessive training during the hyperparameter optimization phase. We demonstrate our workflow for Kohn-Sham density functional theory, the most popular computational method in materials science and chemistry.

The objective of this project was to develop a novel capability to generate synthetic data sets for the purpose of training Machine Learning (ML) algorithms for the detection of malicious activities on satellite systems. The approach experimented with was to a) generate sparse data sets using emulation modeling and b) enlarge the sparse data using Generative Adversarial Networks (GANs). We based our emulation modeling on the Open Source NASA Operational Simulator for Small Satellites (NOS3) developed by the Katherine Johnson Independent Verification and Validation (IV&V) program in West Virginia. Significant new capabilities on NOS3 had to be developed for our data set generation needs. To expand these data sets for the purpose of training ML, we experimented with a) Extreme Learning Machines (ELMs) and b) Wasserstein-GANs (WGAN-GP).

The core function of many neural network algorithms is the dot product, or vector matrix multiply (VMM) operation. Crossbar arrays utilizing resistive memory elements can reduce computational energy in neural algorithms by up to five orders of magnitude compared to conventional CPUs. Moving data between a processor, SRAM, and DRAM dominates energy consumption. By utilizing analog operations to reduce data movement, resistive memory crossbars can enable processing of large amounts of data at lower energy than conventional memory architectures.

For isolated white dwarf (WD) stars, fits to their observed spectra provide the most precise estimates of their effective temperatures and surface gravities. Even so, recent studies have shown that systematic offsets exist between such spectroscopic parameter determinations and those based on broadband photometry. These large discrepancies (10% in Teff, 0.1 M⊙ in mass) provide scientific motivation for reconsidering the atomic physics employed in the model atmospheres of these stars. Recent simulation work of ours suggests that the most important remaining uncertainties in simulation-based calculations of line shapes are the treatment of 1) the electric field distribution and 2) the occupation probability (OP) prescription. We review the work that has been done in these areas and outline possible avenues for progress.

The role of a solid surface for initiating gas-phase reactions is still not well understood. The hydrogen atom (H) is an important intermediate in gas-phase ethane dehydrogenation and is known to interact with surface sites on catalysts. However, direct measurements of H near catalytic surfaces have not yet been reported. Here, we present the first H measurements by laser-induced fluorescence in the gas-phase above catalytic and noncatalytic surfaces. Measurements at temperatures up to 700 °C show H concentrations to be at the highest above inert quartz surfaces compared to stainless steel and a platinum-based catalyst. Additionally, H concentrations above the catalyst decreased rapidly with time on stream. These newly obtained observations are consistent with the recently reported differences in bulk ethane dehydrogenation reactivity of these materials, suggesting H may be a good reporter for dehydrogenation activity.

In global atmospheric modeling, the differences between nonhydrostatic (NH) and hydrostatic (H) dynamical cores are negligible in dry simulations when grid spacing is larger than 10 km. However, recent studies suggest that those differences can be significant at far coarser resolution when moisture is included. To better understand how NH and H differences manifest in global fields, we perform and analyze an ensemble of 28 and 13 km seasonal simulations with the NH and H dynamical cores in the Energy Exascale Earth System Model global atmosphere model, where the differences between H and NH configurations are minimized. A set of idealized rising bubble experiments is also conducted to further investigate the differences. Although NH and H differences are not significant in global statistics and zonal averages, significant differences in precipitation amount and patterns are observed in parts of the tropics. The most prominent differences emerge near India and the Western Pacific in the boreal summer, and the central-southern Indian Ocean and Pacific in the boreal winter. Tropical differences influence surrounding regions through modification of the regional circulation and can propagate to the extratropics, leading to significant temperature and geopotential differences over the middle to high latitudes. While the dry bubble experiments show negligible deviation between H and NH dynamics until grid spacing is below 6.25 km, precipitation amount and vertical velocity are different in the moist case even at 25 km resolution.

Abstract not provided.

This presentation provides details regarding integral experiments at Sandia National Laboratory for fiscal year 2021. The experiments discussed are as follows: IER 230: Characterize the Thermal Capabilities of the 7uPCX; IER 304: Temperature Dependent Critical Benchmarks; IER 305: Critical Experiments with UO₂ Rods and Molybdenum Foils; IER 306: Critical Experiments with UO₂ Rods and Rhodium Foils ; IER 441: Epithermal HEX Lattices with SNL 7uPCX Fuel for Testing Nuclear Data; IER 452: Inversion Point of the Isothermal Reactivity Coefficient; and IER 523: Critical Experiments with ACRR UO₂-BeO Fuel.

We introduce a robust verification tool for computational codes, which we call Stochastic Robust Extrapolation based Error Quantification (StREEQ). Unlike the prevalent Grid Convergence Index (GCI) [1] method, our approach is suitable for both stochastic and deterministic computational codes and is generalizable to any number of discretization variables. Building on ideas introduced in the Robust Verification [2] approach, we estimate the converged solution and orders of convergence with uncertainty using multiple fits of a discretization error model. In contrast to Robust Verification, we perform these fits to many bootstrap samples yielding a larger set of predictions with smoother statistics. Here, bootstrap resampling is performed on the lack-of-fit errors for deterministic code responses, and directly on the noisy data set for stochastic responses. This approach lends a degree of robustness to the overall results, capable of yielding precise verification results for sufficiently resolved data sets, and appropriately expanding the uncertainty when the data set does not support a precise result. For stochastic responses, a credibility assessment is also performed to give the analyst an indication of the trustworthiness of the results. This approach is suitable for both code and solution verification, and is particularly useful for solution verification of high-consequence simulations.

Sandia National Labs has access to unused ACRR fuel, which is unique in its enrichment 35% and material composition BeO. ACRR fuel is available in quantities well above what is needed for experiments. Two experiment concepts have been investigated: UO₂BeO fuel elements and pellets with 7uPCX fuel. The worth of UO₂BeO is large enough to be well above the anticipated experiment uncertainties.

This work presents a new multiscale method for coupling the 3D Maxwell's equations to the 1D telegrapher's equations. While Maxwell's equations are appropriate for modeling complex electromagnetics in arbitrary-geometry domains, simulation cost for many applications (e.g. pulsed power) can be dramatically reduced by representing less complex transmission line regions of the domain with a 1D model. By assuming a transverse electromagnetic (TEM) ansatz for the solution in a transmission line region, we reduce the Maxwell's equations to the telegrapher's equations. We propose a self-consistent finite element formulation of the fully coupled system that uses boundary integrals to couple between the 3D and 1D domains and supports arbitrary unstructured 3D meshes. Additionally, by using a Lagrange multiplier to enforce continuity at the coupling interface, we allow for an absorbing boundary condition to also be applied to non-TEM modes on this boundary. We demonstrate that this feature reduces non-physical reflection and ringing of non-TEM modes off of the coupling boundary. By employing implicit time integration, we ensure a stable coupling, and we introduce an efficient method for solving the resulting linear systems. We demonstrate the accuracy of the new method on two verification problems, a transient O-wave in a rectilinear prism and a steady-state problem in a coaxial geometry, and show the efficiency and weak scalability of our implementation on a cold test of the Z-machine MITL and post-hole convolute.

This presentation discusses activities related to the Nuclear Criticality Safety Program (NCSP) at Sandia National Laboratory in fiscal year 2021. This includes NCSP funding, integral experiment requests, integral experiment spending, highlights, and COVID-19 impacts.

In situ analysis of surfaces during high-flux plasma exposure represents a long-standing challenge in the study of plasma-material interactions. While post-mortem microscopy can provide a detailed picture of structural and compositional changes, in situ techniques can capture the dynamic evolution of the surface. In this study, we demonstrate how spectroscopic ellipsometry can be applied to the real-time characterization of W nanostructure (also known as "fuzz") growth during exposure to low temperature, high-flux He plasmas. Strikingly, over a wide range of sample temperatures and helium fluences, the measured ellipsometric parameters (ψ, Δ) collapse onto a single curve that can be directly correlated with surface morphologies characterized by ex situ helium ion microscopy. The initial variation in the (ψ, Δ) parameters appears to be governed by small changes in surface roughness (<50 nm) produced by helium bubble nucleation and growth, followed by the emergence of 50 nm diameter W tendrils. This basic behavior appears to be reproducible over a wide parameter space, indicating that the spectroscopic ellipsometry may be of general practical use as a diagnostic to study surface morphologies produced by high-flux He implantation in refractory metals. An advantage of the methods outlined here is that they are applicable at low incident ion energies, even below the sputtering threshold. As an example of this application, we apply in situ ellipsometry to examine how W fuzz growth is affected both by varying ion energy and the temperature of the surface.

Here, we utilize electrically detected magnetic resonance (EDMR) measurements to compare high-field stressed, and gamma irradiated Si/SiO₂ metal–oxide–silicon (MOS) structures. We utilize spin-dependent recombination (SDR) EDMR detected using the Fitzgerald and Grove dc I-V approach to compare the effects of high-field electrical stressing and gamma irradiation on defect formation at and near the Si/SiO₂ interface. As anticipated, both greatly increase the concentration of P_b centers (silicon dangling bonds at the interface) densities. The irradiation also generated a significant increase in the dc I-V EDMR response of E' centers (oxygen vacancies in the SiO₂ films), whereas the generation of an E' EDMR response in high-field stressing is much weaker than in the gamma irradiation case. These results likely suggest a difference in their physical distribution resulting from radiation damage and high electric field stressing.

The reactivity of carbonyl oxides has previously been shown to exhibit strong conformer and substituent dependencies. Through a combination of synchrotron-multiplexed photoionization mass spectrometry experiments (298 K and 4 Torr) and high-level theory [CCSD(T)-F12/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ with an added CCSDT(Q) correction], we explore the conformer dependence of the reaction of acetaldehyde oxide (CH3CHOO) with dimethylamine (DMA). The experimental data support the theoretically predicted 1,2-insertion mechanism and the formation of an amine-functionalized hydroperoxide reaction product. Tunable-vacuum ultraviolet photoionization probing of anti- or anti- + syn-CH3CHOO reveals a strong conformer dependence of the title reaction. The rate coefficient of DMA with anti-CH3CHOO is predicted to exceed that for the reaction with syn-CH3CHOO by a factor of ∼34,000, which is attributed to submerged barrier (syn) versus barrierless (anti) mechanisms for energetically downhill reactions.

CPU/GPU heterogeneous compute platforms are an ubiquitous element in computing and a programming model specified for this heterogeneous computing model is important for both performance and programmability. A programming model that exposes the shared, unified, address space between the heterogeneous units is a necessary step in this direction as it removes the burden of explicit data movement from the programmer while maintaining performance. GPU vendors, such as AMD and NVIDIA, have released software-managed runtimes that can provide programmers the illusion of unified CPU and GPU memory by automatically migrating data in and out of the GPU memory. However, this runtime support is not included in GPGPU-Sim, a commonly used framework that models the features of a modern graphics processor that are relevant to non-graphics applications. UVM Smart was developed, which extended GPGPU-Sim 3.x to in- corporate the modeling of on-demand pageing and data migration through the runtime. This report discusses the integration of UVM Smart and GPGPU-Sim 4.0 and the modifications to improve simulation performance and accuracy.

Enhancing the efficiency of second-harmonic generation using all-dielectric metasurfaces to date has mostly focused on electromagnetic engineering of optical modes in the meta-atom. Further advances in nonlinear conversion efficiencies can be gained by engineering the material nonlinearities at the nanoscale, however this cannot be achieved using conventional materials. Semiconductor heterostructures that support resonant nonlinearities using quantum engineered intersubband transitions can provide this new degree of freedom. By simultaneously optimizing the heterostructures and meta-atoms, we experimentally realize an all-dielectric polaritonic metasurface with a maximum second-harmonic generation power conversion factor of 0.5 mW/W2 and power conversion efficiencies of 0.015% at nominal pump intensities of 11 kW/cm2. These conversion efficiencies are higher than the record values reported to date in all-dielectric nonlinear metasurfaces but with 3 orders of magnitude lower pump power. Our results therefore open a new direction for designing efficient nonlinear all-dielectric metasurfaces for new classical and quantum light sources.

A combination of electrodeposition and thermal reduction methods have been utilized for the synthesis of ligand-free FeNiCo alloy nanoparticles through a high-entropy oxide intermediate. These phases are of great interest to the electrocatalysis community, especially when formed by a sustainable chemistry method. This is successfully achieved by first forming a complex five element amorphous FeNiCoCrMn high-entropy oxide (HEO) phase via electrodeposition from a nanodroplet emulsion solution of the metal salt reactants. The amorphous oxide phase is then thermally treated and reduced at 570-600 °C to form the crystalline FeNiCo alloy with a separate CrMnOx cophase. The FeNiCo alloy is fully characterized by scanning transmission electron microscopy and energy-dispersive X-ray spectroscopy elemental analysis and is identified as a face-centered cubic crystal with the lattice constant a = 3.52 Å. The unoptimized, ligand-free FeNiCo NPs activity toward the oxygen evolution reaction is evaluated in alkaline solution and found to have an ∼185 mV more cathodic onset potential than the Pt metal. Beyond being able to synthesize highly crystalline, ligand-free FeNiCo nanoparticles, the demonstrated and relatively simple two-step process is ideal for the synthesis of tailor-made nanoparticles where the desired composition is not easily achieved with classical solution-based chemistries.

The protection systems (circuit breakers, relays, reclosers, and fuses) of the electric grid are the primary component responding to resilience events, ranging from common storms to extreme events. The protective equipment must detect and operate very quickly, generally <0.25 seconds, to remove faults in the system before the system goes unstable or additional equipment is damaged. The burden on protection systems is increasing as the complexity of the grid increases; renewable energy resources, particularly inverter-based resources (IBR) and increasing electrification all contribute to a more complex grid landscape for protection devices. In addition, there are increasing threats from natural disasters, aging infrastructure, and manmade attacks that can cause faults and disturbances in the electric grid. The challenge for the application of AI into power system protection is that events are rare and unpredictable. In order to improve the resiliency of the electric grid, AI has to be able to learn from very little data. During an extreme disaster, it may not be important that the perfect, most optimal action is taken, but AI must be guaranteed to always respond by moving the grid toward a more stable state during unseen events.

There are several engineering applications in which the assumptions of homogenization and scale separation may be violated, in particular, for metallic structures constructed through additive manufacturing. Instead of resorting to direct numerical simulation of the macroscale system with an embedded fine scale, an alternative approach is to use an approximate macroscale constitutive model, but then estimate the model-form error using a posteriori error estimation techniques and subsequently adapt the macroscale model to reduce the error for a given boundary value problem and quantity of interest. Here, we investigate this approach to multiscale analysis in solids with unseparated scales using the example of an additively manufactured metallic structure consisting of a polycrystalline microstructure that is neither periodic nor statistically homogeneous. As a first step to the general nonlinear case, we focus here on linear elasticity in which each grain within the polycrystal is linear elastic but anisotropic.

Electrically detected magnetic resonance and near-zero-field magnetoresistance measurements were used to study atomic-scale traps generated during high-field gate stressing in Si/SiO2 MOSFETs. The defects observed are almost certainly important to time-dependent dielectric breakdown. The measurements were made with spin-dependent recombination current involving defects at and near the Si/SiO2 boundary. The interface traps observed are Pb0 and Pb1 centers, which are silicon dangling bond defects. The ratio of Pb0/Pb1 is dependent on the gate stressing polarity. Electrically detected magnetic resonance measurements also reveal generation of E′ oxide defects near the Si/SiO2 interface. Near-zero-field magnetoresistance measurements made throughout stressing reveal that the local hyperfine environment of the interface traps changes with stressing time; these changes are almost certainly due to the redistribution of hydrogen near the interface.

Graph partitioning has emerged as an area of interest due to its use in various applications in computational research. One way to partition a graph is to solve for the eigenvectors of the corresponding graph Laplacian matrix. This project focuses on the eigensolver LOBPCG and the evaluation of a new preconditioner: Randomized Cholesky Factorization (rchol). This proconditioner was tested for its speed and accuracy against other well-known preconditioners for the method. After experiments were run on several known test matrices, rchol appears to be a better preconditioner for structured matrices. This research was sponsored by National Nuclear Security Administration Minority Serving Institutions Internship Program (NNSA-MSIIP) and completed at host facility Sandia National Laboratories. As such, after discussion of the research project itself, this report contains a brief reflection on experience gained as a result of participating in the NNSA-MSIIP.

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

We present an overview of the magneto-inertial fusion (MIF) concept MagLIF (Magnetized Liner Inertial Fusion) pursued at Sandia National Laboratories and review some of the most prominent results since the initial experiments in 2013. In MagLIF, a centimeter-scale beryllium tube or "liner" is filled with a fusion fuel, axially pre-magnetized, laser pre-heated, and finally imploded using up to 20 MA from the Z machine. All of these elements are necessary to generate a thermonuclear plasma: laser preheating raises the initial temperature of the fuel, the electrical current implodes the liner and quasi-adiabatically compresses the fuel via the Lorentz force, and the axial magnetic field limits thermal conduction from the hot plasma to the cold liner walls during the implosion. MagLIF is the first MIF concept to demonstrate fusion relevant temperatures, significant fusion production (>10^13 primary DD neutron yield), and magnetic trapping of charged fusion particles. On a 60 MA next-generation pulsed-power machine, two-dimensional simulations suggest that MagLIF has the potential to generate multi-MJ yields with significant self-heating, a long-term goal of the US Stockpile Stewardship Program. At currents exceeding 65 MA, the high gains required for fusion energy could be achievable.

Experiments were designed and conducted to investigate the impact that geometric cavities have on the transfer of energy from an embedded explosion to the surface of the physical domain. The experimental domains were fabricated as 3-inch polymer cubes, with varying cavity geometries centered in the cubes. The energy transfer, represented as a shock wave, was generated by the detonation of an exploding bridgewire at the center of the cavity. The shock propagation was tracked by schlieren imaging through the optically accessible polymer. The magnitude of energy transferred to the surface was recorded by an array of pressure sensors. A minimum of five experimental runs were conducted for each cavity geometry and statistical results were developed and compared. Results demonstrated the decoupling effect that geometric cavities produce on the energy field at the surface.

Computational tools to study thermodynamic properties of magnetic materials have, until recently, been limited to phenomenological modeling or to small domain sizes limiting our mechanistic understanding of thermal transport in ferromagnets. Herein, we study the interplay of phonon and magnetic spin contributions to the thermal conductivity in a-iron utilizing non-equilibrium molecular dynamics simulations. It was observed that the magnetic spin contribution to the total thermal conductivity exceeds lattice transport for temperatures up to two-thirds of the Curie temperature after which only strongly coupled magnon-phonon modes become active heat carriers. Characterizations of the phonon and magnon spectra give a detailed insight into the coupling between these heat carriers, and the temperature sensitivity of these coupled systems. Comparisons to both experiments and ab initio data support our inferred electronic thermal conductivity, supporting the coupled molecular dynamics/spin dynamics framework as a viable method to extend the predictive capability for magnetic material properties.

Spectral line-shape models are an important part of understanding high-energy-density (HED) plasmas. Models are needed for calculating opacity of materials and can serve as diagnostics for astrophysical and laboratory plasmas. However, much of the literature on line shapes is directed toward specialists. This perspective makes it difficult for non-specialists to enter the field. We have two broad goals with this topical review. First, we aim to give information so that others in HED physics may better understand the current field. This first goal may help guide future experiments to test different aspects of the theory. Second, we provide an introduction for those who might be interested in line-shape theory, and enough materials to be able to navigate the field and the literature. We give a high-level overview of line broadening process, as well as dive into the formalism, available methods, and approximations.

Electromagnetic (EM) methods are among the original techniques for subsurface characterization in exploration geophysics because of their particular sensitivity to the earth electrical conductivity, a physical property of rocks distinct yet complementary to density, magnetization, and strength. However, this unique ability also makes them sensitive to metallic artifacts - infrastructure such as pipes, cables, and other forms of cultural clutter - the EM footprint of which often far exceeds their diminutive stature when compared to that of bulk rock itself. In the hunt for buried treasure or unexploded ordnance, this is an advantage; in the long-term monitoring of mature oil fields after decades of production, it is quite troublesome indeed. Here we consider the latter through the lens of an evolving energy industry landscape in which the traditional methods of EM characterization for the exploration geophysicist are applied toward emergent problems in well-casing integrity, carbon capture and storage, and overall situational awareness in the oil field. We introduce case studies from these exemplars, showing how signals from metallic artifacts can dominate those from the target itself and impose significant burdens on the requisite simulation complexity. We also show how recent advances in numerical methods mitigate the computational explosivity of infrastructure modeling, providing feasible and real-time analysis tools for the desktop geophysicist. Lastly, we demonstrate through comparison of field data and simulation results that incorporation of infrastructure into the analysis of such geophysical data is, in a growing number of cases, a requisite but now manageable step.