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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Internal Passage Heat Transfer Prediction Using Multiblock Grids and a Kappa-Omega Turbulence Model

Numerical simulations of the three-dimensional flow and heat transfer in a rectangular duct with a 180 C bend were performed. Results are presented for Reynolds numbers of 17,000 and 37,000 and for aspect ratios of 0.5 and I.O. A kappa-omega turbulence model with no reference to distance to a wall is used. Direct comparison between single block and multiblock grid calculations are made. Heat transfer and velocity distributions are compared to available literature with good agreement. The multi-block grid system is seen to produce more accurate results compared to a single-block grid with the same number of cells.

Rigby, David L.↗
NASA NTRS1992

Numerical simulation of turbomachinery flows with advanced turbulence models

A three dimensional full Navier-Stokes (FNS) code is used to simulate complex turbomachinery flows. The code incorporates an explicit multistep scheme and solves a conservative form of the density averaged continuity, momentum, and energy equations. A compressible low Reynolds number form of the k-epsilon turbulence model, and a q-omega model and an algebraic Reynolds stress model have been incorporated in a fully coupled manner to approximate Reynolds stresses. The code is used to predict the viscous flow field in a backswept transonic centrifugal compressor for which laser two focus data is available. The code is also used to simulate the tip clearance flow in a cascade. The code has been extended to include unsteady Euler solutions for predicting the unsteady flow through a cascade due to incoming wakes, simulating rotor-stator interactions.

Lakshminarayana, B.↗
NASA NTRS1993

Supersonic flow calculation using a Reynolds-stress and an eddy thermal diffusivity turbulence model

A second-order model for the velocity field and a two-equation model for the temperature field are used to calculate supersonic boundary layers assuming negligible real gas effects. The modeled equations are formulated on the basis of an incompressible assumption and then extended to supersonic flows by invoking Morkovin's hypothesis, which proposes that compressibility effects are completely accounted for by mean density variations alone. In order to calculate the near-wall flow accurately, correction functions are proposed to render the modeled equations asymptotically consistent with the behavior of the exact equations near a wall and, at the same time, display the proper dependence on the molecular Prandtl number. Thus formulated, the near-wall second order turbulence model for heat transfer is applicable to supersonic flows with different Prandtl numbers. The model is validated against flows with different Prandtl numbers and supersonic flows with free-stream Mach numbers as high as 10 and wall temperature ratios as low as 0.3. Among the flow cases considered, the momentum thickness Reynolds number varies from approximately 4,000 to approximately 21,000. Good correlation with measurements of mean velocity, temperature, and its variance is obtained. Discernible improvements in the law-of-the-wall are observed, especially in the range where the big-law applies.

Sommer, T. P.↗
NASA NTRS1982

The upper atmosphere of Uranus - A critical test of isotropic turbulence models

Observations of the August 15, 1980, Uranus occultation of KM 12, obtained from Cerro Tololo InterAmerican Observatory, European Southern Observatory, and Cerro Las Campanas Observatory, are used to compare the atmospheric structure at points separated by approximately 140 km along the planetary limb. The results reveal striking, but by no means perfect correlation of the light curves, ruling out isotropic turbulence as the cause of the light curve spikes. The atmosphere is strongly layered, and any acceptable turbulence model must accommodate the axial ratios of greater than about 60 which are observed. The mean temperature of the atmosphere is 150 plus or minus 15 K for the region near number density 10 to the 14th per cu cm. Derived temperature variations of vertical scale approximately 130 km and amplitude plus or minus 5 K are in agreement for all stations, and correlated spikes correspond to low-amplitude temperature variations with a vertical scale of several kilometers.

French, R. G.↗
NASA NTRS1991

Elliptic flow computation by low Reynolds number two-equation turbulence models

