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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Effect of Swirl on Turbulent Structures in Supersonic Jets

Direct numerical simulation (DNS) is used to study the mechanism of generation and evolution of turbulence structures in a temporally evolving supersonic swirling round jet and also to examine the resulting acoustic radiations. Fourier spectral expansions are used in the streamwise and azimuthal directions and a 1-D b-spline Galerkin representation is used in the radial direction. Spectral-like accuracy is achieved using this numerical scheme. Direct numerical simulations, using the b-spline spectral method, are carried out starting from mean flow initial conditions which are perturbed by the most unstable linear stability eigenfunctions. It is observed that the initial.helical instability waves evolve into helical vortices which eventually breakdown into smaller scales of turbulence. 'Rib' structures similar to those seen in incompressible mixing layer flow of Rogers and Moser are observed. The jet core breakdown stage exhibits increased acoustic radiations.

Rao, Ram Mohan↗
NASA NTRS1998

Transition and Breakdown to Turbulence in Incompressible Boundary Layers

We have developed a code where the nonlinear terms are treated implicitly. The equations are discretized using the two-point fourth order compact scheme in the y-direction and the backward Euler method in the x-direction. We investigated the transition process in a Blasius boundary layer due to fundamental type breakdown. With 8 modes in the w and 3 planes, we could compute the evolution of disturbances up to Re(x)=910, which is well into the strongly nonlinear region. The transition onset point is located around Re(x)=850. The comparison with the measurements and with the DNS computations are very good up to Re(x)=880.

Balakumar, Ponnampalam↗
NASA NTRS1999

Shared Memory Parallelization of an Implicit ADI-type CFD Code

A parallelization study designed for ADI-type algorithms is presented using the OpenMP specification for shared-memory multiprocessor programming. Details of optimizations specifically addressed to cache-based computer architectures are described and performance measurements for the single and multiprocessor implementation are summarized. The paper demonstrates that optimization of memory access on a cache-based computer architecture controls the performance of the computational algorithm. A hybrid MPI/OpenMP approach is proposed for clusters of shared memory machines to further enhance the parallel performance. The method is applied to develop a new LES/DNS code, named LESTool. A preliminary DNS calculation of a fully developed channel flow at a Reynolds number of 180, Re(sub tau) = 180, has shown good agreement with existing data.

Hauser, Th.↗
NASA NTRS1998

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the Summer Program - No. 7

The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.

Source record↗
NASA NTRS1997

Annual Research Briefs

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Spinks, Debra↗
NASA NTRS1999

Low-Storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations

The derivation of storage explicit Runge-Kutta (ERK) schemes has been performed in the context of integrating the compressible Navier-Stokes equations via direct numerical simulation. Optimization of ERK methods is done across the broad range of properties, such as stability and accuracy efficiency, linear and nonlinear stability, error control reliability, step change stability, and dissipation/dispersion accuracy, subject to varying degrees of memory economization. Following van der Houwen and Wray, 16 ERK pairs are presented using from two to five registers of memory per equation, per grid point and having accuracies from third- to fifth-order. Methods have been assessed using the differential equation testing code DETEST, and with the 1D wave equation. Two of the methods have been applied to the DNS of a compressible jet as well as methane-air and hydrogen-air flames. Derived 3(2) and 4(3) pairs are competitive with existing full-storage methods. Although a substantial efficiency penalty accompanies use of two- and three-register, fifth-order methods, the best contemporary full-storage methods can be pearl), matched while still saving two to three registers of memory.

Kennedy, Chistopher A.↗
NASA NTRS1998

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗
NASA NTRS1999

Large Eddy Simulation in the Computation of Jet Noise

Noise can be predicted by solving Full (time-dependent) Compressible Navier-Stokes Equation (FCNSE) with computational domain. The fluctuating near field of the jet produces propagating pressure waves that produce far-field sound. The fluctuating flow field as a function of time is needed in order to calculate sound from first principles. Noise can be predicted by solving the full, time-dependent, compressible Navier-Stokes equations with the computational domain extended to far field - but this is not feasible as indicated above. At high Reynolds number of technological interest turbulence has large range of scales. Direct numerical simulations (DNS) can not capture the small scales of turbulence. The large scales are more efficient than the small scales in radiating sound. The emphasize is thus on calculating sound radiated by large scales.

Mankbadi, R. R.↗
NASA NTRS1999

Large Eddy Simulations and Turbulence Modeling for Film Cooling

