Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Materials simulation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Interlaced Characterization and Calibration (ICC) for Improved Computational Simulation Credibility

Accurate material characterization and model calibration are pivotal for simulations used for high-consequence engineering decisions. Current characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data is collected for a specific model of interest, (3) provide only mean parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work developed a new paradigm—coined Interlaced Characterization and Calibration (ICC)—which drives forward the state-of-the-art in model calibration by bringing together recent advancements into one improved workflow. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) provides uncertainty metrics on the calibrated model parameters, and (4) incorporates these advances into a quasi real-time feedback loop. The ICC framework was validated synthetically with both low-fidelity and high-fidelity simulations paired with several different elastoplastic material models, and was also demonstrated experimentally with an aluminum 6061 cruciform exemplar specimen. Results showed that the ICC framework—in which Bayesian optimal experimental design actively guided the experiment— resulted in calibrations with similar or better accuracy than predetermined experiments based on subject matter expertise. Moreover, the ICC framework produced a complete model calibration— with quantified uncertainties on model parameters—in 1 week, a 5 - 10× increase in efficiency over traditional approaches. Thus, the ICC paradigm improves both the calibration process and quality, by (1) improving efficiency, which increases agility of solid mechanics modeling and enables utilization of computational simulation (CompSim) at earlier stages of the design cycle and (2) providing quantified, and in some cases reduced, parameter uncertainties, which increases confidence in model predictions and supports credible decision making.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

3D Printing of Cement-Based Materials Using Seawater for Simulated Marine Environments

Global demand for adaptable and rapidly deployable construction solutions in offshore, coastal, and fluvial environments continues to rise, driven by pressing needs to develop energy platforms, improve coastal resilience, and support emergency response in the face of natural disasters. Increased investment in human-made coastal infrastructure, such as piers, support structures for power lines, offshore wind farms, and seawall protection systems, further underscores this trend. This study investigates the development of printable concrete mixtures for underwater environments using seawater as a replacement for freshwater, using a 3D printing syringe-based extrusion system. The effect of seawater addition and the printing medium (in air vs. underwater) was assessed via rheological and mechanical performance characterization. The results indicate rheological properties are favorable for seawater adoption by producing mixtures with higher yield stress and viscosity with the same levels of admixtures used for freshwater. Seawater-based mixtures demonstrated superior dimensional stability compared to freshwater counterparts, maintaining cross-sectional geometry, while compressive strength results showed no statistical differences between in-air and underwater samples. However, flexural strength was significantly influenced by geometry and printing medium. These findings establish critical rheological parameters for printable underwater mixtures and highlight the need for optimized curing strategies and layer bonding techniques to improve interfacial strength in underwater 3D printing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Abstract Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Large radiation back-flux from Monte Carlo simulations of fusion neutron–material interactions

Abstract Fusion power reactors will generate intense neutron fluxes into plasma-facing and structural materials (SMs). Radiation back-fluxes, generated from neutron–material interactions under these fluxes, can dramatically impact the plasma dynamics, e.g. by seeding runaway electrons during disruptions via Compton scattering of background electrons by wall-emitted gamma radiation. Here, we quantify these back-fluxes, including neutrons, gamma rays, and electrons, using Monte Carlo calculations for a range of SM candidates and first wall (FW) thicknesses. The radiation back-flux magnitudes are remarkably large, with neutron and gamma radiation back-fluxes on the same order of magnitude as the incident fusion neutron flux. Electron back-fluxes are two orders of magnitudes lower, but are emitted at sufficiently high energies to impact the sheath and boundary plasma dynamics. Material configuration plays a key role in determining back-flux magnitudes. The SM chiefly determines the neutron back-flux magnitude, while the FW thickness principally attenuates the gamma ray and electron back-fluxes. In addition to prompt back-fluxes, which are emitted immediately after fusion neutrons impact the surface, significant delayed gamma ray and electron back-fluxes arise from nuclear decay processes in the activated materials. These delayed back-flux magnitudes range from 2% to 7% of the prompt back-fluxes, and remain present during transients when fusion no longer occurs. During disruptions, build-up of delayed gamma radiation back-flux represents potential runaway electron seeding mechanisms, posing additional challenges for disruption mitigation in a power reactor compared with non-nuclear plasma operations. This work highlights the impact of these radiation back-fluxes plasma performance and demonstrates the importance of considering back-flux generation in materials selection for fusion power reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Phase-field modeling of crack growth under coupled creep-fatigue

