The damage mechanisms that lead to failure in engineering alloys have been studied extensively, but converting this knowledge into constitutive models that are suitable for engineering-scale analysis remains a challenge. Evolution laws for continuum damage have been developed in the past and have proven effective but suffer from many non-physical assumptions that inhibit the overall accuracy of the model. Further, the assumptions inherent in these existing models prevent them from being applicable to a broad class of materials. At the same time, computational models of fine-scale damage mechanisms continue to advance making it tractable to generate large training data sets through computer simulation. Data-driven machine learning approaches can leverage these data sets to avoid making limiting assumptions, and instead produce models directly from the results of microstructural simulations and/or experiments. Many of these machine learning approaches are rapid and accurate, but they offer little to no insight into the underlying relationships among state variables being discovered. Conversely, genetic programming symbolic regression (GPSR) is a machine learning method that produces analytic expressions relating the state variables, allowing maximal insight and interpretability. To that end, we propose using GPSR as a data-driven method of obtaining microstructurally informed continuum damage models. Data is generated using microstructural simulations of damage evolution, parameterized over microstructural statistics (i.e., pore shape) and nominally applied deformations. Analytic expressions for damage evolution are obtained from the data using GPSR, and these expressions are then utilized within a continuum constitutive model. Overall, this approach is a promising method of automatically obtaining analytic relations describing constitutive phenomena in a material.
Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. We use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton's method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. We show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Furthermore, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.
Capacitance/inductance corrections for grid induced errors for a thin slot models are given for both one and four point testing on a rectangular grid for surface currents surrounding the slot. In addition a formula for translating from one equivalent radius to another is given for the thin-slot transmission line model. Additional formulas useful for this slot modeling are also given.
A method used to solve the problem of water waves on a sloping beach is applied to a thin conducting half plane described by a thin layer impedance boundary condition. The solution for the electric field behavior near the edge is obtained and a simple fit for this behavior is given. This field is used to determine the correction to the impedance per unit length of a conductor due to a sharp edge. The results are applied to the strip conductor. The final appendix also discusses the solution to the dual-sided (impedance surface & perfect conductor surface) half plane problem.
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.
New approaches to preventing and treating infections, particularly of the respiratory tract, are needed. One promising strategy is to reconfigure microbial communities (microbiomes) within the host to improve defense against pathogens. Probiotics and prebiotics for gastrointestinal (GI) infections offer a template for success. We sought to develop comparable countermeasures for respiratory infections. First, we characterized interactions between the airway microbiome and a biodefense-related respiratory pathogen (Burkholderia thailandensis; Bt), using a mouse model of infection. Then, we recovered microbiome constituents from the airway and assessed their ability to re-colonize the airway and protect against respiratory Bt infection. We found that microbiome constituents belonging to Bacillus and related genuses frequently displayed colonization and anti-Bt activity. Comparative growth requirement profiling of these Bacillus strains vs Bt enabled identification of candidate prebiotics. This work serves as proof of concept for airway probiotics, as well as a strong foundation for development of airway prebiotics.
The purpose and scope of the viga span tables project for Rachel Wood Consulting (RWC) is focused on producing tabulated beam span tables for three species of wood vigas commonly used in New Mexico to allow producers, designers and builders to incorporate vigas into their building designs in a prescriptive manner similar to the span tables for sawn lumber incorporated into the International Residential Code (IRC) or the International Log Builders Association (ILBA) publication. The information provided in this report and the associated viga span tables also attempts to address and clarify questions raised by RWC during their review of the 2018 Los Alamos National Laboratory (LANL) New Mexico Small Business Assistance (NMSBA) program report by August Mosimann pertaining to span lengths, loading, deflection calculations, and log grading certification prior to submitting the span tables to the Construction Industries Division (CID) of New Mexico.
Stinchfield, Georgia; Biegler, Lorenz T.; Eslick, John C.; Jacobson, Clas; Miller, David C.; Siirola, John D.; Zamarripa; Zhang, Chen; Zhang, Qi; Laird, Carl D.
This report examines the problem of magnetic penetration of a conductive layer, including nonlinear ferromagnetic layers, excited by an electric current filament. The electric current filament is, for example, a nearby wire excited by a lightning strike. The internal electric field and external magnetic field are determined. Numerical results are compared to various analytical approximations to help understand the physics involved in the penetration.
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.
This document presents the facility-recommended characterization of the neutron, prompt gamma-ray, and delayed gamma-ray radiation fields for the White Sands Missile Range (WSMR) Fast Burst Reactor, also known as molybdenum-alloy Godiva (Molly-G), at the 6-inch and the 24-inch irradiation locations. The neutron, prompt gamma-ray, and delayed gamma-ray energy spectra, uncertainties, and covariance matrices are presented. Code dependent recommended constants are given to facilitate the conversion of various dosimetry readings into radiation metrics desired by experimenters. Representative pulse operations are presented with conversion examples.
Recent efforts at Sandia such as DataSEA are creating search engines that enable analysts to query the institution’s massive archive of simulation and experiment data. The benefit of this work is that analysts will be able to retrieve all historical information about a system component that the institution has amassed over the years and make better-informed decisions in current work. As DataSEA gains momentum, it faces multiple technical challenges relating to capacity storage. From a raw capacity perspective, data producers will rapidly overwhelm the system with massive amounts of data. From an accessibility perspective, analysts will expect to be able to retrieve any portion of the bulk data, from any system on the enterprise network. Sandia’s Institutional Computing is mitigating storage problems at the enterprise level by procuring new capacity storage systems that can be accessed from anywhere on the enterprise network. These systems use the simple storage service, or S3, API for data transfers. While S3 uses objects instead of files, users can access it from their desktops or Sandia’s high-performance computing (HPC) platforms. S3 is particularly well suited for bulk storage in DataSEA, as datasets can be decomposed into object that can be referenced and retrieved individually, as needed by an analyst. In this report we describe our experiences working with S3 storage and provide information about how developers can leverage Sandia’s current systems. We present performance results from two sets of experiments. First, we measure S3 throughput when exchanging data between four different HPC platforms and two different enterprise S3 storage systems on the Sandia Restricted Network (SRN). Second, we measure the performance of S3 when communicating with a custom-built Ceph storage system that was constructed from HPC components. Overall, while S3 storage is significantly slower than traditional HPC storage, it provides significant accessibility benefits that will be valuable for archiving and exploiting historical data. There are multiple opportunities that arise from this work, including enhancing DataSEA to leverage S3 for bulk storage and adding native S3 support to Sandia’s IOSS library.
Some existing approaches to modelling the thermodynamics of moist air make approximations that break thermodynamic consistency, such that the resulting thermodynamics does not obey the first and second laws or has other inconsistencies. Recently, an approach to avoid such inconsistency has been suggested: the use of thermodynamic potentials in terms of their natural variables, from which all thermodynamic quantities and relationships (equations of state) are derived. In this article, we develop this approach for unapproximated moist-air thermodynamics and two widely used approximations: the constant- (Formula presented.) approximation and the dry heat capacities approximation. The (consistent) constant- (Formula presented.) approximation is particularly attractive because it leads to, with the appropriate choice of thermodynamic variable, adiabatic dynamics that depend only on total mass and are independent of the breakdown between water forms. Additionally, a wide variety of material from different sources in the literature on thermodynamics in atmospheric modelling is brought together. It is hoped that this article provides a comprehensive reference for the use of thermodynamic potentials in atmospheric modelling, especially for the three systems considered here.