The objective of the research is to perform Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) for film cooling process, and to evaluate and improve advanced forms of the two equation turbulence models for turbine blade surface flow analysis. The DNS/LES were used to resolve the large eddies within the flow field near the coolant jet location. The work involved code development and applications of the codes developed to the film cooling problems. Five different codes were developed and utilized to perform this research. This report presented a summary of the development of the codes and their applications to analyze the turbulence properties at locations near coolant injection holes.

Acharya, Sumanta↗
NASA NTRS1999

Comparison of AGE and Spectral Methods for the Simulation of Far-Wakes

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

Bisset, D. K.↗
NASA NTRS1999

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

The objective of this numerical/experimental study is to improve the understanding of the effects of gravity on the two-way interaction between dispersed particles (bubbles or liquid droplets) and the carrier turbulent flow. The first phase of the project considers isotropic turbulence. Turbulent homogeneous shear flows laden with droplets/bubbles will be studied in the next phase. The experiments reported here are concerned with the dispersion of liquid droplets by homogeneous turbulence under various gravitational conditions and the effect of these droplets on the evolution of the turbulence of the carrier fluid (air). Direct numerical simulations (DNS) of bubble - laden isotropic decaying turbulence are performed using the two-fluid approach (TF) instead of the Eulerian-Lagrangian approach (EL). The motivation for using the TF formulation is that EL requires considerable computational resources especially for the case of two-way coupling where the instantaneous trajectories of a large number of individual bubbles need to be computed. The TF formulation is developed by spatially averaging the instantaneous equations of the carrier flow and bubble phase over a scale of the order of the Kolmogorov length scale which, in our case, is much larger than the bubble diameter. On that scale, the bubbles are treated as a continuum (without molecular diffusivity) characterized by the bubble phase velocity field and concentration (volume fraction). The bubble concentration, C, is assumed small enough to neglect the bubble-bubble interactions.

Elghobashi, Said↗
NASA NTRS2000

Numerical Simulations of Homogeneous Turbulence Using Lagrangian-Averaged Navier-Stokes Equations

The Lagrangian-averaged Navier-Stokes (LANS) equations are numerically evaluated as a turbulence closure. They are derived from a novel Lagrangian averaging procedure on the space of all volume-preserving maps and can be viewed as a numerical algorithm which removes the energy content from the small scales (smaller than some a priori fixed spatial scale alpha) using a dispersive rather than dissipative mechanism, thus maintaining the crucial features of the large scale flow. We examine the modeling capabilities of the LANS equations for decaying homogeneous turbulence, ascertain their ability to track the energy spectrum of fully resolved direct numerical simulations (DNS), compare the relative energy decay rates, and compare LANS with well-accepted large eddy simulation (LES) models.

Mohseni, Kamran↗
NASA NTRS2000

CVS Decomposition of 3D Homogeneous Turbulence Using Orthogonal Wavelets

This paper compares the filtering used in Coherent Vortex Simulation (CVS) decomposition with an orthogonal wavelet basis, with the Proper Orthogonal Decomposition (POD) or Fourier filtering. Both methods are applied to a field of Direct Numerical Simulation (DNS) data of 3D forced homogeneous isotropic turbulence at microscale Reynolds number R(sub lambda) = 168. We show that, with only 3%N retained modes, CVS filtering separates the coherent vortex tubes from the incoherent background flow. The latter is structureless, has an equipartition energy spectrum, and has a Gaussian velocity probability distribution function (PDF) and an exponential vorticity PDF. On the other hand, the Fourier basis does not extract the coherent vortex tubes cleanly and leaves organized structures in the residual high wavenumber modes whose PDFs are stretched exponentials for both the velocity and the vorticity.

Farge, Marie↗
NASA NTRS2000

CVS Filtering of 3D Turbulent Mixing Layers Using Orthogonal Wavelets

Coherent Vortex Simulation (CVS) filtering has been applied to Direct Numerical Simulation (DNS) data of forced and unforced time-developing turbulent mixing layers. CVS filtering splits the turbulent flow into two orthogonal parts, one corresponding to coherent vortices and the other to incoherent background flow. We have shown that the coherent vortices can be represented by few wavelet modes and that these modes are sufficient to reproduce the vorticity probability distribution function (PDF) and the energy spectrum over the entire inertial range. The remaining incoherent background flow is homogeneous, has small amplitude, and is uncorrelated. These results are compared with those obtained for the same compression rate using large eddy simulation (LES) filtering. In contrast to the incoherent background flow of CVS filtering, the LES subgrid scales have a much larger amplitude and are correlated, which makes their statistical modeling more difficult.

