SEARCH · Engineering Papers

Results for “turbulence methods”

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 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1930

The Measurement of Fluctuations of Air Speed by the Hot-Wire Anemometer

The hot-wire anemometer suggests itself as a promising method for measuring the fluctuating air velocities found in turbulent flow. The only obstacle is the presence of a lag due to the limited energy input which makes even a fairly small wire incapable of following rapid fluctuations with accuracy. This paper gives the theory of the lag and describes an experimental arrangement for compensating for the lag for frequencies up to 100 or more per second when the amplitude of the fluctuation is not too great. An experimental test of the accuracy of compensation and some results obtained with the apparatus in a wind-tunnel air stream are described. While the apparatus is very bulky in its present form, it is believed possible to develop a more portable arrangement. (author)

Dryden, H L↗

NASA NTRS1994

Assimilation by lunar mare basalts: Melting of crustal material and dissolution of anorthite

We discuss techniques for calculating the amount of crustal assimilation possible in lunar magma chambers and dikes based on thermal energy balances, kinetic rates, and simple fluid mechanical constraints. Assuming parent magmas of picritic compositions, we demonstrate the limits on the capacity of such magmas to melt and dissolve wall rock of anorthitic, troctolitic, noritic, and KREEP (quartz monzodiorite) compositions. Significant melting of the plagioclase-rich crustal lithologies requires turbulent convection in the assimilating magma and an efficient method of mixing in the relatively buoyant and viscous new melt. Even when this occurs, the major element chemistry of the picritic magmas will change by less than 1-2 wt %. Diffusion coefficients measured for Al2O3 from an iron-free basalt and an orange glass composition are 10(exp -12) sq m/s at 1340 C and 10(exp -11) sq m/s at 1390 C. These rates are too slow to allow dissolution of plagioclase to significantly affect magma compositions. Picritic magmas can melt significant quantities of KREEP, which suggests that their trace element chemistry may still be affected by assimilation processes; however, mixing viscous melts of KREEP composition with the fluid picritic magmas could be prohibitively difficult. We conclude that only a small part of the total major element chemical variation in the mare basalt and volcanic glass collection is due to assimilation/fractional crystallization processes near the lunar surface. Instead, most of the chemical variation in the lunar basalts and volcanic glasses must result from assimilation at deeper levels or from having distinct source regions in a heterogeneous lunar mantle.

Finnila, A. B.↗

NASA NTRS1996

SCISEAL: A CFD Code for Analysis of Fluid Dynamic Forces in Seals

A 3D CFD code, SCISEAL, has been developed and validated. Its capabilities include cylindrical seals, and it is employed on labyrinth seals, rim seals, and disc cavities. State-of-the-art numerical methods include colocated grids, high-order differencing, and turbulence models which account for wall roughness. SCISEAL computes efficient solutions for complicated flow geometries and seal-specific capabilities (rotor loads, torques, etc.).

Althavale, Mahesh M.↗

NASA NTRS1994

Assimilation by Lunar Mare Basalts: Melting of Crustal Material and Dissolution of Anorthite

We discuss techniques for calculating the amount of crustal assimilation possible in lunar magma chambers and dikes based on thermal energy balances, kinetic rates, and simple fluid mechanical constraints. Assuming parent magmas of picritic compositions, we demonstrate the limits on the capacity of such magmas to melt and dissolve wall rock of anorthitic, troctolitic, noritic, and KREEP (quartz monzodiorite) compositions. Significant melting of the plagioclase-rich crustal lithologies requires turbulent convection in the assimilating magma and an efficient method of mixing in the relatively buoyant and viscous new melt. Even when this occurs, the major element chemistry of the picritic magmas will change by less than 1-2 wt %. Diffusion coefficients measured for Al2O3 from an iron-free basalt and an orange glass composition are 10(exp -12) m(exp 2) s(exp -1) at 1340 C and 10(exp -11) m(exp 2) s(exp -1) at 1390 C. These rates are too slow to allow dissolution of plagioclase to significantly affect magma compositions. Picritic magmas can melt significant quantities of KREEP, which suggests that their trace element chemistry may still be affected by assimilation processes; however, mixing viscous melts of KREEP composition with the fluid picritic magmas could be prohibitively difficult. We conclude that only a small part of the total major element chemical variation in the mare basalt and volcanic glass collection is due to assimilation/fractional crystallization processes near the lunar surface. Instead, most of the chemical variation in the lunar basalts and volcanic glasses must result from assimilation at deeper levels or from having distinct source regions in a heterogeneous lunar mantle.

