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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Superdiffusion resilience in Heisenberg chains with two-dimensional interactions on a quantum processor

Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋⁢𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.

Alagarsamy Manikandan, Keerthi Kumaran [ORNL]↗
DOE OSTI2023

Improved treatment of multi-material cells in thermal radiation transport codes

High-energy-density physics simulations with non-conformal meshes of materials require a multi-material (MM) closure that affects the thermal radiation transport (TRT). We propose a set of novel closures that work for an arbitrary number of materials, both grey and multigroup energy discretizations, and any angular discretization (such as Sn, IMC, or diffusion). For each spatial cell, our closures let each species (ion and electron) of each material have its own temperature, density, and internal energy, but use a single radiation distribution that interacts with all materials within the cell. Our closures maintain energy conservation, do not incur increased computational cost in the TRT solve itself, are compatible with single-material TRT solvers, and do not make any temperature-equilibrium assumptions. Here we test our closures on a wide range of increasingly realistic problems and find them to be robust.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2021

On-Sun Testing of a High-Temperature Solar Receiver's Flux Distribution

Concentrated solar power (CSP) is a promising technology in transitioning to renewable energy because of its abundance in nature and thermal energy storage (TES) capability. Among the four types of available CSP technology, including parabolic trough, linear Fresnel, power tower, and parabolic dishes, a power tower using a central receiver has more potential to generate high-temperature heat in a scale supporting power cycles efficiency and achieve low levelized cost of energy (LCOE). Other than the conventional type of receiver design, the high-absorptive receiver concept developed and presented in this paper is novel in its design approach. The novel receiver design originated from National Renewable Energy Laboratory (NREL) consists of an array of solar flux absorb tubes. The solar absorb tubes require uniform flux distribution and in-depth flux penetration through the three different reflective sections of tubes in a hexagonal shape. To evaluate this unique receiver design and thermal performance, the flux distribution, flux uniformity, and intensity were numerically simulated using ansys fluent and SolTrace modeling program. On-sun testing has been done at NREL high flux solar testing facility, based on the computational analysis.

clean energy↗
DOE OSTI2021

Exceptionally high spallation strength for a high-entropy alloy demonstrated by experiments and simulations

High-entropy alloys are materials with an increasing number of technological applications. Amongst them, the Cantor alloy, FeMnCoCrNi, shows desirable mechanical properties at normal loading conditions. In this study we focus on the performance of the Cantor alloy at the ultra-high deformation rates of shock waves. We study shock-induced spallation using both experiments and atomistic simulations. Experimental loading is achieved using high power laser, with VISAR to obtain velocity profiles and spall strength, followed by transmission electron microscopy of the recovered samples. Molecular Dynamics (MD) simulations of shock-induced spallation are compared with experiments. Both experiments and simulations show a high spall strength which would be beneficial for certain applications, with experiments giving ~8 GPa at ~10 7 s -1 and MD giving almost ~30 GPa at ~10 9 s -1 . The difference between experiments and simulations can be explained by the difference in strain rate. Post-mortem analysis of the experimental samples shows nanotwins near the spall plane, while MD simulations show a highly disordered region giving rise to void nucleation and spall during loading.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Assessment of simulated and observed cavitation-induced erosion damage in Spallation Neutron Source target vessels

Cavitation-induced erosion damage in different Spallation Neutron Source (SNS) target designs are simulated using explicit finite element–based techniques and compared with observations of erosion in targets after operation. The efficacy of the previously developed method, called saturation time, was evaluated using erosion-damaged samples from new target designs. A new metric called maximum bubble size was implemented under the rationale that larger cavitation bubbles will collapse more intensely. The maximum cavitation bubble size over 1 ms of simulated time was calculated based on the Rayleigh–Plesset equation for each element integration point and presented as a contour map at the vessel surface for assessing with erosion observations. SNS targets are now operated with helium gas injection to reduce cavitation damage. A simulation method using a material model for the mixture of mercury and gas bubbles was recently developed and used to account for the effect of small gas bubbles on the structural response of the target vessel. Furthermore, this work compares the new method's results with observed cavitation damage. Maps of the calculated maximum bubble size for targets operated with and without gas injection were compared with photographs of erosion damage observed in SNS targets. The patterns in maximum bubble size maps correlated well with observations of erosion patterns in target vessels after service. Advantages and challenges of the maximum bubble size simulation technique are provided, and differences between results from the previous and the newly proposed metric are discussed.

Jiang, Hao↗
DOE OSTI2024

