Engineering Papers⌕ Search

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.

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

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Non-Diffusive Volume Advection with A High Order Interface Reconstruction Method

We show that non-diffusive volume advection in two-dimensions is achieved with several benchmark problems using a newly developed high-order volume of fluids (VOF) interface reconstruction method. (1) A new VOF interface reconstruction method using circular/corner facets (linear facets are a degenerate case of arcs). We create a circular interface facet in each mixed zone by matching neighbor volume with a hybrid Newton’s-bisection method and the local solution is final. In the general case, the new VOF interface reconstruction has 3rd order accuracy and can be easily made seamless. The new method addresses intrinsic issues with Young’s method such as gaps between interface facets in the case of a curved interface, and inability to define curvature nor identify corners. (2) A non-diffusive volume advection scheme. In an ALE advection step, a well-defined interface can be carried over through a Lagrange step and used to compute volume distribution into a relaxed mesh. Then, an interface reconstruction step is performed to redefine the interface in the relaxed mesh. We must point out that the interface carried over is also a solution of interface reconstruction because all the volume fractions in the relaxed mesh are naturally matched. We provide an interface tracking method compatible with our reconstruction scheme, where it is granted to use the prior info as an initial guess to capture sub-mesh resolution features. As a result, we are able to treat multiple facets inside a single mixed cell and obtain highly accurate, non-diffusive solution for advection problems with rather coarse meshes. We show our solutions for two-dimensional incompressible flows with two materials with a) the X + O diagonal translation; b) the Zalesak rotational test; and c) the single vortex spiral test.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Membrane Development for CO2 Capture from Steel Manufacturing

The U.S. government is targeting a net-zero carbon-emission economy by 2050, offering an exciting opportunity for membrane-based CO2 capture from various industrial point sources. Given that industrial flue gas has low CO2 partial pressures and high volumetric flow rates, high-permeance membranes are needed to make membrane technology economically viable for large-scale deployment. Thin film composite (TFC) membranes are necessary for this implementation because they can provide high permeance by forming a thin selective layer on top of a porous support. This presentation will report the rational design and fabrication of NETL’s highly permeable non-aging TFCs achieved by: synthesizing a high-performance rubbery selective material; developing a high-porosity membrane support; optimizing coating methods to assemble the two materials into scalable membranes; and scaling up membrane supports and TFCs via a roll-to-roll process. The novel rubbery selective material shows mixed-gas CO2 permeability of 930 Barrer and CO2/N2 selectivity of 44, exceeding the 2008 Robeson upper bound. The resulting TFCs yield remarkably high CO2 permeance of 4,500 GPU and CO2/N2 selectivity of 34 at 23ºC. Moreover, the TFCs exhibit not only excellent performance stability (or non-aging behavior) for 1,000 hours in the lab, but also maintain their separation properties in a 700-hour field test at the U.S. DOE’s National Carbon Capture Center using real humid flue gas.

Tran, Thien↗
DOE OSTI2025

Effect of sintering temperature on adhesion of spray-on piezoelectric transducers

Conventionally sol-gel spray-on transducers require a high-temperature (> 700 ◦C) sintering process; however, this process can affect the microstructure of the substrate material. For mechanical elbows and valves utilized for fluid transport in the energy sector, the components are designed to have a specific microstructure, and deviations from these specifications can create weak points in the system. For this reason it is important to investigate how the temperature of the deposition process affects the substrate. This paper investigates the effect of high-temperature and low-temperature (< 150 ◦C) processing conditions on the surface composition of the substrate. Furthermore, the resultant transducers from high- and low-temperature fabrication processes are compared to determine if a low-temperature processing method is feasible. For these studies a sol-gel spray-on process is employed to deposit piezoelectric ceramics onto a stainless-steel 316L substrate. Energy-dispersive X-ray spectroscopy is utilized to determine the composition of the substrate surface before and after transducer deposition. Results indicate that the high-temperature processing conditions may alter the surface composition of the metal due to a diffusion of the metal into the ceramic, which results in a metal surface that is bonded to the ceramic. Furthermore, it is shown that low-temperature processing of spray-on transducers is a viable method for transducer fabrication where the resultant transducers meet the industry minimum requirement of 30 dB signalto-noise ratio. In parallel simulation calculations, finite-element method (FEM) studies were performed to model the adhesive strength of the low-temperature processed transducer to the substrate surface. Comparisons between the simulations and experiments suggest that the bond strength is much greater than the commercial gel bonds and closer to hardened epoxy glue bonds. These results indicate that spray-on transducers fabricated under lowtemperature processing conditions are a viable solution for leave-in-place monitoring of structures.

M. Sinding, Kyle↗
DOE OSTI2024

Residual Stress in Cold Spray SS304L Measured Via Neutron Diffraction and Comparison of Analytical Models to Predict the Residual Stress

