Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solid modeling”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Attosecond synchronization of extreme ultraviolet high harmonics from crystals

The interaction of strong near-infrared (NIR) laser pulses with wide-bandgap dielectrics produces high harmonics in the extreme ultraviolet (XUV) wavelength range. These observations have opened up the possibility of attosecond metrology in solids, which would benefit from a precise measurement of the emission times of individual harmonics with respect to the NIR laser field. In this work, we show that, when high-harmonics are detected from the input surface of a magnesium oxide crystal, a bichromatic probing of the XUV emission shows a clear synchronization largely consistent with a semiclassical model of electron–hole recollisions in bulk solids. On the other hand, the bichromatic spectrogram of harmonics originating from the exit surface of the 200 μm-thick crystal is strongly modified, indicating the influence of laser field distortions during propagation. Our tracking of sub-cycle electron and hole re-collisions at XUV energies is relevant to the development of solid-state sources of attosecond pulses.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Probing atomic physics at ultrahigh pressure using laser-driven implosions

Abstract Spectroscopic measurements of dense plasmas at billions of atmospheres provide tests to our fundamental understanding of how matter behaves at extreme conditions. Developing reliable atomic physics models at these conditions, benchmarked by experimental data, is crucial to an improved understanding of radiation transport in both stars and inertial fusion targets. However, detailed spectroscopic measurements at these conditions are rare, and traditional collisional-radiative equilibrium models, based on isolated-atom calculations and ad hoc continuum lowering models, have proved questionable at and beyond solid density. Here we report time-integrated and time-resolved x-ray spectroscopy measurements at several billion atmospheres using laser-driven implosions of Cu-doped targets. We use the imploding shell and its hot core at stagnation to probe the spectral changes of Cu-doped witness layer. These measurements indicate the necessity and viability of modeling dense plasmas with self-consistent methods like density-functional theory, which impact the accuracy of radiation transport simulations used to describe stellar evolution and the design of inertial fusion targets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Parallel three-dimensional simulations of quasi-static elastoplastic solids

Hypo-elastoplasticity is a flexible framework for modeling the mechanics of many hard materials under small elastic deformation and large plastic deformation. Under typical loading rates, most laboratory tests of these materials happen in the quasi-static limit, but there are few existing numerical methods tailor-made for this physical regime. Here, we extend to three dimensions a recent projection method for simulating quasi-static hypo-elastoplastic materials. The method is based on a mathematical correspondence to the incompressible Navier–Stokes equations, where the projection method of Chorin (1968) is an established numerical technique. We develop and utilize a three-dimensional parallel geometric multigrid solver employed to solve a linear system for the quasi-static projection. Our method is tested through simulation of three-dimensional shear band nucleation and growth, a precursor to failure in many materials. As an example system, we employ a physical model of a bulk metallic glass based on the shear transformation zone theory, but the method can be applied to any elastoplasticity model. We consider several examples of three-dimensional shear banding, and examine shear band formation in physically realistic materials with heterogeneous initial conditions under both simple shear deformation and boundary conditions inspired by friction welding.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A wide-range thermodynamically consistent constitutive model for dynamic loading of partially saturated porous media

SUMMARY This paper presents a novel thermodynamically consistent constitutive model for partially saturated porous rocks across a wide range of conditions. The material states generated behind the shock wave from an explosive source can vary significantly, ranging from crushed and melted rock near the source to a poroelastic medium in the far field. In the model, rock strength is determined by the effective pressure, which is calculated using two independent equations of state: one for the solid rock and another for the pore fluid. The model accounts for shock-induced liquefaction resulting from fluid pressure buildup in the pore spaces near the explosive source. Simultaneously, it describes the increase in wave propagation speed due to elastic pore contraction in both dry and partially saturated rocks. This model is applied to investigate how fluid saturation affects the amplitude and shape of the generated waves, as well as the residual stress surrounding the cavity formed by spherical explosions.

