Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow 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 1,531 records · Page 85

Hypersonic Boundary Layer Instability Over a Corner

A boundary-layer transition study over a compression corner was conducted under a hypersonic flow condition. Due to the discontinuities in boundary layer flow, the full Navier-Stokes equations were solved to simulate the development of disturbance in the boundary layer. A linear stability analysis and PSE method were used to get the initial disturbance for parallel and non-parallel flow respectively. A 2-D code was developed to solve the full Navier-stokes by using WENO(weighted essentially non-oscillating) scheme. The given numerical results show the evolution of the linear disturbance for the most amplified disturbance in supersonic and hypersonic flow over a compression ramp. The nonlinear computations also determined the minimal amplitudes necessary to cause transition at a designed location.

Balakumar, Ponnampalam↗

Computational Study of Laminar Flow Control on a Subsonic Swept Wing Using Discrete Roughness Elements

A combination of parabolized stability equations and secondary instability theory has been applied to a low-speed swept airfoil model with a chord Reynolds number of 7.15 million, with the goals of (i) evaluating this methodology in the context of transition prediction for a known configuration for which roughness based crossflow transition control has been demonstrated under flight conditions and (ii) of analyzing the mechanism of transition delay via the introduction of discrete roughness elements (DRE). Roughness based transition control involves controlled seeding of suitable, subdominant crossflow modes, so as to weaken the growth of naturally occurring, linearly more unstable crossflow modes. Therefore, a synthesis of receptivity, linear and nonlinear growth of stationary crossflow disturbances, and the ensuing development of high frequency secondary instabilities is desirable to understand the experimentally observed transition behavior. With further validation, such higher fidelity prediction methodology could be utilized to assess the potential for crossflow transition control at even higher Reynolds numbers, where experimental data is currently unavailable.

Li, Fei↗

Numerical experiments in boundary-layer stability

Numerical solution of the three-dimensional incompressible Navier-Stokes equations is used to study the instability of a flat-plate boundary layer in a manner analogous to the vibrating-ribbon experiments. Flow-field structures are observed which are very similar to those found in the vibrating-ribbon experiment to which computational initial conditions have been matched. Streamwise periodicity is assumed in the simulation so that the evolution occurs in time, but the events which constitute the instability are so similar to the spatially occurring ones of the laboratory that it seems clear the physical processes involved are the same. A spectral and finite difference numerical algorithm is employed in the simulation.

Wray, A.↗

Effect of 3D Roughness Patch on Instability Amplification in a Supersonic Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of high-speed vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, sinusoidal shape roughness elements on a Mach 3.5 flat plate boundary layer for flow conditions corresponding to the planned conditions of an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at the NASA Langley Research Center. Analysis of convective instabilities in the wake of the roughness patch was reported in a previous paper and the current work extends that analysis to instability amplification across the length of the roughness patch. Quasi-parallel stability analysis of the modified boundary layer flow over the patch indicates two dominant families of unstable disturbances, namely, a group of high frequency modes that amplify in localized regions along the roughness patch and another group of low frequency modes that have smaller peak amplification rates but amplify steadily both above the roughness patch and in the wake region behind it. The results suggest that the amplification factors associated with the high-frequency modes are sufficiently low, at least for the roughness patches considered in this paper, so that these modes are unlikely to have a major influence on the transition process. The amplification of the low-frequency modes within the region of the roughness patch is further quantified via direct numerical simulations. Results confirm the strongly destabilizing influence of the roughness patch on the first mode instabilities, yielding an N-factor increment of ΔN ≈ 3.6 for a roughness patch length of eight wavelengths.

Laminar turbulent transition↗

On the instability of boundary layers on heated flat plates

The stability of a boundary layer on a heated flat plate is investigated in the linear regime. The flow is shown to be unstable to longitudinal vortex structures which in general develop in a nonparallel manner in the streamwise direction. Solutions of the nonparallel equations are obtained numerically at O(1) values of the appropriate stability parameter, i.e., the Grashof number. The particular cases investigated relate to the situations when the instability is induced by localized or distributed wall roughness or nonuniform wall heating. The case when the vortices are induced by freestream disturbances is also considered. The fastest growing mode is found to be governed by a quasi-parallel theory at high wavenumbers. The wavenumber and growth rate of the fastest growing mode are found in closed form. At low wavenumbers the vortex instability is shown to be closely related to Tollmein-Schlichting waves. The effect of wall heating or cooling on the latter type of instability is discussed.

Hall, Philip↗

Boundary-layer stability analysis of Langley Research Center 8-foot LFC experimental data

