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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Asymptotic Behavior of Ensemble-Averaged Linear Disturbances in Homogeneous Shear Flow

In order to expand the predictive capability of single-point turbulence closure models too account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. The calibration is based on homogeneous shear flow where disturbances can be described by rapid distort,ion theory (RDT). The relationship between RDT and linear stability theory is exploit,c d in order to obtain a closed set, of modeled equations. The linear disturbance equations are solved directly so that, the numerical simulation yields a database from which the closure coefficient,s in the ensemble-averaged disturbance equations can he determined.

Thacker, W. D.↗

NASA NTRS2000

Numerical Simulation of Receptivity for a Transition Experiment

The cost of fuel to overcome turbulence induced viscous drag on a commercial airplane constitutes a significant fraction of the operating cost of an airline. Achieving laminar flow and maintaining it over a large portion of the wing can significantly reduce the viscous drag, and hence the cost. Design of such laminar-flow-control wings and their practical operation requires the ability to accurately and reliably predict the transition from laminar to turbulent flow. The transition process begins with the conversion of environmental and surface disturbances into the instability waves of the flow by a process called receptivity. The goal of the current research project has been to improve the prediction of transition through a better understanding of the physics of receptivity. The initial objective of this work was to investigate the specific stability and receptivity characteristics of a particular experimental investigation of boundary layer receptivity at NASA Langley. Some simulation results using direct solutions of the linearized Navier-Stokes equations which modeled this experiment where presented in the 1999 APS DFD meeting. However, based on these initial investigations, it became clear that to cover the vast receptivity parameter space required for a practical transition prediction tool, more efficient methods would be required. Thus, the focus of this research was shifted from modeling this particular experiment to formulating and developing new techniques that could efficiently yet accurately predict receptivity for a wide range of disturbance conditions.

Collis, S. Scott↗

NASA NTRS2020

Verification of Anisotropic Mesh Adaptation for Turbulent Simulations over ONERA M6 Wing

Unstructured anisotropic mesh adaptation is known to be an efficient way to control discretization errors in Computational Fluid Dynamics (CFD) simulations. Method verification is required to provide the confidence for routine use in production analysis. The current work aims at verification of anisotropic mesh adaptation for RANS simulations over the ONERA M6 wing. The present verification study is performed using four different flow solvers, three different implementations of the metric field, and three mesh mechanics packages. Two of the flow solvers use stabilized finite-element discretizations (FUN3D-SFE and GGNS), one uses finite-volume discretization (FUN3D-FV), and the last one uses mixed finite-volume and finite element discretizations (Wolf). The mesh adaptation is based on an error estimator that aims to control the quadratic error term in the linear interpolation of Mach number. Two sets of adaptations were performed; the first one controls the interpolation error in L2 norm and the second one controls the interpolation error in L4 norm. Convergence studies were performed on the forces and the pitching moment using all four solvers, and the results are compared with previously verified convergence studies on fixed (nonadapted) meshes. Both forces and pitching moment on adapted meshes are found to be converging to the fine mesh values faster than those on fixed meshes. In addition to forces and moments, convergence of surface pressure and skin friction coefficients at various measurement locations on the wing are also presented. Adapted-mesh surface pressure distributions agree with the fine fixed mesh pressure distributions. Adapted-mesh skin friction distributions contain high frequency noise with mean values approaching the fixed mesh pressure skin friction distributions.

Aravind Balan↗

NASA NTRS1978

A scheme for computing surface layer turbulent fluxes from mean flow surface observations

A physical model and computational scheme are developed for generating turbulent surface stress, sensible heat flux and humidity flux from mean velocity, temperature and humidity at some fixed height in the atmospheric surface layer, where conditions at this reference level are presumed known from observations or the evolving state of a numerical atmospheric circulation model. The method is based on coupling the Monin-Obukov surface layer similarity profiles which include buoyant stability effects on mean velocity, temperature and humidity to a force-restore formulation for the evolution of surface soil temperature to yield the local values of shear stress, heat flux and surface temperature. A self-contained formulation is presented including parameterizations for solar and infrared radiant fluxes at the surface. Additional parameters needed to implement the scheme are the thermal heat capacity of the soil per unit surface area, surface aerodynamic roughness, latitude, solar declination, surface albedo, surface emissivity and atmospheric transmissivity to solar radiation.

Hoffert, M. I.↗

NASA NTRS1990

Modern developments in shear flow control with swirl

