Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bulk behavior”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Nonaffine Strains Control Ductility of Metallic Glasses

The origin of limited plasticity in metallic glasses is elusive, with no apparent link to their atomic structure. In this study, we propose that the response of the glassy structure to applied stress, not the original structure itself, provides a gauge to predict the degree of plasticity. We carried out high-energy x-ray diffraction on various bulk metallic glasses (BMGs) under uniaxial compression within the elastic limit and evaluated the anisotropic pair distribution function. We show that the extent of local deviation from the affine (uniform) deformation in the elastic regime is strongly correlated with the plastic behavior of BMGs beyond yield, across chemical compositions and sample history. The results suggest that the propensity for collective local atomic rearrangements under stress promotes plasticity.

3-dimensional systems↗
DOE OSTI2020

A Bayesian Approach for Statistical–Physical Bulk Parameterization of Rain Microphysics. Part I: Scheme Description

A novel framework is proposed for the bulk parameterization of rain microphysics: the Bayesian Observationally Constrained Statistical–Physical Scheme (BOSS). It is designed to facilitate direct constraint by observations using Bayesian inference. BOSS combines existing process-level microphysical knowledge with flexible process rate formulations and parameters constrained by observations within a Bayesian framework. Furthermore, using a raindrop size distribution (DSD) normalization method that relates DSD moments to one another via generalized power series, generalized multivariate power expressions are derived for the microphysical process rates as functions of a set of prognostic DSD moments. The scheme is flexible and can utilize any number and combination of prognostic moments and any number of terms in the process rate formulations. This means that both uncertainty in parameter values and structural uncertainty associated with the process rate formulations can be investigated systematically, which is not possible using traditional schemes. In this paper, BOSS is compared to two- and three-moment versions of a traditional bulk rain microphysics scheme (denoted as MORR). It is demonstrated that some process formulations in MORR are analytically equivalent to the generalized power expressions in BOSS using one or two terms, while others are not. BOSS is able to replicate the behavior of MORR in idealized one-dimensional rainshaft tests, but with a much more flexible and systematic design. Part II of this study describes the application of BOSS to derive rain microphysical process rates and posterior parameter distributions in Bayesian experiments using Markov chain Monte Carlo sampling constrained by synthetic observations.

58 GEOSCIENCES↗
DOE OSTI2025

Sound Velocities in Vanadium Reveal Complex Elastic Behavior at High Pressures

Compressional (VP) and shear (VS) wave velocities of polycrystalline vanadium were measured simultaneously up to 11.5 GPa at room temperature using ultrasonic interferometry in a multi-anvil press. Complex softening behavior in VS and resulting shear moduli are discovered, possibly revealing a precursor to the reported phase transition within 30–60 GPa. The current data enables a comprehensive assessment of the elastic and mechanical properties of vanadium at high pressures, including bulk and shear moduli, Young’s modulus, Poisson’s ratio, and Pugh’s ratio. Through fitting to the 3rd-order finite strain equations, the elastic moduli and their pressure derivatives were determined to be K S0 = 151 (2) GPa, G 0 = 46.9 (8) GPa, K$^{’}_{S0}$ = 3.47 (5), and G$^{’}_{0}$ = 0.62 (1). These experimental results allow us to compare with and benchmark the existing Steinberg–Guinan models for extrapolations to extreme pressure and temperature conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Source localization of Jupiter's Io dependent radio emissions

The peak frequencies of the Io-dependent part of the Jovian emissions are compared with the surface gyrofrequency determined from Jovian magnetic models in order to localize the source of Jovian radio emissions. The bulk of the Io-controlled emissions was found to be delayed by up to 70 deg of equatorial longitude from the predicted instantaneous position of the Io flux tube, with the L and S emissions both displaying this same unexpected behavior. It is suggested that the source of these emissions is delayed substantially with respect to Io either as an Alfven-wave delay or because of errors in the magnetic field models.

