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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The Structure and Stability of Laminar Flames

This review paper on the structure and stability of laminar flames considers such phenomena as heterogeneous mixtures, acoustic instabilities, flame balls and related phenomena, radiation effects, the iodate oxidation of arsenous acid and 'liquid flame fronts', approximate kinetic mechanisms and asymptotic approximations, and tribrachial or triple flames. The topics examined here indicate three themes that may play an important role in laminar flame theory in the coming years: microgravity experiments, kinetic modeling, and turbulence modeling. In the discussion of microgravity experiments it is pointed out that access to drop towers, the Space Shuttle and, in due course, the Space Station Freedom will encourage the development of experiments well designed to isolate the fundamental physics of combustion.

Buckmaster, John↗
NASA NTRS1994

Computational/experimental investigation of staged injection into a Mach 2 flow

The stage normal injection of two jets of air into a Mach 2 freestream behind a rearward-facing step has been investigated using laser-induced iodine fluorescence and laser doppler anemometry techniques. A detailed data set has been compiled that includes profiles of pressure, temperature, two components of velocity, and mole fraction of the injectant at 44 locations, plus planar surveys of pressure, temperature, velocity, and mole fraction. A companion numerical study was performed using the SPARK Navier-Stokes code. Good overall agreement was observed for this complex, highly three-dimensional flowfield. Discrepancies were observed in the computed and measured total temperatures, turbulent mixing, and shock strengths. The effect of grid resolution was investigated by calculating solutions on grids of 60,000 points, 200,000 points, and 450,000 points. Differences in the solutions on the two finer grids were small. The effect of turbulence modeling was investigated by calculating solutions with three different algebraic models for the jet turbulence. Overall, the turbulence models were found to have the greatest effect on the numerical solutions, followed by the grid resolution and the injectant Mach number.

Eklund, D. R.↗
NASA NTRSDate not supplied

HIFiRE Flight 2 Unstart Reliability Analysis

The objective of this work was to assess the unstart reliability of the HIFiRE Flight 2 system. To do this, a quantification of margins and uncertainties framework was used for comparing model predictions of the predicted combustion induced shock location to the predicted last stable shock location within the isolator. Uncertainty sources included parametric uncertainty in the flight conditions, the heat release model, and turbulence modeling, as well as model verification errors. Additionally, an estimate of model-form uncertainty was established by comparing the model to measured ground test data. A computationally efficient non-intrusive polynomial chaos approach was used to propagate parametric uncertainty through the computational flu-ids dynamics models of both the ground test configuration and the flight vehicle. Compared to direct-connect ground test data, computational fluid dynamics predictions yielded about two duct heights of model-form uncertainty. This was applied to a prediction of flight vehicle unstart margin at the Mach 6.5 flight condition. Building up all of the computational model uncertainty, including parametric uncertainty, verification errors, and the determined model-form uncertainty, a 95% probability level based confidence ratio, or a ratio of total uncertainty to a statistical margin measure, was found to be 0.31 for the flight system.

Thomas K West↗
NASA NTRS1993

Application of a dynamic subgrid-scale model to turbulent recirculating flows

The dynamic subgrid-scale model of Germano et al. is implemented in a finite volume formulation and applied to the simulation of turbulent flow in a three-dimensional lid-driven cavity at Reynolds number of 7500. The filtering operation is carried out in physical space, and the model coefficient is calculated locally. The computed mean and rms velocities as well as the Reynolds stress are compared with experimental data. It is shown that backscatter from small to large scales is necessary to sustain turbulent fluctuations. The model is being applied to the simulation of turbulent flows in a stratified and rotating environment in complex geometries.

Zang, Y.↗
NASA NTRS1977

A critique of some recent second-order turbulence closure models for compressible boundary layers

Computations based on two recently developed second-order turbulence closure models are compared with a series of boundary-layer experiments and with predictions of these experiments using an algebraic mixing length model. One of the models employs an eddy viscosity, whereas the other evaluates components of the Reynolds stress tensor. For flat plates, the computations are compared with the van Driest skin-friction transformation to assess the handling of compressibility. For boundary layers in pressure gradients, four experiments at Mach 4 and one at Mach 6.7 are used as the bases for comparison. In general, both models represent mean velocities and skin friction reasonably well, but represent the turbulence shear stress less accurately.

