Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow 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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

The HEAO-2 Guest Investigator Program: Non-linear growth of instabilities in line-driven stellar winds

The linear instability of line-driven stellar winds to take proper account of the dynamical effect of scattered radiation were analyzed. It is found that: (1) the drag effect of the mean scattered radiation does greatly reduce the contribution of scattering lines to the instability at the very base of the wind, but the instability growth rate associated with such lines rapidly increases as the flow moves outward from the base, reaching more than 50% of the growth rate for pure absorption lines within a stellar radius of the surface, and eventually reaching 80% of that rate at large radii; (2) perturbations in the scattered radiation field may be important for the propagation of wind disturbances, but they have little effect on the wind instability; and (3) the contribution of strongly shadowed lines to the wind instability is often reduced compared to that of unshadowed lines, but their overall effect is not one of damping in the outer parts of the wind. It is concluded that, even when all scattering effects are taken into account, the bulk of the flow in a line-driven stellar wind is still highly unstable.

Rybicki, G. B.↗
NASA NTRS1996

The Extensional Rheology of Non-Newtonian Materials

The evolution of the transient extensional stresses in dilute and semi-dilute viscoelastic polymer solutions are measured with a filament stretching rheometer of a design similar to that first introduced by Sridhar, et al. The solutions are polystyrene-based (PS) Boger fluids that are stretched at constant strain rates ranging from 0.6 less than or equal to epsilon(0) less than or equal to 4s(exp -1) and to Hencky strains of epsilon greater than 4. The test fluids all strain harden and Trouton ratios exceeding 1000 are obtained at high strains. The experimental data strain hardens at lower strain levels than predicted by bead-spring FENE models. In addition to measuring the transient tensile stress growth, we also monitor the decay of the tensile viscoelastic stress difference in the fluid column following cessation of uniaxial elongation as a function of the total imposed Hencky strain and the strain rate. The extensional stresses initially decay very rapidly upon cessation of uniaxial elongation followed by a slower viscoelastic relaxation, and deviate significantly from FENE relaxation predictions. The relaxation at long times t is greater than or equal to 5 s, is compromised by gravitational draining leading to non-uniform filament profiles. For the most elastic fluids, partial decohension of the fluid filament from the endplates of the rheometer is observed in tests conducted at high strain rates. This elastic instability is initiated near the rigid endplate fixtures of the device and it results in the progressive breakup of the fluid column into individual threads or 'fibrils' with a regular azimuthal spacing. These fibrils elongate and bifurcate as the fluid sample is elongated further. Flow visualization experiments using a modified stretching device show that the instability develops as a consequence of an axisymmetry-breaking meniscus instability in the nonhomogeneous region of highly deformed fluid near the rigid endplate.

Spiegelberg, Stephen H.↗
NASA NTRS2010

Global Instability on Laminar Separation Bubbles-Revisited

In the last 3 years, global linear instability of LSB has been revisited, using state-of-the-art hardware and algorithms. Eigenspectra of LSB flows have been understood and classified in branches of known and newly-discovered eigenmodes. Major achievements: World-largest numerical solutions of global eigenvalue problems are routinely performed. Key aerodynamic phenomena have been explained via critical point theory, applied to our global mode results. Theoretical foundation for control of LSB flows has been laid. Global mode of LSB at the origin of observable phenomena. U-separation on semi-infinite plate. Stall cells on (stalled) airfoil. Receptivity/Sensitivity/AFC feasible (practical?) via: Adjoint EVP solution. Direct/adjoint coupling (the Crete connection). Minor effect of compressibility on global instability in the subsonic compressible regime. Global instability analysis of LSB in realistic supersonic flows apparently quite some way down the horizon.

Theofilis, Vassilis↗
NASA NTRS2019

Effects of Forward-Facing Step Shape on Stationary Crossflow Instability Growth and Breakdown

The effect of forward-facing step shape on stationary-crossflow growth was studied downstream of a two-dimensional step in a swept-wing flow. The step height was set at 81% of the local unperturbed boundary layer thickness at the step location, and the shape was varied by placing two ramps of differing slopes upstream of the step: 30◦ and 45◦. The ramps were successful at delaying transition, with increasing success as the slope was decreased. The ramps significantly altered the mean flow, reducing or eliminating the separation region and crossflow-reversal region downstream of the step, which reduced the amplification of the stationary crossflow modes shortly downstream of the step. The mean flow modification induced by the step also had destabilizing effects that persisted farther downstream as the flow recovered from the near-step region. The downstream growth of the stationary crossflow instabilities was also mitigated by the introduction of the ramps upstream of the step. The presence of harmonics beginning near this region confirms previous findings that there are nonlinear interactions contributing to the growth of the instability.

