Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large-eddy”

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 109 records · Page 6

The transition prediction toolkit: LST, SIT, PSE, DNS, and LES

The e(sup N) method for predicting transition onset is an amplitude ratio criterion that is on the verge of full maturation for three-dimensional, compressible, real gas flows. Many of the components for a more sophisticated, absolute amplitude criterion are now emerging: receptivity theory, secondary instability theory, parabolized stability equations approaches, direct numerical simulation and large-eddy simulation. This paper will provide a description of each of these new theoretical tools and provide indications of their current status.

Zang, Thomas A.↗

A second-order bulk boundary-layer model

Bulk mass-flux models represent the large eddies that are primarily responsible for the turbulent fluxes in the planetary boundary layer as convective circulations, with an associated convective mass flux. In order for such models to be useful, it is necessary to determine the fractional area covered by rising motion in the convective circulations. This fraction can be used as an estimate of the cloud amount, under certain conditions. 'Matching' conditions have been developed that relate the convective mass flux to the ventilation and entrainment mass fluxes. These are based on conservation equations for the scalar means and variances in the entrainment and ventilation layers. Methods are presented to determine both the fractional area covered by rising motion and the convective mass flux. The requirement of variance balance is used to relax the 'well-mixed' assumption. The vertical structures of the mean state and the turbulent fluxes are determined analytically. Several aspects of this simple model's formulation are evaluated using results from large-eddy simulations.

Randall, David A.↗

Direct simulation of the stably stratified turbulent Ekman layer

The Navier-Stokes equations and the Boussinesq approximation were used to compute a 3D time-dependent turbulent flow in the stably stratified Ekman layer over a smooth surface. The simulation data are found to be in very good agreement with atmospheric measurements when nondimensionalized according to Nieuwstadt's local scaling scheme. Results suggest that, when Reynolds number effects are taken into account, the 'constant Froud number' stable layer model (Brost and Wyngaard, 1978) and the 'shearing length' stable layer model (Hunt, 1985) for the dissipitation rate of turbulent kinetic energy are both valid. It is concluded that there is good agreement between the direct numerical simulation results and large-eddy simulation results obtained by Mason and Derbyshire (1990).

Coleman, G. N.↗

Physical processes within the nocturnal stratus-topped boundary layer

There are many physical processes involved in the stratus-topped boundary layer: longwave radiation cooling, entrainment, latent heating, surface heating, solar heating, and drizzling, for example. The manner in which the processes combine to maintain the turbulence within the stratus-topped boundary layer remains an unsolved problem. The large eddy simulation technique is used to examine the first four physical processes mentioned above. First, the contribution of each physical process to the thermodynamic differences between the updraft and downdraft branches of turbulent circulations is examined through a conditional sampling. Second, these mean thermodynamic differences are shown to express well the vertical distributions of heat and moisture fluxes within stratus-topped boundary layers. These provide a method to validate the process partitioning technique. (This technique assumes that the net flux profile can be partitioned into different component-flux profiles according to physical processes, and that each partitioned component flux is linear in height.) In this paper, the heat and moisture fluxes are process partitioned, and each component-flux is found to contribute to the net flux in a way that is consistent with its corresponding process contribution to the mean thermodynamic differences between updrafts and downdrafts. Also, the net flux obtained by summing all component-fluxes agrees very well with that obtained directly from the large-eddy simulations.

Moeng, Chin-Hoh↗

The structure of supersonic jet flow and its radiated sound

Large-eddy simulation of a supersonic jet is presented with emphasis on capturing the unsteady features of the flow pertinent to sound emission. A high-accuracy numerical scheme is used to solve the filtered, unsteady, compressible Navier-Stokes equations while modelling the subgrid-scale turbulence. For random inflow disturbance, the wave-like feature of the large-scale structure is demonstrated. The large-scale structure was then enhanced by imposing harmonic disturbances to the inflow. The limitation of using the full Navier-Stokes equation to calculate the far-field sound is discussed. Application of Lighthill's acoustic analogy is given with the objective of highlighting the difficulties that arise from the non-compactness of the source term.

