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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

On the Goertler instability in hypersonic flows: Sutherland law fluids and real gas effects

The Goertler vortex instability mechanism in a hypersonic boundary layer on a curved wall is investigated. The precise roles of the effects of boundary layer growth, wall cooling, and gas dissociation is clarified in the determination of stability properties. It is first assumed that the fluid is an ideal gas with viscosity given by Sutherland's law. It is shown that when the free stream Mach number M is large, the boundary layer divides into two sublayers: a wall layer of O(M sup 3/2) thickness over which the basic state temperature is O(M squared) and a temperature adjustment layer of O(1) thickness over which the basic state temperature decreases monotonically to its free stream value. Goertler vortices which have wavelengths comparable with the boundary layer thickness are referred to as wall modes. It is shown that their downstream evolution is governed by a set of parabolic partial differential equations and that they have the usual features of Goertler vortices in incompressible boundary layers. As the local wavenumber increases, the neutral Goertler number decreases and the center of vortex activity moves towards the temperature adjustment layer. Goertler vortices with wavenumbers of order one or larger must necessarily be trapped in the temperature adjustment layer and it is this mode which is most dangerous. For this mode, it was found that the leading order term in the Goertler number expansion is independent of the wavenumber and is due to the curvature of the basic state. This term is also the asymptotic limit of the neutral Goertler numbers of the wall mode. To determine the higher order corrections terms in the Goertler number expansion, two wall curvature cases are distinguished. Real gas effects were investigated by assuming that the fluid is an ideal dissociating gas. It was found that both gas dissociation and wall cooling are destabilizing for the mode trapped in the temperature adjustment layer, but for the wall mode trapped near the wall the effect of gas dissociation can be either destabilizing or stabilizing.

Fu, Yibin B.↗
DOE OSTI2022

A mechanism for reduced compression in indirectly driven layered capsule implosions

High-yield implosions on the National Ignition Facility rely on maintaining low entropy in the deuterium–tritium fuel, quantified by its adiabat, in order to efficiently couple energy to the hot spot through high compression of the fuel layer. We present very-high-resolution xRAGE simulation results that study the impacts of interfacial mixing and the jetting of materials due to surface defects, defects on internal interfaces, voids, and engineering features on fuel layer compression. Defects and voids are typically neglected in implosion simulations due to their small size and three-dimensional geometry. Our results showed that supersonic jets of material arise through weak spots in the shell at peak implosion velocity that prevent uniform compression of the fuel layer even when they do not introduce contaminant into the hot spot. This occurs despite maintaining low fuel entropy, since the formation of the weak spots involves nonradial displacement of fuel mass. In contrast, simulations show that fuel–ablator mixing due to interfacial instabilities has a much smaller impact on compression. We show that defects on interior interfaces of plastic capsules decrease compression by 15% to 25% and interfacial mixing between the ablator and fuel decreases compression by less than 1% for implosions with plastic or high-density carbon (HDC) ablators. For low adiabat implosions, the impact of jetting seeded by the support tent can also decrease the compression by 25%. We demonstrate that the inclusion of interior defects in simulations can explain the inferred compression in two fielded plastic capsule implosions and that the inclusion of voids, for which available characterization has large uncertainties, in simulations of HDC capsule implosions has a qualitatively consistent impact. This mechanism offers a potential explanation for persistently overestimated fuel compression in design simulations of layered implosions on the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

New Type of the Interface Evolution in the Richtmyer-Meshkov Instability

We performed systematic theoretical and numerical studies of the nonlinear large-scale coherent dynamics in the Richtmyer-Meshkov instability for fluids with contrast densities. Our simulations modeled the interface dynamics for compressible and viscous uids. For a two-fluid system we observed that in the nonlinear regime of the instability the bubble velocity decays and its surface attens, and the attening is accompanied by slight oscillations. We found the theoretical solution for the system of conservation laws, describing the principal influence of the density ratio on the motion of the nonlinear bubble. The solution has no adjustable parameters, and shows that the attening of the bubble front is a distinct property universal for all values of the density ratio. This property follows from the fact that the RM bubbles decelerate. The theoretical and numerical results validate each other, describe the new type of the bubble front evolution in RMI, and identify the bubble curvature as important and sensitive diagnostic parameter.

