Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Implementation of Detailed Polyethylene Pyrolysis Kinetics into CFD Simulations using Machine Learning

Municipal solid waste (MSW) and waste plastics have received significant attention due to the issues of waste generation and storage, as well as their potential as an energy resource. High-density polyethylene (HDPE) makes up a large portion of plastic waste and has been the subject of several conversion studies. However, the mechanisms associated with converting HDPE through pyrolysis and gasification are extensive and complex making them difficult to implement into high-fidelity computational fluid dynamic (CFD) simulations. For this project, a primary pyrolysis mechanism containing 42 unique species and 737 heterogeneous reactions was used to generate kinetic data over a range of operating conditions. A machine learning (ML) model was developed to replicate the results of the detailed pyrolysis mechanism while significantly increasing the computational efficiency. A deep operator network (DeepONet) architecture was adopted to train the model using time steps relevant to CFD simulations. The ML used physics-based loss functions to ensure mass conservation. The ML model has been deployed in simple MFiX CFD simulations, single particle, and an experimental drop tube reactor, and has shown promising performance compared to the original scheme.

Houston, Ross↗
DOE OSTI2021

Skyshine Calculations for a Large Spent Nuclear Fuel Storage Facility with SCALE 6.2.3

The SCALE code system developed at Oak Ridge National Laboratory includes state-of-the-art capabilities for radiation source term and radiation transport simulations that can be used in numerous applications, including dose rate analyses of complex consolidated interim storage facilities (CISFs). A licensed CISF could be used to store tens of thousands of tonnes of spent nuclear fuel discharged from commercial power reactors using various cask and storage pad designs. A CISF design must comply with the regulatory requirements provided in 10 CFR Part 72, including requirements related to annual dose limits applicable to real individuals located beyond the area controlled by the licensee. Therefore, calculating a dose to the public is a necessary part of the licensing process for the construction of a CISF. These calculations are very challenging because of the complexity of the CISF design and the low magnitude of dose rate at large distances from the facility. This paper describes detailed far-field dose rate calculations performed for a proposed CISF using MAVRIC, the Monte Carlo radiation shielding sequence in SCALE 6.2.3, with automated variance reduction based on discrete ordinates calculations. The method presented in this paper uses a detailed Monte Carlo radiation transport simulation in one step from source to dose rate. A series of independent simulations was made using the complete site geometry (all casks present), but with only one cask containing radiation sources to obtain the dose rate maps produced by each storage cask. The CISF dose rate map was obtained by adding the dose rate maps produced by the independent individual cask simulations. Ample volumes of air and soil extending beyond the location of interest for dose rate calculation were included in the calculation model to properly simulate important radiation attenuation and scattering events that affect far-field dose rates. Furthermore, a comprehensive sensitivity study is included in this paper to illustrate the importance of selecting appropriate air volume, mass density, and composition for CISF skyshine dose rate calculations. Dry soil and soil containing water were analyzed to determine their effects on groundshine radiation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Hardware-Accelerated Ray Tracing of CAD-Based Geometry for Monte Carlo Radiation Transport

Monte Carlo radiation transport (MCRT) methods have been used to simulate radiation environments for many decades by tracking individual particles through a model to accumulate statistical information. MCRT geometry is historically formed using the constructive solid geometry (CSG). Recently, significant work has been performed to support simulations using computer-aided design (CAD)-based tessellated surfaces to support highly complex geometries. Ray tracing acceleration data structures from the rendering and visualization community are applied to accelerate particle tracking in CAD-based models. Despite these efforts, CSG representations provide the superior performance in surface intersection operations during particle flight. Concurrently, pseudo Monte Carlo methods have become prevalent in rendering applications to support more realistic models for scattering media, motivating innovations that are advantageous for MCRT simulations. Finally, the authors’ work extends these innovations by employing Intel’s Embree ray tracing kernel within a geometry toolkit for Monte Carlo to improve the simulation performance using CAD-based models by factors of 1.5 to 2.

42 ENGINEERING↗
DOE OSTI2025

Structure-Based Design of Small-Molecule Inhibitors of Human Interleukin-6

