SEARCH · Engineering Papers

Results for “Material Point Method”

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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2024

High‐throughput combinatorial approach expedites the synthesis of a lead‐free relaxor ferroelectric system

Abstract Developing novel lead‐free ferroelectric materials is crucial for next‐generation microelectronic technologies that are energy efficient and environment friendly. However, materials discovery and property optimization are typically time‐consuming due to the limited throughput of traditional synthesis methods. In this work, we use a high‐throughput combinatorial synthesis approach to fabricate lead‐free ferroelectric superlattices and solid solutions of (Ba 0.7 Ca 0.3 )TiO 3 (BCT) and Ba(Zr 0.2 Ti 0.8 )O 3 (BZT) phases with continuous variation of composition and layer thickness. High‐resolution x‐ray diffraction (XRD) and analytical scanning transmission electron microscopy (STEM) demonstrate high film quality and well‐controlled compositional gradients. Ferroelectric and dielectric property measurements identify the “optimal property point” achieved at the composition of 48BZT–52BCT. Displacement vector maps reveal that ferroelectric domain sizes are tunable by varying {BCT–BZT} N superlattice geometry. This high‐throughput synthesis approach can be applied to many other material systems to expedite new materials discovery and properties optimization, allowing for the exploration of a large area of phase space within a single growth. image

36 MATERIALS SCIENCE↗

DOE OSTI2024

Convergent laser beam shapes: Unveiling the dynamics of Laser-induced elastic waves in composite materials

Overcoming the low signal-to-noise ratio (SNR) in laser ultrasonic testing of composite materials remains a significant challenge. Current efforts focus on enhancing SNR by inserting more energy into the material through temporal and/or spatial modulation of the laser beam. However, potential SNR improvements through wave convergence and wave energy manipulation have been overlooked. This paper addresses this gap by demonstrating the convergence of different wave types to a designated point and by showing the feasibility of directing absorbed laser energy into a specific wave type through spatial modulation of the laser beam. To achieve this, mathematical expressions for the convergent laser beams are derived. Various laser beam profiles are then introduced to the thermoelastic equations and solved using the finite element method. The sample under investigation is a transversely isotropic unidirectional carbon fiber reinforced plastic, characterized by anisotropic thermal expansion coefficients and thermal conductivities. Results reveal pronounced convergence of the intended wave type at the center due to laser beam shaping. This study showcases the ability to direct absorbed laser energy toward a specific wave type through spatial modulation of the laser beam and highlights the role of material anisotropy in energy focusing.

composite materials↗

DOE OSTI2024

STRUCTURAL MODELING TO SUPPORT POST-YIELD ACCEPTANCE CRITERIA FOR SPENT NUCLEAR FUEL CLADDING

Spent nuclear fuel (SNF) is evaluated for structural failure during storage and transportation scenarios. The U.S. Department of Energy’s Spent Fuel and Waste Science and Technology (SFWST) program has sponsored significant research in quantifying mechanical loads on SNF during storage and transportation scenarios using experimental and modeling methods. The SFWST program has also performed significant research on measuring the mechanical behavior of irradiated SNF as defueled cladding segments and cladding with fuel pellets to measure composite behavior. This paper considers some of the key material data from the Sibling Pin testing and uses structural modeling and analysis methods that have been informed by testing to consider post-yield acceptance criteria for SNF cladding structural analysis. Test data published by Oak Ridge National Laboratory (ORNL) and Pacific Northwest National Laboratory (PNNL) are the foundation for informing the material behavior of the models developed in this study. In particular, four-point bend (4PB) tests of fueled and defueled cladding segments provide significant information about the bending failure mode of SNF. ORNL’s 4PB test data is on fueled cladding segments, so the composite behavior of SNF is demonstrated. This paper describes PNNL’s coincident beam model that was developed to approximate the composite behavior of SNF. This paper also presents PNNL’s structural dynamic finite element models of a cask tip-over scenario, which is predicted to cause the strongest mechanical loads on SNF of all postulated storage and transportation scenarios. SNF bending loads predicted in the cask tip-over scenario and cladding acceptance criteria beyond yield are considered, with justification based on the Sibling Pin test data. ASME Boiler and Pressure Vessel code stress intensity limits are also considered. The ultimate goal of this work is to aid in the justification of structural acceptance criteria for SNF cladding beyond the cladding’s irradiated yield strength for use in structural analysis of all storage and transportation scenarios.

Klymyshyn, Nicholas A.↗

DOE OSTI2024

