Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “path average”

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 37 records · Page 2

Thickness-Dependent Cross-Plane Thermal Conductivity Measurements of Exfoliated Hexagonal Boron Nitride

Submicrometer-thick layers of hexagonal boron nitride (hBN) exhibit high in-plane thermal conductivity and useful optical properties, and serve as dielectric encapsulation layers with low electrostatic inhomogeneity for graphene devices. Despite the promising applications of hBN as a heat spreader, the thickness dependence of its cross-plane thermal conductivity is not known, and the cross-plane phonon mean free paths (MFPs) have not been measured. We measure the cross-plane thermal conductivity of hBN flakes exfoliated from bulk crystals. We find that submicrometer thick flakes exhibit thermal conductivities up to 8.1 ± 0.5 W m –1 K –1 at 295 K, which exceeds previously reported bulk values by more than 60%. Surprisingly, the average phonon mean free path is found to be several hundred nanometers at room temperature, a factor of 5 larger than previous predictions. When planar twist interfaces are introduced into the crystal by mechanically stacking multiple thin flakes, the cross-plane thermal conductivity of the stack is found to be a factor of 7 below that of individual flakes with similar total thickness, thus providing strong evidence that phonon scattering at twist boundaries limits the maximum phonon MFPs. These results have important implications for hBN integration in nanoelectronics and improve our understanding of thermal transport in twodimensional materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bounded motion design in the Earth zonal problem using differential algebra based normal form methods

Establishing long-term relative bounded motion between orbits in perturbed dynamics is a key challenge in astrodynamics to enable cluster flight with minimum propellant expenditure. Herein, we present an approach that allows for the design of long-term relative bounded motion considering a zonal gravitational model. Entire sets of orbits are obtained via high-order Taylor expansions of Poincarè return maps about reference fixed points. The high-order normal form algorithm is used to determine a change in expansion variables of the map into normal form space, in which the phase space behavior is circular and can be easily parameterized by action–angle coordinates. The action–angle representation of the normal form coordinates is then used to parameterize the original Poincarè return map and average it over a full phase space revolution by a path integral along the angle parameterization. As a result, the averaged nodal period and drift in the ascending node are obtained, for which the bounded motion conditions are straightforwardly imposed. Sets of highly accurate bounded orbits are obtained, extending over several thousand kilometers, and valid for decades.

79 ASTRONOMY AND ASTROPHYSICS↗

Specifics of Thermal Transport in Graphene Composites: Effect of Lateral Dimensions of Graphene Fillers

In this work, we report on the investigation of thermal transport in noncured silicone composites with graphene fillers of different lateral dimensions. Graphene fillers are comprised of few-layer graphene flakes with lateral sizes in the range from 400 to 1200 nm and the number of atomic planes from 1 to ~100. The distribution of the lateral dimensions and thicknesses of graphene fillers has been determined via atomic force microscopy statistics. It was found that in the examined range of the lateral dimensions, the thermal conductivity of the composites increases with increasing size of the graphene fillers. The observed difference in thermal properties can be related to the average gray phonon mean free path in graphene, which has been estimated to be around ~800 nm at room temperature. The thermal contact resistance of composites with graphene fillers of 1200 nm lateral dimensions was also smaller than that of composites with graphene fillers of 400 nm lateral dimensions. The effects of the filler loading fraction and the filler size on the thermal conductivity of the composites were rationalized within the Kanari model. The obtained results are important for the optimization of graphene fillers for applications in thermal interface materials for heat removal from high-power-density electronics.

36 MATERIALS SCIENCE↗

Demonstration of a compact, multi-joule, diode-pumped Tm:YLF laser

We report the demonstration of a diode-pumped Tm:YLF laser operating at 1.88 µm that produces pulse energies up to 3.88 J in 20 ns. The compact system consists of a Q-switched cavity-dumped oscillator generating 18 mJ pulses, which are then amplified in a 4-pass power amplifier. Energies up to 38.1 J were obtained with long-pulse amplifier operation. These results illustrate the high energy storage and extraction capabilities of diode-pumped Tm:YLF, opening the path to high peak and average power mid-infrared solid-state lasers