Mankbadi, Reda R.↗

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↗

Simulations of curved turbulent boundary layers

The objective of this work is to develop a space-time accurate numerical method for the solution of incompressible Navier-Stokes equations in generalized coordinates. The resulting code is to be used for direct and large-eddy simulation of turbulence in complex geometries. In a previous paper, the system of Navier-Stokes equations in general curvilinear coordinates was solved by a second-order accurate finite-difference scheme. Satisfactory results were obtained for several flows in two and three dimensions. The system of Navier-Stokes for the fluxes are given in Orlandi (1989). The main deficiency of the numerical scheme was the large CPU time required for the solution of the Poisson equation for the 'pressure' field. The point SOR relaxation, in conjunction with a multigrid scheme, was used for the Poisson equation. In some cases, particularly with very fine grids, it was impossible to obtain a divergent-free flow. A preliminary attempt is made to compute the spatially evolving flow of Swearingen & Blackwelder. To reduce the streamwise distance, the inflow was at a distance x = 60 cm from the leading edge.

Orlandi, Paolo↗

Numerical simulations of turbulent thermal convection with differential rotation

The solar nebula, from which the planets in our solar system formed, featured a disk of gas and dust grains in rapid, differential rotation, and at some stage was likely to have been unstable to thermal convection. This situation is suspected by many to lead to significant turbulent Reynolds stress production and angular momentum transport in such systems, and estimates of transport rates have been attempted from unsubstantiated phenomenological models. In order to determine the circumstances and physical conditions under which our own planetary system formed and to explain recent observations of young stellar systems, it is necessary to develop realistic models of heat and angular momentum transport for such flows. Developing an understanding of complicated flows featuring thermal convection, rotation, and shear is also of wide interest in stellar astrophysics and in planetary and terrestrial atmospheric studies. The ultimate objective is to develop workable models based on the numerical simulations for constructing global solar nebula models; viz., relatively simple prescriptions for heat and angular momentum fluxes from given system parameters (e.g., ratios of rotation, shear, and convective lapse rates) are characterized, quantified, and developed. Toward this end, our program has been to attempt to understand the behavior of the direct numerical simulations of Boussinesq convection, which, despite the complexity of the results, is still an overly simplified approximation to the real system and should be more amenable to analysis. These results are also intended to be tested against turbulence models, especially those designed for atmospheric boundary layers, and may provide a basis for subgrid-scale models. In order to make the numerical simulations more realistic with regard to the solar nebula problem, a fully compressible code that will allow incorporation of large density stratifications and realistic thermodynamic and radiative properties is developed. In order to explore the properties of these flows at the very high values of Re found in natural systems and the very low values of Pr found in most astrophysical contexts, we will need to employ large-eddy simulations for which we want to determine the most appropriate subgrid-scale model to incorporate.

Cabot, William H.↗

Similarity states of homogeneous stably-stratified turbulence at infinite Froude number

We present evidence of similarity states which may develop inhomogeneous stably-stratified flows if a dimensionless group in addition to the Reynolds number, the so-called Froude number, is sufficiently large. Here, we define the Froude number as the ratio of the internal wave time-scale to the turbulence time-scale. We examine three different similarity states which may develop depending on the initial conditions of the velocity and density fields. Theoretical arguments and results of large-eddy simulations will be presented. We will conclude this report with some speculative thoughts on similarity states which may develop in stably-stratified turbulence at arbitrary Froude number as well as our future research plans in this area.

Chasnov, Jeffrey R.↗

On the local nature of the energy cascade

