Engineering Papers⌕ Search

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.

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

At least 307 records · Page 17

Numerical Modeling of Turbulent Combustion

The work in numerical modeling is focused on the use of the random vortex method to treat turbulent flow fields associated with combustion while flame fronts are considered as interfaces between reactants and products, propagating with the flow and at the same time advancing in the direction normal to themselves at a prescribed burning speed. The latter is associated with the generation of specific volume (the flame front acting, in effect, as the locus of volumetric sources) to account for the expansion of the flow field due to the exothermicity of the combustion process. The model was applied to the flow in a channel equipped with a rearward facing step. The results obtained revealed the mechanism of the formation of large scale turbulent structure in the wake of the step, while it showed the flame to stabilize on the outer edges of these eddies.

Ghoneim, A. F.↗

Equations and closure methods

Basic differential equations governing compressible turbulent boundary layer flow are reviewed, including conservation of mass and energy, momentum equations derived from Navier-Stokes equations, and equations of state. Closure procedures were broken down into: (1) simple or zeroth-order methods, (2) first-order or mean field closure methods, and (3) second-order or mean turbulence field methods.

Source record↗

Direct simulations of turbulent flow using finite-difference schemes

A high-order accurate finite-difference approach is presented for calculating incompressible turbulent flow. The methods used include a kinetic energy conserving central difference scheme and an upwind difference scheme. The methods are evaluated in test cases for the evolution of small-amplitude disturbances and fully developed turbulent channel flow. It is suggested that the finite-difference approach can be applied to complex geometries more easilty than highly accurate spectral methods. It is concluded that the upwind scheme is a good candidate for direct simulations of turbulent flows over complex geometries.

Rai, Man Mohan↗

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗

Upgraded viscous flow analysis of multielement airfoils

A description of an improved version of the NASA/Lockheed multielement airfoil analysis computer program is presented. The improvements include several major modifications of the aerodynamic model as well as substantial changes of the computer code. The modifications of the aerodynamic model comprise the representation of the boundary layer and wake displacement effects with an equivalent source distribution, the prediction of wake parameters with Green's lag-entrainment method, the calculation of turbulent boundary layer separation with the method of Nash and Hicks, the estimation of the onset of confluent boundary layer separation with a modified form of Goradia's method, and the prediction of profile drag with the formula of Squire and Young. The paper further describes the modifications of the computer program for which the structured approach to computer software development was employed. Important aspects of the structured program development such as the functional decomposition of the aerodynamic theory and its numerical implementation, the analysis of the data flow within the code, and the application of a pseudo code are discussed.

Brune, G. W.↗

Upgraded viscous flow analysis of multi-element airfoils

A description of an improved version of the NASA/Lockheed multi-element airfoil analysis computer program is presented. The improvements include several major modifications of the aerodynamic model as well as substantial changes of the computer code. The modifications of the aerodynamic model comprise the representation of the boundary layer and wake displacement effects with an equivalent source distribution, the prediction of wake parameters with Green's lag-entrainment method, the calculation of turbulent boundary layer separation with the method of Nash and Hicks, the estimation of the onset of confluent boundary layer separation with a modified form of Goradia's method, and the prediction of profile drag with the formula of Squire and Young. The modifications of the computer program for which the structured approach to computer software development was employed are also described. Important aspects of the structured program development such as the functional decomposition of the aerodynamic theory and its numerical implementation, the analysis of the data flow within the code, and the application of a pseudo code are discussed.

Brune, G. W.↗

Spectral multigrid methods for the solution of homogeneous turbulence problems

New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.

Erlebacher, G.↗

Spectral multigrid methods for the solution of homogeneous turbulence problems

New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.

Erlebacher, G.↗

Renormalization group analysis of turbulence

