Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DNS”

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 19 records

DNS and Embedded DNS as Tools for Investigating Unsteady Heat Transfer Phenomena in Turbines

DNS is a powerful tool with high potential for investigating unsteady heat transfer and fluid flow phenomena, in particular for cases involving transition to turbulence and/or large coherent structures. - DNS of idealized configurations related to turbomachinery components is already possible. - For more realistic configurations and the inclusion of more effects, reduction of computational cost is key issue (e.g., hybrid methods). - Approach pursued here: Embedded DNS ( segregated coupling of DNS with LES and/or RANS). - Embedded DNS is an enabling technology for many studies. - Pre-transitional heat transfer and trailing-edge cutback film-cooling are good candidates for (embedded) DNS studies.

vonTerzi, Dominic↗

LES versus DNS: A comparative study

We have performed Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) of forced isotropic turbulence at moderate Reynolds numbers. The subgrid scale model used in the LES is based on an eddy viscosity which adjusts instantaneously the energy spectrum of the LES to that of the DNS. The statistics of the large scales of the DNS (filtered DNS field or fDNS) are compared to that of the LES. We present results for the transfer spectra, the skewness and flatness factors of the velocity components, the PDF's of the angle between the vorticity and the eigenvectors of the rate of strain, and that between the vorticity and the vorticity stretching tensor. The above LES statistics are found to be in good agreement with those measured in the fDNS field. We further observe that in all the numerical measurements, the trend was for the LES field to be more gaussian than the fDNS field. Future research on this point is planned.

Shtilman, L.↗

DNS and LES of a Shear-Free Mixing Layer

The purpose of this work is twofold. First, given the computational resources available today, it is possible to reach, using DNS, higher Reynolds numbers than in Briggs et al.. In the present study, the microscale Reynolds numbers reached in the low- and high-energy homogeneous regions are, respectively, 32 and 69. The results reported earlier can thus be complemented and their robustness in the presence of increased turbulence studied. The second aim of this work is to perform a detailed and documented LES of the shear-free mixing layer. In that respect, the creation of a DNS database at higher Reynolds number is necessary in order to make meaningful LES assessments. From the point of view of LES, the shear-free mixing-layer is interesting since it allows one to test how traditional LES models perform in the presence of an inhomogeneity without having to deal with difficult numerical issues. Indeed, as argued in Briggs et al., it is possible to use a spectral code to study the shear-free mixing layer and one can thus focus on the accuracy of the modelling while avoiding contamination of the results by commutation errors etc. This paper is organized as follows. First we detail the initialization procedure used in the simulation. Since the flow is not statistically stationary, this initialization procedure has a fairly strong influence on the evolution. Although we will focus here on the shear-free mixing layer, the method proposed in the present work can easily be used for other flows with one inhomogeneous direction. The next section of the article is devoted to the description of the DNS. All the relevant parameters are listed and comparison with the Veeravalli & Warhaft experiment is performed. The section on the LES of the shear-free mixing layer follows. A detailed comparison between the filtered DNS data and the LES predictions is presented. It is shown that simple eddy viscosity models perform very well for the present test case, most probably because the flow seems to be almost isotropic in the small-scale range that is not resolved by the LES.

Knaepen, B.↗

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↗

Evaluation of a vortex-based subgrid stress model using DNS databases

The performance of a SubGrid Stress (SGS) model for Large-Eddy Simulation (LES) developed by Misra k Pullin (1996) is studied for forced and decaying isotropic turbulence on a 32(exp 3) grid. The physical viability of the model assumptions are tested using DNS databases. The results from LES of forced turbulence at Taylor Reynolds number R(sub (lambda)) approximately equals 90 are compared with filtered DNS fields. Probability density functions (pdfs) of the subgrid energy transfer, total dissipation, and the stretch of the subgrid vorticity by the resolved velocity-gradient tensor show reasonable agreement with the DNS data. The model is also tested in LES of decaying isotropic turbulence where it correctly predicts the decay rate and energy spectra measured by Comte-Bellot & Corrsin (1971).

Misra, Ashish↗

Hydroacoustic forcing function modeling using DNS database

A wall pressure frequency spectrum model (Blake 1971 ) has been evaluated using databases from Direct Numerical Simulations (DNS) of a turbulent boundary layer (Na & Moin 1996). Good agreement is found for moderate to strong adverse pressure gradient flows in the absence of separation. In the separated flow region, the model underpredicts the directly calculated spectra by an order of magnitude. The discrepancy is attributed to the violation of the model assumptions in that part of the flow domain. DNS computed coherence length scales and the normalized wall pressure cross-spectra are compared with experimental data. The DNS results are consistent with experimental observations.

