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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Effects of eletron heating on the current driven electrostatic ion cyclotron instability and plasma transport processes along auroral field lines

Fluid simulations of the plasma along auroral field lines in the return current region have been performed. It is shown that the onset of electrostatic ion cyclotron (EIC) related anomalous resistivity and the consequent heating of electrons leads to a transverse ion temperature that is much higher than that produced by the current driven EIC instability (CDICI) alone. Two processes are presented for the enhancement of ion heating by anomalous resistivity. The anomalous resistivity associated with the turbulence is limited by electron heating, so that CDICI saturates at transverse temperature that is substantially higher than in the absence of resistivity. It is suggested that this process demonstrates a positive feedback loop in the interaction between CDICI, anomalous resistivity, and parallel large-scale dynamics in the topside ionosphere.

Ganguli, Supriya B.↗

DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

NASA NTRS1985

Magnetic buoyancy instabilities incorporating rotation

Plane layer stability analysis is extended to incorporate the effects of uniform rotation. Detailed studies are made of the interchange, or 'axisymmetric modes' and of the undular or wavelike motions, considering both high and low frequency modes. The force due to rotation on the fluid motions are determined. The types of instability relevant to astrophysical problems, are also discussed with attention given to the occurrence of two distinct modes of instability in a bottom heavy field gradient. Some peculiar stability boundaries, which are due to the absence of a dominating force, are described.

Hughes, D. W.↗

DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗

DOE OSTI2024

Pattern formation in odd viscoelastic fluids

Nonreciprocal interactions fueled by local energy consumption can be found in biological and synthetic active matter at scales where viscoelastic forces are important. Such systems can be described by “odd” viscoelasticity, which assumes fewer material symmetries than traditional theories. Here we study odd viscoelasticity analytically and using lattice Boltzmann simulations. We identify a pattern-forming instability which produces an oscillating array of fluid vortices, and we elucidate which features govern the growth rate, wavelength, and saturation of the vortices. Our observation of pattern formation through odd mechanical response can inform models of biological patterning and guide engineering of odd dynamics in soft active matter systems. Published by the American Physical Society 2024

Floyd, Carlos (ORCID:0000000262707250)↗

DOE OSTI2021

Flow and heat transfer in the wake of a triangular arrangement of spheres

This research work seeks to investigate the influence of spacing and heat transfer on the wake behavior of a triangular arrangement of spheres. Four experimental configurations have been investigated at three Reynolds numbers, Re 1 = 350, Re 2 = 700, and Re 3 = 1050. Two isothermal cases were investigated with spacing between the spheres of zero and one sphere diameter, and two cases were investigated with an applied heat flux at the same spacing conditions. The time resolved particle image velocimetry results revealed various flow phenomena including flow separations, von Karman vortex shedding, and Kelvin–Helmholtz instabilities. The turbulent statistics reveal the effect of proximity and heat transfer on the time averaged values of the wake size, turbulent strengths, and Reynolds shear stress in the wake of each sphere, namely, the laminarization effects from the addition of heat and the suppression of the lead sphere wake from the proximity of the trailing spheres. These results are complemented by the application of proper orthogonal decomposition (POD) to the flow fields, which extracts the coherent structures from the flow. Here the modes that describe the coherent structures are extracted and described in detail, which provide further insight into effects of the experimental conditions on the temporal behavior of the flow. Many of the low order modes are found to be associated in pairs, corresponding to asymmetric structures or advection of a given structure downstream. The capability of POD to produce reduced order models of the flow is then utilized to facilitate vortex identification analysis. A turbulent kinetic energy based mode truncation criteria, which has been found to enhance vortex identification capability, is applied to select the POD modes and temporal coefficients to be used in the reduced order modeling. The reconstructed velocity fields are then analyzed with vortex identification algorithms to extract the vortex cores and boundaries. The combination of these approaches allows the study of the effect of proximity and heat transfer on the vortex characteristics, such as size, strength, and distribution.

42 ENGINEERING↗

DOE OSTI2022

Nonlinear harmonics coupled by parallel wave propagations in a time-dependent plasma flow

In a time-dependent flow, nonlinear harmonics can be excited by coupling between linear waves and flow-induced harmonic waves. Examining the dispersion relations and selection rules for the coupling, we investigate nonlinearly coupled harmonics for waves propagating along the magnetic field line in a magnetized plasma, as well as waves in an unmagnetized plasma. The coupled harmonics in a plasma flow are described by analytic dispersion relations and selection rules. This nonlinear coupling is corroborated by the particle-in-cell simulation. Here, the coupled-harmonics model describes a mechanism for the excitation of nonlinear harmonics from linear waves in a time-dependent flow. The spectral analysis of the dispersion relation provides a useful way to evaluate the spatiotemporal behavior of a plasma flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Observations of heteroclinic bifurcations in resistive magnetohydrodynamic simulations of the plasma response to resonant magnetic perturbations