Human Interleukin-6 (hIL-6) is a pro inflammatory cytokine that binds to its receptor, IL-6Rα followed by binding to gp130 and subsequent dimerization to form a hexamer signaling complex. As a critical inflammation mediator, hIL-6 is associated with a diverse range of diseases and monoclonal antibodies in clinical use that either target IL-6Rα or hIL-6 to inhibit signaling. Here, we perform high-throughput structure-based computational screening using ensemble docking for small-molecule antagonists for which the target conformations were taken from 600 ns long molecular dynamics simulations of the apo protein. Prior knowledge of the contact sites from binary complex studies and experimental work was incorporated into the docking studies. The top 20 scoring ligands from the in silico studies after post analysis were subjected to in vitro functional assays. Among these compounds, the ligand with the second-highest calculated binding affinity experimentally showed an ~84% inhibitory effect on IL6-induced STAT3 reporter activity at 10 μM concentration. This finding may pave the way for designing small-molecule inhibitors of hIL-6 of therapeutic significance.

Human Interleukin-6↗
DOE OSTI2025

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗
DOE OSTI2026

Stratospheric aerosol forcing for CMIP7 – Part 1: optical properties for pre-industrial, historical, and scenario simulations

Stratospheric aerosols, most of which originate from explosive volcanic sulfur emissions into the stratosphere, are a key natural driver of climate variability. They are thus a forcing provided by the Coupled Model Intercomparison Project (CMIP) Climate Forcings Task Team to climate modelling groups participating in phase 7 of CMIP. For the historical period, we provide two datasets covering 1750–2023: (i) a volcanic upper tropospheric-stratospheric sulfur emission dataset, documented in a companion paper; and (ii) a stratospheric sulfate aerosol optical property dataset, documented in this paper. For the satellite era (from 1979 onwards), stratospheric aerosol optical properties are derived from the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) dataset. For the pre-satellite era (1750–1978), optical properties are derived from our volcanic SO 2 emission dataset using a new version of the reduced-complexity volcanic aerosol model Easy Volcanic Aerosol (Height) (EVA_H v2). A background, non-volcanic stratospheric aerosol climatology is derived from the 1998–2001 period with a trend over 1850–1978 accounting for increasing anthropogenic aerosols. A monthly stratospheric aerosol climatology is derived from the 1850–2021 average for both pre-industrial and scenario (future) simulations, with a 9-year ramp over 2022–2030 for scenario simulations to ensure a smooth transition from the historical period. CMIP7 simulations should use version 2.2.1 of our pre-industrial and historical datasets, and version 2.2.2 of our scenario dataset. Our methodology to produce historical aerosol optical properties significantly differs from CMIP6 for the pre-satellite era, and the resulting forcings in turn largely differ. In particular, the CMIP6 dataset was mostly based on the sparse and uncertain pyrheliometer record, which resulted in strongly underrepresented emissions from small-to-moderate magnitude eruptions. The resulting bias is addressed in CMIP7, which is entirely emission-derived in the pre-satellite era and uses more recent ice-core-based volcanic sulfur emission inventories than CMIP6. Our approach results in an overall larger volcanic aerosol forcing for CMIP7, with the 1850–2014 mean mid-visible global mean stratospheric aerosol optical depth (SAOD) in CMIP7 (0.0138) being 29 % higher than in CMIP6 (0.0107). The pre-industrial mean of the same variable is 26 % higher in CMIP7 (0.0135, derived from the historical 1850–2021) than CMIP6 (0.0107, derived from the historical 1850–2014 mean). Using a reduced-complexity climate model, we simulate a global mean surface temperature that is 0.07 °C colder for 1850–1900 when using the CMIP7 dataset instead of CMIP6, whereas 2000–2014 is 0.03 °C warmer in CMIP7. Our dataset also exhibits lower forcing for 1960–1980, resulting in temperatures 0.06 °C warmer when averaged over 1960–1990, a period for which CMIP6 climate models exhibit a cold bias. Given the large uncertainties characterizing the dataset, in particular for the pre-satellite era, we advise against treating the CMIP7 or CMIP6 dataset as uniquely superior for any specific year and highlight the need for further evaluation. We conclude the study by discussing sources of uncertainty for the dataset, future research avenues to improve it, as well as requirements to operationalize the production of the dataset, i.e. extend it and update it on an annual basis instead of every 5–7 years following CMIP cycles.