Passive and active control of swirling turbulent jets is experimentally investigated. Initial swirl distribution is shown to dominate the free jet evolution in the passive mode. Vortex breakdown, a manifestation of high intensity swirl, was achieved at below critical swirl number (S = 0.48) by reducing the vortex core diameter. The response of a swirling turbulent jet to single frequency, plane wave acoustic excitation was shown to depend strongly on the swirl number, excitation Strouhal number, amplitude of the excitation wave, and core turbulence in a low speed cold jet. A 10 percent reduction of the mean centerline velocity at x/D = 9.0 (and a corresponding increase in the shear layer momentum thickness) was achieved by large amplitude internal plane wave acoustic excitation. Helical instability waves of negative azimuthal wave numbers exhibit larger amplification rates than the plane waves in swirling free jets, according to hydrodynamic stability theory. Consequently, an active swirling shear layer control is proposed to include the generation of helical instability waves of arbitrary helicity and the promotion of modal interaction, through multifrequency forcing.

Farokhi, Saeed↗

NASA NTRS1990

Modern developments in shear flow control with swirl

Farokhi, S.↗

NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the nonlocal, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, such as those in the vicinity of roughness elements, steps, gaps, or corners. To compute the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 10^(7) to 10^(8) even in relatively simple flows. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear systems. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE. Additionally, inflow-resolvent analysis based on the HLNSE for discrete-roughness-induced disturbances in the nose tip of a blunt cone at Mach 6 demonstrates the importance of including the disturbance amplification along the near vicinity of the roughness element and separation region.

Boundary Layer Stability↗

NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗

NASA NTRS1990

The application of active controls technology to a generic hypersonic aircraft configuration

Analytical methods are described for the prediction of aerothermoelastic stability of hypersonic aircraft including active control systems. Thermal loads due to aerodynamic heating were applied to the finite element model of the aircraft structure and the thermal effects on flutter were determined. An iterative static aeroelastic trim analysis procedure was developed including thermal effects. And active control technology was assessed for flutter suppression, ride quality improvement, and gust load alleviation to overcome any potential adverse aeroelastic stability or response problems due to aerodynamic heating. A generic hypersonic aircraft configuration was selected which incorporates wing flaps, ailerons, and all moveable fins to be used for active control purposes. The active control system would use onboard sensors in a feedback loop through the aircraft flight control computers to move the surfaces for improved structural dynamic response as the aircraft encounters atmospheric turbulence.

Gilbert, M. G.↗

NASA NTRS2001

Predicting Turbulent Convective Heat Transfer in Fully Developed Duct Flows

The performance of an explicit algebraic stress model (EASM) is assessed in predicting the turbulent flow and forced heat transfer in both straight and wavy ducts, with rectangular, trapezoidal and triangular cross-sections, under fully developed conditions. A comparison of secondary flow patterns. including velocity vectors and velocity and temperature contours, are shown in order to study the effect of waviness on flow dynamics, and comparisons between the hydraulic parameters. Fanning friction factor and Nusselt number, are also presented. In all cases. isothermal conditions are imposed on the duct walls, and the turbulent heat fluxes are modeled using gradient-diffusion type models. The formulation is valid for Reynolds numbers up to 10(exp 5) and this minimizes the need for wall functions that have been used with mixed success in previous studies of complex duct flows. In addition, the present formulation imposes minimal demand on the number of grid points without any convergence or stability problems. Criteria in terms of heat transfer and friction factor needed to choose the optimal wavy duct cross-section for industrial applications among the ones considered are discussed.

Rokni, Masoud↗

NASA NTRS2001

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

Sjoegreen, B.↗

NASA NTRS1972

Unsteady viscous effects in the flow over an oscillating surface

A theoretical model for the interaction of a turbulent boundary layer with an oscillating wavy surface over which a fluid is flowing is developed, with an application to wind-driven water waves and to panel flutter in low supersonic flow. A systematic methodology is developed to obtain the surface pressure distribution by considering separately the effects on the perturbed flow of a mean shear velocity profile, viscous stresses, the turbulent Reynolds stresses, compressibility, and three-dimensionality. The inviscid theory is applied to the wind-water wave problem by specializing to traveling-wave disturbances, and the pressure magnitude and phase shift as a function of the wave phase speed are computed for a logarithmic mean velocity profile and compared with inviscid theory and experiment. The results agree with experimental evidence for the stabilization of the panel motion due to the influence of the unsteady boundary layer.

Lerner, J. I.↗

NASA NTRS1984

Analysis for leakage and rotordynamic coefficients of surface roughened tapered annular gas seals