Here, this study employs neutron diffraction to investigate the relationship between residual stress and coating thickness in cold sprayed 304L austenitic stainless steel. Results show that shot peening predominantly impacts the residual stress profile, leading to substantial in-plane compressive force. The impact of laser heating, a widely used method to alter cold spray's microstructural properties, on the coating's residual stress is also analyzed. The findings indicate that the maximum compressive residual stress in the in-plane component is mainly independent of coating thickness, which suggests that the material properties determine the maximum residual stress. The cold sprayed deposits possessed compressive, nearly biaxial strain and stresses. After laser heating, these stresses were replaced by tensile residual stresses. Two analytical models, the Tsui and Clyne and the Boruah models, for predicting residual stresses are also evaluated, and both models provide reasonable fits to the experimental data. At this point, the deviations between the experimental results and the models are principally caused by the inability of the current models to address plastic deformation and relaxation, and the residual stresses generated by thermal gradients.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Deciphering Fingerprints of Stellar Nucleosynthesis Through Nuclear Reaction Rate Measurements and Isotopic Analyses of Stardust

Microscopic stardust grains, found in trace amounts in primitive meteorites, afford a unique opportunity to study stellar nucleosynthesis in the laboratory by measuring their heavy element isotopic compositions. Applying LLNL’s state-of-the art material characterization capabilities to these micrometer-sized grains yields large new datasets with unprecedented precision. However, the astrophysical models necessary to interpret the isotopic record require improved nuclear physics data to constrain stellar nucleosynthesis conditions. For example, the branch point in the s-process path at 95 Zr controls the relative amounts of heavier nuclides but is poorly understood. We conducted a coordinated effort across disciplines to combine a new determination of the neutron capture cross section of 95 Zr with the laboratory analysis of hundreds of stardust grains to yield a view into stellar interiors. We developed new methods to quantify multiple isotopic systems simultaneously (e.g., Zr, Ba, and W) in stardust grains and improved the determination of the neutron capture cross section of unstable 95 Zr via the surrogate reaction method. This interdisciplinary approach is allowing us to fill gaps in our understanding of the environments in which elements are produced while enhancing LLNL’s experimental and theoretical capabilities relevant to nuclear threat reduction missions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Deciphering Fingerprints of Stellar Nucleosynthesis Through Nuclear Reaction Rate Measurements and Isotopic Analyses of Stardust

Microscopic stardust grains, found in trace amounts in primitive meteorites, afford a unique opportunity to study stellar nucleosynthesis in the laboratory by measuring their heavy element isotopic compositions. Applying LLNL’s state-of-the art material characterization capabilities to these micrometer-sized grains yields large new datasets with unprecedented precision. However, the astrophysical models necessary to interpret the isotopic record require improved nuclear physics data to constrain stellar nucleosynthesis conditions. For example, the branch point in the s-process path at 95 Zr controls the relative amounts of heavier nuclides but is poorly understood. We conducted a coordinated effort across disciplines to combine a new determination of the neutron capture cross section of 95 Zr with the laboratory analysis of hundreds of stardust grains to yield a view into stellar interiors. We developed new methods to quantify multiple isotopic systems simultaneously (e.g., Zr, Ba, and W) in stardust grains and improved the determination of the neutron capture cross section of unstable 95 Zr via the surrogate reaction method. This interdisciplinary approach is allowing us to fill gaps in our understanding of the environments in which elements are produced while enhancing LLNL’s experimental and theoretical capabilities relevant to nuclear threat reduction missions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Ground State Properties of the Wide Band Gap Semiconductor Beryllium Sulfide (BeS)

We report the results from self-consistent calculations of electronic, transport, and bulk properties of beryllium sulfide (BeS) in the zinc-blende phase, and employed an ab-initio local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO). We obtained the ground state properties of zb-BeS with the Bagayoko, Zhao, and Williams (BZW) computational method, as enhanced by Ekuma and Franklin (BZW-EF). Our findings include the electronic energy bands, the total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated band structure clearly shows that zb-BeS has an indirect energy band gap of 5.436 eV, from Γ to a point between Γ and X, for an experimental lattice constant of 4.863 Å. This is in excellent agreement with the experiment, unlike the findings of more than 15 previous density functional theory (DFT) calculations that did not perform the generalized minimization of the energy functional, required by the second DFT theorem, which is inherent to the implementation of our BZW-EF method.

36 MATERIALS SCIENCE↗
DOE OSTI2021

MOSAIC: An Effective FFT-based Numerical Method to Assess Aging Properties of Concrete