Aubier, Monique G.↗
NASA NTRS1995

Penetration of the interplanetary magnetic field B(sub y) magnetosheath plasma into the magnetosphere: Implications for the predominant magnetopause merging site

Magnetosheath plasma peertated into the magnetospere creating the particle cusp, and similarly the interplanetary magnetic field (IMF) B(sub y) component penetrates the magnetopause. We reexamine the phenomenology of such penetration to investigate implications for the magnetopause merging site. Three models are popular: (1) the 'antiparallel' model, in which merging occurs where the local magnetic shear is largest (usually high magnetic latitude); (2) a tilted merging line passing through the subsolar point but extending to very high latitudes; or (3) a tilted merging line passing through the subsolar point in which most merging occurs within a few Earth radii of the equatorial plane and local noon (subsolar merging). It is difficult to distinguish between the first two models, but the third implies some very different predictions. We show that properties of the particle cusp imply that plasma injection into the magnetosphere occurs most often at high magnetic latitudes. In particular, we note the following: (1) The altitude of the merging site inferred from midaltitude cusp ion pitch angle dispersion is typically 8-12 R(sub E). (2) The highest ion energy observable when moving poleward through the cusp drops long before the bulk of the cusp plasma is reached, implying that ions are swimming upstream against the sheath flow shortly after merging. (3) Low-energy ions are less able to enter the winter cusp than the summer cusp. (4) The local time behavior of the cusp as a function of B(sub y) and B(sub z) corroborates predictions of the high-latitude merging models. We also reconsider the penetration of the IMF B(sub y) component onto closed dayside field lines. Our approach, in which closed field lines ove to fill in flux voids created by asymmetric magnetopause flux erosion, shows that strich subsolar merging cannot account for the observations.

Newell, Patrick T.↗
DOE OSTI2022

Flowability of Crumbler Rotary Shear Size-Reduced Granular Biomass: An Experiment-Informed Modeling Study on the Angle of Repose

Biomass has potential as a carbon-neutral alternative to petroleum for chemical and energy products. However, complete replacement of fossil fuel is contingent upon efficient processes to eliminate undesirable characteristics of biomass, e.g., low bulk density, variability, and storage-induced quality problems. Mechanical size reduction via comminution is a processing operation to engineer favorable biomass flowability in handling. Crumbler rotary shear mill has been empirically demonstrated to produce more uniformly shaped particles with higher flowability than hammermilled biomass. This study combines modeling and experimentation to unveil fundamental understandings of the relation between granular particle characteristics and biomass flow behavior, which elucidate underlying mechanisms and guide selection of critical processing parameters. For this purpose, the impact of critical material attributes, including particle size (2–6 mm), particle shape (briquette, chip, clumped-sphere, cube, etc.), and surface roughness, on the angle of repose (AOR) of milled pine chips were investigated using discrete element method (DEM) simulations. Forest Concepts Crumbler rotary shear system is used to produce milled pine particles within the same size range considered in DEM simulations. AOR of different sets of these particles were measured experimentally to benchmark DEM results against experimental data. Specific energy consumption for the comminution of biomass with different particle size and moisture content are measured for technoeconomic analysis. Our results show that the smaller size (2 mm) of pine particle achieves better followability (i.e., smaller AOR) while the energy cost of comminution is significantly higher and bulk density is almost the same as the 6-mm pine particles. For the 2-mm particle size, Crumbles from veneer have better flow properties than Crumbles from chips. Contrarily, no significant difference was observed between the AOR of the two materials for the 6-mm particle size. Furthermore, from DEM simulations, mechanical interlocking between particles was found as a dominant factor in determining AOR of complex-shaped particles such as milled pine, which cannot be accurately captured by using simple particle shapes (e.g., mono-sphere) with a rolling resistance model. Conversely, clumped-sphere model alleviates this limitation without increasing computational cost significantly and can be used for accurate representation of biomass granular particles when simulating free-flow behavior.

