Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic temperature”

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,405 records · Page 78

Impacts of Latent Energy and Snow Fall Speed on a Wintertime Midlatitude Cyclone

This study aimed to understand the impacts of latent energy as well as snow fall speeds on precipitation properties and synoptic-scale storm characteristics of a wintertime midlatitude cyclone. Simulations of a potent winter storm that impacted the Rocky Mountains and northern Great Plains of the United States in February 2017 were performed using the Weather Research and Forecasting model with latent heating or cooling from individual microphysical processes systematically turned off and fall speeds adjusted by ±50%. Results indicated substantial impacts on the microphysical characteristics of the simulated storm to fall speed, cooling from sublimation, and warming from deposition. The impacts of cooling and warming were manifested as differences in accumulated snowfall. Increased (decreased) fall speeds led to smaller (larger) ice crystals and total mass, resulting in offsetting effects on the precipitation flux, and thus minimal impacts on snowfall and large-scale characteristics of the storm. Warming and cooling associated with deposition and sublimation, respectively, impacted the synoptic-scale dynamics, whereby removing warming from deposition resulted in an increased meridional temperature gradient near the jet stream, thus increasing the jet strength and causing it to be more westerly with less curvature aloft. This in turn limited upper-level divergence, creating a weaker surface low and shifting the precipitation shield southward. The opposite occurred with the removal of latent cooling due to sublimation. This study highlights the potential importance of latent energy associated with ice sublimation and deposition and fall speeds in the larger-scale characteristics of winter storms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A high-throughput framework for lattice dynamics

We develop an automated high-throughput workflow for calculating lattice dynamical properties from first principles including those dictated by anharmonicity. The pipeline automatically computes interatomic force constants (IFCs) up to 4th order from perturbed training supercells, and uses the IFCs to calculate lattice thermal conductivity, coefficient of thermal expansion, and vibrational free energy and entropy. It performs phonon renormalization for dynamically unstable compounds to obtain real effective phonon spectra at finite temperatures and calculates the associated free energy corrections. The methods and parameters are chosen to balance computational efficiency and result accuracy, assessed through convergence testing and comparisons with experimental measurements. Deployment of this workflow at a large scale would facilitate materials discovery efforts toward functionalities including thermoelectrics, contact materials, ferroelectrics, aerospace components, as well as general phase diagram construction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Uncovering the hydride ion diffusion pathway in barium hydride via neutron spectroscopy

Solid state materials possessing the ability for fast ionic diffusion of hydrogen have immense appeal for a wide range of energy-related applications. Ionic hydrogen transport research is dominated by proton conductors, but recently a few examples of hydride ion conductors have been observed as well. Barium hydride, BaH 2 , undergoes a structural phase transition around 775 K that leads to an order of magnitude increase in the ionic conductivity. This material provides a prototypical system to understand hydride ion diffusion and how the altered structure produced by the phase transition can have an enormous impact on the diffusion. We employ quasielastic and inelastic neutron scattering to probe the atomic scale diffusion mechanism and vibrational dynamics of hydride ions in both the low- and high-temperature phases. Jump lengths, residence times, diffusion coefficients, and activation energies are extracted and compared to the crystal structure to uncover the diffusion pathways. We find that the hydrogen jump distances, residence times, and energy barriers become reduced following the phase transition, allowing for the efficient conduction of hydride ions through a series of hydrogen jumps of length L = 3.1 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast optical melting of trimer superstructure in layered 1T′-TaTe2

Abstract Quasi-two-dimensional transition-metal dichalcogenides are a key platform for exploring emergent nanoscale phenomena arising from complex interactions. Access to the underlying degrees-of-freedom on their natural time scales motivates the use of advanced ultrafast probes sensitive to self-organised atomic-scale patterns. Here, we report the ultrafast investigation of TaTe 2 , which exhibits unique charge and lattice trimer order characterised by a transition upon cooling from stripe-like chains into a (3 × 3) superstructure of trimer clusters. Utilising MeV-scale ultrafast electron diffraction, we capture the photo-induced TaTe 2 structural dynamics – exposing a rapid ≈ 1.4 ps melting of its low-temperature ordered state followed by recovery via thermalisation into a hot cluster superstructure. Density-functional calculations indicate that the initial quench is triggered by intra-trimer Ta charge transfer which destabilises the clusters, unlike melting of charge density waves in other TaX 2 compounds. Our work paves the way for further exploration and ultimately rapid optical and electronic manipulation of trimer superstructures.

36 MATERIALS SCIENCE↗

Charge carrier dynamics in two-dimensional hybrid perovskites: Dion–Jacobson vs. Ruddlesden–Popper phases