Aubry, Thomas J. [Univ. of Exeter, Devon (United K↗
DOE OSTI2021

Numerical Simulation of Container Breach and Airborne Release of Solids Due to Mechanical Insults

Throughout U.S. Department of Energy (DOE) complexes, safety engineers employ the five-factor formula to calculate the source term (ST) that includes parameters of airborne release fraction (ARF), respirable fraction (RF) and damage ratio (DR). Limited experimental data on fragmentation of solids, such as ceramic pellets (i.e., PuO2), and container breach due to mechanical insults (i.e., drop and forklift impact), can be supplemented by modeling and simulation using high fidelity computational tools to estimate these parameters. This paper presents the use of Sandia National Laboratories' SIERRA solid mechanics (SM) finite element code to investigate the behavior of the widely utilized waste container (such as 7A Drum) subject to a range of free fall impact and puncture scenarios. The resulting behavior of the container is assessed, and the estimates are presented for bounding DRs from calculated breach areas for the various accident conditions considered. This paper also describes a novel multiscale constitutive model recently implemented in SIERRA/SM that simulates the fracture of brittle materials such as PuO 2 and determines ARF during the fracture process. Furthermore, comparisons are made between model predictions and simple bench-top experiments.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Simulating Real Atmospheric Boundary Layers at Gray-Zone Resolutions: How Do Currently Available Turbulence Parameterizations Perform?

Recent computational and modeling advances have led a diverse modeling community to experiment with atmospheric boundary layer (ABL) simulations at subkilometer horizontal scales. Accurately parameterizing turbulence at these scales is a complex problem. The modeling solutions proposed to date are still in the development phase and remain largely unvalidated. This work assesses the performance of methods currently available in the Weather Research and Forecasting (WRF) model to represent ABL turbulence at a gray-zone grid spacing of 333 m. We consider three one-dimensional boundary layer parameterizations (MYNN, YSU and Shin-Hong) and coarse large-eddy simulations (LES). The reference dataset consists of five real-case simulations performed with WRF-LES nested down to 25 m. Results reveal that users should refrain from coarse LES and favor the scale-aware, Shin-Hong parameterization over traditional one-dimensional schemes. Overall, the spread in model performance is large for the cellular convection regime corresponding to the majority of our cases, with coarse LES overestimating turbulent energy across scales and YSU underestimating it and failing to reproduce its horizontal structure. Despite yielding the best results, the Shin-Hong scheme overestimates the effect of grid dependence on turbulent transport, highlighting the outstanding need for improved solutions to seamlessly parameterize turbulence across scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Characterization of Fuel-to-Coolant Heat Transfer During Reactivity-Initiated Accidents Using Tightly Coupled Thermal Hydraulics and Fuel Thermomechanics

The reactivity-initiated accident (RIA) is a complex scenario with several tightly interacting physical phenomena. Accurately predicting fuel behavior during these transients is difficult due to limitations in the modeling of fuel-to-coolant heat transfer. Common approaches to simulate RIAs involve standalone calculations using either a fuel performance code or a thermal-hydraulic code. The complex interdependencies of thermal-hydraulic and fuel mechanical behavior suggest that a tight coupling between these codes may provide more accurate predictions of fuel-to-coolant heat transfer and cladding mechanical response. Here, RELAP5-3D and BISON are coupled in this paper to simulate RIAs, and a sensitivity analysis is performed to rank key thermal properties and two-phase heat transfer parameters relevant for fuel-to-coolant heat transfer and cladding failure mechanisms in UO 2 –Zircaloy-4 systems. Gas gap conductance, film boiling heat transfer uncertainty, pulse width, fuel-specific heat capacity, and cladding-specific heat capacity were identified as important parameters. Variations in figures of merit resulting from changes to pulse width and the material thermal properties indicate that time-dependent heat transfer rates are significant for safety-relevant mechanical parameters due to the time dependence of cladding ductility and pellet-cladding mechanical interaction loading. The results suggest that the thermal-hydraulic factors have a nonnegligible influence on the thermomechanical solution and vice versa. Tight coupling of both sets of physics is recommended to improve prediction of fuel behavior during RIAs. Highlights include the following: 1. The RELAP5-3D thermal-hydraulic code and the BISON fuel performance code are tightly coupled for simulation of RIA transients with energy depositions at the Zircaloy-4 cladding failure threshold. 2. Departure from nucleate boiling occurred for all simulated cases. Due to the ductility of fresh fuel, substantial ballooning occurred in most cases. 3. Gas gap conductance, fuel-specific heat capacity, cladding-specific heat capacity, transient pulse width, and film boiling heat transfer were the dominant thermal factors impacting the safety figures of merit at energy depositions.

Critical Heat Flux (CHF)↗
DOE OSTI2022

Structures of Vanadium-Containing Silicate and Borosilicate Glasses: Vanadium Potential Development and MD Simulations

Vanadium-containing multi-component glasses have been found with widely interests in various applications. However, the complex compositions, thus structures, lead to the challenges in experiments to understand the structural origin of the property change of these glasses. In this work, we developed compatible vanadium-related parameters for a widely used pair wise potential set to enable the simulations of the vanadium-containing multi-component glasses. Various crystal structures and glass compositions (with vanadium in different oxidation states) have been tested using the newly developed parameters, and structural information of cell parameters, pair distribution function, estimated bond distance, and coordination number have been obtained. The results are in good agreement with available experimental data, which indicates the new vanadium parameters can be used to simulate the vanadium-containing multi-component glasses and thus help study those applications with wide interests.

Deng, Lu↗
DOE OSTI2024

AI-assisted detector design for the EIC (AID(2)E)

Artificial Intelligence is poised to transform the design of complex, large-scale detectors like ePIC at the future Electron Ion Collider. Featuring a central detector with additional detecting systems in the far forward and far backward regions, the ePIC experiment incorporates numerous design parameters and objectives, including performance, physics reach, and cost, constrained by mechanical and geometric limits. This project aims to develop a scalable, distributed AI-assisted detector design for the EIC (AID(2)E), employing state-of-the-art multiobjective optimization to tackle complex designs. Supported by the ePIC software stack and using G EANT 4 simulations, our approach benefits from transparent parameterization and advanced AI features. The workflow leverages the PanDA and iDDS systems, used in major experiments such as ATLAS at CERN LHC, the Rubin Observatory, and sPHENIX at RHIC, to manage the compute intensive demands of ePIC detector simulations. Tailored enhancements to the PanDA system focus on usability, scalability, automation, and monitoring. Ultimately, this project aims to establish a robust design capability, apply a distributed AI-assisted workflow to the ePIC detector, and extend its applications to the design of the second detector (Detector-2) in the EIC, as well as to calibration and alignment tasks. Additionally, we are developing advanced data science tools to efficiently navigate the complex, multidimensional trade-offs identified through this optimization process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

General overview of multiphysics modeling and simulation at CEA/DES/IRESNE supporting safety, operating and development of nuclear reactors and facilities

This paper presents a general overview of the innovative multiphysics research/development activities carried on by the CEA/DES/IRESNE institute. Eight main different multiphysics fields addressed by IRESNE have been listed (among others: fuel fabrication, behavior of fuels, of reactors, of corium, radionuclide transfer to environment,...) in relation with the involved coupled physics: neutron transport, fluid mechanics, electro-magnetism, heat transfer, solid mechanics and (physico-)chemistry. The major mission of IRESNE institute is to study current and future nuclear systems integrated into a low carbon energy system. To fulfill this mission, IRESNE's activities rely on the development and the implementation of multiphysics modeling, advanced coupling scheme, data transfer techniques and dedicated multiphysics experimentations. Indeed, such multiphysics simulation approach of nuclear systems is today required to deal with the complex intrinsic features of studied phenomena. This approach aims at taking into account physical phenomena interactions, at each interest modeling scale, in order to reach an accurate and predictive behavior representation of a system component or of the whole system. The final objective is to contribute to the development of digital reactor, core and fuel at CEA, by progressing on applied mathematics, high-performance computing and physical mechanism understanding.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Metrics as tools for bridging climate science and applications

In climate science and applications, the term “metric” is used to describe the distillation of complex, multifaceted evaluations to summarize the overall quality of a model simulation, or other data product, and/or as a means to quantify some response to climate change. Metrics provide insights into the fidelity of processes and outcomes from climate models and can assist with both differentiating models' representation of variables or processes and informing whether models are “fit for purpose.” Metrics can also provide a valuable reference point for co-production of knowledge between climate scientists and climate impact practitioners. Although continued metric developments enable model developers to better understand the impacts of decisions made in the model design process, metrics also have implications for the characterization of uncertainty and facilitating analyses of underlying physical processes. As a result, comprehensive evaluation with multiple metrics enhances usability of climate information by both scientific and stakeholder communities. Here, this paper presents examples of insights gained from the development and appropriate use of metrics, and provides examples of how metrics can be used to engage with stakeholders and inform decision-making.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Computational Investigations of the Reactivity of Metalloporphyrins for Ammonia Oxidation

Density functional theory and molecular dynamics simulations were used to assess the ability of tetraphenylporphyrin (TPP)M complexes (where M = Cr, Mn, Fe, Co, Ni, Mo, Ru, W, and Os) to coordinate and weaken the N-H bonds of ammonia, as well as their reactivity towards N-N bond formation for N 2 generation. Compared to other metalloporphyrins, bis-ammonia complexes (TPP)Mo(NH 3 ) 2 and (TPP)W(NH 3 ) 2 exhibit low and level N-H BDFEs due to a stabilized (TPP)M(NH 3 )(NH) intermediate by multiple metal-ligand bonding. These results resemble those previously obtained for polypyridyl metal complexes, suggesting that broad trends in reactivity towards N-H bond cleavage are more metal-dependent rather than ligand-dependent for a metal in a nitrogen pseudo-octahedral environment. We investigated N-N bond formation via NH 3 nucleophilic attack on M-NH and M-N intermediates, compared to bimolecular coupling of M-NH x intermediates. We evaluated the reactivity of (TPP)Fe(NH 3 ) 2 towards N-N bond formation via a hydrazine pathway, and found amide-amide coupling to form a bridged hydrazido complex to be the most favorable pathway for N-N bond formation. Further investigation of possible N-N bond formation pathways by reaction with NH 3 led us to identify a possible Fe III -•NH species with significant aminyl character that bypasses the nucleophilic attack of NH 3 and that promotes homolytic N-H bond cleavage of ammonia. This reaction forms a Fe-NH 2 moiety and a transient •NH 2 radical that subsequently forms an N-N bond with the Fe-NH 2 moiety to form a (TPP)Fe(NH 3 )(N 2 H 4 ) species. Furthermore, these results indicate the need to evaluate the radical character of imido species and their reactivity towards N-H bond cleavage of ammonia.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Evaluating the grain-scale deformation behavior of a single-phase FCC high entropy alloy using synchrotron high energy diffraction microscopy

We report although the deformation behavior of high-entropy alloys (HEAs) has been extensively studied at the macroscale, many important properties have yet to be explored for these alloys at the microscale, thus hampering accurate prediction of damage and failure processes. Synchrotron high-energy diffraction microscopy (HEDM) and fast-Fourier transform-based crystal plasticity modeling was conducted to investigate the three-dimensional (3D) grain-resolved micromechanical response for approximately 1,900 constituent grains within a single-phase FCC HEA up to 1% applied strain. The evolution of grain-resolved elastic strains, lattice reorientations, and maximum resolved shear stresses (mRSS) were evaluated to quantify elastic, yield, and fully plastic behavior. Overall, the initial critical resolved shear stress (CRSS), determined via in situ HEDM and companion modeling, was found to be > 20% higher than estimated using the classical polycrystalline Taylor factor (M = 3.06). However, a descriptive parameter based on the average grain-resolved Taylor factor ($\overline{M}$) was found to show excellent agreement with plastic yielding behavior observed within HEDM datasets. Noticeable deviations in HEDM lattice reorientations compared to both EVP-FFT simulations and classical predictions for FCC polycrystals were discovered, highlighting the complexity in correlating local lattice reorientations, Taylor, and Schmid factors with plastic response for this material at the grain-scale. Therefore, it is anticipated that the overall trends and parameter identification of 3D grain-resolved properties in this study can serve as an important foundation for continued mesoscale investigation on both well-established and newly developed Cantor-like HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Scanning the RBD-ACE2 molecular interactions in Omicron variant

The emergence of new SARS-CoV-2 variants poses a threat to the human population where it is difficult to assess the severity of a particular variant of the virus. Spike protein and specifically its receptor binding domain (RBD) which makes direct interaction with the ACE2 receptor of the human has shown prominent amino acid substitutions in most of the Variants of Concern. Here, in this work, by using all-atom molecular dynamics simulations we compare the interaction of Wild-type RBD/ACE2 receptor complex with that of the latest Omicron variant of the virus. We observed a very interesting diversification of the charge, dynamics and energetics of the protein complex formed upon mutations. These results would help us in understanding the molecular basis of binding of the Omicron variant with that of SARS-CoV-2 Wild-type.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