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,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Spectral Methods for Computational Fluid Dynamics

As a tool for large-scale computations in fluid dynamics, spectral methods were prophesized in 1944, born in 1954, virtually buried in the mid-1960's, resurrected in 1969, evangalized in the 1970's, and catholicized in the 1980's. The use of spectral methods for meteorological problems was proposed by Blinova in 1944 and the first numerical computations were conducted by Silberman (1954). By the early 1960's computers had achieved sufficient power to permit calculations with hundreds of degrees of freedom. For problems of this size the traditional way of computing the nonlinear terms in spectral methods was expensive compared with finite-difference methods. Consequently, spectral methods fell out of favor. The expense of computing nonlinear terms remained a severe drawback until Orszag (1969) and Eliasen, Machenauer, and Rasmussen (1970) developed the transform methods that still form the backbone of many large-scale spectral computations. The original proselytes of spectral methods were meteorologists involved in global weather modeling and fluid dynamicists investigating isotropic turbulence. The converts who were inspired by the successes of these pioneers remained, for the most part, confined to these and closely related fields throughout the 1970's. During that decade spectral methods appeared to be well-suited only for problems governed by ordinary diSerential eqllations or by partial differential equations with periodic boundary conditions. And, of course, the solution itself needed to be smooth. Some of the obstacles to wider application of spectral methods were: (1) poor resolution of discontinuous solutions; (2) inefficient implementation of implicit methods; and (3) drastic geometric constraints. All of these barriers have undergone some erosion during the 1980's, particularly the latter two. As a result, the applicability and appeal of spectral methods for computational fluid dynamics has broadened considerably. The motivation for the use of spectral methods in numerical calculations stems from the attractive approximation properties of orthogonal polynomial expansions.

Zang, T. A.↗

NASA NTRS1989

An approximate viscous shock layer method for calculating the hypersonic flow over blunt-nosed bodies

The approximate axisymmetric method presented for accurately calculating the surface and flowfield properties of fully viscous hypersonic flow over blunt-nosed bodies incorporates the turbulence model of Cebeci-Smith (1970) and the equilibrium air tables of Hansen (1959). The method is faster than the parabolized Navier-Stokes or viscous shock layer solvers that it could replace for preliminary design determinations. Surface heat transfer and pressure predictions for the present method are comparable with the more accurate viscous shock layer method as well as flight test and wind tunnel data. A starting solution is not required.

Grantz, A. C.↗

NASA NTRS1992

Analysis of an RNG based turbulence model for separated flows

A two-equation turbulence model of the K-epsilon type was recently derived by using Renormalization Group (RNG) methods. It was later reported that this RNG based model yields substantially better predictions than the standard K-epsilon model for turbulent flow over a backward facing step - a standard test case used to benchmark the performance of turbulence models in separated flows. The improvements obtained from the RNG K-epsilon model were attributed to the better treatment of near wall turbulence effects. In contrast to these earlier claims, it is shown in this paper that the original version of the RNG K-epsilon model substantially underpredicts the reattachment point in the backstep problem. This is a deficiency that is traced to the modeling of the production of dissipation term. However, with the most recent improvements in the RNG K-epsilon model, excellent results for the backstep problem are now obtained.

Speziale, C. G.↗

NASA NTRS1980

A comprehensive comparison between experiment and prediction for a transonic turbulent separated flow

Attempts to predict surface pressure distributions on lifting surfaces have been relatively unsuccessful in the transonic regime when the shock wave is of sufficient strength to produce an extensive region of turbulent separated flow. For these conditions, the viscous flow behavior must be accurately described even to obtain reasonable predictions of surface pressure. The present paper addresses this problem. Detailed comparisons between prediction and experiment are made for a transonic, turbulent boundary-layer separation (freestream Mach number = 0.875) for which the turbulent flow properties (including the turbulent Reynolds stress) had been measured by the laser velocimeter technique from upstream of the separated region through reattachment. The flow was generated on an axisymmetric 'bump' model designed to simulate the flow on an airfoil at transonic conditions. The numerical methods used in the comparisons include the solution of the time-dependent, mass-averaged Navier-Stokes equations, and the solution of the compressible boundary-layer equations by the inverse method. Solutions were obtained for the well established Cebeci-Smith algebraic turbulence model and the more recently developed Wilcox-Rubesin two-equation turbulence model.