Fusion Energy Research at Idaho National Laboratory: Experimentation and Simulation to Support Safety and Rapid Technology Development

Research into fusion energy is growing rapidly, responding to a call for sustainable sources of energy to replace fossil fuels and mitigate climate change. Within the United States, at least, researchers are also responding to the “Bold Decadal Vision” proposed by the White House, seeking to have a commercially relevant fusion pilot plant deployed within a decade. Before this can become a reality, many Fusion Science & Technology (FS&T) gaps remain. For over 45 years, Idaho National Laboratory has been at the forefront of addressing these FS&T gaps in the context of fusion safety and technology via the operation of world-leading experimental facilities within the Safety and Tritium Applied Research (STAR) Facility. Here, INL focuses on the tritium fuel cycle, conceptual system design studies, risk assessment, waste management, and materials safety. Modeling and Simulation (M&S) has also been a component of this portfolio of research, but, early on, focused on individual systems. Since 2019, active development and research on integrated whole device modeling tools based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework has been undertaken. This has culminated in a MOOSE-based version of the Tritium Migration and Analysis Program (TMAP), an INL code historically focused on tritium permeation and trapping within fusion systems. More recently, INL Laboratory Directed Research and Development funds have been used to create the Fusion ENergy Integrated multiphys-X (FENIX) code focused on scrape-off layer plasma physics and the first wall of a magnetically confined fusion device. This talk will focus on an overview of INL activities in the FS&T research area, with a particular focus on recent M&S activities and results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Ferrimagnetic spin waves in honeycomb and triangular layers of Mn{sub 3}Si{sub 2}Te{sub 6}.

A detailed analysis of the ferrimagnetic ground state of Mn 3 Si 2 Te 6 has been performed using inelastic neutron scattering. Although the proposed valence of the nominal Mn 2+ ions would have quenched orbital angular momentum, a significant exchange anisotropy exists in Mn 3 Si 2 Te 6 . This apparent exchange anisotropy is a manifestation of a weak spin-orbit coupling in the layered material. We employ a detailed simulation of the spin-wave spectrum coupling traditional refinement of dispersion parameters to image analysis techniques, while including Monte Carlo simulations of the instrumental resolution to accurately identify the exchange couplings to the third nearest neighbor. An independent validation of our results is made by comparing our final Hamiltonian to heat capacity measurements.

Sala, G.↗
DOE OSTI2023

A Sparse Distributed Gigascale Resolution Material Point Method

In this paper, we present a four-layer distributed simulation system and its adaptation to the Material Point Method (MPM). The system is built upon a performance portable C++ programming model targeting major High-Performance-Computing (HPC) platforms. A key ingredient of our system is a hierarchical block-tile-cell sparse grid data structure that is distributable to an arbitrary number of Message Passing Interface (MPI) ranks. We additionally propose strategies for efficient dynamic load balance optimization to maximize the efficiency of MPI tasks. Our simulation pipeline can easily switch among backend programming models, including OpenMP and CUDA, and can be effortlessly dispatched onto supercomputers and the cloud. Finally, we construct benchmark experiments and ablation studies on supercomputers and consumer workstations in a local network to evaluate the scalability and load balancing criteria. We demonstrate massively parallel, highly scalable, and gigascale resolution MPM simulations of up to 1.01 billion particles for less than 323.25 seconds per frame with 8 OpenSSH-connected workstations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Jan-Dec 2021 IC Project: w20_earthquake (Yearly Report)

The simulations conducted under this IC project used an in-house code developed at Los Alamos National Laboratory called HOSS (Hybrid Optimization Software Suite) and its extended version called HOSS-FSIS (fluid-structure interaction solver). HOSS has been under development in LANL’s Geophysics Group since 2004 and was used to simulate the signal propagation in pieces of material with pre-existing defects and also to simulate the near source effects of above ground nuclear explosions. One of the key feature of HOSS is that it is superior to pure finite element method and discrete element method implementations because it allows for a natural transition from continuum to discontinuum when the material breaks or fractures.

58 GEOSCIENCES↗
DOE OSTI2021

Multiwavelength anomalous X-ray diffraction for combined imaging of atomic displacement and strain

