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 613 records · Page 34

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbomachinery

A rapid quasi three-dimensional analysis was developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations which are solved b a two-stage Runge-Kutta scheme made efficient by use of vectorization, a variable time-step, and a flux-based multigrid scheme, are described. A stability analysis is presented for the two-stage. Results for a flat-plate model problem show the applicability of the method to axial, radial, and rotating geometries. Results for a centrifugal impeller and a radial diffuser show that the quasi three-dimensional viscous analysis can be a practical design tool.

Chima, R. V.↗

Convective heat flow in space cryogenics plugs - Critical and moderate He II heat flux densities

Plug flow rates of entropy, heat and normal fluid in phase separators and in zero net mass flow systems are, to some extent, quite similar. A simplified analysis of critical conditions is presented in agreement with data trends. A critical temperature gradient arises on the basis of the He II two-fluid model at the stability limit constraining the thermohydrodynamics of the system. Thus, the question of critical thermodynamic fluctuations associated with nucleation versus the possibility of critical gradients in externally imposed parameters is answered in favor of the latter route toward turbulence. Furthermore, a similarity equation is presented which incorporates size dependent rates for moderate heat flow densities observed in experiments.

Yuan, S. W. K.↗

Does ion tearing exist?

Theoretical models of reconnection in the geomagnetic tail and its relationship to the collisionless tearing mode are examined analytically. The focus is on attempts to remove the stabilizing effects of electron compressibility in the quasi-neutral sheet. It is shown that this goal can be achieved by introducing spatial diffusion of electrons across the magnetic-flux surfaces, but not by wave-turbulent pitch-angle scattering (Coroniti, 1980) or nonadiabatic stochastic first-invariant diffusion (Buechner and Zelenyi, 1988). Hence the ion tearing mode does not exist, but the electron tearing mode does.

Pellat, R.↗

Alleviation of side force on tangent-ogive forebodies using passive porosity

An experimental investigation to determine the effectiveness of porosity for alleviating side forces on forebodies was conducted in the NASA Langley Research Center 7 x 10 ft High-Speed Wind Tunnel. The study consisted of a comparison of experimental force, moment, and surface pressure results obtained on a fineness ratio 5.0 tangent-ogive porous forebody with 0.020 in. hole diameter and 22 percent porosity with results obtained on a solid forebody. The forebodies were tested with cylindrical afterbodies. The solid forebody was tested with transition grit to simulate fully turbulent conditions and without transition grit to simulate free transition conditions. The extent of porosity on the forebody was varied to determine the extent of porosity needed to alleviate side forces. Static longitudinal and lateral-directional stability and surface pressure data were obtained at Mach numbers of 0.2, 0.5, and 0.8. The angle of attack range was from 5 to 45 deg and roll angles from -90 to 180 deg. The solid forebody exhibited large asymmetries at moderate to high angles of attack causing large side forces and yawing moments. The porous forebody exhibited no significant side forces or yawing moments at any angle of attack tested.

Bauer, Steven X. S.↗

Computational fluid mechanics

Two papers are included in this progress report. In the first, the compressible Navier-Stokes equations have been used to compute leading edge receptivity of boundary layers over parabolic cylinders. Natural receptivity at the leading edge was simulated and Tollmien-Schlichting waves were observed to develop in response to an acoustic disturbance, applied through the farfield boundary conditions. To facilitate comparison with previous work, all computations were carried out at a free stream Mach number of 0.3. The spatial and temporal behavior of the flowfields are calculated through the use of finite volume algorithms and Runge-Kutta integration. The results are dominated by strong decay of the Tollmien-Schlichting wave due to the presence of the mean flow favorable pressure gradient. The effects of numerical dissipation, forcing frequency, and nose radius are studied. The Strouhal number is shown to have the greatest effect on the unsteady results. In the second paper, a transition model for low-speed flows, previously developed by Young et al., which incorporates first-mode (Tollmien-Schlichting) disturbance information from linear stability theory has been extended to high-speed flow by incorporating the effects of second mode disturbances. The transition model is incorporated into a Reynolds-averaged Navier-Stokes solver with a one-equation turbulence model. Results using a variable turbulent Prandtl number approach demonstrate that the current model accurately reproduces available experimental data for first and second-mode dominated transitional flows. The performance of the present model shows significant improvement over previous transition modeling attempts.