Two-dimensional (2D) halide perovskites are promising materials for environmentally stable next-generation optoelectronic device applications. Besides the widely investigated Ruddlesden–Popper (RP) phases, Dion–Jacobson (DJ) phases are attracting considerable attention due to their rapid emergence as efficient solar cell materials. However, there is very little atomistic understanding of the charge carrier dynamics under ambient conditions for these DJ-phases, limiting the possibilities to tune their optoelectronic performances through compositional engineering routes. In this study, by combining nonadiabatic molecular dynamics with time-domain density functional theory methods at room temperature, we compare the dominant non-radiative carrier recombination and dephasing processes in RP and DJ monolayered lead halide perovskites. Our systematic study demonstrates that performance-limiting nonradiative carrier recombination processes greatly depends on the electron–phonon interactions induced by structural fluctuations and instantaneous charge localization in these materials. The stiffer interlayer packing due to the presence of single spacer dications, which separates the lead iodide slabs, reduces the thermal fluctuations in the DJ-phase to a greater extent than that in the RP-phase 2D-perovskites. Specific electronic coupling between the closely spaced lead iodide layers enhances the delocalization of band-edge charge densities in DJ-phase systems. Compared to the RP-phase, reduced inelastic electron–phonon scattering in DJ-phase perovskites significantly limits intrinsic non-radiative recombination processes. The consequent enhancement in the photogenerated charge carrier lifetime makes DJ-phase perovskites potentially suitable for various optoelectronic devices. The computational insights gained from this study allow us to outline a set of robust design principles for DJ-phase perovskites to strategically tune their optoelectronic properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast THz emission spectroscopy of spin currents in the metamagnet FeRh

Heterostructures of ferromagnetic (FM) and noble metal (NM) thin films have recently attracted considerable interest as viable platforms for the ultrafast generation, control, and transduction of light-induced spin currents. In such systems, an ultrafast laser can generate a transient spin current in the FM layer, which is then converted to a charge current at the FM/NM interface due to strong spin–orbit coupling in the NM layer. Whether such conversion can happen in a single material and how the resulting spin current can be quantified are open questions under active study. Here, we report ultrafast THz emission from spin–charge conversion in a bare FeRh thin film without any NM layer. Our results highlight that the magnetic material by itself can enable spin–charge conversion in the same order as that in a FM/NM heterostructure. We further propose a simple model to estimate the light-induced spin current in FeRh across its metamagnetic phase transition temperature. Our findings have implications for the study of the ultrafast dynamics of magnetic order in quantum materials using THz emission spectroscopy.

36 MATERIALS SCIENCE↗

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗

Mixed stochastic-deterministic density functional theoretic decomposition of Kubo–Greenwood conductivities in the projector augmented wave formalism

Pairing the accuracy of the Kohn–Sham density-functional framework with the efficiency of a stochastic algorithmic approach, mixed stochastic-deterministic density functional theory (mDFT) achieves a favorable computational scaling with system sizes and electronic temperatures. We employ the recently developed mDFT formalism to investigate the dynamic charge-transport properties of systems in the warm dense matter regime. The optical conductivity spectra are computed for single- and multi- component mixtures of carbon, hydrogen, and beryllium using two complementary approaches: Kubo–Greenwood in the mDFT picture and real-time time-dependent mDFT. We further devise a decomposition of the Onsager coefficients leading up to the Kubo–Greenwood spectra to exhibit contributions from the deterministic, stochastic, and mixed electronic state transitions at different incident photon energies.

36 MATERIALS SCIENCE↗

Particle acceleration during classical phase transitions on a spherical lattice

Abstract When compressed, certain lattices undergo phase transitions that may allow nuclei to gain significant kinetic energy. To explore the dynamics of this phenomenon, we develop a methodology to study Coulomb coupled N -body systems constrained to a sphere, as in the Thomson problem. We initialize N total Boron nuclei as point particles on the surface of the sphere, allowing them to equilibrate via Coulomb scattering with a viscous damping term. To simulate a phase transition, we remove N rm particles, forcing the system to rearrange into a new equilibrium. With this model, we consider the Thomson problem as a dynamical system, providing a framework to explore how non-zero temperature affects structural imperfections in Thomson minima. We develop a scaling relation for the average peak kinetic energy attained by a single particle as a function of N and N rm . For certain values of N , we find an order of magnitude energy gain when increasing N rm from 1 to 6. The model may help to design a lattice that maximizes the energy output.

Bachmann, Aidan M. (ORCID:0000000294706404)↗

Normal mode description of phases of matter: Application to heat capacity

