Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid instabilities”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Emergence of Unstable Focused Flow Induced by Variable‐Density Flows in Vertical Fractures

Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.

Environmental Sciences & Ecology↗
NASA NTRS1994

Effect of chemical heat release in a temporally evolving mixing layer

Two-dimensional numerical simulations of a temporally evolving mixing layer with an exothermic infinitely fast diffusion flame between two unmixed reactants have been carried out in the limit of zero Mach number to study the effect of the heat release on the early stages of the evolution of the flow. Attention has been directed to relatively large values of the oxidizer-to-fuel mass stoichiometric ratio typical of hydrocarbon flames, and initial vorticity distributions thicker than the temperature and species distributions have been chosen to mimic the situation at the outlet of a jet. The results show that, during the stages of the evolution covered by the present simulations, enhancement of combustion occurs by local stretching of the flame without much augmentation of its area. The rate of product generation depends strongly on the initial conditions, which suggests the possibility of controlling the combustion by acting on the flow. Rollup and vortex amalgamation still occur in these reacting flows but are very much affected by the production of new vorticity by baroclinic torques. These torques lead to counter rotating vortex pairs around the flame and, more importantly, in thin layers of light fluid that leave the vicinity of the flame when the Kelvin-Helmholtz instability begins to develop. Propelled by the vortex pairs, these layers wind around, split on reaching high pressure regions, and originate new vortex pairs in a process that ends up building large-scale vortices with a vorticity distribution more complex than for a constant density fluid.

Higuera, F. J.↗
DOE OSTI2020

Confinement Controls the Bend Instability of Three-Dimensional Active Liquid Crystals

Spontaneous growth of long-wavelength deformations is a defining feature of active liquid crystals. We investigate the effect of confinement on the instability of 3D active liquid crystals in the isotropic phase composed of extensile microtubule bundles and kinesin molecular motors. When shear aligned, such fluids exhibit finite-wavelength self-amplifying bend deformations. By systematically changing the channel size we elucidate how the instability wavelength and its growth rate depend on the channel dimensions. Experimental findings are qualitatively consistent with a minimal hydrodynamic model, where the fastest growing deformation is set by a balance of active driving and elastic relaxation. Furthermore, our results demonstrate that confinement determines the structure and dynamics of active fluids on all experimentally accessible length scales.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1988

Radiative instabilities in sheared magnetic field

The structure and growth rate of the radiative instability in a sheared magnetic field B have been calculated analytically using the Braginskii fluid equations. In a shear layer, temperature and density perturbations are linked by the propagation of sound waves parallel to the local magnetic field. As a consequence, density clumping or condensation plays an important role in driving the instability. Parallel thermal conduction localizes the mode to a narrow layer where K(parallel) is small and stabilizes short wavelengths k larger-than(c) where k(c) depends on the local radiation and conduction rates. Thermal coupling to ions also limits the width of the unstable spectrum. It is shown that a broad spectrum of modes is typically unstable in tokamak edge plasmas and it is argued that this instability is sufficiently robust to drive the large-amplitude density fluctuations often measured there.

Drake, J. F.↗
NASA NTRS1985

Theoretical and experimental design studies for the Atmospheric General Circulation Experiment

The major criterion for the Atmospheric General Circulation Experiment (AGCE) design is that it be possible to realize strong baroclinic instability in the spherical configuration chosen. A configuration was selected in which a hemispherical shell of fluid is subjected to latitudinal temperature gradients on its spherical boundaries and the latitudinal boundaries are insulators. Work in the laboratory with a cylindrical version of this configuration revealed more instabilities than baroclinic instability. Since researchers fully expect these additional instabilities to appear in the spherical configuration also, they decided to continue the laboratory cylindrical annulus studies. Four flow regimes were identified: an axisymmetric Hadley circulation, boundary layer convection, baroclinic waves and deep thermal convection. Regime diagrams were prepared.

Fowlis, W. W.↗
NASA NTRS1994

Local shear instabilities in weakly ionized, weakly magnetized disks

