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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A Numerical Investigation of Turbulent Flow in Noncircular Ducts

Since turbulent duct flows play an important role in engineering, continuous efforts to investigate this problem have been made. Until recently, these investigations were limited mostly to experiments and various semi-analytic methods. Owing to the improvement of turbulence modeling, the prevailing method of prediction is now mainly numerical. The majority of these studies deals, however, with turbulent flows in circular or square ducts and only limited information is available for straight noncircular ducts. In view of this situation, we propose to conduct a numerical investigation of turbulent flow in a class of ducts, whose cross sections vary from a circle to a near square. Turbulent flow in a noncircular duct is characterized by the presence of secondary flow for which a more refined turbulence model than the k-epsilon equations is required. In order to show that the calculated results are credible, various modes of verification were used to examine the results for a selected configuration, including an accuracy check by a scaling law and observing the decay of secondary flow as the cross section changes from a noncircular shape to a circle. After this was done, computations were performed for other configurations and with different Reynolds numbers from which wall shear stresses and friction factors are plotted.

Kao, Hsiao C.↗
NASA NTRS2016

Hybrid Reynolds-Averaged/Large Eddy Simulation of the Flow in a Model SCRamjet Cavity Flameholder

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. Experimental data available for this configuration include velocity statistics obtained from particle image velocimetry. Several turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged/large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This e ort was undertaken to not only assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community, but to also begin to understand how this capability can best be used to augment standard Reynolds-averaged simulations. The numerical errors were quantified for the steady-state simulations, and at least qualitatively assessed for the scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results displayed a high degree of variability when comparing the flameholder fuel distributions obtained from each turbulence model. This prompted the consideration of applying the higher-fidelity scale-resolving simulations as a surrogate "truth" model to calibrate the Reynolds-averaged closures in a non-reacting setting prior to their use for the combusting simulations. In general, the Reynolds-averaged velocity profile predictions at the lowest fueling level matched the particle imaging measurements almost as well as was observed for the non-reacting condition. However, the velocity field predictions proved to be more sensitive to the flameholder fueling rate than was indicated in the measurements.

Baurle, R. A.↗
NASA NTRS1993

Calculations of turbulent separated flows

A numerical study of incompressible turbulent separated flows is carried out by using two-equation turbulence models of the K-epsilon type. On the basis of realizability analysis, a new formulation of the eddy-viscosity is proposed which ensures the positiveness of turbulent normal stresses - a realizability condition that most existing two-equation turbulence models are unable to satisfy. The present model is applied to calculate two backward-facing step flows. Calculations with the standard K-epsilon model and a recently developed RNG-based K-epsilon model are also made for comparison. The calculations are performed with a finite-volume method. A second-order accurate differencing scheme and sufficiently fine grids are used to ensure the numerical accuracy of solutions. The calculated results are compared with the experimental data for both mean and turbulent quantities. The comparison shows that the present model performs quite well for separated flows.

Zhu, J.↗
NASA NTRS1992

CFD modeling of turbulent flows around the SSME main injector assembly using porosity formulation

Hot gas turbulent flow distribution around the main injector assembly of the Space Shuttle Main Engine (SSME) and LOX flow distribution through the LOX posts have a great effect on the combustion phenomena inside the main combustion chamber. In order to design a CFD model to be an effective engineering analysis tool with good computational turn-around time (especially for 3-D flow problems) and still maintain good accuracy in describing the flow features, the concept of porosity was employed to describe the effects of blockage and drag force due to the presence of the LOX posts in the turbulent flow field around the main injector assembly of the SSME. Two-dimensional numerical studies were conducted to identify the drag coefficients of the flows, both through tube banks and round the shielded posts, over a wide range of Reynolds numbers. Empirical, analytical expressions of the drag coefficients as a function of local flow Reynolds number were then deduced. The porosity model was applied to the turbulent flow around the main injector assembly of the SSME, and analyses were performed. The 3-D CFD analysis was divided into three parts: LOX dome, hot gas injector assembly, and hydrogen cavity. The numerical results indicate that the mixture ratio at the downstream of injector face was close to stoichiometric around baffle elements.

Cheng, Gary C.↗
NASA NTRS2001

Isolating Curvature Effects in Computing Wall-Bounded Turbulent Flows