The X-ray Bragg coherent diffractive imaging (CDI) technique assumes that the structure factor holds constant over the measured crystal. This approximation breaks down for materials exhibiting variations in the unit-cell configuration, such as piezo- and ferroelectrics. In that case, the strain field cannot be reliably determined from the reconstruction because the lattice deformation and the structure factor contribute concomitantly. Proposed here is a solution to this problem achieved by combining Bragg CDI and the multiwavelength anomalous diffraction approach that measures a Friedel pair of reflections at two different photon energies near an absorption edge. Comparing the obtained reconstructions with a parametric model that includes calculating the scattering amplitude as a function of wavelength and the unit-cell configuration, the contributions of the lattice deformation and the structure factor are separated. Simulations of the ferroelectric material BaTiO 3 demonstrate the possibility of simultaneous probing of the strain and displacement of the Ti atoms. The proposed method opens up an opportunity to apply coherent X-ray diffraction for nanoscale-resolved 3D mapping of polarization domains in micro- and nanocrystals.

lattice deformation↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2023

Computationally Guided Design of Multiple Impurities Tolerant Cathode for Solid Oxide Fuel Cell Applications

Solid oxide fuel cells (SOFCs) have been recognized as promising candidates in the demand of reducing the dependence on fossil fuels, which is one of the solutions to the threat of global warming. It was reported that the direct chemical-to-electrical energy conversion in SOFCs comes with negligible emission of greenhouse gases and up to 60% working efficiency. SOFC has unique advantages over traditional power generation technologies and other types of fuel cells due to the use of relatively inexpensive materials, low sensitivity to poisoning impurities in fuels, relatively high conversion efficiency, and fuel flexibility (like carbon-based fuels). However, during the SOFC operation, the cathode is exposed to intrinsic and extrinsic airborne contaminants at high temperatures (600–900°C), leading to their reaction with electrode materials, irreversible chemical and structural changes, and electrochemical performance degradation. of which CO2, SO2, Cr6+ and H2O are known to be concerns. Thus, developing a comprehensive understanding of the multiple impurities poisoning of SOFC cathodes in the presence of single or multiple impurities is necessary to address the long-term degradation of the SOFC cathode systems at operating conditions. An integrated computational and experimental approach has been utilized in the current work aiming at the comprehensive investigation of the cathode materials under single or multiple impurities poisoning effects, understanding the multiple impurities poisoning mechanism(s), and improving the long-term durability of the SOFC cathodes to develop better performing impurity-tolerant cathode materials. We first examined the agreement between the simulation approach and the experimental observations at various atmospheric and temperature conditions to confirm the reliability of CALPHAD for poisoning simulations in these cathodes. Then, using this approach, we investigated the accelerated testing protocol, which has become the standard method for conducting impurity poisoning experiments involving these cathode materials. We found our simulation approach is able to predict in which systems, i.e., cathode material and treatment environment, the accelerated testing protocols reflect actual operating conditions. Finally, in-depth simulations were done for these three candidate cathodes under different experimental conditions and suggested directions toward alternative cathode materials that have superior impurity tolerance.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Irradiation Testing of Additively Manufactured Materials for LWR Applications: Alloy 718 and 316L Stainless Steel

The objective of this NSUF Project is to assess the changes in irradiated additively manufactured (AM) material properties as compared to non-irradiated material. Type 316L stainless steel and Alloy 718 samples were produced using Direct Metal Laser Melting (DMLM) fabrication. Materials produced from this fabrication method have several potential applications within the nuclear industry as reactor internal repair parts, fuel debris resistant filters, or fuel spacers within existing light water reactors (LWRs). AM materials have been shown to achieve equivalent mechanical behavior in simulated reactor environments as compared to wrought materials, but have significantly more flexibility when it comes to unique design features. The increased component design flexibility makes these AM materials an attractive choice for both current LWR applications as well as for small modular reactor (SMR) designs. Prior to use of these materials in reactor fleet operation, the industry as a whole must evaluate the effects of irradiation on their material properties. Standard 0.4 inch thick Compact Tension specimens and SSJ3 type tensile bars were neutron irradiated at the Advanced Test Reactor to ~1 dpa for the purpose of performing a variety of mechanical tests in a range of simulated environments applicable to LWRs. For the ductile austenitic Type 316L stainless steel, the irradiated data will be used to confirm that the AM process produces materials with properties that are equivalent to wrought materials under testing conditions applicable to LWR operation. Transmission electron microscopy analysis was also performed in order to understand microstructural and microchemical changes induced in each material in response to neutron irradiation. If possible, data collected from these AM 316L samples will be used to remove fluence limits from specifications of ASME code cases for this alloy, which will give vendors much more flexibility in building future components.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

