Search NASA⌕ Search

SEARCH · Search NASA

Results for “Transport coefficients”

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 217 records · Page 12

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗

Higher-form symmetry and chiral transport in real-time Abelian lattice gauge theory

We study classical lattice simulations of theories of electrodynamics coupled to charged matter at finite temperature, interpreting them using the higher-form symmetry formulation of magnetohydrodynamics (MHD). We compute transport coefficients using classical Kubo formulas on the lattice and show that the properties of the simulated plasma are in complete agreement with the predictions from effective field theories. In particular, the higher-form formulation allows us to understand from hydrodynamic considerations the relaxation rate of axial charge in the chiral plasma observed in previous simulations. A key point is that the resistivity of the plasma – defined in terms of Kubo formulas for the electric field in the 1-form formulation of MHD – remains a well-defined and predictive quantity at strong electromagnetic coupling. However, the Kubo formulas used to define the conventional conductivity vanish at low frequencies due to electrodynamic fluctuations, and thus the concept of the conductivity of a gauged electric current must be interpreted with care.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamic Structure of Magnetized High Energy Density Plasmas (Final Report)

This report summarizes results of project DE-SC0022202 ``Dynamic Structure of Magnetized High Energy Density Plasmas,'' which was originally funded for the period 08/01/2021 - 7/31/2022, but was no-cost extended to 7/31/2023. The primary objectives of the work, as stated in the original proposal, were ``develop a theory to describe a novel state of magnetized high energy density plasmas, benchmark the theory using first-principles molecular dynamics simulations, and apply the results to develop a means to interpret x-ray scattering measurements in future experiments.'' These were the goals of the original project proposal, which was for a 3 year project. The project was ultimately funded for 1 year, and so the work completed focused on the goals described for only the first year of the proposed work. This included further developing a kinetic theory for strongly magnetized plasmas, applying it to compute fluid-scale (MHD) transport coefficients, and to test the predictions using molecular dynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The right conditions for high-precision dynamic temperature and heat capacity measurement via pyrometry and conductivity

The pursuit of accurate bulk temperature T under extreme conditions has been a long-standing goal of the high pressure science community, complicated by a lack of data to inform models. To reach these extremely high-pressure, high-temperature (high P − T) conditions, a combination of dynamic and heated static experiments (e.g., diamond or gem anvil cel experiments) are used. For example, in a diamond anvil cell (DAC) experiment, a sample placed in the DAC is first pressurized. Following pressurization, the sample T is increased either by heating the entire DAC (usually using resistive heating, and limited to ∼1000K) or by applying intense laser power to the sample surfaces. In a dynamic experiment, the process of pressurizing the sample also heats it. In the case of shock physics experiments, such heating is substantial, easily reaching thousands of Kelvin; in our work we have seen T ∼17000K. Most methods of measuring temperature at ambient are not compatible with experiments under these high-pressure, high-temperature conditions: thermocouples break, melt, or have conductivity properties that differ from ambient where they are calibrated; thermometers would melt; both are too slow. As a result most methods are based on non-contact techniques such as x-ray diffraction broadening, neutron scattering, or optical methods. Of these, optical methods using the visible and near-infrared region of the spectrum are the most commonly used as the sources and detectors are readily available. In the case of optical methods the optical depth, and therefore the measurement location, is limited to the surface. When a window or anvil material is used, heat flows from the sample into the window/anvil. Likewise, if the sample undergoes a change in thermodynamic state, such as expansion upon release, different T may be expected. As a result, the surface or apparent temperature T app measurement will differ from the bulk or interior temperature that is desired. This surface measurement must be related to the bulk measurement using thermal transport models and material models. While it is tempting to conclude that one should just use x-ray methods that directly probe the interior, even these methods have been shown to depend on thermal transport and material models. Regardless of the method used to create the high P − T condition, therefore, we must understand the role of thermal transport and material models upon our interpretation of the T measurement, as well as the errors and uncertainties associated with the choice of models used in the analysis. This is a substantial area of research and this paper is by no means a complete survey of the relevant sources of uncertainty. For example, we have yet to begin to address alternate transport models in a detailed manner (e.g., Tan-Ahrens), or the many models that use additional layers to approximate melting, turbulence, or epitaxial phenomena). Likewise, we have not explored the impact upon uncertainty of thermal models that use temperature-dependent thermal transport coefficients, or the wide range of material models that can be applied. Instead, this paper focuses on using one simple model, the Urtiew-Grover model, to understand the sources of error in T measurement so that we may identify how best to focus future research efforts to return the best improvements and avoid working on over-optimizing a single type of measurement. To this end, we work through some of the best and worst case scenarios for T measurement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quarkonium Polarization Kinetic Equation from Open Quantum Systems and Effective Field Theories

