Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluids”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Vortex splitting in two-dimensional fluids and non-neutral electron plasmas with smooth vorticity profiles

Initially, elliptical, quasi-two-dimensional (2D) fluid vortices can split into multiple pieces if the aspect ratio is sufficiently large due to the growth and saturation of perturbations known as Love modes on the vortex edge. Presented here are experiments and numerical simulations, showing that the aspect ratio threshold for vortex splitting is significantly higher for vortices with realistic, smooth edges than that predicted by a simple “vortex patch” model, where the vorticity is treated as piecewise constant inside a deformable boundary. The experiments are conducted by exploiting the E × B drift dynamics of collisionless, pure electron plasmas in a Penning–Malmberg trap, which closely model 2D vortex dynamics due to an isomorphism between the Drift–Poisson equations describing the plasmas and the Euler equations describing ideal fluids. The simulations use a particle-in-cell method to model the evolution of a set of point vortices. The aspect ratio splitting threshold ranges up to about twice as large as the vortex patch prediction and depends on the edge vorticity gradient. This is thought to be due to spatial Landau damping, which decreases the vortex aspect ratio over time and, thus, stabilizes the Love modes. Near the threshold, asymmetric splitting events are observed in which one of the split products contains much less circulation than the other. These results are relevant to a wide range of quasi-2D fluid systems, including geophysical fluids, astrophysical disks, and drift-wave eddies in tokamak plasmas.

Physics↗
DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The smallest fluid on earth

High energy heavy ion collisions create quark gluon plasmas that behave like almost perfect fluids. Very similar features to those that led to this insight have also been observed in experimental data from collisions of small systems, involving protons or other light nuclei. Here we describe recent developments aimed at understanding whether, and if so how, systems that produce relatively few particles (orders of magnitude less than in typical heavy ion collisions) and are only one to a few times the size of a proton, can behave like fluids. This involves a deeper understanding of fluid dynamics and its applicability, improvements of our understanding of the initial geometry of the collisions by considering fluctuations of the proton shape, as well as advancements in the calculation of initial state effects within an effective theory of quantum chromodynamics, which can affect the observables that are used to study fluid behavior. We further address open questions and discuss future directions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Relativistic analogue of the Newtonian fluid energy equation with nucleosynthesis

In Newtonian fluid dynamics simulations in which composition has been tracked by a nuclear reaction network, energy generation due to composition changes has generally been handled as a separate source term in the energy equation. Here, a relativistic equation in conservative form for total fluid energy, obtained from the spacetime divergence of the stress-energy tensor, in principle encompasses such energy generation; but it is not explicitly manifest. An alternative relativistic energy equation in conservative form—in which the nuclear energy generation appears explicitly, and that reduces directly to the Newtonian internal+kinetic energy in the appropriate limit—emerges naturally and self-consistently from the difference of the equation for total fluid energy and the equation for baryon number conservation multiplied by the average baryon mass m, when m is expressed in terms of contributions from the nuclear species in the fluid, and allowed to be mutable.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Field theories and fluids for an interacting dark sector

We consider the relationship between fluid models of an interacting dark sector and the field theoretical models that underlie such descriptions. This question is particularly important in light of suggestions that such interactions may help alleviate a number of current tensions between different cosmological datasets. We construct consistent field theory models for an interacting dark sector that behave exactly like the coupled fluid ones, even at the level of linear perturbations, and can be trusted deep in the nonlinear regime. As a specific example, we focus on the case of a Dirac, Born-Infeld (DBI) field conformally coupled to a quintessence field. We show that the fluid linear regime breaks before the field gradients become large; this means that the field theory is valid inside a large region of the fluid nonlinear regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Fluid flow through packings of elastic shells

Fluid transport in porous materials is commonly studied in geological samples (soil, sediments, etc.) or idealized systems, but the fluid flow through compacted granular materials, consisting of substantially strained granules, remains relatively unexplored. As a step toward filling this gap, we study a model of liquid transport in packings of deformable elastic shells using finite-element and lattice-Boltzmann methods. In this work, we find that the fluid flow abruptly vanishes as the porosity of the material falls below a critical value, and the flow obstruction exhibits features of a percolation transition. We further show that the fluid flow can be captured by a simplified permeability model in which the complex porous material is replaced by a collection of disordered capillaries, which are distributed and shaped by the percolation transition. To that end, we numerically explore the divergence of hydraulic tortuosity $τ_H$ and the decrease of a hydraulic radius Rh as the percolation threshold is approached. We interpret our results in terms of scaling predictions derived from the percolation theory applied to random packings of spheres.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Universal molecular-kinetic scaling relation for slip of a simple fluid at a solid boundary

Using the observation that slip in simple fluids at low and moderate shear rates is a thermally activated process driven by the shear stress in the fluid close to the solid boundary, we develop a molecular-kinetic model for simple fluid slip at solid boundaries. Here, the proposed model, which is in the form of a universal scaling relation that connects slip and shear rate, reduces to the well-known Navier slip condition under low shear conditions, providing a direct connection between molecular parameters and the slip length. Molecular-dynamics simulations are in very good agreement with the predicted dependence of slip on system parameters, including the temperature and fluid-solid interaction strength. Connections between our model and previous work, as well as simulation and experimental results, are explored and discussed.