Hassan, H. A.↗

Numerical simulation of supersonic boundary layer transition

The present contribution reviews some of the recent progress obtained at our group in the direct numerical simulation (DNS) of compressible boundary layer transition. Elements of the different simulation approaches and numerical techniques employed are surveyed. Temporal and spatial simulations, as well as comparisons with results obtained from Parabolized Stability Equations, are discussed. DNS results are given for flat plate boundary layers in the Mach number range 1.6 to 4.5. A temporal DNS at Mach 4.5 has been continued through breakdown all the way to the turbulent stage. In addition results obtained with a recently developed extended temporal DNS approach are presented, which takes into account some nonparallel effects of a growing boundary layer. Results from this approach are quite close to those of spatial DNS, while preserving the efficiency of the temporal DNS.

Guo, Y.↗

Novel Epoxy Particulate Composites for Mitigation of Insect Residue Adhesion on Future Aircraft Surfaces

Drag is reduced significantly for airflow over surfaces when laminar flow can be maintained over greater chord lengths, the distance from the leading edge of an airfoil.1 However, surface imperfections, such as chipped paint, scratches, and events that change topography on a microscopic scale can introduce airflow instabilities resulting in premature transition to turbulent flow.1 Although many of these surface imperfections can be avoided with proper maintenance, advanced materials, and advanced manufacturing practices, topographical surface anomalies arising during flight from insect impacts cannot be controlled and can influence laminar flow stability. Practical solutions to this operational challenge need to be developed for future aircraft to have full advantage of laminar flow designs that improve fuel efficiency.2 Researchers have investigated various methods to mitigate insect residue adhesion for decades.3 Although several techniques have demonstrated efficacy including mechanical scrapers, active liquid discharge systems, and sacrificial paper coatings, they have not been commercially implemented due to increased manufacturing and operational complexity, environmental impact, and weight penalties. Coatings offer a simple route for passive insect residue adhesion prevention without many of the challenges associated with maintenance of laminar flow.4 In our previous work, we determined that most commercially available materials were not effective at insect residue adhesion.5 We also identified improvements when both surface energy could be controlled by surface modifying agents and the topography could be altered through the use of micron-sized and nanometer-sized filler materials.6 In this work, these general principles were applied to an epoxy system to evaluate the behavior of the surface modifying agent, a fluorinated alkyl ether oligomer, on surface energy and insect residue adhesion properties.

Wohl, Christopher J.↗

Beam-plasma instability in the presence of low-frequency turbulence

General equations are derived for a linear beam-plasma instability in the presence of low-frequency turbulence. Within a 'quasi-linear' statistical approximation, these equations contain Langmuir wave scattering, diffusion, resonant and nonresonant anomalous absorption, and a 'plasma laser' effect. It is proposed that naturally occurring density irregularities in the solar wind may stabilize the beam-unstable Langmuir waves which occur during type III solar emissions.

Goldman, M. V.↗

Incompressible viscous flow simulations of the NFAC wind tunnel

The capabilities of an existing 3-D incompressible Navier-Stokes flow solver, INS3D, are extended and improved to solve turbulent flows through the incorporation of zero- and two-equation turbulence models. The two-equation model equations are solved in their high Reynolds number form and utilize wall functions in the treatment of solid wall boundary conditions. The implicit approximate factorization scheme is modified to improve the stability of the two-equation solver. Applications to the 3-D viscous flow inside the 80 by 120 feet open return wind tunnel of the National Full Scale Aerodynamics Complex (NFAC) are discussed and described.

Champney, Joelle Milene↗

Stabilization Of Combustion Of Sprayed Fuel

