SEARCH · Engineering Papers
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.
Nonsteady flow conditions in subsonic section of Laval nozzle to investigate nonlinear transverse combustion instability in liquid propellant rocket motors
Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing and the mode selected for forcing corresponds to the most amplified secondary instability mode that, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. An inlet boundary condition is carefully designed to allow for accurate injection of instability wave modes and minimize acoustic reflections at numerical boundaries. Nonlinear parabolized stability equation (PSE) predictions compare well with the DNS in terms of modal amplitudes and modal shape during the strongly nonlinear phase of the secondary instability mode. During the transition process, the skin friction coefficient rises rather rapidly and the wall-shear distribution shows a sawtooth pattern that is analogous to the previously documented surface flow visualizations of transition due to stationary crossflow instability. Fully turbulent features are observed in the downstream region of the flow.
As the result of global, non-parallel flow stability analysis the single value of the disturbance growth-rate and respective frequency is obtained. This complex value characterizes the stability of the whole flow configuration and is not referred to any particular flow pattern. The global analysis assures that all the flow elements (wake, boundary and shear layer) are taken into account. The physical phenomena connected with the wake instability are properly reproduced by the global analysis. This enhances the investigations of instability of any 2-D flows, including ones in which the boundary layer instability effects are known to be of dominating importance. Assuming fully 2-D disturbance form, the global linear stability problem is formulated. The system of partial differential equations is solved for the eigenvalues and eigenvectors. The equations, written in the pure stream function formulation, are discretized via FDM using a curvilinear coordinate system. The complex eigenvalues and corresponding eigenvectors are evaluated by an iterative method. The investigations performed for various Reynolds numbers emphasize that the wake instability develops into the Karman vortex street. This phenomenon is shown to be connected with the first mode obtained from the non-parallel flow stability analysis. The higher modes are reflecting different physical phenomena as for example Tollmien-Schlichting waves, originating in the boundary layer and having the tendency to emerge as instabilities for the growing Reynolds number. The investigations are carried out for a circular cylinder, oblong ellipsis and airfoil. It is shown that the onset of the wake instability, the waves in the boundary layer, the shear layer instability are different solutions of the same eigenvalue problem, formulated using the non-parallel theory. The analysis offers large potential possibilities as the generalization of methods used till now for the stability analysis.
This paper (the first in a series) focuses on using active-control methods to maintain laminar flow in a region of the flow in which the natural instabilities, if left unattended, lead to turbulent flow. The authors review previous studies that examine wave cancellation (currently the most prominent method) and solve the unsteady, nonlinear Navier-Stokes equations to evaluate this method of controlling instabilities. It is definitively shown that instabilities are controlled by the linear summation of waves (i.e., wave cancellation). Although a mathematically complete method for controlling arbitrary instabilities has been developed, the review, duplication, and physical explanation of previous studies are important steps for providing an independent verification of those studies, for establishing a framework for the work which will involve automated transition control, and for detailing the phenomena by-which the automated studies can be used to expand knowledge of flow control.
A spanwise core instability of the primary vortices in a plane mixing layer has been identified as a viable mechanism for the generation of streamwise vortical structures far downstream of the flow partition. The core instability is excited by a time-harmonic wave train with spanwise phase variation, synthesized by a mosaic of individually controlled surface heaters flush mounted on the flow partition. The flow is visualized in planes of constant cross-stream elevation by means of a schlieren system. As a result of the forcing, the primary vortices undergo spanwise deformations that induce secondary vortical structures, the shape and strength of which depends on the magnitude of the prescribed spanwise phase variation. The appearance of small-scale structures within the large coherent vortices suggests that core instability is an important contributor to mixing.
Laminar-turbulent transition in shear flows is still an enigma in the area of fluid mechanics. The conventional explanation of the phenomenon is based on the instability of the shear flow with respect to infinitesimal disturbances. The conventional hydrodynamic stability theory deals with the analysis of normal modes that might be unstable. The latter circumstance is accompanied by an exponential growth of the disturbances that might lead to laminar-turbulent transition. Nevertheless, in many cases, the transition scenario bypasses the exponential growth stage associated with the normal modes. This type of transition is called bypass transition. An understanding of the phenomenon has eluded us to this day. One possibility is that bypass transition is associated with so-called algebraic (non-modal) growth of disturbances in shear flows. In the present work, an analysis of the optimal disturbances/streamwise vortices associated with the transient growth mechanism is performed for boundary layers in the presence of a streamwise pressure gradient. The theory will provide the optimal spacing of the control elements in the spanwise direction and their placement in the streamwise direction.
A steady shear flow is generated in an axisymmetric liquid bridge surrounded by a passive gas through the imposition of a temperature gradient along the bridge and the driving of the motion by thermocapillarity. This dynamic state is susceptible to convective instabilities, which can lead to propagating hydrothermal waves that feed on the underlying temperature gradients. The convective instabilities of axisymmetric return flow are presently given as functions of the liquid's Prandtl number and the surface Biot number of the interface.
The flow about a slender body of revolution placed at incidence to an oncoming stream is numerically investigated for angles of attack ranging from 20 to 80 degrees and a Reynolds number of 200,000, based on freestream conditions and maximum body diameter. Navier-Stokes computations are applied to investigate the phenomena governing the onset of vortex asymmetry. Time-accurate solutions are given for an ogive-cylinder body with and without a space-fixed time-invariant disturbance added near the tip. At a 40-degree angle of attack the flow is found to be steady but becomes highly asymmetric when the disturbance is added. The level of asymmetry is dependent upon the size and location of the disturbance. It is suggested that the origin of the asymmetry is a convective-type instability of the originally symmetric flow. For higher angles of attack (60 and 80 degrees) the flow around the cylindrical part becomes unsteady, and vortex shedding is observed. It is suggested that the origin of flow unsteadiness and vortex shedding in the wake is an absolute-type instability of the originally steady flow.
The intent of the original proposal and the work carried out under the grant was to explore the instabilities of turbulent free shear flows of an incompressible fluid. The purpose for doing so was that coherent instabilities are susceptible to control via upstream forcing and because at sufficient finite amplitudes, the coherent modes enter into nonlinear interactions with the mean motion and thus, through such nonlinear interactions, the control of mixing and of mean flow spreading can be effected through the control of the coherent structures.
Parabolized stability equations (PSE) are used in a variational approach to study the optimal, non-modal disturbance growth in a Mach 3 at plate boundary layer and a Mach 6 circular cone boundary layer. As noted in previous works, the optimal initial disturbances correspond to steady counter-rotating streamwise vortices, which subsequently lead to the formation of streamwise-elongated structures, i.e., streaks, via a lift-up effect. The nonlinear evolution of the linearly optimal stationary perturbations is computed using the nonlinear plane-marching PSE for stationary perturbations. A fully implicit marching technique is used to facilitate the computation of nonlinear streaks with large amplitudes. To assess the effect of the finite-amplitude streaks on transition, the linear form of plane- marching PSE is used to investigate the instability of the boundary layer flow modified by spanwise periodic streaks. The onset of bypass transition is estimated by using an N- factor criterion based on the amplification of the streak instabilities. Results show that, for both flow configurations of interest, streaks of sufficiently large amplitude can lead to significantly earlier onset of transition than that in an unperturbed boundary layer without any streaks.
Results of computer code time dependent solutions of the two dimensional compressible Navier-Stokes equations and the results of independent experiments are compared to verify the Mach number range for instabilities in the transonic flow field about a 14 percent thick biconvex airfoil at an angle of attack of 0 deg and a Reynolds number of 7 million. The experiments were conducted in a transonic, slotted wall wind tunnel. The computer code included an algebraic eddy viscosity turbulence model developed for steady flows, and all computations were made using free flight boundary conditions. All of the features documented experimentally for both steady and unsteady flows were predicted qualitatively; even with the above simplifications, the predictions were, on the whole, in good quantitative agreement with experiment. In particular, predicted time histories of shock wave position, surface pressures, lift, and pitching moment were found to be in very good agreement with experiment for an unsteady flow. Depending upon the free stream Mach number for steady flows, the surface pressure downstream of the shock wave or the shock wave location was not well predicted.
It has long been recognized that flow in the melt can have a profound influence on the dynamics of a solidifying interface and hence the quality of the solid material. In particular, flow affects the heat and mass transfer, and causes spatial and temporal variations in the flow and melt composition. This results in a crystal with nonuniform physical properties. Flow can be generated by buoyancy, expansion or contraction upon phase change, and thermo-soluto capillary effects. In general, these flows can not be avoided and can have an adverse effect on the stability of the crystal structures. This motivates crystal growth experiments in a microgravity environment, where buoyancy-driven convection is significantly suppressed. However, transient accelerations (g-jitter) caused by the acceleration of the spacecraft can affect the melt, while convection generated from the effects other than buoyancy remain important. Rather than bemoan the presence of convection as a source of interfacial instability, Hurle in the 1960s suggested that flow in the melt, either forced or natural convection, might be used to stabilize the interface. Delves considered the imposition of both a parabolic velocity profile and a Blasius boundary layer flow over the interface. He concluded that fast stirring could stabilize the interface to perturbations whose wave vector is in the direction of the fluid velocity. Forth and Wheeler considered the effect of the asymptotic suction boundary layer profile. They showed that the effect of the shear flow was to generate travelling waves parallel to the flow with a speed proportional to the Reynolds number. There have been few quantitative, experimental works reporting on the coupling effect of fluid flow and morphological instabilities. Huang studied plane Couette flow over cells and dendrites. It was found that this flow could greatly enhance the planar stability and even induce the cell-planar transition. A rotating impeller was buried inside the sample cell, driven by an outside rotating magnet, in order to generate the flow. However, it appears that this was not a well-controlled flow and may also have been unsteady. In the present experimental study, we want to study how a forced parallel shear flow in a Hele-Shaw cell interacts with the directionally solidifying crystal interface. The comparison of experimental data show that the parallel shear flow in a Hele-Shaw cell has a strong stabilizing effect on the planar interface by damping the existing initial perturbations. The flow also shows a stabilizing effect on the cellular interface by slightly reducing the exponential growth rate of cells. The left-right symmetry of cells is broken by the flow with cells tilting toward the incoming flow direction. The tilting angle increases with the velocity ratio. The experimental results are explained through the parallel flow effect on lateral solute transport. The phenomenon of cells tilting against the flow is consistent with the numerical result of Dantzig and Chao.
A spectral algorithm for simulating three-dimensional, incompressible, parallel shear flows is described. It applies to the channel, to the parallel boundary layer, and to other shear flows with one wall-bounded and two periodic directions. Simulations performed for the channel identify the one nonlinear interaction which is crucial to the secondary instability. Moreover, it is shown that subharmonic instability ensues for Reynolds numbers as low as 1500. Simulations performed for the heated water boundary layer reveal the operation of the secondary instability in this controlled flow.
The results of numerical model calculations are compared with space laboratory experiments for rotating, thermally driven flow in a hemispherical shell, and various flow regimes are described. A spherically symmetric body force, analogous to gravity, is imposed in the radial direction through the use of a dielectric fluid and an electrostatic potential difference across the gap. The spherical boundaries are maintained at constant temperature profiles (with the pole and the inner sphere being the warmer surfaces), and the equatorial wall is an insulator. Typical parameter combinations result in highly nonlinear, but laminar, flow. For weak enough buoyant forcing, the flow is axisymmetric. It consists of a single meridional cell, rising in warm latitudes (in this case, near the pole) and sinking in cool latitudes, with prograde flow in the equatorial region and near the inner hemisphere for other latitudes, and retrograde motion elsewhere. For fast rotation, the first transition due to the instability of this simple flow is also axisymmetric, consisting of rings of convection in the polar region, which propagate poleward. The first nonaxisymmetric convection occurs at Rayleigh numbers which increase with rotation rate. The form of the convection near the transition also depends upon the rate of rotation. Selected flow patterns near the transition as well as those beyond it are studied numerically. For those cases where there exist laboratory experiments with which to compare, the numerical and experimental results agree very well.
Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing derived from a two-dimensional, partial-differential-equation based eigenvalue computation; and the mode selected for forcing corresponds to the most amplified secondary instability mode which, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. Both the growth of the secondary instability wave and the resulting onset of laminar-turbulent transition are captured within the DNS computations. The growth of the secondary instability wave in the DNS solution compares well with linear secondary instability theory when the amplitude is small; the linear growth is followed by a region of reduced growth resulting from nonlinear effects before an explosive onset of laminar breakdown to turbulence. The peak fluctuations are concentrated near the boundary layer edge during the initial stage of transition, but rapidly propagates towards the surface during the process of laminar breakdown. Both time-averaged statistics and flow visualization based on the DNS reveal a sawtooth transition pattern that is analogous to previously documented surface flow visualizations of transition due to stationary crossflow instability. The memory of the stationary crossflow vortex is found to persist through the transition zone and well beyond the location of the maximum skin friction.
Boundary-layer transition and relaminarization may have a critical effect on the flow development about multielement high-lift systems of subsonic transport aircraft with swept wings. The purpose of this paper is a study of transition phenomena in the leading-edge region of the various elements of a high-lift system. The flow phenomena studied include transition of the attachment-line flow, relaminarization, and crossflow instability. The calculations are based on pressure distributions measured in flight on the NASA Transport Systems Research Vehicle (Boeing 737-100) at a wing station where the flow approximated infinite swept wing conditions. The results indicate that significant regions of laminar flow can exist on all flap elements in flight. In future flight experiments the extent of these regions will be measured, and the transition mechanisms and the effect of laminar flow on the high-lift characteristics of the multi-element system will be further explored.
This paper reports a series of wind tunnel tests simulating the near-field behavior of the Space Shuttle Orbiter Boundary Layer Transition Detailed Test Objective (BLT DTO) flight experiment. Hypersonic flow over a flat plate with an attached BLT DTO-shaped trip was tested in a Mach 10 wind tunnel. The sharp-leading-edge flat plate was oriented at an angle of 20 degrees with respect to the freestream flow, resulting in post-shock edge Mach number of approximately 4. The flowfield was visualized using nitric oxide (NO) planar laser-induced fluorescence (PLIF). Flow visualizations were performed at 10 Hz using a wide-field of view and high-resolution NO PLIF system. A lower spatial resolution and smaller field of view NO PLIF system visualized the flow at 500 kHz, which was fast enough to resolve unsteady flow features. At the lowest Reynolds number studied, the flow was observed to be laminar and mostly steady. At the highest Reynolds number, flow visualizations showed streak instabilities generated immediately downstream of the trip. These instabilities transitioned to unsteady periodic and spatially irregular structures downstream. Quantitative surface heating imagery was obtained using the Temperature Sensitive Paint (TSP) technique. Comparisons between the PLIF flow visualizations and TSP heating measurements show a strong correlation between flow patterns and surface heating trends.
The stability of the flow of a viscous incompressible fluid over a curved compliant wall to longitudinal Goertler vortices is investigated. The compliant wall is modeled by a particularly simple equation relating the induced wall displacement to the pressure in the overlying fluid. Attention is restricted to the large Goertler number regime; this regime being appropriate to the most unstable Goertler mode. The effect of wall compliance on this most unstable mode is investigated.