Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence modeling”

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 1,261 records · Page 70

CFD validation for aerodynamic flows - Challenge for the '90's

The process of computational fluid dynamics (CFD) validation is described from two perspectives, numerical error validation and physical model validation. Errors associated with time and space discretization, grid refinement, numerical dissipation, and level of solution convergence, are all possible sources of numerical error. The second aspect of CFD code validation is associated with the governing equations and the physical models, e.g., chemistry and turbulence models. This type of validation requires comparisons with carefully conducted experiments and is the primary aspect of the validation process discussed in this paper. Examples of validation are shown for selected Reynolds-averaged Navier-Stokes codes for conditions from transonic through hypersonic. Since code applications are becoming more complex, it no longer suffices to use data from surface or integral quantities alone to provide the required validation. Flow field surveys and experimental boundary condition measurements are emerging as critical data that must be obtained for CFD code validation. Progress in the overall level of validation is shown to be improving. However, future challenges remain, some of which are enumerated.

Marvin, Joseph G.↗

Symmetry plane model for turbulent flows with vortex generators

An approximate procedure is proposed for predicting the performance of counterrotating vortex-generator installations in incompressible flow. An inviscid calculation that includes the motion of the vortices is used to obtain crossflow velocities at the boundary-layer edge as a function of initial position, spacing, and strength of the vortices, and local values of the spanwise gradient are then folded into an integral turbulent-boundary layer procedure applied in the plane of symmetry. Special attention is paid to the consistency of the approximations and equations used. The two-dimensional aerodynamics of vortex generator installations on a NACA 0016 airfoil at angle-of-attack are estimated in this manner, and the results compared with experiments carried out with a 30-cm chord wing mounted in a 2.4 x 3.6-m cross-section wind tunnel and tested at chord Reynolds numbers of 0.7 and 1.4 x 10 to the 6th. Agreement in the separation location is found for these complex flows for a range of conditions.

Arnaud, Gilles L.↗

Ashfall dispersal for the 16 September 1986, eruption of Lascar, Chile, calculated by a turbulent diffusion model

A two dimensional model for the dispersal of volcanic ash, modified from that of Suzuki (1983), is applied to the September 16, 1986 eruption of Lascar, in northern Chile. The result is an isopach map predicting ground concentrations of ash. The map shows two local maxima along distinct dispersal axes and predicts that the deposit was detached from the source. Predicted concentrations are in agreement with a single reported thickness measurement.

Glaze, Lori S.↗

Unstructured mesh algorithms for aerodynamic calculations

The use of unstructured mesh techniques for solving complex aerodynamic flows is discussed. The principle advantages of unstructured mesh strategies, as they relate to complex geometries, adaptive meshing capabilities, and parallel processing are emphasized. The various aspects required for the efficient and accurate solution of aerodynamic flows are addressed. These include mesh generation, mesh adaptivity, solution algorithms, convergence acceleration, and turbulence modeling. Computations of viscous turbulent two-dimensional flows and inviscid three-dimensional flows about complex configurations are demonstrated. Remaining obstacles and directions for future research are also outlined.

Mavriplis, D. J.↗

Unstructured mesh algorithms for aerodynamic calculations

The use of unstructured mesh techniques for solving complex aerodynamic flows is discussed. The principle advantages of unstructured mesh strategies, as they relate to complex geometries, adaptive meshing capabilities, and parallel processing are emphasized. The various aspects required for the efficient and accurate solution of aerodynamic flows are addressed. These include mesh generation, mesh adaptivity, solution algorithms, convergence acceleration, and turbulence modeling. Computations of viscous turbulent two-dimensional flows and inviscid three-dimensional flows about complex configurations are demonstrated. Remaining obstacles and directions for future research are also outlined.

Mavriplis, D. J.↗

Overview of aerothermodynamic loads definition study

The objective of the Aerothermodynamic Loads Definition Study is to develop methods of accurately predicting the operating environment in advanced Earth-to-Orbit (ETO) propulsion systems, such as the Space Shuttle Main Engine (SSME) powerhead. Development of time averaged and time dependent three dimensional viscous computer codes as well as experimental verification and engine diagnostic testing are considered to be essential in achieving that objective. Time-averaged, nonsteady, and transient operating loads must all be well defined in order to accurately predict powerhead life. Described here is work in unsteady heat flow analysis, improved modeling of preburner flow, turbulence modeling for turbomachinery, computation of three dimensional flow with heat transfer, and unsteady viscous multi-blade row turbine analysis.

Gaugler, Raymond E.↗

CFD Modeling of Turbulent Duct Flows for Coolant Channel Analysis