Vorobiev, Oleg [Lawrence Livermore National Labora↗
DOE OSTI2024

Integrated modeling of boron powder injection for real-time plasma-facing component conditioning

An integrated modeling framework for investigating the application of solid boron (B) powder injection for real-time surface conditioning of plasma-facing components (PFCs) in tokamak environments is presented. Utilizing the DIII-D impurity powder dropper (IPD) setup, this study simulates B powder injection scenarios ranging from milligrams to tens of milligrams per second, corresponding to boron flux rates of 10 20 –10 21 B/s in standard L-mode conditions. The comprehensive modeling approach combines EMC3-EIRENE for simulating the deuterium plasma background and the Dust Injection Simulator (DIS) for the ablation and transport of the boron powder particles. EMC3 trace impurity fluid modeling results show substantial boron transport to the inboard lower divertor, predominantly influenced by the main ion plasma flow. The dependency on powder particle size (5-250 µm) was found to be insignificant for the scenario considered. The effects of erosion and redeposition were considered to reconcile the discrepancies with experimental observations, which saw substantial deposition on the outer divertor plasma-facing components. For this purpose, the WallDYN3D code was updated to include boron sources within the plasma domain and integrated into the modeling framework. The mixed-material migration modeling shows evolving boron deposition patterns, suggesting the formation of mixed B-C layers or predominantly B coverage depending on the powder mass flow rate. While the modeling outcomes at lower B injection rates tend to align with DIII-D experimental observations, the prediction of near-pure boron layers at higher rates has yet to be experimentally verified in the carbon environment of the DIII-D tokamak. The extensive reach of boron layers found in the modeling suggests the need for modeling that encompasses the entire wall geometry for more accurate experimental correlations. This integrated approach sets a precedent for analyzing and applying real-time in-situ boron coating techniques in advanced tokamak scenarios, potentially extendable to ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comparison of methods to derive infrared optical constants for organic materials for optical modeling

Reflectance spectroscopy, especially at infrared wavelengths, is often used for contact, standoff, and remote sensing of solid materials. The reflectance spectra of solids, however, are complex, relying on many factors, even for the same material. Such phenomena can be modeled if the optical constants as a function of wavelength or wavenumber, n(?) and k(?), are known. Methods to measure the optical constants of solids, however, are challenging, particularly for powdered materials. For powdered materials, a pressed pellet of the neat material is often used when a crystalline specimen is not available. In this work, three techniques, including ellipsometry, single-angle reflectance and KBr transmission spectroscopy, are applied to an organic material, acetaminophen, for comparison, and the effects on the modeled spectra are shown.

ellipsometry, infrared (IR) spectroscopy, FTIR, hy↗
DOE OSTI2021

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗
DOE OSTI2025

SOC Microstructural Property Estimator

This pre-trained ML model is a tool that uses basic compositional parameters for porous solid oxide cell (SOC) electrodes - the phase fractions and mean particle/pore diameters – as inputs and uses them to estimate additional electrochemical performance parameters: active (i.e., connected) TPB density, all tortuosity factors, and phase pair specific interfacial areas. The electrode is assumed to be composed of two solid phases and a pore phase. The property calculations are performed using neural network regression models trained on a large bank of synthetic electrode microstructural data that NETL has generated using the program DREAM3D (that bank is also hosted on EDX: https://edx.netl.doe.gov/dataset/soc-synthetic-microstructure-bank). This means the generated parameters are based on training from actual measured properties from 3D microstructures, not estimated from geometric simplifications. This tool was developed and is intended to replace percolation theory calculations in models that use hypothetical electrode properties. An example use case would be running SOC performance simulations across a parametric sweep of electrode designs (e.g., varying phase fractions and particle sizes) and assessing how it impacts the electrochemical performance of the SOC. Within the parameter space of the training data (statistics of that parameter space is provided in the readme file), this model achieves sub-5% mean absolute percent errors, an order of magnitude less error than percolation theory across the same parameter space. However, be aware that this tool was developed with parametric simulations in mind, and users are encouraged to assess accuracy for their own specific use case rather than taking accuracy metrics at face value. More info, including a usage guide, is in the included readme file. This tool should be cited with the DOI number provided.

Electrode Microstructure↗
DOE OSTI2025

Pressure–Modulated Luminescence Enhancement and Quenching in a Hydrogen–Bonded Organic Framework

Light emission in the solid state is central for illumination, sensing, and imaging applications. Unlike luminescence in dilute solutions, where the excited states are unimolecular in nature, intermolecular interaction plays a significant role in the quantum yield of solid-state luminophores, manifested as competing aggregation-caused quenching (ACQ) and aggregation-induced enhancement (AIE). Both effects are extensively studied in various systems; however, it remains unclear how their competition depends on molecular conformation and intermolecular stacking. Here the direct observation of pressure-modulated AIE-ACQ competition in a crystalline hydrogen-bonded organic framework (HOF) is reported. Using in situ spectroscopies and computational modeling, the intramolecular vibration and intermolecular π–π stacking directly responsible for the non-radiative decay of the excited state are identified. The extent of these two contributions is modulated by hydrostatic pressure and guest molecules in the HOF pores. Furthermore, this work demonstrates a physically neat model system to understand and control solid-state luminescence, and a potential material platform for piezoluminescent sensing.

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

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗
DOE OSTI2024

Stoichiometric Lanthanide Compounds with Diglycolamides: A Synthetic Approach toward Understanding Rare-Earth Speciation in Solution

A fundamental understanding of coordination chemistry across the lanthanide series is essential for explaining the chemical behavior of rare-earth metals in complex liquid–liquid extraction processes. Probing the exact bonding between the extractant and the metal is sometimes done through the synthesis of solid-state compounds that can serve as models for metal speciation in solution. In the case of diglycolamide (DGA), a commonly used neutral diamide extractant, extensive studies identify the stepwise formation of 1:1 [Ln(DGA)(H 2 O) 6 ] 3+ , 1:2 [Ln(DGA) 2 (H 2 O) 3 ] 3+ , and 1:3 [Ln(DGA) 3 ] 3+ complexes in solution. The crystallographic reports, however, exclusively show a 1:3 [Ln(DGA) 3 ] 3+ moiety in the solid state, while the structures of 1:1 and 1:2 complexes are yet to be isolated and comprehensively studied. In this work, we report the synthesis and characterization of three new families of stoichiometric N,N,N′,N’ -tetramethyldiglycolamide (TMDGA) compounds: [Ln(TMDGA)(H 2 O) 5 Cl]Cl 2 ·2H 2 O, [Ln(TMDGA) 2 (H 2 O) 3 ]Cl 3 ·3H 2 O, and [Ln(TMDGA) 3 ]Cl 3 ·7H 2 O, where Ln = Nd, Eu, Gd, with Ln:TMDGA ratios of 1:1, 1:2, and 1:3, respectively. Further, the compounds have been analyzed using vibrational spectroscopy (both Raman and FT-IR), as well as variable temperature fluorescence spectroscopy. A spectrophotometric titration of Eu(III) was performed to confirm the presence of the [Eu(TMDGA) n ] 3+ ( n = 1, 2, and 3) species in solution and to compare the individual solid-state emission spectra to their respective analogues in solution.

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

The Microscopic Mechanism of Lithiation and Delithiation in the Ag/C Buffer Layer for Anode‐Free Solid‐State Batteries

Abstract Lithium metal solid‐state batteries (LMSSBs) have demonstrated their high energy density and cycling performance at high current densities in an anode‐free architecture, featuring a thin Ag/C composite buffer layer (BL) between the current collector (CC) and the solid electrolyte (SE). This study explains the microscopic mechanism of the Ag/C BL by using first‐principles atomistic and continuum modeling. It is shown that Ag effectively acts as a homogeneous solid‐solution beyond AgLi 2.32 and maintains a positive potential even at AgLi 25 during lithiation. Key factors underlying the working of the Ag/C BL include lower interfacial resistance at the BL/CC interface than at the BL/SE interface, leading to predominant Li deposition on BL/CC, and substantial Ag–Li volume expansion during lithiation. This, combined with stronger BL/SE adhesion, causes BL/SE separation and Ag–Li extrusion toward the CC side. During delithiation, Ag re‐precipitates as nanoparticles uniformly on the CC, with its positive lithiation potential homogenizing Li currents in subsequent cycles. Other metals are less effective due to their relatively large overpotential, premature lithiation termination, and limited volume expansions hindering movement toward the CC. The study aids the BL design, focusing on metal choice and optimization material and microstructural properties, such as the Li‐ion conductivity and interfacial resistance.

25 ENERGY STORAGE↗
DOE OSTI2023

Direct imaging of shock wave splitting in diamond at Mbar pressure

Understanding the behavior of matter at extreme pressures of the order of a megabar (Mbar) is essential to gain insight into various physical phenomena at macroscales—the formation of planets, young stars, and the cores of super-Earths, and at microscales—damage to ceramic materials and high-pressure plastic transformation and phase transitions in solids. Under dynamic compression of solids up to Mbar pressures, even a solid with high strength exhibits plastic properties, causing the induced shock wave to split in two: an elastic precursor and a plastic shock wave. This phenomenon is described by theoretical models based on indirect measurements of material response. The advent of x-ray free-electron lasers (XFELs) has made it possible to use their ultrashort pulses for direct observations of the propagation of shock waves in solid materials by the method of phase-contrast radiography. However, there is still a lack of comprehensive data for verification of theoretical models of different solids. Here, we present the results of an experiment in which the evolution of the coupled elastic–plastic wave structure in diamond was directly observed and studied with submicrometer spatial resolution, using the unique capabilities of the x-ray free-electron laser (XFEL). The direct measurements allowed, for the first time, the fitting and validation of the 2D failure model for diamond in the range of several Mbar. Our experimental approach opens new possibilities for the direct verification and construction of equations of state of matter in the ultra-high-stress range, which are relevant to solving a variety of problems in high-energy-density physics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Continuous recovery of phosphoric acid and Rare-Earths containing particles from phosphoric acid sludge using a decanter centrifuge

Recovery of rare earth elements (REEs) from various industrial and natural streams currently draws significant attention in efforts to meet the demands of the manufacturing industry. Among many industrial byproducts and waste streams, phosphoric acid sludge could be one of the most economically feasible resources for the recovery of REEs because solid particles in the sludge contain relatively concentrated REEs, up to 3,000 ppm, while the liquid component of the sludge is valuable phosphoric acid (P 2 O 5 ) that can be recovered and returned to the main product. Due to high viscosity and large solids content (e.g., 30–40 %), however, this byproduct stream requires multistep separation and purification processes. In this study, a single-step process involving a continuous-flow decanter centrifuge (CFDC) was employed to investigate its feasibility for continuous solid/liquid separation from real phosphoric acid sludge. High centrifugal forces generated from up to 1500 G gravity acceleration separate solid particles from the sludge, generating a liquid-rich stream and a solids-rich stream at the exit of the CFDC. A single pass of phosphoric-acid sludge through the CFDC yielded 95 % liquid recovery and 90 % recovery of REEs-containing solids from 20 to 34 wt% solids-containing sludge. Finally, a reduced order model developed for the CFDC operation showed good agreement with experimental data, and preliminary technoeconomic analysis revealed potential process feasibility.

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

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