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 541 records · Page 30

Large Eddy Simulations of Severe Convection Induced Turbulence

Convective storms can pose a serious risk to aviation operations since they are often accompanied by turbulence, heavy rain, hail, icing, lightning, strong winds, and poor visibility. They can cause major delays in air traffic due to the re-routing of flights, and by disrupting operations at the airports in the vicinity of the storm system. In this study, the Terminal Area Simulation System is used to simulate five different convective events ranging from a mesoscale convective complex to isolated storms. The occurrence of convection induced turbulence is analyzed from these simulations. The validation of model results with the radar data and other observations is reported and an aircraft-centric turbulence hazard metric calculated for each case is discussed. The turbulence analysis showed that large pockets of significant turbulence hazard can be found in regions of low radar reflectivity. Moderate and severe turbulence was often found in building cumulus turrets and overshooting tops.

Ahmad, Nash'at↗

Modeling of compressible turbulent shear flows

Despite all the recent developments in computer technologies and numerical algorithms, full numerical simulations of turbulent flows are feasible only at moderate Reynolds numbers and for flows with relatively simple geometries. The main goal of this research is to develop new second order moment closures for compressible turbulence. It has been shown that the models based on the extension of those developed originally for incompressible flows fail to adequately predict turbulent flows at high Mach numbers. In this attempt, the compressibility effects are explicitly considered. A successful development of these models that directly takes into account the compressibility effects may have a range of technological implications in the design of supersonic and hypersonic vehicles.

Liou, William W.↗

Turbulent Convection: Is 2D a good proxy of 3D?

Several authors have recently carried out 2D simulations of turbulent convection for both solar and massive stars. Fitting the 2D results with the MLT, they obtain that alpha(sub MLT) greater than 1 specifically, 1.4 less than alpha(sub MLT) less than 1.8. The authors further suggest that this methodology could be used to calibrate the MLT used in stellar evolutionary codes. We suggest the opposite viewpoint: the 2D results show that MLT is internally inconsistent because the resulting alpha(sub MLT) greater than 1 violates the MLT basic assumption that alpha(sub MLT) less than 1. When the 2D results are fitted with the CM model, alpha(sub CMT) less than 1, in accord with the basic tenet of the model. On the other hand, since both MLT and CM are local models, they should be replaced by the next generation of non-local, time dependent turbulence models which we discuss in some detail.

Canuto, V. M.↗

Numerical Simulation of High-Speed Turbulent Reacting Flows

The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES) methods for computational analysis of high-speed reacting turbulent flows. We have just completed the first year of Phase 3 of this research.

Jaberi, F. A.↗

Analysis of discretization errors in LES

