Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gradient flow”

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

Thermal relaxation and the complete set of second-order transport coefficients for the unitary Fermi gas from kinetic theory

We compute the complete set of second-order transport coefficients of the unitary Fermi gas, a dilute gas of spin-1/2 particles interacting via an 𝑠 -wave interaction tuned to infinite scattering length. The calculation is based on kinetic theory and the Chapman-Enskog method at second order in the Knudsen expansion. We take into account the exact two-body collision integral. We extend previous results on second-order coefficients related to shear stress by including terms related to heat flow and gradients of the fugacity. We confirm that the thermal relaxation time is given by the simple estimate 𝜏 𝜅 = 𝜅⁢𝑚/(𝑐 𝑃 ⁢𝑇) even if the full collision kernel is taken into account. Furthermore, 𝜅 is the thermal conductivity, 𝑚 is the mass of the particles, 𝑐𝑃 is the specific heat at constant pressure, and 𝑇 is the temperature.

Kinetic theory↗
DOE OSTI2021

Coupling between Alfven wave and Kelvin-Helmholtz waves in the low latitude boundary layer

The Kelvin-Helmholtz (KH) instability of magnetohydrodynamic surface waves at the low latitude boundary layer is examined using both an eigenfrequency analysis and a time-dependent wave simulation. The analysis includes the effects of sheared flow and Alfven velocity gradient. When the magnetosheath flows are perpendicular to the ambient magnetic field direction, unstable KH waves that propagate obliquely to the sheared flow direction occur at the sheared flow surface when the Alfv\'en Mach number is higher than an instability threshold. Including a shear transition layer between the magnetosphere and magnetosheath leads to secondary KH waves (driven by the sheared flow) that are coupled to the resonant surface Alfven wave. There are remarkable differences between the primary and the secondary KH waves including wave frequency, the growth rate, and the ratio between transverse and the compressional component. The secondary KH wave energy is concentrated near the shear Alfven wave frequency at the magnetosheath with a lower frequency than the primary KH waves. Although the growth rate of the secondary KH waves is lower than the primary KH waves, the threshold condition is lower, so it is expected that these types of waves will dominate at lower Mach number. Because the transverse component of the secondary KH waves is stronger than the primary KH waves, more efficient wave energy transfer from the boundary layer to the inner magnetosphere is also predicted.

Alfven wave↗
DOE OSTI2022

Coupling Between Alfvén Wave and Kelvin–Helmholtz Waves in the Low Latitude Boundary Layer

The Kelvin–Helmholtz (KH) instability of magnetohydrodynamic surface waves at the low latitude boundary layer is examined using both an eigenfrequency analysis and a time-dependent wave simulation. The analysis includes the effects of sheared flow and Alfvén velocity gradient. When the magnetosheath flows are perpendicular to the ambient magnetic field direction, unstable KH waves that propagate obliquely to the sheared flow direction occur at the sheared flow surface when the Alfvén Mach number is higher than an instability threshold. Including a shear transition layer between the magnetosphere and magnetosheath leads to secondary KH waves (driven by the sheared flow) that are coupled to the resonant surface Alfvén wave. There are remarkable differences between the primary and the secondary KH waves, including wave frequency, the growth rate, and the ratio between the transverse and compressional components. The secondary KH wave energy is concentrated near the shear Alfvén wave frequency at the magnetosheath with a lower frequency than the primary KH waves. Although the growth rate of the secondary KH waves is lower than the primary KH waves, the threshold condition is lower, so it is expected that these types of waves will dominate at a lower Mach number. Because the transverse component of the secondary KH waves is stronger than that of the primary KH waves, more efficient wave energy transfer from the boundary layer to the inner magnetosphere is also predicted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Electron Injection via Modified Diffusive Shock Acceleration in High-Mach-number Collisionless Shocks