Rubesin, M. W.↗
NASA NTRS1973

A stochastic model of turbulent mixing with chemical reaction: Nitric oxide formulation in a plug-flow burner

A stochastic model of turbulent mixing was developed for a reactor in which mixing is represented by n-body fluid particle interactions. The model was used to justify the assumption (made in previous investigations of the role of turbulent mixing on burner generated thermal nitric oxide and carbon monoxide emissions) that for a simple plug flow reactor, composition nonuniformities can be described by a Gaussian distribution function in the local fuel:air equivalence ratio. Recent extensions of this stochastic model to include the combined effects of turbulent mixing and secondary air entrainment on thermal generation of nitric oxide in gas turbine combustors are discussed. Finally, rate limited upper and lower bounds of the nitric oxide produced by thermal fixation of molecular nitrogen and oxidation of organically bound fuel nitrogen are estimated on the basis of the stochastic model for a plug flow burner; these are compared with experimental measurements obtained using a laboratory burner operated over a wide range of test conditions; good agreement is obtained.

Flagan, R. C.↗
NASA NTRS2015

OpenFOAM Simulations of Atmospheric-Entry Capsules in the Subsonic Regime

The open-source Computational Fluid Dynamics software OpenFOAM is gaining wider acceptance in industry and academia for incompressible flow simulations. To date, there has been relatively little utilization of OpenFOAM for compressible external aerodynamic applications. The numerous turbulence models available in OpenFOAM makes it an attractive option for evaluating alternate Reynolds-Averaged Navier-Stokes (RANS) turbulent models to assess separated flow on atmospheric entry vehicles in the subsonic regime, where traditional turbulent models show reduced accuracy. This paper presents simulations of an axisymmetric capsule geometry at subsonic conditions using an OpenFOAM compressible flow solver. These results are compared with results from the NASA CFD code OVERFLOW and experimental data. These OpenFOAM simulations serve as a basis to explore OpenFOAMs extended turbulence models on compressible separated flows such as found on entry capsules.

Nikaido, Ben E.↗
NASA NTRS1991

Turbulence closure modeling near rigid boundaries

The near-wall region plays an essential role in turbulent boundary layers: it is a region of high shear; the peak rate of production and peak intensity of turbulence occurs there; and the peak rate of dissipation occurs right at the wall. Nevertheless, this region has received less attention from modelers than have more nearly homogeneous flows. One reason for this is that when the boundary layer is near equilibrium, experimental data can be used to prescribe the flow in the wall layer. Another reason is that most turbulence models are developed under assumptions of near homogeneity. This is a poor approximation in the wall region. A single-point moment closure model for the strongly non-homogeneous A turbulent flow near a rigid boundary is developed.

Durbin, Paul A.↗
NASA NTRS2018

Numerical Study of Turbulent Separation Bubbles with Varying Pressure Gradient and Reynolds Number

A family of cases each containing a small separation bubble is treated by direct numerical simulation (DNS), varying two parameters: the severity of the pressure gradients, generated by suction and blowing across the opposite boundary, and the Reynolds number. Each flow contains a well-developed entry region with essentially zero pressure gradient, and all are adjusted to have the same value for the momentum thickness, extrapolated from the entry region to the centre of the separation bubble. Combined with fully defined boundary conditions this will make comparisons with other simulations and turbulence models rigorous; we present results for a set of eight Reynolds-averaged Navier–Stokes turbulence models. Even though the largest Reynolds number is approximately 5.5 times higher than in a similar DNS study we presented in 1997, the models have difficulties matching the DNS skin friction very closely even in the zero pressure gradient, which complicates their assessment. In the rest of the domain, the separation location per se is not particularly difficult to predict, and the most definite disagreement between DNS and models is near reattachment. Curiously, the better models tend to cluster together in their predictions of pressure and skin friction even when they deviate from the DNS, although their eddy-viscosity levels are widely different in the outer region near the bubble (or they do not rely on an eddy viscosity). Stratford’s square-root law is satisfied by the velocity profiles, both at separation and reattachment. The Reynolds-number range covers a factor of two, with the Reynolds number based on the extrapolated momentum thickness equal to approximately 1500 and 3000. This allows tentative estimates of the improvements that even higher values will bring to the model comparisons. The solutions are used to assess models through pressure, skin friction and other measures; the flow fields are also used to produce effective eddy-viscosity targets for the models, thus guiding turbulence-modelling work in each region of the flow.