47 OTHER INSTRUMENTATION↗

Induced microearthquakes predict permeability creation in the brittle crust

Predicting the evolution of permeability accurately during stimulation at the reservoir scale and at the resolution of individual fractures is essential to characterize the fluid transport and the reactive/heat-transfer characteristics of reservoirs where stress exerts significant control. Here, we develop a hybrid machine learning (ML) model to visualize in situ permeability evolution for an intermediate-scale (~10 m) hydraulic stimulation experiment. This model includes an ML model that was trained using the well history of flow rate and wellhead pressure and MEQ data from the first three stimulation episodes to predict average permeability from the statistical features of the MEQs alone for later episodes. Moreover, a physics-inspired model is integrated to estimate in situ fracture permeability spatially. This method relates fracture permeability to fracture dilation and scales dilation to the equivalent MEQ magnitude, according to laboratory observations. The seismic data are then applied to define incremental changes in permeability in both space and time. Our results confirm the excellent agreement between the ground truth and model-predicted permeability evolution. The resulting permeability map defines and quantifies flow paths in the reservoir with the averaged permeability comparing favorably with the ground truth of permeability.

induced seismicity↗

Phonon gas model for thermal conductivity of dense, strongly interacting liquids

Developing predictive thermal property models for liquids based on microscopic principles has been elusive. The difficulty is that liquids have gas-like and solid-like attributes that are at odds when considering the frameworks of microscopic models: Models for gases are simple due to randomness and low density, whereas models for crystalline solids rely on symmetry and long-range order for easier calculation. The short-range order in liquids does, however, provide structure to neighboring molecules similar to amorphous solids, and there have been recent advances indicating that collective vibrational modes store heat in liquids. Models combining Debye approximations from solid-state physics and Frenkel’s theory of liquids can accurately predict the heat capacity of liquids. Phonon-like dispersions in liquids have also been widely observed in neutron scattering experiments. These developments motivate us to propose a model where high-frequency vibrational modes, which travel at the speed of sound and have a mean free path on the order of the average intermolecular distance, conduct heat in liquids. We use this liquid phonon gas model to calculate the thermal conductivity of liquids with varying intermolecular interaction energies from strongest to weakest—Coulomb, hydrogen-bonding, Keesom, and London dispersion energy. Generally, the model is more accurate as the intermolecular interaction energy and density of liquids increase. The calculated thermal conductivity of Coulombic-bound molten sodium nitrate and hydrogen-bonded water is within 1.46% and 2.98% of the experimentally measured values, respectively, across their entire temperature ranges. Further modal analysis of the velocity and the mean free path of collective vibrations could establish the liquid phonon gas model as an accurate model for weakly interacting liquids as well.

36 MATERIALS SCIENCE↗

Data-Driven Transition Path Analysis Yields a Statistical Understanding of Sudden Stratospheric Warming Events in an Idealized Model

Abstract Atmospheric regime transitions are highly impactful as drivers of extreme weather events, but pose two formidable modeling challenges: predicting the next event (weather forecasting) and characterizing the statistics of events of a given severity (the risk climatology). Each event has a different duration and spatial structure, making it hard to define an objective “average event.” We argue here that transition path theory (TPT), a stochastic process framework, is an appropriate tool for the task. We demonstrate TPT’s capacities on a wave–mean flow model of sudden stratospheric warmings (SSWs) developed by Holton and Mass, which is idealized enough for transparent TPT analysis but complex enough to demonstrate computational scalability. Whereas a recent article (Finkel et al. 2021) studied near-term SSW predictability, the present article uses TPT to link predictability to long-term SSW frequency. This requires not only forecasting forward in time from an initial condition, but also backward in time to assess the probability of the initial conditions themselves. TPT enables one to condition the dynamics on the regime transition occurring, and thus visualize its physical drivers with a vector field called the reactive current . The reactive current shows that before an SSW, dissipation and stochastic forcing drive a slow decay of vortex strength at lower altitudes. The response of upper-level winds is late and sudden, occurring only after the transition is almost complete from a probabilistic point of view. This case study demonstrates that TPT quantities, visualized in a space of physically meaningful variables, can help one understand the dynamics of regime transitions.

