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 757 records · Page 42

Implicit Large-Eddy Simulation of Single-Injector Cooling Flow

The turbulent heat flux vector plays a critical role in aero-propulsion flow. Currently, traditional RANS based methods do not currently have reliable closures which can accurately predict the turbulent transport of heat. Scale Resolving Simulations can provide critical insight into this fundamental fluid mechanism. Large-eddy Simulations (LES) were performed on the single-hole film-cooling experimental configuration from the recent experiments conducted at NASA Glenn Research Center (GRC)[1]. The experiment of Wernet et al. obtained mean and fluctuating temperatures which can be used to validate LES results predicting the turbulent heat flux vector. The FDL3DI flow solver was used to complete simulations of two blowing ratios of the THX - III dataset. A grid sensitivity study was also performed. Initial LES results show promising agreement between the simulation and experimental measurements.

High-Order methods↗

High speed turbulent reacting flows: DNS and LES

Work on understanding the mechanisms of mixing and reaction in high speed turbulent reacting flows was continued. Efforts, in particular, were concentrated on taking advantage of modern computational methods to simulate high speed turbulent flows. In doing so, two methodologies were used: large eddy simulations (LES) and direct numerical simulations (DNS). In the work related with LES the objective is to study the behavior of the probability density functions (pdfs) of scalar properties within the subgrid in reacting turbulent flows. The data base obtained by DNS for a detailed study of the pdf characteristics within the subgrid was used. Simulations are performed for flows under various initializations to include the effects of compressibility on mixing and chemical reactions. In the work related with DNS, a two-dimensional temporally developing high speed mixing layer under the influence of a second-order non-equilibrium chemical reaction of the type A + B yields products + heat was considered. Simulations were performed with different magnitudes of the convective Mach numbers and with different chemical kinetic parameters for the purpose of examining the isolated effects of the compressibility and the heat released by the chemical reactions on the structure of the layer. A full compressible code was developed and utilized, so that the coupling between mixing and chemical reactions is captured in a realistic manner.

Givi, Peyman↗

Nonnormality of Sonic Boom Loudness Metrics in the Turbulent Atmospheric Boundary Layer at Large Lateral Distances from the Flight Path

Atmospheric boundary layer (ABL) turbulence causes variability of the sonic boom waveform at the ground. Recent numerical investigations of sonic boom propagation through kinematic velocity fluctuations indicate that loudness metric distributions are positively skewed relative to a normal distribution. This skewness depends on the propagation distance and turbulence intensity. Propagation simulations of N-waves and shaped booms through inhomogeneous ABL turbulence are presented. Meteorological conditions are varied to examine different daytime ABL conditions and their effect on sonic boom loudness distributions. Two outcomes are observed: 1) the loudness metric distributions become increasingly positively skewed as the propagation distance through the ABL increases, and 2) the distributions become increasingly positively skewed at the same lateral distance from the flight path as the convection level of the daytime ABL is increased. Thus, results indicate that ground level measurements of sonic boom loudness from flight tests performed at large lateral distances from the flight path may not be normally distributed, due to turbulence present in the ABL. (This research is supported by the Commercial Supersonic Technology Project of the National Aeronautics and Space Administration under Grant No. 80NSSC19K1685.)

sonic boom↗

Direct simulation of particle dispersion in a decaying isotropic turbulence

Results of a numerical investigation of the dispersion of solid particles in decaying isotropic turbulence are presented. The 3D time-dependent velocity field of a homogeneous nonstationary turbulence is computed using the method of direct numerical simulation (DNS). The dispersion characteristics of three different solid particles (corn, copper, and glass) injected in the flow are obtained by integrating the complete equation of particle motion along the instantaneous trajectories of 22-cubed particles for each particle type, and then by performing ensemble averaging. Good agreement was achieved between the present DNS results and the measured time development of the mean-square displacement of the particles. Questions of how and why the dispersion statistics of a solid particle differ from those of its corresponding fluid point and surrounding fluid and what influences inertia and gravity have on these statistics are also discussed.

Elghobashi, S.↗

Assessment of UMUSCL Scheme for DNS of Turbulent Flows

