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 757 records · Page 42

Assessment of RANS-based Transition Models based on Experimental Data of the Common Research Model with Natural Laminar Flow

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework often rely upon transition correlations that were derived from a limited number of low-speed experiments and these models often fail to account for all of the relevant transition mechanisms and/or the variation in those mechanisms with respect to changes in the significant flow parameters. Available data from a recent experiment on the Common Research Model with Natural Laminar Flow (CRM-NLF) in the National Transonic Facility at the NASA Langley Research Center are used to assess the current transition modeling capability in NASA's OVERFLOW 2.3b solver for a swept wing configuration with nonzero taper and transonic cruise conditions. Specifically, the OVERFLOW solutions are used to evaluate the accuracy and robustness of the transport-equation-based transition models. Results highlight that the Spalart-Allmaras-based amplification factor transport (AFT-2017b) equation model and Menter’s shear-stress transport equation (SST2003)-based Langtry-Menter transition models (either with or without the modeling of crossflow transition) significantly underpredict the reported extent of laminar flow region over the entire span of the wing, irrespective of which instability mechanism(s) is expected to dominate the onset of the transition process. We show that the transition correlations underlying these models fail to account for the stabilizing effect of compressibility on the Tollmien-Schlichting transition, which is likely to be a major contributor to the underprediction of the laminar flow region on the CRM-NLF. The SST-2003-based Langtry-Menter model also appears to inaccurately predict the chordwise pressure variation along the majority of the wing span at all the flow conditions studied herein, due to how the turbulence intensity levels were enforced in the computations and how that was interfering with the functioning of the underlying turbulence model within the boundary layer. The AFT and Langtry-Menter models appear to be sensitive to the level of the freestream turbulence intensity, but the degree of sensitivity varies across the models.

CFD modeling↗

A numerical method for the prediction of high-speed boundary-layer transition using linear theory

A method is described of estimating the location of transition in an arbitrary laminar boundary layer on the basis of linear stability theory. After an examination of experimental evidence for the relation between linear stability theory and transition, a discussion is given of the three essential elements of a transition calculation: (1) the interaction of the external disturbances with the boundary layer; (2) the growth of the disturbances in the boundary layer; and (3) a transition criterion. The computer program which carried out these three calculations is described. The program is first tested by calculating the effect of free-stream turbulence on the transition of the Blasius boundary layer, and is then applied to the problem of transition in a supersonic wind tunnel. The effects of unit Reynolds number and Mach number on the transition of an insulated flat-plate boundary layer are calculated on the basis of experimental data on the intensity and spectrum of free-stream disturbances. Reasonable agreement with experiment is obtained in the Mach number range from 2 to 4.5.

Mack, L. M.↗

The Venus SAGE Atmospheric Structure Investigation

Experiment goals and objectives are: a) To accurately define the state properties as a function of altitude from below the 10(exp -4) mb level (approx.150 km) to 92 bars (surface); b) To measure the stability of the atmosphere, and identify convective layers and stable layers, where they exist; c) To detect cloud levels from changes in the lapse rate at their boundaries; d) To provide state properties within the cloud levels, and thus provide supplementary information on cloud composition; e) To search for and characterize wave structure within the atmosphere; f) To search for and measure the intensity and scale of turbulence; g) To measure descent and surface wind speed and direction; h) To provide Lander altitude and attitude during decent for descent imaging analysis; and i) To provide a back-up landing sensor.

Colaprete, Anthony↗

Ensemble Diffraction Measurements of Spray Combustion in a Novel Vitiated Coflow Turbulent Jet Flame Burner