A class of topological magnetic island bifurcations that has not previously been observed in toroidal plasmas is described. Increasing an externally applied three-dimensional magnetic field in resistive magnetohydrodynamic simulations results in the asymmetric elongation of resonant island flux surfaces followed by a sequence of heteroclinic bifurcations. Here, these bifurcations produce new sets of hyperbolic-elliptic fixed points as predicted by the Poincaré-Birkoff fixed point theorem. Field line calculations verify that the new fixed points do not connect to those of the prebifurcated islands as required for heteroclinic bifurcations on a torus with winding numbers composed of common integer factors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Exploration of Rayleigh-Taylor seeding mechanisms in laser- and pulsed-power-driven inertial fusion

The ubiquitous Rayleigh-Taylor instability (RTI) has been observed in fluids and plasmas in a wide range of parameter regimes from astrophysical to laboratory plasmas. In high-energy-density (HED) laboratory plasmas, such as laser-driven inertial confinement fusion (ICF) plasmas, the RTI can have a detrimental impact on achieving fusion ignition through the generation of hydrodynamic mix. Hence, an understanding of the seeding mechanisms that produce the RTI and identifying ways to mitigate hydrodynamic mix is of critical importance. Appropriately aligned magnetic fields have been known to stabilize short-wavelength RTI. A number of numerical and experimental studies have demonstrated the benefits of using imposed magnetic fields in laser-driven ICF to achieve higher ion temeratures, higher neutron yields, and a relaxation of the ignition criteria. This work addresses seeding mechanisms in laser-driven implosions to understand critical early-stage physics that ultimately leads to substantial growth of the RTI along with mechanisms for mitigation of this growth. Surface perturbations due to machining tolerances and single-feature seeds (for example due to fill tubes) can produce substantial RTI growth in the ignition-relevant high-convergence ratio targets for laser- driven ICF implosions. RTI growth from single-feature seeding in laser-based implosions has the potential to be mitigated through appropriately aligned externally applied magnetic fields and the goal of this work has been to quantify that numerically and experimentally. Resistive-magnetohydrodynamic (MHD) simu- lations are used to study the seeding and evolution of the RTI leveraging previous support from the DOE HEDLP program. Furthermore, experimental data has been obtained and applied towards code validation of unmagnetized and magnetized evolution of single-feature seeded RTI growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Simulations and models for the Richtmyer– Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P∝ k m , where k is the wavenumber. The hydrodynamic condition is the same as in the $\theta$-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The $\theta$-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ~1.86 shock and gamma-law = 5/3 fluids with Atwood number A = 0.5. The bubble amplitude h B is found to grow in two stages. Initially, h B undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, h B grows as a power law $~t^{\theta}$ when $\langle{k}\rangle{|h_B|}$ > O(1), where $\langle{k}\rangle$ is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1990

Interfacial dynamics of two liquids under an oscillating gravitational field

The evolution of two miscible liquids meeting at an initially sharp interface inside a cavity under microgravity g-jitter conditions is studied numerically. The response of the interface and kinematics of the flowfield to various g-jitter accelerations and aspect ratio variations is shown. The interface region acts like a vortex source sheet, and it can be unstable to Kelvin-Helmholtz and Rayleigh-Taylor instabilities. The vortices produced along the interface can serve as a stirring mechanism to promote local mixing. Below the critical Stokes-Reynolds number, the destabilization of the interface results in deformation into wavy structures. In some parameter regions, these structures oscillate in time; in others, they are quasi-steady. Above the critical Stokes-Reynolds number, 'chaotic' instability results, and the interface breaks into concentration pockets. The morphology of the initial breakup is similar to observed Rayleigh-Taylor instability. Subsequent mixing of the two fluids after the breakup of the interface is then very rapid.

Duval, Walter M. B.↗

NASA NTRS1991

Nonaxisymmetric instabilities in rotating shear flows - Internal gravity modes in stratified media and analogies with plane flows

The role of the internal gravity modes in mediating the growth of nonaxisymmetric instabilities is investigated by studying the instability of stratified incompressible differentially rotating fluid cylinders to global nonaxisymmetric modes. The results indicate that, in addition to a modified version of the well-known principal branch mediated by surface modes of the system (analogous to f-modes of stars), there exist unstable branches of the dispersion relation mediated by internal gravity modes of the system (similar to the g-modes of stars). These branches arise due to the interaction between the g-modes. It is shown that the maximum growth rate on one of the new branches can sometimes equal or exceed that on the principal branch, thus modifying the principal branch.

Ghosh, P.↗

NASA NTRS1992

Anomalous ion mixing within a Kelvin-Helmholtz vortex in a collisionless plasma

Anomalously fast ion mixing is observed in a hybrid code simulation (particle ions and a massless electron fluid) of the Kelvin-Helmholtz instability in a collisionless plasma. While the traditional view predicts that the ion mixing occurs on a time scale longer than the roll-up time T(v) of the Kelvin-Helmholtz vortex by a factor of (velocity shear scale length)/(ion Larmor radius), the simulation results show that the mixing occurs within a time scale comparable to T(v). It is shown that this fast mixing is due to the scattering of ions by unsteady electromagnetic fields.

Terasawa, T.↗

NASA NTRS2004

Torque Balances on the Taylor Cylinders in the Geomagnetic Data Assimilation

