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

Direct Measure of Electrode Spatial Heterogeneity: Influence of Processing Conditions on Anode Architecture and Performance

The spatial (in)homogeneity of aqueous processed silicon electrodes using standard poly(acrylic acid)-based binders and slurry preparation conditions is demonstrated in this work. X-ray nanotomography shows segregation of materials into submicron-thick layers depending on the mixing method and starting binder molecular weights. Using a dispersant, or in situ production of dispersant from the cleavage of the binder into smaller molecular weight species, increases the resulting lateral homogeneity while drastically decreasing the vertical homogeneity as a result of sedimentation and separation due to gravitational forces. This data explains some of the variability in the literature with respect to silicon electrode performance and demonstrates two potential ways to improve slurry-based electrode fabrications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Characterization of pitch carbon coating properties affecting the electrochemical behavior of silicon nanoparticle lithium-ion battery anodes

Silicon is an exciting material for next-generation lithium-ion battery anodes, due to its high theoretical capacity and availability, but its widespread implementation has been limited by extensive volume changes during cycling that causes mechanical damage and limits cycle life. One approach to mitigating these deleterious effects while still maintaining silicon's benefits is the addition of a pitch carbon coating to nanometer-sized silicon particles. Here, in this study, we present characterization of pitch-coated silicon electrodes without annealing and annealed at 700 °C and 1000 °C to determine the optimum temperature treatment, as well as pitch-only and silicon-only electrodes with and without annealing to elucidate the impacts of each component on early electrochemical behavior. Using Fourier transform infrared and Raman spectroscopies, atomic force and scanning spreading resistance microscopy, and electrochemical analysis, we find 700 °C to be the optimal annealing temperature, as amorphous carbon is created, which improves conductivity prior to cycling and facilitates ion storage during cycling that increases capacity. We also present in situ Raman spectroscopy data demonstrating the heterogeneous aging that occurs across the electrode surface during cycling.

25 ENERGY STORAGE↗
DOE OSTI2020

The influence of packing structure and interparticle forces on ultrasound transmission in granular media

Ultrasound propagation through externally stressed, disordered granular materials was experimentally and numerically investigated. Experiments employed piezoelectric transducers to excite and detect longitudinal ultrasound waves of various frequencies traveling through randomly packed sapphire spheres subjected to uniaxial compression. The experiments featured in situ X-ray tomography and diffraction measurements of contact fabric, particle kinematics, average per-particle stress tensors, and interparticle forces. The experimentally measured packing configuration and inferred interparticle forces at different sample stresses were used to construct spring networks characterized by Hessian and damping matrices. The ultrasound responses of these network were simulated to investigate the origins of wave velocity, acoustic paths, dispersion, and attenuation. Results revealed that both packing structure and interparticle force heterogeneity played an important role in controlling wave velocity and dispersion, while packing structure alone quantitatively explained most of the observed wave attenuation. This research provides insight into time- and frequency-domain features of wave propagation in randomly packed granular materials, shedding light on the fundamental mechanisms controlling wave velocities, dispersion, and attenuation in such systems.

36 MATERIALS SCIENCE↗
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↗
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↗