Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear stability analysis”

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 181 records · Page 10

Electrohydrodynamic Stability of a Liquid Bridge: The "ALEX" Experiment

To provide insight into the roles of electrical forces, experiments on the stability of a liquid bridge were carried out during the 1996 Life And Microgravity Science Mission on the space shuttle Columbia. In terrestrial laboratories a Plateau configuration (where the bridge is surrounded by a matched density liquid) is necessary to avoid deformation due to buoyancy. This complicates the electrical boundary conditions, since charge is transported across the liquid-liquid interface. In the microgravity environment, a cylindrical bridge can be deployed in a gas which considerably simplifies the boundary condition. Nevertheless, to provide a tie-in to terrestrial experiments, two-phase experiments were carried out. The agreement with previous work was excellent. Then several experiments were conducted with a bridge deployed in a dielectric gas, SF6. In experiments with steady fields, it was found that the bridge was less stable than predicted by a linearized stability analysis using the Taylor-Melcher leaky dielectric model.

Burcham, C. L.↗

Noise Sources in a Low-Reynolds-Number Turbulent Jet at Mach 0.9

The mechanisms of sound generation in a Mach 0.9, Reynolds number 3600 turbulent jet are investigated by direct numerical simulation. Details of the numerical method are briefly outlined and results are validated against an experiment at the same flow conditions. Lighthill's theory is used to define a nominal acoustic source in the jet, and a numerical solution of Lighthill's equation is compared to the simulation to verify the computational procedures. The acoustic source is Fourier transformed in the axial coordinate and time and then filtered in order to identify and separate components capable of radiating to the far field. This procedure indicates that the peak radiating component of the source is coincident with neither the peak of the full unfiltered source nor that of the turbulent kinetic energy. The phase velocities of significant components range from approximately 5% to 50% of the ambient sound speed which calls into question the commonly made assumption that the noise sources convect at a single velocity. Space-time correlations demonstrate that the sources are not acoustically compact in the streamwise direction and that the portion of the source that radiates at angles greater than 45 deg. is stationary. Filtering non-radiating wavenumber components of the source at single frequencies reveals that a simple modulated wave forms for the source, as might be predicted by linear stability analysis. At small angles from the jet axis the noise from these modes is highly directional, better described by an exponential than a standard Doppler factor.

Freund, Jonathan B.↗

Boundary Layer Transition on Slender Cones in Conventional and Low Disturbance Mach 6 Wind Tunnels

An experimental investigation was conducted on a 5-degree half-angle cone and a 5-degree half-angle flared cone in a conventional Mach 6 wind tunnel to examine the effects of facility noise on boundary layer transition. The influence of tunnel noise was inferred by comparing transition onset locations determined from the present test to that previously obtained in a Mach 6 low disturbance quiet tunnel. Together, the two sets of experiments are believed to represent the first direct comparison of transition onset between a conventional and a low disturbance wind tunnel using a common test model and transition detection technique. In the present conventional hypersonic tunnel experiment, separate measurements of heat transfer and adiabatic wall temperatures were obtained on the conical models at small angles of attack over a range of Reynolds numbers, which resulted in laminar, transitional, and turbulent flow. Smooth model turbulent heating distributions are compared to that obtained with transition forced via discrete surface roughness. The model nosetip radius was varied to examine the effects of bluntness on transition onset. Despite wall to total temperature differences between the transient heating measurements and the adiabatic wall temperature measurement, the two methods for determining sharp cone transition onset generally yielded equivalent locations. In the 'noisy' mode of the hypersonic low disturbance tunnel, transition onset occurred earlier than that measured in the conventional hypersonic tunnel, suggesting higher levels of freestream acoustic radiation relative to the conventional tunnel. At comparable freestream conditions, the transition onset Reynolds number under low disturbance conditions was a factor of 1.3 greater than that measured on flared cone in the LaRC conventional hypersonic tunnel and a factor of 1.6 greater that the flared cone run in the low disturbance tunnel run 'noisy'. Navier-Stokes mean flow computations and linear stability analysis were conducted to assess the experimental results and have indicated N factors associated with sharp flared cone transition onset to be approximately a factor of 2 lower than that inferred from the corresponding low disturbance tunnel measurements.

Horvath, Thomas J.↗

Rayleigh-Benard Instability in a Vertical Cylinder with a Rotating Magnetic Field

