Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

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 847 records · Page 47

Flow measurements in a model swirl combustor

Flow measurements in a model combustor composed of two confined coaxial swirling jets under noncombusting conditions are presented. Mean flow results are obtained for five different flow conditions to determine the effect of outer swirl on the recirculation zone (which is used for flame stabilization under combustion conditions). As the outer swirl magnitude is first decreased from maximum counter-swirl to zero and then increased to give co-swirl conditions, the size and the reverse flow velocity in the recirculation zone diminish. Detailed time mean and fluctuating flow measurements are obtained for a co-swirl and a counter-swirl condition with a directional pitot probe and hot-wire anenometry. For these two cases, recirculation zone occurs only with counter-swirl, near the exit of the inner jet. The recirculation zone is in the form of a one celled toroidal vortex having very low swirl velocity. Axial development of fluctuation levels, energy dissipation rates and turbulence length scales are described for the two flow conditions. Spectral analysis reveals periodic oscillations in both flows. The oscillations originate from the inner jet and up to 4 harmonics are observed. The fundamental frequencies are comparable to the rotational frequencies of the inner jets at the exit which are approximately under solid body rotation. The significance of the present results for the combustion process is discussed.

Vu, B. T.↗

On the stability of the ionopause of Venus

The stability of the Venus ionopause is examined in light of the importance of gravitation and curvature. Using a one-fluid approximation for the equation of motion of the plasma, and ignoring the effects of neutrals, a dispersion relation is obtained that includes the effects of the magnetic field, sheared plasma flow, buoyancy, centrifugal force and magnetic tension due to boundary curvature. It is found that buoyancy acts to neutralize the flute instability. As expected, the Kelvin-Helmholtz mode is the dominant instability over most of the dayside ionopause. The expected growth times of this mode are short in comparison with the wave-convection time over the boundary; the waves can grow and saturate quickly, producing a turbulent boundary that may affect electrodynamic coupling between the solar wind and ionospheric plasmas.

Elphic, R. C.↗

Analysis of rotordynamic coefficients of helically-grooved turbulent annular seals

An analysis for helically-grooved turbulent annular seals is developed to predict leakage and dynamic coefficients, as related to rotordynamics. The grooved surface pattern is formulated as an inhomogeneous directivity in surface shear stress. The zeroth-order equations define the steady-state leakage and the circumferential velocity development due to wall shear for a centered rotor position. The first-order equations define perturbations in the pressure and axial and circumferential velocity fields due to small motion of the rotor about the centered position. Numerical results are presented for proposed grooved seals in the High Pressure Oxygen Turbopump (HPTOP) of the Space Shuttle Main Engine (SSME) and for a water-pump application. The results show that an optimum helix angle exists from a rotordynamic stability viewpoint. Further, a properly designed helically-grooved stator is predicted to have pronounced stability advantages over other currently used seals.

Kim, C.-H.↗

Reynolds stress closure in jet flows using wave models

Numerical methods were developed that will form the computational part of the turbulence closure scheme. A wave model was developed for the two-dimensional shear layer. This configuration is being used as a test case for the closure schemes. Various numerical schemes were examined to give efficient solutions of the Rayleigh equation for this geometry. These include both spectral and finite difference methods. Secondly, numerical methods are under development to solve the non-separable Rayleigh equation. This solution is required for the closure scheme in more complex geometries. A model problem was used to assist in the algorithm development. Two-dimensional spectral methods and a hybrid spectral/finite difference technique were developed. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique were also examined. These solutions are forming the basis of a model for the shock structure in jets of arbitrary geometry.

Morris, P. J.↗

Receptivity of flat-plate boundary layer in a non-uniform free stream (vorticity normal to the plate)

Recent progress in both the linear and nonlinear aspects of stability theory has highlighted the importance of the receptivity problem. One of the most unclear aspects of receptivity study is the receptivity of boundary-layer flow normal to vortical disturbances. Some experimental and theoretical results permit the proposition that quasi-steady outer-flow vortical disturbances may trigger by-pass transition. In present work such interaction is investigated for vorticity normal to a leading edge. The interest in these types of vortical disturbances arise from theoretical work, where it was shown that small sinusoidal variations of upstream velocity along the spanwise direction can produce significant variations in the boundary-layer profile. In the experimental part of this work, such non-uniform flow was created and the laminar-turbulent transition in this flow was investigated. The experiment was carried out in a low-turbulence direct-flow wind tunnel T-361 at the Central Aerohydrodynamic Institute (TsAGI). The non-uniform flow was produced by laminar or turbulent wakes behind a wire placed normal to the plate upstream of the leading edge. The theoretical part of the work is devoted to studying the unstable disturbance evolution in a boundary layer with strongly non-uniform velocity profiles similar to that produced by outer-flow vorticity. Specifically, the Tollmien-Schlichting wave development in the boundary layer flow with spanwise variations of velocity is investigated.