Building Stock Models for Embodied Carbon Emissions—A Review of a Nascent Field

Building stock modeling emerges as a critical tool in the strategic reduction of embodied carbon emissions, which is pivotal in reshaping the evolving construction sector. This review provides an overall view of modern methodologies in building stock modeling, homing in on the nuances of embodied carbon analysis in construction. Examining 23 seminal papers, our study delineates two primary modeling paradigms—top-down and bottom-up—each further compartmentalized into five innovative methods. This study points out the challenges of data scarcity and computational demands, advocating for methodological advancements that promise to refine the precision of building stock models. A groundbreaking trend in recent research is the incorporation of machine learning algorithms, which have demonstrated remarkable capacity, improving stock classification accuracy by 25% and urban material quantification by 40%. Furthermore, the application of remote sensing has revolutionized data acquisition, enhancing data richness by a factor of five. This review offers a critical examination of current practices and charts a course toward an environmentally prudent future. It underscores the transformative impact of building stock modeling in driving ecological stewardship in the construction industry, positioning it as a cornerstone in the quest for sustainability and its significant contribution toward the grand vision of an eco-efficient built environment.

Hu, Ming (ORCID:0000000325831161)↗

DOE OSTI2023

Defects go green: using defects in nanomaterials for renewable energy and environmental sustainability

Induction of point defects in nanomaterials can bestow upon them entirely new physics or augment their pre-existing physical properties, thereby expanding their potential use in green energy technology. Predicting structure-property relationships for defects a priori is challenging, and developing methods for precise control of defect type, density, or structural distribution during synthesis is an even more formidable task. Hence, tuning the defect structure to tailor nanomaterials for enhanced device performance remains an underutilized tool in materials design. We review here the state of nanomaterial design through the lens of computational prediction of defect properties for green energy technology, and synthesis methods to control defect formation for optimal performance. We illustrate the efficacy of defect-focused approaches for refining nanomaterial physics by describing several specific applications where these techniques hold potential. Most notably, we focus on quantum dots for reabsorption-free solar windows and net-zero emission buildings, oxide cathodes for high energy density lithium-ion batteries and electric vehicles, and transition metal dichalcogenides for electrocatalytic green hydrogen production and carbon-free fuels.

14 SOLAR ENERGY↗

DOE OSTI2025

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE↗

DOE OSTI2022

On the Prediction of Uniaxial Tensile Behavior Beyond the Yield Point of Wrought and Additively Manufactured Ti-6Al-4V

In this paper, phenomenological relationships are presented that permit the prediction of the plastic regime of stress–strain curves using a limited number of parameters. These relationships were obtained from both conventional (wrought + β annealed) and additively manufactured (i.e., “3D printed”) Ti-6Al-4V. Three different methods of additive manufacturing have been exploited to produce the materials, including large-volume electron beam additive manufacturing, large-volume laser hot wire additive manufacturing, and small-volume selective laser melting. The general fundamental expressions are independent not only of the additive manufacturing process, but also of a wide variety of post-deposition heat treatments, however the coefficients are specific to material states. Thus, this work demonstrates that it is possible to predict not only the ultimate tensile strength, but also the full true stress, true strain curves, if certain parameters of the material are known. In general, the prediction of ultimate tensile strength are within 5% of the experimentally measured values across all additive manufacturing variants and subsequent heat treatments. The absolute values of ultimate tensile strength range from ~ 910 MPa to ~ 1170 MPa for the single alloy Ti-6Al-4V. Data representing 113 explicit samples are included in this work.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗

DOE OSTI2021

Revisiting trends in the exchange current for hydrogen evolution

Nørskov and collaborators proposed a simple kinetic model to explain the volcano relation for the hydrogen evolution reaction on transition metal surfaces such that j 0 = k 0 f(ΔG H ) where j 0 is the exchange current density, f(ΔG H ) is a function of the hydrogen adsorption free energy ΔG H as computed from density functional theory, and k 0 is a universal rate constant. Herein, focusing on the hydrogen evolution reaction in acidic medium, we revisit the original experimental data and find that the fidelity of this kinetic model can be significantly improved by invoking metal-dependence on k 0 such that the logarithm of k 0 linearly depends on the absolute value of ΔG H . Here, we further confirm this relationship using additional experimental data points obtained from a critical review of the available literature. Our analyses show that the new model decreases the discrepancy between calculated and experimental exchange current density values by up to four orders of magnitude. Furthermore, we show the model can be further improved using machine learning and statistical inference methods that integrate additional material properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

