Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence simulation”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Verifying Implementation of the Dryden Turbulence Model and MIL-F-8785 Gust Gradient

Turbulence modeling in human-in-the-loop simulation is important to assessing aircraft handling qualities and pilot performance and to provide additional realism for pilot training. In the simulation community, the Dryden turbulence spectra is a popular choice for modeling the linear turbulent gusts because its rational form is efficiently reproduced by passing white noise through linear filters. The MIL-F-8785 gust gradients similarly use additional linear filters to model the gradient of the turbulent gust over the wing, and it represents the gust gradients as perturbations to the air-relative rotational rates. The Cockpit Motion Facility at NASA Langley Research Center (LaRC) models continuous random turbulence using the Dryden one-dimensional spectra and MIL-F-8785 gust gradient. The facility recently reviewed and updated its verification of these models as part of an initiative to improve motion cueing under turbulence. This exercise introduced improved methods for verifying the turbulence models and led to rediscovery of model assumptions that informed improvements to implementation.

Madden, Michael M.↗
NASA NTRS1993

High Reynolds number calculations using the dynamic subgrid-scale stress model

The dynamic subgrid-scale eddy viscosity model was used in the large-eddy simulation of the turbulent flow in a plane channel for Reynolds numbers based on friction velocity and channel half-width ranging between 200 and 2000, a range including values significantly higher than in previous simulations. The computed wall stress, mean velocity, and Reynolds stress profiles compare very well with experimental and direct simulation data. Comparison of higher moments is also satisfactory. Although the grid in the near-wall region is fairly coarse, the results are quite accurate: the turbulent kinetic energy peaks at y(+) of about 12, and the near-wall behavior of the resolved stresses is captured accurately. The model coefficient is o(0.001) in the buffer layer and beyond, where the cutoff wave numbers are in the decaying region of the spectra; in the near-wall region, the cutoff wave numbers are nearer the energy-containing range, and the resolved turbulent stresses become a constant fraction of the resolved stresses. This feature is responsible for the correct near-wall behavior of the model coefficient.

Piomelli, Ugo↗
NASA NTRS2005

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Wang, Chi R.↗
NASA NTRS1984

Interactive computer modeling of combustion chemistry and coalescence-dispersion modeling of turbulent combustion

An interactive computer code for simulation of a high-intensity turbulent combustor as a single point inhomogeneous stirred reactor was developed from an existing batch processing computer code CDPSR. The interactive CDPSR code was used as a guide for interpretation and direction of DOE-sponsored companion experiments utilizing Xenon tracer with optical laser diagnostic techniques to experimentally determine the appropriate mixing frequency, and for validation of CDPSR as a mixing-chemistry model for a laboratory jet-stirred reactor. The coalescence-dispersion model for finite rate mixing was incorporated into an existing interactive code AVCO-MARK I, to enable simulation of a combustor as a modular array of stirred flow and plug flow elements, each having a prescribed finite mixing frequency, or axial distribution of mixing frequency, as appropriate. Further increase the speed and reliability of the batch kinetics integrator code CREKID was increased by rewriting in vectorized form for execution on a vector or parallel processor, and by incorporating numerical techniques which enhance execution speed by permitting specification of a very low accuracy tolerance.

Pratt, D. T.↗
NASA NTRS1993

Summary of the GASP code application and evaluation effort for scramjet combustor flowfields

Numerical simulations of 3D turbulent mixing and reacting flows have been systematically evaluated by comparison with experimental data and other numerical solutions to provide confidence in the General Aerodynamic Simulation Program (GASP) code version 1.3. Specifically, the GASP flow solver has been used to model high speed flow through scram jet combustors. Unit injector problems evaluated include (1) UVA staged sonic normal injection of air into Mach 2 air stream, (2) the VPI Mach 1.7 low angled (15 deg) flush wall helium injection into an unconfined Mach 6 'cold air' stream, and (3) the HYPULSE angled (30 deg) flush wall hydrogen injection into a Mach 6, high enthalpy (flight Mach 17 simulation) nitrogen and air test gases. Results of the numerical simulation for the first two test cases have been compared with detailed in-stream measurements and with SPARK CFD solutions. The results for the last two cases have been compared with measured wall pressure and heat flux data and with SPARK CFD solutions. The present GASP solutions compare favorably with both the experimental data and the SPARK solutions.

Srinivasan, Shivakumar↗
NASA NTRS1994

Computation and modeling of aero-thermal fields in turbine cascades and strongly curved ducts