All numerical simulations of turbulence (DNS or LES) involve some discretization errors. The integrity of such simulations therefore depend on our ability to quantify and control such errors. In the classical literature on analysis of errors in partial differential equations, one typically studies simple linear equations (such as the wave equation or Laplace's equation). The qualitative insight gained from studying such simple situations is then used to design numerical methods for more complex problems such as the Navier-Stokes equations. Though such an approach may seem reasonable as a first approximation, it should be recognized that strongly nonlinear problems, such as turbulence, have a feature that is absent in linear problems. This feature is the simultaneous presence of a continuum of space and time scales. Thus, in an analysis of errors in the one dimensional wave equation, one may, without loss of generality, rescale the equations so that the dependent variable is always of order unity. This is not possible in the turbulence problem since the amplitudes of the Fourier modes of the velocity field have a continuous distribution. The objective of the present research is to provide some quantitative measures of numerical errors in such situations. Though the focus of this work is LES, the methods introduced here can be just as easily applied to DNS. Errors due to discretization of the time-variable are neglected for the purpose of this analysis.

Ghosal, Sandip↗

Development of an adaptive failure detection and identification system for detecting aircraft control element failures

A methodology for designing a failure detection and identification (FDI) system to detect and isolate control element failures in aircraft control systems is reviewed. An FDI system design for a modified B-737 aircraft resulting from this methodology is also reviewed, and the results of evaluating this system via simulation are presented. The FDI system performed well in a no-turbulence environment, but it experienced an unacceptable number of false alarms in atmospheric turbulence. An adaptive FDI system, which adjusts thresholds and other system parameters based on the estimated turbulence level, was developed and evaluated. The adaptive system performed well over all turbulence levels simulated, reliably detecting all but the smallest magnitude partially-missing-surface failures.

Bundick, W. Thomas↗

Comparison of CFD Predictions with Shuttle Global Flight Thermal Imagery and Discrete Surface Measurements

Surface temperature measurements from the STS-119 boundary-layer transition experiment on the space shuttle orbiter Discovery provide a rare opportunity to assess turbulent CFD models at hypersonic flight conditions. This flight data was acquired by on-board thermocouples and by infrared images taken off-board by the Hypersonic Thermodynamic Infrared Measurements (HYTHIRM) team, and is suitable for hypersonic CFD turbulence assessment between Mach 6 and 14. The primary assessment is for the Baldwin-Lomax and Cebeci-Smith algebraic turbulence models in the DPLR and LAURA CFD codes, respectively. A secondary assessment is made of the Shear-Stress Transport (SST) two-equation turbulence model in the DPLR code. Based upon surface temperature comparisons at eleven thermocouple locations, the algebraic-model turbulent CFD results average 4% lower than the measurements for Mach numbers less than 11. For Mach numbers greater than 11, the algebraic-model turbulent CFD results average 5% higher than the three available thermocouple measurements. Surface temperature predictions from the two SST cases were consistently 3 4% higher than the algebraic-model results. The thermocouple temperatures exhibit a change in trend with Mach number at about Mach 11; this trend is not reflected in the CFD results. Because the temperature trends from the turbulent CFD simulations and the flight data diverge above Mach 11, extrapolation of the turbulent CFD accuracy to higher Mach numbers is not recommended.

Wood, William A.↗

Numerical Dissipation and Wrong Propagation Speed of Discontinuities for Stiff Source Terms

In compressible turbulent combustion/nonequilibrium flows, the constructions of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and (b) obtaining correct propagation speed of discontinuities for stiff reacting terms on coarse grids share one important ingredient - minimization of numerical dissipation while maintaining numerical stability. Here coarse grids means standard mesh density requirement for accurate simulation of typical non-reacting flows. This dual requirement to achieve both numerical stability and accuracy with zero or minimal use of numerical dissipation is most often conflicting for existing schemes that were designed for non-reacting flows. The goal of this paper is to relate numerical dissipations that are inherited in a selected set of high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities for two representative stiff detonation wave problems.

Yee, H. C.↗

Numerical Dissipation and Wrong Propagation Speed of Discontinuities for Stiff Source Terms

In compressible turbulent combustion/nonequilibrium flows, the constructions of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and (b) obtaining correct propagation speed of discontinuities for stiff reacting terms on coarse grids share one important ingredient - minimization of numerical dissipation while maintaining numerical stability. Here coarse grids means standard mesh density requirement for accurate simulation of typical non-reacting flows. This dual requirement to achieve both numerical stability and accuracy with zero or minimal use of numerical dissipation is most often conflicting for existing schemes that were designed for non-reacting flows. The goal of this paper is to relate numerical dissipations that are inherited in a selected set of high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities as a function of stiffness of the source term and the grid spacing.

Yee, H. C.↗

Direct numerical simulation of reacting flows

The objectives of this work are: (1) to extend the technique of direct numerical simulations to turbulent, chemically reacting flows, (2) to test the validity of the method by comparing computational results with laboratory data, and (3) to use the simulations to gain a better understanding of the effects of turbulence on chemical reactions. The effects of both the large scale structure and the smaller scale turbulence on the overall reaction rates are addressed. The relationship between infinite reaction rate and finite reaction rate chemistry is compared with some of the results of calculations with existing theories and laboratory data. The direct numerical simulation method involves the numerical solution of the detailed evolution of the complex turbulent velocity and concentration fields. Using very efficient numerical methods (e.g., pseudospectral methods), the fully nonlinear (possibly low pass filtered) equations of motion are solved and no closure assumptions or turbulence models are used. Statistical data are obtained by performing spatial, temporal, and/or ensemble averages over the computed flow fields.

Riley, J. J.↗

Turbulent flow in a 180 deg bend: Modeling and computations

A low Reynolds number k-epsilon turbulence model was presented which yields accurate predictions of the kinetic energy near the wall. The model is validated with the experimental channel flow data of Kreplin and Eckelmann. The predictions are also compared with earlier results from direct simulation of turbulent channel flow. The model is especially useful for internal flows where the inflow boundary condition of epsilon is not easily prescribed. The model partly derives from some observations based on earlier direct simulation results of near-wall turbulence. The low Reynolds number turbulence model together with an existing curvature correction appropriate to spinning cylinder flows was used to simulate the flow in a U-bend with the same radius of curvature as the Space Shuttle Main Engine (SSME) Turn-Around Duct (TAD). The present computations indicate a space varying curvature correction parameter as opposed to a constant parameter as used in the spinning cylinder flows. Comparison with limited available experimental data is made. The comparison is favorable, but detailed experimental data is needed to further improve the curvature model.

Kaul, Upender K.↗

Collaborative testing of turbulence models

A review is given of an ongoing international project, in which data from experiments on, and simulations of, turbulent flows are distributed to developers of (time-averaged) engineering turbulence models. The predictions of each model are sent to the organizers and redistributed to all the modelers, plus some experimentalists and other experts (total approx. 120), for comment. The 'reaction time' of modelers has proved to be much longer than anticipated, partly because the comparisons with data have prompted many modelers to improve their models or numerics.

Bradshaw, Peter↗

Turbulence in a Global Magnetohydrodynamic Simulation of the Earth's Magnetosphere during Northward and Southward Interplanetary Magnetic Field

We report the results of MHD simulations of Earth's magnetosphere for idealized steady solar wind plasma and interplanetary magnetic field (IMF) conditions. The simulations feature purely northward and southward magnetic fields and were designed to study turbulence in the magnetotail plasma sheet. We found that the power spectral densities (PSDs) for both northward and southward IMF had the characteristics of turbulent flow. In both cases, the PSDs showed the three scale ranges expected from theory: the energy-containing scale, the inertial range, and the dissipative range. The results were generally consistent with in-situ observations and theoretical predictions. While the two cases studied, northward and southward IMF, had some similar characteristics, there were significant differences as well. For southward IMF, localized reconnection was the main energy source for the turbulence. For northward IMF, remnant reconnection contributed to driving the turbulence. Boundary waves may also have contributed. In both cases, the PSD slopes had spatial distributions in the dissipative range that reflected the pattern of resistive dissipation. For southward IMF there was a trend toward steeper slopes in the dissipative range with distance down the tail. For northward IMF there was a marked dusk-dawn asymmetry with steeper slopes on the dusk side of the tail. The inertial scale PSDs had a dusk-dawn symmetry during the northward IMF interval with steeper slopes on the dawn side. This asymmetry was not found in the distribution of inertial range slopes for southward IMF. The inertial range PSD slopes were clustered around values close to the theoretical expectation for both northward and southward IMF. In the dissipative range, however, the slopes were broadly distributed and the median values were significantly different, consistent with a different distribution of resistivity.

magnetohydrodynamics (MHD)↗

Evaluation of the Use of Second Generation Wavelets in the Coherent Vortex Simulation Approach

The objective of this study is to investigate the use of the second generation bi-orthogonal wavelet transform for the field decomposition in the Coherent Vortex Simulation of turbulent flows. The performances of the bi-orthogonal second generation wavelet transform and the orthogonal wavelet transform using Daubechies wavelets with the same number of vanishing moments are compared in a priori tests using a spectral direct numerical simulation (DNS) database of isotropic turbulence fields: 256(exp 3) and 512(exp 3) DNS of forced homogeneous turbulence (Re(sub lambda) = 168) and 256(exp 3) and 512(exp 3) DNS of decaying homogeneous turbulence (Re(sub lambda) = 55). It is found that bi-orthogonal second generation wavelets can be used for coherent vortex extraction. The results of a priori tests indicate that second generation wavelets have better compression and the residual field is closer to Gaussian. However, it was found that the use of second generation wavelets results in an integral length scale for the incoherent part that is larger than that derived from orthogonal wavelets. A way of dealing with this difficulty is suggested.

Goldstein, D. E.↗

Modeling Laminar-to-Turbulent Transition in the Panel Test Facility Arcjet

A computational fluid dynamics (CFD) laminar-to-turbulence transition model was developed for the NASA Ames Research Center’s 20 MW Panel Test Facility (PTF). Surface pressure, heat flux to a water-cooled plate, and surface temperature on a tile plate coated with reaction-cured glass were measured across several conditions in the facility and compared with laminar and fully-turbulent CFD simulations. The potential for bypass transition in the PTF nozzle was assessed via application of the Langtry-Menter four-equation transitional shear-stress transport (SST) model. Results from the transition model were inconsistent with measurements. Flow interaction with a boundary conditioning plate feature inside the nozzle was also investigated as a potential source of laminar-to-turbulent transition using two turbulence models with specified transition locations. The Baldwin-Lomax turbulence model was configured to simulate a transition at the upstream edge of the boundary conditioning plate and produced results consistent with the surface pressure measurements but not the cold-wall heat flux. Finally, the SST turbulence model was calibrated to transition at the upstream edge of the boundary conditioning plate and produced results consistent with both the surface pressure and cold-wall heat flux measurements. The SST-based model demonstrated reasonable agreement with surface temperature measurements on the reaction-cured glass tile, albeit with some discrepancies.

Computational Fluid Dynamics↗

Transport of passive scalars in a turbulent channel flow

A direct numerical simulation of a turbulent channel flow with three passive scalars at different molecular Prandtl numbers is performed. Computed statistics including the turbulent Prandtl numbers are compared with existing experimental data. The computed fields are also examined to investigate the spatial structure of the scalar fields. The scalar fields are highly correlated with the streamwise velocity; the correlation coefficient between the temperature and the streamwise velocity is as high as 0.95 in the wall region. The joint probability distributions between the temperature and velocity fluctuations are also examined; they suggest that it might be possible to model the scalar fluxes in the wall region in a manner similar to the Reynolds stresses.

Kim, John↗