Understanding thermodynamics in liquids at the atomic level is challenging because of strong atomic interactions and lack of spatial symmetry. Recent prior theoretical works have focused on describing heat capacity of liquids in terms of phonon-like excitations but often rely on fitting factors and assumptions. In this work, we propose characterizing various phases in terms of instantaneous normal modes (INMs) of structural snapshots from molecular dynamics simulations of single-element systems over wide ranges of temperature and pressure. We use the INMs to build a mode-level microscopic description of heat capacity and demonstrate that heat capacity of liquids can be described by a combination of both solidlike and gaslike degrees of freedom, leading to a more unified framework to fundamentally describe heat capacity of all three phases of matter: solid, liquid, and gas. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Distributionally Robust Bilevel Optimization Model for Distribution Network With Demand Response Under Uncertain Renewables Using Wasserstein Metrics

Here, we consider a distribution network integrating demand response (DR) participants in the presence of uncertain renewable suppliers and outdoor temperatures. A bilevel optimization model is proposed to capture the intricate dynamics between price-incentivized DR participants and distribution system operations, including energy procurement and active/reactive power flows. The model is formulated as a distributional robust bilevel optimization using Wasserstein metrics. We show favorable data-driven properties including out-of-sample guarantee and asymptotic consistency. Furthermore, we present a tractable mixed-integer linear programming reformulation and characterize the worst-case distribution. Computational experiments are conducted on a modified 33-bus system. Our findings underscore the efficacy of the pricing strategies derived from the proposed bilevel optimization model. These strategies not only effectively manage DR participants' behavior but also bring equity considerations among households with various characteristics to light. The results contribute to a deeper understanding of the interplay between distribution system operators and DR participants.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Flux Detection Probe to Quantify Dynamic Groundwater-Surface Water Exchange in the Hyporheic Zone

A new probe was designed to quantify groundwater-surface water exchange in the hyporheic zone under dynamic stage condition. Current methods focus on either vertical pore water velocity or Darcy flux measurements. Both parameters must be understood to evaluate residence time and mass flux of constituents. Furthermore, most instruments are not well suited for monitoring instantaneous velocity or flux under dynamic exchange conditions. For this reason, the flux detection probe (FDP) was designed that employs electrogeophysical measurements to estimate in situ sediment porosity, which can be used to convert pore water velocity to Darcy flux. Dynamic pore water velocity is obtained by monitoring fluid conductivity and temperature along the FDP probe. Pressure sensors deployed at the top and bottom of the probe provide the additional information necessary to estimate vertical permeability. This study focuses on the use of a geophysical method to estimate pore water velocity, porosity, and permeability within a controlled soil column where simulated river water displaces simulated groundwater. The difference between probe derived and theoretical pore water velocity using natural tracers such as electrical conductivity and temperature was -4.9 and 3.9% for downward flow and 1.1 and 12.8% for upward flow, respectively. The difference in porosity calculated from mass and volume packed in the soil column and probe measure porosity ranged between -3.2% and 1.5%. Also, the calculated hydraulic conductivity differed from probe derived values by -8.9%.

54 ENVIRONMENTAL SCIENCES↗

Direct Numerical Simulation of the Flow through a Randomly Packed Pebble Bed

The proposition for molten salt and high-temperature gas-cooled reactors has increased the focus on the dynamics and physics in randomly packed pebble beds. Research is being conducted on the validity of these designs as a possible contestant for the fourth-generation nuclear power systems. A detailed understanding of the coolant flow behavior is required in order to ensure proper cooling of the reactor core during normal and accident conditions. In order to increase the understanding of the flow through these complex geometries and enhance the accuracy of lower-fidelity modeling, high-fidelity approaches such as direct numerical simulation (DNS) can be utilized. Nek5000, a spectral-element computational fluid dynamics (CFD) code, was used to develop DNS fluid flow data. The flow domain consisted of 147 pebbles enclosed by a bounding wall. In the work presented, the Reynolds numbers ranged from 430 to 1050 based on the pebble diameter and inlet velocity. Characteristics of the flow domain such as volume averaged porosity, axial porosity, and radial porosity were studied and compared with correlations available in the literature. Friction factors from the DNS results for all Reynolds numbers were compared with correlations in the literature. The first- and second-order statistics show good agreement with the available experimental data. Turbulence length scales were analyzed in the flow. Reynolds stress anisotropy was characterized by utilizing invariant analysis. Overall, the results of the analysis in this study provide deeper understanding of the flow behavior and the effect of the wall in packed beds.

Yildiz, Mustafa Alper↗

Massively distributed fiber strain sensing using Brillouin lasing