Abarzhi, S. I.↗
DOE OSTI2022

Faster ablative Kelvin–Helmholtz instability growth in a magnetic field

Shear flows along a plasma interface will quickly grow unstable due to the Kelvin–Helmholtz instability. If there is a concurrent temperature gradient across the interface, higher modes are stabilized by the thermal diffusion. These ablative effects must be considered in, for example, jet features in inertial confinement fusion hot-spots, or plasma plumes in young supernovae. We show that magnetization of the plasma can greatly affect the instability, even if magnetic pressure is small. This is because electrons are localized by their gyromotion, reducing the heat flux and material ablation. We use a two-dimensional numerical extended-magnetohydrodynamics approach to assess this effect for dense fusion conditions. In comparison with the unmagnetized case, self-generated Biermann fields make only a minor difference to growth rates. However, simulations with a large 50 kT external field found that the growth rate of the least stable mode increased by 40%. In conclusion, this has implications for mix processes in Z-pinches and magnetized inertial confinement fusion concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗
DOE OSTI2022

Experiments on the single-mode Richtmyer–Meshkov instability with reshock at high energy densities

Here, the hydrodynamic instability growth of a reshocked single-mode interface between high energy density fluids is studied. A laser-driven shock wave is used to drive an initially solid, sinusoidal interface between a dense plastic (1.43 g/cc) and a light foam (≈ 0.110 g/cc). After the interface has grown to a nonlinear state where the amplitude is of order of the wavelength, it is reshocked. The reshock compresses the nonlinear perturbation, which then grows at about twice the rate. While the pre-reshock growth rate is sensitive to the initial amplitude and wavelength of the perturbation, the post-reshock growth rate is comparatively insensitive to the initial condition. Qualitatively, we observe that the perturbations are less coherent after reshock, consistent with the idea that having a reshock accelerates the transition to turbulence. We find that some memory of the initial condition remains, even after reshock at late time: it appears if the initial perturbations have large enough wavelengths, and the flow structure of size comparable to the initial wavelength persists through reshock. Our results agree with design simulations and are consistent with the phenomenology of reshock studies in conventional gaseous shock tubes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2005

Problems in Microgravity Fluid Mechanics: G-Jitter Convection

This is the final report on our NASA grant, Problems in Microgravity Fluid Mechanics NAG3-2513: 12/14/2000 - 11/30/2003, extended through 11/30/2004. This grant was made to Stanford University and then transferred to the University of California at Santa Barbara when the PI relocated there in January 2001. Our main activity has been to conduct both experimental and theoretical studies of instabilities in fluids that are relevant to the microgravity environment, i.e. those that do not involve the action of buoyancy due to a steady gravitational field. Full details of the work accomplished under this grant are given below. Our work has focused on: (i) Theoretical and computational studies of the effect of g-jitter on instabilities of convective states where the convection is driven by forces other than buoyancy (ii) Experimental studies of instabilities during displacements of miscible fluid pairs in tubes, with a focus on the degree to which these mimic those found in immiscible fluids. (iii) Theoretical and experimental studies of the effect of time dependent electrohydrodynamic forces on chaotic advection in drops immersed in a second dielectric liquid. Our objectives are to acquire insight and understanding into microgravity fluid mechanics problems that bear on either fundamental issues or applications in fluid physics. We are interested in the response of fluids to either a fluctuating acceleration environment or to forces other than gravity that cause fluid mixing and convection. We have been active in several general areas.

Homsy, G. M.↗
NASA NTRS1990

Stability of an oscillated fluid with a uniform density gradient

Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Jacqmin, David↗
DOE OSTI2022

Shock-driven dispersal of a corrugated finite-thickness particle layer