An experimental investigation is presented of a novel vitiated coflow spray flame burner. The vitiated coflow emulates the recirculation region of most combustors, such as gas turbines or furnaces; additionally, since the vitiated gases are coflowing, the burner allows exploration of the chemistry of recirculation without the corresponding fluid mechanics of recirculation. As such, this burner allows for chemical kinetic model development without obscurations caused by fluid mechanics. The burner consists of a central fuel jet (droplet or gaseous) surrounded by the oxygen rich combustion products of a lean premixed flame that is stabilized on a perforated, brass plate. The design presented allows for the reacting coflow to span a large range of temperatures and oxygen concentrations. Several experiments measuring the relationships between mixture stoichiometry and flame temperature are used to map out the operating ranges of the coflow burner. These include temperatures as low 300 C to stoichiometric and oxygen concentrations from 18 percent to zero. This is achieved by stabilizing hydrogen-air premixed flames on a perforated plate. Furthermore, all of the CO2 generated is from the jet combustion. Thus, a probe sample of NO(sub X) and CO2 yields uniquely an emission index, as is commonly done in gas turbine engine exhaust research. The ability to adjust the oxygen content of the coflow allows us to steadily increase the coflow temperature surrounding the jet. At some temperature, the jet ignites far downstream from the injector tube. Further increases in the coflow temperature results in autoignition occurring closer to the nozzle. Examples are given of methane jetting into a coflow that is lean, stoichiometric, and even rich. Furthermore, an air jet with a rich coflow produced a normal looking flame that is actually 'inverted' (air on the inside, surrounded by fuel). In the special case of spray injection, we demonstrate the efficacy of this novel burner with a methanol spray in a vitiated coflow. As a proof of concept, an ensemble light diffraction (ELD) optical instrument was used to conduct preliminary measurements of droplet size distribution and liquid volume fraction.

Cabra, R.↗

Dynamics of face and annular seals with two-phase flow

A detailed study was made of face and annular seals under conditions where boiling, i.e., phase change of the leaking fluid, occurs within the seal. Many seals operate in this mode because of flashing due to pressure drop and/or heat input from frictional heating. High pressure, water pumps, industrial chemical pumps, and cryogenic pumps are mentioned as a few of many applications. The initial motivation was the LOX-GOX seals for the space shuttle main engine, but the study was expanded to include any face or annular seal where boiling occurs. Some of the distinctive behavior characteristics of two-phase seals were discussed, particularly their axial stability. While two-phase seals probably exhibit instability to disturbances of other degrees of freedom such as wobble, etc., under certain conditions, such analyses are too complex to be treated at present. Since an all liquid seal (with parallel faces) has a neutral axial stiffness curve, and is stabilized axially by convergent coning, other degrees of freedom stability analyses are necessary. However, the axial stability behavior of the two-phase seal is always a consideration no matter how well the seal is aligned and regardless of the speed. Hence, axial stability is thought of as the primary design consideration for two-phase seals and indeed the stability behavior under sub-cooling variations probably overshadows other concerns. The main thrust was the dynamic analysis of axial motion of two-phase face seals, principally the determination of axial stiffness, and the steady behavior of two-phase annular seals. The main conclusions are that seals with two-phase flow may be unstable if improperly balanced. Detailed theoretical analyses of low (laminar) and high (turbulent) leakage seals are presented along with computer codes, parametric studies, and in particular a simplified PC based code that allows for rapid performance prediction. A simplified combined computer code for the performance prediction over the laminar and turbulent ranges of a two-phase seal is described and documented. The analyses, results, and computer codes are summarized.

Hughes, William F.↗

Rotordynamic Instability Problems in High-Performance Turbomachinery 1996