Direct numerical simulation (DNS) is performed using an unstructured node-based finite-volume methodology for canonical turbulent flow configurations to assess the effect of the upwinding parameter (κ) in the unstructured monotonic upstream scheme for conservation laws (UMUSCL) reconstruction scheme. The NASA FUN3D solver is used, and four turbulent validation cases are considered: periodic Taylor-Green vortex, (2) flow through a plane channel, (3) flow through a channel with a constriction, and (4) flow over a flat plate. The results are presented for κ = 0.5, 0.9 and 1 on hexahedral grids, with the last value corresponding to a zero-numerical dissipation scheme. Detailed comparisons with available reference data indicate that for the smallest length-scales simulated, κ = 0.5 is too dissipative while κ = 1 could lead to an energy build-up. Overall, κ =0.9 yields satisfactory results in terms of accuracy and robustness for all the cases considered in this study.

Direct Numerical Simulation↗

Direct and large-eddy simulations of three-dimensional compressible Navier-Stokes turbulence

This paper reports results from the numerical implementation and testing of the compressible large eddy-simulation (LES). Relevant quantities from 32-cubed coarse grid LES solutions are compared with results generated from direct numerical simulations (DNS) of three-dimensional compressible turbulence that have been run both with sufficient resolution, at 96 cubed. The 32 cubed LES results overall agree well with their 96 cubed DNS counterparts. Moreover, the new DNS results confirm several recent conclusions about compressible turbulence that have been based primarily on two-dimensional simulations.

Zang, T. A.↗

Helicity fluctuations and turbulent energy production in rotating and non-rotating pipes

Finite-difference second-order accurate direct simulation of a turbulent pipe has been used to investigate how the turbulence production and dissipation change when a solid body rotation is applied. It is shown that when the helicity increases, the dissipation is reduced. It is asserted that to have a drag reduction the external action should be such as to disrupt the symmetry of right- and left-handed helical structures. In this study the Navier-Stokes equations in rotational form permit the turbulent energy production to be split into a part related to the energy cascade from large to small scales and into a part related to the convection by large scales. The full simulation data have shown the latter is greater than the former in the wall region and that, on the contrary, these two terms balance each other in the central region. From the pdf of the former, it has been shown how the vortical structures are changed in the wall region by the background radiation and how they are related to the changes in the energy production.

Orlandi, P.↗

The analysis and modeling of dilatational terms in compressible turbulence

It is shown that the dilatational terms that need to be modeled in compressible turbulence include not only the pressure-dilatation term but also another term - the compressible dissipation. The nature of these dilatational terms in homogeneous turbulence is explored by asymptotic analysis of the compressible Navier-Stokes equations. A non-dimensional parameter which characterizes some compressible effects in moderate Mach number, homogeneous turbulence is identified. Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, and their results are found to be in agreement with the theoretical analysis. A model for the compressible dissipation is proposed; the model is based on the asymptotic analysis and the direct numerical simulations. This model is calibrated with reference to the DNS results regarding the influence of compressibility on the decay rate of isotropic turbulence. An application of the proposed model to the compressible mixing layer has shown that the model is able to predict the dramatically reduced growth rate of the compressible mixing layer.

Sarkar, S.↗

The analysis and modelling of dilatational terms in compressible turbulence

It is shown that the dilatational terms that need to be modeled in compressible turbulence include not only the pressure-dilatation term but also another term - the compressible dissipation. The nature of these dilatational terms in homogeneous turbulence is explored by asymptotic analysis of the compressible Navier-Stokes equations. A non-dimensional parameter which characterizes some compressible effects in moderate Mach number, homogeneous turbulence is identified. Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, and their results are found to be in agreement with the theoretical analysis. A model for the compressible dissipation is proposed; the model is based on the asymptotic analysis and the direct numerical simulations. This model is calibrated with reference to the DNS results regarding the influence of compressibility on the decay rate of isotropic turbulence. An application of the proposed model to the compressible mixing layer has shown that the model is able to predict the dramatically reduced growth rate of the compressible mixing layer.

Sarkar, S.↗

An implicit solution of the three-dimensional Navier-Stokes equations for an airfoil spanning a wind tunnel

An implicit finite-difference algorithm is developed for the numerical solution of the incompressible three dimensional Navier-Stokes equations in the non-conservative primitive-variable formulation. The flow field about an airfoil spanning a wind-tunnel is computed. The coordinate system is generated by an extension of the two dimensional body-fitted coordinate generation techniques of Thompson, as well as that of Sorenson, into three dimensions. Two dimensional grids are stacked along a spanwise coordinate defined by a simple analytical function. A Poisson pressure equation for advancing the pressure in time is arrived at by performing a divergence operation on the momentum equations. The pressure at each time-step is calculated on the assumption that continuity be unconditionally satisfied. An eddy viscosity coefficient, computed according to the algebraic turbulence formulation of Baldwin and Lomax, simulates the effects of turbulence.