The ability of collisionless shocks to efficiently accelerate nonthermal electrons via diffusive shock acceleration (DSA) is thought to require an injection mechanism capable of preaccelerating electrons to high enough energy where they can start crossing the shock front potential. We propose, and show via fully kinetic plasma simulations, that in high-Mach-number shocks electrons can be effectively injected by scattering in kinetic-scale magnetic turbulence produced near the shock transition by the ion Weibel, or current filamentation, instability. We describe this process as a modified DSA mechanism where initially thermal electrons experience the flow velocity gradient in the shock transition and are accelerated via a first-order Fermi process as they scatter back and forth. The electron energization rate, diffusion coefficient, and acceleration time obtained in the model are consistent with particle-in-cell simulations and with the results of recent laboratory experiments where nonthermal electron acceleration was observed. This injection model represents a natural extension of DSA and could account for electron injection in high-Mach-number astrophysical shocks, such as those associated with young supernova remnants and accretion shocks in galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Analysis of Thermal Physical Processes and Peculiarities of Modification of a Steel Surface by Pulsed Plasma Flows

The temperature gradient and melting depth of the surface of iron-based alloys under the action of high-temperature pulsed plasma beams have been estimated and compared with the experimental data. It has been shown that the steel surface melts under the action of a pulsed plasma flow at an energy density of 15–20 J/cm{sup 2}; the heat front propagates to a depth of up to 20 μm, the thickness of the molten layer being less than 10 μm. The characteristic size of the microstructure formed as a result of thermal perturbation is estimated at 17 nm. The formation of new phases with crystallite sizes that vary in the range of 16–270 nm is demonstrated experimentally. It has been shown that the formation of nanosize crystalline structure and modified near-surface region are the main factors responsible for strengthening steels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

IMPROVED PERFORMANCES IN SUBSONIC FLOWS OF AN SPH SCHEME WITH GRADIENTS ESTIMATED USING AN INTEGRAL APPROACH

In this paper, we present results from a series of hydrodynamical tests aimed at validating the performance of a smoothed particle hydrodynamics (SPH) formulation in which gradients are derived from an integral approach. We specifically investigate the code behavior with subsonic flows, where it is well known that zeroth-order inconsistencies present in standard SPH make it particularly problematic to correctly model the fluid dynamics. In particular, we consider the Gresho–Chan vortex problem, the growth of Kelvin–Helmholtz instabilities, the statistics of driven subsonic turbulence and the cold Keplerian disk problem. We compare simulation results for the different tests with those obtained, for the same initial conditions, using standard SPH. We also compare the results with the corresponding ones obtained previously with other numerical methods, such as codes based on a moving-mesh scheme or Godunov-type Lagrangian meshless methods. We quantify code performances by introducing error norms and spectral properties of the particle distribution, in a way similar to what was done in other works. We find that the new SPH formulation exhibits strongly reduced gradient errors and outperforms standard SPH in all of the tests considered. In fact, in terms of accuracy, we find good agreement between the simulation results of the new scheme and those produced using other recently proposed numerical schemes. These findings suggest that the proposed method can be successfully applied for many astrophysical problems in which the presence of subsonic flows previously limited the use of SPH, with the new scheme now being competitive in these regimes with other numerical methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

O 2 based resonantly ionized photoemission thermometry analysis of supersonic flows