The first rotordynamics workshop proceedings emphasized a feeling of uncertainty in predicting the stability of characteristics of high-performance turbomachinery. In the second workshop proceedings these uncertainties were reduced through programs established to systematically resolve problems, with emphasis on experimental validation of the forces that influence rotordynamics. In the third proceedings many programs for predicting or measuring forces and force coefficients in high-performance turbomachinery produced results. Data became available for designing new machines with enhanced stability characteristics or for upgrading existing machines. In the fourth proceedings there emerged trends towards a more unified view of rotordynamic instability problems and several encouraging new analytical developments. The fifth workshop supported the continuing trend toward a unified view with several new developments in the design and manufacture of new turbomachineries with enhanced stability characteristics along with new data and associated numerical/theoretical results. The sixth workshop report provided field experience and experimental results, and expanded the use of computational and control techniques with integration of damper, bearing, and eccentric seal operation results. The seventh workshop report provided field experiences, numerical, theoretical, and experimental results and control methods for seals, bearings, and dampers with some attention given to variable thermophysical properties and turbulence measurements, and introduction of two-phase flow results. In the present workshop, active magnetic bearings (AMB's) evolve into a new method of measuring rotordynamic coefficients with discussions on honeycomb seals, drop of magnetically supported rotors, seals, bearings and dampers with new data being reported. The intent of the workshop and this proceedings is to provide a continuing impetus for an understanding and resolution of these problems.

Source record↗

Subcritical transition to turbulence in planar shear flows

The two-dimensional steady and time dependent properties of plane Poiseuille and plane Couette flows are analyzed using iterative techniques and full numerical simulation of the Navier-Stokes equations. It is shown that the finite-amplitude two-dimensional states investigated are strongly unstable to very small three-dimensional perturbations. It is also shown, through full numerical simulation, that this explosive secondary instability can explain the subcritical transitions that occur in real flows. Finally, it is shown that the three-dimensional instability can be analyzed by a linear stability analysis of a two-dimensional flow consisting of the basic parallel flow and a steady (or quasi-steady) finite-amplitude two-dimensional cellular motion.

Orszag, S. A.↗

Studies of heat source driven natural convection

Natural convection energy transport in a horizontal layer of internally heated fluid with a zero heat flux lower boundary, and an isothermal upper boundary, has been studied. Quantitative information on the time-mean temperature distribution and the fluctuating component of temperature about the mean temperature in steady turbulent convection are obtained from a small thermocouple inserted into the layer through the upper bounding plate. Data are also presented on the development of temperature at several vertical positions when the layer is subject to both a sudden increase and to a sudden decrease in power input. For changes of power input from zero to a value corresponding to a Rayleigh number much greater than the critical linear stability theory value, a slight hysteresis in temperature profiles near the upper boundary is observed between the heat-up and cool-down modes.

Kulacki, F. A.↗

On the theory of solitary Rossby waves

The evolution of long, finite amplitude Rossby waves in a horizontally sheared zonal current is studied. The wave evolution is described by the Korteweg-de Vries equation or the modified Korteweg-de Vries equation depending on the atmospheric stratification. In either case, the cross-stream modal structure of these waves is given by the long-wave limit of the neutral eigensolutions of the barotropic stability equation. Both non-singular and singular eigensolutions are considered and the appropriate analysis is developed to yield a uniformly valid description of the motion in the critical-layer region where the wave speed matches the flow velocity. The analysis demonstrates that coherent, propagating, eddy structures can exist in stable shear flows and that these eddies have peculiar interaction properties quite distinct from the traditional views of turbulent motion.

Redekopp, L. G.↗

Magnetospheric Multiscale (MMS) Mission Attitude Ground System Design

This paper describes the attitude ground system (AGS) design to be used for support of the Magnetospheric MultiScale (MMS) mission. The AGS exists as one component of the mission operations control center. It has responsibility for validating the onboard attitude and accelerometer bias estimates, calibrating the attitude sensors and the spacecraft inertia tensor, and generating a definitive attitude history for use by the science teams. NASA's Goddard Space Flight Center (GSFC) in Greenbelt, Maryland is responsible for developing the MMS spacecraft, for the overall management of the MMS mission, and for mission operations. MMS is scheduled for launch in 2014 for a planned two-year mission. The MMS mission consists of four identical spacecraft flying in a tetrahedral formation in an eccentric Earth orbit. The relatively tight formation, ranging from 10 to 400 km, will provide coordinated observations giving insight into small-scale magnetic field reconnection processes. By varying the size of the tetrahedron and the orbital semi-major axis and eccentricity, and making use of the changing solar phase, this geometry allows for the study of both bow shock and magnetotail plasma physics, including acceleration, reconnection, and turbulence. The mission divides into two phases for science; these phases will have orbit dimensions of 1.2 x 12 Earth radii in the first phase and 1.2x25 Earth radii in the second in order to study the dayside magnetopause and the nightside magnetotail, respectively. The orbital periods are roughly one day and three days for the two mission phases. Each of the four MMS spacecraft will be spin stabilized at 3 revolutions per minute (rpm), with the spin axis oriented near the ecliptic north pole but tipped approximately 2.5 deg towards the Sun line. The main body of each spacecraft will be an eight-sided platform with diameter of 3.4 m and height of 1.2 m. Several booms are attached to this central core: two axial booms of 14.9 m length, two radial magnetometer booms of 5 m length, and four radial wire booms of 60 m length. Attitude and orbit control will use a set of axial and radial thrusters. A four-head star tracker and a slit-type digital Sun sensor (DSS) provide input for attitude determination. In addition, an accelerometer will be used for closed-loop orbit maneuver control. The primary AGS product will be a daily definitive attitude history. Due to power limitations, the star tracker and accelerometer data will not be available at all times. However, tracker data from at least 10 percent of each orbit and continuous DSS data will be provided. An extended Kalman filter (EKF) will be used to estimate the three-axis attitude (i.e., spin axis orientation and spin phase) and rotation rate for all times when the tracker data is valid. For other times, the attitude is generated by assuming a constant angular momentum vector in the inertial frame. The DSS sun pulse will provide a timing signal to maintain an accurate spin phase. There will be times when the Sun is occulted and DSS data is not available. If this occurs at the start or end of a definitive attitude product, then the spin phase will be extrapolated using the mean rate determined by the EKF.

Sedlak, Joseph E.↗

Observations of mesoscale vertical velocities around frontal zones

Vertical velocity and reflectivity data obtained with a VHF Doppler radar over a 15-day period in October and November of 1981 are analyzed. Standard radiosonde data and surface observations were used to locate two occluded fronts, two warm fronts, and a cold front that passed the radar site. These fronts are also evident in the radar reflectivity data. Most studies of the vertical circulation patterns associated with mososcale systems have used precipitation and cloud formations as tracers. Unlike other observational techniques, the VHF radar permits the continuous measurement of the three-dimensional air velocity vector in time and height from a fixed location. With the beam in a vertically pointing position, signals are scattered from turbulent variations in the refractive index with half the scale of the radar wavelength and by regions with sudden changes in the refractive index associated with horizontally stratified layers. Generally, the strongest echoes occur at the maximum in the vertical gradient of refractivity, usually at the base of a temperature inversion, such as the tropopause. VHF radars can also be used to locate atmospheric fronts, which are characterized by static stability, large horizontal temperature gradients, large vorticities, and vertical wind shears. These radars can provide the velocity field data needed to study wave motions associated with fronts and compare the actual vertical circulation to theoretical predictions.

Dennis, T. S.↗

On recent advances and future research directions for computational fluid dynamics

This paper highlights some recent accomplishments regarding CFD numerical algorithm constructions for generation of discrete approximate solutions to classes of Reynolds-averaged Navier-Stokes equations. Following an overview of turbulent closure modeling, and development of appropriate conservation law systems, a Taylor weak-statement semi-discrete approximate solution algorithm is developed. Various forms for completion to the final linear algebra statement are cited, as are a range of candidate numerical linear algebra solution procedures. This development sequence emphasizes the key building blocks of a CFD RNS algorithm, including solution trial and test spaces, integration procedure and added numerical stability mechanisms. A range of numerical results are discussed focusing on key topics guiding future research directions.

Baker, A. J.↗

Multigrid Acceleration of Time-Accurate Navier-Stokes Calculations

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is fully implicit in time and is unconditionally stable (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability. This is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. A multigrid-multiblock, steady-state, three-dimensional Navier-Stokes solver, TLNS3D, was modified to iteratively invert the equations at each physical time step. The implementation of this procedure in TLNS3D is discussed. The implications of applying several popular turbulence models to unsteady flow are also considered. Numerical results are presented to show the application of the scheme to various two-dimensional turbulent flows. The results of a three-dimensional laminar flow calculation are also given.

Melson, N. Duane↗

High-Order Simulation of Induced Disturbance in a Mach 6 Boundary Layer

Simulations of the early evolution of induced unstable disturbances in a Mach 6 hypersonic boundary layer are presented for the purposes of validating a high-order discontinuous Galerkin Navier-Stokes simulation code. The simulations are modeled after experiments performed in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. In these experiments, a forcing mechanism was employed to induce reproducible disturbances in a hypersonic boundary layer providing for the controlled study of the growth and breakdown of these disturbances into turbulent spots. Simulations revealed that the form and strength of the excitation can greatly influence the growth of the disturbance. In particular, at large forcing amplitude, the simulated forcing produces large advecting transients that appear to enhance the growth of the wave packet, relative to that of low amplitude forcing. Simulation results with large amplitude forcing agree well with experimental results while results from low amplitude forcing agree with linear stability theory.

Atkins, Harold L.↗

Transition Analysis for Isolated Trips on BOLT-II Wind-Tunnel and Flight Configuration

The BOLT-II configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. The present paper is related to the NASA-led roughness experiment on the secondary side of the BOLT-II configuration. First, computations are performed to aid mean-flow based correlations for the transition measurements acquired during an experimental campaign in the Ludwieg Tube facility at the United States Air Force Academy. To that end, both laminar and turbulent flow computations are performed for selected conditions of the experiment with a 1:3 scale model of the slightly shorter BOLT configuration. The variation in well-known correlation parameters for transition due to isolated roughness elements with respect to element position, height, and flow Reynolds number is delineated. Furthermore, the wake instability theory is used to examine the amplification of unstable disturbances behind the centerline trip at the design conditions of the BOLT-II flight experiment and for its dynamically scaled version on the 1:3 scale model of the wind tunnel experiment. To our knowledge, this study represents the first analysis of the wake instabilities in the presence of a prominent entropy layer outside the boundary layer. At the nominal flight design condition of Re ∞ ) = 5.44 and Re ∞ = 2.5 x 10 6 /m, a diamond planform (“pizza-box”) trip with k/δ = 0.70 and planform-halfwidth-to-height ratio of b/k = 3.0 at X/L = 0.254 induces a strong distortion of the boundary-layer flow in the vicinity of the symmetry plane, such that the wake structure includes a centerline streak produced by the flow retardation behind the trip along with multiple additional, prominent vorticity structures representing the necklace vortex system. Stability results reveal that each of these streaks can support a strong amplification of wake instabilities within a frequency band of 95 kHz–255 kHz, with a peak N-factor of greater than 15 beyond X/L ≈ 0.42. The dynamically scaled wind tunnel configuration with Re ∞ = 20.25 x 10 6 /m is also found to support N > 15 at similar axial stations.

BOLT↗

Prediction of an axisymmetric combusting flow

A numerical model for turbulent, recirculating combusting flow is developed and applied to a research combustor. The model is based on the time-averaged Navier-Stokes equations with k-epsilon turbulence closure and the compositional fluctuations at each point are given probabilistically in terms of the mixture fraction. The probability density function is derived from transport equations for its first two moments along with an assumption regarding its shape. The resulting equations are solved using a standard line relaxation algorithm. It is found that the predictions of the model are in good agreement with data from an axisymmetric, bluff-body stabilized research combustor, while the major discrepancies are similar to those found in isothermal flow comparisons. The peculiar features of this flow which contribute to the errors are examined. The agreement between the theory and data deteriorates as the central jet velocity is increased which indicates an enhanced role for unsteady effects.

Correa, S. M.↗

Propagation of high-energy laser beams through the earth's atmosphere II; Proceedings of the Meeting, Los Angeles, CA, Jan. 21-23, 1991

Consideration is given to turbulence at the inner scale, modeling turbulent transport in laser beam propagation, variable wind direction effects on thermal blooming correction, realistic wind effects on turbulence and thermal blooming compensation, wide bandwidth spectral measurements of atmospheric tilt turbulence, remote alignment of adaptive optical systems with far-field optimization, focusing infrared laser beams on targets in space without using adaptive optics, and a simplex optimization method for adaptive optics system alignment. Consideration is also given to ground-to-space multiline propagation at 1.3 micron, a path integral approach to thermal blooming, functional reconstruction predictions of uplink whole beam Strehl ratios in the presence of thermal blooming, and stability analysis of semidiscrete schemes for thermal blooming computation.

Ulrich, Peter B.↗

Numerical simulation of two-dimensional spatially-developing mixing layers

Two-dimensional, incompressible, spatially developing mixing layer simulations are performed at Re = 10(exp 2) and 10(exp 4) with two classes of perturbations applied at the inlet boundary; combinations of discrete modes from linear stability theory, and a broad spectrum of modes derived from experimentally measured velocity spectra. The effect of the type and strength of inlet perturbations on vortex dynamics and time-averaged properties are explored. Two-point spatial velocity and autocorrelations are used to estimate the size and lifetime of the resulting coherent structures and to explore possible feedback effects. The computed time-averaged properties such as mean velocity profiles, turbulent statistics, and spread rates show good agreement with experimentally measured values. It is shown that by forcing with a broad spectrum of modes derived from an experimental energy spectrum many experimentally observed phenomena can be reproduced by a 2-D simulation. The strength of the forcing merely affected the length required for the dominant coherent structures to become fully-developed. Thus intensities comparable to those of the background turbulence in many wind tunnel experiments produced the same results, given sufficient simulation length.

Wilson, R. V.↗