In this presentation we report on our continuing effort in geomagnetic data assimilation, aiming at understanding and predicting geomagnetic secular variation on decadal time scales. In particular, we focus on the effect of the torque balances on the cylindrical surfaces in the core co-axial with the Earth's rotation axis (the Taylor cylinders) on the time evolution of assimilated solutions. We use our MoSST core dynamics,model and observed geomagnetic field at the Earth's surface derived via Comprehensive Field Model (CFM) for the geomagnetic data assimilation. In our earlier studies, a model solution is selected randomly from our numerical database. It is then assimilated with the observations such that the poloidal field possesses the same field tomography on the core-mantel boundary (CMB) continued downward from surface observations. This tomography change is assumed to be effective through out the outer core. While this approach allows rapid convergence between model solutions and the observations, it also generates sevee numerical instabilities: the delicate balance between weak fluid inertia and the magnetic torques on the Taylor cylinders are completely altered. Consequently, the assimilated solution diverges quickly (in approximately 10% of the magnetic free-decay time in the core). To improve the assimilation, we propose a partial penetration of the assimilation from the CMB: The full-scale modification at the CMB decreases linearly and vanish at an interior radius r(sub a). We shall examine from our assimilation tests possible relationships between the convergence rate of the model solutions to observations and the cut-off radius r(sub a). A better assimilation shall serve our nudging tests in near future.

Kuang, Weijia↗

NASA NTRS2000

Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Skarda, J. Raymond Lee↗

NASA NTRS2002

Phase Transition in Dusty Plasmas: A Microphysical Description

Dust grains immersed in plasma discharges acquire a large negative charge and settle into a dust cloud at the edge of the sheath. In this region, the plasma ions stream toward the electrode at a velocity u approx. cs=(T(sub e)/m(sub i))(exp 1/2). Experimentally at sufficiently high gas pressure P, the random kinetic energy of the grains is damped by gas friction, and the grains are strongly coupled and self-organize into a crystalline configuration. For lower pressures despite the dissipation of grain kinetic energy to gas friction, the dust grains reach a steady-state kinetic temperature T(sub d) which is much larger than the temperature of any other component in the plasma. T(sub d) is so large that the dust acts like a fluid. We have used the dynamically shielded dust (DSD) model to simulate these physical processes. We find that the known experimental features are nicely reproduced in the simulations, and that additional features are revealed. In the figure we plot the variation of T(sub d) as P is continuously varied in a DSD code run. A marked difference is evident between the critical pressure P(sub m) for the melting transition as P is decreased, and the critical pressure P(sub c) for the condensation transition as P is increased. For P(sub m) is less than P is less than P(sub c), mixed phase states are seen. This hysteresis occurs because the instability which triggers melting is different from the instability that heats the dust in the fluid phase and inhibits freezing. At low pressure, the dust is subject to a two-stream instability with the ions. This instability is responsible for the high temperature of the dust at low pressure. The basic physics underlying the melting transition has been elucidated in a series of papers. We are developing a first-principles analytic approach to the melting transition, which embodies the same physics that is present in the DSD code.

Joyce, Glenn↗

NASA NTRS2019

Operational Stability Limits in Rotating Detonation Engine Numerical Simulations

An instability is described which arises in computational fluid dynamic (CFD) simulations of semi-idealized rotating detonation engines (RDE) configured with a throat at the exit. Its existence is verified by examining output from two independently developed CFD codes simulating the same configuration and producing solutions that agree well. The instability is shown to be thermo-acoustic in that a spatial integral of the product of pressure and heat release fluctuations develops a regular oscillation which grows in time. The instability can become severe enough to cause detonation failure. Its onset is shown to be closely linked to the size of the exit throat and the size of the inlet restriction; both parameters that strongly influence RDE performance. It is shown that the instability places a cap on ideal RDEperformance, but that an optimized exhaust throat and inlet restriction combination still yields substantial pressure gain. Other parametric sensitivities are also examined in terms ofinstability growth. These include axial length, inlet manifold pressure, and air-fuel ratio.

Paxson, Daniel E.↗

NASA NTRS1996

Electric Field Induced Interfacial Instabilities

The study of the interface in a charge-free, nonpolar, critical and near-critical binary fluid in the presence of an externally applied electric field is presented. At sufficiently large fields, the interface between the two phases of the binary fluid should become unstable and exhibit an undulation with a predefined wavelength on the order of the capillary length. As the critical point is approached, this wavelength is reduced, potentially approaching length-scales such as the correlation length or critical nucleation radius. At this point the critical properties of the system may be affected. In zero gravity, the interface is unstable at all long wavelengths in the presence of a field applied across it. It is conjectured that this will cause the binary fluid to break up into domains small enough to be outside the instability condition. The resulting pattern formation, and the effects on the critical properties as the domains approach the correlation length are of acute interest. With direct observation, laser light scattering, and interferometry, the phenomena can be probed to gain further understanding of interfacial instabilities and the pattern formation which results, and dimensional crossover in critical systems as the critical fluctuations in a particular direction are suppressed by external forces.

Kusner, Robert E.↗