The flow over the zero-pressure-gradient So-Mellor convex curved wall is simulated using the Navier-Stokes equations. An inviscid effective outer wall shape, undocumented in the experiment, is obtained by using an adjoint optimization method with the desired pressure distribution on the inner wall as the cost function. Using this wall shape with a Navier-Stokes method, the abilities of various turbulence models to simulate the effects of curvature without the complicating factor of streamwise pressure gradient can be evaluated. The one-equation Spalart-Allmaras turbulence model overpredicts eddy viscosity, and its boundary layer profiles are too full. A curvature-corrected version of this model improves results, which are sensitive to the choice of a particular constant. An explicit algebraic stress model does a reasonable job predicting this flow field. However, results can be slightly improved by modifying the assumption on anisotropy equilibrium in the model's derivation. The resulting curvature-corrected explicit algebraic stress model possesses no heuristic functions or additional constants. It lowers slightly the computed skin friction coefficient and the turbulent stress levels for this case (in better agreement with experiment), but the effect on computed velocity profiles is very small.

Rumsey, Christopher L.↗
NASA NTRS1980

Conceptual model of turbulent flameholding for scramjet combustors

New concepts and approaches to scramjet combustor design are presented. Blowoff was from failure of the recirculation-zone (RZ) flame to reach the dividing streamline (DS) at the rear stagnation zone. Increased turbulent exchange across the DS helped flameholding due to forward movement of the flame anchor point inside the RZ. Modeling of the blowoff phenomenon was based on a mass conservation concept involving the traverse of a flame element across the RZ and a flow element along the DS. The scale required to achieve flameholding, predicted by the model, showed a strong adverse effect of low pressure and low fuel equivalence ratio, moderate effect of flight Mach number, and little effect of temperature recovery factor. Possible effects of finite rate chemistry on flameholding and flamespreading in scramjets are discussed and recommendations for approaches to engine combustor design as well as for needed research to reduce uncertainties in the concepts are made.

Huber, P. W.↗
NASA NTRS2016

Computational and Experimental Study of Supersonic Nozzle Flow and Aft-Deck Interactions

NASA has been conducting research into reducing sonic boom and changing FAA regulations to allow for supersonic commercial transport over land in the United States. This particular study looks at a plume passing through a shock generated from an aft deck on a nacelle; the aft deck is meant to represent the trailing edge of a wing. NASA Langley Research Center USM3D CFD code results are compared to the experimental data taken at the NASA Glenn Research Center 1-foot by 1-foot Supersonic Wind Tunnel. This study included examining two turbulence models along with different volume sourcing methods for grid generation. The results show that using the k-epsilon turbulence model within USM3D produced shock signatures that closely follow the experimental data at a variety of nozzle pressure ratio settings.

Bruce, Walter E., IV↗
NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the Volume of Fluid (VOF) and Sharp Interface multiphase 2D axisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. Turbulence modeling; turbulence damping at the interface; and using constant vs. temperature dependent fluid properties were shown to have the most profound influence on predicted tank pressures and temperatures. Current study indicates that including phase change at the interface into the computational model is crucial for accurate prediction of the tank self-pressurization process. The effect of accommodation coefficient was also studied. Tank pressure values predicted by the VOF, and Sharp Interface models are within 4% of the experimental ones.

Computational Fluid Dynamics↗
DOE OSTI2026

High-Fidelity Simulation Aerodynamics Dataset of NACA 0012 and 0021 Airfoils

