Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials 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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

COMPUTER-AIDED LATTICE DESIGN AND ADVANCED MODELING FOR THE DEVELOPMENT OF LIGHTWEIGHT STRUCTURAL MATERIALS

The notion of a “fission battery” conveys a vision focused on realizing very simple “plug-and-play” nuclear systems that can be integrated into a variety of applications requiring affordable, reliable energy in the form of electricity and/or heat and function without operations and maintenance staff. Fission batteries require lightweight structural materials to increase their mobility, and the lightweight materials must demonstrate structural resilience under various conditions. The objective of this work is to develop lattice structured lightweight structural material featuring a good combination of mechanical properties using advanced modeling and simulation together with an advanced additive manufacturing technique such as laser powder bed fusion. The preliminary results show that different lattice structures and types can be successfully meshed using nTopology software, and the lattice structure data can be successfully transformed to Multiphysics Object-Oriented Simulation (MOOSE) Environment input. Finite Element Analysis (FEA) displays that, at macro/engineering scale simulation, the weight saving design has an obvious effect on tensile behavior such as effective elastic modulus and yield stress. The novel approaches of this work are (1) development of lattice structures for improved mechanical properties using advanced simulation and modeling techniques; and (2) model predictions of the mechanical properties (e.g., strength and stress distribution) of macroscopic materials in order to preliminarily select a lattice structure for additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of Lightweight Structural Materials with Improved Properties for Fission Batteries

The notion of a “fission battery” conveys a vision focused on realizing very simple “plug-and-play” nuclear systems that can be integrated into a variety of applications requiring affordable, reliable energy in the form of electricity and/or heat and function without operations and maintenance staff. Fission batteries require lightweight structural materials to increase their mobility, and the lightweight materials must demonstrate structural resilience under various conditions. The objective of this work is to develop lightweight structural material featuring a good combination of mechanical properties using advanced modeling and simulation together with an advanced additive manufacturing technique such as laser powder bed fusion. The preliminary results show that different lattice structures and types can be successfully meshed using nTopology software, and the lattice structure data can be successfully transformed to Multiphysics Object-Oriented Simulation Environment (MOOSE) input. Finite Element Analysis (FEA) displays that, at macro/engineering scale simulation, the weight saving design has an obvious effect on tensile behavior such as effective elastic modulus and yield stress. The novel approaches of this work are (1) development of optimized lattice structures for improved mechanical properties using advanced simulation and modeling techniques; and (2) model predictions of the mechanical properties (e.g., strength and stress distribution) of macroscopic materials in order to preliminarily select a lattice structure for additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Extended Finite Element Based Approach in Additive Manufacturing Modeling for Optimizing Highly Complex Manifold in Protonic Ceramic Electrochemical Cells

The objective of this project is to develop a novel numerical approach based on the extended finite element method (XFEM) for modeling moving boundaries of material deposited in additive manufacturing (AM) process with improved accuracy and reduced computational cost. One of the major challenges in simulating AM processes is that the boundaries of the computational domain need to change as material is deposited. Previous approaches have achieved this by activating new finite elements, but this requires a high level of mesh refinement to capture small movement of the boundary. XFEM allows the solution boundary to move independently of the mesh, permitting smooth representation of the evolution of the boundary as material is deposited. As a result, this new method will significantly improve the simulation accuracy while reducing the computational cost compared with existing approaches. Furthermore, this project will bring significant improvements for the design and AM process optimization in an iterative fashion among numerical simuation, parameter optimization, and experimental validation. Specifically, this project will support the protonic ceramic electrochemical cells (PCEC) stack development at INL by providing insights into the density, residual stress, thermal and mechanical properties of the PCEC manifold and interconnect, which is critical in enhancing the system lifetime and reducing the overall cost. Upon the success of proposed development and validation of the proposed approach, simulation will be applied to determine an optimal set of AM process parameters for the PCEC manifold production. This project will consolidate Idaho National Laboratory (INL)’s simulation capabilities provided by the MOOSE framework and the Valhalla AM simulation application, with encouraging expansion to emerging PCEC applications and AM technologies at INL and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Exploring the Structural, Dynamic, and Functional Properties of Metal–Organic Frameworks through Molecular Modeling