Crack growth under coupled creep-fatigue strongly influences service life of metallic components operating at high temperatures. However, there lack numerical models that enable direct simulations of crack growth and plasticity development under coupled creep-fatigue. Here, in this study, a unified phase-field model is developed for simulating crack growth under fatigue, creep, and coupled creep-fatigue. This model is able to reproduce the Paris laws for both brittle and ductile materials in pure fatigue mode, which shows that the degradation rate due to cyclic stress in ductile materials is generally larger than that in brittle materials. Fatigue simulations with different yield strengths demonstrate that plasticity facilitates crack growth in typical ductile materials. High-throughput simulations are performed under coupled creep-fatigue conditions, with varying hold-stress levels and hold-time. By analyzing crack growth rate under different hold-time, a creep-fatigue interaction term is obtained. Simulation results demonstrate that the creep-fatigue interaction is caused by smaller stress gradient and enhanced degradation rate ahead of the crack tips under coupled creep-fatigue, compared to the situation of pure fatigue. This work reveals the origin of different Paris-law exponents in brittle and ductile materials and demonstrates how the interplay between creep and fatigue affects crack growth.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2020

Towards simulation of the dynamics of materials on quantum computers

A highly anticipated application for quantum computers is as a universal simulator of quantum many-body systems, as was conjectured by Richard Feynman in the 1980s. The last decade has witnessed the growing success of quantum computing for simulating static properties of quantum systems, i.e., the ground state energy of small molecules. However, it remains a challenge to simulate quantum many-body dynamics on current-to-near-future noisy intermediate-scale quantum computers. In this work, we demonstrate successful simulation of nontrivial quantum dynamics on IBM's Q16 Melbourne quantum processor and Rigetti's Aspen quantum processor; namely, ultrafast control of emergent magnetism by THz radiation in an atomically-thin two-dimensional material. The full code and step-by-step tutorials for performing such simulations are included to lower the barrier to access for future research on these two quantum computers. As such, this work lays a foundation for the promising study of a wide variety of quantum dynamics on near-future quantum computers, including dynamic localization of Floquet states and topological protection of qubits in noisy environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The joint automated repository for various integrated simulations (JARVIS) for data-driven materials design

The Joint Automated Repository for Various Integrated Simulations (JARVIS) is an integrated infrastructure to accelerate materials discovery and design using density functional theory (DFT), classical force-fields (FF), and machine learning (ML) techniques. JARVIS is motivated by the Materials Genome Initiative (MGI) principles of developing open-access databases and tools to reduce the cost and development time of materials discovery, optimization, and deployment. The major features of JARVIS are: JARVIS-DFT, JARVIS-FF, JARVIS-ML, and JARVIS-tools. To date, JARVIS consists of ≈40,000 materials and ≈1 million calculated properties in JARVIS-DFT, ≈500 materials and ≈110 force-fields in JARVIS-FF, and ≈25 ML models for material-property predictions in JARVIS-ML, all of which are continuously expanding. JARVIS-tools provides scripts and workflows for running and analyzing various simulations. We compare our computational data to experiments or high-fidelity computational methods wherever applicable to evaluate error/uncertainty in predictions. In addition to the existing workflows, the infrastructure can support a wide variety of other technologically important applications as part of the data-driven materials design paradigm. The JARVIS datasets and tools are publicly available at the website: https://jarvis.nist.gov .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Development and Advanced Process Simulation for Additive Manufacturing with Fiber-Reinforced Thermoplastics (Final Technical Report)