A research area emerging in the multiphase flow community is the study of shock-driven multiphase instability (SDMI), a gas–particle analog of the traditional fluid-fluid Richtmyer–Meshkov instability (RMI). In this work, we study the interaction of planar air shocks with corrugated glass particle curtains through the use of numerical simulations with an Eulerian–Lagrangian approach. One objective of this study is to compare the simulated particle curtains to a comparable set of shock tube experiments performed to analyze traditional RMI of a gas curtain. The simulations are set to match the experimental shock Mach numbers and perturbation wavelengths (3.6 and 7.2 mm) while also matching the Atwood number of the experiments to the multiphase Atwood number of the simulations. Varying particle diameters are tested in the simulations to explore the impact of particle diameter on the evolution of the particle curtain. This simulation setup allows for a one-to-one comparison between RMI and SDMI under comparable conditions while also allowing for a separate study into the validity of the use of the multiphase Atwood number to compare the single-phase and multiphase instabilities. In particular, we show that the comparison depends on the diameter of the particles (thus, dependent on the Stokes number of the flow). A second objective of this study is to analyze the effect of the initial particle volume fraction on the evolution of the curtain and the behavior of the instability. This is done through analyzing the effect of the multiphase terms of the vorticity evolution equation on the vorticity deposition in SDMI. Also discussed is the effect of the particle diameter on the multiphase generation terms as well as in the baroclinic vorticity generation term in SDMI as the shock passes over the curtain.

42 ENGINEERING↗
NASA NTRS1990

Buoyancy-driven instability in a vertical cylinder: Binary fluids with Soret effect. I - General theory and stationary stability results

The buoyancy-driven instability of a monocomponent or binary fluid completely contained in a vertical circular cylinder is investigated, including the influence of the Soret effect for the binary mixture. The Boussinesq approximation is used, and the resulting linear stability problem is solved using a Galerkin technique. The analysis considers fluid mixtures ranging from gases to liquid metals. The flow structure is found to depend strongly on both the cylinder aspect ratio and the magnitude of the Soret effect. The predicted stability limits are shown to agree closely with experimental observations.

Hardin, G. R.↗
NASA NTRS1980

Extracting energy from natural flow

Three concepts for extracting energy from wind, waterflow, and tides utilize flow instability to generate usable energy. Proposed converters respond to vortex excitation motion, galloping or plunging motion, and flutter. Fluid-flow instability is more efficient in developing lift than is direct flow.

Delionback, L. M.↗
DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gravothermal phase transition, black holes and space dimensionality

In the framework of gravothermal evolution of an ideal monatomic fluid, I examine the dynamical instability of the fluid sphere in (N+1) dimensions by exploiting Chandrasekhar’s criterion to each quasi-static equilibrium along the sequence of the evolution. Once the instability is triggered, it would probably collapse into a black hole if no other interaction halts the process. From this viewpoint, the privilege of (3+1)-dimensional spacetime is manifest, as it is the marginal dimensionality in which the ideal monatomic fluid is stable but not too stable. Moreover, it is the unique dimensionality that allows stable hydrostatic equilibrium with positive cosmological constant. While all higher dimensional (N > 3) spheres are genuinely unstable. In contrast, in (2+1)-dimensional spacetime it is too stable either in the context of Newton’s theory of Gravity or Einstein’s General Relativity. It is well-known that the role of negative cosmological constant is crucial to have the Banados-Teitelboim-Zanelli (BTZ) black hole solution and the equilibrium configurations of a fluid disk. Owing to the negativeness of the cosmological constant, there is no unstable configuration for a homogeneous fluid disk to collapse into a naked singularity, which supports the Cosmic Censorship Conjecture. However, BTZ holes of mass M BTZ > 0 could emerge from collapsing fluid disks. The implications of spacetime dimensionality are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1996

Current Results and Proposed Activities in Microgravity Fluid Dynamics