Eppink, Jenna L.↗
NASA NTRS1993

Role of secondary instability theory and parabolized stability equations in transition modeling

In modeling the laminar-turbulent transition region, the designer depends largely on benchmark data from experiments and/or direct numerical simulations that are usually extremely expensive. An understanding of the evolution of the Reynolds stresses, turbulent kinetic energy, and quantifies in the transport equations like the dissipation and production is essential in the modeling process. The secondary instability theory and the parabolized stability equations method are used to calculate these quantities, which are then compared with corresponding quantities calculated from available direct numerical simulation data for the incompressible boundary-layer flow of laminar-turbulent transition conditions. The potential of the secondary instability theory and the parabolized stability equations approach in predicting these quantities is discussed; results indicate that inexpensive data that are useful for transition modeling in the early stages of the transition region can be provided by these tools.

El-Hady, Nabil M.↗
NASA NTRS2003

An Accretion Model for the Growth of Black Hole in Quasars

A possible accretion model associated with the ionization instability of quasar disks is proposed to address the growth of the central black hole harbored in the host galaxy. The evolution of quasars in cosmic time is assumed to change from a highly active state to a quiescent state triggered by the S-shaped ionization instability of the quasar accretion disk. For a given external mass transfer rate ionization instability can modify accretion rate in the disk and separates the accretion flows of the disk into three different phases like a S-shape. We suggest that the bright quasars observed today are those quasars with disks in the upper branch of S-shaped instability and the dormant quasars are the system in the lower branch. The disk is assumed to evolve as ADIOS configuration in the lower branch. The mass ratio between black hole and its host galactic bulge is a nature consequence of ADIOS. Our model also demonstrates that a seed black hole 2 x 10(exp 6) solar masses similar to those found in spiral galaxies today is needed to produce a black hole with a final mass 2 x 10(exp 8) solar masses.

Lu, Ye↗
NASA NTRS2003

The Growth of Central Black Hole and the Ionization Instability of Quasar Disk

A possible accretion model associated with the ionization instability of quasar disks is proposed to address the growth of the central black hole harbored in the host galaxy. The evolution of quasars in cosmic time is assumed to change from a highly active state to a quiescent state triggered by the S-shaped ionization instability of the quasar accretion disk. For a given external mass transfer rate supplied by the quasar host galaxy, ionization instability can modify accretion rate in the disk and separates the accretion flows of the disk into three different phases, like a S-shape. We suggest that the bright quasars observed today are those quasars with disks in the upper branch of S-shaped instability, and the faint or 'dormant' quasars are simply the system in the lower branch. The middle branch is the transition state which is unstable. We assume the quasar disk evolves according to the advection-dominated inflow-outflow solutions (ADIOS) configuration in the stable lower branch of S-shaped instability, and Eddington accretion rate is used to constrain the accretion rate in each phase. The mass ratio between black hole and its host galactic bulge is a nature consequence of ADIOS. Our model also demonstrates that a seed black hole (BH) similar to those found in spiral galaxies today is needed to produce a BH with a final mass 2 x 10(exp 8) solar mases.

Lu, Ye↗
NASA NTRSDate not supplied

On A Stabilization of the Ingard-Myers Impedance Boundary Condition

It has been well-known that the Ingard-Myers impedance condition, while simple to apply, is subject to the hydrodynamic Kelvin-Helmholtz instability due to its use of a vortex sheet in modeling the flow at the liner boundary. Recently, in the development of a time domain boundary element method for acoustic scattering by treated surfaces, it was found that by neglecting a certain second-order spatial derivative term in the Ingard-Myers formulation, the hydrodynamic instability can be avoided. The present paper aims to provide further analysis of this modified condition, hereby referred to as Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC). It will be shown, based on the dispersion relations of linear waves, that the instability intrinsic to the Ingard-Myers condition is eliminated in the proposed new formulation. Quantitative assessments on the accuracy of TIMIBC for scattering of acoustic waves by lined surfaces will be carried out, and its effectiveness will be demonstrated by numerical examples. Specifically, the accuracy is assessed by comparing solutions obtained by the Ingard-Myers condition with that by the proposed TIMIBC where theoretical reflection coefficients at a lined surface are derived for cases of plane and spherical incident waves. It is found that the TIMIBC provides a good approximation to the original Ingard-Myers condition for flows of low to mid subsonic Mach numbers. Time domain implementations of TIMIBC are also discussed and illustrated with a numerical example using a finite difference scheme. As many studies have shown the Ingard-Myers condition to be the correct limit of boundary layer thickness going to zero, the proposed TIMIBC can offer a practical solution for overcoming the intrinsic instability associated with the Ingard-Myers condition.

