Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TAYLOR INSTABILITY”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Modeling of Turbulence Effect on Liquid Jet Atomization

Recent studies indicate that turbulence behaviors within a liquid jet have considerable effect on the atomization process. Such turbulent flow phenomena are encountered in most practical applications of common liquid spray devices. This research aims to model the effects of turbulence occurring inside a cylindrical liquid jet to its atomization process. The two widely used atomization models Kelvin-Helmholtz (KH) instability of Reitz and the Taylor analogy breakup (TAB) of O'Rourke and Amsden portraying primary liquid jet disintegration and secondary droplet breakup, respectively, are examined. Additional terms are formulated and appropriately implemented into these two models to account for the turbulence effect. Results for the flow conditions examined in this study indicate that the turbulence terms are significant in comparison with other terms in the models. In the primary breakup regime, the turbulent liquid jet tends to break up into large drops while its intact core is slightly shorter than those without turbulence. In contrast, the secondary droplet breakup with the inside liquid turbulence consideration produces smaller drops. Computational results indicate that the proposed models provide predictions that agree reasonably well with available measured data.

Trinh, H. P.↗
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↗
NASA NTRS1984

Evolution equations for Taylor vortices in the small-gap limit

The centrifugal instability of the viscous fluid flow between concentric circular cylinders in the small gap limit is considered. The amplitude of the Taylor vortex is allowed to depend on a slow time variable, a slow axial variable, and the polar angle. It is shown that the amplitude of the vortex cannot in general be described by a single amplitude equation. However, if the axial variations are periodic a single amplitude equation can be derived. In the absence of any slow axial variations it is shown that a Taylor vortex remains stable to wavy vortex perturbations. Furthermore, in this situation, stable nonaxisymmeric modes can occur but do not bifurcate from the Taylor vortex state. The stability of these modes is shown to be governed by a modified form of the Eckhaus criterion. Previously announced in STAR as N84-10501

Hall, P.↗
DOE OSTI2020

Supernova Hydrodynamics: A Lab-scale Study of the Blast-driven Instability Using High-speed Diagnostics

A novel experimental approach to study the blast-driven instability at a nondiffuse, gaseous interface with a density gradient is presented. Under Euler similarity, this approach enables study of dissipative-scale hydrodynamics relevant to many astrophysical and laboratory high energy density phenomena in a well-resolved manner. The instability is initiated by passing a Taylor–Sedov blast wave originating from a controlled detonation through a perturbed and stably stratified interface between two gases. The facility and driving blast wave are characterized to obtain repeatable conditions and capture large ensembles of time-resolved Mie scattering imaging that show consistent hydrodynamic development. We analyze the instability evolution between different gas pairs to demonstrate the wide range of development and turbulent behavior that may occur between different supernova layers. The mean evolution of the hydrodynamic instability is compared to a buoyancy–drag model that is frequently used to estimate perturbation growth in supernova mixing research. We propose a time delay to this model in order to reproduce the measured instability behavior and demonstrate model robustness in handling flows driven by a time-varying acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗
NASA NTRS1987

On the onset of three-dimensionality and time-dependence in the Goertler vortex problem

The instability of large amplitude Goertler vortices in a growing boundary layer is discussed in the fully nonlinear regime. It is shown that a three-dimensional breakdown to a flow with wavy vortex boundaries similar to that which occurs in the Taylor vortex problem takes place. However, the instability is confined to the thin shear layers which have been shown to trap the region of vortex activity. The disturbance eigenfunctions decay exponentially away from the center of these layers so that the upper and lower shear layers can support independent modes of instability. The structure of the instability, in particular its location and speed of downstream propagation, is found to be entirely consistent with recent experimental results. Furthermore, it is shown that the upper and lower layers support wavy vortex instabilities with quite different frequencies. This result is again consistent with the available experimental observations.

Hall, Philip↗
DOE OSTI2023

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Reduction of acoustic disturbances in the test section of supersonic wind tunnels by laminarizing their nozzle and test section wall boundary layers by means of suction

The feasibility of quiet, suction laminarized, high Reynolds number (Re) supersonic wind tunnel nozzles was studied. According to nozzle wall boundary layer development and stability studies, relatively weak area suction can prevent amplified nozzle wall TS (Tollmien-Schlichting) boundary layer oscillations. Stronger suction is needed in and shortly upstream of the supersonic concave curvature nozzle area to avoid transition due to amplified TG (Taylor-Goertler) vortices. To control TG instability, moderately rapid and slow expansion nozzles require smaller total suction rates than rapid expansion nozzles, at the cost of larger nozzle length Re and increased TS disturbances. Test section mean flow irregularities can be minimized with suction through longitudinal or highly swept slots (swept behind local Mach cone) as well as finely perforated surfaces. Longitudinal slot suction is optimized when the suction-induced crossflow velocity increases linearly with surface distance from the slot attachment line toward the slot (through suitable slot geometry). Suction in supersonic blowdown tunnels may be operated by one or several individual vacuum spheres.

Pfenninger, W.↗
NASA NTRS1989

On the onset of three-dimensionality and time-dependence in Goertler vortices