Johnson, D. A.↗

NASA NTRS1993

Large eddy simulation and direct numerical simulation of high speed turbulent reacting flows

The objective of this research is to make use of Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first phase of this research conducted within the past three years have been directed in several issues pertaining to intricate physics of turbulent reacting flows. In our previous 5 semi-annual reports submitted to NASA LaRC, as well as several technical papers in archival journals, the results of our investigations have been fully described. In this progress report which is different in format as compared to our previous documents, we focus only on the issue of LES. The reason for doing so is that LES is the primary issue of interest to our Technical Monitor and that our other findings were needed to support the activities conducted under this prime issue. The outcomes of our related investigations, nevertheless, are included in the appendices accompanying this report. The relevance of the materials in these appendices are, therefore, discussed only briefly within the body of the report. Here, results are presented of a priori and a posterior analyses for validity assessments of assumed Probability Density Function (PDF) methods as potential subgrid scale (SGS) closures for LES of turbulent reacting flows. Simple non-premixed reacting systems involving an isothermal reaction of the type A + B yields Products under both chemical equilibrium and non-equilibrium conditions are considered. A priori analyses are conducted of a homogeneous box flow, and a spatially developing planar mixing layer to investigate the performance of the Pearson Family of PDF's as SGS models. A posteriori analyses are conducted of the mixing layer using a hybrid one-equation Smagorinsky/PDF SGS closure. The Smagorinsky closure augmented by the solution of the subgrid turbulent kinetic energy (TKE) equation is employed to account for hydrodynamic fluctuations, and the PDF is employed for modeling the effects of scalar fluctuations. The implementation of the model requires the knowledge of the local values of the first two SGS moments. These are provided by additional modeled transport equations. In both a priori and a posteriori analyses, the predicted results are appraised by comparison with subgrid averaged results generated by DNS. Based on these results, the paths to be followed in future investigations are identified.

Adumitroaie, V.↗

NASA NTRS1994

Nonlinear decay of magnetic helicity in magnetohydrodynamic turbulence with a mean magnetic field

We show that the magnetic helicity associated with fluctuations in homogeneous incompressible magnetohydrodynamic (MHD) turbulence with a mean magnetic field decays in time because of nonlinear processes. Evidence is obtained numerically, by use of both dissipative and nondissipative spectral method simulations. The described effect stands in contrast to expectations based on studies of MHD turbulence without an applied mean field, in which magnetic helicity is transferred nonlinearly to long wavelengths and is preserved in time because of selective decay when dissipation is present. The process of nonlinear decay is described in terms of a generalized ideal helicity invariant and is characterized by the transient production of a mean induced electric field aligned with the applied magnetic field, an effect reminiscent of the alpha effect in dynamo theory. A simple phenomenological model for the decay process is proposed.

Stribling, Troy↗

NASA NTRS1987

Numerical optimization of conical flow waveriders including detailed viscous effects

A family of optimized hypersonic waveriders is generated and studied wherein detailed viscous effects are included within the optimization process itself. This is in contrast to previous optimized waverider work, wherein purely inviscid flow is used to obtain the waverider shapes. For the present waveriders, the undersurface is a streamsurface of an inviscid conical flowfield, the upper surface is a streamsurface of the inviscid flow over a tapered cylinder (calculated by the axisymmetric method of characteristics), and the viscous effects are treated by integral solutions of the boundary layer equations. Transition from laminar to turbulent flow is included within the viscous calculations. The optimization is carried out using a nonlinear simplex method. The resulting family of viscous hypersonic waveriders yields predicted high values of lift/drag, high enough to break the L/D barrier based on experience with other hypersonic configurations. Moreover, the numerical optimization process for the viscous waveriders results in distinctly different shapes compared to previous work with inviscid-designed waveriders. Also, the fine details of the viscous solution, such as how the shear stress is distributed over the surface, and the location of transition, are crucial to the details of the resulting waverider geometry. Finally, the moment coefficient variations and heat transfer distributions associated with the viscous optimized waveriders are studied.