Schneider, Kai↗
NASA NTRS2000

A New Method for Accurate Treatment of Flow Equations in Cylindrical Coordinates Using Series Expansions

The motivation of this work is the ongoing effort at the Center for Turbulence Research (CTR) to use large eddy simulation (LES) techniques to calculate the noise radiated by jet engines. The focus on engine exhaust noise reduction is motivated by the fact that a significant reduction has been achieved over the last decade on the other main sources of acoustic emissions of jet engines, such as the fan and turbomachinery noise, which gives increased priority to jet noise. To be able to propose methods to reduce the jet noise based on results of numerical simulations, one first has to be able to accurately predict the spatio-temporal distribution of the noise sources in the jet. Though a great deal of understanding of the fundamental turbulence mechanisms in high-speed jets was obtained from direct numerical simulations (DNS) at low Reynolds numbers, LES seems to be the only realistic available tool to obtain the necessary near-field information that is required to estimate the acoustic radiation of the turbulent compressible engine exhaust jets. The quality of jet-noise predictions is determined by the accuracy of the numerical method that has to capture the wide range of pressure fluctuations associated with the turbulence in the jet and with the resulting radiated noise, and by the boundary condition treatment and the quality of the mesh. Higher Reynolds numbers and coarser grids put in turn a higher burden on the robustness and accuracy of the numerical method used in this kind of jet LES simulations. As these calculations are often done in cylindrical coordinates, one of the most important requirements for the numerical method is to provide a flow solution that is not contaminated by numerical artifacts. The coordinate singularity is known to be a source of such artifacts. In the present work we use 6th order Pade schemes in the non-periodic directions to discretize the full compressible flow equations. It turns out that the quality of jet-noise predictions using these schemes is especially sensitive to the type of equation treatment at the singularity axis. The objective of this work is to develop a generally applicable numerical method for treating the singularities present at the polar axis, which is particularly suitable for highly accurate finite-differences schemes (e.g., Pade schemes) on non-staggered grids. The main idea is to reinterpret the regularity conditions developed in the context of pseudo-spectral methods. A set of exact equations at the singularity axis is derived using the appropriate series expansions for the variables in the original set of equations. The present treatment of the equations preserves the same level of accuracy as for the interior scheme. We also want to point out the wider utility of the method, proposed here in the context of compressible flow equations, as its extension for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical coordinates with finite-differences schemes of various level of accuracy is straightforward. The robustness and accuracy of the proposed technique is assessed by comparing results from simulations of laminar forced-jets and turbulent compressible jets using LES with similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r = 0.

Constantinescu, G.S.↗
NASA NTRS2000

Perspectives on the Future of CFD

This viewgraph presentation gives an overview of the future of computational fluid dynamics (CFD), which in the past has pioneered the field of flow simulation. Over time CFD has progressed as computing power. Numerical methods have been advanced as CPU and memory capacity increases. Complex configurations are routinely computed now and direct numerical simulations (DNS) and large eddy simulations (LES) are used to study turbulence. As the computing resources changed to parallel and distributed platforms, computer science aspects such as scalability (algorithmic and implementation) and portability and transparent codings have advanced. Examples of potential future (or current) challenges include risk assessment, limitations of the heuristic model, and the development of CFD and information technology (IT) tools.

Kwak, Dochan↗
NASA NTRS2000

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗
NASA NTRS1997

Direct Numerical Simulation of Supersonic Turbulent Boundary Layers

The talk will present some initial results from the direct numerical simulation (DNS) of compressible turbulent boundary layers. We solve numerically the compressible Navier-Stokes equations using a method based on Spalart's transformation for the incompressible turbulent boundary layer. This allows the spatially developing boundary layer to be transformed to a calculation with periodic boundary conditions in the streamwise and spanwise directions. The equations are solved using Fourier expansions in the horizontal directions and B-splines in the wall-normal direction. The first simulation is at Mach 2.5 with a momentum thickness Reynolds number based on wall viscosity of R(sub theta(sup 1)) = 825. We are examining the physics of the compressible boundary layer using turbulence statistics and budget equations. The turbulence statistics include: rms (root mean square) and mean profiles, energy spectra, and two-point correlations. It is found that there are large density gradients which require significantly more resolution than the incompressible case.

Guarini, Stephen↗