Zawadzki, I.↗

A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing

Currently, the modeling of supercritical-pressure flows through Large Eddy Simulation (LES) uses models derived for atmospheric-pressure flows. Those atmospheric-pressure flows do not exhibit the particularities of high densitygradient magnitude features observed both in experiments and simulations of supercritical-pressure flows in the case of two species mixing. To assess whether the current LES modeling is appropriate and if found not appropriate to propose higher-fidelity models, a LES a posteriori study has been conducted for a mixing layer that initially contains different species in the lower and upper streams, and where the initial pressure is larger than the critical pressure of either species. An initially-imposed vorticity perturbation promotes roll-up and a double pairing of four initial span-wise vortices into an ultimate vortex that reaches a transitional state. The LES equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equation set uses transport properties depending on the thermodynamic variables. Unlike all LES models to date, the differential equations contain, additional to the subgrid scale (SGS) fluxes, a new SGS term that is a pressure correction in the momentum equation. This additional term results from filtering of Direct Numerical Simulation (DNS) equations, and represents the gradient of the difference between the filtered pressure and the pressure computed from the filtered flow field. A previous a priori analysis, using a DNS database for the same configuration, found this term to be of leading order in the momentum equation, a fact traced to the existence of high-densitygradient magnitude regions that populated the entire flow; in the study, models were proposed for the SGS fluxes as well as this new term. In the present study, the previously proposed constantcoefficient SGS-flux models of the a priori investigation are tested a posteriori in LES, devoid of or including, the SGS pressure correction term. The present pressure-correction model is different from, and more accurate as well as less computationally intensive than that of the a priori study. The constant-coefficient SGS-flux models encompass the Smagorinsky (SMC), in conjunction with the Yoshizawa (YO) model for the trace, the Gradient (GRC) and the Scale Similarity (SSC) models, all exercised with the a priori study constant coefficients calibrated at the transitional state. The LES comparison is performed with the filtered- and-coarsened (FC) DNS, which represents an ideal LES solution. Expectably, an LES model devoid of SGS terms is shown to be considerably inferior to models containing SGS effects. Among models containing SGS effects, those including the pressure-correction term are substantially superior to those devoid of it. The sensitivity of the predictions to the initial conditions and grid size are also investigated. Thus, it has been discovered that, additional to the atmospheric-pressure models currently used, a new model is necessary to simulate supercritical-pressure flows. This model depends on the thermodynamic characteristics of the chemical species involved.

Bellan, Josette↗

DNS, Enstrophy Balance, and the Dissipation Equation in a Separated Turbulent Channel Flow

The turbulent flows through a plane channel and a channel with a constriction (2-D hill) are numerically simulated using DNS and RANS calculations. The Navier-Stokes equations in the DNS are solved using a higher order kinetic energy preserving central schemes and a fifth order accurate upwind biased WENO scheme for the space discretization. RANS calculations are performed using the NASA code CFL3D with the komega SST two-equation model and a full Reynolds stress model. Using DNS, the magnitudes of different terms that appear in the enstrophy equation are evaluated. The results show that the dissipation and the diffusion terms reach large values at the wall. All the vortex stretching terms have similar magnitudes within the buffer region. Beyond that the triple correlation among the vorticity and strain rate fluctuations becomes the important kinematic term in the enstrophy equation. This term is balanced by the viscous dissipation. In the separated flow, the triple correlation term and the viscous dissipation term peak locally and balance each other near the separated shear layer region. These findings concur with the analysis of Tennekes and Lumley, confirming that the energy transfer terms associated with the small-scale dissipation and the fluctuations of the vortex stretching essentially cancel each other, leaving an equation for the dissipation that is governed by the large-scale motion.

Balakumar, Ponnampalam↗

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↗

DNS of Flow in a Low-Pressure Turbine Cascade Using a Discontinuous-Galerkin Spectral-Element Method

A new computational capability under development for accurate and efficient high-fidelity direct numerical simulation (DNS) and large eddy simulation (LES) of turbomachinery is described. This capability is based on an entropy-stable Discontinuous-Galerkin spectral-element approach that extends to arbitrarily high orders of spatial and temporal accuracy and is implemented in a computationally efficient manner on a modern high performance computer architecture. A validation study using this method to perform DNS of flow in a low-pressure turbine airfoil cascade are presented. Preliminary results indicate that the method captures the main features of the flow. Discrepancies between the predicted results and the experiments are likely due to the effects of freestream turbulence not being included in the simulation and will be addressed in the final paper.

Low-Pressure Turbine Cascade↗

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗

Compressible Turbulent Channel Flows: DNS Results and Modeling

The present paper addresses some topical issues in modeling compressible turbulent shear flows. The work is based on direct numerical simulation of two supersonic fully developed channel flows between very cold isothermal walls. Detailed decomposition and analysis of terms appearing in the momentum and energy equations are presented. The simulation results are used to provide insights into differences between conventional time-and Favre-averaging of the mean-flow and turbulent quantities. Study of the turbulence energy budget for the two cases shows that the compressibility effects due to turbulent density and pressure fluctuations are insignificant. In particular, the dilatational dissipation and the mean product of the pressure and dilatation fluctuations are very small, contrary to the results of simulations for sheared homogeneous compressible turbulence and to recent proposals for models for general compressible turbulent flows. This provides a possible explanation of why the Van Driest density-weighted transformation is so successful in correlating compressible boundary layer data. Finally, it is found that the DNS data do not support the strong Reynolds analogy. A more general representation of the analogy is analysed and shown to match the DNS data very well.

Huang, P. G.↗

The Use of DNS in Turbulence Modeling

The use of Direct numerical simulations (DNS) data in developing and testing turbulence models is reviewed. The data is used to test turbulence models at all levels: algebraic, one-equation, two-equation and full Reynolds stress models were tested. Particular examples on the development of models for the dissipation rate equation are presented. Homogeneous flows are used to test new scaling arguments for the various terms in the dissipation rate equation. The channel flow data is used to develop modifications to the equation model that take into account near-wall effects. DNS of compressible flows under mean compression are used in testing new compressible modifications to the two-equation models.

Mansour, Nagi N.↗

Detailed Comparison of DNS to PSE for Oblique Breakdown at Mach 3

A pair of oblique waves at low amplitudes is introduced in a supersonic flat-plate boundary layer. Their downstream development and the concomitant process of laminar to turbulent transition is then investigated numerically using Direct Numerical Simulations (DNS) and Parabolized Stability Equations (PSE). This abstract is the last part of an extensive study of the complete transition process initiated by oblique breakdown at Mach 3. In contrast to the previous simulations, the symmetry condition in the spanwise direction is removed for the simulation presented in this abstract. By removing the symmetry condition, we are able to confirm that the flow is indeed symmetric over the entire computational domain. Asymmetric modes grow in the streamwise direction but reach only small amplitude values at the outflow. Furthermore, this abstract discusses new time-averaged data from our previous simulation CASE 3 and compares PSE data obtained from NASA's LASTRAC code to DNS results.

Mayer, Christian S. J.↗

DNS of Laminar-Turbulent Transition in Swept-Wing Boundary Layers

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing and the mode selected for forcing corresponds to the most amplified secondary instability mode that, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. An inlet boundary condition is carefully designed to allow for accurate injection of instability wave modes and minimize acoustic reflections at numerical boundaries. Nonlinear parabolized stability equation (PSE) predictions compare well with the DNS in terms of modal amplitudes and modal shape during the strongly nonlinear phase of the secondary instability mode. During the transition process, the skin friction coefficient rises rather rapidly and the wall-shear distribution shows a sawtooth pattern that is analogous to the previously documented surface flow visualizations of transition due to stationary crossflow instability. Fully turbulent features are observed in the downstream region of the flow.

Duan, L.↗

DNS of Low-Pressure Turbine Cascade Flows with Elevated Inflow Turbulence Using a Discontinuous-Galerkin Spectral-Element Method

Recent progress towards developing a new computational capability for accurate and efficient high-fidelity direct numerical simulation (DNS) and large-eddy simulation (LES) of turbomachinery is described. This capability is based on an entropy- stable Discontinuous-Galerkin spectral-element approach that extends to arbitrarily high orders of spatial and temporal accuracy, and is implemented in a computationally efficient manner on a modern high performance computer architecture. An inflow turbulence generation procedure based on a linear forcing approach has been incorporated in this framework and DNS conducted to study the effect of inflow turbulence on the suction- side separation bubble in low-pressure turbine (LPT) cascades. The T106 series of airfoil cascades in both lightly (T106A) and highly loaded (T106C) configurations at exit isentropic Reynolds numbers of 60,000 and 80,000, respectively, are considered. The numerical simulations are performed using 8th-order accurate spatial and 4th-order accurate temporal discretization. The changes in separation bubble topology due to elevated inflow turbulence is captured by the present method and the physical mechanisms leading to the changes are explained. The present results are in good agreement with prior numerical simulations but some expected discrepancies with the experimental data for the T106C case are noted and discussed.

Garai, Anirban↗