The local nature of the energy cascade in space and time is studied using direct numerical simulation of decaying and forced isotropic turbulence. To examine the concept that large scales evolve into smaller ones, we compute the Lagrangian correlation coefficient between local kinetic energy at different scales. This correlation is found to peak at a Lagrangian time-delay that increases with scale separation. The results show that, on average, the flow of energy to smaller scales is predominantly local in physical space and that the view of eddies decaying into smaller ones while transferring their kinetic energy appears to be, on average, quite realistic. To examine the spectral characteristics of the cascade under unsteady conditions, a pulse of large-scale energy is added to the large-eddy simulation of forced isotropic turbulence. As time progresses, the evolution of this pulse through bands of increasing wavenumbers is studied.

Meneveau, C.↗

Parallel spatial direct numerical simulations on the Intel iPSC/860 hypercube

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube is documented. The direct numerical simulation approach is used to compute spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows. The feasibility of using the PSDNS on the hypercube to perform transition studies is examined. The results indicate that the direct numerical simulation approach can effectively be parallelized on a distributed-memory parallel machine. By increasing the number of processors nearly ideal linear speedups are achieved with nonoptimized routines; slower than linear speedups are achieved with optimized (machine dependent library) routines. This slower than linear speedup results because the Fast Fourier Transform (FFT) routine dominates the computational cost and because the routine indicates less than ideal speedups. However with the machine-dependent routines the total computational cost decreases by a factor of 4 to 5 compared with standard FORTRAN routines. The computational cost increases linearly with spanwise wall-normal and streamwise grid refinements. The hypercube with 32 processors was estimated to require approximately twice the amount of Cray supercomputer single processor time to complete a comparable simulation; however it is estimated that a subgrid-scale model which reduces the required number of grid points and becomes a large-eddy simulation (PSLES) would reduce the computational cost and memory requirements by a factor of 10 over the PSDNS. This PSLES implementation would enable transition simulations on the hypercube at a reasonable computational cost.

Joslin, Ronald D.↗

Computational aero-acoustics as a branch of turbulence research

Computational aero-acoustics (CAA) is concerned with calculations of the aerodynamically generated sound source and its propagation. This paper presents a background to CAA, prediction of the time-dependent sound source using large-eddy simulations, and sound propagation to the far-field.

Mankbadi, Reda R.↗

Local isotropy in distorted turbulent boundary layers at high Reynolds number

This is a report on the continuation of our experimental investigations of the hypothesis of local isotropy in shear flows. This hypothesis, which states that at sufficiently high Reynolds numbers the small-scale structures of turbulent motions are independent of large-scale structures and mean deformations, has been used in theoretical studies of turbulence and computational methods such as large-eddy simulation. Since Kolmogorov proposed his theory, there have been many experiments, conducted in wakes, jets, mixing layers, a tidal channel, and atmospheric and laboratory boundary layers, in which attempts have been made to verify - or refute - the local-isotropy hypothesis. However, a review of the literature over the last five decades indicated that, despite all these experiments in shear flows, there was no consensus in the scientific community regarding this hypothesis, and, therefore, it seemed worthwhile to undertake a fresh experimental investigation into this question.

Saddoughi, Seyed G.↗

Control of low-speed turbulent separated flow over a backward-facing ramp