Time Domain BEM↗
NASA NTRSDate not supplied

On A Stabilization of the Ingard-Myers Impedance Boundary Condition

It has been well-known that the Ingard-Myers impedance condition, while simple to apply, is subject to the hydrodynamic Kelvin-Helmholtz instability due to its use of a vortex sheet in modeling the flow at the liner boundary. Recently, in the development of a time domain boundary element method for acoustic scattering by treated surfaces, it was found that by neglecting a certain second-order spatial derivative term in the Ingard-Myers formulation, the hydrodynamic instability can be avoided. The present paper aims to provide further analysis of this modified condition, hereby referred to as Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC). It will be shown, based on the dispersion relations of linear waves, that the instability intrinsic to the Ingard-Myers condition is eliminated in the proposed new formulation. Quantitative assessments on the accuracy of TIMIBC for scattering of acoustic waves by lined surfaces will be carried out, and its effectiveness will be demonstrated by numerical examples. Specifically, the accuracy is assessed by comparing solutions obtained by the Ingard-Myers condition with that by the proposed TIMIBC where theoretical reflection coefficients at a lined surface are derived for cases of plane and spherical incident waves. It is found that the TIMIBC provides a good approximation to the original Ingard-Myers condition for flows of low to mid subsonic Mach numbers. Time domain implementations of TIMIBC are also discussed and illustrated with a numerical example using a finite difference scheme. As many studies have shown the Ingard-Myers condition to be the correct limit of boundary layer thickness going to zero, the proposed TIMIBC can offer a practical solution for overcoming the intrinsic instability associated with the Ingard-Myers condition.

Time Domain BEM↗
NASA NTRS1985

Instabilities in line-driven stellar winds. II - Effect of scattering

An earlier analysis (Owocki and Rybicki) of the linear instability of line-driven stellar winds is extended to take proper account of the dynamical effect of scattered radiation. The principal findings are as follows: (1) the drag effect of the mean scattered radiation does indeed greatly reduce the contribution of scattering lines to the instability at the very base of the wind, but the instability growth rate associated with such lines rapidly increases as the flow moves outward from the base, reaching more than 50 percent of the growth rate for pure absorption lines within a stellar radius of the surface, and eventually reaching 80 percent of that rate at large radii; (2) perturbations in the scattered radiation field may be important for the propagation of wind disturbances, but they have little effect on the wind instability; (3) the contribution of a strongly shadowed line to the wind instability is often reduced compared to that of an unshadowed line, but its effect is not one of damping in the outer parts of the wind. The primary conclusion derived from these results is thus that, even when all scattering effects are taken into account, the bulk of the flow in a line-driven stellar wind is still highly unstable.

Owocki, S. P.↗
NASA NTRS1990

The role of magnetic fields in cluster cooling flows

An investigation is made of the dynamical effects of the intracluster magnetic field, whose radial inflow and shear can produce a dramatic increase in the field's strength while rendering it more radial, with cooling flows. It is found that field reconnection is the most likely dominant-loss mechanism, so that buoyancy effects are probably not important. Attention is given to the effect of the magnetic field on thermal instabilities. The most important observable effect of the magnetic field in cooling flows will probably be very strong Faraday rotation of the polarization of radio sources within or behind the cooling flow.

Soker, Noam↗
NASA NTRS1982

Secondary effects in combustion instabilities leading to flashback

The secondary effects in turbulent combustion instabilities leading to flashback are investigated, including those due to buoyancy and contraction at the combustor outlet. Experiments were conducted in an oblong, rectangular cross-section combustion tunnel, where the effects of a bluff-body flames holder were generated by a rear-facing step behind a streamlined inlet nozzle. The results of experiments leading to flashback with the step mounted at the bottom of the combustion chamber were compared to those of experiments in which it was located at the top. Irrespective of the flow obstructions introduced downstream, the critical equivalence ratio for flashback was consistently lower with the step at the bottom, indicating that buoyancy was enhancing the growth of the recirculation zone that pushed the flame upstream and caused flashback. The contraction at the end of the combustion chamber had a promoting influence on the process of vortex pairing, re-enforcing the influence of the trailing vortices over that of the recirculation vortex system, and thereby curbing the tendency to flashback. Provided that the flow velocity was low, however, the characteristic features of combustion instabilities leading to flashback in the absence of contraction could still be established in its presence.