Meteorology & Atmospheric Sciences↗

Thermal Properties of Sodium Borosilicate Glasses as a Function of Sulfur Content

SO3 additions, up to 3.0 wt%, were systematically investigated for effects on the physical properties of sodium borosilicate glass melted in air, with a sulfur-free composition of 50SiO2-10Al2O3-12B2O3-21Na2O-7CaO (wt%). Solubility measurements, using electron microscopy chemical analysis, determined the maximum loading to be ~1.5 wt% SO3. It was found that sulfur (here as sulfate) in this glass increased the glass transition temperature, thermal diffusivity, heat capacity, and thermal conductivity, and decreased the mass density. Structural analysis, performed with Raman spectroscopy, indicated that the borosilicate network polymerized with sulfur additions, presumably due to Na2O being required to charge compensate the ionic SO42- additions, thus becoming unavailable to form non-bridging oxygen in the silicate network. It is postulated that this increased crosslinking of the borosilicate backbone lead to a structure with higher dimensionality and average bond energy. This increased the mean free paths and vibration frequency of the phonons, which resulted in the observed increase in thermal properties.

Thermal properties, Borosilicate glass, Sulfur sol↗

Operator dynamics in Floquet many-body systems

We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading. Published by the American Physical Society 2025

Yoshimura, Takato (ORCID:0000000309159846)↗

Simulation Inputs for the METS-R Simulator

This dataset served as the input for the METS-R simulator. Data include the historical and predicted demand, cache of transit scheduling results, cache of candidate paths for routing, and link-level average speed and corresponding standard deviations.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Exploring the Effects of Emplacement Conditions on Explosion P/S Ratios across Local to Regional Distances

High-frequency (~> 2 Hz) seismic P/S amplitude ratios are well-established as a discriminant to distinguish between natural earthquakes and underground explosions at regional distances (~200–1500 km). As research shifts toward identifying lower-yield events, work has begun to investigate the potential of this discriminant for use at local distances (<200 km), in which initial results raise questions about its effectiveness. Here, we utilize data from several chemical explosion experiment series at the Nevada National Security Site in southern Nevada in the United States to study explosion Pg/Lg ratios across the range of local to regional distances. Additionall, the experiments are conducted over differing emplacement conditions, with contrasting geologies and a variety of yields and depths of burial, including surface explosions. We first establish the similarities of Pg/Lg ratios from chemical explosions to those from historic nuclear tests and conclude that, as previous data have suggested, chemical explosion ratios are good proxies for nuclear tests. We then examine Pg/Lg ratios from the new experiment series as functions of distance, yield, depth of burial, and scaled depth of burial (SDOB). At far-local and regional distances, we observe consistently higher ratios from hard-rock explosions compared to ones in a weaker dry alluvium medium, consistent with prior regional distance results. No other trends with yield, depth of burial, or SDOB are strongly evident. Scatter in the observed ratios is very high, particularly at the shortest event-to-station distances, suggesting that small-scale path effects play a significant role. On average, the local distance explosion Pg/Lg ratios show remarkable consistency across all the variations in emplacement. Explosion source models will need to reproduce these results.

58 GEOSCIENCES↗

Predicting electrical conductivity in Cu/Nb composites: A combined model-experiment study