Coleman, G. N.↗
NASA NTRS2017

Grid Convergence for Three Dimensional Benchmark Turbulent Flows

Grid convergence studies are performed to establish reference solutions for benchmark three dimensional turbulent flows in support of the ongoing turbulence model verification and validation e ort at the Turbulence Modeling Resource website curated by NASA. The bench- mark cases are a subsonic flow around a hemisphere cylinder and a transonic flow around the ONERA M6 wing with a sharp trailing edge. The study applies widely-used computational fluid dynamics codes developed and supported at the NASA Langley Research Center: FUN3D, USM3D, and CFL3D. Reference steady-state solutions are computed for the Reynolds-Averaged Navier-Stokes equations with the Spalart-Allmaras turbulence model on families of consistently-refined grids composed of different types of cells. Coarse-to- ne and code-to-code solution variation is described in detail.

Diskin, Boris↗
DOE OSTI2025

Adaptive Grid Redistribution for a 1D Model of Turbulence and Clouds

In global atmospheric models, resolving stratocumulus (Sc) in the vertical is computationally expensive. However, Sc appear only under special meteorological conditions. Therefore, there is motivation to refine the vertical grid levels adaptively. In order to facilitate the possibility of parallelization on graphical processing units, our grid adaptation method prescribes the number of vertical levels a priori. Then grid levels are relocated toward altitude ranges in need of refinement. Because the method relocates existing grid levels, rather than adding extra levels, there is a risk of creating regions with overly coarse grid spacing, that is, voids in the grid mesh. To prevent such voids from forming, a simple method is developed to impose a maximum grid spacing. To decide where to place enhanced resolution, the authors develop an empirical mesh refinement criterion. It refines grid spacing near the ground, near strong temperature gradients, and within clouds. Our grid adaptation method is implemented in a single-column model and evaluated on four test cases: decaying stratocumulus, developing shallow cumulus, a quasi-stationary stratocumulus deck, and the diurnal cycle of a dry boundary layer. In the stratocumulus cases, mesh refinement leads to improvements in both the time evolution of fields and their time averages. The other two cases show smaller differences.