Finnila, A. B.↗

NASA NTRS1999

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

NASA NTRS1994

Progress In Incompressible Pump Flow Calculations

Steady and unsteady flows for propulsion systems are efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard kappa - epsilon model equations. The resulting computer code is applied to the flow analysis inside an advanced rocket pump impeller in steadily rotating reference frames. Numerical results are compared with experimental measurements. The effects of exit and shroud cavities with the leak-age flow are investigated. Time-accurate incompressible Navier-Stokes formulation with the overlapped grid scheme capability was evaluated by using MIT flapping foil experiment. The grid dependency, turbulence model effects, and the effect of order of differencing were investigated. Numerical results were compared against experimental data. The resulting procedure were applied to unsteady flapping foil calculations. Two upstream NACA 0025 foils perform high frequency synchronized motion and generate unsteady flow conditions to the downstream larger stationary foil. Comparison between unsteady experimental data and numerical results from two different moving boundary procedures will be presented.

Kiris, Cetin↗

NASA NTRS1999

An Elasticity-Based Mesh Scheme Applied to the Computation of Unsteady Three-Dimensional Spoiler and Aeroelastic Problems

This paper presents a modification of the spring analogy scheme which uses axial linear spring stiffness with selective spring stiffening/relaxation. An alternate approach to solving the geometric conservation law is taken which eliminates the need for storage of metric Jacobians at previous time steps. Efficiency and verification are illustrated with several unsteady 2-D airfoil Euler computations. The method is next applied to the computation of the turbulent flow about a 2-D airfoil and wing with two and three- dimensional moving spoiler surfaces, and the results compared with Benchmark Active Controls Technology (BACT) experimental data. The aeroelastic response at low dynamic pressure of an airfoil to a single large scale oscillation of a spoiler surface is computed. This study confirms that it is possible to achieve accurate solutions with a very large time step for aeroelastic problems using the fluid solver and aeroelastic integrator as discussed in this paper.

Bartels, Robert E.↗

NASA NTRS2003

Sound Sources Identified in High-Speed Jets by Correlating Flow Density Fluctuations With Far-Field Noise

Noise sources in high-speed jets were identified by directly correlating flow density fluctuation (cause) to far-field sound pressure fluctuation (effect). The experimental study was performed in a nozzle facility at the NASA Glenn Research Center in support of NASA s initiative to reduce the noise emitted by commercial airplanes. Previous efforts to use this correlation method have failed because the tools for measuring jet turbulence were intrusive. In the present experiment, a molecular Rayleigh-scattering technique was used that depended on laser light scattering by gas molecules in air. The technique allowed accurate measurement of air density fluctuations from different points in the plume. The study was conducted in shock-free, unheated jets of Mach numbers 0.95, 1.4, and 1.8. The turbulent motion, as evident from density fluctuation spectra was remarkably similar in all three jets, whereas the noise sources were significantly different. The correlation study was conducted by keeping a microphone at a fixed location (at the peak noise emission angle of 30 to the jet axis and 50 nozzle diameters away) while moving the laser probe volume from point to point in the flow. The following figure shows maps of the nondimensional coherence value measured at different Strouhal frequencies ([frequency diameter]/jet speed) in the supersonic Mach 1.8 and subsonic Mach 0.95 jets. The higher the coherence, the stronger the source was.

Panda, Jayanta↗

NASA NTRS1987

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.↗

NASA NTRS1997

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↗

NASA NTRS1975

Cross-correlation methods for studying near- and farfield noise characteristics of flow-surface interactions

Systematic methods based on cross-correlation techniques are presented for experimental studies of near- and farfield noise characteristics in airflow-surface interaction problems. Analyses show that, in near- and farfields, the noise characteristics due to the surface contribution of fluctuating pressures and velocities and due to the volume contribution of the turbulence in the flow can be determined separately. Both farfield noise intensities and nearfield acoustic energy fluxes can be expressed in terms of appropriate cross correlations. These correlations can be obtained by making microphone measurements in the farfield, in the nearfield, and on the surface. Examples of the applications to the noise field generated by impinging jets, by surface blowing jets, and by turbulent flow over trailing edges are presented. Advantages of the present method over conventional methods are also discussed.

