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 883 records · Page 49

CMOTT Turbulence Module for NPARC

This is a user's manual of the CMOTT turbulence module, version 2.0, developed for the NPARC code. The module is written in a self-contained manner so that the user can use any turbulence model in the module without concern as to how it is implemented and solved. Three two-equation turbulence models have been built into the module: Chien, Shih-Lumley and CMOTT models, and all of them have both the low Reynolds number and wall function options. Unlike Chien's model, both the Shih-Lumley and CMOTT models do not involve the dimensionless wall distance y(sup +) in the low Reynolds number approach, an advantage for separated flow calculations. The Van Driest transformation is used so that the wall functions can be applied to both incompressible and compressible flows. The manual gives the details of the turbulence models used and their numerical implementation. It also gives two application examples, one for subsonic and the other for transonic flow, for demonstration. The module can be easily linked to the NPARC code for practical applications.

Zhu, J.↗

Testing and Implementation of Advanced Reynolds Stress Models

A research program was proposed for the testing and implementation of advanced turbulence models for non-equilibrium turbulent flows of aerodynamic importance that are of interest to NASA. Turbulence models that are being developed in connection with the Office of Naval Research ARI in Non-equilibrium are provided for implementation and testing in aerodynamic flows at NASA Langley Research Center. Close interactions were established with researchers at Nasa Langley RC and refinements to the models were made based on the results of these tests. The models that have been considered include two-equation models with an anisotropic eddy viscosity as well as full second-order closures. Three types of non-equilibrium corrections to the models have been considered in connection with the ARI on Nonequilibrium Turbulence: conducted for ONR.

Speziale, Charles G.↗

Validation of the NPARC code for nozzle afterbody flows at transonic speeds

The NPARC code, a Reynolds-averaged full Navier-Stokes code, was validated for nozzle afterbody (boatail) flow fields at transonic speeds. The flow fields about three geometries were studied: an axisymmetric nozzle with attached flow; an axisymmetric nozzle with separated flow: and a two-dimensional (rectangular) nozzle with separated flow. Three turbulence models, Baldwin-Lomax, Baldwin-Barth, and Chien k-epsilon, were used to determine the effect of turbulence model selection on the flow field solution. Static pressure distributions on the nozzle surfaces and pitot pressure measurements in the exhaust plume were examined. Results from the NPARC code compared very well with experimental data for all cases. For attached flow fields, the effect of the turbulence models showed no discernable differences. The Baldwin-Barth model yielded better results than either the Chien k-epsilon or the Baldwin-Lomax model for separated flow fields.

Debonis, James R.↗

Flame-Vortex Interactions in Microgravity to Improve Models of Turbulent Combustion

A unique flame-vortex interaction experiment is being operated in microgravity in order to obtain fundamental data to assess the Theory of Flame Stretch which will be used to improve models of turbulent combustion. The experiment provides visual images of the physical process by which an individual eddy in a turbulent flow increases the flame surface area, changes the local flame propagation speed, and can extinguish the reaction. The high quality microgravity images provide benchmark data that are free from buoyancy effects. Results are used to assess Direct Numerical Simulations of Dr. K. Kailasanath at NRL, which were run for the same conditions.

Driscoll, James F.↗

A One‐Dimensional Model for Turbulent Mixing in the Benthic Biolayer of Stream and Coastal Sediments

Abstract In this paper, we develop and validate a rigorous modeling framework, based on Duhamel's Theorem, for the unsteady one‐dimensional vertical transport of a solute across a flat sediment‐water interface (SWI) and through the benthic biolayer of a turbulent stream. The modeling framework is novel in capturing the two‐way coupling between evolving solute concentrations above and below the SWI and in allowing for a depth‐varying diffusivity. Three diffusivity profiles within the sediment (constant, exponentially decaying, and a hybrid model) are evaluated against an extensive set of previously published laboratory measurements of turbulent mass transfer across the SWI. The exponential diffusivity profile best represents experimental observations and its reference diffusivity scales with the permeability Reynolds number, a dimensionless measure of turbulence at the SWI. The depth over which turbulence‐enhanced diffusivity decays is of the order of centimeters and comparable to the thickness of the benthic biolayer. Thus, turbulent mixing across the SWI may serve as a universal transport mechanism, supplying the nutrient and energy fluxes needed to sustain microbial growth, and nutrient processing, in the benthic biolayer of stream and coastal sediments.