Vaneveld, L.↗
NASA NTRS1990

Effects of shock on the stability of hypersonic boundary layers

A set of linearized shock boundary conditions is derived, which is then imposed at the shock to account for the interaction of the shock wave with the boundary/shock layer instability wave; these boundary conditions are used to study the effect of shock on hypersonic boundary layer stability under the assumption of quasi-parallel flow. The result show that the shock has little effect on the boundary layer instability (subsonic first and second mode disturbances) when the shock is located outside the boundary layer edge. When the shock is located near the boundary layer edge, it exerts a stabilizing influence on the first and second modes. The shock also induces unstable supersonic modes with oscillatory structure in the shock layer, but these modes grow slower than the subsonic modes.

Chang, Chau-Lyan↗
NASA NTRS1992

Weakly nonlinear models for turbulent mixing in a plane mixing layer

New closure models for turbulent free shear flows are presented in this paper. They are based on a weakly nonlinear theory with a description of the dominant large-scale structures as instability waves. Two models are presented that describe the evolution of the free shear flows in terms of the time-averaged mean flow and the dominant large-scale turbulent structure. The local characteristics of the large-scale motions are described using linear theory. Their amplitude is determined from an energy integral analysis. The models have been applied to the study of an incompressible mixing layer. For both models, predictions of the mean flow developed are made. In the second model, predictions of the time-dependent motion of the large-scale structures in the mixing layer are made. The predictions show good agreement with experimental observations.

Liou, William W.↗
DOE OSTI2024

Fluid-Dynamic Mechanisms Underlying Wind Turbine Wake Control with Strouhal-Timed Actuation

A reduction in wake effects in large wind farms through wake-aware control has considerable potential to improve farm efficiency. This work examines the success of several emerging, empirically derived control methods that modify wind turbine wakes (i.e., the pulse method, helix method, and related methods) based on Strouhal numbers on the O(0.3). Drawing on previous work in the literature for jet and bluff-body flows, the analyses leverage the normal-mode representation of wake instabilities to characterize the large-scale wake meandering observed in actuated wakes. Idealized large-eddy simulations (LES) using an actuator-line representation of the turbine blades indicate that the n=0 and ±1 modes, which correspond to the pulse and helix forcing strategies, respectively, have faster initial growth rates than higher-order modes, suggesting these lower-order modes are more appropriate for wake control. Exciting these lower-order modes with periodic pitching of the blades produces increased modal growth, higher entrainment into the wake, and faster wake recovery. Modal energy gain and the entrainment rate both increase with streamwise distance from the rotor until the intermediate wake. This suggests that the wake meandering dynamics, which share close ties with the relatively well-characterized meandering dynamics in jet and bluff-body flows, are an essential component of the success of wind turbine wake control methods. A spatial linear stability analysis is also performed on the wake flows and yields insights on the modal evolution. In the context of the normal-mode representation of wake instabilities, these findings represent the first literature examining the characteristics of the wake meandering stemming from intentional Strouhal-timed wake actuation, and they help guide the ongoing work to understand the fluid-dynamic origins of the success of the pulse, helix, and related methods.

17 WIND ENERGY↗
DOE OSTI2020

The Role of Multiscale Interaction in Tropical Cyclogenesis and Its Predictability in Near-Global Aquaplanet Cloud-Resolving Simulations

Tropical cyclogenesis (TCG) is a multiscale process that involves interactions between large-scale circulation and small-scale convection. A near-global aquaplanet cloud-resolving model (NGAqua) with 4-km horizontal grid spacing that produces tropical cyclones (TCs) is used to investigate TCG and its predictability. This study analyzes an ensemble of three 20-day NGAqua simulations, with initial white-noise perturbations of low-level humidity. TCs develop spontaneously from the northern edge of the intertropical convergence zone (ITCZ), where large-scale flows and tropical convection provide necessary conditions for barotropic instability. Zonal bands of positive low-level absolute vorticity organize into cyclonic vortices, some of which develop into TCs. A new algorithm is developed to track the cyclonic vortices. A vortex-following framework analysis of the low-level vorticity budget shows that vertical stretching of absolute vorticity due to convective heating contributes positively to the vorticity spinup of the TCs. A case study and composite analyses suggest that sufficient humidity is key for convective development. TCG in these three NGAqua simulations undergoes the same series of interactions. The locations of cyclonic vortices are broadly predetermined by planetary-scale circulation and humidity patterns associated with ITCZ breakdown, which are predictable up to 10 days. Whether and when the cyclonic vortices become TCs depend on the somewhat more random feedback between convection and vorticity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