This paper presents a linear stability analysis for the Rayleigh-Benard convection in a finite-length, vertical cylinder with a rotating magnetic field. The vertical wall of the cylinder is adiabatic, and planar top and bottom walls are isothermal with a higher temperature at the bottom. The stabilizing effects of the rotating magnetic field are studied for four values of the Prandtl number. Results for one Prandtl number are compared to previously published experimental results.

Walker, J. S.↗

Effects of Gravity on Bubble Formation in an Annular Jet

The effects of gravity on the bubble formation in an annular jet were studied. The experiments were conducted in the 2.2-second drop tower at the NASA Glenn Research Center. Terrestrial gravity experiments were conducted at the Fluid Dynamics Research Laboratory at the University of Oklahoma. Stainless steel tubing with inner diameters of 1/8" (gas inner annulus) and 5/16" (liquid outer annulus) served as the injector. A rectangular test section, 6" x 6" x 14" tall, made out of half-inch thick Lexan was used. Images of the annular jet were acquired using a high-speed camera. The effects of gravity and varying liquid and gas flow rates on bubble size, wavelength, and breakup length were documented. In general, the bubble diameter was found to be larger in terrestrial gravity than in microgravity for varying Weber numbers (0.05 - 0.16 and 5 - 11) and liquid flow rates (1.5 ft/s - 3.0 ft/s). The wavelength was found to be larger in terrestrial gravity than in microgravity, but remained constant for varying Weber numbers. For low Weber numbers (0.05 - 0.16), the breakup length in microgravity was significantly higher than in terrestrial gravity. Comparison with linear stability analysis showed estimated bubble sizes within 9% of experimental bubble sizes. Bubble size compared to other terrestrial gravity experiments with same flow conditions showed distinct differences in bubble size, which displayed the importance of injector geometry on bubble formation.

Koepp, R. A.↗

An Evaluation of the Roll-Rate Stabilization System of the Sidewinder Missile at Mach Numbers from 0.9 to 2.3

A linear stability analysis and flight-test investigation has been performed on a rolleron-type roll-rate stabilization system for a canard-type missile configuration through a Mach number range from 0.9 to 2.3. This type damper provides roll damping by the action of gyro-actuated uncoupled wing-tip ailerons. A dynamic roll instability predicted by the analysis was confirmed by flight testing and was subsequently eliminated by the introduction of control-surface damping about the rolleron hinge line. The control-surface damping was provided by an orifice-type damper contained within the control surface. Steady-state rolling velocities were at all times less than 1 radian per second between the Mach numbers of 0.9 to 2.3 on the configurations tested. No adverse longitudinal effects were experienced in flight because of the tendency of the free-floating rollerons to couple into the pitching motion at the low angles of attack and disturbance levels investigated herein after the introduction of control-surface damping.

Nason, Martin L.↗

Transition Within Leeward Plane of Axisymmetric Bodies at Incidence in Supersonic Flow

Boundary layer transition along the leeward symmetry plane of axisymmetric bodies at nonzero angle of attack in supersonic flow was investigated experimentally and numerically as part of joint research between the Japan Aerospace Exploration Agency (JAXA) and National Aeronautics and Space Administration (NASA). Transition over four axisymmetric bodies (namely, Sears-Haack body, semi-Sears-Haack body, straight cone and flared cone) with different axial pressure gradients was measured in two different facilities with different unit Reynolds numbers. The semi-Sears-Haack body and flared cone were designed at JAXA to broaden the range of axial pressure distributions. For a body shape with an adverse pressure gradient (i.e., flared cone), the experimentally measured transition patterns show an earlier transition location along the leeward symmetry plane in comparison with the neighboring azimuthal locations. For nearly zero pressure gradient (i.e.,straight cone), this feature is only observed at the larger unit Reynolds number. Later transition along the leeward plane was observed for the remaining two body shapes with a favorable pressure gradient. The observed transition patterns are only partially consistent with the numerical predictions based on linear stability analysis. Additional measurements are used in conjunction with the stability computations to explore the phenomenon of leeward line transition and the underlying transition mechanism in further detail.

Tokugawa, Naoko↗

Crossflow Instability on a Wedge-Cone at Mach 3.5

As a follow-on activity to the HyBoLT flight experiment, a six degree half angle wedge-cone model at zero angle of attack has been employed to experimentally and computationally study the boundary layer crossflow instability at Mach 3.5 under low disturbance freestream conditions. Computed meanflow and linear stability analysis results are presented along with corresponding experimental Pitot probe data. Using a model-mounted probe survey apparatus, data acquired to date show a well defined stationary crossflow vortex pattern on the flat wedge surface. This effort paves the way for additional detailed, calibrated flow field measurements of the crossflow instability, both stationary and traveling modes, and transition-to-turbulence under quiet flow conditions as a means of validating existing stability theory and providing a foundation for dynamic flight instrumentation development.