A detailed comparison of ten low-Reynolds-number k-epsilon models is carried out. The flow solver, based on an implicit approximate factorization method, is designed for incompressible, steady two-dimensional flows. The conservation of mass is enforced by the artificial compressibility approach and the computational domain is discretized using centered finite differences. The turbulence model predictions of the flow past a hill are compared with experiments at Re = 10 exp 6. The effects of the grid spacing together with the numerical efficiency of the various formulations are investigated. The results show that the models provide a satisfactory prediction of the flow field in the presence of a favorable pressure gradient, while the accuracy rapidly deteriorates when a strong adverse pressure gradient is encountered. A newly proposed model form that does not explicitly depend on the wall distance seems promising for application to complex geometries.

Michelassi, V.↗
NASA NTRS1980

Simple turbulence models and their application to boundary layer separation

Measurements in the boundary layer and wake of a stalled airfoil are presented in two coordinate systems, one aligned with the airfoil chord, the other being conventional boundary layer coordinates. The NACA 4412 airfoil is studied at a single angle of attack corresponding to maximum lift, the Reynolds number based on chord being 1.5 x 10 to the 6th power. Turbulent boundary layer separation occurred at the 85 percent chord position. The two-dimensionality of the flow was documented and the momentum integral equation studied to illustrate the importance of turbulence contributions as separation is approached. The assumptions of simple eddy-viscosity and mixing-length turbulence models are checked directly against experiment. Curvature effects are found to be important as separation is approached.

Wadcock, A. J.↗
DOE OSTI2022

Dynamic calibration of differential equations using machine learning, with application to turbulence models

We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Verification Test Suite for Spalart-Allmaras QCR2000 Turbulence Model

The paper presents three benchmark cases for verification of Reynolds-averaged Navier-Stokes solvers. The verification studies focus on a one-equation Spalart-Allmaras model, SA-[neg]-QCR2000, that uses a version of quadratic constitutive relations.The benchmark cases are a two-dimensional subsonic flow around a Joukowski airfoil, a three-dimensional supersonic flow through a square duct, and a three-dimensional flow over a wing-fuselage configuration. The turbulence-model formulation, geometry, flow conditions, grids, and the expected output are described in detail. Reference solutions computed by several established codes are shown

Computational Aerodynamics↗
NASA NTRS1991

Computation of supersonic and low subsonic cascade flows using an explicit Navier-Stokes technique and the kappa-epsilon turbulence model

A fully explicit two-dimensional flow solver, based on a four-stage Runge-Kutta scheme, was developed and utilized to predict two-dimensional viscous flow through turbomachinery cascades for which experimental data is available. The formulation is applied to the density averaged Navier-Stokes equations. Several features of the technique improve the ability of the code to predict high Reynolds number flows on highly stretched grids. These include a low Reynolds number compressible form of the k-epsilon turbulence model, anisotropic scaling of artificial dissipation terms and locally varying timestep evaluation based on hyperbolic and parabolic stability considerations. Comparisons between computation and experiment are presented for both a supersonic and a low-subsonic compressor cascade. These results indicate that the code is capable of predicting steady two-dimensional viscous cascade flows over a wide range of Mach numbers in reasonable computational times.

Kunz, R. F.↗
NASA NTRS2010

Effect of Turbulence Modeling on an Excited Jet

The flow dynamics in a high-speed jet are dominated by unsteady turbulent flow structures in the plume. Jet excitation seeks to control these flow structures through the natural instabilities present in the initial shear layer of the jet. Understanding and optimizing the excitation input, for jet noise reduction or plume mixing enhancement, requires many trials that may be done experimentally or computationally at a significant cost savings. Numerical simulations, which model various parts of the unsteady dynamics to reduce the computational expense of the simulation, must adequately capture the unsteady flow dynamics in the excited jet for the results are to be used. Four CFD methods are considered for use in an excited jet problem, including two turbulence models with an Unsteady Reynolds Averaged Navier-Stokes (URANS) solver, one Large Eddy Simulation (LES) solver, and one URANS/LES hybrid method. Each method is used to simulate a simplified excited jet and the results are evaluated based on the flow data, computation time, and numerical stability. The knowledge gained about the effect of turbulence modeling and CFD methods from these basic simulations will guide and assist future three-dimensional (3-D) simulations that will be used to understand and optimize a realistic excited jet for a particular application.

