Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transport phenomena”

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

Quantum dynamics simulation of the advection-diffusion equation

The advection-diffusion equation is simulated via several quantum algorithms. Three formulations are considered: (1) Trotterization, (2) variational quantum time evolution (VarQTE), and (3) adaptive variational quantum dynamics simulation (AVQDS). These schemes were originally developed for the Hamiltonian simulation of many-body quantum systems. The finite-difference discretized operator of the transport equation is formulated as a Hamiltonian and solved without the need for ancillary qubits. Computations are conducted on a quantum simulator (IBM Qiskit Aer) and a superconducting quantum hardware (IBM Fez). The former emulates the latter without the noise. The actual hardware implementation experiences significant noise. The results of the quantum simulator are compared with data from direct numerical simulation (DNS) with infidelities of the order 10 −5 . In the quantum simulator, Trotterization is observed to have the lowest infidelity and is suitable for fault-tolerant computation. The AVQDS algorithm requires the lowest gate count and circuit depth. The VarQTE algorithm is the next best in terms of gate counts, but the number of its optimization variables is directly proportional to the number of qubits. Due to current hardware limitations, Trotterization cannot be implemented, as it has an overwhelmingly large number of operations. Meanwhile, AVQDS and VarQTE can be executed at the hardware level. These algorithms present a new paradigm for computational transport phenomena on quantum computers.