Several modifications of nozzle spraying liquid propellant into combustion chamber proposed to stabilize combustion. Proposed changes in design of nozzle alter flow field in combustion chamber, according to fluid-mechanical principles, in such way as to suppress oscillations. Various alternative nozzle configurations include bluff or toroidal body to generate turbulence. Other features helping to suppress oscillations include down-stream recirculation zones, baffles, and damping cavities. Similar modifications help to suppress oscillations in industrial combustion chambers and in commercial and domestic oil-burning furnaces.

Voecks, Gerald E.↗

Combined Bluntness and Roughness Effects on Cones at Hypersonic Speeds

This computational study investigates the effects of discrete roughness elements on a blunt cone at zero degrees angle of attack in a Mach 6 flow. Motivation was provided by experiments conducted in the Air Force Research Laboratory Mach 6 High Reynolds Number facility on a 7-degree half-angle cone with a roughness array located at 45 degrees from the apex on two nose tips of different blutness but equivalent roughness Reynolds number. Transition was only affected on the blunter cone, indicating that the transition onset is associated with the combined effects of bluntness and roughness. The present study investigates the 15.24 mm nose radius, 420 azimuthal wavenumber case via Navier-Stokes computations of the laminar base flow and instability analysis. Plane-marching parabolized stability equations (PSE) and inflow-resolvent analysis based on the three-dimensional, harmonic linearized Navier-Stokes equations (HLNSE) are used to calculate the amplification of disturbances along the roughness wake as well as over the roughness nearfield. Results show that the roughness shape can have a great impact on the characteristics of the most amplified wake instabilities. For the experimental configuration with cubic roughness elements of 15 μ m height, the flow is marginally unstable. For prismatic elements of 20 μ m height, the PSE predicts a logarithmic disturbance amplification ratio of N = 5.6 along its wake, but this ratio increases to N = 9.6 when the amplification over the roughness and separation regions is included in the inflow-resolvent analysis.

Boundary-layer transition↗

Physics and Modeling of Quasi-Steady and Forced Shock Trains

The dual-mode scramjet is a propulsive concept for contemporary hypersonic vehicles promising flight cruise speeds upwards of Mach 5. Despite decades of development and a recent resurgence in interest, such engines are not routinely deployed for commercial or military purposes. Instead, they have been confined to the realm of ground-test research and limited one-off flight test missions. One of the key stumbling blocks on the critical path to deployment of a dual-mode scramjet is the optimal design and stable operation of the isolator, a typically constant area duct that separates the scramjet inlet from the combustor. During ramjet mode at Mach numbers from 3 to 5 and during the transition to scramjet mode at higher Mach numbers, the isolator contains the pseudo-shock. Necessary to meet the back-pressure in the combustor, the pseudo-shock is typically divided into a shock dominated shock train region followed by a turbulence driven mixing region. Fluctuations in combustor pressure due to vehicle maneuvering, throttling, or combustion instabilities drive the shock train to shift streamwise within the isolator. In response to sharp pressure fluctuations, the shock train is unable to stabilize within the isolator, and the scramjet unstarts by ejecting the shock train out the inlet. This work addresses three key gaps in understanding of the isolator pseudo-shock. The phenomenon is studied in three experimental facilities to investigate the quasi-steady and unsteady fluid dynamics of the shock train in fueled and un-fueled test conditions. The strong link between boundary layer morphology and shock train structure is studied in detail and used as a tool to evaluate the state-of-the-art definitions and models of the pseudo-shock. Identifying the limitations of these models, a new process driven flux-conserved pseudo-shock model is proposed. This provides a computationally efficient model grounded in the fluid dynamic mechanisms of the pseudo-shock with better predictive capability than the current state-of-the-art across all facilities tested. The unsteady behavior of the shock train is further investigated to identify and understand two distinct forced dynamic modes: mechanically induced, low-rate forcing and high-rate combustion driven unsteadiness that sheds upstream propagating shocks. The proposed process driven model is found to be effective in the the limiting case of high-rate combustion driven unsteadiness by evaluating the model in a shock stationary reference frame. Further analysis identifies the approximately rigid chain motion of the shock train structure to mechanically induced forcing. The time delay between the onset of forced dynamics is found to saturate to a constant value and maximum shock train speed along the isolator with respect to the perturbation rate of change. The result of this work provides dual-mode scramjet designers with improved tools and physical understanding to develop more efficient, robust, and capable propulsion systems for hypersonic vehicles.

