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,441 records · Page 80

Whirl/whip demonstration

Fluid flow in bearings and seals, set in motion by shaft rotation, generates dynamic forces which may result in a well recognized instability known as whirl and whip. These are lateral, forward precessional, self excited, subsynchronous vibrations in which the amplitude may vary from very small to nearly the limit of the bearing or seal clearances. Oil whirl in lubricated bearings, in particular, typically occurs at somewhat less than half rotative speed. As the rotative speed increases, the frequency relationship remains constant until the whirl frequency approaches the first balance resonance. Now the whirl is smoothly replaced by whip at a nearly constant frequency asymptotically approaching first balance resonance, independent of increasing rotative speed. Changes in bearing/seal radial loading can permit, prevent, or eliminate this instability. The oil whirl/whip rig demonstrates the effects of fluid dynamic forces generated by the rotating shaft. At low rotative speeds, this produces changes of the journal static equilibrium position within the bearing. The demonstrator shows the relationship between any load direction and the average journal equilibrium position. At higher rotative speeds, the instability threshold is observed as a function of unidirectional radial load, unbalance, and rotor configuration.

Grissom, R.↗

Singularities of the Euler equation and hydrodynamic stability

Equations governing the motion of a specific class of singularities of the Euler equation in the extended complex spatial domain are derived. Under some assumptions, it is shown how this motion is dictated by the smooth part of the complex velocity at a singular point in the unphysical domain. These results are used to relate the motion of complex singularities to the stability of steady solutions of the Euler equation. A sufficient condition for instability is conjectured. Several examples are presented to demonstrate the efficacy of this sufficient condition which include the class of elliptical flows and the Kelvin-Stuart Cat's Eye.

Tanveer, S.↗

Singularities of the Euler equation and hydrodynamic stability

Equations governing the motion of a specific class of singularities of the Euler equation in the extended complex spatial domain are derived. Under some assumptions, it is shown how this motion is dictated by the smooth part of the complex velocity at a singular point in the unphysical domain. These results are used to relate the motion of complex singularities to the stability of steady solutions of the Euler equation. A sufficient condition for instability is conjectured. Several examples are presented to demonstrate the efficacy of this sufficient condition which include the class of elliptical flows and the Kelvin-Stuart Cat's Eye.

Tanveer, S.↗

Higher-order accurate Osher schemes with application to compressible boundary layer stability

Two fourth order accurate Osher schemes are presented which maintain higher order accuracy on nonuniform grids. They use either a conservative finite difference or finite volume discretization. Both methods are successfully used for direct numerical simulations of flat plate boundary layer instability at different Mach numbers. Results of growth rates of Tollmien-Schlichting waves compare well with direct simulations of incompressible flow and for compressible flow with results obtained by solving the parabolic stability equations.

Van Der Vegt, J. J. W.↗

Higher-order accurate Osher schemes with application to compressible boundary layer stability

Two fourth order accurate Osher schemes are presented which maintain higher order accuracy on nonuniform grids. They use either a conservative finite difference or finite volume discretization. Both methods are successfully used for direct numerical simulations of flat plate boundary layer instability at different Mach numbers. Results of growth rates of Tollmien-Schlichting waves compare well with direct simulations of incompressible flow and for compressible flow with results obtained by solving the parabolic stability equations.

Vandervegt, J. J. W.↗

Boundary-layer receptivity due to distributed surface imperfections of a deterministic or random nature

Acoustic receptivity of a Blasius boundary layer in the presence of distributed, 2D surface irregularities is investigated analytically. It is shown that, out of the entire spatial spectrum of the surface irregularities, only a narrow band of Fourier components can lead to an efficient conversion of the acoustic input at any given frequency to an unstable eigenmode of the boundary-layer flow. Random irregularities which are spatially homogeneous in nature, and also possess a continuous spectrum in the streamwise direction, lead to instability amplitudes that are intermediate to those caused by the periodic and isolated irregularities, respectively, being, in fact, of the same order as the geometrical mean of the amplitudes in the latter two cases. A physical explanation for these asymptotic scalings is given, in addition to providing an analytical expression for the expected value of the instability amplitude for an ensemble of statistically irregular surfaces with random phase distributions.

