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 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRS1990

A direct-inverse method for transonic and separated flows about airfoils

A direct-inverse technique and computer program called TAMSEP that can be used for the analysis of the flow about airfoils at subsonic and low transonic freestream velocities is presented. The method is based upon a direct-inverse nonconservative full potential inviscid method, a Thwaites laminar boundary layer technique, and the Barnwell turbulent momentum integral scheme; and it is formulated using Cartesian coordinates. Since the method utilizes inverse boundary conditions in regions of separated flow, it is suitable for predicting the flow field about airfoils having trailing edge separated flow under high lift conditions. Comparisons with experimental data indicate that the method should be a useful tool for applied aerodynamic analyses.

Carlson, Leland A.↗

NASA NTRS1975

Interaction between a shock wave and a turbulent boundary layer in transonic flow

Interaction between a shock wave and an unseparated turbulent boundary layer is considered. The method of matched asymptotic expansions is used, with solutions valid in the double limit as Reynolds number tends to infinity and Mach number tends to unity. The shock is weak enough that interaction effects can be considered as perturbations to the undisturbed flow; the case considered is that where the sonic line is near the outer edge of the boundary layer. It is shown that, with order estimates for Reynolds stress perturbations, the induced wall pressure distribution can be calculated using only the two outer interaction regions, independent of a specific closure condition and that this solution is in fact a turbulent free interaction solution. A detailed analysis of the inner regions, for which an eddy viscosity model for the Reynolds shear stress is used, provides a description of the variations in velocity, temperature and density near and at the wall.

Adamson, T. C., Jr.↗

NASA NTRSDate not supplied

Split-Step Simulations of Sonic Boom Propagation Beyond the Lateral Cutoff in a Turbulent Atmosphere

Recent flight tests during the Quiet Supersonic Flights 2018 (QSF18) study reported sonic booms heard outside of the primary carpet region. In the absence of turbulence, the lateral cutoff region separates the primary sonic boom carpet from the shadow zone, where the sonic boom signal experiences significant attenuation. However, when turbulence is present in the atmospheric boundary layer (ABL), additional scattering of the sonic boom to the shadow zone region occurs. A method is presented for simulating sonic boom propagation in a turbulent atmospheric boundary layer beyond the lateral cutoff region into the shadow zone. A split-step method is used to integrate a partially one-way equation for the acoustic pressure. Inhomogeneous turbulence, representative of the ABL, is generated in the computational domain with a Fourier synthesis approach. Distributions of several loudness metrics in the shadow zone region for a sonic boom N-wave and a shaped boom are examined. Increasing both turbulence root-mean-square velocity and integral length scale are found to increase the average loudness of booms in the shadow zone. (This research is supported by the Commercial Supersonic Technology Project of the National Aeronautics and Space Administration under Grant No. 80NSSC19K1685.)

sonic boom↗

NASA NTRS1996

High Reynolds number analysis of flat plate and separated afterbody flow using non-linear turbulence models

The ability of the three-dimensional Navier-Stokes method, PAB3D, to simulate the effect of Reynolds number variation using non-linear explicit algebraic Reynolds stress turbulence modeling was assessed. Subsonic flat plate boundary-layer flow parameters such as normalized velocity distributions, local and average skin friction, and shape factor were compared with DNS calculations and classical theory at various local Reynolds numbers up to 180 million. Additionally, surface pressure coefficient distributions and integrated drag predictions on an axisymmetric nozzle afterbody were compared with experimental data from 10 to 130 million Reynolds number. The high Reynolds data was obtained from the NASA Langley 0.3m Transonic Cryogenic Tunnel. There was generally good agreement of surface static pressure coefficients between the CFD and measurement. The change in pressure coefficient distributions with varying Reynolds number was similar to the experimental data trends, though slightly over-predicting the effect. The computational sensitivity of viscous modeling and turbulence modeling are shown. Integrated afterbody pressure drag was typically slightly lower than the experimental data. The change in afterbody pressure drag with Reynolds number was small both experimentally and computationally, even though the shape of the distribution was somewhat modified with Reynolds number.

Carlson, John R.↗

NASA NTRS1974

Aerodynamic heating on 3-D bodies including the effects of entropy-layer swallowing

A relatively simple method was developed previously (authors, 1973) for calculating laminar, transitional, and turbulent heating rates on three-dimensional bodies in hypersonic flows. This method was shown to yield reasonably accurate results for laminar heating on blunted circular and elliptical cones and an earlier version of the space shuttle vehicle. As the boundary layer along the surface grows, more and more of the inviscid-flow mass is entrained into the boundary layer, and the streamlines which passed through the nearly normal portion of the bow shock wave are 'swallowed' by the boundary layer. This phenomenon is often referred to as entropy-layer or streamline swallowing, and it can have a significant effect on the calculated heating rates. An approximate, yet simple, method for including the effects of entropy-layer swallowing in the heating-rate calculations is given.