Kogan, M. N.↗

Convergence Acceleration of Runge-Kutta Schemes for Solving the Navier-Stokes Equations

The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 can be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. This RK/implicit scheme is used as a smoother for multigrid. Fourier analysis is applied to determine damping properties. Numerical dissipation operators based on the Roe scheme, a matrix dissipation, and the CUSP scheme are considered in evaluating the RK/implicit scheme. In addition, the effect of the number of RK stages is examined. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. Turbulent flows over an airfoil and wing at subsonic and transonic conditions are computed. The effects of the cell aspect ratio on convergence are investigated for Reynolds numbers between 5:7 x 10(exp 6) and 100 x 10(exp 6). It is demonstrated that the implicit preconditioner can reduce the computational time of a well-tuned standard RK scheme by a factor between four and ten.

Swanson, Roy C., Jr.↗

A Flight Investigation of the Effects of Varied Lateral Damping on the Effectiveness of a Fighter Airplane as a Gun Platform

Flight tests were conducted at three conditions of damping and two atmospheric turbulence levels with fixed-reticle, gyro computing, and telescopic gunsights. Results of simulated strafing runs indicate that the gun-line dispersion could be expected to be decreased about 7 percent by increased lateral damping and to be increased about 85 percent by decreased damping and about 40 percent by rough air of the type encountered. Use of the telescopic sight indicated a 20-percent decrease in gun-line dispersion.

OPERATING PROBLEMS, PHYSIOLOGICAL↗

Modeling of Wake-vortex Aircraft Encounters

There are more people passing through the world's airports today than at any other time in history. With this increase in civil transport, airports are becoming capacity limited. In order to increase capacity and thus meet the demands of the flying public, the number of runways and number of flights per runway must be increased. In response to the demand, the National Aeronautics and Space Administration (NASA), in conjunction with the Federal Aviation Administration (FAA), airport operators, and the airline industry are taking steps to increase airport capacity without jeopardizing safety. Increasing the production per runway increases the likelihood that an aircraft will encounter the trailing wake-vortex of another aircraft. The hazard of a wake-vortex encounter is that heavy load aircraft can produce high intensity wake turbulence, through the development of its wing-tip vortices. A smaller aircraft following in the wake of the heavy load aircraft will experience redistribution of its aerodynamic load. This creates a safety hazard for the smaller aircraft. Understanding this load redistribution is of great importance, particularly during landing and take-off. In this research wake-vortex effects on an encountering 10% scale model of the B737-100 aircraft are modeled using both strip theory and vortex-lattice modeling methods. The models are then compared to wind tunnel data that was taken in the 30ft x 60ft wind tunnel at NASA Langley Research Center (LaRC). Comparisons are made to determine if the models will have acceptable accuracy when parts of the geometry are removed, such as the horizontal stabilizer and the vertical tail. A sensitivity analysis was also performed to observe how accurately the models could match the experimental data if there was a 10% error in the circulation strength. It was determined that both models show accurate results when the wing, horizontal stabilizer, and vertical tail were a part of the geometry. When the horizontal stabilizer and vertical tail were removed there were difficulties modeling the sideforce coefficient and pitching moment. With the removal of only the vertical tail unacceptable errors occurred when modeling the sideforce coefficient and yawing moment. Lift could not be modeled with either the full geometry or the reduced geometry attempts.

Smith, Sonya T.↗

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Zhang, Hai↗

Simulations of Boundary-Layer Transition