Choudhari, Meelan↗

Testing the frozen flow approximation

We investigate the accuracy of the frozen-flow approximation (FFA), recently proposed by Matarrese, et al. (1992), for following the nonlinear evolution of cosmological density fluctuations under gravitational instability. We compare a number of statistics between results of the FFA and n-body simulations, including those used by Melott, Pellman & Shandarin (1993) to test the Zel'dovich approximation. The FFA performs reasonably well in a statistical sense, e.g. in reproducing the counts-in-cell distribution, at small scales, but it does poorly in the crosscorrelation with n-body which means it is generally not moving mass to the right place, especially in models with high small-scale power.

Lucchin, Francesco↗

Increased Stability in Crystal Growth Kinetics in Response to Bulk Transport Enhancement

High-resolution interferometry on facets of lysozyme crystals growing under steady external conditions revealed intrinsic kinetics instabilities that occur as fluctuations in the spacing of growth steps. The amplitude of these fluctuations is reduced by forced solution flow. This increased stability is interpreted in terms of the ratio of the characteristic rates of lattice incorporation and bulk transport. Apparently a steady state is approached as an asymptote of a continuous family of unsteady states.

Vekilov, Peter G.↗

Lattice Boltzmann Method for 3-D Flows with Curved Boundary

In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.

Mei, Renwei↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

Noise Radiation by Instability Waves in Coaxial Jets

In this paper predictions are made for the noise radiation from supersonic coaxial jets. The noise in the downstream arc of a supersonic jet is dominated by highly directional radiation from the supersonically convecting large scale structures in the jet mixing layer. Since the mean flow is not described easily in terms of simple analytic functions, a numerical prediction is made for its development. The compressible Reynolds-averaged boundary layer equations in cylindrical polar coordinates are solved. A mixing length turbulence model is used. Empirical correlations are developed the effects of velocity and temperature ratio and Mach number. Both normal and inverted velocity profiles are considered. Comparisons with measurements for both single and coaxial jets show good agreement. The large scale structures are modeled as instability waves. The noise radiation generated by the instability waves is determined by a matching between the inner instability wave solution and the outer acoustic solution. Predictions are made for the differences between the noise radiated by coaxial jets with different operating conditions and a single equivalent jet with the same exit area, thrust, and mass-flow.

Dahl, Milo D.↗

Winds at the Mars 2020 Landing Site: 1. Near-Surface Wind Patterns at Jezero Crater

This is the first part of a two-part paper. NASA's Mars 2020 Perseverance rover measured winds on the Jezero crater floor close to the delta of an ancient river. A mostly repeatable diurnal cycle was observed and presented two regimes: (a) a convective regime, from dawn to sunset, with average easterly to southeasterly winds, during which maximum wind speeds were measured, and (b) a nighttime regime with westerly-northwesterly winds followed by a relatively calm period with highly variable wind directions as a function of sol and time of night. The timing and magnitude of the observed regimes are consistent with primary control by regional and local slope flows. Data suggest that the surface circulation at Jezero region in northern spring and summer is highly unaffected by large-scale circulation except during particular periods in the diurnal cycle or generally during dust storms, which is supported by MarsWRF model simulations. Consequently, the observed seasonal variability was weak. However, sol-to-sol and seasonal variability were measured, most of it during certain nighttime periods. Traveling waves consistent with baroclinic instability were clearly observed in surface winds at L_s ∼ 75°. The early MY36/2022A regional dust storm at L_s ∼ 153° disturbed the wind patterns with changes suggesting enhanced tidal flows. After sunset, the dust storm also produced detectable gravity wave activity, increasing the mixing in the nighttime planetary boundary layer during storm conditions. Inferred wind directions from dust devil movies strongly suggest that prevailing winds continued to be slope-driven during the late summer, fall and early winter seasons.

D. Viúdez-Moreiras↗

Shear stabilization of the capillary breakup of a cylindrical interface