Dejarnette, F. R.↗

NASA NTRS1982

Numerical modelling of turbulent flow in a combustion tunnel

A numerical technique is presented for the analysis of turbulent flow associated with combustion. The technique uses Chorin's random vortex method (rvm), an algorithm capable of tracing the action of elementary turbulent eddies and their cumulative effects without imposing any restriction upon their motion. In the past, the rvm has been used with success to treat nonreacting turbulent flows, revealing in particular the mechanics of large-scale flow patterns, the so-called coherent structures. Introduced here is a flame propagation algorithm, also developed by Chorin, in conjunction with volume sources modelling the mechanical effects of the exothermic process of combustion. As an illustration of its use, the technique is applied to flow in a combustion tunnel where the flame is stabilized by a back-facing step. Solutions for both nonreacting and reacting flow fields are obtained which satisfactorily describe the essential features of turbulent combustion in a lean propane-air mixture that were observed in the laboratory by means of high speed Schlieren photography.

Ghoniem, A. F.↗

NASA NTRS2003

Comparison of Artificial Compressibility Methods

Various artificial compressibility methods for calculating three-dimensional, steady and unsteady, laminar and turbulent, incompressible Navier-Stokes equations are compared in this work. Each method is described in detail along with appropriate physical and numerical boundary conditions. Analysis of well-posedness and numerical solutions to test problems for each method are provided. A comparison based on convergence behavior, accuracy, stability and robustness is used to establish the relative positive and negative characteristics of each method.

Kiris, Cetin↗

NASA NTRS1989

Turbulent and directed plasma motions in solar flares

An improved method for fitting asymmetric soft X-ray line profiles from solar flares is presented. A two-component model is used where one component represents the total emission from directed upflow plasma and the other the emission from the plasma at rest. Unlike previous methods, the width of the moving component is independent from that of the stationary component. Time variations of flare plasma characteristics (i.e., temperature, emission measure of moving and stationary plasma, upflow and turbulent velocities) are derived from the Ca XIX and Fe XXV spectra recorded by the Bent Crystal Spectrometer on the Solar Maximum Mission. The fitting technique provides a statistical estimation for the uncertainties in the fitting parameters. The relationship between the directed and turbulent motions has been studied, and a correlation of the random and directed motions has been found in some flares with intensive plasma upflows. Mean temperatures of the upflowing and stationary plasmas are compared for the first time from ratios of calcium to iron X-ray line intensities. Finally, evidence for turbulent motions and the possibility of plasma upflow late into the decay phase is presented and discussed.

Fludra, A.↗

DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1992

A PDF closure model for compressible turbulent chemically reacting flows

The objective of the proposed research project was the analysis of single point closures based on probability density function (pdf) and characteristic functions and the development of a prediction method for the joint velocity-scalar pdf in turbulent reacting flows. Turbulent flows of boundary layer type and stagnation point flows with and without chemical reactions were be calculated as principal applications. Pdf methods for compressible reacting flows were developed and tested in comparison with available experimental data. The research work carried in this project was concentrated on the closure of pdf equations for incompressible and compressible turbulent flows with and without chemical reactions.

Kollmann, W.↗

NASA NTRS1993

Fluid dynamics and convective heat transfer in impinging jets through implementation of a high resolution liquid crystal technique

A combined convective heat transfer and fluid dynamics investigation in a turbulent round jet impinging on a flat surface is presented. The experimental study uses a high resolution liquid crystal technique for the determination of the convective heat transfer coefficients on the impingement plate. The heat transfer experiments are performed using a transient heat transfer method. The mean flow and the character of turbulent flow in the free jet is presented through five hole probe and hot wire measurements, respectively. The flow field character of the region near the impingement plate plays an important role in the amount of convective heat transfer. Detailed surveys obtained from five hole probe and hot wire measurements are provided. An extensive validation of the liquid crystal based heat transfer method against a conventional technique is also presented. After a complete documentation of the mean and turbulent flow field, the convective heat transfer coefficient distributions on the impingement plate are presented. The near wall of the impingement plate and the free jet region is treated separately. The current heat transfer distributions are compared to other studies available from the literature. The present paper contains complete sets of information on the three dimensional mean flow, turbulent velocity fluctuations, and convective heat transfer to the plate. The experiments also prove that the present nonintrusive heat transfer method is highly effective in obtaining high resolution heat transfer maps with a heat transfer coefficient uncertainty of 5.7 percent.