Brown, Clifford A.↗
DOE OSTI2024

Modeling turbulent impurity transport in the SOL of DIII-D with a reduced model

A novel impurity transport model that approximates SOL turbulence as a fluctuating poloidal electric field is shown to be an acceptable replacement for the traditional approach of assigning an arbitrary radial diffusion coefficient to the impurity ions. The model is implemented in the DIVIMP impurity transport code and applied to an L-Mode tungsten divertor experiment on DIII-D. The poloidal electric field is represented as fluctuating between ±1000 V m –1 based on previous measurements. The resulting intermittent v r = E θ × B T transport causes ions to transport both into the core as well as into the far-SOL. Simultaneous agreement with estimates of the W density just inside the separatrix as well as in the far-SOL is obtained (n W ~ 10 14 m –3 and n W ~ 10 12 m –3 , respectively). Prompt re-deposition of the W ions was necessary to obtain agreement (f redep ~ 99%). Here, we conclude that simulating impurity transport using a physics-based approximation for turbulence in the SOL, versus arbitrarily assigning diffusion coefficients, may enable better reactor scale predictions of core impurity contamination.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Implementation of a kappa-epsilon turbulence model to RPLUS3D code

The RPLUS3D code has been developed at the NASA Lewis Research Center to support the National Aerospace Plane (NASP) project. The code has the ability to solve three dimensional flowfields with finite rate combustion of hydrogen and air. The combustion process of the hydrogen-air system are simulated by an 18 reaction path, 8 species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though a fully implicit numerical algorithm is used. A k-epsilon turbulence model has recently been incorporated into the code. Initial validations have been conducted for a flow over a flat plate. Results of the validation studies are shown. Some difficulties in implementing the k-epsilon equations to the code are also discussed.

Chitsomboon, Tawit↗
NASA NTRS1991

Turbulence modeling in supersonic combusting flows

To support the National Aerospace Plane project, the RPLUS3D CFD code has been developed at NASA Lewis. The code has the ability to solve three-dimensional flowfields with finite rate combustion of hydrogen and air. The combustion processes of the hydrogen-air system are simulated by an 18-reaction path, 8-species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though it uses a fully implicit numerical algorithm. A k-epsilon (two equation) turbulence model is incorporated into the RPLUS3D code.

Chitsomboon, Tawit↗
NASA NTRS1989

Modeling of turbulence for compression corner flows and internal flows

A compressible Reynolds averaged Navier-Stokes code is applied to examine the performances of basic turbulence models for unseparated and separated flows. The turbulence models considered are a zero-equation model and several two-equation models, including a new k-epsilon model with an eddy viscosity damping function depending upon the Reynolds number and the distance from the wall. The turbulence models are assessed using the following experimental flows: a channel flow, a backward facing step, and a two-dimensional compression corner at Mach 2.8. Discussions and comparisons of the experimental and numerical results are given. Mesh refinement is shown to have a significant beneficial effect on free shear layer flow results obtained with two-equation turbulence models. The results illustrate the usefulness of several turbulence models as a design tool for fluid engineering systems.

Champney, J. M.↗
NASA NTRS2010

Assessment of Hybrid RANS/LES Turbulence Models for Aeroacoustics Applications