Recent measurements of polarization phenomena in relativistic heavy ion collisions have aroused a great interest in understanding dynamical spin evolution of the QCD matter. In particular, the spin alignment signature of J/ψ has been recently observed in Pb-Pb collisions at LHC, which may infer nontrivial spin transport of quarkonia in quark gluon plasmas. Motivated by this, we study the spin-dependent in-medium dynamics of quarkonia by using the potential nonrelativistic QCD (pNRQCD) and the open quantum system framework. By applying the Markovian approximation and Wigner transformation, we systematically derive the Boltzmann transport equation for vector quarkonia with polarization dependence in the quantum optical limit. As opposed to the previous study for the spin-independent case where the collision terms depend on chromoelectric correlators, the new kinetic equation incorporates gauge invariant correlators of chromomagnetic fields that determine the recombination and dissociation terms with polarization dependence at the order we are working in the multipole expansion. In the quantum Brownian motion limit, the Lindblad equation with new transport coefficients defined in terms of the chromomagnetic field correlators have also been derived. Our formalism is generic and valid for both weakly-coupled and strongly-coupled quark gluon plasmas. It may be further applied to study spin alignment of vector quarkonia in heavy ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Heavy Flavor Physics at the sPHENIX Experiment

The sPHENIX experiment is a state-of-the-art jet and heavy flavor physics detector, which successfully recorded its first Au + Au collision data at 200 GeV at the Relativistic Heavy Ion Collider (RHIC). sPHENIX will provide heavy flavor physics measurements at RHIC, covering an unexplored kinematic region and unprecedented precision, to probe the parton energy loss mechanism, parton transport coefficients in quark–gluon plasma, and the hadronization process under various medium conditions. At the center of sPHENIX, the monolithic active pixel sensor (MAPS)-based VerTeX detector (MVTX) is a high-precision silicon pixel detector. The MVTX provides excellent position resolution and the capability of operating in continuous streaming readout mode, allowing precise vertex determination and recording a large data sample, both of which are particularly crucial for heavy flavor physics measurements. In this work, we will show the general performance of heavy-flavor hadron reconstruction. In addition, we will discuss the commissioning experience with sPHENIX. Finally, we will provide the projection of b-hadron and jet observables and discuss the estimated constraints on theoretical models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling Neutron Star Matter in the Age of Multimessenger Astrophysics

Abstract The interpretation of available and forthcoming data obtained from multimessenger astrophysical observations—potentially providing unprecedented access to neutron star properties—will require the development of novel, accurate theoretical models of dense matter. Of great importance, in this context, will be the capability for devising a description of thermal effects applicable to the study of quantities other than the equation of state, such as the transport coefficients and the neutrino mean free path in the nuclear medium. A formalism based on correlated basis states and the cluster expansion technique has been previously employed to derive a well-behaved effective interaction—suitable for use in standard perturbation theory—from a state-of-the-art nuclear Hamiltonian, including phenomenological two- and three-nucleon potentials. Here, we provide a comprehensive and self-contained account of the extension of this approach to the treatment of finite-temperature effects, and report the results of numerical calculations of a number of properties of nuclear matter with arbitrary neutron excess and temperature up to 50 MeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Bottomonium suppression in pNRQCD and open quantum system approach