In order to soften the effects of rub, the smooth stators of turbine gas seals are sometimes replaced by a honeycomb surface. This deliberately roughened stator and smooth rotor combination retards the seal leakage and may lead to enhanced rotor stability. However, many factors determine the rotordynamic coefficients and little is known as to the effectiveness of these honeycomb seals under various changes in the independent seal parameters. An analytical-computational method to solve for the rotordynamic coefficients of this type of compressible-flow seal is developed. The governing equations for surface roughned tapered annular gas seals are based on a modified Hirs' turbulent bulk flow model. A perturbation analysis is employed to develop zeroth and first-order perturbation equations. These equations are numerically integrated to solve for the leakage, pressure, density, and velocity for small motion of the shaft about the centered position. The resulting pressure distribution is then integrated to find the corresponding rotor-dynamic coefficients.

Nelson, C. C.↗

NASA NTRS1993

Numerical simulation of laminar-turbulent transition in a spatially-developing flat plate wake

Laminar-turbulent transition of an incompressible flat-plate wake is investigated by direct numerical integration of the Navier-Stokes equations. For the numerical integration, a combination of finite-difference and spectral methods along with an ADI/Crank-Nicolson/Adams-Bashforth time integration scheme is employed. Subject to 2D forcing, the wake exhibited a rapidly-growing fundamental disturbance that quickly saturated. This saturation was due partly to the stabilizing effect of the mean flow distortion. Downstream of the saturation point, disturbance energy was concentrated in the fundamental disturbance, the second harmonic, and the mean flow distortion component. At large amplitude levels, a Karman vortex street formed. Variations in the 2D forcing level did not alter the qualitative behavior of the disturbances. Simulations of 3D breakdown indicates that the presence of large-amplitude, 2D disturbances tends to initially suppress small-amplitude 3D disturbance growth. Following this initial suppression, a resumption of 3D growth is observed that may have been due to a secondary instability mechanism. For high levels of 3D disturbance energy, lambda-shaped vortical structures formed between adjacent Karman vortices.

Dratler, D. I.↗

NASA NTRS1993

The Structure and Stability of Laminar Flames

This review paper on the structure and stability of laminar flames considers such phenomena as heterogeneous mixtures, acoustic instabilities, flame balls and related phenomena, radiation effects, the iodate oxidation of arsenous acid and 'liquid flame fronts', approximate kinetic mechanisms and asymptotic approximations, and tribrachial or triple flames. The topics examined here indicate three themes that may play an important role in laminar flame theory in the coming years: microgravity experiments, kinetic modeling, and turbulence modeling. In the discussion of microgravity experiments it is pointed out that access to drop towers, the Space Shuttle and, in due course, the Space Station Freedom will encourage the development of experiments well designed to isolate the fundamental physics of combustion.

Buckmaster, John↗

NASA NTRS2002

Simultaneous Raman-Rayleigh-LIF Measurements and Numerical Modeling Results of a Lifted H2/N2 Turbulent Jet Flame in a Vitiated Coflow

An experimental and numerical investigation is presented of a H2/N2 turbulent jet flame burner that has a novel vitiated coflow. 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 for exploration of recirculation chemistry without the corresponding fluid mechanics of recirculation. Thus the vitiated coflow burner design facilitates the development of chemical kinetic combustion models without the added complexity of recirculation fluid mechanics. Scalar measurements are reported for a turbulent jet flame of H2/N2 in a coflow of combustion products from a lean ((empty set) = 0.25) H2/Air flame. The combination of laser-induced fluorescence, Rayleigh scattering, and Raman scattering is used to obtain simultaneous measurements of the temperature, major species, as well as OH and NO. Laminar flame calculation with equal diffusivity do agree when the premixing and preheating that occurs prior to flame stabilization is accounted for in the boundary conditions. Also presented is an exploratory pdf model that predicts the flame's axial profiles fairly well, but does not accurately predict the lift-off height.

Cabra, R.↗

NASA NTRS2011

ARTEMIS Mission Overview: From Concept to Operations

ARTEMIS (Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun) repurposed two spacecraft to extend their useful science (Angelopoulos, 2010) by moving them via lunar gravity assists from elliptical Earth orbits to L1 and L2 Earth-Moon libration orbits and then to lunar orbits by exploiting the Earth-Moon-Sun dynamical environment. This paper describes the complete design from conceptual plans using weak stability transfer options and lunar gravity assist to the implementation and operational support of the Earth-Moon libration and lunar orbits. The two spacecraft of the ARTEMIS mission will have just entered lunar orbit at this paper's presentation.

Folta, David↗

NASA NTRS2011

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 l.2xl2 Earth radii in the first phase and l.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.↗