Grant, Stanley B.↗

Reynolds number and pressure gradient effects on compressible turbulent boundary layers

A detailed investigation of attached supersonic turbulent boundary layers over an extensive range of Reynolds numbers (12 x 10 to the 6th to 314 x 10 to the 6th) is presented. Experimental measurements were obtained for adverse pressure gradients ranging in magnitude from those of previous investigations to those approaching separation. The measurements include mean values of surface pressure and skin-friction, mean-flow profiles, and profiles of the three turbulent velocity fluctuation components and turbulent shear stress. Numerical solutions, employing three turbulence models of various degrees of complexity have been compared with the details of the measured flow fields. Generally, it was found that the more sophisticated turbulence models are superior to a mixing length model for predicting the Reynolds number and pressure gradient effects. However, some details of the turbulent fluctuations as well as the exact Reynolds number trends indicated by the data were not accurately predicted with any of the turbulence models considered.

Acharya, M.↗

Parametrics on 2D Navier-Stokes analysis of a Mach 2.68 bifurcated rectangular mixed-compression inlet

The supersonic diffuser of a Mach 2.68 bifurcated, rectangular, mixed-compression inlet was analyzed using a two-dimensional (2D) Navier-Stokes flow solver. Parametric studies were performed on turbulence models, computational grids and bleed models. The computer flowfield was substantially different from the original inviscid design, due to interactions of shocks, boundary layers, and bleed. Good agreement with experimental data was obtained in many aspects. Many of the discrepancies were thought to originate primarily from 3D effects. Therefore, a balance should be struck between expending resources on a high fidelity 2D simulation, and the inherent limitations of 2D analysis. The solutions were fairly insensitive to turbulence models, grids and bleed models. Overall, the k-e turbulence model, and the bleed models based on unchoked bleed hole discharge coefficients or uniform velocity are recommended. The 2D Navier-Stokes methods appear to be a useful tool for the design and analysis of supersonic inlets, by providing a higher fidelity simulation of the inlet flowfield than inviscid methods, in a reasonable turnaround time.

Mizukami, M.↗

Assessment of Turbulent Shock-Boundary Layer Interaction Computations Using the OVERFLOW Code

The performance of two popular turbulence models, the Spalart-Allmaras model and Menter s SST model, and one relatively new model, Olsen & Coakley s Lag model, are evaluated using the OVERFLOWcode. Turbulent shock-boundary layer interaction predictions are evaluated with three different experimental datasets: a series of 2D compression ramps at Mach 2.87, a series of 2D compression ramps at Mach 2.94, and an axisymmetric coneflare at Mach 11. The experimental datasets include flows with no separation, moderate separation, and significant separation, and use several different experimental measurement techniques (including laser doppler velocimetry (LDV), pitot-probe measurement, inclined hot-wire probe measurement, preston tube skin friction measurement, and surface pressure measurement). Additionally, the OVERFLOW solutions are compared to the solutions of a second CFD code, DPLR. The predictions for weak shock-boundary layer interactions are in reasonable agreement with the experimental data. For strong shock-boundary layer interactions, all of the turbulence models overpredict the separation size and fail to predict the correct skin friction recovery distribution. In most cases, surface pressure predictions show too much upstream influence, however including the tunnel side-wall boundary layers in the computation improves the separation predictions.

Oliver, A. B.↗

Computation of Flow Over a Drag Prediction Workshop Wing/Body Transport Configuration Using CFL3D