For incompressible benchmark flows, we have demonstrated the capability of the parabolized stability equations (PSE) to simulate the transition process in excellent agreement with microscopic experiments and direct Navier-Stokes simulations at modest computational cost. Encouraged by these results, we have developed the PSE methodology of three-dimensional boundary-layers in general curvilinear coordinates for the range from low to hypersonic speeds, and for both linear and nonlinear problems. For given initial and boundary conditions, the approach permits simulations from receptivity through linear and secondary instabilities into the late stages of transition where significant changes in skin friction and heat transfer coefficients occur. We have performed transition simulations for a variety of two- and three-dimensional similarity solutions and for realistic flows over swept wings at subsonic and supersonic speeds, the pressure ans suction side of turbine blades at low and medium turbulence levels, and over a blunt cone at Mach number Ma = 8. We present selected results for different transition mechanisms with emphasis on the late stage of transition and the evolution of wall-shear stress and heat transfer.

Herbert, Thorwald↗

Mid-Fidelity Computational Fluid Dynamics Analysis of the Elytron 4S UAV Concept

The Elytron 4S Unmanned Aerial Vehicle (UAV) Concept was developed to combine the advantages of fixed- and rotary-wing technology. The 4S Concept is a box-wing configuration with rotors mounted on a centrally located tiltwing. The UAV is intended to be capable of both conventional takeoff and landing (CTOL) and vertical takeoff and landing (VTOL), and is envisioned to excel in UAV performance because of the combined efficiency of fixed-wing aircraft and the hover and VTOL capabilities of regular drones or quadcopters. A mid-fidelity Unsteady Reynolds-Averaged Navier-Stokes (URANS) approach using Rotorcraft CFD (RotCFD) is performed to analyze and characterize the performance of the aircraft. The flow field is coupled with a rotor model based on blade-element momentum theory to model the 4S UAV rotors. Turbulence is modeled using a realizable k- turbulence model with special wall function. The code is used to generate aerodynamic forces and moments on the body at cruise conditions, and during VTOL. The results and their uncertainties are characterized, and an angle- of-attack and sideslip sweep are computed, both with and without rotors on. Simulations are compared with the wind tunnel tests in the 7- by 10-Foot U.S. Army Wind Tunnel at NASA Ames Research Center, performed in 2017. Results show promising comparison with experimental data, despite a late change in rotor size and rudder size of the physical model that cause the expected deviations from the simulation. A slight change in the net thrust value, when rotors are modeled, is observed because of the rotor diameter increase on the physical model. A noticeable difference in the directional stability was observed because of the increased rudder surface and added strakes. These changes were implemented to improve on the design as simulated, which is observed in the results. The simulation results paved the way to the first successful flight of the UAV Concept.

Koning, Witold J. K.↗

Design Exploration of a Transonic Cruise Slotted Airfoil

A knowledge-based aerodynamic design method, CDISC, has been leveraged for the computational design exploration of the cruise slotted airfoil as a drag-saving technology for single aisle, transonic transport aircraft. Aerodynamic predictions were generated using the NASA USM3D Reynolds-averaged Navier-Stokes flow solver, and laminar flow assessments were conducted using the NASA BLSTA3D boundary layer profile solver paired with the LASTRAC stability analysis and transition prediction software. The CDISC design method was used to parametrically vary several slotted-airfoil design variables to understand their impact on aerodynamic performance. In addition to establishing best practices for slotted airfoil design, it was observed that the historically noted skin-friction penalty for cruise slotted airfoils was attributed to the low-Reynolds-number boundary layer of the flap. At cruise, this skin-friction penalty negated any potential decrease in pressure drag enabled by the cruise slotted airfoil architecture. This observation led to the development of an Aft-Laminar Multi-element Airfoil concept that uses airfoil shaping to promote natural laminar flow on the flap. At cruise conditions, the concept is predicted to offer a 2.7% reduction in sectional drag relative to a fully turbulent supercritical airfoil. Off-design analyses showed that laminar flow could be maintained with near-cruise variations in angle of attack and Mach number, resulting in sustained aerodynamic efficiency improvements and a delay in the drag rise Mach number. A low-speed, high-lift analysis at takeoff conditions predicted a maximum lift coefficient of 2.4 could be achieved by use of a variable-camber leading edge on the main element and simple deflection of the flap. These computationally-predicted benefits have motivated future research toward an aft-laminar cruise slotted wing design as an incremental step toward the application of natural laminar flow technology on transonic commercial transports.

Brett R Hiller↗

Piloted Evaluation of Modernized Limited Authority Control Laws in the NASA-Ames Vertical Motion Simulator (VMS)