This review spotlights the role of atomic-level modeling in research on metal-organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster-based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high-throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Efficient fully-coherent quantum signal processing algorithms for real-time dynamics simulation

Simulating the unitary dynamics of a quantum system is a fundamental problem of quantum mechanics, in which quantum computers are believed to have significant advantage over their classical counterparts. One prominent such instance is the simulation of electronic dynamics, which plays an essential role in chemical reactions, non-equilibrium dynamics, and material design. These systems are time-dependent, which requires that the corresponding simulation algorithm can be successfully concatenated with itself over different time intervals to reproduce the overall coherent quantum dynamics of the system. In this paper, we quantify such simulation algorithms by the property of being fully-coherent: the algorithm succeeds with arbitrarily high success probability 1 − δ while only requiring a single copy of the initial state. Here we subsequently develop fully-coherent simulation algorithms based on quantum signal processing (QSP), including a novel algorithm that circumvents the use of amplitude amplification while also achieving a query complexity additive in time t, ln(1/δ), and ln(1/ϵ) for error tolerance ϵ: $Θ‖\mathscr{H}‖|t|+ln(1/ϵ)+ln(1/δ)$. Furthermore, we numerically analyze these algorithms by applying them to the simulation of the spin dynamics of the Heisenberg model and the correlated electronic dynamics of an H2 molecule. Since any electronic Hamiltonian can be mapped to a spin Hamiltonian, our algorithm can efficiently simulate time-dependent ab initio electronic dynamics in the circuit model of quantum computation. Accordingly, it is also our hope that the present work serves as a bridge between QSP-based quantum algorithms and chemical dynamics, stimulating a cross-fertilization between these exciting fields.

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

Quench Simulation and Validation for EIC Direct Wind Magnets

The Electron-Ion Collider (EIC), a new facility to be built in the United States at the U.S. Department of Energy’s Brookhaven National Laboratory in collaboration with Thomas Jefferson National Accelerator Facility-will be the only new particle collider in the U.S. for the next 10∼20 years. The EIC project will investigate the origin of mass and building blocks of the Universe. It will accelerate and collide beams of electrons with protons and nuclei. Unlike any collider that has ever existed before, both EIC’s particle beams will be polarized, meaning their spins will be aligned in a controlled way to allow for precision measurements. Different types of new superconducting magnets will be required near the interaction region (IR) of EIC. A 6-around-1 NbTi cable will be used for some of those new superconducting magnets. Given the relatively smaller conductors used for those direct wind magnets, the overall engineering current density will be 2–3 times higher in comparison to traditional collared superconducting magnets in other accelerator magnets. This is technically challenging and potentially represents high risks for quench protection. Here, in this paper we use different tools and methods for quench simulation for several EIC direct wind magnets. We also present results of a sensitivity study, exploring different material properties and Jc curves. This paper also discusses and validates the quench simulation results with measurement results from a small-scale serpentine magnet tested at BNL in October 2024.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Enabling Efficient Water Splitting with Advanced Materials Designed for High pH Membrane Interface

This project was focused on developing the durable, high-performance materials and interfaces for advanced water splitting, enabling a clear pathway for achieving <$2/Kg H2 (on scale) with efficiency of 43 kWh/kg H 2 via anion exchange membrane (AEM)-based electrolysis. We aimed to advance these final goals via an improved fundamental understanding of both hydrogen and oxygen evolution reactions (HER/OER) leading to novel platinum group metal (PGM)-free catalyst materials in conjunction with critical improvements in membrane and ionomers and gas evolution electrodes with corresponding characterization and testing. Northeastern University (NU) lead this effort focusing on catalyst development and characterization (both in situ and ex situ) while project partners lead improvements in ionomer and membrane materials and will aid in the development of specialized electrode and membrane electrode assemblies. In addition, close collaboration occured with the HydroGEN Energy Materials Network (EMN) National Laboratory consortium including efforts related to use of advanced ionomers, durability protocols and validation of electrolyzer materials (e.g. NREL), multiscale modeling and computation (e.g. LBNL), and molecular dynamics (MD) simulations of the membrane catalyst interface (e.g. SNL). The interactions with HydroGEN included exchange of data and materials as needed to facilitate project success.