SCALE 6.3 Validation: Radiation Shielding

Safe and reliable use of scientific and engineering computer codes requires validation for the types of applications in which they will be used. An example in the nuclear reactor engineering and licensing field is radiation transport employed in shielding analyses. The validity of computer codes for shielding applications is demonstrated in this report for SCALE version 6.3.0. Representative benchmarks corresponding to shielding analyses are selected for the validation study. Typical measurement results analyzed from these benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma dose rates, neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between measurement and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, the Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) radiation transport computer code with built-in variance reduction methods distributed with the SCALE computer code system agrees well with the measurement results. In selected benchmarks, MAVRIC is also compared to Monte Carlo N- Particle® (MCNP® ) 1 calculations. Both computer codes generally agree well within the estimated uncertainties. With the release of SCALE 6.3.0, Shift was integrated as an alternative transport solver in MAVRIC, denoted MAVRIC-Shift. Although the traditional MAVRIC using Monaco was used primarily in this validation study, many results have also been generated using MAVRIC-Shift. Agreement between MAVRIC-Monaco and MAVRIC-Shift is generally very good. The benchmarks presented in this report were obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook, the Shielding Integral Benchmark Archive & Database, and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to expand the shielding validation suite.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2025

A Modified Model Dielectric Function for Analyzing Optical Spectra of InGaN Nanofilms on Sapphire Substrates

Due to a lower InN bandgap energy Eg~0.7 eV, InxGa1−xN/Sapphire epifilms are considered valuable in the development of low-dimensional heterostructure-based photonic devices. Adjusting the composition x and thickness d in epitaxially grown films has offered many possibilities of light emission across a wide spectral range, from ultraviolet through visible into near-infrared regions. Optical properties have played important roles in making semiconductor materials useful in electro-optic applications. Despite the efforts to grow InxGa1−xN/Sapphire samples, no x- and d-dependent optical studies exist for ultrathin films. Many researchers have used computationally intensive methods to study the electronic band structures Ejk→, and subsequently derive optical properties. By including inter-band transitions at critical points from Ejk→, we have developed a semiempirical approach to comprehend the optical characteristics of InN, GaN and InxGa1−xN. Refractive indices of InxGa1−xN and sapphire substrate are meticulously integrated into a transfer matrix method to simulate d- and x-dependent reflectivity RE and transmission TE spectra of nanostructured InxGa1−xN/Sapphire epifilms. Analyses of RE and TE have offered accurate x-dependent shifts of energy gaps for InxGa1−xN (x = 0.5, 0.7) in excellent agreement with the experimental data.

Talwar, Devki N. (ORCID:0000000202723195)↗

DOE OSTI2022

Ultra-rapid, physics-based development pathway for reactor-relevant RF antenna materials

This paper presents a rapid, atomistically-informed, experimental development pathway for fusion reactor-relevant radio frequency (RF) antenna materials in the Cu-Cr-(Nb,Al,Zr) composition system, with the goal of improving upon GRCop-84. RF antennas in a tokamak fusion reactor will face a unique set of challenges as both structural and functional materials. The desired material must simultaneously achieve and maintain high electrical conductivity, high strength, high thermal conductivity, resist high temperatures, possess low nuclear activation, and incur low damage due to neutron bombardment. The GRCop-84 alloy serves as a starting point for iterative improvement, with the desire to reduce or eliminate Nb from the material to minimize nuclear activation. The rapid development pathway makes use of a multi-target combinatorial thick film sputtering process to produce full ternary phase diagrams on a Si wafer substrate. Transient grating spectroscopy (TGS), a laser-ultrasonic method, will determine spatially-varying thermo-elastic properties, while four terminal electrical conductivity measurements will map out the best per- forming regions of the sample for in-depth study at larger length scales. High energy proton and self-ion irradiation emulates the effects of neutron damage on the thermal/electric properties. With rapid turnaround time (∼days) in terms of mapping radiation damage-induced material property changes in the full ternary system, these techniques allow rapid iteration towards an optimal material, testing hundreds of nearby compositions in the time it took to test one. Focused testing of larger, single composition samples (produced in an arc furnace or by laser sintering) provides data on structural and high power RF properties, and validates our thick-film based workflow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum

Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.

58 GEOSCIENCES↗

DOE OSTI2025

Emulation of radiation transport in 3D stochastic media using 1D planar Monte Carlo stochastic media radiation transport algorithms