The UH-60 BLACK HAWK was designed in the 1970s, when the US Army primarily operated during the day in good visual conditions. Subsequently, the introduction of night-vision goggles increased the BLACK HAWK'S mission effectiveness, but the accident rate also increased. The increased accident rate is strongly tied to increased pilot workload as a result of a degradation in visual cues. Over twenty years of research in helicopter flight control and handling qualities has shown that these degraded handling qualities can be recovered by modifying the response type of the helicopter in low speed flight. Sikorsky Aircraft Corporation initiated a project under the National Rotorcraft Technology Center (NRTC) to develop modern flight control laws while utilizing the existing partial authority Stability Augmentation System (SAS) of the BLACK HAWK. This effort resulted in a set of Modernized Control Laws (MCLAWS) that incorporate rate command and attitude command response types. Sikorsky and the US Army Aeroflightdynamics Directorate (AFDD) conducted a piloted simulation on the NASA-Ames Vertical h4otion Simulator, to assess potential handling qualities and to reduce the risk of subsequent implementation and flight test of these modern control laws on AFDD's EH-60L helicopter. The simulation showed that Attitude Command Attitude Hold control laws in pitch and roll improve handling qualities in the low speed flight regime. These improvements are consistent across a range of mission task elements and for both good and degraded visual environments. The MCLAWS perform better than the baseline UH-60A control laws in the presence of wind and turbulence. Finally, while the improved handling qualities in the pitch and roll axis allow the pilot to pay more attention to the vertical axis and hence altitude performance also improves, it is clear from pilot comments and altitude excursions that the addition of an Altitude Hold function would further reduce workload and improve overall handling qualities of the aircraft.

Sahasrabudhe, Vineet↗

Rotordynamic coefficients for stepped labyrinth gas seals

The basic equations are derived for compressible flow in a stepped labyrinth gas seal. The flow is assumed to be completely turbulent in the circumferential direction where the friction factor is determined by the Blasius relation. Linearized zeroth and first-order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth-order pressure distribution is found by satisfying the leakage equation while the circumferential velocity distribution is determined by satisfying the momentum equations. The first order equations are solved by a separation of variables solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients. The results of this analysis are presented in the form of a parametric study, since there are no known experimental data for the rotordynamic coefficients of stepped labyrinth gas seals. The parametric study investigates the relative rotordynamic stability of convergent, straight and divergent stepped labyrinth gas seals. The results show that, generally, the divergent seal is more stable, rotordynamically, than the straight or convergent seals. The results also show that the teeth-on-stator seals are not always more stable, rotordynamically, then the teeth-on-rotor seals as was shown by experiment by Childs and Scharrer (1986b) for a 15 tooth seal.

Scharrer, Joseph K.↗

Numerical investigation of laminar-turbulent transition in a flat plate wake

Lamina-turbulent transition of high-deficit flat plate wakes is investigated by direct numerical simulations using the complete Navier Stokes equations. The simulations are based on a spatial model so that both the base flow and the disturbance flow can develop in the downstream direction. The Navier Stokes equations are used in a vorticity-velocity form and are solved using a combination of finite difference and spectral approximations. Fourier series are used in the spanwise direction. Second-order finite-differences are used to approximate the spatial derivatives in the streamwise and transverse directions. For the temporal discretion, a combination of ADI, Crank-Nicolson, and Adams-Bashforth methods is employed. The discretized velocity equations are solved using fast Helmholtz solvers. Code validation is accomplished by comparison of the numerical results to both linear stability and to experiments. Calculations of two- and/or three-dimensional sinuous and mode disturbances in the wake of flat plate are undertaken. For calculations of two-dimensional disturbances, the wake is forced at an amplitude level so that nonlinear disturbance development may be observed. In addition, the forcing amplitude is varied in order to determine its effect on the disturbance behavior. To investigate the onset of three-dimensionality, the wake is forced with a small-amplitude three-dimensional disturbance and a larger amplitude two-dimensional disturbance. The two-dimensional forcing amplitude is varied in order to determine its influence on the three-dimensional flow field.

Fasel, Hermann F.↗

The influence of bulges on boundary-layer instability