Moitra, A.↗

Large Eddy Simulation of Aircraft Wake Vortices: Atmospheric Turbulence Effects

Crow instability can develop in most atmospheric turbulence levels, however, the ring vortices may not form in extremely strong turbulence cases due to strong dissipation of the vortices. It appears that strong turbulence tends to accelerate the occurrences of Crow instability. The wavelength of the most unstable mode is estimated to be about 5b(sub 0), which is less than the theoretical value of 8.6b(sub 0) (Crow, 1970) and may be due to limited domain size and highly nonlinear turbulent flow characteristics. Three-dimensional turbulence can decay wake vortices more rapidly. Axial velocity may be developed by vertical distortion of a vortex pair due to Crow instability or large turbulent eddy motion. More experiments with various non-dimensional turbulence levels are necessary to get useful statistics of wake vortex behavior due to turbulence. Need to investigate larger turbulence length scale effects by enlarging domain size or using grid nesting.

Han, Jongil↗

DNS of Flows over Periodic Hills using a Discontinuous-Galerkin Spectral-Element Method

Direct numerical simulation (DNS) of turbulent compressible flows is performed using a higher-order space-time discontinuous-Galerkin finite-element method. The numerical scheme is validated by performing DNS of the evolution of the Taylor-Green vortex and turbulent flow in a channel. The higher-order method is shown to provide increased accuracy relative to low-order methods at a given number of degrees of freedom. The turbulent flow over a periodic array of hills in a channel is simulated at Reynolds number 10,595 using an 8th-order scheme in space and a 4th-order scheme in time. These results are validated against previous large eddy simulation (LES) results. A preliminary analysis provides insight into how these detailed simulations can be used to improve Reynoldsaveraged Navier-Stokes (RANS) modeling

DNS↗

Turbulence Effects on Shaped Booms: Finite Impulse Response Filter Development

Numerical simulations of propagation through turbulent atmospheres can quantify effects on ground waveforms, but such simulations are computationally expensive. To enable quick turnaround analyses as required by NASA’s Quesst Mission, updating the N-wave filtering approach developed by researchers at The Pennsylvania State University to include shaped booms is proposed as an alternative method for estimating turbulence effects on acoustic metrics more quickly. Beginning with a nearfield pressure cylinder modeled after the on-design X-59 configuration, a database of propagation results at 45 turbulence conditions was compiled using nonlinear turbulence propagation modeling code (KZKFourier) and used as input to a process for generating finite impulse response (FIR) filters. Ground waveforms distorted by turbulence were selected to represent mean and mean±standard deviation levels for six metrics, and corresponding FIR filters were generated through a matrix deconvolution process. In order to evaluate how well the FIR filters perform, additional KZKFourier verification cases were devised with different input conditions, and results used as a benchmark. Convolution of shaped boom waveforms modeled using nonturbulent propagation simulations with the new FIR filters showed better agreement on average with KZKFourier statistical results than the N-wave-based FIR filters.

turbulence↗

On the measurement of lateral velocity derivatives in turbulent flows

Direct numerical simulation data for the lateral velocity derivative delta(u)/delta(y) at the centerline of a fully developed turbulent channel flow provide reasonable support for Wyngaard's analysis of the error involved in measuring this quantity using parallel hot wires. Numerical data in the wall region of the channel flow also provide a useful indication of how to select the separation between the wires. Justification for this choice is obtained by comparing several measured statistics of delta(u)/delta(y) with the corresponding numerical data.

Antonia, R. A.↗

Modeling the pressure-dilatation correlation

It is generally accepted that pressure dilatation, which is an additional compressibility term in turbulence transport equations, may be important for high speed flows. Recent direct simulations of homogeneous shear turbulence have given concrete evidence that the pressure dilatation is important insofar that it contributes to the reduced growth of turbulent kinetic energy due to compressibility effects. The problem of modeling pressure dilatation is addressed. A component of the pressure dilatation is isolated which exhibits temporal oscillations and, using direct numerical simulations of homogeneous shear turbulence and isotropic turbulence, show that it has a negligible contribution to the evolution of turbulent kinetic energy. Then, an analysis for the case of homogeneous turbulence is performed to obtain a model for the nonoscillatory pressure dilatation. This model algebraically relates the pressure dilatation to quantities traditionally obtained in incompressible turbulence closures. The model is validated by direct comparison with the pressure dilatation data obtained from the simulations.

Sarkar, S.↗