Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heterogeneous materials”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Microstructure modeling of nuclear structural materials: Recent progress and future directions

Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spatial Control of Substitutional Dopants in Hexagonal Monolayer WS 2 : The Effect of Edge Termination

The ability to control the density and spatial distribution of substitutional dopants in semiconductors is crucial for achieving desired physicochemical properties. Substitutional doping with adjustable doping levels has been previously demonstrated in 2D transition metal dichalcogenides (TMDs); however, the spatial control of dopant distribution remains an open field. In this work, edge termination is demonstrated as an important characteristic of 2D TMD monocrystals that affects the distribution of substitutional dopants. Particularly, in chemical vapor deposition (CVD)-grown monolayer WS 2 , it is found that a higher density of transition metal dopants is always incorporated in sulfur-terminated domains when compared to tungsten-terminated domains. Here two representative examples demonstrate this spatial distribution control, including hexagonal iron- and vanadium-doped WS 2 monolayers. Density functional theory (DFT) calculations are further performed, indicating that the edge-dependent dopant distribution is due to a strong binding of tungsten atoms at tungsten-zigzag edges, resulting in the formation of open sites at sulfur-zigzag edges that enable preferential dopant incorporation. Based on these results, it is envisioned that edge termination in crystalline TMD monolayers can be utilized as a novel and effective knob for engineering the spatial distribution of substitutional dopants, leading to in-plane hetero-/multi-junctions that display fascinating electronic, optoelectronic, and magnetic properties.

2D transition metal dichalcogenides↗
DOE OSTI2020

Developing Scaling Relationships for Molecular Electrocatalysis through Studies of Fe-Porphyrin-Catalyzed O 2 Reduction

The oxygen reduction reaction (ORR) is a multiproton/multielectron transformation in which dioxygen (O 2 ) is reduced to water or hydrogen peroxide and serves as the cathode reaction in most fuel cells. The ORR (O 2 + 4e – + 4H + → 2H 2 O) involves up to nine substrates and thus requires navigating a complicated reaction landscape, typically with several high-energy intermediates. Many catalysts can perform this reaction, though few operate with fast rates and at low overpotentials (close to the thermodynamic potential). Attempts to optimize these parameters, both in homogeneous and heterogeneous electrocatalytic systems, have focused on modifying catalyst design and understanding kinetic/thermodynamic relationships between catalytic intermediates. One such method for analyzing and predicting catalyst reactivity and efficiency has been the development of “molecular scaling relationships”. In this work, we share our experience deriving and utilizing molecular scaling relationships for soluble, iron-porphyrin-catalyzed O 2 reduction in organic solvents. These relationships correlate turnover frequencies (TOF max ) and effective overpotentials (η eff ), properties uniquely defined for homogeneous catalysts. Following a general introduction of scaling relationships for both homogeneous and heterogeneous electrocatalysis, we describe the components of such scaling relationships: (i) the overall thermochemistry of the reaction and (ii) the rate and rate law of the catalyzed reaction. We then show how connecting these thermodynamic and kinetic parameters reveals multiple molecular scaling relationships for iron-porphyrin-catalyzed O 2 reduction. For example, the log(TOF max ) responds steeply to changes in ηeff that result from different catalyst reduction potentials (18.5 decades in TOF max /V in η eff ) but much less dramatically to changes in η eff that arise from varying the pKa of the acid buffer (5.1 decades in TOF max /V in η eff ). Thus, a single scaling relationship is not always sufficient for describing molecular electrocatalysis. This is particularly evident when the catalyst identity and reaction conditions are coupled. Using these multiple scaling relationships, we demonstrate that the metrics of turnover frequency and effective overpotential can be predictably tuned to achieve faster rates at lowered overpotentials. This Account uses a collection of related stories describing our research on soluble iron-porphyrin-catalyzed ORR to show how molecular scaling relationships can be derived and used for any electrocatalytic reaction. Such scaling relationships are powerful tools that connect the thermochemistry, mechanism, and rate law for a catalytic system. We hope that this collection shows the utility and simplicity of the molecular scaling approach for understanding catalysis, for enabling direct comparisons between catalyst systems, and for optimizing catalytic processes.

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

Kinetic description of site ensembles on catalytic surfaces