Advanced turbulence models are crucial for accurate prediction of rocket engine flows, due to existence of very large extra strain rates, such as strong streamline curvature. Numerical simulation of the turbulent flow in a strongly curved turn-around duct (TAD) has been carried out with a Reynolds stress model (RSM), an algebraic Reynolds stress model (ARSM) and a kappa-epsilon model. the RSM model and the ARSM model are found to capture the turbulence damping due to the convex curvature, but underpredict the turbulence enhancement caused by the concave curvature. To capture the concave curvature effects, it is necessary to modify the epsilon-equation. The modification of episilon-equation suggested by Launder, et.al, provides the correct trend, but over-corrects the curvature effects. A comparative study of two modes of transition in gas turbine, the by-pass transition and the separation-induced transition, has been carried out with several low-Reynolds-number (LRN) kappa-epsilon models. Effects of blade surface pressure gradient, freestream turbulence and Reynolds number on the blade boundary layer development, and particularly the inception of transition are examined in detail. The present study indicates that the turbine blade transition, in the presence of high freestream turbulence, is predicted well with LRN kappa-epsilon models employed.

Luo, J.↗
NASA NTRS1996

Turbulence Modeling and Computation of Turbine Aerodynamics and Heat Transfer

The objective of the present research is to develop improved turbulence models for the computation of complex flows through turbomachinery passages, including the effects of streamline curvature, heat transfer and secondary flows. Advanced turbulence models are crucial for accurate prediction of rocket engine flows, due to existance of very large extra strain rates, such as strong streamline curvature. Numerical simulation of the turbulent flows in strongly curved ducts, including two 180-deg ducts, one 90-deg duct and a strongly concave curved turbulent boundary layer have been carried out with Reynolds stress models (RSM) and algebraic Reynolds stress models (ARSM). An improved near-wall pressure-strain correlation has been developed for capturing the anisotropy of turbulence in the concave region. A comparative study of two modes of transition in gas turbine, the by-pass transition and the separation-induced transition, has been carried out with several representative low-Reynolds number (LRN) k-epsilon models. Effects of blade surface pressure gradient, freestream turbulence and Reynolds number on the blade boundary layer development, and particularly the inception of transition are examined in detail. The present study indicates that the turbine blade transition, in the presence of high freestream turbulence, is predicted well with LRN k-epsilon models employed. The three-dimensional Navier-Stokes procedure developed by the present authors has been used to compute the three-dimensional viscous flow through the turbine nozzle passage of a single stage turbine. A low Reynolds number k-epsilon model and a zonal k-epsilon/ARSM (algebraic Reynolds stress model) are utilized for turbulence closure. An assessment of the performance of the turbulence models has been carried out. The two models are found to provide similar predictions for the mean flow parameters, although slight improvement in the prediction of some secondary flow quantities has been obtained by the ARSM model. It's found that the wake profiles inside the endwall boundary layers are predicted better than those near the mid-span.

Lakshminarayana, B.↗
DOE OSTI2020

The Role of Outflows, Radiation Pressure, and Magnetic Fields in Massive Star Formation

Stellar feedback in the form of radiation pressure and magnetically driven collimated outflows may limit the maximum mass that a star can achieve and affect the star formation efficiency of massive prestellar cores. Here we present a series of 3D adaptive mesh refinement radiation–magnetohydrodynamic simulations of the collapse of initially turbulent, massive prestellar cores. Our simulations include radiative feedback from both the direct stellar and dust-reprocessed radiation fields, and collimated outflow feedback from the accreting stars. We find that protostellar outflows punch holes in the dusty circumstellar gas along the star’s polar directions, thereby increasing the size of optically thin regions through which radiation can escape. Precession of the outflows as the star’s spin axis changes due to the turbulent accretion flow further broadens the outflow, and causes more material to be entrained. Additionally, the presence of magnetic fields in the entrained material leads to broader entrained outflows that escape the core. We compare the injected and entrained outflow properties and find that the entrained outflow mass is a factor of ∼3 larger than the injected mass and the momentum and energy contained in the entrained material are ∼25% and ∼5% of the injected momentum and energy, respectively. As a result, we find that, when one includes both outflows and radiation pressure, the former are a much more effective and important feedback mechanism, even for massive stars with significant radiative outputs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1996

Evaluation of Subgrid-Scale Models for Large Eddy Simulation of Compressible Flows