09 BIOMASS FUELS↗
DOE OSTI2023

Chemical and Structural Evolution of AgCu Catalysts in Electrochemical CO 2 Reduction

Silver–copper (AgCu) bimetallic catalysts hold great potential for electrochemical carbon dioxide reduction reaction (CO 2 RR), which is a promising way to realize the goal of carbon neutrality. Although a wide variety of AgCu catalysts have been developed so far, it is relatively less explored how these AgCu catalysts evolve during CO 2 RR. The absence of insights into their stability makes the dynamic catalytic sites elusive and hampers the design of AgCu catalysts in a rational manner. Here, in this study, we synthesized intermixed and phase-separated AgCu nanoparticles on carbon paper electrodes and investigated their evolution behavior in CO 2 RR. Our time-sequential electron microscopy and elemental mapping studies show that Cu possesses high mobility in AgCu under CO 2 RR conditions, which can leach out from the catalysts by migrating to the bimetallic catalyst surface, detaching from the catalysts, and agglomerating as new particles. Besides, Ag and Cu manifest a trend to phase-separate into Cu-rich and Ag-rich grains, regardless of the starting catalyst structure. The composition of the Cu-rich and Ag-rich grains diverges during the reaction and eventually approaches thermodynamic values, i.e., Ag 0.88 Cu 0.12 and Ag 0.05 Cu 0.95 . The separation between Ag and Cu has been observed in the bulk and on the surface of the catalysts, highlighting the importance of AgCu phase boundaries for CO 2 RR. In addition, an operando high-energy-resolution X-ray absorption spectroscopy study confirms the metallic state of Cu in AgCu as the catalytically active sites during CO 2 RR. Taken together, this work provides a comprehensive understanding of the chemical and structural evolution behavior of AgCu catalysts in CO 2 RR.

Chen, Peng-Cheng↗
DOE OSTI2026

Machine learning-based interatomic potential development and phase transition analysis of ferroelectric hafnium dioxide

The ferroelectric phase (𝑃⁢𝑐⁢𝑎⁢2 1 , which is in orthorhombic symmetry) of hafnium dioxide (HfO 2 ) has gained much attention due to its potential applications in nanoelectronics and advanced memory devices. However, its complex phase behavior under external stimuli, such as pressure and temperature, remains a subject of intense investigation. This study focuses on developing a machine learning-based interatomic potential (MLIP) that is trained with data from density-functional theory (DFT) calculations to simulate phase transitions and mechanical properties of HfO 2 . The developed MLIP predicts lattice parameters, equations of state, bulk and shear moduli, and elastic constants that closely align with DFT predictions for several phases and at various pressures. Once validated, the MLIP is used to investigate the phase transitions of ferroelectric HfO 2 (𝑃⁢𝑐⁢𝑎⁢2 1 ) under both isobaric and constant stress conditions at elevated temperatures ranging from 200 to 2500 K. We used several complementary methods, including local symmetry identification, radial distribution function, and x-ray diffraction characterization, to identify interesting phase transitions among several competitive hafnia phases predicted from our simulations. The suggested methods uniformly reveal that under pure deviatoric condition, the system favors a transition from the orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase to a tetragonal (𝑃⁢4 2 /𝑛⁢𝑚⁢𝑐) phase, whereas a zero stress condition drives the system from the 𝑃⁢𝑐⁢𝑎⁢2 1 phase to another orthorhombic (𝑃⁢𝑏⁢𝑐⁢𝑛) phase. These findings provide crucial insights into stress and temperature-induced phase behavior of hafnia, guiding future experimental and theoretical studies for optimizing hafnia-based ferroelectric devices.

Ferroelectric HfO2↗
DOE OSTI2024

Magnetic structure and properties of the compositionally complex perovskite (Y 0.2 La 0.2 Pr 0.2 Nd 0.2 Tb 0.2 )MnO 3