08 HYDROGEN↗
DOE OSTI2025

Elasto-viscoplastic fast Fourier transform modeling framework for assessing microstructural effects on stress intensity factors characterizing fracture toughness

A large-strain elasto-viscoplastic fast Fourier transform (LS-EVPFFT) model with non-periodic (NP) velocity-based boundary conditions is adapted to simulate the sensitivity of stress intensity factors on microstructure for 304L stainless steel. The material was characterized via electron backscattered diffraction (EBSD) serial-sectioning to obtain a measured 3-D microstructural cell to perform simulations. The NP-LS-EVPFFT model, including the simulation setup and boundary conditions, was verified using a crystal plasticity finite element (CPFE) model. To this end, the generation of meshes of notched specimens was developed, which involved creating Python scripts for mesh “cutting” in Abaqus, and Sculpt scripts in Cubit for meshing of the measured microstructural cell processed with DREAM.3D. The complexity of the mesh preparation highlighted the advantages of the FFT-based model, which circumvents the mesh generation process. Given the efficiency of the FFT-based model, statistical distribution of stress intensity factors in function of crystal orientation at the crack tip, grain structure, and crystallographic texture surrounding the crack tip were predicted. Further, the distributions reveal about 10% variation of stress intensity factors with microstructure with the most significant sensitivity found to be the crystal orientation at the crack tip. The methodology developed in this work is discussed as a practical simulation tool for predicting the sensitivity of stress intensity factors on microstructural variability in metallic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accurate frozen core approximation for all-electron density-functional theory

We implement and benchmark the frozen core approximation, a technique commonly adopted in electronic structure theory to reduce the computational cost by means of mathematically fixing the chemically inactive core electron states. The accuracy and efficiency of this approach are well controlled by a single parameter, the number of frozen orbitals. Explicit corrections for the frozen core orbitals and the unfrozen valence orbitals are introduced, safeguarding against seemingly minor numerical deviations from the assumed orthonormality conditions of the basis functions. A speedup of over twofold can be achieved for the diagonalization step in all-electron density-functional theory simulations containing heavy elements, without any accuracy degradation in terms of the electron density, total energy, and atomic forces. This is demonstrated in a benchmark study covering 103 materials across the Periodic Table and a large-scale simulation of CsPbBr3 with 2560 atoms. Our study provides a rigorous benchmark of the precision of the frozen core approximation (sub-meV per atom for frozen core orbitals below −200 eV) for a wide range of test cases and for chemical elements ranging from Li to Po. The algorithms discussed here are implemented in the open-source Electronic Structure Infrastructure software package.

Yu, Victor Wen-zhe↗
DOE OSTI2022

A mechanism for reduced compression in indirectly driven layered capsule implosions

High-yield implosions on the National Ignition Facility rely on maintaining low entropy in the deuterium–tritium fuel, quantified by its adiabat, in order to efficiently couple energy to the hot spot through high compression of the fuel layer. We present very-high-resolution xRAGE simulation results that study the impacts of interfacial mixing and the jetting of materials due to surface defects, defects on internal interfaces, voids, and engineering features on fuel layer compression. Defects and voids are typically neglected in implosion simulations due to their small size and three-dimensional geometry. Our results showed that supersonic jets of material arise through weak spots in the shell at peak implosion velocity that prevent uniform compression of the fuel layer even when they do not introduce contaminant into the hot spot. This occurs despite maintaining low fuel entropy, since the formation of the weak spots involves nonradial displacement of fuel mass. In contrast, simulations show that fuel–ablator mixing due to interfacial instabilities has a much smaller impact on compression. We show that defects on interior interfaces of plastic capsules decrease compression by 15% to 25% and interfacial mixing between the ablator and fuel decreases compression by less than 1% for implosions with plastic or high-density carbon (HDC) ablators. For low adiabat implosions, the impact of jetting seeded by the support tent can also decrease the compression by 25%. We demonstrate that the inclusion of interior defects in simulations can explain the inferred compression in two fielded plastic capsule implosions and that the inclusion of voids, for which available characterization has large uncertainties, in simulations of HDC capsule implosions has a qualitatively consistent impact. This mechanism offers a potential explanation for persistently overestimated fuel compression in design simulations of layered implosions on the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Pressure-driven density match nucleates metastable r8 phases from amorphous Si and Ge