42 ENGINEERING↗
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 OSTI2016

Density Fluctuations of Hard-Sphere Fluids in Narrow Confinement

Spatial confinement induces microscopic ordering of fluids, which in turn alters many of their dynamic and thermodynamic properties. However, the isothermal compressibility has hitherto been largely overlooked in the literature, despite its obvious connection to the underlying microscopic structure and density fluctuations in confined geometries. Here, we address this issue by probing density profiles and structure factors of hard-sphere fluids in various narrow slits, using x-ray scattering from colloid-filled nanofluidic containers and integral-equation-based statistical mechanics at the level of pair distributions for inhomogeneous fluids. Most importantly, we demonstrate that density fluctuations and isothermal compressibilities in confined fluids can be obtained experimentally from the long-wavelength limit of the structure factor, providing a formally exact and experimentally accessible connection between microscopic structure and macroscopic, thermodynamic properties. Our approach will thus, for example, allow direct experimental verification of theoretically predicted enhanced density fluctuations in liquids near solvophobic interfaces.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Asynchronous in situ connected-components analysis for complex fluid flows

Fluid flow in porous media is at the heart of problems such as groundwater contamination and carbon sequestration, and presents an important challenge for scientific computing. For this class of problem, large three-dimensional simulations are performed to advance scientific inquiry. On massively parallel computing systems, the volume of data generated by such approaches can become a productivity bottleneck if the raw data generated from the simulation is analyzed in a post-processing step. We present a physics-based framework for in situ data reduction that is theoretically grounded in multiscale averaging theory. We show how task parallelism can be exploited to concurrently perform a variety of analysis tasks with data-dependent costs, including the generation of iso-surfaces, morphological analyses, and connected components analysis. A task management framework is constructed to launch asynchronous analysis threads, manage dependencies between different tasks, promote data locality and hide the impact of data transfers. The framework is applied to analyze GPU-based simulations of two-fluid flow in porous media, generating a set of averaged measures that represents the overall system behavior. We demonstrate how the approach can be applied to perform physically-consistent averaging over fluid subregions using connected components analysis. Simulations performed on Oak Ridge National Lab s Titan supercomputer are profiled to demonstrate the performance of the associated multi-threaded in situ analysis for typical production simulation of two-fluid flow.

Mcclure, James↗
DOE OSTI2023

Relationships between fluid mixing, biodiversity, and chemosynthetic primary productivity in Yellowstone hot springs

The factors that influence biodiversity and productivity of hydrothermal ecosystems are not well understood. Here we investigate the relationship between fluid mixing, biodiversity, and chemosynthetic primary productivity in three co-localized hot springs (RSW, RSN, and RSE) in Yellowstone National Park that have different geochemistry. All three springs are sourced by reduced hydrothermal fluid, but RSE and RSN receive input of vapour phase gas and oxidized groundwaters, with input of both being substantially higher in RSN. Metagenomic sequencing revealed that communities in RSN were more biodiverse than those of RSE and RSW in all dimensions evaluated. Microcosm activity assays indicate that rates of dissolved inorganic carbon (DIC) uptake were also higher in RSN than in RSE and RSW. Together, these results suggest that increased mixing of reduced volcanic fluid with oxidized fluids generates additional niche space capable of supporting increasingly biodiverse communities that are more productive. These results provide insight into the factors that generate and maintain chemosynthetic biodiversity in hydrothermal systems and that influence the distribution, abundance, and diversity of microbial life in communities supported by chemosynthesis. These factors may also extend to other ecosystems not supported by photosynthesis, including the vast subterranean biosphere and biospheres beneath ice sheets and glaciers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Electrically controlled interlayer trion fluid in electron-hole bilayers

Here, the combination of repulsive and attractive Coulomb interactions in a quantum electron-hole (e-h) fluid can produce correlated phases of multiparticle charge complexes, such as excitons, trions, and biexcitons. We report an experimental realization of an electrically controlled interlayer trion fluid in van der Waals heterostructures. In strongly coupled e-h bilayers, electrons and holes spontaneously form three-particle trion bound states. The interlayer trions can assume 1e-2h and 2e-1h configurations. We show that the two holes in 1e-2h trions form a spin-singlet with a spin gap of approximately one milli–electron volt. By electrostatic gating, the equilibrium state can be continuously tuned into an exciton fluid, a trion fluid, an exciton-trion mixture, or a trion-charge mixture. Our work demonstrates a platform to study correlated phases of tunable Bose-Fermi mixtures.