The generation of high magnetic fields requires materials with high electric conductivity and good strength properties. Cu/Nb composites are considered to be good candidates for this purpose. In this work, we aim to predict, from theory, the dependence of electric conductivity on the microstructure, most notably on the layer thickness and grain sizes. We also conducted experiments to calibrate and validate our simulations. Bimetal interfaces and grain boundaries are confirmed to have the largest impact on conductivity in this composite material. In this approach, a distribution of the layer thickness is accounted for in order to better model the experimentally observed microstructure. Because layer thicknesses below the mean free path of Cu significantly degrade the conductivity, an average layer thickness larger than expected may be needed to meet conductivity requirements in order to minimize these smaller layers in the distribution. We also investigate the effect of variations in volume fraction of Nb and temperature on the material’s conductivity.

36 MATERIALS SCIENCE↗

From anti-Arrhenius to Arrhenius behavior in a dislocation-obstacle bypass: Atomistic simulations and theoretical investigation

Dislocations are the primary carriers of plasticity in metallic materials. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate of an edge dislocation overcoming the long-range elastic interaction with a prismatic loop in tungsten. Close to the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius-like behavior. Here, we have computed the activation entropy analytically along the minimum energy path following Schoeck’s method, which captures the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI), another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck’s formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck’s and PAFI compared to molecular dynamics.

36 MATERIALS SCIENCE↗

Measurement-induced purification in large- N hybrid Brownian circuits

Competition between unitary dynamics that scrambles quantum information nonlocally and local measurements that probe and collapse the quantum state can result in a measurement-induced entanglement phase transition. Here we study this phenomenon in an analytically tractable all-to-all Brownian hybrid circuit model composed of qubits. The system is initially entangled with an equal sized reference, and the subsequent hybrid system dynamics either partially preserves or totally destroys this entanglement depending on the measurement rate. Our approach can access a variety of entropic observables which are distinguished by the averaging procedure, and for concreteness we focus on a particular purity quantity for which the averaging is particularly simple. We represent the purity as a path integral coupling four replicas with twisted boundary conditions. Saddle-point analysis reveals a second-order phase transition corresponding to replica permutation symmetry breaking below a critical measurement rate. The transition is mean-field-like and we characterize the critical properties near the transition in terms of a simple Ising field theory in 0 + 1 dimensions. In addition to studying the purity of the entire system, we study subsystem purities and relate these results to manifestations of quantum error correction in the model. Here, we also comment on the experimental feasibility for simulating this averaged purity, and corroborate our results with exact diagonalization for modest system sizes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermal transport in Yb-based 1-2-20 materials

Abstract Lattice thermal-transport properties of heavy-fermion YbT 2 Zn 20 -based (T = Co, Rh, and Ir) 1-2-20 compounds are calculated with the single-mode relaxation-time approximation and the full solution of the linearized phonon Boltzmann transport equation from first-principles anharmonic phonon calculations. We predict low lattice thermal conductivity κ L with a maximum value of ∼ 24.16 · K at 85.14 K, 137.29 · K at 9.67 K, and 23.55 · K at 67.50 K and ∼ 12.02, 40.04, and 10.30 · K at room temperature for YbCo 2 Zn 20 , YbRh 2 Zn 20 , and YbIr 2 Zn 20 , respectively. Based on the analysis of the cumulative κ L as a function of the phonon mean free path and the frequency-dependent joint density of states at various temperatures, we attribute the low κ L to a rattling mode with an average characteristic rattling frequency that enhanced phonon scattering processes, which reduced the phonon mean free paths, suppressed the phonon lifetime, and enhanced the probability of three-phonon scattering events. The predicted low κ L , especially in YbCo 2 Zn 20 and YbIr 2 Zn 20 , makes them promising candidate materials for thermoelectric applications and thermal management.

Physics↗

Grain boundary strain as a determinant of localized sink efficiency