An analytical study of linear-amplifying instabilities of a laminar boundary layer as found in the experimental data of the LaRC/8-foot laminar-flow control (LFC) experiment was completed and the results are presented. The LFC airfoil used for this experiment was a swept, supercritical design which removed suction air through spanwise slots. The amplification of small disturbances by linear processes on a swept surface such as this can be due to either Tollmien-Schlichting (TS) and/or crossflow (CF) mechanisms. This study consists of the examination of these two instabilities by both the commonly used incompressible (SALLY and MARIA) analysis and the more involved compressible (COSAL) analysis. A wide range of experimental test conditions with variations in Mach number, Reynolds number, and suction distributions were available for this study. Experimentally determined transition locations were found from thin-film techniques and were used to correlate the n-factors at transition for the range of test cases.

Berry, Scott↗

Pick-Up Ion Instabilities at Planetary Magnetospheres

This effort involved the analysis of low frequency waves as observed by the Galileo spacecraft near the Galilean moon, Io. Io is a significant source of material, especially SO2, and various products of dissociation, and further these atoms and molecules are readily ionized. The initial velocity of the ions is essentially that of the neutral species, i.e., the Keplerian velocity. The plasma, on the other hand is co-rotating, and there is a differential flow of the order 57 km/s between the plasma and the neutral particles. Thus pick-up ion instabilities are Rely to occur within the Jovian magnetosphere. Indeed, magnetometer observations from the Galileo spacecraft clearly show ion cyclotron waves that have been identified with a large variety of plasma species, such as O+, S++ (which has the same gyro-frequency as O+), S+, and SO2+. Typically, however, the dominant frequency is near the SO2+ gyro-frequency. The research effort was originally planned to be a team effort between Robert J. Strangeway as the Principal Investigator, and Debbie Huddleston, who was an Assistant Research Geophysicist at UCLA. Unfortunately, Dr. Huddleston took a position within Industry. The effort was therefore descoped, and Dr. Strangeway instead pursued a collaboration with Dr. Xochitl Blanco-Cano, of the Instituto de Geofisica, Universidad Nacional Autonoma de Mexico. This has proved to be a productive collaboration, with several papers and publications arising out of the effort. The magnetic field oscillations near lo generally fall into two types: ion cyclotron waves, with frequencies near an ion gyro-frequency, and lower frequency mirror-mode waves. The ion cyclotron waves are mainly transverse, and frequently propagate along the ambient magnetic field. The mirror-mode waves are compressional waves, and they have essentially zero frequency in the plasma rest frame. One of the purposes of our investigation is to understand what controls the types of wave modes that occur, since both wave modes can be drive unstable by the pressure anisotropy associated with the pick up ions. The pick ion velocity is perpendicular to the ambient magnetic field, and is generally much larger than the thermal velocity, at least initially. At its simplest, we found that the ion cyclotron waves are controlled by the parameters of the species in gyro-resonance with the wave. Thus, while the growth rates for the lower mass (higher gyro-frequency) pick-up ions are generally larger, we found that the heavier SO2+ ion cyclotron waves are generally preferred. This is because one of the effects of the wave instability is to diffuse the ions in pitch angle and energy. The lower mass ions therefore consist of both a ring of recently created pick-up ions, and a thermal background. This thermal background quenches the ion cyclotron instability. SO2+ is different, however. Being a molecule, the species can also dissociate. Our analysis suggests that the dissociation acts on a time scale comparable to or faster than the velocity space diffusion time scale. There are consequently no thermal SO2+ ions to quench the instability. We have also investigated the mirror-mode. This mode can at times grow more rapidly than the individual ion cyclotron waves. This is mainly because the mirror-mode can grow off the pressure anisotropy of the individual species with which the waves are in resonance. Lastly, as part of this effort we have begun to investigate the instability for obliquely propagating modes. Galileo observations show that at times the ion cyclotron waves are significantly elliptically polarized, and further the wave vector is at a large angle to the field, significant growth can occur for oblique propagation. Depending on the group velocity of the waves, it is possible that obliquely propagating modes have higher advective growth rate, but this has yet to be determined.

Strangeway, Robert J.↗

On flow phenomena that emit X-rays in hot star winds

The X-ray emission from O stars may originate from instability-generated shocks in their stellar winds. Previous numerical simulations that assumed the wind to be isothermal could only draw limited conclusions concerning this emission. Calculations that include the energy transfer in the wind are presented. They confirm that up to a few stellar radii, radiative cooling is efficient which implies that the shock cooling zones are short in comparison with the dynamical lengths. At larger radii, the wind structure changes as all the shocks are destroyed due to a broadening of their cooling zones. The following flow phenomena are discussed as possible X-ray sources: inner radiative shock; shock merging, outer adiabatic shocks, remaining hot gas from shock destruction, and an outer corona.

Feldmeier, A.↗

Supersonic jet noise generated by large scale instabilities

The role of large scale wavelike structures as the major mechanism for supersonic jet noise emission is examined. With the use of aerodynamic and acoustic data for low Reynolds number, supersonic jets at and below 70 thousand comparisons are made with flow fluctuation and acoustic measurements in high Reynolds number, supersonic jets. These comparisons show that a similar physical mechanism governs the generation of sound emitted in he principal noise direction. These experimental data are further compared with a linear instability theory whose prediction for the axial location of peak wave amplitude agrees satisfactorily with measured phased averaged flow fluctuation data in the low Reynolds number jets. The agreement between theory and experiment in the high Reynolds number flow differs as to the axial location for peak flow fluctuations and predicts an apparent origin for sound emission far upstream of the measured acoustic data.

