Assessment of Scale-Resolving Turbulence Modeling Approaches for Axisymmetric Supersonic Wake Prediction.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
Abstract not provided.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.
Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.
This paper presents a numerical study of the effects of expansion ratio on two-dimensional separating and reattaching flows in plane backward-facing step geometries with parallel walls. Closure of the Reynolds equations was achieved by four different turbulence models: k-epsilon, 'modified' k-epsilon, algebraic stress, and 'modified' algebraic stress models. The k-epsilon model relates the Reynolds stresses to the mean rate of strain through the definition of an isotropic turbulent viscosity. The more advanced algebraic stress model calculates the stresses from implicit algebraic relationships containing the stresses themselves, the mean rate of strain, and the turbulent kinetic energy and its dissipation rate. 'Modified' versions of the models employ a new dissipation rate equation whose production term was made more sensitive to streamwise curvature effects. A new nonequilibrium wall function treatment proposed by Chieng and Launder (1980) was also incorporated into each model.
Model for statistically isotropic homogeneous turbulence in incompressible fluid, representing turbulence as superposition of individual vortex sheets
The background information in establishing several proposed atmospheric turbulence models for use on motion based aircraft simulators was documented. A specific model was proposed which, in addition to varying turbulence intensity (rms velocity), varies the atmospheric turbulence scale length to achieve compatibility with real atmospheric turbulence. With a suitable combination of scale length and intensity distribution, the model will simulate various atmospheric conditions characterized by altitude, stability, and terrain. The model is mechanized to be included in a flight simulator experiment in order to determine to what extent the pilots are sensitive to changes in atmospheric conditions and the realism of the model. The following topics were covered: literature survey, presently used techniques, proposed model, and simulation details.
The turbulent kinetic energy transport equation, which is the basis of turbulence models, is investigated for homogeneous, compressible turbulence using direct numerical simulations performed at CTR. It is shown that the partition between dilatational and solenoidal modes is very sensitive to initial conditions for isotropic decaying turbulence but not for sheared flows. The importance of the dilatational dissipation and of the pressure-dilatation term is evidenced from simulations and a transport equation is proposed to evaluate the pressure-dilatation term evolution. This transport equation seems to work well for sheared flows but does not account for initial condition sensitivity in isotropic decay. An improved model is proposed.
The defects of turbulence models are summarized and the importance of so-called nongradient diffusion in turbulent fluxes is discussed. The mathematical theory of the flow of two interpenetrating continua is reviewed, and the mathematical formulation of the two fluid model is outlined. Results from plane wake, axisymmetric jet, and combustion studies are shown.
Two-dimensional, integral boundary-layer methods of the entrainment and dissipation type, with closure relationships which account for turbulent separated flows through the velocity and turbulence models, are developed. Kinematic shape factor relations used in the two methods are evaluated by comparison with recently obtained experimental data for separated flows. Prescribed-displacement-thickness computations are made with both methods and compared to experiment for a series of separated flows, including low speed and transonic shock-induced separation. Effects of upstream history on the turbulent shear stress development are shown to be important for the accurate prediction of the separated flows considered and the experimental trends are reasonably modeled by a rate equation derived from the turbulent kinetic energy equation.
The atmospheric turbulence velocities seen by nonrotating aircraft components and rotating blades can be substantially different. The differences are due to the spatial motion of the rotor blades, which move fore and aft through the gust waves. Body-fixed atmospheric turbulence refers to the actual atmospheric turbulence experienced by a point fixed on a nonrotating aircraft component such as the aircraft's center of gravity or the rotor hub, while blade-fixed atmospheric turbulence refers to the atmospheric turbulence experienced by an element of the rotating rotor blade. An example is presented, which, though overly simplified, shows important differences between blade- and body-fixed rotorcraft atmospheric turbulence models. All of the information necessary to develop the dynamic equations describing the atmospheric turbulence velocity field experienced by an aircraft is contained in the atmospheric turbulence velocity correlation matrix. It is for this reason that a generalized formulation of the correlation matrix describing atmospheric turbulence that a rotating blade encounters is developed. From this correlation matrix, earlier treated cases restricted to a rotor flying straight and level directly into the mean wind can be recovered as special cases.
This article details an effort to improve the understanding and prediction of turbulent flow inside a droplet of molten metal levitated in an electromagnetic field. It is shown that the flow field in a test case, a nickel droplet levitated under microgravity conditions, is in the transitional regime between laminar and turbulent flow. Past research efforts have used laminar, enhanced viscosity, and k-epsilon turbulence models to describe the flow. The method highlighted in our study is the renormalization group (RNG) algorithm. We show that an accurate description of the turbulent eddy viscosity is critical in order to obtain realistic velocity fields, and that the turbulent eddy viscosity cannot be uniform in levitated droplets. The RNG method does not impose isotropic length or time scales on the flow field, thus allowing such nonuniform features to be captured. A number of other materials processing applications exhibit similarly complex flow characteristics, such as highly recirculating, transitional, and free surface flows, for which this modeling approach may prove useful.
The paper considers the mixing and subsequent combustion within turbulent reacting shear layers. A computer program was developed to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method integrates the Reynolds averaged Navier-Stokes equations using a finite volume approach while advancing the solution forward in time using a Runge-Kutta scheme. Three separate flowfields are investigated and it is found that no single turbulence model considered could accurately predict the degree of mixing for all three cases.
Boundary layer mass removal (bleed) through spanwise bands of holes on a surface is used to prevent or control separation and to stabilize the normal shock in supersonic inlets. The addition of a transport equation lag relationship for eddy viscosity to the rough wall algebraic turbulence model of Cebeci and Chang was found to improve agreement between predicted and measured mean velocity distributions downstream of a bleed band. The model was demonstrated for a range of bleed configurations, bleed rates, and local free stream Mach numbers. In addition, the model was applied to the boundary layer development over acoustic lining materials for the inlets and nozzles of commercial aircraft. The model was found to yield accurate results for integral boundary layer properties unless there was a strong adverse pressure gradient.
JeNo (Version 1.0) is a Fortran90 computer code that calculates the far-field sound spectral density produced by axisymmetric, unheated jets at a user specified observer location and frequency range. The user must provide a structured computational grid and a mean flow solution from a Reynolds-Averaged Navier Stokes (RANS) code as input. Turbulence kinetic energy and its dissipation rate from a k-epsilon or k-omega turbulence model must also be provided. JeNo is a research code, and as such, its development is ongoing. The goal is to create a code that is able to accurately compute far-field sound pressure levels for jets at all observer angles and all operating conditions. In order to achieve this goal, current theories must be combined with the best practices in numerical modeling, all of which must be validated by experiment. Since the acoustic predictions from JeNo are based on the mean flow solutions from a RANS code, quality predictions depend on accurate aerodynamic input.This is why acoustic source modeling, turbulence modeling, together with the development of advanced measurement systems are the leading areas of research in jet noise research at NASA Glenn Research Center.
A new model dissipation rate equation is proposed based on the dynamic equation of the mean-square vorticity fluctuation for large Reynolds number turbulence. The advantage of working with the vorticity fluctuation equation is that the physical meanings of the terms in this equation are more clear than those in the dissipation rate equation. Hence, the model development based on the vorticity fluctuation equation is more straightforward. The resulting form of the model equation is consistent with the spectral energy cascade analysis introduced by Lumley. The proposed model dissipation rate equation is numerically well behaved and can be applied to any level of turbulence modeling. It is applied to a realizable eddy viscosity model. Flows that are examined include: rotating homogeneous shear flows; free shear flows; a channel flow and flat plate boundary layers with and without pressure gradients; and backward facing step separated flows. In most cases, the present model predictions show considerable improvement over the standard kappa-epsilon model.