We extend the analysis of axisymmetric magnetic shear instabilities from ideal magnetohydrodynamic (MHD) flows to weakly ionized plasmas with coupling between ions and neutrals caused by collisions, ionization, and recombination. As part of the analysis, we derive the single-fluid MHD dispersion relation without invoking the Boussinesq approximation. This work expands the range of applications of these instabilities from fully ionized accretion disks to molecular disks in galaxies and, with somewhat more uncertainty, to protostellar disks. Instability generally requires the angular velocity to decrease outward, the magnetic field strengths to be subthermal, and the ions and neutrals to be sufficiently well coupled. If ionization and recombination processes can be neglected on an orbital timescale, adequate coupling is achieved when the collision frequency of a given neutral with the ions exceeds the local epicyclic freqency. When ionization equilibrium is maintained on an orbital timescale, a new feature is present in the disk dynamics: in contrast to a single-fluid system, subthermal azimuthal fields can affect the axisymmetric stability of weakly ionized two-fluid systems. We discuss the underlying causes for this behavior. Azimuthal fields tend to be stabilizing under these circumstances, and good coupling between the neutrals and ions requires the collision frequency to exceed the epicyclic frequency by a potentially large secant factor related to the magnetic field geometry. When the instability is present, subthermal azimuthal fields may also reduce the growth rate unless the collision frequency is high, but this is important only if the field strengths are very subthermal and/or the azimuthal field is the dominant field component. We briefly discuss our results in the context of the Galactic center circumnuclear disk, and suggest that the shear instability might be present there, and be responsible for the observed turbulent motions.

Blaes, Omer M.↗
DOE OSTI2021

Measurement of Deformation Heterogeneity During Shale Swelling Using Digital Image Correlation

Abstract Rock-fluid interactions in shale formations are one of the main sources of wellbore instability issues and inadequate stimulation performance. For better planning of fracturing operations and optimizing production from these reservoirs, it is important to understand the mechanisms behind these interactions. These issues are especially prevalent in clay-rich shales, such as Tuscaloosa Marine Shale, which is the subject of this study. Conventional techniques to quantify these shale-fluid interactions comprise of measuring swelling in powdered rock grains or measurement of deformation in the whole pieces of the core using a linear variable differential transformer and strain gages. However, the contribution from individual laminae to overall deformation cannot be evaluated using these methods. In this study, we developed an experimental setup to evaluate the spatial deformation in shale during interaction with water using digital image correlation (DIC). Deformation of two shale samples, with 34 wt% to 51 wt% clay content, was studied. White paint was used to generate a random speckle pattern on the specimen and then immersed in deionized water. The deformation process was captured using a digital camera and images were analyzed using DIC to quantify the deformation. The implementation of the DIC technique enables the visualization and quantification of spatial deformation in the specimen during interacting with water. The results show the localization of large strains in select laminations. The results provide a better understanding of shale deformation when interacting with water in comparison to traditional measurements that can provide only an average strain value.

Energy & Fuels↗
NASA NTRS1998

Finite Element Modeling of Magnetically-Damped Convection during Solidification

A fully 3-D, transient finite element model is developed to represent the magnetic damping effects on complex fluid flow, heat transfer and electromagnetic field distributions in a Sn- 35.5%Pb melt undergoing unidirectional solidification. The model is developed based on our in- house finite element code for the fluid flow, heat transfer and electromagnetic field calculations. The numerical model is tested against numerical and experimental results for water as reported in literature. Various numerical simulations are carried out for the melt convection and temperature distribution with and without the presence of a transverse magnetic field. Numerical results show that magnetic damping can be effectively applied to stabilize melt flow, reduce turbulence and flow levels in the melt and over a certain threshold value a higher magnetic field resulted in a greater reduction in velocity. Also, for the study of melt flow instability, a long enough running time is needed to ensure the final fluid flow recirculation pattern. Moreover, numerical results suggest that there seems to exist a threshold value of applied magnetic field, above which magnetic damping becomes possible and below which the 0 convection in the melt is actually enhanced.

deGroh, H. C.↗
DOE OSTI2026

Linear simulations of a wide pedestal quiescent H-mode plasma with the extended-MHD code NIMROD

We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

LES of a Jet Excited by the Localized Arc Filament Plasma Actuators

The fluid dynamics of a high-speed jet are governed by the instability waves that form in the free-shear boundary layer of the jet. Jet excitation manipulates the growth and saturation of particular instability waves to control the unsteady flow structures that characterize the energy cascade in the jet.The results may include jet noise mitigation or a reduction in the infrared signature of the jet. The Localized Arc Filament Plasma Actuators (LAFPA) have demonstrated the ability to excite a high-speed jets in laboratory experiments. Extending and optimizing this excitation technology, however, is a complex process that will require many tests and trials. Computational simulations can play an important role in understanding and optimizing this actuator technology for real-world applications. Previous research has focused on developing a suitable actuator model and coupling it with the appropriate computational fluid dynamics (CFD) methods using two-dimensional spatial flow approximations. This work is now extended to three-dimensions (3-D) in space. The actuator model is adapted to a series of discrete actuators and a 3-D LES simulation of an excited jet is run. The results are used to study the fluid dynamics near the actuator and in the jet plume.