Significance The kinetic assessment of catalytic reactions is essential to the design and application of thermochemical processes useful in petrochemical synthesis, energy conversion, and environmental remediation. We demonstrate that the ubiquitously used Langmuir–Hinshelwood formalism is incomplete and fails to correctly describe even simple catalytic reactions (e.g., A + A → A 2 ). An accurate kinetic description requires explicit rate equations for multisite configurations (e.g., A*–A* pairs) to account for the nonrandom clustering/isolation of surface species. We rigorously describe clustering/isolation phenomena by derivation of coverage-dependent stoichiometric coefficients, which capture the distinct influence each elementary step exerts on each catalytic site ensemble—thereby overcoming limitations of the mean-field (i.e., Hinshelwood) assumption in describing reactions of mobile Langmuirian adsorbates on catalytic surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗
DOE OSTI2020

Anatomical Fraction segmentation in the biomass bales

The U.S. Department of Energy Bioenergy Technologies Office (BETO) is committed to the development of sustainable, nationwide, commercial biofuel production to displace petroleum-derived fuels, increase domestic energy production, and encourage the creation of a domestic bioenergy and bioproducts industry. Operating a commercial-scale bioenergy operation requires significant technological advancements for determining biomass feedstock quality during the preprocessing stage. Usage of non-food feedstocks- e.g., corn-stover, pine residue, and forest residue - at the biorefineries reduces strain on the food supply chain. But the non-food feedstock heterogeneity- physical size, shape, and chemical composition-poses a significant challenge during milling, conveyance, feeding, and biofuel conversion processes. 3D X-Ray imaging (CT) is capable of distinguishing material features by detecting density and compositional differences. But the inherent complexity introduced during harvesting and bailing makes the reconstruction and interpretation of baled biomass materials from x-ray data time consuming, laborious, and expensive. Feedstocks are low-density materials that produce small contrast differences in the x-ray images. This paper focuses on using the shape and texture properties, using 3D image processing techniques like 3D skeletonization, directional statistics to characterize and extract volumetric content of the different tissue samples in the biomass bales.

09 BIOMASS FUELS↗
DOE OSTI2025

Hyperspectral x-ray imaging mapping capabilities for nuclear forensics

Nuclear forensics relies on the integration of complementary signatures to constrain the origins and history of materials. Outcomes benefit from the timeliness and precision of the disparate methods that form typical analysis chains. Sample forms are often either minute in quantity or contain signatures like morphology or composition heterogeneity encoded on a microscale, so many analysis techniques focus on resolving signatures on ever-smaller length scales. The new hyperspectral x-ray imaging (HXI) instrument developed at Los Alamos National Laboratory seeks to improve the information available from scanning electron microscopy (SEM) x-ray spectrum analysis through superior spectral energy resolution vs. typical energy dispersive spectroscopy (EDS) systems in common use in nuclear forensics and other microanalysis fields. Based on arrays of transition-edge sensor (TES) microcalorimeter detectors, this instrument achieves a typical energy resolution of 7 eV full-width at half-maximum (FWHM) at 2 keV, opening new possibilities in trace element detection/analysis and chemical state determination through spectral shape shifts. We present here some of the first applications of the HXI instrument to actinide samples and discuss potential maturation of this nascent technology for future analysis pipelines.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2008

Mixed Mode Fracture of Plasma Sprayed Thermal Barrier Coatings: Effects of Anisotropy and Heterogeneity

The combined mode I-mode II fracture behavior of anisotropic ZrO2-8wt%Y2O3 thermal barrier coatings was determined in asymmetric flexure loading at both ambient and elevated temperatures. A fracture envelope of KI versus KII was determined for the coating material at ambient and elevated temperatures. Propagation angles of fracture as a function of KI/KII were also determined. The mixed-mode fracture behavior of the microsplat coating material was modeled using Finite Element approach to account for anisotropy and micro cracked structures, and predicted in terms of fracture envelope and propagation angle using mixed-mode fracture theories.

Zhu, Dongming↗
DOE OSTI2022

Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard

Here we perform numerical simulations of sequences of earthquake and aseismic slip on planar rate and state faults separating dissimilar material within the 2-D plane strain approximation. We resolve all stages of the earthquake cycle from aseismic slip to fast ruptures while incorporating full inertia effects during seismic event propagation. We show that bimaterial coupling results in favorable nucleation site and subsequent asymmetric rupture propagation. We demonstrate that increasing the material contrast enhances this asymmetry leading to higher slip rates and normal stress drops in the preferred rupture propagation direction. The normal stress drop, induced by the bimaterial effect, leads to strong dynamic weakening of the fault and may destabilize the creeping region on a heterogeneous rate and state fault, resulting in extended rupture propagation. Such rupture penetration into creeping patches may lead to more frequent opening of earthquake gates, causing increased seismic hazard. Furthermore, bimaterial coupling may lead to irregular seismicity pattern in terms of event length, peak slip rates, and hypocenter location, depending on the properties of the creeping patches bordering the seismogenically active part of the fault. Our results highlight robust characteristics of bimaterial interfaces that persist over long sequence of events and suggest the need for further exploration of the role of material contrast in earthquake physics and models of seismic hazard.

58 GEOSCIENCES↗
DOE OSTI2021

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Response of Sulfonated Polystyrene Melts to Nonlinear Elongation Flows

Ionizable polymers form dynamic networks with domains controlled by two distinct energy scales, ionic interactions and van der Waals forces; both evolve under elongational flows during their processing into viable materials. A molecular level insight of their nonlinear response, paramount to controlling their structure, is attained by fully atomistic molecular dynamics simulations of a model ionizable polymer, polystyrene sulfonate. As a function of increasing elongational flow rate, the systems display an initial elastic response, followed by an ionic fraction-dependent strain hardening, stress overshoot, and eventually strain-thinning. As the sulfonation fraction increases, the chain elongation becomes more heterogeneous. Finally, flow-driven ionic assembly dynamics that continuously break and reform control the response of the system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of non-uniform void distributions on the yielding of metals

High-throughput (several thousand) calculations have been carried out to investigate the yield behavior of porous materials with randomly distributed pores, porosity levels over four orders of magnitude and up to a hundred pores per simulation box. To this end, a Galerkin based fast Fourier transform (FFT) formulation was enhanced to deal with high phase contrast materials. In addition, GPU parallelization was employed in solving the governing equation for strain fluctuations using a Krylov iterative solver. Emphasis is laid on the conditions under which percolation of plastically non-deforming zones through the porous network emerge, a regime termed unhomogeneous yielding. By way of contrast, the regime where the plastic strain fluctuations (associated with the heterogeneous void-matrix aggregate) fall below the percolation threshold is defined as homogeneous yielding. Here, we find that nonuniform pore distributions only affect unhomogeneous yielding and have a universal softening effect. The extent of this distribution softening is analyzed as a function of porosity, cell size and number of realizations. Whether the uncovered universal distribution softening has direct implications on failure resistance of porous materials is discussed.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Probing the Effect of Electrode Thermodynamics on Reaction Heterogeneity in Thick Battery Electrodes

Thick electrodes present a viable strategy for enhancing energy density and reducing manufacturing costs of lithium-ion batteries. However, reaction heterogeneity during cycling compromises their rate capability and cycle life. While this nonuniformity is commonly attributed to sluggish charge transport, it is demonstrated here that the thermodynamic properties of the electrode material play an equally critical role. Through combined X-ray fluorescence microscopy and absorption near-edge structure spectroscopy, reaction distributions in LiFePO 4 (LFP) and LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC) thick electrodes with matched porosity and tortuosity are compared. LFP electrodes develop pronounced depth-oriented state-of-charge (SOC) gradients that worsen with increasing discharge rates, whereas NMC maintains much more uniform SOC distributions under such conditions. This difference originates from their distinct SOC dependence of equilibrium potentials and is quantifiable through a dimensionless “reaction uniformity” number. Intriguingly, LFP thick electrodes also exhibit lateral SOC variations that strengthen during slow discharge. In conclusion, the enhanced reaction uniformity in NMC correlates with better active material utilization and slower capacity fade than LFP, highlighting electrode thermodynamics as a key design consideration for thick electrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Peridynamic Models for Random Media Found by Coarse Graining

Using coarse graining, the upscaled mechanical properties of a solid with small scale heterogeneities are derived. The method maps internal forces at the small scale onto peridynamic bond forces in the coarse grained mesh. These upscaled bond forces are used to calibrate a peridynamic material model with position-dependent parameters. These parameters incorporate mesoscale variations in the statistics of the small scale system. The upscaled peridynamic model can have a much coarser discretization than the original small scale model, allowing larger scale simulations to be performed efficiently. The convergence properties of the method are investigated for representative random microstructures. In conclusion, a bond breakage criterion for the upscaled peridynamic material model is also demonstrated.

36 MATERIALS SCIENCE↗