We construct a family of embedded pairs for optimal explicit strong stability preserving Runge–Kutta methods of order 2≤p≤4 to be used to obtain numerical solution of spatially discretized hyperbolic PDEs. In this construction, the goals include non-defective property, large stability region, and small error values as defined in Dekker and Verwer (1984) and Kennedy et al. (2000). The new family of embedded pairs offer the ability for strong stability preserving (SSP) methods to adapt by varying the step-size. Through several numerical experiments, we assess the overall effectiveness in terms of work versus precision while also taking into consideration accuracy and stability.
To detect a specific radio-frequency (rf) magnetic field, rf optically pumped magnetometers (OPMs) require a static magnetic field to set the Larmor frequency of the atoms equal to the frequency of interest. However, unshielded and variable magnetic field environments (e.g., an rf OPM on a moving platform) pose a problem for rf OPM operation. Here, we demonstrate the use of a natural-abundance rubidium vapor to make a comagnetometer to address this challenge. Our implementation builds upon the simultaneous application of several OPM techniques within the same vapor cell. First, we use a modified implementation of an OPM variometer based on 87Rb to detect and actively cancel unwanted external fields at frequencies 60Hz using active feedback to a set of field control coils. We exploit this stabilized field environment to implement a high-sensitivity rf magnetometer using 85Rb. Using this approach, we demonstrate the ability to measure rf fields with a sensitivity of approximately 9fTHz-1/2 inside a magnetic shield in the presence of an applied field of approximately 20μT along three mutually orthogonal directions. This demonstration opens up a path toward completely unshielded operation of a high-sensitivity rf OPM.
Tang, Yanfei; Mclaughlan, John E.; Grest, Gary S.; Cheng, Shengfeng
A method of simulating the drying process of a soft matter solution with an implicit solvent model by moving the liquid-vapor interface is applied to various solution films and droplets. For a solution of a polymer and nanoparticles, we observe “polymer-on-top” stratification, similar to that found previously with an explicit solvent model. Furthermore, “polymer-on-top” is found even when the nanoparticle size is smaller than the radius of gyration of the polymer chains. For a suspension droplet of a bidisperse mixture of nanoparticles, we show that core-shell clusters of nanoparticles can be obtained via the “small-on-outside” stratification mechanism at fast evaporation rates. “Large-on-outside” stratification and uniform particle distribution are also observed when the evaporation rate is reduced. Polymeric particles with various morphologies, including Janus spheres, core-shell particles, and patchy particles, are produced from drying droplets of polymer solutions by combining fast evaporation with a controlled interaction between the polymers and the liquid-vapor interface. Our results validate the applicability of the moving interface method to a wide range of drying systems. The limitations of the method are pointed out and cautions are provided to potential practitioners on cases where the method might fail.
The V31 containment vessel was procured by the US Army Recovered Chemical Material Directorate (RCMD) as a third - generation EDS containment vessel. It is the fifth EDS vessel to be fabricated under Code Case 2564 of the 2019 ASME Boiler and Pressure Vessel Code, which provides rules for the design of impulsively loaded vessels. The explosive rating for the vessel, based on the code case, is twenty-four (24) pounds TNT - equivalent for up to 1092 detonations. This report documents the results of explosive tests that were performed on the vessel at Sandia National Laboratories in Albuquerque, New Mexico to qualify the vessel for field operations use. There were three design basis configurations for qualification testing. Qualification test (1) consisted of a simulated M55 rocket motor and warhead assembly of 24lbs of Composition C-4 (30 lb TNT equivalent). This test was considered the maximum load case, based on modeling and simulation methods performed by Sandia prior to the vessel design phase. Qualification test (2) consisted of a regular, right circular cylinder, unitary charge, located central to the vessel interior of 19.2 lb of Composition C-4 (24 lb TNT equivalent). Qualification test (3) consisted of a 12-pack of regular, right circular cylinders of 2 lb each, distributed evenly inside the vessel (totaling 19.2 lb of C-4, or 24 lb TNT equivalent). All vessel acceptance criteria were met.
Conversion cathode materials are gaining interest for secondary batteries due to their high theoretical energy and power density. However, practical application as a secondary battery material is currently limited by practical issues such as poor cyclability. To better understand these materials, we have developed a pseudo-two-dimensional model for conversion cathodes. We apply this model to FeS2 – a material that undergoes intercalation followed by conversion during discharge. The model is derived from the half-cell Doyle–Fuller–Newman model with additional loss terms added to reflect the converted shell resistance as the reaction progresses. We also account for polydisperse active material particles by incorporating a variable active surface area and effective particle radius. Using the model, we show that the leading loss mechanisms for FeS2 are associated with solid-state diffusion and electrical transport limitations through the converted shell material. The polydisperse simulations are also compared to a monodisperse system, and we show that polydispersity has very little effect on the intercalation behavior yet leads to capacity loss during the conversion reaction. We provide the code as an open-source Python Battery Mathematical Modeling (PyBaMM) model that can be used to identify performance limitations for other conversion cathode materials.
Radiation-induced segregation is a phenomenon commonly observed in many alloys which consists of the redistribution of elements (solute or interstitial impurities) under irradiation. The onset and development of radiation-induced segregation can only occur when a sufficient flux of defects is sustained and defect sinks are present. Irradiation dose, dose rate, and particle types all affect defect flux. In this work, we employ a phase-field model to examine the effects of dose, dose rate, and type of incident particles on radiation-induced segregation behavior in a model binary alloy. The phase-field model takes into account the formation and evolution of point defects as well as defect clusters, the diffusion and clustering of alloy species, the presence of additional extrinsic defect sinks in the form of dislocations, and two different methods of radiation-damage insertion, which are intended to simulate either light-ion/electron irradiation via Frenkel pairs or heavy-ion irradiation in the form of cascades. Our results show a dose-rate and particle-type dependence on the amount of solute segregation. We show that the material systems exposed to higher dose rates are less subjected to solute segregation at equivalent doses. We also show that such dose-rate-dependence behavior is due to a delay of the incubation dose at which radiation-induced segregation effectively starts. Particle type and the presence of dislocations can accentuate this behavior. Our model predictions correlate with many experimental observations made over the years on radiation-induced segregation providing credence to the simulation results. The methodology presented in this study allows for a first-order prediction of the dose rate at which proxy irradiation experiments could be performed to approximate radiation-induced segregation behaviors seen in targeted irradiation conditions.
Efficiently embedding and/or integrating mechanistic information with data-driven models is essential if it is desired to simultaneously take advantage of both engineering principles and data-science. The opportunity for hybridization occurs in many scenarios, such as the development of a faster model of an accurate high-fidelity computer model; the correction of a mechanistic model that does not fully-capture the physical phenomena of the system; or the integration of a data-driven component approximating an unknown correlation within a mechanistic model. At the same time, different techniques have been proposed and applied in different literatures to achieve this hybridization, such as hybrid modeling, physics-informed Machine Learning (ML) and model calibration. In this paper we review the methods, challenges, applications and algorithms of these three research areas and discuss them in the context of the different hybridization scenarios. Moreover, we provide a comprehensive comparison of the hybridization techniques with respect to their differences and similarities, as well as advantages and limitations and future perspectives. Finally, we apply and illustrate hybrid modeling, physics-informed ML and model calibration via a chemical reactor case study.
This work details the development of a concentrating solar power (CSP) and thermal (CST) library archive. This work included digitization of one-of-a-kind documents that could be degraded or destroyed over time. Sandia National Laboratories (SNL) National Solar Thermal Test Facility (NSTTF) and Sandia's Technical Library departments collaborated to establish and maintain the first and only digital collection in the world of Concentrating Solar Power (CSP) related historical documents. These date back to the CSP program inception here at Sandia in the early 1970's thru to the present.
Significant variety is observed in spherical crystal x-ray imager (SCXI) data for the stagnated fuel–liner system created in Magnetized Liner Inertial Fusion (MagLIF) experiments conducted at the Sandia National Laboratories Z-facility. As a result, image analysis tasks involving, e.g., region-of-interest selection (i.e. segmentation), background subtraction and image registration have generally required tedious manual treatment leading to increased risk of irreproducibility, lack of uncertainty quantification and smaller-scale studies using only a fraction of available data. We present a convolutional neural network (CNN)-based pipeline to automate much of the image processing workflow. This tool enabled batch preprocessing of an ensemble of Nscans = 139 SCXI images across Nexp = 67 different experiments for subsequent study. The pipeline begins by segmenting images into the stagnated fuel and background using a CNN trained on synthetic images generated from a geometric model of a physical three-dimensional plasma. The resulting segmentation allows for a rules-based registration. Our approach flexibly handles rarely occurring artifacts through minimal user input and avoids the need for extensive hand labelling and augmentation of our experimental dataset that would be needed to train an end-to-end pipeline. We also fit background pixels using low-degree polynomials, and perform a statistical assessment of the background and noise properties over the entire image database. Our results provide a guide for choices made in statistical inference models using stagnation image data and can be applied in the generation of synthetic datasets with realistic choices of noise statistics and background models used for machine learning tasks in MagLIF data analysis. We anticipate that the method may be readily extended to automate other MagLIF stagnation imaging applications.