Bowcutt, Kevin G.↗

NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

NASA NTRS1988

Some numerical and physical aspects of unsteady Navier-Stokes computations over airfoils using dynamic meshes

An upwind-biased implicit approximate factorization algorithm is applied to several unsteady flows on dynamic meshes. The thin-layer form of the compressible Navier-Stokes equations is used to solve both laminar and turbulent flows over airfoils pitching about the quarter chord. Numerical aspects of the solutions are investigated, including grid and time step effects. Two methods for determining fluxes - flux-vector splitting and flux-difference splitting - are compared. Flux-difference splitting predicts results more accurately than flux-vector splitting on a coarse mesh, but both methods agree on a fine mesh. Physical aspects of the computations are also examined. An equilibrium turbulent boundary layer model computes generally better unsteady results in comparison with experiment than a nonequilibrium model for the transonic case analyzed. Also, the size and location of the primary shed vortex for an airfoil pitching up at a constant rate is calculated in good agreement with experiment for two pitch rates.

Rumsey, Christopher L.↗

NASA NTRS1977

Instrumentation techniques for studying heterogeneous combustion

Velocity measurements in flames and the simultaneous measurement of droplet size and velocity with the aid of laser anemometry are considered along with methods of particle size measurement and techniques based on laser Raman spectroscopy. Attention is given to high-speed photography and holography, computer compensated thermocouples, and the determination of turbulence characteristics. The employment of suction probes is also discussed, taking into account sampling methods, the effects of temperature change in sampling probes on particle formation, and the chemical analysis of samples.

Chigier, N. A.↗

NASA NTRSDate not supplied

Numerical Methods for Jet Noise Predictions Using the Generalized Acoustic Analogy

Numerical methods used in several implementations of the Generalized Acoustic Analogy for the prediction of turbulent jet noise are discussed. The equation governing sound propagation in the analogy formulation is presented along with analytical simplifications to reduce the computational burden of obtaining solutions. Numerical methods for the resulting boundary-value problem for two classes of jet flows, axisymmetric and non-axisymmetric mean flows, are discussed in detail. Example results from problems of the two classed of jets are presented.

jet noise↗

NASA NTRS1985

Semiempirical method of determining flow coefficients for pitot rake mass flow rate measurements

Flow coefficients applicable to area-weighted pitot rake mass flow rate measurements are presented for fully developed, turbulent flow in an annulus. A turbulent velocity profile is generated semiempirically for a given annulus hub-to-tip radius ratio and integrated numerically to determine the ideal mass flow rate. The calculated velocities at each probe location are then summed, and the flow rate as indicated by the rake is obtained. The flow coefficient to be used with the particular rake geometry is subsequently obtained by dividing the ideal flow rate by the rake-indicated flow rate. Flow coefficients ranged from 0.903 for one probe placed at a radius dividing two equal areas to 0.984 for a 10-probe area-weighted rake. Flow coefficients were not a strong function of annulus hub-to-tip radius ratio for rakes with three or more probes. The semiempirical method used to generate the turbulent velocity profiles is described in detail.

Trefny, C. J.↗

NASA NTRS1995

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗

NASA NTRS1979

Computation of supersonic viscous flows over ogive-cylinders at angle of attack