Beeler, George B.↗

Assessment of Transition Modeling Capability in OVERFLOW with Emphasis on Swept-Wing Configurations

In preparation for comparisons with data obtained from the recently concluded experiments in the National Transonic Facility at the NASA Langley Research Center on the new common research model with natural laminar flow (CRM-NLF), an assessment of the current transition modeling capability in the NASA OVERFLOW 2.2o code has been carried out. A combination of the available experimental data and linear stability analysis is used to evaluate the accuracy and robustness of these models for selected swept-wing type configurations, with significant crossflow. An additional goal for this work involves providing a comparative assessment of the relevant transition models in the context of a single flow solver and identifying model limitations as well as the potential for future improvements that would help strengthen the physical basis of such transition models. Included in this investigation is an assessment of the sensitivities of the underlying transition models to grid resolution (wall-normal, as well as streamwise and spanwise) and the values of extra input parameters such as the level of surface roughness and freestream turbulence variables. The flow configurations targeted in this assessment include the NASA NLF(2)-0415 swept-wing configuration, the sickle-shaped wing introduced by the Technical University of Braunschweig, and the wing-body configuration of the CRM model from the fourth and fifth AIAA CFD Drag Prediction Workshops.

Freestream velocity↗

Laminar to Turbulence Transition in Boundary Layers due to Tripping Devices

Direct numerical simulations (DNS)and linear stability analyses were performed to identify the mechanisms by whichtwo-and three-dimensional trips promote early transition in subsonic boundary layers over a flat plate. Investigations were conducted for threeboundary-layer transition trip configurations:a two-dimensional smooth hump,a two-dimensional rectangular-shapedrod and a three-dimensional zig-zagtape. The smooth hump was used by Park et al.1in the linear and nonlinear parabolized stability (PSE) computations. Reynolds numbersbased on the tripheights are 710, 683,and 300 for the two-and three-dimensional trips, respectively. Simulations withthe smooth hump showed that the flow becomes stable upstream and strongly unstable downstream of the hump. Simulations with the two-dimensional rectangular trip showed that the flow separates at the top of the trip and forms a long recirculation zonedownstream. However, the flow remains steady in the simulation without any external disturbances. Linear stability analysis showed that the inflectional profile in the recirculation zone is strongly unstable to frequencies in the range of 100-300 Hz. Simulations with the external disturbances forced at the frequency of 200 Hz showed that disturbances grow strongly inside the recirculation zone as observed in the Klebanoff2 experiment. Three-dimensional simulations with secondary disturbances of fundamental, subharmonic, and oblique types produced turbulent flows shortlydownstream. The three-dimensional zig-zag trip with a small roughness Reynolds number lead to early breakdown to a turbulent flow.

Stability↗

Automatic Boundary-Layer Adaptation of Structured Grids in VULCAN-CFD

In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.

boundary layer transition↗

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

Boundary layer transition↗

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗

Liapunov functions for non-linear difference equation stability analysis.

Liapunov functions to determine the stability of non-linear autonomous difference equations can be developed through the use of auxiliary exact difference equations. For this purpose definitions are introduced for the gradient of an implicit function of a discrete variable, a principal sum, a definite sum and an exact difference equation, and a theorem for exactness of a difference form is proved. Examples illustrate the procedure.

Park, K. E.↗

Nonparallel stability of three-dimensional compressible boundary layers. Part 1: Stability analysis

A compressible linear stability theory is presented for nonparallel three-dimensional boundary-layer flows, taking into account the normal velocity component as well as the streamwise and spanwise variations of the basic flow. The method of multiple scales is used to account for the nonparallelism of the basic flow, and equations are derived for the spatial evolution of the disturbance amplitude and wavenumber. The numerical procedure for obtaining the solution of the nonparallel problem is outlined.

El-Hady, N. M.↗

Towards Stability Analysis of Jump Linear Systems with State-Dependent and Stochastic Switching

This paper analyzes the stability of hierarchical jump linear systems where the supervisor is driven by a Markovian stochastic process and by the values of the supervised jump linear system s states. The stability framework for this class of systems is developed over infinite and finite time horizons. The framework is then used to derive sufficient stability conditions for a specific class of hybrid jump linear systems with performance supervision. New sufficient stochastic stability conditions for discrete-time jump linear systems are also presented.

Tejada, Arturo↗