The pressure–temperature phase behavior of covalent disordered solids such as amorphous silicon and germanium is complex. Questions remain on possible glass transitions, on polyamorphism via amorphous–amorphous transitions, on connections with liquid–liquid transitions, on structure-behavior relationships, and on their potential as precursor for novel methods for material discovery. Here we demonstrate experimentally the nucleation of a metastable, four-fold coordinated rhombohedral r8 phase from pure amorphous silicon and germanium upon room temperature compression at pressures below 10 GPa. Accompanying theory reveals a strong pressure-driven distortion of the bond angle transforming the starting tetrahedral low-density amorphous network to a distorted four-fold coordinated medium-density state. This state is of lower density than metallic high-density networks, resembles the crystalline r8 phase and initiates its nucleation. Our finding shows that polyamorphism is not the only possible transformation mode for these amorphous solids and that instead nucleation of interesting functional phases at potentially useful pressures is possible. Such novel access modes to metastable structures are critical for future exploitability and could be useful for other tetrahedral materials including carbon, where the related (bc8) post-diamond phase remains elusive. Our observed density match between an amorphous and a metastable crystalline phase clearly allows for a new phase transition pathway, while corresponding theory demonstrates how carefully validated atomistic simulations can guide prediction, discovery and synthesis of novel material structures.

Materials discovery↗
DOE OSTI2021

Computational study of Li3BO3 and Li3BN2 I: Electrolyte properties of pure and doped crystals

Both Li3BO3 and Li3BN2 materials have promising properties for use in all-solid-state batteries and other technologies dependent on electrolytes with significant ionic conductivity. As the first of a two-part study, this paper reports the analysis of detailed simulations of Li ion diffusion in the monoclinic forms of these materials. Using both NEB and MD methods, it is clear that Li ion migration via vacancy mechanisms provides the most efficient ion transport in each material. While the results suggest that interstitial defects in these materials do not play a direct role in Li ion migration, their relative stability seems to enhance vacancy production via the formation of Frenkel-type defects. This may partially explain why the Li ion conductivities computed from MD simulations of samples initially containing a single Li ion vacancy are in reasonable agreement with measured values of this work for Li3BO3 and those reported in the literature for poorly crystalline samples of both materials. The possibility of increasing vacancy concentrations by substitutional doping (F for O in Li3BO3 and C for B in Li3BN2) is also examined, finding simulated conductivities comparable to those of the ideal vacancy model.

Li, Yan↗
DOE OSTI2025

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting​

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Liquid-Vapor Coexistence and Critical Point of Mg 2 SiO 4 From Ab Initio Simulations

Hypervelocity impact-driven vaporization is characteristic of late-stage planet formation. Yet the behavior and properties of liquid-vapor mixtures of planetary materials of interest are typically unknown. Multiphase equations of state used in hydrodynamic simulations of planet impacts therefore lack reliable data for this important phenomenon. Here, we present the first constraints on the liquid-vapor critical point and coexistence phase boundary of Mg2SiO4 computed from ab initio molecular dynamics simulations. We found that the vapor is depleted in magnesium and enriched in silica and oxygen, while the coexisting liquid is enriched in magnesium and depleted in oxygen, from which we infer vaporization is incongruent. The critical point was estimated from an equation of state fit to the data. The results are in line with recent calculations of MgSiO3 and together confirm that extant multiphase equation of state (EOS) models used in planetary accretion modeling significantly underestimate the amount of supercritical material postimpact.

58 GEOSCIENCES↗
DOE OSTI2022

Gradient Architecture Design in Scalable Porous Battery Electrodes