Brillouin based distributed fiber sensors present a unique set of characteristics amongst fiber sensing architectures. They are able to measure absolute strain and temperature over long distances, with high spatial resolution, and very large dynamic range in off-the-shelf fiber. However, Brillouin sensors traditionally provide only modest sensitivity due to the weak dependence of the Brillouin frequency on strain and the high signal to noise ratio required to identify the resonance’s peak frequency to within a small fraction of its linewidth. Recently, we introduced a technique which substantially improves the precision of Brillouin fiber sensors by exciting a series of lasing modes in a fiber loop cavity that experience Brillouin amplification at discrete locations in the fiber. The narrow-linewidth and high intensity of the lasing modes enabled ultra-low noise Brillouin sensors with large dynamic range. However, our initial demonstration was only modestly distributed: measuring strain at 40, non-contiguous positions along a 400 m fiber. In this work, we greatly extend this methodology to enable fully distributed sensing at 1000 contiguous locations along 3.5 km of fiber—an order of magnitude increase in sensor count and range. This highly-multiplexed Brillouin fiber laser sensor provides a strain noise as low as 34 nɛ/√Hz and we analyze the limiting factors in this approach.

47 OTHER INSTRUMENTATION↗

LLNL FESP Theory Highlights: November 2024

This talk presented a refined understanding of plasma-facing thermionic emitter cathodes by examining the limitations of the conventional space charge-limited model, which overlooks ion trapping dynamics within the virtual cathode. While conventional theory distinguishes between temperature-limited and space charge-limited regimes based on electron current constraints, recent findings reveal that trapped ions can form a quasi-neutral "second plasma," expanding upstream and interacting with the primary plasma. These effects are amplified in non-planar geometries, such as cylindrical or spherical sheaths, where trapped ions create a positive feedback loop with the upstream plasma. To address these phenomena, a new "Aid-and-Compete" model was developed, offering insights into current enhancement and sheath mode transitions relevant for plasma devices, including emissive cathodes in electric propulsion, fusion applications, thin film deposition systems, and plasma diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Filament Formation and Melt Spinning of Coal-Based Mesophase Pitch for Carbon Fiber Production

High-performance carbon fibers excel as high specific strength and modulus materials, and are utilized in composites applications ranging across aerospace, automotive, energy, infrastructure, and sports equipment sectors. Compared to PAN-based carbon fiber, mesophase pitch-based carbon fiber has lower tensile strength but provides a higher modulus, higher thermal conductivity, and a lower cost potential. However, challenges in stable, continuous melt spinning processing are a serious limiting factor. Unlike typical melt spinning of long, linear chain polymers, mesophase pitch is a comprised of relatively shorter polycyclic aromatic hydrocarbons which form a liquid crystal. In its ‘green’ or ‘as-spun’ state, the fibers are very fragile. Moreover, its temperature of processing is quite high, often approaching 400ºC. Complex flow dynamics combined with short length and time scale heat transfer render its melt spinning a formidable processing challenge. Improved understanding of the root causes for nascent filament breakage and spinning instabilities are needed. This work aims to determine fundamental phenomena that govern mesophase pitch filament formation by comparison of required draw force and uninterrupted spinning minutes with applied draw down ratio. The effect of perturbations on spinning stability are analyzed using capillary rheometry, microscopy of elongating filaments, and filtration.

01 COAL, LIGNITE, AND PEAT↗

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

Measurements and Prediction of Ash Deposition in a Cyclone-Fired Boiler Operating under Variable Load Conditions

Measurements of ash deposition rates were made between the secondary superheater and reheater sections of a 450 MW cyclone-fired lignite boiler as the operational load varied from 33 to 100%. Significant reductions in deposition rates with a decrease in operational load were observed. To uncover the causative mechanisms behind these observations, operational data from the power plant were used to carry out computational fluid dynamic (CFD) simulations of the boiler. After ascertaining that the gas temperatures and velocities at various sections within the boiler were being represented adequately, decoupled simulations of the ash deposition process on the deposit probe were carried out using a finely resolved boundary layer mesh. Fly ash particle size distribution (PSD) and its concentration for the decoupled calculations were determined from stand-alone cyclone barrel simulations. The ash partitioning (mass %) between the fly ash and slag was found to be ~50:50, which was in line with previous field observations, and it did not vary significantly across different cyclone loads. The predicted PSD of the deposit ash was concentrated in the size range 10–30 microns, which was in agreement with cross-sectional images of the deposit obtained from the measurements. At lower loads, sharp variations in the deposition rates were predicted in the gas temperature range 950–1150 K. The particle kinetic energy—particle viscosity-based capture methodology utilized in this study in conjunction with appropriate ash compositions, ash viscosity models and gas temperature estimates can help estimate slagging propensities at different loads reasonably well in these systems.

Krishnamoorthy, Gautham (ORCID:0000000255205092)↗