As the nuclear fleet in the United States ages and subsequent license renewal applications grow, the prediction of concrete durability at extended operation becomes more important. To address this issue, a Fast-Fourier Transform (FFT) method is utilized to simulate aging-related degradation of concrete within the Microstructure Oriented Scientific Analysis of Irradiated Concrete (MOSAIC) software. MOSAIC utilizes compositional phase maps to simulate damage from radiation-induced volumetric expansion (RIVE), applied force, creep, and thermal expansion. This compositional detail allows each mineral in the microstructure to be assigned specific material properties, allowing the simulation to be as accurate and representative as possible. The principal goal of MOSAIC is to simulate the effects of nonlinear aging mechanisms occurring in nuclear concrete on the macroscopic mechanical properties, using only the aggregate microstructure compositional information as a starting point. Here, several realistic example simulations are shown to demonstrate the utility and uniqueness of the MOSAIC software.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Uncertainty quantification for neural network potential foundation models

Abstract For neural network potentials (NNPs) to gain widespread use, researchers must be able to trust model outputs. However, the blackbox nature of neural networks and their inherent stochasticity are often deterrents, especially for foundation models trained over broad swaths of chemical space. Uncertainty information provided at the time of prediction can help reduce aversion to NNPs. In this work, we detail two uncertainty quantification (UQ) methods. Readout ensembling, by finetuning the readout layers of an ensemble of foundation models, provides information about model uncertainty, while quantile regression, by replacing point predictions with distributional predictions, provides information about uncertainty within the underlying training data. We demonstrate our approach with the MACE-MP-0 model, applying UQ to the foundation model and a series of finetuned models. The uncertainties produced by the readout ensemble and quantile methods are demonstrated to be distinct measures by which the quality of the NNP output can be judged.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stable diagonal stripes in the t-J model at $\bar{n}$ h = 1/8 doping from fPEPS calculations

We investigate the two-dimensional t–J model at a hole doping of $\bar{n}$ h= 1 / 8 using recently developed high accuracy fermionic projected entangled pair states method. By applying the stochastic gradient descent method combined with the Monte Carlo sampling technique, we obtain the ground state hole energy Ehole = -1.621 for J/t = 0.4. We show that the ground state has stable diagonal stripes instead of vertical stripes with a width of 4 unit cells, and stripe filling ρ l = 0.5. We further show that the long-range superconductivity order is suppressed at this point.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Adaptively driven X-ray diffraction guided by machine learning for autonomous phase identification

Machine learning (ML) has become a valuable tool to assist and improve materials characterization, enabling automated interpretation of experimental results with techniques such as X-ray diffraction (XRD) and electron microscopy. Because ML models are fast once trained, there is a key opportunity to bring interpretation in-line with experiments and make on-the-fly decisions to achieve optimal measurement effectiveness, which creates broad opportunities for rapid learning and information extraction from experiments. Here, we demonstrate such a capability with the development of autonomous and adaptive XRD. By coupling an ML algorithm with a physical diffractometer, this method integrates diffraction and analysis such that early experimental information is leveraged to steer measurements toward features that improve the confidence of a model trained to identify crystalline phases. We validate the effectiveness of an adaptive approach by showing that ML-driven XRD can accurately detect trace amounts of materials in multi-phase mixtures with short measurement times. The improved speed of phase detection also enables in situ identification of short-lived intermediate phases formed during solid-state reactions using a standard in-house diffractometer. Our findings showcase the advantages of in-line ML for materials characterization and point to the possibility of more general approaches for adaptive experimentation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Charged point defect benchmark of Hybrid and GGA-PBE

Data used by the publication "High-throughput calculations of charged point defect properties with semi-local density functional theory - performance benchmarks for materials screening applications." This work presented an in-depth benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We considered three different a-posteriori correction sets for semi-local calculations, implemented in a fully automated workflow, and consider the qualitative and quantitative differences for four different categories of defect information: thermodynamic transition levels, formation energies, fermi levels, and dopability limits. We highlighted the type of qualitative information about point defect properties that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy that can be achieved by these approaches.

Broberg, Danny↗
DOE OSTI2026

Grid Operator Analytics and Assessment Tools for Inverter- Based Resources Dominated Grid (GOAAT-IBR) Project Update

This presentation provides an update on the OPTIMA GOAAT project, with emphasis on the cloud-native data platform developed in-house to ingest, manage, and operationalize high-resolution power system data. Since our last NASPI presentation, accessible via OSTI ID #2671437, the project team advanced the design and deployment of a scalable architecture capable of handling both synchronized and non-synchronized streams, including PMU, point-on-wave (POW), COMTRADE, and SCADA data. These materials review the project status, recent progress, and key lessons learned. The core of the presentation examines the architecture and engineering of our cloud-native ingestion and data management platform. We then explain how pipelines were designed to collect, normalize, time-align, store, and serve heterogeneous data at scale. We will discuss design choices such as data models, streaming versus batch ingestion, storage tiers, and interoperability with analytics applications. Practical experiences with cloud-native technologies were shared during the event, including benefits, limitations, and integration challenges in a utility environment, along with methods used to improve performance, reduce latency, and optimize resource usage. The presentation also showcases user interface designs and visualization tools that convert raw measurements and analytics results into intuitive, actionable insights for operators and engineers. During the presentation examples were provided demonstrating how visualization, event views, and summarized analytics enhance situational awareness and support operational decision-making. These use cases illustrate how a well-designed data infrastructure can bridge the gap between high-volume measurements and practical grid operations.