The relative performance and flow phenomena associated with several devices for controlling turbulent separated flow were investigated at low speeds. Relative performance of the devices was examined for flow over a curved, backward-facing ramp in a wind tunnel, and the flow phenomena were examined in a water tunnel using dye-flow visualization. Surface static pressure measurements and oil-flow visualization results from the wind tunnel tests indicated that transverse grooves, longitudinal grooves, submerged vortex generators, vortex generator jets (VGJ's), Viets' fluidic flappers, elongated arches at positive angle of attack, and large-eddy breakup devices (LEBU's) at positive angle of attack placed near the baseline separation location reduce flow separation and increase pressure recovery. Spanwise cylinders reduce flow separation but decrease pressure recovery downstream. Riblets, passive porous surfaces, swept grooves, Helmholtz resonators, and arches and LEBU's with angle of attack less than or = 0 degrees had no significant effect in reducing the extent of the separation region. Wall-cooling computations indicated that separation delay on a partially-cooled ramp is nearly the same as on a fully-cooled ramp, while minimizing the frictional drag increase associated with the wall cooling process. Dry-flow visualization tests in the water tunnel indicated that wishbone vortex generators in the forward orientation shed horseshoe vortices; wishbone vortex generators oriented in the reverse direction and doublet vortex generators shed streamwise counterrotating vortices; a spanewise cylinder located near the wall and LEBU's at angle of attack = -10 degrees produced eddies or transverse vortices which rotated with the same sign as the mean vorticity in a turbulent boundary layer; and the most effective VGJ's produced streamwise co-rotating vortices. Comparative wind-tunnel test results indicated that transferring momentum from the outer region of a turbulent boundary layer through the action of embedded streamwise vortices is more effective than by transverse vortices for the separation control application studied herein.

Lin, John C.↗

Lumley's PODT definition of large eddies and a trio of numerical procedures

Lumley's 1967 Moscow paper provided, for the first time, a completely rational definition of the physically-useful term 'large eddy', popular for a half-century. The numerical procedures based upon his results are: (1) PODT (Proper Orthogonal Decomposition Theorem), which extracts the Large Eddy structure of stochastic processes from physical or computer simulation two-point covariances, and 2) LEIM (Large-Eddy Interaction Model), a predictive scheme for the dynamical large eddies based upon higher order turbulence modeling. Earlier Lumley's work (1964) forms the basis for the final member of the triad of numerical procedures: this predicts the global neutral modes of turbulence which have surprising agreement with both structural eigenmodes and those obtained from the dynamical equations. The ultimate goal of improved engineering design tools for turbulence may be near at hand, partly due to the power and storage of 'supermicrocomputer' workstations finally becoming adequate for the demanding numerics of these procedures.

Payne, Fred R.↗

Identification of modeled ocean plumes in Greenland gyre ERS-1 SAR data

Oceanic convective plumes modeled with a thermobaric large-eddy simulation and driven by conditions similar to those of the Greenland Sea are compared to observations from ERS-1 Synthetic Aperture Radar (SAR) data from the Greenland Sea for the winter of 1992. In both form and size the two representations are seen to compare favorably. The plume-filled area of the SAR image occupies a region about 20 km by 90 km at the ice edge of the open water in 'Nordbukta', the large seasonal ice retreat, in the 'Odden' ice protuberance in the southern Greenland gyre. In the SAR data the plumes appear to be ice covered while the convective-return areas are open.

Carsey, Frank D.↗

Micro Blowing Simulations Using a Coupled Finite-Volume Lattice-Boltzman n L ES Approach

Three dimensional large-eddy simulations (LES) of single and multiple jet-in-cross-flow (JICF) are conducted using the 19-bit Lattice Boltzmann Equation (LBE) method coupled with a conventional finite-volume (FV) scheme. In this coupled LBE-FV approach, the LBE-LES is employed to simulate the flow inside the jet nozzles while the FV-LES is used to simulate the crossflow. The key application area is the use of this technique is to study the micro blowing technique (MBT) for drag control similar to the recent experiments at NASA/GRC. It is necessary to resolve the flow inside the micro-blowing and suction holes with high resolution without being restricted by the FV time-step restriction. The coupled LBE-FV-LES approach achieves this objectives in a computationally efficient manner. A single jet in crossflow case is used for validation purpose and the results are compared with experimental data and full LBE-LES simulation. Good agreement with data is obtained. Subsequently, MBT over a flat plate with porosity of 25% is simulated using 9 jets in a compressible cross flow at a Mach number of 0.4. It is shown that MBT suppresses the near-wall vortices and reduces the skin friction by up to 50 percent. This is in good agreement with experimental data.

Menon, S.↗

The NASA Turbulent Heat Flux Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds- averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