The design of modern liquid rocket engines requires the analysis of chamber coolant channels to maximize the heat transfer while minimizing the coolant flow. Coolant channels often do not remain at a constant cross section or at uniform curvature. New designs require higher aspect ratio coolant channels than previously used. To broaden the analysis capability and to complement standard analysis tools an investigation on the accuracy of CFD predictions for coolant channel flow has been initiated. Validation of CFD capabilities for coolant channel analysis will enhance the capabilities for optimizing design parameters without resorting to extensive experimental testing. The eventual goal is to use CFD to determine the flow fields of unique coolant channel designs and therefore determine critical heat transfer coefficients. In this presentation the accuracy of a particular CFD code is evaluated for turbulent flows. The first part of the presentation is a comparison of numerical results to existing cold flow data for square curved ducts (NASA CR-3367, 'Measurements of Laminar and Turbulent Flow in a Curved Duct with Thin Inlet Boundary Layers'). The results of this comparison show good agreement with the relatively coarse experimental data. The second part of the presentation compares two cases of higher aspect ratio channels (AR=2.5,10) to show changes in axial and secondary flow strength. These cases match experimental work presently in progress and will be used for future validation. The comparison shows increased secondary flow strength of the higher aspect ratio case due to the change in radius of curvature. The presentation includes a test case with a heated wall to demonstrate the program's capability. The presentation concludes with an outline of the procedure used to validate the CFD code for future design analysis.

Ungewitter, Ronald J.↗

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

Effect of non-Poisson samples on turbulence spectra from laser velocimetry

Spectral analysis of laser velocimetry (LV) data plays an important role in characterizing a turbulent flow and in estimating the associated turbulence scales, which can be helpful in validating theoretical and numerical turbulence models. The determination of turbulence scales is critically dependent on the accuracy of the spectral estimates. Spectral estimations from 'individual realization' laser velocimetry data are typically based on the assumption of a Poisson sampling process. What this Note has demonstrated is that the sampling distribution must be considered before spectral estimates are used to infer turbulence scales.

Sree, Dave↗

Separated and Recovering Turbulent Boundary Layer Flow Behind a Backward Facing Step For Different Reynolds Numbers

Experimental results for a two-dimensional separated turbulent boundary layer behind a backward facing step for five different Reynolds numbers are reported. Results are presented in the form of tables, graphs and a floppy disk for an easy access of the data. Reynolds number based on the step height was varied by changing the reference velocity upstream of the step, U(sub o), and the step height, h. Hot-wire measurement techniques were used to measure three Reynolds stresses and four triple-velocity correlations. In addition, surface pressure and skin friction coefficients were measured. All hot-wire measurements were acquired in a measuring domain which excluded recirculating flow region due to the directional insensitivity of hot-wires. The downstream extent of the domain from the step was 51 h for the largest and I 14h for the smallest step height. This significant downstream length permitted extensive study of the flow recovery. Prediction of perturbed flows and their recovery is particularly attractive for popular turbulence models since variations of turbulence length and time scales and flow interactions in different regions are generally inadequately predicted. The data indicate that the flow in the free shear layer region behaves like the plane mixing layer up to about 2/3 of the mean reattachment length when the flow interaction with the wall commences the flow recovery to that of an ordinary turbulent boundary layer structure. These changes of the flow do not occur abruptly with the change of boundary conditions. A reattachment region represents a transitional region where the flow undergoes the most dramatic adjustments to the new boundary conditions. Large eddies, created in the upstream free-shear layer region, are being torn, recirculated, reentrained back into the main stream interacting with the incoming flow structure. It is foreseeable that it is quite difficult to describe the physics of this region in a rational and quantitative manner other than statistical. Downstream of the reattachment point the flow recovers at different rates near the wall, in the newly developing internal boundary layer, and in the outer part of the flow. It appears that Reynolds stresses do not fully recover up to the longest recovery length of 114 h.

Jovic, Srba↗

A Comparative Study of Three Methodologies for Modeling Dynamic Stall

During the past two decades, there has been an increased reliance on the use of computational fluid dynamics methods for modeling rotors in high speed forward flight. Computational methods are being developed for modeling the shock induced loads on the advancing side, first-principles based modeling of the trailing wake evolution, and for retreating blade stall. The retreating blade dynamic stall problem has received particular attention, because the large variations in lift and pitching moments encountered in dynamic stall can lead to blade vibrations and pitch link fatigue. Restricting to aerodynamics, the numerical prediction of dynamic stall is still a complex and challenging CFD problem, that, even in two dimensions at low speed, gathers the major difficulties of aerodynamics, such as the grid resolution requirements for the viscous phenomena at leading-edge bubbles or in mixing-layers, the bias of the numerical viscosity, and the major difficulties of the physical modeling, such as the turbulence models, the transition models, whose both determinant influences, already present in static maximal-lift or stall computations, are emphasized by the dynamic aspect of the phenomena.

Sankar, L.↗

Modeling of Turbulence Generated Noise in Jets