JAXA↗

A comparison of airborne eddy correlation and bulk aerodynamic methods for ocean-air turbulent fluxes during cold-air outbreaks

The viscous interfacial-sublayer model of Liu et al. (1979) is used to derive four bulk schemes (LKB, FG, D, and DB), with the flux-profile relationships of Lie et al., Francey and Garratt (1981), Dyer (1974), and Dyer and Bradley (1982). These schemes, with stability-dependent transfer coefficients, are tested against the eddy-correlation fluxes measured at the 50-m flight level above the western Atlantic Ocean during cold-air outbreaks. The bulk fluxes of momentum (tau), sensible heat (H), and latent heat (E) are found to increase with various von Karman constants. The dependence of transfer coefficients on wind speeds and roughness lengths is discussed. The transfer coefficients for tau and E agree excellently between LKB and FG. The ratio of the coefficent for H of LKB to that of FG, increasing with decreasing stability, is very sensitive to stability at low winds, but approaches the neutral value of 1.25 at high winds.

Chou, Shu-Hsien↗

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary-layer transition↗

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary layer transition↗

Theoretical studies of solar oscillations

Possible sources for the excitation of the solar 5 minute oscillations were investigated and a linear non-adiabatic stability code was applied to a preliminary study of the solar g-modes with periods near 160 minutes. Although no definitive conclusions concerning the excitation of these modes were reached, the excitation of the 5 minute oscillations by turbulent stresses in the convection zone remains a viable possibility. Theoretical calculations do not offer much support for the identification of the 160 minute global solar oscillation (reported by several independent observers) as a solar g-mode. A significant advance was made in attempting to reconcile mixing-length theory with the results of the calculations of linearly unstable normal modes. Calculations show that in a convective envelope prepared according to mixing length theory, the only linearly unstable modes are those which correspond to the turbulent eddies which are the basic element of the heuristic mixing length theory.

Goldreich, P.↗

The drive for Aircraft Energy Efficiency

NASA's Aircraft Energy Efficiency (ACEE) program, which began in 1976, has mounted a development effort in four major transport aircraft technology fields: laminar flow systems, advanced aerodynamics, flight controls, and composite structures. ACEE has explored two basic methods for achieving drag-reducing boundary layer laminarization: the use of suction through the wing structure (via slots or perforations) to remove boundary layer turbulence, and the encouragement of natural laminar flow maintenance through refined design practices. Wind tunnel tests have been conducted for wide bodied aircraft equipped with high aspect ratio supercritical wings and winglets. Maneuver load control and pitch-active stability augmentation control systems reduce fuel consumption by reducing the drag associated with high aircraft stability margins. Composite structures yield lighter airframes that in turn call for smaller wing and empennage areas, reducing induced drag for a given payload. In combination, all four areas of development are expected to yield a fuel consumption reduction of 40 percent.

James, R. L., Jr.↗

Transition flight experiments on a swept wing with suction

Flight boundary-layer transition experiments were conducted on a 30-degree swept wing with a perforated leading-edge suction panel. The transition location on the panel was changed by systematically varying the location and amount of suction. Transition from laminar to turbulent flow was due to leading-edge turbulence contamination or crossflow disturbance growth and/or Tollmien-Schlichting disturbance growth, depending on flight condition and suction variation. Amplification factor correlations with transition location were made for various suction configurations using a state-of-the-art linear stability theory which accounts for body and streamline curvature and compressibility.

Maddalon, D. V.↗