Aminifar, Farrokh↗
DOE OSTI2022

Estimates of Quantum Tunneling Effects for Hydrogen Diffusion in PuO 2

We detail the estimation of activation energies and quantum nuclear vibrational tunneling effects for hydrogen diffusion in PuO 2 based on Density Functional Theory calculations and a quantum double well approximation. We find that results are relatively insensitive to choice of exchange correlation functional. In addition, the representation of spin in the system and use of an extended Hubbard U correction has only a small effect on hydrogen point defect formation energies when the PuO 2 lattice is held fixed at the experimental density. We then compute approximate activation energies for transitions between hydrogen interstitial sites seeded by a semi-empirical quantum model and determine the quantum tunneling enhancement relative to classical kinetic rates. Our model indicates that diffusion rates in H/PuO 2 systems could be enhanced by more than one order of magnitude at ambient conditions and that these effects persist at high temperature. The method we propose here can be used as a fast screening tool for assessing possible quantum nuclear vibrational effects in any number of condensed phase materials and surfaces, where hydrogen hopping tends to follow well defined minimum energy pathways.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Decoding defect statistics from diffractograms via machine learning

Abstract Diffraction techniques can powerfully and nondestructively probe materials while maintaining high resolution in both space and time. Unfortunately, these characterizations have been limited and sometimes even erroneous due to the difficulty of decoding the desired material information from features of the diffractograms. Currently, these features are identified non-comprehensively via human intuition, so the resulting models can only predict a subset of the available structural information. In the present work we show (i) how to compute machine-identified features that fully summarize a diffractogram and (ii) how to employ machine learning to reliably connect these features to an expanded set of structural statistics. To exemplify this framework, we assessed virtual electron diffractograms generated from atomistic simulations of irradiated copper. When based on machine-identified features rather than human-identified features, our machine-learning model not only predicted one-point statistics (i.e. density) but also a two-point statistic (i.e. spatial distribution) of the defect population. Hence, this work demonstrates that machine-learning models that input machine-identified features significantly advance the state of the art for accurately and robustly decoding diffractograms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dynamic Tensile Characterization of Thin-Sheet Brittle Metallic Materials

Refractory metals are favorable materials in applications where high strength and ductility are needed at elevated temperatures. In some cases, operating temperatures may be near the melting point of the material. However, as temperature drops, refractory metals typically undergo a significant mechanical response change - ductile-to-brittle transition. These materials may be subjected to high strain rate loading at an ambient temperature state, such as an impact or crash. Knowledge of the high rate material properties are essential for design as well as simulation of impact events. The high rate stress-strain behavior of brittle metallic materials at ambient temperature is rarely studied because of experimental challenges, particularly when failure is involved. Failure typically occurs within the non-gage section of the material, which invalidates any collected stress-strain information. Here, a method to determine a specimen geometry which will produce failures in the gage section is presented. Pure tungsten in thin-sheet form was used as a trial material to select a specimen geometry for high rate Kolsky tension bar experiments. A finite element simulation was conducted to derive a strain correction for more accurate results. The room temperature stress-strain behavior of pure tungsten at a strain rate of 24 s –1 is presented. The outcome of this experimental technique can be applied to other brittle materials for dynamic tensile characterization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Combinatorial Exploration and Mapping of Phase Transformation in a Ni–Ti–Co Thin Film Library

Combinatorial synthesis and high-throughput characterization of a Ni–Ti–Co thin film materials library are reported for exploration of reversible martensitic transformation. The library was prepared by magnetron co-sputtering, annealed in vacuum at 500 °C without atmospheric exposure, and evaluated for shape memory behavior as an indicator of transformation. Composition, structure, and transformation behavior of the 177 pads in the library were characterized using high-throughput wavelength dispersive spectroscopy (WDS), X-ray photoelectron spectroscopy (XPS), X-ray diffraction (XRD), and four-point probe temperature-dependent resistance (R(T)) measurements. A new, expanded composition space having phase transformation with low thermal hysteresis and Co > 10 at. % is found. Unsupervised machine learning methods of hierarchical clustering were employed to streamline data processing of the large XRD and XPS data sets. Through cluster analysis of XRD data, we identified and mapped the constituent structural phases. Finally, composition–structure–property maps for the ternary system are made to correlate the functional properties to the local microstructure and composition of the Ni–Ti–Co thin film library.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