The objective is to understand and extend a recent theory of turbulence based on dynamic renormalization group (RNG) techniques. The application of RNG methods to hydrodynamic turbulence was explored most extensively by Yakhot and Orszag (1986). An eddy viscosity was calculated which was consistent with the Kolmogorov inertial range by systematic elimination of the small scales in the flow. Further, assumed smallness of the nonlinear terms in the redefined equations for the large scales results in predictions for important flow constants such as the Kolmogorov constant. It is emphasized that no adjustable parameters are needed. The parameterization of the small scales in a self-consistent manner has important implications for sub-grid modeling.

Smith, Leslie M.↗

Some practical turbulence modeling options for Reynolds-averaged full Navier-Stokes calculations of three-dimensional flows

New turbulence modeling options recently implemented for the 3-D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Validation test cases include the incompressible and compressible flat plate turbulent boundary layers, turbulent developing S-duct flow, and glancing shock wave/turbulent boundary layer interaction. Good agreement is obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(sup +) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. The test cases show that the highly optimized one-and two-equation turbulence models can be used in routine 3-D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.

Bui, Trong T.↗

CFD Sensitivity Analysis of a Modern Civil Transport Near Buffet-Onset Conditions

A computational fluid dynamics (CFD) sensitivity analysis is conducted for a modern civil transport at several conditions ranging from mostly attached flow to flow with substantial separation. Two different Navier-Stokes computer codes and four different turbulence models are utilized, and results are compared both to wind tunnel data at flight Reynolds number and flight data. In-depth CFD sensitivities to grid, code, spatial differencing method, aeroelastic shape, and turbulence model are described for conditions near buffet onset (a condition at which significant separation exists). In summary, given a grid of sufficient density for a given aeroelastic wing shape, the combined approximate error band in CFD at conditions near buffet onset due to code, spatial differencing method, and turbulence model is: 6% in lift, 7% in drag, and 16% in moment. The biggest two contributers to this uncertainty are turbulence model and code. Computed results agree well with wind tunnel surface pressure measurements both for an overspeed 'cruise' case as well as a case with small trailing edge separation. At and beyond buffet onset, computed results agree well over the inner half of the wing, but shock location is predicted too far aft at some of the outboard stations. Lift, drag, and moment curves are predicted in good agreement with experimental results from the wind tunnel.

Rumsey, Christopher L.↗

Recent advances in PDF modeling of turbulent reacting flows

This viewgraph presentation concludes that a Monte Carlo probability density function (PDF) solution successfully couples with an existing finite volume code; PDF solution method applied to turbulent reacting flows shows good agreement with data; and PDF methods must be run on parallel machines for practical use.

Leonard, Andrew D.↗

Scaling and modeling of three-dimensional, end-wall, turbulent boundary layers

The method of matched asymptotic expansion was employed to identify the various subregions in three dimensional, turbomachinery end wall turbulent boundary layers, and to determine the proper scaling of these regions. The two parts of the b.l. investigated are the 3D pressure driven part over the endwall, and the 3D part located at the blade/end wall juncture. Models are proposed for the 3d law of the wall and law of the wake. These models and the data of van den Berg and Elsenaar and of Mueller are compared and show good agreement between models and experiments.

Goldberg, U. C.↗

Indirect measurement of the thermal-acoustic efficiency spectrum of a long turbulent burner

A new method is described for deducing the thermal-acoustic efficiency spectrum (defined as the fraction of combustion heat release converted to acoustic energy at a given frequency) of a long turbulent burner from the sound spectrum measured in the far field. The method, which is based on a one-dimensional model of the unsteady flow in the burner, is applied to a tubular diffusion-flame hydrogen burner whose length is large compared to its diameter. The results for thermal powers ranging from 4.5 to 22.3 kW show that the thermal-acoustic efficiency is relatively insensitive to the burner power level, decreasing from a value of around 0.0001 at 150 Hz with a slope of about 20 dB per decade. Evidence is presented indicating that acoustic agitation of the flame below 500 Hz, especially in the neighborhood of the resonant frequencies of the burner, is a significant acoustic source.

Mahan, J. R.↗