Alipanah, Hirad [Univ. of Pittsburgh, PA (United S↗
DOE OSTI2020

Anomalous phonon-mode dependence in polarized Raman spectroscopy of the topological Weyl semimetal TaP

Topological Weyl semimetals (WSMs) have attracted widespread interest due to the chiral Weyl fermions and surface Fermi arcs that enable unique optical and transport phenomena. In this work, we present angle-resolved Raman spectroscopy of TaP, a prototypical noncentrosymmetric WSM, for five excitation wavelengths ranging from 364 to 785 nm. The Raman-active modes, $A_1$, $B^1_1$, and $B^2_1$ modes, exhibit two main unique features beyond the conventional Raman theory. First, the relative intensities of Raman-active modes change as a function of the excitation wavelength. Second, angle-resolved polarized Raman spectra show systematic deviation from the Raman tensor theory. In particular, the $B^1_1$ mode is absent for 633-nm excitation, whereas the $B^2_1$ mode shows an unusual twofold symmetry instead of a fourfold symmetry for 488-, 532-, and 633-nm excitations. Finally, these unconventional phenomena are attributed to the interference effect in the Raman process owing to the existence of multiple carrier pockets with almost the same energy but different symmetries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Competition between magnetic order and charge localization in Na 2 IrO 3 thin crystal devices thin crystal devices

Spin orbit assisted Mott insulators such as sodium iridate (Na 2 IrO 3 ) have been an important subject of study in recent years. In these materials, the interplay of electronic correlations, spin-orbit coupling, crystal field effects, and a honeycomb arrangement of ions bring exciting ground states, predicted in the frame of the Kitaev model. The insulating character of Na 2 IrO 3 has hampered its integration to an electronic device, desirable for applications, such as the manipulation of quasiparticles interesting for topological quantum computing. Here we show through electronic transport measurements supported by angle-resolved photoemission spectroscopy (ARPES) experiments, that electronic transport in Na 2 IrO 3 is ruled by variable range hopping and it is strongly dependent on the magnetic ordering transition known for bulk Na 2 IrO 3 , as well as on external electric fields. Furthermore, electronic transport measurements allow us to deduce a value for the localization length and the density of states in our Na 2 IrO 3 thin crystal devices, and offer an alternative approach to study insulating 2D-materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Magnetotransport properties of the topological nodal-line semimetal CaCdSn

Topological nodal-line semimetals support protected band crossings, which form nodal lines or nodal loops between the valence and conduction bands and exhibit novel transport phenomena. Here we address the topological state of the nodal-line semimetal candidate material CaCdSn, and report magnetotransport properties of its single crystals grown by the self-flux method. In this work, our first-principles calculations show that the electronic structure of CaCdSn harbors a single nodal loop around the Γ point in the absence of spin-orbit coupling effects. The nodal crossings in CaCdSn are found to lie above the Fermi level and yield a Fermi surface that consists of both electron and hole pockets. CaCdSn exhibits high mobility (μ ≈ 3.44 × 10 4 cm 2 V –1 s –1 ) and displays a field-induced metal-semiconductor-like crossover with a plateau in resistivity at low temperature. We observe an extremely large and quasilinear nonsaturating transverse as well as longitudinal magnetoresistance (MR) at low temperatures (≈ 7.44 × 10 3 % and ≈ 1.71 × 10 3 %, respectively, at 4 K). We also briefly discuss possible reasons behind such a large quasilinear magnetoresistance and its connection with the nontrivial band structure of CaCdSn.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

c -axis transport in UTe 2 : Evidence of three-dimensional conductivity component

Here we study the temperature dependence of electrical resistivity for currents directed along all crystallographic axes of the spin-triplet superconductor UTe 2 . We focus particularly on an accurate determination of the resistivity along the c axis (ρ c ) by using a generalized Montgomery technique that allows extraction of crystallographic resistivity components from a single sample. In contrast to expectations from the observed highly anisotropic band structure, our measurement of the absolute values of resistivities in all current directions reveals a surprisingly nearly isotropic transport behavior at temperatures above Kondo coherence, with ρ c ~ ρ b ~ 2 ρ a , that evolves to reveal qualitatively distinct behaviors on cooling. The temperature dependence of ρ c exhibits a peak at a temperature much lower than the onset of Kondo coherence observed in ρ a and ρ b , consistent with features in magnetotransport and magnetization that point to a magnetic origin. A comparison to the temperature-dependent evolution of the scattering rate observed in angle-resolved photoemission spectroscopy experiments provides important insights into the underlying electronic structure necessary for building a microscopic model of superconductivity in UTe 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2020

Magnetic Damping Modulation in IrMn 3 / Ni 80 Fe 20 via the Magnetic Spin Hall Effect

Noncollinear antiferromagnets can have additional spin Hall effects due to the net chirality of their magnetic spin structure, which provides for more complex spin-transport phenomena compared to ordinary nonmagnetic materials. Here we investigated how ferromagnetic resonance of permalloy (Ni 80 Fe 20 ) is modulated by spin Hall effects in adjacent epitaxial IrMn 3 films. We observe a large dc modulation of the ferromagnetic resonance linewidth for currents applied along the [001] IrMn 3 direction. Furthermore, this very strong angular dependence of spin-orbit torques from dc currents through the bilayers can be explained by the magnetic spin Hall effect where IrMn 3 provides novel pathways for modulating magnetization dynamics electrically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Record-quality GaAs two-dimensional hole systems

The complex band structure, large spin-orbit induced band splitting, and heavy effective mass of two-dimensional (2D) hole systems hosted in GaAs quantum wells render them rich platforms to study many-body physics and ballistic transport phenomena. Here we report ultra-high-quality (001) GaAs 2D hole systems, fabricated using molecular beam epitaxy and modulation doping, with mobility values as high as 5.8 × 10 6 cm 2 /(V s) at a hole density of p = 1.3 × 10 11 /cm 2 , implying a mean free path of ≃ 27 μm. In the low-temperature magnetoresistance trace of this sample, we observe high-order fractional quantum Hall states up to the Landau level filling ν = 12/25 near ν = 1/2. Furthermore, we see a deep minimum develop at ν = 1/5 in the magnetoresistance of a sample with a much lower hole density of p = 4.0 × 10 10 /cm 2 where we measure a mobility of 3.6 × 10 6 cm2/(V s). Finally, these improvements in sample quality were achieved by the reduction of residual impurities both in the GaAs channel and in the AlGaAs barrier material, as well as optimization in the design of the sample structure.

2-dimensional systems↗
DOE OSTI2022

A room temperature polar magnetic metal

The emergence of long-range magnetic order in noncentrosymmetric compounds has stimulated interest in the possibility of exotic spin transport phenomena and topologically protected spin textures for applications in next-generation spintronics. Polar magnets, with broken symmetries of spatial inversion and time reversal, usually host chiral spin textures. This work reports on a wurtzite-structure polar magnetic metal, identified as AA'-stacked (Fe 0.5 Co 0.5 ) 5 GeTe 2 , which exhibits a Néel-type skyrmion lattice as well as a Rashba-Edelstein effect at room temperature. Atomic resolution imaging of the structure reveals a structural transition as a function of Co-substitution, leading to the emergence of the polar phase at 50% Co. We report this discovery reveals an unprecedented layered polar magnetic system for investigating intriguing spin topologies, and it ushers in a promising new framework for spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Carrier properties of Bi(111) grown on mica and Si(111)

A comparison is presented between properties of Bi(111) films grown by Stranski-Krastanov epitaxy on Si(111) 7 × 7 and by van der Waals epitaxy on mica substrates. Thin film morphologies and electronic transport properties of Bi(111) films of variable thickness are investigated for each growth method. Atomic force micrographs for Bi(111) films on mica reveal clearly defined triangular regions consisting of layered steps with height 0.4 nm, corresponding to the Bi(111) bilayer height. Variable-temperature electronic transport measurements show the existence of a quantum confinement-induced energy gap in the film interiors, resulting in a semimetal-to-semiconductor transition. Magnetotransport analysis in a three-carrier model including metallic electrons in surface states and electrons and holes in the films' interiors provides a detailed study of densities, mobilities, and mean-free paths of the three carrier types. Further, improved electronic transport properties are found in Bi(111) films of higher thickness on mica compared to on Si(111), a likely result of the largely strain-free van der Waals epitaxial growth.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultrafast nonequilibrium dynamics and high-harmonic generation in two-dimensional quantum spin Hall materials

For this work, we develop the theoretical framework of nonequilibrium ultrafast photonics in monolayer quantum spin Hall insulators supporting a multitude of topological states. In these materials, ubiquitous strong light-matter interactions in the femtosecond scale lead to nonadiabatic quantum dynamics, resulting in topology-dependent nonlinear optoelectronic transport phenomena. We investigate the mechanism driving topological Dirac fermions interacting with strong ultrashort light pulses and uncover various experimentally accessible physical quantities that encode fingerprints of the quantum material's topological electronic state from the high-harmonic generated spectrum. Our work sets the theoretical cornerstones to realize the full potential of time-resolved harmonic spectroscopy for understanding nonequilibrium processes in quantum topological systems and identifying topological invariants in two-dimensional quantum spin Hall solid state systems.

2-dimensional systems↗
DOE OSTI2025

Potential quantum advantage for simulation of fluid dynamics

Numerical simulation of turbulent fluid dynamics needs to either parametrize turbulence—which introduces large uncertainties—or explicitly resolve the smallest scales—which is prohibitively expensive. Here, we provide evidence through analytic bounds and numerical studies that a potential quantum speedup can be achieved to simulate fluid dynamics using quantum computing. Specifically, we provide a lattice Boltzmann formulation of fluid dynamics for which we give evidence that low-order Carleman linearization is much more accurate than previously believed for these systems. This is achieved via a combination of reformulating the Navier-Stokes nonlinearity (u·$\triangledown$u) to lattice-Boltzmann nonlinearity (u 2 ) and accurately linearizing the dynamical equations, which effectively trades nonlinearity for additional degrees of freedom that add negligible expense in the quantum solver. Based on this, we apply a quantum algorithm for simulating the Carleman-linearized lattice Boltzmann equation and provide evidence that its cost scales logarithmically with system size compared with polynomial scaling in the best known classical algorithms. In this paper, we suggest that a quantum advantage may exist for simulating fluid dynamics, paving the way for simulating nonlinear multiscale transport phenomena in a wide range of disciplines using quantum computing.

42 ENGINEERING↗
DOE OSTI2025

Heat conduction in magnetic insulators via hybridization of acoustic phonons and spin-flip excitations

We present a comprehensive study on the longitudinal magnetothermal transport in a paramagnetic effective spin-1/2 magnetic insulator CsYbSe 2 , by introducing a minimal model requiring only Zeeman splitting and magnetoelastic coupling. We use it to argue that hybridized excitations—formed from acoustic phonons and localized spin-flip-excitations across the Zeeman gap of the crystal electric field ground doublet—are responsible for a nonmonotonic field dependence of longitudinal thermal conductivity. Beyond highlighting a starring role for phonons, our results raise the prospect of universal magnetothermal transport phenomena in paramagnetic insulators that originate from simple features shared across many systems. Published by the American Physical Society 2025

Pocs, Christopher A.↗
DOE OSTI2025

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗
DOE OSTI2022

Real-time deep-learning inversion of seismic full waveform data for CO 2 saturation and uncertainty in geological carbon storage monitoring

Deep-learning inversion has recently drawn attention in geological carbon storage research due to its potential of imaging and monitoring carbon storage in real time, significantly improving efficiency and safety of carbon storage operations. We present a deep-learning full waveform inversion method that after the neural network has been trained can image CO 2 saturation and its uncertainty in real time. Our deep-learning inversion method is based on the U-Net architecture with the neural network trained on pairs of synthetic seismic data and CO 2 saturation models. Accordingly, our training establishes a mapping relationship between seismic data and CO 2 saturation models and once fully trained directly estimates CO 2 saturation as a function of subsurface location. We further quantify uncertainties of CO 2 saturation estimates using the Monte Carlo dropout method and a bootstrap aggregating method. For this proof-of-concept study, the CO 2 training models and data are derived from the Kimberlina 1.2 model, a hypothetical 3D geological carbon storage model that is constructed based on various geological and hydrological data from the Southern San Joaquin Basin, California. We perform deep-learning inversion experiments using noise-free and noisy training and test data sets and compare the results. Our modelling experiments show that (1) the deep-learning inversion can estimate 2D distributions of CO 2 fairly well even in the presence of Gaussian random noise and (2) both CO 2 saturation imaging and uncertainty quantification can be done in real time. Our results suggest that the deep-learning inversion method can serve as a robust real-time monitoring tool for geological carbon storage and/or other time-varying reservoir/aquifer properties that result from injection, extraction, and/or other subsurface transport phenomena.

58 GEOSCIENCES↗
DOE OSTI2020

Mesoscale Effects of Composition and Calendering in Lithium-Ion Battery Composite Electrodes

Macrohomogeneous battery models are widely used to predict battery performance, necessarily relying on effective electrode properties, such as specific surface area, tortuosity, and electrical conductivity. While these properties are typically estimated using ideal effective medium theories, in practice they exhibit highly non-ideal behaviors arising from their complex mesostructures. In this paper, we computationally reconstruct electrodes from X-ray computed tomography of 16 nickel–manganese–cobalt-oxide electrodes, manufactured using various material recipes and calendering pressures. Due to imaging limitations, a synthetic conductive binder domain (CBD) consisting of binder and conductive carbon is added to the reconstructions using a binder bridge algorithm. Reconstructed particle surface areas are significantly smaller than standard approximations predicted, as the majority of the particle surface area is covered by CBD, affecting electrochemical reaction availability. Finite element effective property simulations are performed on 320 large electrode subdomains to analyze trends and heterogeneity across the electrodes. Significant anisotropy of up to 27% in tortuosity and 47% in effective conductivity is observed. Electrical conductivity increases up to 7.5× with particle lithiation. We compare the results to traditional Bruggeman approximations and offer improved alternatives for use in cell-scale modeling, with Bruggeman exponents ranging from 1.62 to 1.72 rather than the theoretical value of 1.5. We also conclude that the CBD phase alone, rather than the entire solid phase, should be used to estimate effective electronic conductivity. This study provides insight into mesoscale transport phenomena and results in improved effective property approximations founded on realistic, image-based morphologies.

25 ENERGY STORAGE↗
DOE OSTI2024

Coupled Multiphysics Modeling of Lithium-Ion Batteries for Automotive Crashworthiness Applications

Considerable advances have been made in battery safety models, but achieving predictive accuracy across a wide range of conditions continues to be challenging. Interactions between dynamically evolving mechanical, electrical, and thermal state variables make model prediction difficult during mechanical abuse scenarios. In this study, we develop a physics-based modeling approach that allows for choosing between different mechanical and electrochemical models depending on the required level of analysis. We demonstrate the use of this approach to connect cell-level abuse response to electrode-level and particle-level transport phenomena. A pseudo-two-dimensional model and simplified single-particle models are calibrated to electrical-thermal cycling data and applied to mechanically induced short-circuit scenarios to understand how the choice of electrochemical model affects the model prediction under abuse scenarios. These models are implemented using user-defined subroutines on ls-dyna finite element software and can be coupled with existing automotive crash safety models.

analysis and design of components↗