Carstensen, Steffen [Univ. of Wisconsin, Milwaukee↗
NASA NTRS1991

Low Reynolds number two-equation modeling of turbulent flows

A k-epsilon model that accounts for viscous and wall effects is presented. The proposed formulation does not contain the local wall distance thereby making very simple the application to complex geometries. The formulation is based on an existing k-epsilon model that proved to fit very well with the results of direct numerical simulation. The new form is compared with nine different two-equation models and with direct numerical simulation for a fully developed channel flow at Re = 3300. The simple flow configuration allows a comparison free from numerical inaccuracies. The computed results prove that few of the considered forms exhibit a satisfactory agreement with the channel flow data. The model shows an improvement with respect to the existing formulations.

Michelassi, V.↗
NASA NTRS1992

Computation of Turbulent Recirculating Flow in Channels, and for Impingement Cooling

Fully elliptic forms of the transport equations have been solved numerically for two flow configurations. The first is turbulent flow in a channel with transverse rectangular ribs, and the second is impingement cooling of a plane surface. Both flows are relevant to proposed designs for active cooling of hypersonic vehicles using supercritical hydrogen as the coolant. Flow downstream of an abrupt pipe expansion and of a backward-facing step were also solved with various near-wall turbulence models as benchmark problems. A simple form of periodicity boundary condition was used for the channel flow with transverse rectangular ribs. The effects of various parameters on heat transfer in channel flow with transverse ribs and in impingement cooling were investigated using the Yap modified Jones and Launder low Reynolds number k-epsilon turbulence model. For the channel flow, predictions were in adequate agreement with experiment for constant property flow, with the results for friction superior to those for heat transfer. For impingement cooling, the agreement with experiment was generally good, but the results suggest that improved modelling of the dissipation rate of turbulence kinetic energy is required in order to obtain improved heat transfer prediction, especially near the stagnation point. The k-epsilon turbulence model was used to predict the mean flow and heat transfer for constant and variable property flows. The effect of variable properties for channel flow was investigated using the same turbulence model, but comparison with experiment yielded no clear conclusions. Also, the wall function method was modified for use in the variable properties flow with a non-adiabatic surface, and an empirical model is suggested to correctly account for the behavior of the viscous sublayer with heating.

Chang, Byong Hoon↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Full-wave model for the lower hybrid wave electric field vector with synthetic turbulence on Alcator C-Mod

The effects of synthetic turbulence on lower hybrid (LH) wave field magnitude and polarization are studied on Alcator C-Mod. Three synthetic turbulence models are evaluated to assess the impact of filaments, holes, and periodic fluctuations on the LH wave electric fields. For all three cases, the synthetic turbulence can greatly impact the LH wave magnitude and polarization. Back scattering, interference and edge/SOL absorption of the LH wave is observed. This impact is shown to depend on the wavelength and amplitude of the fluctuations. The strongest effect appears for fluctuations at high amplitudes and wavelengths comparable to the LH wavelength. This synthetic turbulence model can be used as inputs into a synthetic diagnostic to calculate the lower hybrid wave field magnitude and direction measured by dynamic Stark effect spectroscopy. It will be shown that fluctuations at these high amplitudes and wavelengths may explain the experimentally measured results. Here, the effect of fluctuations on modifying LH polarization are also shown to be anti-correlated with SOL power losses, similar to recently observed experimental trends on Alcator C-Mod. These SOL losses may have detrimental impact on LHCD efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Simulation as an aid to phenomenological modeling

Certain difficulties arise in connection with the development of models for the prediction of turbulent flow characteristics. The present paper is mainly concerned with an approach which provides an aid in such development work by making it possible to test Reynolds average turbulence models. This approach is based on full and large eddy simulations of turbulent flows. Attention is given to details regarding Reynolds average turbulence models, mixing length models, two equation (k-epsilon) models, Reynolds stress models, and passive scalars. It is shown that a great deal about modeling can be learned from simulations of turbulent flows. There is reason to believe that these results will lead to improvements in the quality of predictions of turbulent flows.

Ferziger, J. H.↗
NASA NTRS2007

Flow Control Predictions using URANS Modeling: A Parametric Study

A computational study was performed for steady and oscillatory flow control over a hump model with flow separation to assess how well the steady and unsteady Reynolds-averaged Navier-Stokes equations predict trends due to Reynolds number, control magnitude, and control frequency. As demonstrated previously, the hump model case is useful because it clearly demonstrates a failing in all known turbulence models: they under-predict the turbulent shear stress in the separated region and consequently reattachment occurs too far downstream. In spite of this known failing, three different turbulence models were employed to determine if trends can be captured even though absolute levels are not. The three turbulence models behaved similarly. Overall they showed very similar trends as experiment for steady suction, but only agreed qualitatively with some of the trends for oscillatory control.

Rumsey, Christopher L.↗
DOE OSTI2026

Simulation of Channel Flow with Square Ribs for Blanket First-Wall Cooling: Geometry-Specific Tuning of k-ω Model Using Adjoint Method

Cooling of the plasma-facing first wall is challenging in the design of blanket components because of the high heat flux (on the order of 𝑀𝑊/𝑚2) from the plasma, especially when a low thermal mass medium like helium is chosen as the coolant. Therefore, heat transfer enhancement in which the convective heat transfer rate is augmented by the addition of turbulence-promoting structures becomes a key initiative for providing sufficient cooling capability with helium. Previously, computational fluid dynamics simulations had been performed on pipe flows with different transverse and longitudinal ribbed geometries at Oak Ridge National Laboratory to compare the enhancement performance among different ribbed geometries. Rib shape morphing had been conducted to obtain an optimized rib profile. In the work presented here, the adjoint method is adopted in the ANSYS Fluent solver for turbulence model augmentation, and the Generalized k-ω (GEKO) turbulence model is employed because of its ability of tuning the turbulence model. The Nusselt number and pressure drop obtained from the channel flow with bottom ribbed wall experiments are used as the targets. Sensitivity analysis provides information as guidance to improve the turbulence model accuracy. The augmented GEKO model is tuned for the studied ribbed channel geometry and flow conditions, providing improved predictive accuracy within this context. Extension to other configurations offers potential but may require additional tuning and validation.

Xu, Tracy [ORNL] (ORCID:0009000193700887)↗