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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Comments on a proposed standard wind hazard environment and its use in real-time aircraft simulations

A set of mean-wind profiles and associated turbulence models have been proposed by the FAA as a standard wind hazard environment to be used in piloted training and research simulators. The 'wind hazard package' is described along with a 'strawman' method of implementation into the real-time simulation of NASA's Terminal Configured Vehicle (TCV). Simulation results of both piloted and automatic landings have raised questions about the validity of the wind hazard environment, and the use of a lumped parameter quasi-steady flow aerodynamic model in this environment. Areas of discussion will include turbulence models, span and area averaging filters, gust penetration effects, and the need for the addition of terms in the moment equations to represent unsteady flow characteristics.

Dieudonne, J. E.↗
NASA NTRS1993

GGOT total pressure loss control concept evaluation

Total pressure loss is one of the most important parameters in the design of a turbine. This parameter effects not only the turbine performance, but consequently the engine power balance and engine performance. Computational Fluid Dynamics (CFD) can be an effective tool in predicting turbine total pressure loss, and also for performing sensitivity studies to achieve an optimal design with respect to pressure loss. In the present study, the AEROVISC code was used to predict the total pressure loss in the Turbine Technology Team Gas Generator Oxidizer Turbine (GGOT). The objectives in this study are two-fold. It is first necessary to determine an optimal methodology in predicting total pressure loss. The type of grid, grid density and distribution are parameters which may affect the loss prediction. Also, the effect of using a standard K-epsilon turbulence model with wall functions versus a two-layer turbulence model needs to be investigated. The use of grid embedding to resolve areas with high flow gradients needs to be explored. The second objective of the study is to apply the optimal methodology toward evaluating different tip leakage control concepts.

Blumenthal, R. F.↗
NASA NTRS2004

Computations of Wall Distances Based on Differential Equations

The use of differential equations such as Eikonal, Hamilton-Jacobi and Poisson for the economical calculation of the nearest wall distance d, which is needed by some turbulence models, is explored. Modifications that could palliate some turbulence-modeling anomalies are also discussed. Economy is of especial value for deforming/adaptive grid problems. For these, ideally, d is repeatedly computed. It is shown that the Eikonal and Hamilton-Jacobi equations can be easy to implement when written in implicit (or iterated) advection and advection-diffusion equation analogous forms, respectively. These, like the Poisson Laplacian term, are commonly occurring in CFD solvers, allowing the re-use of efficient algorithms and code components. The use of the NASA CFL3D CFD program to solve the implicit Eikonal and Hamilton-Jacobi equations is explored. The re-formulated d equations are easy to implement, and are found to have robust convergence. For accurate Eikonal solutions, upwind metric differences are required. The Poisson approach is also found effective, and easiest to implement. Modified distances are not found to affect global outputs such as lift and drag significantly, at least in common situations such as airfoil flows.

Tucker, Paul G.↗
NASA NTRS1984

A convective model for turbulent mixing in rotating convection zones

The effects of rotation are included in an analytical model for the convective motions in a plane-parallel layer of an ideal fluid. The turbulent stress tensor, formed by taking products and averages of the various velocity components, is calculated for an arbitrary eddy size and shape. Heuristic formulae presented for determining the size and shape of the dominant eddy then give a fully specified stress tensor. Applications for this stress tensor in problems of stellar internal dynamics, heat flow, scalar diffusion, and dynamo theory are suggested. The resultant stresses tend to produce differential rotation profiles with rapidly rotating equators and interiors. The dynamo activity associated with these convective motions tends to occur near the lower boundary of the convection zone.

Hathaway, D. H.↗
NASA NTRS1990

On modelling the pressure-strain correlations in wall bounded flows

Turbulence models for the pressure-strain term of the Reynolds-stress equations in the vicinity of a moving wall are evaluated for a high Reynolds number flow using decaying grid turbulence as a model problem. The data of Thomas and Hancock are used as a base for evaluating the different turbulence models. In particular, the Rotta model for return-to-isotropy is evaluated both in its inclusion into the Reynolds-stress equation model and in comparison to a nonlinear model advanced by Sarkar and Speziale. Further, models for the wall correction to the transfer term advanced by Launder et al., Shir, and Shih and Lumley are compared. Initial data using the decaying grid turbulence experiment as a base suggests that the coefficients proposed for these models are high perhaps by as much as an order of magnitude. The Shih and Lumley model which satisfies realizability constraints, in particular, seems to hold promise in adequately modeling the Reynolds stress components of this flow. Extensions of this work are to include testing the homogeneous transfer model by Shih and Lumley and the testing of the wall transfer models using their proposed coefficients and the coefficients chosen from this work in a flow with mean shear component.

Peltier, L. J.↗
NASA NTRS1990

A Monte Carlo simulation technique for low-altitude, wind-shear turbulence