Brown, Clifford A.↗
NASA NTRS2016

First-Principles Modeling of IAT-Driven Anomalous Resistivity in Hollow Cathode Discharges I: Theory

Self-consistent expressions are derived for the impact of ion acoustic turbulence (IAT) on the fluid properties of electrons and ions in the plume of a hollow cathode. This instability is believed to have a dominant role in the plasmadynamics in this region. The random phase approximation is applied to the perturbed Boltzmann equation to determine the effective drag and heating of the plasma species induced by the IAT. These driving terms are then incorporated into momentum and energy fluid equations for the ions and electrons. A fluid-like relation governing the evolution of the total IAT wave energy density is also presented. The physical significance of the IAT-based terms is discussed in the context of energy and momentum conservation, and a path forward is outlined for implementing the effects of IAT in the fluid-based model in Part II.

Jorns, Benjamin A.↗
NASA NTRS1973

Ion streaming instabilities with application to collisionless shock wave structure

The electromagnetic dispersion relation for two counterstreaming ion beams of arbitrary relative strength flowing parallel to a dc magnetic field is derived. The beams flow through a stationary electron background and the dispersion relation in the fluid approximation is unaffected by the electron thermal pressure. Magnetic effects on the ion beams are included, but the electrons are treated as a magnetized fluid. The dispersion relation is solved with a zero net current condition applied and the regions of instability in the k-U space (U is the relative velocity between the two ion beams) are presented. These results are extensions of Kovner's analysis for weak beams. The parameters are then chosen to be applicable for parallel shocks. It is found that unstable waves with zero group velocity in the shock frame can exist near the leading edge of the shock for upstream Alfven Mach numbers greater than 5.5.

Golden, K. I.↗
DOE OSTI2020

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗
DOE OSTI2025

A fluctuating hydrodynamics model for nanoscale surfactant-laden interfaces

A multispecies diffuse interface model is formulated in a fluctuating hydrodynamics framework for the purpose of simulating surfactant interfaces at the nanoscale. The model generalizes previous work to ternary mixtures, employing a Cahn-Hilliard free energy density combined with incompressible, isothermal fluctuating hydrodynamics where dissipative fluxes include both deterministic and stochastic terms. The intermolecular parameters in the free energy are chosen such that one species acts as a partially miscible surfactant. From Laplace pressure measurements, we show that in this model the surface tension decreases linearly with surfactant concentration, leading to Marangoni convection for interfaces with concentration gradients. In the capillary wave spectrum for interfaces with and without surfactant, we find that for the former, the spectrum deviates significantly from classical capillary wave theory, presumably due to Gibbs elasticity. In non-equilibrium simulations of the Rayleigh-Plateau instability, deterministic simulations showed that the surfactant delays pinching of a fluid cylinder into droplets. However, stochastic simulations indicate that thermal fluctuations disrupt the surfactant's stabilizing effect. Similarly, the spreading of a patch of surfactant, driven by Marangoni convection, was found to be partially suppressed by thermal fluctuations.

Capillary waves↗
NASA NTRS1973

The reconnection rate of magnetic fields.

Rapid reconnection of magnetic lines of force occurs whenever two regions of opposite field are pressed together without sufficient fluid pressure between to keep them apart. It has been suggested that rapid reconnection is the basis for solar flares, and that in the magnetopause and the neutral sheet of the geomagnetic tail it plays an essential role in geomagnetic activity. It is pointed out that the escape of fluid from between the fields is the result of magnetohydrodynamic exchange instabilities, leading to reconnection rates which are generally some fraction of the Alfven speed.

Parker, E. N.↗
NASA NTRS1975

Multifluids description of dynamics of upper atmosphere

A multifluids model to investigate ionospheric dynamics was established on kinetic theory. Its resultant equations are used to examine the following dynamic problems in the gamma region of 80-2000 Km of the ionosphere: (1) propagation of acoustic modes in the 500-2,000 Km of the ionosphere (two fluid model); (2) the relation between the cross field plasma drift instabilities and type I and type II ionospheric irregularities; and (3) time dependent neutral wind structure and horizontal pressure gradient.

Wu, S. T.↗