A numerically calculated Green's function is used to predict jet noise spectrum and its far-field directivity. A linearized form of Lilley's equation governs the non-causal Green s function of interest, with the non-linear terms on the right hand side identified as the source. In this paper, contributions from the so-called self- and shear-noise source terms will be discussed. A Reynolds-averaged Navier-Stokes solution yields the required mean flow as well as time- and length scales of a noise-generating turbulent eddy. A non-compact source, with exponential temporal and spatial functions, is used to describe the turbulence velocity correlation tensors. It is shown that while an exact non-causal Green's function accurately predicts the observed shift in the location of the spectrum peak with angle as well as the angularity of sound at moderate Mach numbers, at high subsonic and supersonic acoustic Mach numbers the polar directivity of radiated sound is not entirely captured by this Green's function. Results presented for Mach 0.5 and 0.9 isothermal jets, as well as a Mach 0.8 hot jet conclude that near the peak radiation angle a different source/Green's function convolution integral may be required in order to capture the peak observed directivity of jet noise.

Khavaran, Abbas↗

Grid Convergence for Turbulent Flows(Invited)

A detailed grid convergence study has been conducted to establish accurate reference solutions corresponding to the one-equation linear eddy-viscosity Spalart-Allmaras turbulence model for two dimensional turbulent flows around the NACA 0012 airfoil and a flat plate. The study involved three widely used codes, CFL3D (NASA), FUN3D (NASA), and TAU (DLR), and families of uniformly refined structured grids that differ in the grid density patterns. Solutions computed by different codes on different grid families appear to converge to the same continuous limit, but exhibit different convergence characteristics. The grid resolution in the vicinity of geometric singularities, such as a sharp trailing edge, is found to be the major factor affecting accuracy and convergence of discrete solutions, more prominent than differences in discretization schemes and/or grid elements. The results reported for these relatively simple turbulent flows demonstrate that CFL3D, FUN3D, and TAU solutions are very accurate on the finest grids used in the study, but even those grids are not sufficient to conclusively establish an asymptotic convergence order.

Diskin, Boris↗

Comparison of High-Fidelity Computational Tools for Wing Design of a Distributed Electric Propulsion Aircraft

A variety of tools, from fundamental to high order, have been used to better understand applications of distributed electric propulsion to aid the wing and propulsion system design of the Leading Edge Asynchronous Propulsion Technology (LEAPTech) project and the X-57 Maxwell airplane. Three high-fidelity, Navier-Stokes computational fluid dynamics codes used during the project with results presented here are FUN3D, STAR-CCM+, and OVERFLOW. These codes employ various turbulence models to predict fully turbulent and transitional flow. Results from these codes are compared for two distributed electric propulsion configurations: the wing tested at NASA Armstrong on the Hybrid-Electric Integrated Systems Testbed truck, and the wing designed for the X-57 Maxwell airplane. Results from these computational tools for the high-lift wing tested on the Hybrid-Electric Integrated Systems Testbed truck and the X-57 high-lift wing presented compare reasonably well. The goal of the X-57 wing and distributed electric propulsion system design achieving or exceeding the required 𝐶 (sub L) = 3.95 for stall speed was confirmed with all of the computational codes.

Deere, Karen A.↗

Convergence Analysis of Turbulent Flow Solutions

Data from the "Turbulence Modeling Resource" website for turbulent flow over an NACA-0012 airfoil is analyzed to determine the convergence behavior of three second-order CFD (Computational Fluid Dynamics) codes: CFL3D (Computational Fluids Lab 3 Dimensional flow solver), FUN3D (Fully Unstructured Navier-stokes flow solver), and TAU (German Aerospace Center (DLR) 2 dimensional code for unstructured hybrid grids solving the Reynolds-Averaged Navier-Stokes equations or the Euler equations). The convergence of both integrated properties and pointwise data are examined. Several different methods for estimating errors and computing convergence rates are compared. A high-order extension to the Richardson extrapolation is developed that improves the accuracy of the mesh limit values and provides a quantitative estimate of the threshold of the asymptotic regime. The coefficient of total drag exhibits second-order convergence for all three codes, and convergence is monotone over a sequence of 7 grids. Other force coefficients are not so well behaved. The convergence rates of the viscous component of drag on the three nest grids ranges from 3:0 for CFL3D to 1:0 for FUN3D. The three codes are converging to similar but not identical solutions. The largest differences between the codes are in the coefficient of lift for which the difference between CFL3D and FUN3D is greater than 10 (sup minus 4). The best agreement occurs in the viscous component of drag, which is the only force component for which all three codes are converging towards each other at a rate of second-order. The agreement between the two unstructured grid codes is good with all properties except lift converging towards common values at a rate of second-order. No one code was universally better than the other. The TAU code has the lowest error in total drag, FUN3D has the lowest error in lift, and CFL3D has the lowest error in the viscous component of drag.

Atkins, Harold L.↗