Local disturbances caused by a spanwise surface corrugation affect the position of the boundary-layer transition, and so the drag, of an object. This premature transition from laminar to turbulent flow is often associated with a separation of the laminar boundary-layer from its surface. Also the roughness-induced separation bubble provides an important link between the pressure and velocity fluctuations in the environment and the development of the disturbance in the laminar boundary-layer, i.e., the receptivity problem. To investigate the influence of a laminar separation bubble on boundary-layer instability, a separated flow generated by a velocity gradient over a flat plate was analyzed by direct numerical simulation using finite-difference solutions of the Navier-Stokes equations. The bubble acts as a strong amplifier of the instability waves and a highly nonlinear flow field is shown to develop downstream of the bubble. Consequently, the results of the direct numerical simulation differ noticeably from those of the classical linear stability theory proving the fact that the nonparallel effects together with the nonlinear interactions are crucial to this flow development. In the present paper, the effect of physical perturbations such as humps and hollows on boundary-layer instability is analyzed. This problem has been considered theoretically by several researchers (e.g., Nayfeh et al., 1987 and 1990; Cebeci et al., 1988). They used linear stability theory in their approach which does not include the nonparallel nor the nonlinear effects. Therefore, to account for these important effects in studying flow over humps and hollows the direct simulation technique is being implemented in generalized coordinates.

Elli, S.↗

Measurements of Aerodynamic Heat Transfer and Boundary-Layer Transition on a 15 deg. Cone in Free Flight at Supersonic Mach Numbers up to 5.2

Measurements of aerodynamic heat transfer have been made at several stations on the 15 deg total-angle conical nose of a rocket-propelled model in free flight at Mach numbers up to 5.2. Data are presented for a range of local Mach number just outside the boundary layer from 1.40 to 4.65 and a range of local Reynolds number from 3.8 x 10(exp 6) to 46.5 x 10(exp 6), based on length from the nose tip to a measurement station. Laminar, transitional, and turbulent heat-transfer coefficients were measured. The laminar data were in agreement with laminar theory for cones, and the turbulent data agreed well with turbulent theory for cones using Reynolds number based on length from the nose tip. At a nearly constant ratio of wall to local static temperature of 1.2 the Reynolds number of transition increased from 14 x 10(exp 6) to 30 x 10(exp 6) as Mach number increased from 1.4 to 2.9 and then decreased to 17 x 10(exp 6) as Mach number increased to 3.7. At Mach numbers near 3.5, transition Reynolds numbers appeared to be independent of skin temperature at skin temperatures very cold with respect to adiabatic wall temperature. The transition Reynolds number was 17.7 x 10(exp 6) at a condition of Mach number and ratio of wall to local static temperature near that for which three-dimensional disturbance theory has been evaluated and has predicted laminar boundary-layer stability to very high Reynolds numbers (approximately 10(exp 12)).

Rumsey, Charles B.↗

Analysis of Fuel Vaporization, Fuel-Air Mixing, and Combustion in Integrated Mixer-Flame Holders

Requirements to limit pollutant emissions from the gas turbine engines for the future High-Speed Civil Transport (HSCT) have led to consideration of various low-emission combustor concepts. One such concept is the Integrated Mixer-Flame Holder (IMFH). This report describes a series of IMFH analyses performed with KIVA-II, a multi-dimensional CFD code for problems involving sprays, turbulence, and combustion. To meet the needs of this study, KIVA-II's boundary condition and chemistry treatments are modified. The study itself examines the relationships between fuel vaporization, fuel-air mixing, and combustion. Parameters being considered include: mixer tube diameter, mixer tube length, mixer tube geometry (converging-diverging versus straight walls), air inlet velocity, air inlet swirl angle, secondary air injection (dilution holes), fuel injection velocity, fuel injection angle, number of fuel injection ports, fuel spray cone angle, and fuel droplet size. Cases are run with and without combustion to examine the variations in fuel-air mixing and potential for flashback due to the above parameters. The degree of fuel-air mixing is judged by comparing average, minimum, and maximum fuel/air ratios at the exit of the mixer tube, while flame stability is monitored by following the location of the flame front as the solution progresses from ignition to steady state. Results indicate that fuel-air mixing can be enhanced by a variety of means, the best being a combination of air inlet swirl and a converging-diverging mixer tube geometry. With the IMFH configuration utilized in the present study, flashback becomes more common as the mixer tube diameter is increased and is instigated by disturbances associated with the dilution hole flow.

Deur, J. M.↗