The repository contains time series data of pressure, viscous and moment forces for the NACA 0012 and 0021 airfoils at angles of attack (AOA), 5, 17, 30, 45, 60 and 90 degrees. The data was generated during the high-fidelity simulations of the airfoils using the open-source CFD code, Nalu-Wind (https://github.com/Exawind/nalu-wind). The research objective was to study the effects of mesh resolution and turbulence model on the three-dimensional deep stall aerodynamics of the airfoils, which has been published in the Journal of Turbulence (https://doi.org/10.1080/14685248.2023.2225141). Several meshes of the O-grid type were generated using the modified Pointwise Glyph script of Carrigan (https://github.com/pointwise/AirfoilMesh) by varying the wall normal and spanwise resolutions. Note that, the geometry was extruded by four chord lengths in the spanwise direction. In the paper, we showed that Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS/LES turbulence model with a minimum of 24 cells per chord length in the spanwise direction is necessary to correctly predict the loads in the deep stall regime. Under the NACA 0012 airfoil directory, the data for several combinations of wall-normal and spanwise resolutions are presented in each of the AOA subdirectories. The corresponding experimental data for the lift and drag polars can be found in the "exp_data" subdirectory. The directory for the NACA 0021 airfoil consists of CFD force data and airfoil polars for Reynolds numbers 2.7e5 and 2.0e6.

17 WIND ENERGY↗
NASA NTRS1973

Some problems of the calculation of three-dimensional boundary layer flows on general configurations

An accurate solution of the three-dimensional boundary layer equations over general configurations such as those encountered in aircraft and space shuttle design requires a very efficient, fast, and accurate numerical method with suitable turbulence models for the Reynolds stresses. The efficiency, speed, and accuracy of a three-dimensional numerical method together with the turbulence models for the Reynolds stresses are examined. The numerical method is the implicit two-point finite difference approach (Box Method) developed by Keller and applied to the boundary layer equations by Keller and Cebeci. In addition, a study of some of the problems that may arise in the solution of these equations for three-dimensional boundary layer flows over general configurations.

Cebeci, T.↗
NASA NTRS1990

Numerical study of twin-jet impingement upwash flow

Two horizontally spaced jets impinging normally on a flat surface create a fountain upwash flow due to the collision of the radially flowing wall jets. This fountain flow is of importance to the dynamics and propulsion of STOVL aircraft. The fountain flow influences the lift forces on the aircraft and the ingestion of hot gases and debris by the engine inlet. In this paper, a multigrid based finite-difference numerical procedure has been applied to solve the equations governing this three-dimensional flow. The standard k-epsilon turbulence model has been used. Comparisons with experimental data reveal that while the mean velocities are predicted with reasonable accuracy, the turbulent kinetic energies are seriously in error. The reasons for this discrepancy could be the intense unsteadiness and large-scale structures of the flow in the near-wall region, which cannot be captured well by any Reynolds-averaged turbulence model.

Pegues, W. J.↗
NASA NTRS1990

A Detailed Investigation of Staged Normal Injection into a Mach 2 Flow

A study of the staged injection of two jets of air behind a rearward facing step into a Mach 2 flow was performed using the SPARK 3-D Navier-Stokes code. Calculated mole fraction distributions were compared with an extensive set of planar mole fraction measurements made with a laser induced iodine fluorescence technique. A statistical measure, the standard deviation, was used to help assess agreement between calculation and experiment. Overall, good agreement was found between calculated and measured values. Generally, agreement was better in the far field of the injectors. The effect of grid resolution was investigated by calculating solutions on grids of 60,000, 200,000, and 450,000 points. Differences in the solutions on the two finer grids were small. However, the mole fraction distributions were distinguishable. The effect of turbulence modeling was investigated by employing three different algebraic models for the jet turbulence: the Baldwin-Lomax model, the Prandtl mixing length model, and the Eggers mixing length model. Overall, the Eggers mixing length model was found to be superior for this case. Finally, the effect of the jet exit conditions was examined. A recently proposed Mach number distribution at the jet exit was found to slightly improve agreement between measurement and calculation.

Eklund, Dean R.↗
NASA NTRS1996

Numerical simulation of jet aerodynamics using the three-dimensional Navier-Stokes code PAB3D

This report presents a unified method for subsonic and supersonic jet analysis using the three-dimensional Navier-Stokes code PAB3D. The Navier-Stokes code was used to obtain solutions for axisymmetric jets with on-design operating conditions at Mach numbers ranging from 0.6 to 3.0, supersonic jets containing weak shocks and Mach disks, and supersonic jets with nonaxisymmetric nozzle exit geometries. This report discusses computational methods, code implementation, computed results, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions. The Navier-Stokes method using the standard Jones-Launder two-equation kappa-epsilon turbulence model can accurately predict jet flow, and such predictions are made without any modification to the published constants for the turbulence model.

Pao, S. Paul↗
NASA NTRS2009

Enhancement of USM3D Unstructured Flow Solver for High-Speed High-Temperature Shear Flows

Large temperature and pressure fluctuations have a profound effect on turbulence development in transonic and supersonic jets. For high-speed, high-temperature jet flows, standard turbulence models lack the ability to predict the observed mixing rate of a shear layer. Several proposals to address this deficiency have been advanced in the literature to modify the turbulence transport equations in a variety of ways. In the present study, some of the most proven and simple modifications to two-equation turbulence models have been selected and implemented in NASA's USM3D tetrahedral Navier-Stokes flow solver. The modifications include the addition of compressibility correction and pressure dilatation terms in the turbulence transport equations for high-speed flows, and the addition of a simple modification to the Boussinesq's closure model coefficient for high-temperature jets. The efficacy of the extended models is demonstrated by comparison with experimental data for two supersonic axisymmetric jet test cases at design pressure ratio.

Pandya, Mohagna J.↗
NASA NTRS2001

Unsteady Analysis of Separated Aerodynamic Flows Using an Unstructured Multigrid Algorithm

An implicit method for the computation of unsteady flows on unstructured grids is presented. The resulting nonlinear system of equations is solved at each time step using an agglomeration multigrid procedure. The method allows for arbitrarily large time steps and is efficient in terms of computational effort and storage. Validation of the code using a one-equation turbulence model is performed for the well-known case of flow over a cylinder. A Detached Eddy Simulation model is also implemented and its performance compared to the one equation Spalart-Allmaras Reynolds Averaged Navier-Stokes (RANS) turbulence model. Validation cases using DES and RANS include flow over a sphere and flow over a NACA 0012 wing including massive stall regimes. The project was driven by the ultimate goal of computing separated flows of aerodynamic interest, such as massive stall or flows over complex non-streamlined geometries.

Pelaez, Juan↗
NASA NTRS2005

Computational Study of Separating Flow in a Planar Subsonic Diffuser

A computational study of the separated flow through a 2-D asymmetric subsonic diffuser has been performed. The Wind Computational Fluid Dynamics code is used to predict the separation and reattachment behavior for an incompressible diffuser flow. The diffuser inlet flow is a two-dimensional, turbulent, and fully-developed channel flow with a Reynolds number of 20,000 based on the centerline velocity and the channel height. Wind solutions computed with the Menter SST, Chien k-epsilon, Spalart-Allmaras and Explicit Algebraic Reynolds Stress turbulence models are compared with experimentally measured velocity profiles and skin friction along the upper and lower walls. In addition to the turbulence model study, the effects of grid resolution and use of wall functions were investigated. The grid studies varied the number of grid points across the diffuser and varied the initial wall spacing from y(sup +) = 0.2 to 60. The wall function study assessed the applicability of wall functions for analysis of separated flow. The SST and Explicit Algebraic Stress models provide the best agreement with experimental data, and it is recommended wall functions should only be used with a high level of caution.

DalBello, Teryn↗
NASA NTRS2015

Hybrid Reynolds-Averaged/Large Eddy Simulation of a Cavity Flameholder; Assessment of Modeling Sensitivities

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. The cases simulated corresponded to those used to examine this flowfield experimentally using particle image velocimetry. A variety of turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged / large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This effort was undertaken to formally assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community. The numerical errors were quantified for both the steady-state and scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results showed a high degree of variability when comparing the predictions obtained from each turbulence model, with the non-linear eddy viscosity model (an explicit algebraic stress model) providing the most accurate prediction of the measured values. The hybrid Reynolds-averaged/large eddy simulation results were carefully scrutinized to ensure that even the coarsest grid had an acceptable level of resolution for large eddy simulation, and that the time-averaged statistics were acceptably accurate. The autocorrelation and its Fourier transform were the primary tools used for this assessment. The statistics extracted from the hybrid simulation strategy proved to be more accurate than the Reynolds-averaged results obtained using the linear eddy viscosity models. However, there was no predictive improvement noted over the results obtained from the explicit Reynolds stress model. Fortunately, the numerical error assessment at most of the axial stations used to compare with measurements clearly indicated that the scale-resolving simulations were improving (i.e. approaching the measured values) as the grid was refined. Hence, unlike a Reynolds-averaged simulation, the hybrid approach provides a mechanism to the end-user for reducing model-form errors.

Baurle, R. A.↗
NASA NTRS1993

Numerical analysis of base flowfield at high altitude for a four-engine clustered nozzle configuration

The objective of this study is to benchmark a four-engine clustered nozzle base flowfield with a computational fluid dynamics (CFD) model. The CFD model is a pressure based, viscous flow formulation. An adaptive upwind scheme is employed for the spatial discretization. The upwind scheme is based on second and fourth order central differencing with adaptive artificial dissipation. Qualitative base flow features such as the reverse jet, wall jet, recompression shock, and plume-plume impingement have been captured. The computed quantitative flow properties such as the radial base pressure distribution, model centerline Mach number and static pressure variation, and base pressure characteristic curve agreed reasonably well with those of the measurement. Parametric study on the effect of grid resolution, turbulence model, inlet boundary condition and difference scheme on convective terms has been performed. The results showed that grid resolution and turbulence model are two primary factors that influence the accuracy of the base flowfield prediction.

Wang, Ten-See↗