DuPont, Local Motors, and the Composites Manufacturing and Simulation Center (CMSC) at Purdue worked together to advance process simulation capabilities for Additive Manufacturing (AM) and develop printable short fiber reinforced semi-crystalline engineering thermoplastics suitable for vehicle components. DuPont provided one candidate thermoplastic material, which was evaluated against the following three criteria: 1) Printability using extrusion additive manufacturing processes. 2) Suitability for vehicle applications based on mechanical coupon screening tests. 3) Predictive capability of process simulations. Indiana Manufacturing Institute provided modeling and simulation technology and access to Purdue’s lab scale AM printer. Local Motors provided mechanical property requirements and printed demonstration prototype components using their Thermwood Large Scale Additive Manufacturing (LSAM) printer. To meet the performance and cost requirements of the automotive industry, high fiber content semi-crystalline engineering materials are required. Printing experience and property data for a 50% by weight short glass fiber reinforced material were gathered. Further, material properties required for simulating the printing process with the glass-fiber reinforced material were characterized. Knowledge was acquired with regards to the processability and performance of this thermoplastic material. Predictions of part deformation and crystallinity level were verified against experimental observations. Any significant results other than verification of simulation, e.g., novel material identified here would be impactful.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulation and Modeling of Optical Properties of U, Th, Pb, and Co Nanoparticles of Interest to Nuclear Security Using Finite Element Analysis

In this work, the optical characteristics of uranium (U), lead (Pb), cobalt (Co), and thorium (Th) nanoparticles are fashioned and simulated employing the finite element analysis (FEA) approach concerning multiple particle sizes. Applying finite element analysis, it was found that the simulated absorption peaks of electronic excitations of nuclear nanoparticles are red-shifted from 365 nm to 555 nm for U; from 355 nm to 550 nm for Pb; from 415 nm to 610 nm for Co; and from 350 nm to 540 nm for Th, comparing expanding particle sizes from 60 nm to 100 nm (except for Co, which varied from 70 nm to 100 nm). The FEA-simulated optical band gap energies and far-field radiation patterns were also obtained for nuclear materials. The simulation approach in this research enables the prediction of optical properties and design of nuclear materials before manufacture for nuclear security applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Influence of fuel ash on the recoverability of copper from the spent material of chemical looping combustion

Chemical looping combustion (CLC) is an advanced energy production technology that inherently separates the CO 2 generated. The oxygen needed for fuel conversion is supplied by the oxygen carrier. When using copper-based oxygen carriers, oxygen carrier attrition is a challenge to the process economics, unless the attrited copper that is mixed with fuel ash in the spent material can be recovered by, for instance, acid leaching. However, fuel ash may transform the attrited copper into certain complexes of low leachability (e.g., CuAl 2 O 4 ), hindering the copper recovery. In this study, the copper content in the spent material was first evaluated. The effects of ash composition and copper content on the occurrence mode of copper in spent material were predicted by thermodynamic modeling. Redox cycles of simulated spent materials based on various coal ashes were performed to reveal the characteristics of CuAl 2 O 4 /CuFe 2 O 4 formation. The CuO-Al 2 O 3 reaction was investigated under the temperatures and solid residence times experienced in CLC. The materials were analyzed by X-ray diffraction. The results suggest that the fuel ash with a high Al/Ca ratio may lower the recoverability of copper when the temperature is above 925 °C. The slow formation of CuFe 2 O 4 is not expected to affect the copper recovery.

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

Python Urban Deployment Model (PyUDM) v1.0.0

The Python Urban Deployment Model (PyUDM) is a simulation tool used to investigate networks of static radiation detectors in urban environments. PyUDM simulates traffic, stationary NaI gamma-ray detectors, and moving radioactive sources on simulated vehicles. The analyzed output of the simulations contain valuable insights into the performance of different configurations of urban radiological detector configurations such as their ability to detect sources moving through the environment. To accurately simulate radioactive material moving through urban environments, PyUDM combines Monte Carlo simulation tools, publicly available map and traffic data, and measured gamma-ray spectra from urban environments.

Rofors, Emil↗