Kim, K.↗

NASA NTRS1998

Assessment of an Euler-Interacting Boundary Layer Method Using High Reynolds Number Transonic Flight Data

Flight-measured high Reynolds number turbulent-flow pressure distributions on a transport wing in transonic flow are compared to unstructured-grid calculations to assess the predictive ability of a three-dimensional Euler code (USM3D) coupled to an interacting boundary layer module. The two experimental pressure distributions selected for comparative analysis with the calculations are complex and turbulent but typical of an advanced technology laminar flow wing. An advancing front method (VGRID) was used to generate several tetrahedral grids for each test case. Initial calculations left considerable room for improvement in accuracy. Studies were then made of experimental errors, transition location, viscous effects, nacelle flow modeling, number and placement of spanwise boundary layer stations, and grid resolution. The most significant improvements in the accuracy of the calculations were gained by improvement of the nacelle flow model and by refinement of the computational grid. Final calculations yield results in close agreement with the experiment. Indications are that further grid refinement would produce additional improvement but would require more computer memory than is available. The appendix data compare the experimental attachment line location with calculations for different grid sizes. Good agreement is obtained between the experimental and calculated attachment line locations.

Bonhaus, Daryl L.↗

NASA NTRS1951

On the turbulent friction layer for rising pressure

Among the information presented are included displacement, momentum, and kinetic energy thicknesses, shearing stress distributions across boundary layer, and surface friction coefficients. The Gruschwitz method and its modifications are examined and tested. An energy theorem for the turbulent boundary layer is introduced and discussed but does not lead to a method for the prediction of the behavior of the turbulent boundary layer because relations for the shearing stress and the surface friction are lacking.

SKIN FRICTION - PLATES↗

DOE OSTI2023

Drift-ordered fluid vorticity equation with energy consistency

Although drift-ordered fluid models are widely applied in tokamak edge turbulence simulations, the models used are acknowledged not to conserve energy or even electrical charge. The present paper aims to remove many of the existing pitfalls in drift-fluid models, however, with the objective of finding a solution simple enough to be implemented in numerical applications. Our main result is an improved version of the drift-Braginskii equations involving a generalized vorticity function. In the new drift-Braginskii system, the quasi-neutrality condition translates into a transport equation for a generalized vorticity, expressed in conservation form, and related to the total mass-weighted circulation. It is found that kinetic energy conservation can be achieved if the polarization flow is defined recursively. In conclusion, the resulting model conserves the kinetic energy associated with E×B and diamagnetic flows and retains the associated perpendicular kinetic energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

Implicit Large-Eddy Simulation of Single-Injector Cooling Flow

The turbulent heat flux vector plays a critical role in aero-propulsion flow. Currently, traditional RANS based methods do not currently have reliable closures which can accurately predict the turbulent transport of heat. Scale Resolving Simulations can provide critical insight into this fundamental fluid mechanism. Large-eddy Simulations (LES) were performed on the single-hole film-cooling experimental configuration from the recent experiments conducted at NASA Glenn Research Center (GRC)[1]. The experiment of Wernet et al. obtained mean and fluctuating temperatures which can be used to validate LES results predicting the turbulent heat flux vector. The FDL3DI flow solver was used to complete simulations of two blowing ratios of the THX - III dataset. A grid sensitivity study was also performed. Initial LES results show promising agreement between the simulation and experimental measurements.

High-Order methods↗

NASA NTRS1982

Numerical simulation of wall-bounded turbulent shear flows

Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included.

Moin, P.↗

NASA NTRS1982

Numerical simulation of wall-bounded turbulent shear flows

Moin, P.↗

NASA NTRS1992

Turbulence modeling for hypersonic flight

The objective of the present work is to develop, verify, and incorporate two equation turbulence models which account for the effect of compressibility at high speeds into a three dimensional Reynolds averaged Navier-Stokes code and to provide documented model descriptions and numerical procedures so that they can be implemented into the National Aerospace Plane (NASP) codes. A summary of accomplishments is listed: (1) Four codes have been tested and evaluated against a flat plate boundary layer flow and an external supersonic flow; (2) a code named RANS was chosen because of its speed, accuracy, and versatility; (3) the code was extended from thin boundary layer to full Navier-Stokes; (4) the K-omega two equation turbulence model has been implemented into the base code; (5) a 24 degree laminar compression corner flow has been simulated and compared to other numerical simulations; and (6) work is in progress in writing the numerical method of the base code including the turbulence model.

Bardina, Jorge E.↗