By employing the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory within an open quantum system framework, we derive a Lindblad equation governing the evolution of the heavy-quarkonium reduced density matrix, accurate to next-to-leading order (NLO) in the ratio of the state's binding energy to the medium's temperature [1]. The derived NLO Lindblad equation provides a more reliable description of heavy-quarkonium evolution in the quark-gluon plasma at low temperatures compared to the leading-order truncation. For phenomenological applications, we numerically solve this equation using the quantum trajectories algorithm. By averaging over Monte Carlo-sampled quantum jumps, we obtain solutions without truncation in the angular momentum quantum number of the considered states. Our analysis highlights the importance of quantum jumps in the nonequilibrium evolution of bottomonium states within the quark-gluon plasma [2]. Additionally, we demonstrate that the quantum regeneration of singlet states from octet configurations is essential to explain experimental observations of bottomonium suppression. The heavy-quarkonium transport coefficients used in our study align with recent lattice QCD determinations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Diffusional growth of cloud droplets in homogeneous isotropic turbulence: DNS, scaled-up DNS, and stochastic model

This paper presents a novel methodology to use direct numerical simulation (DNS) to study the impact of isotropic homogeneous turbulence on the condensational growth of cloud droplets. As shown by previous DNS studies, the impact of turbulence increases with the computational domain size, that is, with the Reynolds number, because larger eddies generate higher and longer-lasting supersaturation fluctuations that affect growth of individual cloud droplets. The traditional DNS can only simulate a limited range of scales because of the excessive computational cost that comes from resolving all scales involved, that is, from large scales at which the turbulent kinetic energy (TKE) is introduced down to the Kolmogorov microscale, and from following every single droplet. The novel approach is referred to as the “scaled-up DNS”. The scaling up is done in two parts, first by increasing both the computational domain and the Kolmogorov microscale and second by using super-droplets instead of real droplets. To ensure proper dissipation of TKE and scalar variance at small scales, molecular transport coefficients are appropriately scaled up with the grid length. For the scaled-up domains, say, meters and tens of meters, one needs to follow billions of real droplets. This is not computationally feasible, and so-called super-droplets are applied in scaled-up DNS simulations. Each super-droplet represents an ensemble of identical real droplets, and the number of real droplets represented by a super-droplet is referred to as the multiplicity attribute. After simple tests showing the validity of the methodology, scaled-up DNS simulations are conducted for five domains, the largest of 64 3 m 3 volume using a DNS of 256 3 grid points and various multiplicities. All simulations are carried out with vanishing mean vertical velocity and with no mean supersaturation, similarly to past DNS studies. As expected, the supersaturation fluctuations as well as the spread in droplet size distribution increase with the domain size, with the droplet radius variance increasing in time t as t 1∕2 as identified in previous DNS studies. Scaled-up simulations with different multiplicities document numerical convergence of the scaled-up solutions. Finally, we compare the scaled-up DNS results with a simple stochastic model that calculates supersaturation fluctuations based on the vertical velocity fluctuations updated using the Langevin equation. Overall, the results document similar scaling to previous small-domain DNS simulations and support the notion that the stochastic subgrid-scale model is a valuable tool for the multi-scale simulation of droplet spectral evolution applying a large-eddy simulation model.

54 ENVIRONMENTAL SCIENCES↗

Conservation laws and effective hadronization models

Hadronization models based on local string-breaking dynamics are typically Markovian by construction, yet the physical ensemble of final states is shaped by global constraints that couple the entire fragmentation trajectory. Recasting hadronization as a conditioned stochastic diffusion process provides a precise mathematical resolution to this tension. In particular, this language reveals explicitly that constraints stemming from conservation laws induce non-Markovian correlations between otherwise independent fragmentation steps, and that these correlations can be absorbed exactly into a renormalization of the local dynamics through a Doob $h$-transform. We develop this formalism for a $q\bar{q}$ string in the chiral limit, where the longitudinal-transverse factorization of the Lund kernel becomes exact, enabling systematic power counting and clean ultraviolet (UV)/infrared (IR) separation. The dynamics organize naturally into a tower of effective theories distinguished by the remaining string mass, spanning a UV fixed point with scale-invariant transport coefficients, an intermediate regime where transverse phase space induces controlled running, and an IR boundary layer where non-local effects enter at leading order. The tower exhibits genuine Wilsonian structure, including $β$-functions, anomalous dimensions, and systematic matching conditions. The resulting framework achieves a clean factorization of universal microscopic fragmentation dynamics from infrared constraint effects, and opens new directions for both the theoretical analysis and practical simulation of hadronization.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗

Kinetic-energy dissipation and fluctuations in strongly damped heavy-ion collisions within the stochastic mean-field approach

Microscopic mean-field approaches have been successful in describing the most probable reaction outcomes in low-energy heavy-ion reactions. However, those approaches are known to severely underestimate dispersions of observables around the average values that has limited their applicability. Recently it has been shown that a quantal transport approach based on the stochastic mean-field (SMF) theory significantly improves the description, while its application has been limited so far to fragment mass and charge dispersions. In this work, we extend the quantal transport approach based on the SMF theory for relative kinetic energy dissipation and angular momentum transfer in low-energy heavy-ion reactions. Based on the SMF concept, analytical expressions are derived for the radial and tangential friction and associated diffusion coefficients. Those quantal transport coefficients are calculated microscopically in terms of single-particle orbitals within the time-dependent Hartree-Fock (TDHF) approach. As the first application of the proposed formalism, we consider the radial linear momentum dispersion, neglecting the coupling between radial and angular momenta. We analyze the total kinetic energy (TKE) distribution of binary reaction products in the 136 Xe + 208 Pb reaction at E c.m. =526 MeV and compare with experimental data. From time evolution of single-particle orbitals in TDHF, the radial diffusion coefficient is computed on a microscopic basis, while a phenomenological treatment is introduced for the radial friction coefficient. By solving the quantal diffusion equation for the radial linear momentum, the dispersion of the radial linear momentum is obtained, from which one can construct the TKE distribution. We find that the calculations provide a good description of the TKE distribution for strongly damped events with large energy losses, TKEL ≳ 150 MeV. However, the calculations underestimate the TKE distribution for smaller energy losses. Further studies are needed to improve the technical details of calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fully generalized, turbulent trace impurity transport with Gkeyll and Flan in the DIII-D far-SOL

The Monte Carlo trace impurity turbulent transport code Flan is introduced for the first time. Flan follows impurities in a turbulent background plasma from Gkeyll using the Lorentz force to resolve the full particle gyro orbit. Collisions are handled using the Nanbu collision algorithm (Nanbu 1997 Phys. Rev. E 55 4642–52), and ionization/recombination is handled via ADAS coupling. The far-SOL of a generic DIII-D L-mode is simulated with and without the collision model to show how collisions affect radial tungsten transport. Anomalous diffusion coefficient (D r ) and pinch velocity (v p ) profiles are extracted from fits to the results. With collisions, D r and v p are between 0–1.0 m 2 s −1 and −100–100 m s −1 , respectively. Without collisions, D r and v p are between 0–0.3 m 2 s −1 and −50–50 m s −1 , respectively. Exponential fits to the radial W density profiles and experimental data from W deposition along a collector probe are in good agreement, demonstrating Flan as a useful interpretive modeling tool. Additional simulations show that impurity transport away from the wall increases with atomic number, though it is not clear why. Flan has the potential to better interpret existing data and improve reactor scale predictions of core contamination because the underlying physics model is very general and does not rely on arbitrary user-defined transport coefficients.

DIII-D↗

Novel method for accurately estimating membrane transport properties and mass transfer coefficients in reverse osmosis

Here, we present a simple and robust method to simultaneously characterize the water and salt permeability (A, B) of reverse osmosis (RO) membranes and mass transfer coefficient (k) in membrane modules. The proposed methodology comprises a set of RO experiments performed at different operating pressures or stages. The measured water and salt fluxes in each stage are simultaneously fitted to the RO transport equations by performing a non-linear regression, using A, B, and k as regression parameters. We first perform a systematic accuracy analysis of the proposed method across the full operational range of RO. The assessment shows that the method accuracy is substantially higher than current methods and increases with number of experimental stages and driving forces. This assessment is used to inform the design of an experimental protocol that minimizes errors in estimated A, B, and k. We then evaluate two commercial RO membranes following the new protocol. For both membranes, A and B parameters decrease by 17% and 15% from the dilute solution to seawater concentrations, whereas the k parameter remains constant. Our study demonstrates that the proposed method, informed by data-driven experimental designs, provides a new approach for accurately characterizing transport phenomenon in membrane processes with feeds of less than 100 g/L total dissolved solids.