Characterization of the thermal gradients within supersonic and hypersonic flows is essential for understanding transition, turbulence, and aerodynamic heating. Developments in novel, impactful non-intrusive techniques are key for enabling flow characterizations of sufficient detail that provide experimental validation datasets for computational simulations. In this work, Resonantly Ionized Photoemission Thermometry (RIPT) signals are directly imaged using an ICCD camera to realize the techniques 1D measurement capability for the first time. The direct imaging scheme presented for oxygen-based RIPT (O 2 RIPT) uses the previously established calibration data to direct excite various resonant rotational peaks within the S-branch of the C 3 Π, ( v = 2) ← X 3 Σ( v ′ = 0) absorption band of O 2 . The efficient ionization of O 2 liberates electrons that induce electron avalanche ionization of local N 2 molecules generating N 2 + , which primarily deexcites via photoemissions of the first negative band of N 2 + ( B 2 Σ <#comment/> u + − <#comment/> X 2 Σ <#comment/> g + ) . When sufficient lasing energy is used, the ionization region and subsequent photoemission signal is achieved along a 1D line thus, if directly imaged can allow for gas temperature assignments along said line; demonstrated here of up to five centimeters in length. The temperature gradients present within the ensuing shock train of a supersonic under expanded free jet serves as a basis of characterization for this new RIPT imaging scheme. The O 2 RIPT results are extensively compared and validated against well-known and established techniques (i.e., CARS and CFD). The direct imaging capability fully realizes the technique’s fundamental potential and is expected to be the standard of implementation going forward. The direct imaging capability can play instrumental roles in future scientific studies that rely upon acute characterization of thermal gradients within a medium that cannot be easily resolved by a point. Furthermore, the removal of the spectrometer greatly reduces the cost, complexity, and optical alignment associated with prior RIPT measurements.

42 ENGINEERING↗
DOE OSTI2016

A Radiotracer Assessment of Pulsatile Flow

Positron Emission Particle Tracking (PEPT) is a Lagrangian measurement technique used to examine flow behavior. In PEPT, tracer particles are labeled with a positron-emitting radionuclide. The emitted positron annihilates with a nearby electron, producing two back-to-back gamma rays of 511 keV. These gamma rays are detected using a Positron Emission Tomography (PET) scanner. The PET scanner monitors the detection sites and corresponding time in time steps on the order of milliseconds. The detection sites of the two gamma rays are connected via lines of response (LORs), and with multiple LORs a particle location is determined. Many methods have been explored in the past to track tracer particles in the field of view (FOV) of the scanner. The initial method for single particle PEPT was established by Parker et al. This method was refined to track multiple particles by Yang et al. Another method was developed by Bickell, et al based on the 'line density' of LORs in the scanner's FOV. This method was modified by Wiggins, et al to locate and track multiple particles in the FOV without a priori information about the particles. Pulsatile flow is seen in human blood flow and hydraulic and mechanical systems. The pressure gradient of the flow has been studied in an elastic tube with varying viscosity and stiffness by Anssari-Benam et al. Velocity and Reynold stresses have been examined in a rigid tube by Adamec et al. In this research, a new pulsatile flow system is constructed to analyze 3D velocity fields in Lagrangian specifications with capability of conversion to Eulerian specification. The absolute pressure across the test section is also examined. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Limit cycle oscillations in the zonal-flow-catalyzed interactions of ion-temperature-gradient turbulence

Limit-cycle oscillations are studied for ion temperature gradient turbulence, which, in the absence of large diamagnetic (mean) shear flows, saturates through energy transfer from unstable modes to large-scale stable modes via zonal-flow intermediary modes. Oscillations of zonal flow and turbulence levels are strongly constrained by the reactive, largely non-dissipative character of the zonal flows. Since existing predator–prey models for observed oscillations in experiments do not include energy transfer through zonal flows to stable modes, low-order fluid models with this physics are constructed and investigated. A simple three-wave truncation produces low-amplitude zonal flows that slowly oscillate around a zero mean, with turbulence oscillations between coupled wavenumbers that exceed linear frequencies by orders of magnitude. This inconsistency with experimental observations is caused by the weak non-linear drive of zonal flows in three-wave systems and the lack of multiple-wavenumber turbulent interactions. A more comprehensive model that preserves multiple wavenumber interactions within the context of conservative zonal-flow-mediated energy transfer to stable modes accurately reflects observed dynamics when the phase between stable and unstable modes is occasionally randomized.

Li, P. -Y. (ORCID:0000000295254171)↗
DOE OSTI2022

Zonal shear layer collapse and the power scaling of the density limit: old L-H wine in new bottles

Edge shear layer collapse causes edge cooling and aggravates radiative effects. This paper details on the microscopic dynamics of the emergence of power (Q) scaling of density limit (DL) from the shear layer collapse transport bifurcation scenario. The analysis is based on a novel 4-field model, which evolves turbulence energy, zonal flow energy, temperature gradient and density, including the neoclassical screening of zonal flow response. Bifurcation analysis yields power scaling of critical density for shear layer collapse as $n_{crit}\sim Q^{1/3}$. The favorable Q scaling of the DL emerges from the fact that the shear layer strength increases with Q, thus preventing shear layer collapse. This in turn reduces particle transport and improves particle confinement. RMP induced ambient stochastic fields degrade the shear layer by inducing decoherence in the Reynolds stress. As a result the particle transport increases and particle confinement degrades. This leads to the emergence of unfavorable stochastic field intensity ($b_{st}^{2}$) scaling of the critical density as $n_{crit}\sim(1+b_{st}^{2})^{-5/3}$. All fields, including zonal flow shear, exhibit hysteresis when the power (Q) is ramped cyclically across the bifurcation point. The hysteresis is due to dynamical delay in bifurcation on account of critical slowing down. Furthermore, the dynamical hysteresis here is fundamentally different from the hysteresis associated with the existence of bi-stable states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A configurational force for adaptive re-meshing of gradient-enhanced poromechanics problems with history-dependent variables

Here, we introduce a mesh-adaption framework that employs a multi-physical configurational force and Lie algebra to capture multiphysical responses of fluid-infiltrating geological materials while maintaining the efficiency of the computational models. To resolve sharp gradients of both displacement and pore pressure, we introduce an energy-estimate-free re-meshing criterion by extending the configurational force theory to consider the energy dissipation due to the fluid diffusion and the gradient-dependent plastic flow. To establish new equilibria after remeshing, the local tensorial history-dependent variables at the integration points are first decomposed into spectral forms. Then, the principal values and directions are projected onto smooth fields interpolated by the basis function of the finite element space via the Lie-algebra mapping. Our numerical results indicate that this Lie algebra operator in general leads to a new trial state closer to the equilibrium than the ones obtained from the tensor component mapping approach. A new configurational force for dissipative fluid-infiltrating porous materials that exhibit gradient-dependent plastic flow is introduced such that the remeshing may accommodate the need to resolve the sharp pressure gradient as well as the strain localization. The predicted responses are found to be not influenced by the mesh size due to the micromorphic regularization, while the adaptive meshing enables us to capture the width of deformation bands without the necessity of employing fine mesh everywhere in the domain.

42 ENGINEERING↗
DOE OSTI2018

Distributed gradient descent for solving optimal power flow in radial networks

Node controllers and power distribution networks in accordance with embodiments of the invention enable distributed power control. One embodiment includes a node controller comprising a memory containing: a plurality of node operating parameters; and a plurality of node operating parameters describing operating parameters for a set of at least one node selected from the group consisting of at least one downstream node and at least one upstream node; wherein the processor is configured by the node controller application to: receive and store in memory a plurality of coordinator parameters describing operating parameters of a node coordinator by the network interface; and calculate updated node operating parameters using an iterative gradient projection process to determine updated node parameters using node operating parameters that describe operating parameters of node and operating parameters of the set of at least one node, where each iteration is determined by the coordinator parameters.

Gan, Lingwen↗
DOE OSTI2018

The plastic flow machining: A new SPD process for producing metal sheets with gradient structures

Highlights: • A new process named Plastic Flow Machining is presented to produce thin sheets by an extrusion operation from bulk material. • The produced sheets present ultra-fine grain structures, good mechanical properties and high formability. • A deformation gradient was found across the thickness of the sheet regarding plastic strain, microstructure and texture. - Abstract: A new severe plastic deformation process is presented – the plastic flow machining – an extrusion process during which a sheet is formed through a lateral opening of the die in a single operation. The sheet is heavily deformed presenting an ultra-fine grain structure. Experimental results are presented for commercially pure aluminum concerning microstructures, textures and mechanical properties. A gradient structure was found across the thickness of the sheet concerning plastic strain, grain size and texture. Only simple shear textures were seen. The material strength increased by a factor of three with significant ductility and with relatively high Lankford parameter.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A micromorphically regularized Cam-clay model for capturing size-dependent anisotropy of geomaterials

Here we introduce a regularized anisotropic modified Cam-clay (MCC) model which captures the size-dependent anisotropic elastoplastic responses for clay, mudstone, shales, and sedimentary rock. By homogenizing the multiscale anisotropic effects induced by clay particle aggregate, clusters, peds, micro-fabric, and mineral contact across length scales, we introduce two distinctive anisotropic mechanisms for the MCC model at the material point and mesoscale levels. We first employ a mapping that links the anisotropic stress state to a fictitious isotropic principal stress-space to introduce anisotropy at the material point scale. Then, the mesoscale anisotropy is introduced via an anisotropic regularization mechanism. This anisotropic regularization mechanism is triggered by introducing gradient-dependence of the internal variables through a penalty method such that the resultant gradient-enhanced plastic flow may exhibit anisotropic responses non-coaxial to the stress gradient of the yield function. The influence of the size-dependent anisotropy on the formation of the shear band and the macroscopic responses of the effective media are analyzed in 2D and 3D numerical examples.

42 ENGINEERING↗
DOE OSTI2022

Observing flow of He II with unsupervised machine learning

Abstract Time dependent observations of point-to-point correlations of the velocity vector field (structure functions) are necessary to model and understand fluid flow around complex objects. Using thermal gradients, we observed fluid flow by recording fluorescence of $${\text{He}}_{2}^{*}$$ He 2 ∗ excimers produced by neutron capture throughout a ~ cm 3 volume. Because the photon emitted by an excited excimer is unlikely to be recorded by the camera, the techniques of particle tracking (PTV) and particle imaging (PIV) velocimetry cannot be applied to extract information from the fluorescence of individual excimers. Therefore, we applied an unsupervised machine learning algorithm to identify light from ensembles of excimers (clusters) and then tracked the centroids of the clusters using a particle displacement determination algorithm developed for PTV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗
DOE OSTI2024

Effect of Wettability Gradient on the Scale Formation in Falling Film Flow

Mitigation of scale formation and performance degradation remains a vital challenge for falling film evaporators in various industries. In this work, an experimental study of falling film flow on a horizontal tube is conducted to investigate the effects of wettability gradients on thermal, hydraulic, and fouling behavior. It is revealed that certain hydrophobic coating patterns, such as strip, ring, and grid patterns, lead to unwetted heat transfer area, which results in decreased heat transfer compared to fully wetted plain tube. By adjusting the geometry and position of the wettability gradient, the hybrid coating demonstrates improved heat transfer performance. Based on the characteristics of horizontal tube falling film flow, impinging jet, thin film flow, and liquid retention at the tube bottom, a hybrid coating pattern is developed to improve surface wetting and mitigate the scaling coverage. It is revealed that scale deposition is regulated by wettability gradient. Crystals tend to be dense and compact in hydrophilic areas, while they appear scattered or even absent in hydrophobic regions, depending on the dimension of the hydrophobic area. While at the hydrophilic/hydrophobic boundary, a noticeable scale thickness step is observed, which raises the potential for self-cleaning. Furthermore, the balance of minimization of scaling layer coverage and maximization of wetting area requires an optimal design in coating dimensions, for which a systemic study of both flow dynamics and fouling characteristics on the falling film is necessary in the future.

42 ENGINEERING↗