A case is made for including anisotropy in a Monte Carlo flight simulation scheme of low-altitude wind-shear turbulence by means of power spectral density. This study attempts to eliminate all flight simulation-induced deficiencies in the basic turbulence model. A full-scale low-altitude wind-shear turbulence simulation scheme is proposed with particular emphasis on low cost and practicality for near-ground flight. The power spectral density statistic is used to highlight the need for realistic estimates of energy transfer associated with low-altitude wind-shear turbulence. The simulation of a particular anisotropic turbulence model is shown to be a relatively simple extension from that of traditional isotropic (Dryden) turbulence.

Bowles, Roland L.↗
NASA NTRSDate not supplied

Loci/STREAM Sharp Interface Model Validation of Cryogenic Propellant Tank Self-Pressurization

This paper documents a model validation effort for the prediction of cryogenic propellant tank selfpressurization physics using the Loci/STREAM sharp interface CFD model developed at NASA MSFC. Cryogenic fluid management (CFM) applications for long duration spaceflight missions require accurate modeling of heat transfer and thermodynamics. Therefore, detailed validation of those physics is important to ensure simulation accuracy. In the past, the MSFC Propulsion Fluid Dynamics Branch has validated and used the Loci/STREAM Volume of Fluid (VOF) module extensively for propellant slosh applications, which are dominated by fluid motion. For cases where the interface gas/liquid interface is essentially static, the Loci/STREAM sharp interface was developed for more rapid and efficient CFM analysis. In the present study, Loci/STREAM sharp interface simulations were compared with the NASA K-site self-pressurization experiment of a flight weight, partially full, LH2 tank in normal gravity. This study investigated several parametric effects and made a number of important observations on the performance of the Loci/STREAM sharp interface method: - First, a homogeneous model was derived to verify the CFD model. The homogeneous model assumes both gas and liquid phases are uniform and are at saturation temperature. The predicted pressure rise rates at 3.5 W/m2 and 2.0 W/m2 with a 49% fill level are all lower than the experimental measurement. This is expected due to the observed temperature stratification near the gas-liquid interface during testing. - The k-w and k-e turbulence models predict comparable pressurization rates, which are much lower than the experiment measurement. On the other hand, the pressure values and pressurization rates from both models are almost the same as those of the homogeneous model. The current turbulence model is too diffusive and is not capable of capturing the thermal stratification in the liquid and ullage. - When using the laminar flow model, the pressurization rate is found to be very close to the experimental value when a quasi-steady flow field is allowed to develop by continuously venting the tank prior to pressurization as was done during testing. On the other hand, without initial venting, the rate is not constant. Additionally, the development time to reach a quasi-steady flow field differed from the experiment due to the lack of turbulent mixing away from the gas-liquid interface that is expected in reality. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 15% of the experimental value for 49% fill level at 3.5 W/m 2 and is within 18% for 2.0 W/m 2 heat loading. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 18% of the experimental value for a 29% fill level and within 15% for an 83% fill level.

H. Q. Yang↗
NASA NTRSDate not supplied

Computational Fluid Dynamics Simulation of Methane Slosh and Drain Experiments

NASA possesses a wealth of historical cryogenic experiments that provide valuable insights into the design, troubleshooting, and understanding involved in the complex fluid and thermodynamics of managing cryogenic propellants. One approach to leveraging these historical datasets is by simulating these experiments to validate the accuracy of simulation environments and constituent models. This study focused on simulating a selection of the K site test series for both static and sloshing pressurized liquid methane draining experiments conducted at NASA in the 1970’s, utilizing computational fluid dynamics. The simulations were performed in the ANSYS FLUENT environment using the Volume of Fluid (VOF) numerical approach. A k-omega turbulence model was used with interfacial turbulence damping, and accurately predicted the amount of pressurant needed to maintain the required tank pressure throughout the static drain. The static simulation predicted the temperature stratification in the ullage observed at the end of the drain. During the methane expulsion with sloshing test, many features were successfully captured using the k-omega turbulence model with interfacial turbulence damping included. The rate of phase change and liquid temperature was overpredicted compared to the experimental measurements. The overprediction may be attributed to uncertainties in the vessel geometry, methane pressurant temperature and composition, and methodological differences in how the sloshing frequency was adjusted during the expulsion.

Cryogenic Propellants↗
NASA NTRS1986

Computational and experimental study of turbulent flows in 180 degree bends

A theoretical and experimental investigation of turbulent flows inside a 180-deg strongly curved channel is included in the present study. A 10 by 100 centimeter water channel with a 10 centimeter radius turn was employed for the experimental study. Data measured are velocity profiles and distributions of longitudinal turbulence intensity. Numerical investigations are performed by solving the incompressible Navier-Stokes equations and turbulence models using a finite difference flow solver in generalized curvilinear coordinate systems. A widely used two-equation turbulence model and its modified version, which includes the effect of strong curvature on the turbulence length scale through a relation of algebraic Reynolds stress model, are employed to predict the turbulence quantities inside the flow field. Improvements on the predictions are obtained using the present approach.