A Drag Prediction Workshop was held in conjunction with the 19th AIAA Applied Aerodynamics Conference in June 2001. The purpose of the workshop was to assess the prediction of drag by computational methods for a wing/body configuration (DLR-F4) representative of subsonic transport aircraft. This report details computed results submitted to this workshop using the Reynolds-averaged Navier-Stokes code CFL3D. Two supplied grids were used: a point-matched 1-to-1 multi-block grid, and an overset multi-block grid. The 1-to-1 grid, generally of much poorer quality and with less streamwise resolution than the overset grid, is found to be too coarse to adequately resolve the surface pressures. However, the global forces and moments are nonetheless similar to those computed using the overset grid. The effect of three different turbulence models is assessed using the 1-to-1 grid. Surface pressures are very similar overall, and the drag variation due to turbulence model is 18 drag counts. Most of this drag variation is in the friction component, and is attributed in part to insufficient grid resolution of the 1-to-1 grid. The misnomer of 'fully turbulent' computations is discussed; comparisons are made using different transition locations and their effects on the global forces and moments are quantified. Finally, the effect of two different versions of a widely used one-equation turbulence model is explored.

Rumsey, Christopher L.↗

Summary of Data from the Fifth AIAA CFD Drag Prediction Workshop

Results from the Fifth AIAA CFD Drag Prediction Workshop (DPW-V) are presented. As with past workshops, numerical calculations are performed using industry-relevant geometry, methodology, and test cases. This workshop focused on force/moment predictions for the NASA Common Research Model wing-body configuration, including a grid refinement study and an optional buffet study. The grid refinement study used a common grid sequence derived from a multiblock topology structured grid. Six levels of refinement were created resulting in grids ranging from 0.64x10(exp 6) to 138x10(exp 6) hexahedra - a much larger range than is typically seen. The grids were then transformed into structured overset and hexahedral, prismatic, tetrahedral, and hybrid unstructured formats all using the same basic cloud of points. This unique collection of grids was designed to isolate the effects of grid type and solution algorithm by using identical point distributions. This study showed reduced scatter and standard deviation from previous workshops. The second test case studied buffet onset at M=0.85 using the Medium grid (5.1x106 nodes) from the above described sequence. The prescribed alpha sweep used finely spaced intervals through the zone where wing separation was expected to begin. Some solutions exhibited a large side of body separation bubble that was not observed in the wind tunnel results. An optional third case used three sets of geometry, grids, and conditions from the Turbulence Model Resource website prepared by the Turbulence Model Benchmarking Working Group. These simple cases were intended to help identify potential differences in turbulence model implementation. Although a few outliers and issues affecting consistency were identified, the majority of participants produced consistent results.

Levy, David W.↗

The present state and future directions of PDF methods

The objectives of the workshop are presented in viewgraph format, as is this entire article. The objectives are to discuss the present status and the future direction of various levels of engineering turbulence modeling related to Computational Fluid Dynamics (CFD) computations for propulsion; to assure that combustion is an essential part of propulsion; and to discuss Probability Density Function (PDF) methods for turbulent combustion. Essential to the integration of turbulent combustion models is the development of turbulent model, chemical kinetics, and numerical method. Some turbulent combustion models typically used in industry are the k-epsilon turbulent model, the equilibrium/mixing limited combustion, and the finite volume codes.

Pope, S. B.↗

A diagonally inverted LU implicit multigrid scheme for the 3-D Navier-Stokes equations and a two equation model of turbulence

An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows.

Yokota, Jeffrey W.↗

A diagonally inverted LU implicit multigrid scheme for the 3-D Navier-Stokes equations and a two equation model of turbulence

An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows.

Yokota, Jeffrey W.↗

Multivalued velocity field model of turbulence

The considered study is based on new theoretical concepts regarding a post-instability model of a fluid discussed by Zak (1980). The model permits the completion of the governing equations of turbulence by introducing multivalued fields of velocities. Attention is given to the mechanism of energy dissipation, the characteristic wave propagation, a simplified model, the formation of turbulence around stagnation points, the formulation of boundary conditions, and the mechanism of turbulence formation. The mechanism of turbulence formation can be understood as propagation of initial discontinuities from the boundaries into a flow with the characteristic velocity which is defined by the normal (to the boundary) velocity components. These components emerge at the boundary as a result of jumps in the tangential components due to the continuity equation.