A cylindrical interface containing a viscous liquid set into axial motion is subject to a capillary and to a surface-wave instability. Clues from previous studies suggest that, even though both mechanisms separately are destabilizing, under certain circumstances their mutual interaction can lead to a stable interface; shear can stabilize capillary breakup. Here, an axial flow through an annular cross section bounded on the inside by a rigid rod and on the outside by a deformable interface is considered. The competition between the two mechanisms is studied through the temporal growth of infinitesimal axisymmetric and nonaxisymmetric disturbances. This examination of temporal stability shows that, indeed, for geometries corresponding to thin annular layers both instabilities can be completely suppressed (disturbances of all wavelengths decay).

Russo, Mathew J.↗

Turbulent swirling jets with excitation

An existing cold-jet facility at NASA Lewis Research Center was modified to produce swirling flows with controllable initial tangential velocity distribution. Two extreme swirl profiles, i.e., one with solid-body rotation and the other predominated by a free-vortex distribution, were produced at identical swirl number of 0.48. Mean centerline velocity decay characteristics of the solid-body rotation jet flow exhibited classical decay features of a swirling jet with S - 0.48 reported in the literature. However, the predominantly free-vortex distribution case was on the verge of vortex breakdown, a phenomenon associated with the rotating flows of significantly higher swirl numbers, i.e., S sub crit greater than or equal to 0.06. This remarkable result leads to the conclusion that the integrated swirl effect, reflected in the swirl number, is inadequate in describing the mean swirling jet behavior in the near field. The relative size (i.e., diameter) of the vortex core emerging from the nozzle and the corresponding tangential velocity distribution are also controlling factors. Excitability of swirling jets is also investigated by exciting a flow with a swirl number of 0.35 by plane acoustic waves at a constant sound pressure level and at various frequencies. It is observed that the cold swirling jet is excitable by plane waves, and that the instability waves grow about 50 percent less in peak r.m.s. amplitude and saturate further upstream compared to corresponding waves in a jet without swirl having the same axial mass flux. The preferred Strouhal number based on the mass-averaged axial velocity and nozzle exit diameter for both swirling and nonswirling flows is 0.4.

Taghavi, Rahmat↗

Small-scale motions observed by aircraft in the tropical lower stratosphere - Evidence for mixing and its relationship to large-scale flows

Measurements of temperature and ozone from instrumented aircraft in the tropical lower stratosphere show the presence of small-scale disturbances generated by underlying cumulus convection and Kelvin-Helmholtz instability. The disturbances associated with underlying convection have peak-to-peak vertical parcel excursions of about 300 m. Flying conditions were smooth, suggesting an ensemble of gravity waves and little or no turbulent mixing. It is nevertheless possible that these waves break at other altitudes, leading to turbulent mixing and net fluxes of vertically stratified tracers. Disturbances attributed to KH instability implied vertical parcel excursions of 300-400 m. The disturbances coincided with rough flying conditions, suggesting turbulent mixing. A linear stability analysis of the atmospheric basic state defined by high-resolution radiosondes shows fastest growing waves with horizontal wavelengths of 1.4-1.8 km, consistent with the aircraft observations. The strong shears responsible for the KH instability are due to large-scale waves propagating into a region of small intrinsic frequency. Radiosonde observations show that the zonal length scale of these waves is about 1000 km.

Pfister, Leonhard↗

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.↗

Steady, subsonic, lifting surface theory for wings with swept, partial span, trailing edge control surfaces

A method for computing the lifting pressure distribution on a wing with partial span, swept control surfaces is presented. This method is valid within the framework of linearized, steady, potential flow theory and consists of using conventional lifting surface theory in conjuction with a flap pressure mode. The cause of a numerical instability that can occur during the quadrature of the flap pressure mode is discussed, and an efficient technique to eliminate the instability is derived. This technique is valid for both the flap pressure mode and regular pressure modes and could be used to improve existing lifting surface methods. Examples of the use of the flap pressure mode and comparisons among this method, other theoretical methods, and experiments are given. Discrepancies with experiment are indicated and candidate causes are presented. It is concluded that the method can lead to an efficient and accurate solution of the mathematical problem when a partial span, trailing edge flap is involved.

Medan, R. T.↗