Large configurational disorder in compositionally complex ceramics can lead to unique functional properties that deviate from traditional rules of alloy mixing. In recent years, compositionally complex oxides (CCOs) have shown intriguing magnetic behavior including long-range order, enhanced magnetic exchange couplings, and mixed phase magnetic structures. This work focuses on how large local spin disorder affects magnetic ordering in a CCO. Specifically, we investigated the A-site alloyed perovskite, (Y 0.2 La 0.2 Pr 0.2 Nd 0.2 Tb 0.2 )MnO 3 , or (5A)MnO 3 , using a combination of bulk magnetometry, synchrotron X-ray diffraction, and temperature-dependent neutron diffraction. The five A-site ions have an average spin and ionic radius nearly equal to that of Nd 3+ ions, which minimizes structural distortions and allows for an understanding of the local spin disorder effects through a direct comparison with NdMnO 3 . Our magnetometry data show that (5A)MnO 3 exhibits two distinct phase transitions associated with the A-site and B-site sublattices, as seen in NdMnO 3 , as well as the presence of domain pinning and exchange bias at low temperature, suggesting a mixed phase magnetic ground state, as seen in other magnetic CCOs. Neutron powder diffraction shows clear long-range antiferromagnetic ordering below 67 K and refines to a Pn'ma' magnetic structure at low temperature, in excellent agreement with the well-studied behavior of NdMnO 3 . The two most notable differences in (5A)MnO 3 magnetism apparent from our data are a slight suppression of the B-site ordering temperature, which is explained by a smaller Mn–O–Mn bond angle in (5A)MnO 3 than NdMnO 3 , and the presence of a magnetic susceptibility transition above the B-site ordering, which could indicate the formation of a cluster glass but requires further study. Finally, this work demonstrates a general method of isolated investigation of size and spin disorder in CCOs and motivates future work using local structure probes to better understand the effects of nanoscale clustering and local spin disorder in magnetic CCOs.

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

Redshift space power spectrum beyond Einstein-de Sitter kernels

Abstract We develop a framework to compute the redshift space power spectrum (PS), with kernels beyond Einstein-de Sitter (EdS), that can be applied to a wide variety of generalized cosmologies. We build upon a formalism that was recently employed for standard cosmology in Chen, Vlah & White (2020), and utilize an expansion of the density-weighted velocity moment generating function that explicitly separates the magnitude of the k-modes and their angle to the line-of-sight direction dependencies. We compute the PS for matter and biased tracers to 1-loop Perturbation Theory (PT) and show that the expansion has a correct infrared and ultraviolet behavior, free of unwanted divergences. We also add Effective Field Theory (EFT) counterterms, necessary to account for small-scale contributions to PT, and employ an IR-resummation prescription to properly model the smearing of the BAO due to large scale bulk flows within Standard-PT. To demonstrate the applicability of our formalism, we apply it on the ΛCDM and the Hu-Sawicki f(R) models, and compare our numerical results against the elephant suite of N-body simulations, finding very good agreement up to k = 0.27 Mpc -1 h at z = 0.5 for the first three non-vanishing Legendre multipoles of the PS. To our knowledge, the model presented in this work is the most accurate theoretical EFT-PT for modified gravity to date, being the only one that accounts for beyond linear local biasing in redshift-space. Hence, we argue our RSD modeling is a promising tool to construct theoretical templates in order to test deviations from ΛCDM using real data obtained from the next stage of cosmological surveys such as DESI and LSST.

Astronomy & Astrophysics↗
DOE OSTI2021

Critical Barriers to Successful Implementation of Earth-Abundant, Mn-Rich Cathodes for Vehicle Applications and Beyond: A Detailed Study of Low SOC Impedance