Controlling substrate temperature with infrared heating to improve mechanical properties of large-scale printed parts

Additively manufactured parts made with polymer extrusion techniques can be 50-75% weaker in the z-direction (across layers) than in the x- and y-directions. This has been attributed to poor mobility of polymer chains and a low degree of entanglement across a cold deposition interface. This is particularly a challenge when printing large-scale parts, such as with the Big Area Additive Manufacturing (BAAM) system, because layer times can exceed several minutes. The current work presents a method for controlling the temperature of the substrate material on the BAAM just prior to deposition using infrared heating lamps. Additionally, long layer times were simulated by actively cooling the material following deposition of each layer. The effect of substrate temperature on the z-direction mechanical properties of 20% carbon fiber reinforced acrylonitrile butadiene styrene (ABS) was measured for an initial temperature ranging from 50 °C to 150 °C and a preheated temperature ranging from 150 °C to 220 °C. Infrared preheating proved to be very effective when applied to substrates that had cooled considerably, almost doubling the tensile strength and increasing the fracture toughness by a factor of 7x.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Utilizing Site Disorder in the Development of New Energy-Relevant Semiconductors

We report controlling site disorder in ternary and multinary compounds enables tuning optical and electronic properties at fixed lattice constants and stoichiometries, moving beyond many of the challenges facing binary alloy systems. Here, we consider possible enhancements to energy-related applications through the integration of disorder-tunable materials in devices such as light-emitting diodes, photonics, photovoltaics, photocatalytic materials, batteries, and thermoelectrics. However, challenges remain in controlling and characterizing disorder. Focusing primarily on II–IV–V 2 materials, we identify three metrics for experimentally characterizing cation site disorder. Complementary to these experiments, we discuss simulation methods to understand disordered materials. Nonidealities, such as off-stoichiometry and oxygen incorporation, can occur while synthesizing metastable disordered materials. While nonidealities may seem undesirable, we describe how if harnessed they could provide another knob for tuning disorder and subsequently properties. To illustrate the effects of disorder on device-relevant properties, we provide case examples of disordered materials and their potential in device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Breaking the mold: Overcoming the time constraints of molecular dynamics on general-purpose hardware

The evolution of molecular dynamics (MD) simulations has been intimately linked to that of computing hardware. For decades following the creation of MD, simulations have improved with computing power along the three principal dimensions of accuracy, atom count (spatial scale), and duration (temporal scale). Since the mid-2000s, computer platforms have, however, failed to provide strong scaling for MD, as scale-out central processing unit (CPU) and graphics processing unit (GPU) platforms that provide substantial increases to spatial scale do not lead to proportional increases in temporal scale. Important scientific problems therefore remained inaccessible to direct simulation, prompting the development of increasingly sophisticated algorithms that present significant complexity, accuracy, and efficiency challenges. While bespoke MD-only hardware solutions have provided a path to longer timescales for specific physical systems, their impact on the broader community has been mitigated by their limited adaptability to new methods and potentials. In this work, we show that a novel computing architecture, the Cerebras wafer scale engine, completely alters the scaling path by delivering unprecedentedly high simulation rates up to 1.144 M steps/s for 200 000 atoms whose interactions are described by an embedded atom method potential. This enables direct simulations of the evolution of materials using general-purpose programmable hardware over millisecond timescales, dramatically increasing the space of direct MD simulations that can be carried out. In this paper, we provide an overview of advances in MD over the last 60 years and present our recent result in the context of historical MD performance trends.

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

Mixed Material Scintillator Systems for Neutron Detection

Scintillators are a fundamental technology that enables radiation detection. Their use has historically enabled the detection of incoming gamma radiation, with recent advances in plastic scintillator development enabling detection and discrimination of neutron and gamma radiation through pulse shape discrimination. However, scintillators have so far consisted of a single material, typically one primary dye and one secondary dye embedded in a plastic matrix. Based on the simulation studies described previously, mixed material scintillator systems with multiple secondary dyes arranged in alternating, periodic micro-sized regions opens the door for scintillators with a wide variety of different capability, such as directional sensing, positional reconstruction, particle identification, and pulse shape discrimination. However, fabrication of a mixed material scintillator system by conventional manufacturing methods is difficult to manufacture and scale, creating a barrier to their wide spread implementation. On the other hand, additive manufacturing additive manufacturing presents itself as the optimal method for fabricating mixed material scintillator systems based on its flexibility, ability to fabricate fine features, and potential for multimaterial fabrication.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