Zak, M.↗

Navier-Stokes computation of compressible turbulent flows with a second order closure, part 1

A second order closure turbulence model for compressible flows is developed and implemented in a 2D Reynolds-averaged Navier-Stokes solver. From the beginning where a kappa-epsilon turbulence model was implemented in the bidiagonal implicit method of MACCORMACK (referred to as the MAC3 code) to the final stage of implementing a full second order closure in the efficient line Gauss-Seidel algorithm, numerous work was done, individually and collectively. Besides the collaboration itself, the final product of this work is a second order closure derived from the Launder, Reece, and Rodi model to account for near wall effects, which has been called FRAME model, which stands for FRench-AMerican-Effort. During the reporting period, two different problems were worked out. The first was to provide Ames researchers with a reliable compressible boundary layer code including a wide collection of turbulence models for quick testing of new terms, both in two equations and in second order closure (LRR and FRAME). The second topic was to complete the implementation of the FRAME model in the MAC5 code. The work related to these two different contributions is reported. dilatation in presence of stron shocks. This work, which has been conducted during a work at the Center for Turbulence Research with Zeman aimed also to cros-check earlier assumptions by Rubesin and Vandromme.

Haminh, Hieu↗

Segregation parameters and pair-exchange mixing models for turbulent nonpremixed flames

The progress of chemical reactions in nonpremixed turbulent flows depends on the coexistence of reactants, which are brought together by mixing. The degree of mixing can strongly influence the chemical reactions and it can be quantified by segregation parameters. In this paper, the relevance of segregation parameters to turbulent mixing and chemical reactions is explored. An analysis of the pair-exchange mixing models is performed and an explanation is given for the peculiar behavior of such models in homogeneous turbulence. The nature of segregation parameters in a H2/Ar-air nonpremixed jet flame is investigated. The results show that Monte Carlo simulation with the modified Curl's mixing model predicts segregation parameters in close agreement with the experimental values, providing an indirect validation for the theoretical model.

Chen, J.-Y.↗

The use of transition region characteristics to improve the numerical simulation of heat transfer in bypass transitional flows

A method is presented for improving the numerical prediction of bypass transition heat transfer on a flat plate in a high-disturbance environment with zero or favorable pressure gradient. The method utilizes low Reynolds number k-epsilon turbulence models in combination with the characteristic parameters of the transition region. The parameters representing the characteristics of the transition region used are the intermittency, transition length and turbulent spot properties. An analysis is made of the transition length in terms of turbulent spot variables. The nondimensional spot formation rate, required for the prediction of the transition length, is shown by the analysis to be a function of the spot spreading angle, the dimensionless spot velocity ratio and the dimensionless spot area ratio. The intermittency form of the k-epsilon equations were derived from conditionally averaged equations which have been shown to be an improvement over global-time-averaged equations for the numerical calculation of the transition region. The numerical predictions are in general good agreement with the experimental data and indicate the potential use of the method in accelerating flows. Turbulence models of the k-epsilon type are known to underpredict the transition length. The present work demonstrates how incorporating transition region characteristics improves the ability of two-equation turbulence models to simulate bypass transition for flat plates with potential application to turbine vanes and blades.

Simon, Frederick F.↗

Further improvements of a new model for turbulent convection in stars

The effects of including a variable molecular weight and of using the newest opacities of Rogers and Iglesias (1991) as inputs to a recent model by Canuto and Mazzitelli (1991) for stellar turbulent convection are studied. Solar evolutionary tracks are used to conclude that the the original model for turbulence with mixing length Lambda = z, Giuli's variable Q unequal to 1 and the new opacities yields a fit to solar T(eff) within 0.5 percent. A formulation of Lambda is proposed that extends the purely nonlocal Lambda = z expression to include local effects. A new expression for Lambda is obtained which generalizes both the mixing length theory (MLT) phenomenological expression for Lambda as well as the model Lambda = z. It is argued that the MLT should now be abandoned.

Canuto, V. M.↗