36 MATERIALS SCIENCE↗

Structural, magnetic and electrical transport properties of YMnAl

The intermetallic compound YMnAl, crystallizing in the cubic C-15 (MgCu 2 -type) structure, was prepared using arc-melting and annealing. Rietveld refinement of the XRD pattern yields a lattice parameter of a = 7.816(1) Å⁠. The system exhibits antiferromagnetic order with a transition temperature of approximately 38 K, along with indications of a spin-glass–like phase. The effective magnetic moment, estimated from Curie–Weiss fitting of the susceptibility versus temperature curve, is 6.06 μ B /f.u. and the Curie–Weiss temperature of −79 K confirming the antiferromagnetic nature. Isothermal magnetization measurements below and above the Néel temperature show no saturation, with a maximum magnetization of 0.09 μ B /f.u observed at an applied field of 9 T. The temperature dependence of resistivity exhibits a negative temperature coefficient, characteristic of dirty metals. Hall effect measurements indicate the absence of an anomalous Hall contribution. The Hall coefficient RH, carrier concentration n, and mobility μ estimated from R XY (H) are -1.05×10 -8 m 3 /C, 5.92×10 26 m -3 and 1.05×10 -3 m 2 /V.s, respectively. Our first-principles calculations indicate that YMnAl in the orthorhombic (Imma) structure has the lowest energy; however, the substitutional disorder occurring during experimental synthesis may favor the formation of the cubic Laves phase.

Crystallography↗

Efficacy of the Method of Four Coefficients to Determine Charge-Carrier Scattering

The investigation of the electronic properties of semiconductors from transport measurements (i.e., resistivity, Hall, and Seebeck coefficient measurements) is challenging due to the averaging of charge-carrier dynamics inherent in such measurements. Here, we investigate the incorporation of a fourth measurement of electronic transport, the Nernst coefficient, into the analysis, termed the method of four coefficients. This approach yields the Fermi level, effective mass, scattering exponent, and relaxation time prefactor. We begin with a review of the underlying mathematics and investigate the mapping between the four-dimensional material property and transport coefficient spaces. We then investigate how the traditional single-parabolic band method yields a single, potentially incorrect point on the solution subspace. This uncertainty can be resolved through Nernst coefficient measurements and we map the span of the ensuing subspace. We conclude with an investigation of how sensitive the analysis of transport coefficients is to experimental error for different sample types.

charge carriers↗

Studying QGP with flow: A theory overview

I review recent developments in the phenomenological study of the quark-gluon plasma (QGP) transport properties based on a personal selection of results that were presented at Quark Matter 2019. Additionally, the constraints on the temperature dependence of QGP shear and bulk viscosity are summarized. I discuss new theory advancements towards more realistic 3D dynamical simulations of heavy-ion collisions at finite baryon density. The challenges and opportunities of applying hydrodynamics to small collision systems are highlighted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying selective solvent transport under an electric field in mixed-solvent electrolytes

Electrolytes in lithium-ion batteries comprise solvent mixtures, but analysis of ion transport is always based on treating the solvents as a single-entity. We combine electrophoretic NMR (eNMR) measurements and molecular dynamics (MD) simulations to quantify electric-field-induced transport in a concentrated solution containing LiPF 6 salt dissolved in an ethylene carbonate/ethyl methyl carbonate (EC/EMC) mixture. The selective transport of EC relative to EMC is reflected in the difference between two transference numbers, defined as the fraction of current carried by cations relative to the velocity of each solvent species. This difference arises from the preferential solvation of cations by EC and its dynamic consequences. The simulations reveal the presence of a large variety of transient solvent-containing clusters which migrate at different velocities. Rigorous averaging over different solvation environments is essential for comparing simulated and measured transference numbers. Our study emphasizes the necessity of acknowledging the presence of four species in mixed-solvent electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