Pan, Y. S.↗

NASA NTRS1991

Numerical investigation of an internal layer in turbulent flow over a curved hill

A numerical investigation of incompressible and compressible turbulent flows over strongly curved surfaces is presented. The turbulent flow equations are solved by a pressure based Navier-Stokes equations solver. In the method, the conservation of mass equation is replaced by a pressure correction equation applicable for both compressible and incompressible flows. The turbulence is described by a multiple time scale turbulence model supplemented with a near-wall turbulence model. The numerical results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. The development of the internal layer is attributed to the enormous mean flow strain rate caused by the streamline curvature. The external boundary layer flow responds rather slowly to the streamline curvature. Thus, the turbulence field of the forward corner of the curved hill is characterized by two turbulence fields interacting with each other. The turbulence intensity of the internal layer is much stronger than that of the external boundary layer, so that the development of a new boundary layer in the downstream region of the curved hill depends mostly on the internal layer. These numerical results are in good agreement with the measured data, and show that the turbulence model can resolve the turbulence field subjected to the strong streamline curvature.

Kim, S.-W.↗

NASA NTRS1990

A review of instability and noise propagation in supersonic flows

Originally analytical and numerical models were to be developed for noise production in supersonic jets, wakes and free shear layers. While the effort was concentrated initially on these aspects, other topics were also pursued, most were of interest to the Jet Noise Group of the Aeroacoustics Branch. An overview is given of subjects reviewed and the investigations that were carried out. A significant effort was devoted to numerically predicting the flow field of a turbulent supersonic wall jet. This information is necessary for computing the pressure in the far field. The wall jet was selected because it represents a generic flow that can be associated with plug nozzle in supersonic engines. It combines the characteristic of a boundary layer with that of a free shear flow. The spatially evolving flow obtained using Dash's code would form the input for the stability analysis program. This analysis would determine the large scale instability wave within the flow. The far field pressure can be computed from the shape of the evolving large scale structure by asymptotic methods. Flow characteristics obtained from a program that analyses the turbulent downstream supersonic flow in a nozzle are described and compared with experimental results.

Daudpota, Q. Isa↗

NASA NTRS1975

The simulation of turbulent boundary layer separation on multi-element infinite swept wings

An improved method of simulating the effect of flow separation on potential flow over two-dimensional and infinite swept multi-element wings was developed. Two types of flow separation are considered: at the wing trailing edge, and at the cove region between the main wing and the flap. Potential flow is assumed outside the boundary layer and the separated wake. A previously obtained computer program (VIP) is modified by matching the outflow distribution for the potential flow such that the potential flow part of the real flow is duplicated and the results agree with the experimental observations. Results are given for the NACA 2412 airfoil. Iterative drag calculations are presented.

Dvorak, F. A.↗

NASA NTRS2020

Implicit Large-Eddy Simulations of Hot and Cold Supersonic Jets in Loci-CHEM

Turbulence is a critical factor in most aero- and propulsion-related flows. For many years, the industry-standard approach to incorporating the effect of turbulence in computational fluid dynamics (CFD)calculations has been the Reynolds Averaged Navier Stokes (RANS) modeling approach. RANS implicitly time-averages turbulent motion, and models its effect on the mixing of species, momentum, and energy in a flow field. The approach is economical and often allows a CFD calculation to proceed to steady-state. However, it makes several simplifying assumptions, and the results will only be as good as the RANS turbulence model. Meanwhile, the last two decades have witnessed rapid and explosive growth in computational power, the development of CFD codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. CFD engineering calculations of 108– 109 finite-volume cells have become routine. Thus, large eddy simulation (LES), which attempts to directly resolve the unsteady motion of the largest scales of turbulence, is increasingly of interest for many fluids engineering problems. Additionally, hybrid RANS-LES approaches such as Spalart’s Detached Eddy Simulation (DES), which attempt to marry the best strengths of both methods, are in active use and development. Because it seeks to accurately resolve the unsteady motion of the largest scales of turbulence, LES places certain requirements on the numerical method used to discretize and advance the Navier Stokes equations. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide range of length scales.