The parabolic Navier-Stokes (PNS) marching finite-difference method is applied to 3-D viscous flow over pointed ogive-cylinders, and to turbulent flow over a cone. Ogive computations were performed using the new technique recently reported by Vigneron, Rakich, and Tannehill. Comparison is made with experiment and inviscid computations. The present results show that this method, which neglects part of the pressure gradient in the x-momentum equation, is nevertheless valid for flows with a strong favorable pressure gradient. In addition, turbulent separated flow over a cone has been computed using the older PNS code due to Lubard and Helliwell. It is found that one must freeze the turbulent eddy-viscosity model upstream of 3-D separation to get agreement with experiment.

Rakich, J. V.↗

NASA NTRS1947

Some Recent Contributions to the Study of Transition and Turbulent Boundary Layers

The first part of this paper reviews the present state of the problem of the instability of laminar boundary layers which has formed an important part of the general lectures by von Karman at the first and fourth Congresses and by Taylor at the fifth Congress. This problem may now be considered as essentially solved as the result of work completed since 1938. When the velocity fluctuations of the free-stream flow are less than 0.1 percent of the mean speed, instability occurs as described by the well-known Tollmien-Schlichting theory. The Tollmien-Schlichting waves were first observed experimentally by Schubauer and Skramstad in 1940. They devised methods of introducing controlled small disturbances and obtained measured values of frequency, damping, and wave length at various Reynolds numbers which agreed well with the theoretical results. Their experimental results were confirmed by Liepmann. Much theoretical work was done in Germany in extending the Tol1mien-Schlichting theory to other boundary conditions, in particular to flow along a porous wall to which suction is applied for removing part of the boundary layer. The second part of this paper summarizes the present state of knowledge of the mechanics of turbulent boundary layers, and of the methods now being used for fundamental studies of the turbulent fluctuations in turbulent boundary layers. A brief review is given of the semi-empirical method of approach as developed by Buri, Gruschwitz, Fediaevsky, and Kalikhman. In recent years the National Advisory.Commsittee for Aeronautics has sponsored a detailed study at the National Bureau of Standards of the turbulent fluctuations in a turbulent boundary layer under adverse pressure gradient sufficient to produce separation. The aims of this investigation and its present status are described.

Dryden, Hugh L↗

NASA NTRS1992

Development of an algebraic turbulence model for analysis of propulsion flows

A simple turbulence model that will be applicable to propulsion flows having both wall bounded and unbounded regions was developed and installed within the PARC Navier-Stokes code by linking two existing algebraic turbulence models. The first is the Modified Mixing Length (MML) model which is optimized for wall bounded flows. The second is the Thomas model, the standard algebraic turbulence model in PARC which has been used to calculate both bounded and unbounded turbulent flows but was optimized for the latter. This paper discusses both models and the method employed to link them into one model (referred to as the MMLT model). The PARC code with the MMLT model was applied to two dimensional turbulent flows over a flat plate and over a backward facing step to validate and optimize the model and to compare its predictions to those obtained with the three turbulence models already available in PARC.

Georgiadis, N. J.↗

NASA NTRS1992

Development of an algebraic turbulence model for analysis of propulsion flows

Georgiadis, N. J.↗

NASA NTRS2000

CVS Decomposition of 3D Homogeneous Turbulence Using Orthogonal Wavelets

This paper compares the filtering used in Coherent Vortex Simulation (CVS) decomposition with an orthogonal wavelet basis, with the Proper Orthogonal Decomposition (POD) or Fourier filtering. Both methods are applied to a field of Direct Numerical Simulation (DNS) data of 3D forced homogeneous isotropic turbulence at microscale Reynolds number R(sub lambda) = 168. We show that, with only 3%N retained modes, CVS filtering separates the coherent vortex tubes from the incoherent background flow. The latter is structureless, has an equipartition energy spectrum, and has a Gaussian velocity probability distribution function (PDF) and an exponential vorticity PDF. On the other hand, the Fourier basis does not extract the coherent vortex tubes cleanly and leaves organized structures in the residual high wavenumber modes whose PDFs are stretched exponentials for both the velocity and the vorticity.

Farge, Marie↗