We report that earth-abundant, sustainable cathode chemistries based on Mn are increasingly attractive for enabling a broader portfolio of cathode oxides. The well-known class of lithium- and manganese-rich cathodes still represents very viable options. However, despite the important work on understanding the mechanisms of voltage fade, hysteresis, and oxygen activity, relatively little attention has been given to understanding the impedance characteristics of these electrodes. In particular, an anomalous rise in area specific impedance at lower states of charge, as well as overall impedance rise and surface damage due to electrolyte interactions, represent critical barriers to implementation. This work presents a comprehensive study of impedance behavior in cobalt-free, lithium- and manganese-rich electrodes. The use of a robust surface treatment allows for long-term behavior to be probed in the absence of surface damage, capacity loss, and impedance rise due to electrolyte interactions. The anomalous rise in impedance could not be correlated to surface changes, or surface phase formation, but could be directly correlated with the bulk processes of voltage fade and voltage hysteresis. The activated material can be explained as a percolating network of higher-voltage, layered-type sites having facile Li diffusion. Interspersed throughout this network are lower-voltage, disordered sites that represent a significant barrier.

25 ENERGY STORAGE↗
DOE OSTI2025

Evaluating the Impact of Managed EV Charging for Reliable Operation of Bulk Power Systems with High Non-Dispatchable Generation

The growth of electric vehicles (EVs) and variable-generation (VG) sources introduces new challenges for power-system operations. This study introduces a modeling framework and evaluates five EV charging strategies under projected 2040 grid conditions in the Evergy service territory with high non-dispatchable generation. Using realistic EV behavior and generation models, their impacts on system peak demand, ramp rate, and reserve capacity are evaluated. Results show that only the peak-avoidance strategy effectively reduces system peak demand, while decentralized strategies-particularly cost based dynamic charging-can exacerbate peaks due to synchronized user behavior. However, ramp-rate minimization strategy significantly reduce the stress on dispatchable generation achieving the lowest maximum absolute ramp rate (MARR) (56.81 MW) and lowest reserve requirement (2.39 GW). In contrast, unmanaged and TOU random strategies increase the stress on dispatchable generation sources with increased MARR and reserve requirements. These findings highlight the importance of coordinated, system-aware managed charging strategies to ensure reliable and affordable grid operation in the presence of EVs and VG sources.

14 - SOLAR ENERGY↗
DOE OSTI2023

A constitutive model for elastomers tailored by ionic bonds and entanglements

Over the past decade or two, the concept has emerged of using multiple types of weak interactions simultaneously to enhance the mechanical properties of elastomers. These weak interactions include physical entanglements, hydrogen bonds, metal-coordination bonds, dynamic covalent bonds, and ionic bonds. The combination of entanglements and ionic bonding has been minimally explored and is particularly exciting because of the broad application space for polyelectrolytes. In this work, a constitutive model framework is developed to describe the response of elastomers with both ionic bonds and entanglements. Here we formulate a micromechanical model that couples together chain stretching, ionic bond slipping, and entanglement evolution. The ionic bonds provide toughness by enabling plastic deformation in comparison to covalently crosslinked material and add strength compared to a linear polymer. Evolution of the entanglement density is taken as a key mechanism that can govern stiffness, toughness, and self-recovery in elastomers. The model is used to match bulk polyelectrolytes with different fractions of ionic components under a variety of loading histories. The variations in material parameters are then used to help understand the relative importance of different governing mechanisms in the bulk polymers. We show that the theoretical framework can explain our experimental uniaxial tensile experimental results for polyelectrolytes. This model can help to design better material with high stiffness and toughness. We expect that our model can be extended to explain the mechanical behavior of other polyelectrolytes and other soft materials with a wide range of dynamic bonds.

42 ENGINEERING↗
DOE OSTI2026

Response to the comments on creep-induced elemental redistribution at grain boundaries of 304L stainless steel – An experimental evidence for diffusional creep mechanisms