The opportunity to achieve radiation tolerance in crystalline materials hinges on understanding the structure and response of grain boundary sinks to irradiation. A common descriptor of grain boundary efficiency as a defect sink is the denuded zone, which is a defect free zone adjacent to the grain boundary dictated by its ability to absorb radiation induced defects. This descriptor is often used at the mesoscale, which requires an averaging of absorption events. In this paper, we resolve the defect sink efficiency as a function of interfacial strain with respect to grain boundary character, and correlate high levels of grain boundary strain to an enhanced absorption efficiency. Here, we also introduce a key relationship between localized strain in proximity with the grain boundary sink and the variation absorption efficiency associated with these regions, revealing the pitfalls of averaging absorption events along a grain boundary, and presenting a path forward toward improved models for denuded zones and localized grain boundary absorption phenomena.

36 MATERIALS SCIENCE↗

Design-Technology Co-Optimization for NVM-based Neuromorphic Processing Elements

An emerging use-case of machine learning (ML) is to train a model on a high-performance system and deploy the trained model on energy-constrained embedded systems. Neuromorphic hardware platforms, which operate on principles of the biological brain, can significantly lower the energy overhead of a machine learning inference task, making these platforms an attractive solution for embedded ML systems. In this paper, we present a design-technology tradeoff analysis to implement such inference tasks on the processing elements (PEs) of a Non-Volatile Memory (NVM)-based neuromorphic hardware. Through detailed circuit-level simulations at scaled process technology nodes, we show the negative impact of technology scaling on the information-processing latency, which impacts the quality-of-service (QoS) of an embedded ML system. At a finer granularity, the latency inside a PE depends on 1) the delay introduced by parasitic components on its current paths, and 2) the varying delay to sense different resistance states of its NVM cells. Based on these two observations, we make the following three contributions. First, on the technology front, we propose an optimization scheme where the NVM resistance state that takes the longest time to sense is set on current paths having the least delay, and vice versa, reducing the average PE latency, which improves the QoS. Second, on the architecture front, we introduce isolation transistors within each PE to partition it into regions that can be individually power-gated, reducing both latency and energy. Finally, on the system-software front, we propose a mechanism to leverage the proposed technological and architectural enhancements when implementing a machine-learning inference task on neuromorphic PEs of the hardware. Evaluations with a recent neuromorphic hardware architecture show that our proposed design-technology co-optimization approach improves both performance and energy efficiency of machine-learning inference tasks without incurring high cost-per-bit.

42 ENGINEERING↗

Drizzle, Turbulence, and Density Currents Below Post Cold Frontal Open Cellular Marine Stratocumulus Clouds

Nine cases of postcold frontal marine stratocumulus clouds exhibiting open cellular mesoscale organization are analyzed to characterize the drizzle, turbulence, and density currents below them. Data collected by the vertically pointing Doppler cloud radar and multiple lidars part of the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site are used in these analyses. A total of 251 drizzle shafts passed over the site with 76 shafts sampled for more than 15 min by the vertically pointing instruments. On average the drizzle shafts were ~18 km wide with below-cloud drizzle water path of ~40 g m −2 , cloud base rain rate of 7.78 mm day −1 , and cloud base drizzle modal diameter of 393 μm. The widths of the drizzle shafts did not exhibit any relationship with the cloud base rain rate, below-cloud drizzle water path, and the total water removed from the cloud. The downdrafts in the lowest 500 m within the drizzle cells strengthened with increasing cloud base rain rates, while the updrafts within drizzle cells did not exhibit this behavior. On average, the surface air density, pressure, water vapor mixing ratio, and wind speed along the background wind during the drizzle shaft were higher by 1.88 g m −3 , 8.60 Pa, 0.24 g kg −1 , and 1.22 m s −1 , respectively, from their background values. On the other hand, the surface air temperature was lower by 0.43 K during the drizzle shaft than its background value. The coincident measurements highlight the drizzle-turbulence-surface coupling in these cloud systems.

59 BASIC BIOLOGICAL SCIENCES↗