A subset of stochastic media radiation transport problems involves those in which radiation is incident on a thin slab of stochastic material. Particle tracking in 3D for such problems is expensive, and 1D planar models lack accuracy because they only allow the material to change in one dimension. Therefore, we propose dimensional emulation, which through a slight modification allows existing 1D planar geometry stochastic media radiation transport models to reproduce results from the equivalent 3D models by allowing the material to change in all three dimensions, reproducing the fidelity of the 3D model for the low computational cost of the 1D planar model. In this work, we apply dimensional emulation to three Monte Carlo stochastic media radiation transport models: Chord Length Sampling (CLS), the Local Realization Preserving method (LRP), and a variant of Conditional Point Sampling (CoPS). For a common Markovian benchmark set, the 3D emulation variants of these algorithms are numerically verified to reproduce the results of the 3D variants within statistics while running 1.3 to 2 times faster in the implementation within Sandia National Laboratories open-source research code PlaybookMC. The 3D emulation variants are also shown to yield a 72%–92% reduction in error for the thin slab problems in comparison to the 1D benchmark. As a result, the 3D emulation variant of CLS and CoPS-1 are shown to reproduce 3D CLS results that were used to approximate results for a 3D spherical inclusion geometry benchmark set.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Measurement of crack length in width tapered beam experiments

The width tapered beam method for measuring fracture toughness has been contemplated for use in measuring the adhesion of photovoltaic materials, is promoted by it being a viable method for use directly on modules rather than to engineered test specimens. However, precise determination of crack length is often difficult as failure often involves tendril formation, cavitation voids, and other plastic deformation at the crack front in addition to obfuscation by opaque materials. Because fracture energy varies as the square of crack length, imprecise crack tip length results in large measurement uncertainties. In this work, we develop a method to estimate the crack tip length using the relationship between the applied force and the displacement of the load frame. The crack tip is determined for every data point and can be defined from a calibration curve related to the amount of energy put into the beam, or from a theoretically idealized estimation. Once a calibration curve is set up and the equations programmed, the fracture energy as a function of position on the newly formed crack interface can be accurately estimated. The test additionally provides much easier analysis and insight into the joint failure through more precise knowledge of fracture toughness at each de-bonded area.

42 ENGINEERING↗

DOE OSTI2026

Integrated Process Model Utilization and Development for Inertial Fusion Energy

This CRADA between LLNS and Longview Fusion Energy Systems utilized a modernized LLNL Integrated Process Model (IPM) to assess technoeconomic viability and power plant configurations for Inertial Fusion Energy (IFE). During the collaboration period, a streamlined IPM model was produced by LLNL, which consolidated a combination of 3 models under one file. Additionally, obsolete parameters were removed and parameters available for trade space analysis were organized in a user-friendly fashion on the front interface. Gain scaling curves were discussed and analyzed in the framework of the IPM. Additional updates were made to account for cost scaling to today's dollars using an average inflation rate. Heat transfer material costing was updated using a ground up approach with public vendor data. Further levelized cost of electricity (LCOE) methods and models were discussed and reviewed for applicability to the IPM. Prior published literature on optical scaling was also discussed. Moreover, as part of this CRADA, the Participant developed a fusion technology development roadmap and point designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Integrated Process Model Utilization and Development for Inertial Fusion Energy

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Synthesis of Layered Gold Tellurides AuSbTe and Au 2 Te 3 and Their Semiconducting and Metallic Behavior

Previous studies on natural samples of pampaloite (AuSbTe) revealed the crystal structure of a potentially cleavable and/or exfoliable material, while studies on natural and synthetic montbrayite (Sb-containing Au 2 Te 3 ) claimed various chemical compositions for this low-symmetry compound. Few investigations of synthetic samples have been reported for both materials, leaving much of their chemical, thermal, and electronic characteristics unknown. Here, in this study, we investigate the stability, electronic properties, and synthesis of the gold antimony tellurides AuSbTe and Au 1.9 Sb 0.46 Te 2.64 (montbrayite). Differential thermal analysis and in situ powder X-ray diffraction revealed that AuSbTe is incongruently melting, while Au 1.9 Sb 0.46 Te 2.64 is congruently melting. Calculations of the band structures and four-point resistivity measurements showed that AuSbTe is a semiconductor and Au 1.9 Sb 0.46 Te 2.64 a metal. Various synthesis attempts confirmed the limited stable chemical composition of Au 1.9 Sb 0.46 Te 2.64 , identified successful methods to synthesize both compounds, and highlighted the challenges associated with single-crystal synthesis of AuSbTe.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