The Institute for Problems in Mechanics' Laboratory work in mathematical and physical modelling of fluid mechanics develops models, methods, and software for analysis of fluid flow, instability analysis, direct numerical modelling and semi-empirical models of turbulence, as well as experimental research and verification of these models and their applications in technological fluid dynamics, microgravity fluid mechanics, geophysics, and a number of engineering problems. This paper presents an overview of the results in microgravity fluid dynamics research during the last two years. Nonlinear problems of weakly compressible and compressible fluid flows are discussed.

Polezhaev, V. I.↗
DOE OSTI2023

Reduced mixing in inertial confinement fusion with early-time interface acceleration

In inertial confinement fusion (ICF) implosions, the interface between the cryogenic DT fuel and the ablator is unstable to shock acceleration (the Richtmyer-Meshkov instability, RM) and constant acceleration (Rayleigh-Taylor instability, RT). Instability growth at this interface can reduce the final compression, limiting fusion burnup. If the constant acceleration is in the direction of the lighter material (negative Atwood number), the RT instability produces oscillatory motion that can stabilize against RM growth. Theory and simulations suggest this scenario occurred at early times in some ICF experiments on the National Ignition Facility, possibly explaining their favorable performance compared to one-dimensional simulations. In conclusion, this characteristic is being included in newer, lower adiabat designs, seeking to improve compression while minimizing ablator mixing into the fuel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of outer-cylinder rotation on the radially heated Taylor–Couette flow

A Taylor–Couette setup with radial heating is considered where a Boussinesq fluid is sheared in the annular region between two concentric, independently rotating cylinders maintained at different temperatures. Linear stability analysis is performed to determine the Taylor number for the onset of instability. Two radius ratios corresponding to wide and thin gaps with several rotation rate ratios are considered. The rotation of the outer cylinder is found to have a general stabilizing effect on the stability threshold as compared to pure inner-cylinder rotation, with a few exceptions. The radial heating sets up an axial flow which breaks the reflection symmetry of isothermal Taylor–Couette flow in the axial coordinate. This symmetry breaking separates linear stability thresholds, and we find the fastest growing modes with both positive and negative azimuthal numbers for different parameters. Another important finding of the current study is the discovery of unstable modes in the Rayleigh-stable regime. Furthermore, closed disconnected neutral curves (CDNCs) are observed for both wide and thin gaps which can separate from or merge into open neutral–stability curves. Alternatively, CDNCs can also morph into open neutral stability curves as the rotation rate ratio is changed. CDNCs are observed to be sensitive to changes in control parameters, and their appearance/disappearance is shown to induce discontinuous jumps in the critical Taylor number. Finally, for both wide and thin gaps, the fastest-growing modes found in the pure corotation case are shown to have their origins in the instability islands at smaller values of rotation rate ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Vortex-sheet modeling of hydrodynamic instabilities produced by an oblique shock interacting with a perturbed interface in the HED regime

In this work, we consider hydrodynamic instabilities produced by the interaction of an oblique shock with a perturbed material interface under high-energy-density (HED) conditions. During this interaction, a baroclinic torque is generated along the interface due to the misalignment between the density and pressure gradients, thus leading to perturbation growth. Our objective is to understand the competition between the impulsive acceleration due to the normal component of the shock velocity, which drives the Richtmyer–Meshkov instability, and the shear flow across the interface due to the tangential component of the shock velocity, which drives the Kelvin–Helmholtz instability, as well as its relation to perturbation growth. Since the vorticity resulting from the shock-interface interaction is confined to the interface, we describe the perturbation growth using a two-dimensional vortex-sheet model. We demonstrate the ability of the vortex-sheet model to reproduce roll-up dynamics for non-zero Atwood numbers by comparing to past laser-driven HED experiments. We determine the dependence of the interface dynamics on the tilt angle and propose a time scaling for the perturbation growth at early time. Eventually, this scaling will serve as a platform for the design of future experiments. This study is the first attempt to incorporate into a vortex-sheet model the time-dependent interface decompression and the deceleration (as well as the corresponding Rayleigh–Taylor instability) arising from laser turn-off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