Chen, Y. S.↗
NASA NTRS2011

Toward Supersonic Retropropulsion CFD Validation

This paper begins the process of verifying and validating computational fluid dynamics (CFD) codes for supersonic retropropulsive flows. Four CFD codes (DPLR, FUN3D, OVERFLOW, and US3D) are used to perform various numerical and physical modeling studies toward the goal of comparing predictions with a wind tunnel experiment specifically designed to support CFD validation. Numerical studies run the gamut in rigor from code-to-code comparisons to observed order-of-accuracy tests. Results indicate that this complex flowfield, involving time-dependent shocks and vortex shedding, design order of accuracy is not clearly evident. Also explored is the extent of physical modeling necessary to predict the salient flowfield features found in high-speed Schlieren images and surface pressure measurements taken during the validation experiment. Physical modeling studies include geometric items such as wind tunnel wall and sting mount interference, as well as turbulence modeling that ranges from a RANS (Reynolds-Averaged Navier-Stokes) 2-equation model to DES (Detached Eddy Simulation) models. These studies indicate that tunnel wall interference is minimal for the cases investigated; model mounting hardware effects are confined to the aft end of the model; and sparse grid resolution and turbulence modeling can damp or entirely dissipate the unsteadiness of this self-excited flow.

Kleb, Bil↗
NASA NTRS1985

Turbulent flow model for vapor collection efficiency of a high-purity silicon reactor

In this study a mathematical model and a computer code based on this model was developed to allow prediction of the product distribution in chemical reactors for converting gaseous silicon compounds to condensed-phase silicon. Specifically, the model formulated describes the silicon vapor separation/collection from the developing turbulent flow stream within reactors of the Westinghouse type. Migration of the silicon vapor to the reactor walls was described by the parametric solutions presented here, in order to reduce the experimentation necessary in the design of such reactors. Calculations relating to the collection efficiencies of such reactors are presented as a function of the reactor throughflow and distance along its length.

Srivastava, R.↗
NASA NTRS1995

Efficiency and Accuracy of Time-Accurate Turbulent Navier-Stokes Computations

The accuracy and efficiency of two types of subiterations in both explicit and implicit Navier-Stokes codes are explored for unsteady laminar circular-cylinder flow and unsteady turbulent flow over an 18-percent-thick circular-arc (biconvex) airfoil. Grid and time-step studies are used to assess the numerical accuracy of the methods. Nonsubiterative time-stepping schemes and schemes with physical time subiterations are subject to time-step limitations in practice that are removed by pseudo time sub-iterations. Computations for the circular-arc airfoil indicate that a one-equation turbulence model predicts the unsteady separated flow better than an algebraic turbulence model; also, the hysteresis with Mach number of the self-excited unsteadiness due to shock and boundary-layer separation is well predicted.

Rumsey, Christopher L.↗
DOE OSTI2021

Assessment of Turbulent Combustion Models for Simulating Prechamber Ignition in a Natural Gas Engine

Lean combustion is a promising strategy to increase thermal efficiency in an internal combustion engine, by exploiting a favorable specific heat ratio of the fresh mixture while simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition and slow flame propagation at fuel-lean conditions lead to large cycle-to-cycle variability and limit the high-efficiency engine operating range. Prechamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling a faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling standpoint to date science base and available simulation tools are inadequate to properly understand and predict the combustion processes in prechamber ignited engines. In this paper, conceptually different Reynolds-averaged Navier–Stokes (RANS) combustion models widely adopted in the engine modeling community are used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a prechamber spark-ignition system. Here, a flamelet-based turbulent combustion model, i.e., G-equation, and a multizone well-stirred reactor model are employed for this modeling study. Simulation results are compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two evaluated approaches.

42 ENGINEERING↗
NASA NTRSDate not supplied

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗
NASA NTRSDate not supplied

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗
NASA NTRSDate not supplied

1st AIAA CFD Transition Modeling and Prediction Workshop: OVERFLOW Results for the SST and Langtry-Menter Models

This paper reports the results of CFD simulations of the test cases for the AIAA 1stCFDTransition Modeling and Prediction Workshop held at the AIAA Scitech Forum on 21-22January, 2021 using the NASA solver OVERFLOW. The Langtry-Menter transition model was used in combination with the Shear Stress Transport (SST) turbulence model. The testcases include the ERCOFTAC zero-pressure-gradient flat plate cases T3A and T3B with by-pass transition, the natural laminar flow NLF(1)-0416 airfoil, the DLR inclined 6:1 prolate spheroid, and the NASA Common Research Model with wing design optimized for natural laminar flow transition (CRM-NLF). OVERFLOW load and transition predictions obtained using best practices are compared with experiments and prior simulations. Additional results are provided for modified turbulence model boundary specifications to improve alignment to reference data and assess model sensitivity.

ARMD↗