Sandström and He reanalyzed our 304L creep data (Acta Materialia 121137, 2025) at 10 MPa/750 °C and proposed that: (i) primary dislocation creep dominates based on the ϕ -model; (ii) measured creep rates deviate from published 304H data and Coble creep predictions using literature grain boundary (GB) diffusivity; and (iii) creep-induced elemental redistribution (CIER) could be explained by dislocation emission/absorption (DEAA) and GB creep (GBC), as opposed to diffusional creep. Here, in this response, we propose that: (i) the creep behavior at 15 MPa/750 °C—for which we demonstrated CIER—is only partially captured by the ϕ -model; (ii) direct comparison between 304L and 304H creep data is inappropriate due to carbides in 304H, and Coble creep predictions estimated from CALPHAD-based bulk diffusivity span experimental rates; and (iii) CIER exhibits distinct fingerprint for diffusional creep characterized by anisotropy and stress orientation dependence whereas neither DEAA nor GBC have been mechanistically substantiated.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

Tribochemical changes of alloy 800HT under sliding contact at elevated temperature in impure helium environment

This work investigates the chemical changes of the surface of Alloy 800HT during sliding contact under helium atmosphere in the presence of impurities at 750 °C, which is of interest for the nuclear industry. The alloy samples are exposed to helium (He) atmosphere at high temperature for long duration of 500 h to identify the chemical and topographical changes of the alloy surface. The oxidation rate constant is reported based on the weight change data. For the tribo-experimental samples, the surfaces inside and outside the wear track are studied thoroughly using techniques such as secondary electron microscopy, energy dispersive spectroscopy, x-ray diffraction, Raman spectroscopy and secondary ion mass spectroscopy. The formation and evolution of surface oxides are discussed in He atmosphere, compared with air. The presence of Fe-rich thin layer of oxide on top of a thicker Cr-rich oxide and a mixed spinel-based inner layer on the bulk is found. The persistent presence of thicker Fe3O4 in the wear track in high temperature air proved to be the reason of lower friction property. Pre-oxidation of the samples in air or He does not improve but worsen the tribological behavior in He atmosphere.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Octane Requirements of Lean Mixed-Mode Combustion in a Direct-Injection Spark-Ignition Engine

Here, this study investigates the octane requirements of a hybrid flame propagation and controlled autoignition mode referred to as mixed-mode combustion (MMC), which allows for strong control over combustion parameters via a spark-initiated deflagration phase. Due to the throughput limitations associated with both experiments and 3-D computational fluid dynamics calculations, a hybrid 0-D and 1-D modeling methodology was developed, supported by experimental validation data. This modeling approach relied on 1-D, two-zone engine simulations to predict bulk in-cylinder thermodynamic conditions over a range of engine speeds, compression ratios, intake pressures, trapped residual levels, fueling rates, and spark timings. Those predictions were then transferred to a 0-D chemical kinetic model, which was used to evaluate the autoignition behavior of fuels when subjected to temperature–pressure trajectories of interest. Finally, the predicted autoignition phasings were screened relative to the progress of the modeled deflagration-based combustion in order to determine if an operating condition was feasible or infeasible due to knock or stability limits. The combined modeling and experimental results reveal that MMC has an octane requirement similar to modern stoichiometric spark-ignition engines in that fuels with high research octane number (RON) and high octane sensitivity (S) enable higher loads. Experimental trends with varying RON and S were well predicted by the model for 1000 and 1400 rpm, confirming its utility in identifying the compatibility of a fuel’s autoignition behavior with an engine configuration and operating strategy. However, the model was not effective in predicting (nor designed to predict) operability limits due to cycle-to-cycle variations, which experimentally inhibited operation of some fuels at 2000 rpm. Putting the operable limits and efficiency from MMC in the context of a state-of-the-art engine, the MMC showed superior efficiencies over the range investigated, demonstrating the potential to further improve fuel economy.

33 ADVANCED PROPULSION SYSTEMS↗