Seiner, J. M.↗

Convection and segregation in directional solidification of dilute and non-dilute binary alloys - Effects of ampoule and furnace design

A Petrov-Galerkin/finite-element method is used to analyze the effect of furnace configuration and ampoule design on the temperature field, the convection in the melt, the shape of the melt-solid interface, and the segregation of solute in the crystal, in the directional solidification of several dilute and nondilute binary semiconductor alloys. The vertical Bridgman-Stockbarger system leads to a two-cell flow structure, and the constant gradient furnace has only a single cell near the interface. Radial temperature gradients interact with the solute field to cause the previously predicted sideways diffusive instability, and addition of a solute that is less dense than the bulk melt and that is incorporated upon solidification decreases the intensity of the flow near the interface by increasing the melt density there. The present results have application to understanding of the effects on solute segregation of microgavity solidification and of applied magnetic fields.

Adornato, Peter M.↗

Does the ballooning instability trigger substorms in the near-Earth magnetotail?

The stability of the near-Earth magnetotail against ballooning (or configurational) instability is examined in the framework of the MHD approximation. It is emphasized that a change in plasma pressure induced by a meriodional electric field drift delta u(sub n) is an important factor that determines the stability. We have to consider two ways in which plasma pressure changes, that is, a convective change -delta u(sub n) grad(P(sub 0)), where P(sub 0) is background plasma pressure, and plasma expansion/compression -P(sub 0) dive (delta u(sub n)). Since delta u(sub n) is perpendicular to the magnetic field and its magnitude is inversely proportional to the magnetic field strength, delta u(sub n) diverges/converges in usual tail magnetic field configurations. For the instability, the convective change must overwhelm the effects of the plasma expansion/compression. However, near the equator in the near-Earth tail, the latter may overcompensate for the former. We describe the ballooning instability in terms of a coupling between the Alfven and slow magnetosonic waves in an inhomogeneous plasma and derive instability conditions. The result shows that the excessive curvature stabilizes, rather than destabilizes, perturbations. It is also found that the field-aligned flow stabilizes perturbations, as well as the field-aligned current. We infer that under quiet conditions, the plasma pressure gradient in the near-Earth tail is not sharp enough to trigger the instability. The plasma sheet is expected to become more stable during the substorm growth phase because of an increase in the field line curvature associated with the plasma sheet thinning. In the region closer to the Earth, including the ring current, the plasma pressure gradient may be localized in a limited range of the radial distance during the growth phase. However, recently reported plasma and magnetic field parameters before substorm onsets do not provide very convincing evidence that the ballooning instability is the triggering mechanism of substorms.

Ohtani, Shin-Ichi↗

Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

The regulation of electron heat transport in high-β, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as β$^{-1}_{e}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of β$_e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected β$^{-1}_{e}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90° scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of Nose Bluntness on Hypersonic Boundary-Layer Receptivity and Stability Over Cones

The receptivity to freestream acoustic disturbances and the stability properties of hypersonic boundary layers are numerically investigated for boundary-layer flows over a 5 straight cone at a freestream Mach number of 6.0. To compute the shock and the interaction of the shock with the instability waves, the Navier-Stokes equations in axisymmetric coordinates were solved. In the governing equations, inviscid and viscous flux vectors are discretized using a fifth-order accurate weighted-essentially-non-oscillatory scheme. A third-order accurate total-variation-diminishing Runge-Kutta scheme is employed for time integration. After the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. The appearance of instability waves near the nose region and the receptivity of the boundary layer with respect to slow mode acoustic waves are investigated. Computations confirm the stabilizing effect of nose bluntness and the role of the entropy layer in the delay of boundary-layer transition. The current solutions, compared with experimental observations and other computational results, exhibit good agreement.

Kara, Kursat↗

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Functional Renormalization Group Study of Superconductivity in Rhombohedral Trilayer Graphene

We employ a functional renormalization group approach to ascertain the pairing mechanism and symmetry of the superconducting phase observed in rhombohedral trilayer graphene. Superconductivity in this system occurs in a regime of carrier density and displacement field with a weakly distorted annular Fermi sea. We find that repulsive Coulomb interactions can induce electron pairing on the Fermi surface by taking advantage of momentum-space structure associated with the finite width of the Fermi sea annulus. Furthermore, the degeneracy between spin-singlet and spin-triplet pairing is lifted by valley-exchange interactions that strengthen under the RG flow and develop nontrivial momentum-space structure. We find that the leading pairing instability is d-wave-like and spin-singlet, and that the theoretical phase diagram versus carrier density and displacement field agrees qualitatively with experiment.

36 MATERIALS SCIENCE↗