Qi, Ruishi [University of California, Berkeley, CA↗
DOE OSTI2015

Stochastic Reductions for Inertial Fluid-Structure Interactions Subject to Thermal Fluctuations

We explore the dynamics of elastic microstructures within a fluid that are subjected to thermal fluctuations. We perform analysis to obtain systematically simplified descriptions of the mechanics in the limiting regimes when (i) the coupling forces that transfer momentum between the fluid and microstructures are strong, (ii) the mass of the microstructures is small relative to the displaced mass of the fluid, and (iii) the response to stresses results in hydrodynamics that relax rapidly to a quasi-steady-state relative to the motions of the microstructure. Furthermore, we derive effective equations using a singular perturbation analysis of the backward Kolmogorov equations of the stochastic process. Our continuum mechanics description is based on the stochastic Eulerian--Lagrangian method which provides a framework for approximation of the fluid-structure interactions when subject to thermal fluctuations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗
DOE OSTI2023

Hamiltonian structure of 2D fluid dynamics with broken parity

Isotropic fluids in two spatial dimensions can break parity symmetry and sustain transverse stresses which do not lead to dissipation. Corresponding transport coefficients include odd viscosity, odd torque, and odd pressure. We consider an isotropic Galilean invariant fluid dynamics in the adiabatic regime with momentum and particle density conservation. We find conditions on transport coefficients that correspond to dissipationless and separately to Hamiltonian fluid dynamics. The restriction on the transport coefficients will help identify what kind of hydrodynamics can be obtained by coarse-graining a microscopic Hamiltonian system. Interestingly, not all parity-breaking transport coefficients lead to energy conservation and, generally, the fluid dynamics is energy conserving but not Hamiltonian. We show how this dynamics can be realized by imposing a nonholonomic constraint on the Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

User's Guide of TOUGH2-EGS. A Coupled Geomechanical and Reactive Geochemical Simulator for Fluid and Heat Flow in Enhanced Geothermal Systems Version 1.0

TOUGH2-EGS is a numerical simulation program coupling geomechanics and chemical reactions for fluid and heat flows in porous media and fractured reservoirs of enhanced geothermal systems. The simulator includes the fully-coupled geomechanical (THM) module, the fully-coupled geochemical (THC) module, and the sequentially coupled reactive geochemistry (THMC) module. The fully-coupled flow-geomechanics model is developed from the linear elastic theory for the thermo-poro-elastic system and is formulated with the mean normal stress as well as pore pressure and temperature. The chemical reaction is sequentially coupled after solution of flow equations, which provides the flow velocity and phase saturation for the solute transport calculation at each time step. In addition, reservoir rock properties, such as porosity and permeability, are subjected to change due to rock deformation and chemical reactions. The relationships between rock properties and geomechanical and chemical effects from poro-elasticity theories and empirical correlations are incorporated into the simulator. This report provides the user with detailed information on both mathematical models and instructions for using TOUGH2-EGS for THM, THC or THMC simulations. The mathematical models include the fluid and heat flow equations, geomechanical equation, reactive geochemistry equations, and discretization methods. Although TOUGH2-EGS has the capability for simulating fluid and heat flows coupled with both geomechanical and chemical effects, it is up to the users to select the specific coupling process, such as THM, THC, or THMC in a simulation. There are several example problems illustrating the applications of this program. These example problems are described in details and their input data are presented. The results demonstrate that this program can be used for field-scale geothermal reservoir simulation with fluid and heat flow, geomechanical effect, and chemical reaction in porous and fractured media.

15 GEOTHERMAL ENERGY↗
DOE OSTI2013

User's guide of TOUGH2-EGS-MP: A Massively Parallel Simulator with Coupled Geomechanics for Fluid and Heat Flow in Enhanced Geothermal Systems VERSION 1.0

TOUGH2-EGS-MP is a parallel numerical simulation program coupling geomechanics with fluid and heat flow in fractured and porous media, and is applicable for simulation of enhanced geothermal systems (EGS). TOUGH2-EGS-MP is based on the TOUGH2-MP code, the massively parallel version of TOUGH2. In TOUGH2-EGS-MP, the fully-coupled flow-geomechanics model is developed from linear elastic theory for thermo-poro-elastic systems and is formulated in terms of mean normal stress as well as pore pressure and temperature. Reservoir rock properties such as porosity and permeability depend on rock deformation, and the relationships between these two, obtained from poro-elasticity theories and empirical correlations, are incorporated into the simulation. This report provides the user with detailed information on the TOUGH2-EGS-MP mathematical model and instructions for using it for Thermal-Hydrological-Mechanical (THM) simulations. The mathematical model includes the fluid and heat flow equations, geomechanical equation, and discretization of those equations. In addition, the parallel aspects of the code, such as domain partitioning and communication between processors, are also included. Although TOUGH2-EGS-MP has the capability for simulating fluid and heat flows coupled with geomechanical effects, it is up to the user to select the specific coupling process, such as THM or only TH, in a simulation. There are several example problems illustrating applications of this program. These example problems are described in detail and their input data are presented. Their results demonstrate that this program can be used for field-scale geothermal reservoir simulation in porous and fractured media with fluid and heat flow coupled with geomechanical effects.

15 GEOTHERMAL ENERGY↗