The secondary instability of large-amplitude Goertler vortices in a growing boundary layer is discussed in the fully nonlinear regime. It is shown that the three-dimensional breakdown to a flow with wavy vortex boundaries, similar to that which occurs in the Taylor vortex problem takes place. However, the instability is confined to the thin shear layers which were shown by Hall & Lakin (1988) to trap the region of vortex activity. The disturbance eigenfunctions decay exponentially away from the center of these layers, so that the upper and lower shear layers can support independent modes of instability. The structure of the instability, in particular its location and speed of downstream propagation, is found to be entirely consistent with recent experimental results. Furthermore, it is shown that the upper and lower layers support wavy vortex instabilities with quite different frequencies. This result is again consistent with the available experimental observations.

Hall, Philip↗
NASA NTRS2002

Bridgman Growth of Germanium Crystals in a Rotating Magnetic Field

A series of (100)-oriented gallium-doped germanium crystals have been grown by the Bridgman method and under the influence of a rotating magnetic field (RMF). The RMF has a marked affect on the interface shape, changing it from concave to nearly flat. The onset of time-dependent flow instabilities occurs when the critical magnetic Taylor number is exceeded, and this can be observed by noting the appearance of striations in the grown crystals. The critical magnetic Taylor number is a sensitive function of the aspect ratio and, as the crystal grows under a constant applied magnetic field, the induced striations change from nonperiodic to periodic, undergo a period-doubling transition, and then cease to exist. Also, by pulsing the RMF on and off, it is shown that intentional interface demarcations can be introduced.

Volz, M. P.↗
NASA NTRS2004

Bridgman Growth of Germanium Crystals in a Rotating Magnetic Field

A series of (100)-oriented gallium-doped germanium crystals have been grown by the Bridgman method and under the influence of a rotating magnetic field (RMF). Time-dependent flow instabilities occur when the critical magnetic Taylor number (Tm(sup c) is exceeded, and this can be observed by noting the appearance of striations in the grown crystals. The experimental data indicate that Tm(sup c) increases as the aspect ratio of the melt decreases. Modeling calculations predicting Tm(sup c) as a function of aspect ratio are in reasonable agreement with the experimental data. The RMF has a marked affect on the interface shape, changing it from concave to nearly flat as the RMF strength is increased. Also, by pulsing the RMF on and off, it is shown that intentional interface demarcations can be introduced.

Volz, M. P.↗
NASA NTRS2005

Bridgman Growth of Germanium Crystals in a Rotating Magnetic Field

A series of (100)-oriented gallium-doped germanium crystals has been grown by the vertical Bridgman method and under the influence of a rotating magnetic field (RMF). Time-dependent flow instabilities occur when the critical magnetic Taylor number (Tm(sup c)) is exceeded, and this can be observed by noting the appearance of striations in the grown crystals. Tm(sup c) decreases as the aspect ratio of the melt increases, and approaches the theoretical limit expected for an infinite cylinder. Intentional interface demarcations are introduced by pulsing the RMF on and off The RMF has a marked affect on the interface shape, changing it from concave to nearly flat as the RMF strength is increased.

Volz, M. P.↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “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,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.

42 ENGINEERING↗
NASA NTRS1991

Short-time Lyapunov exponent analysis and the transition to chaos in Taylor-Couette flow

The physical mechanism driving the weakly chaotic Taylor-Couette flow is investigated using the short-time Liapunov exponent analysis. In this procedure, the transition from quasi-periodicity to chaos is studied using direct numerical 3D simulations of axially periodic Taylor-Couette flow, and a partial Liapunov exponent spectrum for the flow is computed by simultaneously advancing the full solution and a set of perturbations. It is shown that the short-time Liapunov exponent analysis yields more information on the exponents and dimension than that obtained from the common Liapunov exponent calculations. Results show that the chaotic state studied here is caused by a Kelvin-Helmholtz-type instability of the outflow boundary jet of Taylor vortices.

Vastano, John A.↗
NASA NTRS1995

Interstellar clouds in high-speed, supersonic flows: Two-dimensional simulations

We present a series of gasdynamical simulations of the interaction of a dense, cool interstellar cloud with a high-speed, supersonic wind that confines and accelerates the embedded cloud. Our goal is to attempt to determine if such clouds can survive various potentially disruptive instabilities, that occur at their peripheries, long enough to be accelerated to speeds which are comparable to the wind velocity. These simulations are performed using two-dimensional, Eulerian gas dynamics on both an axisymmetric (about the cloud axis) and 'slab' geometric grid. The spatial and temporal resolutions of the simulations are varied over a wide range to investigate the effects of small-scale instabilities on the overall acceleration of clouds and the development of large-scale, disruptive instabilities. Also, we study the effects of wind/cloud Mach number variations by changing the wind speed constant at about 12 km/s (which corresponds to a cloud temperature of 10,000 K). The current simulations track the evolution of clouds as they are accelerated to speeds approximately 4-5 times greater than their internal sound speeds. Furthermore, the models with the highest resolution were extended far beyond quasi-linear Rayleigh-Taylor growth times reaching 6-7 Rayleigh-Taylor growth times for the largest scale instabilities before being terminated because of the accumulation of errors at the rear grid boundary.

Schiano, A. V. R.↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