The objective of this project was to evaluate and develop subgrid-scale (SGS) turbulence models for large eddy simulations (LES) of compressible flows. During the first phase of the project results from LES using the dynamic SGS model were compared to those of direct numerical simulations (DNS) of compressible homogeneous turbulence. The second phase of the project involved implementing the dynamic SGS model in a NASA code for simulating supersonic flow over a flat-plate. The model has been successfully coded and a series of simulations has been completed. One of the major findings of the work is that numerical errors associated with the finite differencing scheme used in the code can overwhelm the SGS model and adversely affect the LES results. Attached to this overview are three submitted papers: 'Evaluation of the Dynamic Model for Simulations of Compressible Decaying Isotropic Turbulence'; 'The effect of the formulation of nonlinear terms on aliasing errors in spectral methods'; and 'Large-Eddy Simulation of a Spatially Evolving Compressible Boundary Layer Flow'.

Blaisdell, Gregory A.↗
NASA NTRS2003

Large-Eddy Simulation of Conductive Flows at Low Magnetic Reynolds Number

In this paper we study the LES method with dynamic procedure in the context of conductive flows subject to an applied external magnetic field at low magnetic Reynolds number R(sub m). These kind of flows are encountered in many industrial applications. For example, in the steel industry, applied magnetic fields can be used to damp turbulence in the casting process. In nuclear fusion devices (Tokamaks), liquid-lithium flows are used as coolant blankets and interact with the surrounding magnetic field that drives and confines the fusion plasma. Also, in experimental facilities investigating the dynamo effect, the flow consists of liquid-sodium for which the Prandtl number and, as a consequence, the magnetic Reynolds number is low. Our attention is focused here on the case of homogeneous (initially isotropic) decaying turbulence. The numerical simulations performed mimic the thought experiment described in Moffatt in which an initially homogeneous isotropic conductive flow is suddenly subjected to an applied magnetic field and freely decays without any forcing. Note that this flow was first studied numerically by Schumann. It is well known that in that case, extra damping of turbulence occurs due to the Joule effect and that the flow tends to become progressively independent of the coordinate along the direction of the magnetic field. Our comparison of filtered direct numerical simulation (DNS) predictions and LES predictions show that the dynamic Smagorinsky model enables one to capture successfully the flow with LES, and that it automatically incorporates the effect of the magnetic field on the turbulence. Our paper is organized as follows. In the next section we summarize the LES approach in the case of MHD turbulence at low R(sub m) and recall the definition of the dynamic Smagorinsky model. In Sec. 3 we describe the parameters of the numerical experiments performed and the code used. Section 4 is devoted to the comparison of filtered DNS results and LES results. Conclusions are presented in Sec. 5.

Knaepen, B.↗
NASA NTRS2001

The Numerical Analysis of a Turbulent Compressible Jet

A numerical method to simulate high Reynolds number jet flows was formulated and applied to gain a better understanding of the flow physics. Large-eddy simulation was chosen as the most promising approach to model the turbulent structures due to its compromise between accuracy and computational expense. The filtered Navier-Stokes equations were developed including a total energy form of the energy equation. Subgrid scale models for the momentum and energy equations were adapted from compressible forms of Smagorinsky's original model. The effect of using disparate temporal and spatial accuracy in a numerical scheme was discovered through one-dimensional model problems and a new uniformly fourth-order accurate numerical method was developed. Results from two- and three-dimensional validation exercises show that the code accurately reproduces both viscous and inviscid flows. Numerous axisymmetric jet simulations were performed to investigate the effect of grid resolution, numerical scheme, exit boundary conditions and subgrid scale modeling on the solution and the results were used to guide the three-dimensional calculations. Three-dimensional calculations of a Mach 1.4 jet showed that this LES simulation accurately captures the physics of the turbulent flow. The agreement with experimental data was relatively good and is much better than results in the current literature. Turbulent intensities indicate that the turbulent structures at this level of modeling are not isotropic and this information could lend itself to the development of improved subgrid scale models for LES and turbulence models for RANS simulations. A two point correlation technique was used to quantify the turbulent structures. Two point space correlations were used to obtain a measure of the integral length scale, which proved to be approximately 1/2 D(sub j). Two point space-time correlations were used to obtain the convection velocity for the turbulent structures. This velocity ranged from 0.57 to 0.71 U(sub j).

DeBonis, James R.↗
DOE OSTI2020

Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry

Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code XGC, in order to study density-pump out and electron heat confinement.The RMP field is imported from the extended magneto-hydrodynamics (MHD) code M3D-C1, taking into account the linear two-fluid plasma response.With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate, and suppression of the electron heat flux in the steepest part of the edge pedestal.In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.