C I Morris↗

NASA NTRSDate not supplied

Sparse Linear Algebra Toolkit for Computational Aerodynamics

Finding solutions to sparse linear systems of equations is an essential step in Computational Engineering applications of interest to NASA. Linear systems of equations are composed and solved in almost every computational engineering application. The characteristics of linear systems vary greatly from one application to another. Accordingly, there are a wide variety of methods for the solution of linear systems of equations. The operations and methods prepared by the authors are focused on linear systems of interest to NASA, primarily those associated with Computational Fluid Dynamics (CFD), Aeroelasticity, and Aeroacoustics. The Sparse Linear Algebra Toolkit (SLAT) is a coordinated collection of software featuring operations, methods, and data structures that are useful when solving sparse linear systems of equations on modern computer architectures. The implemented operations and methods are designed and tuned for parallelism in shared memory, in distributed memory, and across the hybrid combination of distributed-shared memory. The toolkit includes novel methods and implementations for modern architectures and facilitates development of new approaches for meeting NASA’s evolving computational engineering challenges using evolving computer architectures that are not available in vendor libraries. In this paper, significant features and interfaces within SLAT are presented and verified for simulations performed with NASA’s CFD solver, FUN3D. The runtime and scaling performance of the Generalized Minimum Residual (GMRES) method implemented in SLAT is analyzed for the linear subproblems within the solution of turbulent Navier-Stokes equations employed in the simulation of high-lift configurations. Prior to this work, the SPARSKIT GMRES implementation was the only Krylov subspace method available within FUN3D. A strong scaling study shows the SLAT GMRES implementation facilitates accurate Reynolds-averaged Navier-Stokes CFD solutions between 15% and 56% faster than the SPARSKIT GMRES implementation.

Stephen L Wood↗

NASA NTRS1983

Experimental and theoretical study of combustion jet ignition

A combustion jet ignition system was developed to generate turbulent jets of combustion products containing free radicals and to discharge them as ignition sources into a combustible medium. In order to understand the ignition and the inflammation processes caused by combustion jets, the studies of the fluid mechanical properties of turbulent jets with and without combustion were conducted theoretically and experimentally. Experiments using a specially designed igniter, with a prechamber to build up and control the stagnation pressure upstream of the orifice, were conducted to investigate the formation processes of turbulent jets of combustion products. The penetration speed of combustion jets has been found to be constant initially and then decreases monotonically as turbulent jets of combustion products travel closer to the wall. This initial penetration speed to combustion jets is proportional to the initial stagnation pressure upstream of the orifice for the same stoichiometric mixture. Computer simulations by Chorin's Random Vortex Method implemented with the flame propagation algorithm for the theoretical model of turbulent jets with and without combustion were performed to study the turbulent jet flow field. In the formation processes of the turbulent jets, the large-scale eddy structure of turbulence, the so-called coherent structure, dominates the entrainment and mixing processes. The large-scale eddy structure of turbulent jets in this study is constructed by a series of vortex pairs, which are organized in the form of a staggered array of vortex clouds generating local recirculation flow patterns.

Chen, D. Y.↗

NASA NTRS1975

Turbulent flow between a rotating disk and a parallel wall

A scalar effective viscosity method is used to calculate the incompressible flow in the gap between two infinite, parallel smooth disks, one of them rotating and the other stationary. The flow is calculated for two values of the gap Reynolds number omega(s squared)/nu, namely 296 and 2852, where omega is the disk angular velocity, s the gap width, and nu the kinematic viscosity of the fluid; and over a range of radial Reynolds number omega(r squared)/nu from 0 to 10 million, where r is the radius from the axis of rotation. A numerical shooting technique is used to match the two boundary-layer flows arising on the disk and wall. Gradual transition is assumed for radial Reynolds numbers from 160,000 to 250,000. The results are in good agreement with Daily and Nece's measurements of velocity profiles and disk friction drag for enclosed rotating disks, and at sufficiently large radius the local skin friction approaches that of turbulent Couette flow.

Cooper, P.↗