Because it has been demonstrated to be effective toward faster ion diffusion inside the pore space, low-tortuosity porous architecture has become the focus in thick electrode designs, and other possibilities are rarely investigated. Here, to advance current understanding in the structure-affected electrochemistry and to broaden horizons for thick electrode designs, we present a gradient electrode design, where porous channels are vertically aligned with smaller openings on one end and larger openings on the other. With its 3D morphology carefully visualized by Raman mapping, the electrochemical properties between opposite orientations of the gradient electrodes are compared, and faster energy storage kinetics is found in larger openings and more concentrated active material near the separator. As further verified by simulation, this study on gradient electrode design deepens the knowledge of structure-related electrochemistry and brings perspectives in high-energy battery electrode designs.

25 ENERGY STORAGE↗
DOE OSTI2021

BISON Simulation Development for ALD Coated Particles (Progress Report, FY21)

Argonne has on ongoing effort to perform atomic layer deposition coatings on micron-scale fuel particles. Initial results showed cracking of the coating layer above a specific coating thickness, which motived the development of a BISON model for the coated particle system to help explain the behavior. This report describes the initial development of the BISON model, the materials models used, and the conditions used in the simulation. A 2D model has been developed, with sensitivity studies performed on several key parameters. Based on the 2D model results, and 3D model was also developed, with results from all calculations described. First principles calculations were also performed on the fuel/coating interface to help describe the observed behavior. Potential future activities are also described.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Cold-Sprayed NMC622 Composite as a Cathode for Lithium-Ion Batteries

The growing demand for high-energy, low-cost lithium-ion batteries (LIBs) to power electric vehicles (EVs) necessitates advances in both materials and manufacturing processes. Conventional cathode fabrication methods, such as slurry casting and drying, are energy-intensive and pose challenges for scalability and environmental compliance. In this study, we propose cold-spray (CS) deposition as a solvent-free approach for fabricating LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) composite cathodes. Powder blends of NMC622, poly(vinylidene fluoride) (PVDF), and carbon black (CB) are directly deposited onto stainless steel and Inconel substrates under varied gas temperatures, pressures (and thus velocities), and standoff distances. The effects of temperature on the deposit morphology, coating density, and volume are systematically investigated. Computational fluid dynamics simulations reveal that increasing the gas temperature enhances the particle velocity, narrows the spray angle, and reduces the mass concentration radially at the nozzle outlet. CS deposition results in a dense cathode microstructure, accompanied by fracture in polycrystalline NMC622 particles. X-ray diffraction analysis further verifies that there are no phase changes during the deposition process. The electrochemical performance of the cold-sprayed cathodes reveals an initial capacity of approximately 96 mAh g –1 for single-crystal NMC622 and 167 mAh g –1 for polycrystalline NMC622. While these values are modest compared to state-of-the-art slurry-cast cathodes, which typically exhibit 180–200 mAh g –1 under optimized conditions. The results demonstrate a competitive performance given the solvent-free nature of the CS process and compare favorably with tape-cast samples made from identical feedstock. In conclusion, he CS process enables the formation of dense, binder-integrated cathode coatings without the need for solvent processing, offering a promising pathway for scalable, energy-efficient dry electrode manufacturing of next-generation LIBs.

Batteries↗
DOE OSTI2024

Physical properties of asteroid Dimorphos as derived from the DART impact

On 26 September 2022, NASA’s Double Asteroid Redirection Test (DART) mission successfully impacted Dimorphos, the natural satellite of the binary near-Earth asteroid (65803) Didymos. Numerical simulations of the impact provide a means to find the surface material properties and structures of the target that are consistent with the observed momentum deflection efficiency, ejecta cone geometry and ejected mass. Our simulation that best matches the observations indicates that Dimorphos is weak, with a cohesive strength of less than a few pascals, like asteroids (162173) Ryugu and (101955) Bennu. We find that the bulk density of Dimorphos ρ B is lower than ~2,400 kg m – 3 and that it has a low volume fraction of boulders (≲40 vol%) on the surface and in the shallow subsurface, which are consistent with data measured by the DART experiment. These findings suggest that Dimorphos is a rubble pile that might have formed through rotational mass shedding and reaccumulation from Didymos. Our simulations indicate that the DART impact caused global deformation and resurfacing of Dimorphos. ESA’s upcoming Hera mission may find a reshaped asteroid rather than a well-defined crater.

79 ASTRONOMY AND ASTROPHYSICS↗