DIII-D↗
DOE OSTI2023

Direct measurement of the electron turbulence-broadening edge transport barrier to facilitate core–edge integration in tokamak fusion plasmas

Abstract The integration of a high-performance core and a dissipative divertor, or the so-called ‘core–edge integration,’ has been widely identified as a critical gap in the design of future fusion reactors. In this letter, we report, for the first time, direct experimental evidence of electron turbulence at the DIII-D H-mode pedestal that correlates with the broadening of the pedestal and thus facilitates core–edge integration. In agreement with gyrokinetic simulations, this electron turbulence is enhanced by high η e ( η e = L n / L T e , where L n is the density scale length and L T e is the electron temperature scale length), which is due to a strong shift between the density and temperature pedestal profiles associated with a closed divertor. The modeled turbulence drives significant heat transport with a lower pressure gradient that may broaden the pedestal to a greater degree than the empirical and theoretically predicted pedestal width scalings. Such a wide pedestal, coupled with a closed divertor, enables us to achieve a good core–edge scenario that integrates a high-temperature low-collisionality pedestal (pedestal top temperature T e,ped > 0.8 keV and a pedestal top collisionality ν * ped < 1) under detached divertor conditions. This paves a new path toward solving the core–edge integration issue in future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Turbulent flow in a channel with a wall with progressive waves

The effects of prescribed wall motion on turbulent channel flow were examined with the objective of understanding the drag mechanisms and obtaining a possible means for drag reduction. A computer program for large-eddy simulation (LES) of turbulent flow in a channel was adapted to treat the case of time-varying wall motion using linearized boundary conditions. The code was applied to flow in a channel with prescribed wall motion on one wall. It was found that the results can be explained by a model based on a Stokes-like layer near the moving wall. This provides a theory which can be used to guide the selection of parameters. Considering both viscous and pressure drag components, it was found that net drag reduction is possible under certain conditions of wall motion. However, the work required to drive the wall motion can be greater than the reduced through-flow energy requirement. The major effects of the wall motion are restricted to the viscous sublayer of the turbulent flow.

Kuhn, G. D.↗
NASA NTRS1992

The structure of intense vorticity in homogeneous isotropic turbulence

The structure of the intense vorticity regions is studied in numerically simulated homogeneous, isotropic, equilibrium turbulent flow fields at four different Reynolds numbers in the range Re(sub lambda) = 36-171. In accordance with previous investigators, this vorticity is found to be organized in coherent, cylindrical or ribbon-like, vortices ('worms'). A statistical study suggests that they are just especially intense features of the background, O(omega'), vorticity. Their radii scale with the Kolmogorov microscale and their lengths with the integral scale of the flow. An interesting observation is that the Reynolds number based on the circulation of the intense vortices, gamma/nu, increases monotonically with Re(sub lambda), raising the question of the stability of the structures in the limit of Re(sub lambda) approaching infinity. One and two-dimensional statistics of vorticity and strain are presented; they are non-gaussian, and the behavior of their tails depends strongly on the Reynolds number. There is no evidence of convergence to a limiting distribution in our range of Re(sub lambda), even though the energy spectra and the energy dissipation rate show good asymptotic properties in the higher Reynolds number cases. Evidence is presented to show that worms are natural features of the flow and that they do not depend on the particular forcing scheme.

Jimenez, J.↗
NASA NTRS1994

Validation of two-equation turbulence models for propulsion flowfields

The objective of the study is to assess the capability of two-equation turbulence models for simulating propulsion-related flowfields. The standard kappa-epsilon model with Chien's low Reynolds number formulation for near-wall effects is used as the baseline turbulence model. Several experimental test cases, representative of rocket combustor internal flowfields, are used to catalog the performance of the baseline model. Specific flowfields considered here include recirculating flow behind a backstep, mixing between coaxial jets and planar shear layers. Since turbulence solutions are notoriously dependent on grid and numerical methodology, the effects of grid refinement and artificial dissipation on numerical accuracy are studied. In the latter instance, computational results obtained with several central-differenced and upwind-based formulations are compared. Based on these results, improved turbulence modes such as enhanced kappa-epsilon models as well as other two-equation formulations (e.g., kappa-omega) are being studied. In addition, validation of swirling and reacting flowfields are also currently underway.

Deshpande, Manish↗