Predicting the noise from aircraft with exposed landing gear remains a challenging problem for the aeroacoustics community. Although computational fluid dynamics (CFD) has shown promise as a technique that could produce high-fidelity flow solutions, generating grids that can resolve the pertinent physics around complex configurations can be very challenging. Structured grids are often impractical for such configurations. Unstructured grids offer a path forward for simulating complex configurations. However, few unstructured grid codes have been thoroughly tested for unsteady flow problems in the manner needed for aeroacoustic prediction. A widely used unstructured grid code, FUN3D, is examined for resolving the near field in unsteady flow problems. Although the ultimate goal is to compute the flow around complex geometries such as the landing gear, simpler problems that include some of the relevant physics, and are easily amenable to the structured grid approaches are used for testing the unstructured grid approach. The test cases chosen for this study correspond to the experimental work on single and tandem cylinders conducted in the Basic Aerodynamic Research Tunnel (BART) and the Quiet Flow Facility (QFF) at NASA Langley Research Center. These configurations offer an excellent opportunity to assess the performance of hybrid RANS/LES turbulence models that transition from RANS in unresolved regions near solid bodies to LES in the outer flow field. Several of these models have been implemented and tested in both structured and unstructured grid codes to evaluate their dependence on the solver and mesh type. Comparison of FUN3D solutions with experimental data and numerical solutions from a structured grid flow solver are found to be encouraging.

Vatsa, Veer N.↗
NASA NTRS1990

One- and two-dimensional Stirling machine simulation using experimentally generated flow turbulence models

Investigations of one- and two-dimensional (1- or 2-D) simulations of Stirling machines centered around experimental data generated by the U. of Minnesota Mechanical Engineering Test Rig (METR) are covered. This rig was used to investigate oscillating flows about a zero mean with emphasis on laminar/turbulent flow transitions in tubes. The Space Power Demonstrator Engine (SPDE) and in particular, its heater, were the subjects of the simulations. The heater was treated as a 1- or 2-D entity in an otherwise 1-D system. The 2-D flow effects impacted the transient flow predictions in the heater itself but did not have a major impact on overall system performance. Information propagation effects may be a significant issue in the simulation (if not the performance) of high-frequency, high-pressure Stirling machines. This was investigated further by comparing a simulation against an experimentally validated analytic solution for the fluid dynamics of a transmission line. The applicability of the pressure-linking algorithm for compressible flows may be limited by characteristic number (defined as flow path information traverses per cycle); this warrants further study. Lastly the METR was simulated in 1- and 2-D. A two-parameter k-w foldback function turbulence model was developed and tested against a limited set of METR experimental data.

Goldberg, Louis F.↗
NASA NTRS1996

Transonic Turbulent Flow Predictions With Two-Equation Turbulence Models

Solutions of the Favre-averaged Navier-Stokes equations for two well-documented transonic turbulent flows are compared in detail with existing experimental data. While the boundary layer in the first case remains attached, a region of extensive flow separation has been observed in the second case. Two recently developed k-epsilon, two-equation, eddy-viscosity models are used to model the turbulence field. These models satisfy the realizability constraints of the Reynolds stresses. Comparisons with the measurements are made for the wall pressure distribution, the mean streamwise velocity profiles, and turbulent quantities. Reasonably good agreement is obtained with the experimental data.

Liou, William W.↗
NASA NTRS1991

Analysis of two-equation turbulence models for recirculating flows

The two-equation kappa-epsilon model is used to analyze turbulent separated flow past a backward-facing step. It is shown that if the model constraints are modified to be consistent with the accepted energy decay rate for isotropic turbulence, the dominant features of the flow field, namely the size of the separation bubble and the streamwise component of the mean velocity, can be accurately predicted. In addition, except in the vicinity of the step, very good predictions for the turbulent shear stress, the wall pressure, and the wall shear stress are obtained. The model is also shown to provide good predictions for the turbulence intensity in the region downstream of the reattachment point. Estimated long time growth rates for the turbulent kinetic energy and dissipation rate of homogeneous shear flow are utilized to develop an optimal set of constants for the two equation kappa-epsilon model. The physical implications of the model performance are also discussed.

Thangam, S.↗