SEARCH · Search NASA
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.
Large Eddy simulation (LES) models are presented and evaluated on a database obtained from Direct Numerical Simulation (DNS) of a three-dimensional temporal mixing layer with evaporating drops.
This technology development originated from the need to assess the debris threat resulting from soil material erosion induced by landing spacecraft rocket plume impingement on extraterrestrial planetary surfaces. The impact of soil debris was observed to be highly detrimental during NASA s Apollo lunar missions and will pose a threat for any future landings on the Moon, Mars, and other exploration targets. The innovation developed under this program provides a simulation tool that combines modeling of the diverse disciplines of rocket plume impingement gas dynamics, granular soil material liberation, and soil debris particle kinetics into one unified simulation system. The Unified Flow Solver (UFS) developed by CFDRC enabled the efficient, seamless simulation of mixed continuum and rarefied rocket plume flow utilizing a novel direct numerical simulation technique of the Boltzmann gas dynamics equation. The characteristics of the soil granular material response and modeling of the erosion and liberation processes were enabled through novel first principle-based granular mechanics models developed by the University of Florida specifically for the highly irregularly shaped and cohesive lunar regolith material. These tools were integrated into a unique simulation system that accounts for all relevant physics aspects: (1) Modeling of spacecraft rocket plume impingement flow under lunar vacuum environment resulting in a mixed continuum and rarefied flow; (2) Modeling of lunar soil characteristics to capture soil-specific effects of particle size and shape composition, soil layer cohesion and granular flow physics; and (3) Accurate tracking of soil-borne debris particles beginning with aerodynamically driven motion inside the plume to purely ballistic motion in lunar far field conditions.
Motivated by large eddy simulation (LES) modeling of supercritical turbulent flows, transitional states of databases obtained from direct numerical simulations (DNS) of binary-species supercritical temporal mixing layers were examined to understand the subgrid-scale dissipation, and its variation with filter size. Examination of the DSN-scale domain- averaged dissipation confirms previous findings that, out of the three modes of viscous, temperature and species-mass dissipation, the species-mass dissipation is the main contributor to the total dissipation. The results revealed that the percentage of species-mass by total dissipation is nearly invariant across species systems and initial conditions. This dominance of the species-mass dissipation is due to high-density-gradient magnitude (HDGM) regions populating the flow under the supercritical conditions of the simulations; such regions have also been observed in fully turbulent supercritical flows. The domain average being the result of both the local values and the extent of the HDGM regions, the expectations were that the response to filtering would vary with these flow characteristics. All filtering here is performed in the dissipation range of the Kolmogorov spectrum, at filter sizes from 4 to 16 times the DNS grid spacing. The small-scale (subgrid scale, SGS) dissipation was found by subtracting the filtered-field dissipation from the DNS-field dissipation. In contrast to the DNS dissipation, the SGS dissipation is not necessarily positive; negative values indicate backscatter. Backscatter was shown to be spatially widespread in all modes of dissipation and in the total dissipation (25 to 60 percent of the domain). The maximum magnitude of the negative subgrid- scale dissipation was as much as 17 percent of the maximum positive subgrid- scale dissipation, indicating that, not only is backscatter spatially widespread in these flows, but it is considerable in magnitude and cannot be ignored for the purposes of LES modeling. The Smagorinsky model, for example, is unsuited for modeling SGS fluxes in the LES because it cannot render backscatter. With increased filter size, there is only a modest decrease in the spatial extent of backscatter. The implication is that even at large LES grid spacing, the issue of backscatter and related SGS-flux modeling decisions are unavoidable. As a fraction of the total dissipation, the small-scale dissipation is between 10 and 30 percent of the total dissipation for a filter size that is four times the DNS grid spacing, with all OH cases bunched at 10 percent, and the HN cases spanning 24 30 percent. A scale similarity was found in that the domain-average proportion of each small-scale dissipation mode, with respect to the total small-scale dissipation, is very similar to equivalent results at the DNS scale. With increasing filter size, the proportion of the small-scale dissipation in the dissipation increases substantially, although not quite proportionally. When the filter size increases by four-fold, 52 percent for all OH runs, and 70 percent for HN runs, of the dissipation is contained in the subgrid-scale portion with virtually no dependence on the initial conditions of the DNS. The indications from the dissipation analysis are that modeling efforts in LES of thermodynamically supercritical flows should be focused primarily on mass-flux effects, with temperature and viscous effects being secondary. The analysis also reveals a physical justification for scale-similarity type models, although the suitability of these will need to be confirmed in a posteriori studies.
Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The favorable pressure gradient is induced by the streamwise curvature of the two-dimensional, planar, convex measurement surface used during experiments at the Texas A&M University. The DNS data shows good comparison with the measured velocity profiles, strain rates, and some, but not all, of the Reynolds-stress components. The discrepancies between the predicted and the measured wall-normal as well as shear stress components are primarily attributed to the lower than actual values inferred from typical PIV measurements of turbulent boundary layers. The DNS data shows a zero or slightly negative Reynolds shear stress in the outer part of the boundary-layer, which is indicative of the decaying turbulent motion under a strong favorable pressure gradient. The DNS data is also compared with the results of RANS computations based on commonly used zero, one, and two equation eddy-viscosity models. The RANS models yield reasonable comparisons with the DNS-based skin friction under zero and weak pressure gradients, but significant discrepancies under a strong pressure gradient. The k-w SST model provided the best overall predictions of skin friction, except in the region where the flow transitions from a favorable to an adverse pressure gradient. While the RANS models examined herein also give good predictions of the Reynolds shear stress under a sufficiently weak pressure gradient, none of those models are able to appropriately capture the reduction in the Reynolds stresses when the flow was subjected to a strong pressure gradient. An a priori assessment of the turbulent heat-flux prediction based on the assumption of a constant turbulent Prandtl number with the DNS data shows that while the constant turbulent Prandtl number model is effective in predicting the wall-normal component of turbulent heat flux, it does not capture the turbulent heat transfer in the streamwise direction for all the pressure gradient cases. The failure of the constant turbulent Prandtl number model highlights a requirement for more advanced models of the turbulent heat flux.
The current status of numerical simulation techniques for the transition to turbulence in incompressible channel and boundary-layer flows is surveyed, and typical results are presented graphically. The focus is on direct numerical simulations based on the full nonlinear time-dependent Navier-Stokes equations without empirical closure assumptions for prescribed initial and boundary conditions. Topics addressed include the vibrating ribbon problem, space and time discretization, initial and boundary conditions, alternative methods based on the triple-deck approximation, two-dimensional channel and boundary-layer flows, three-dimensional boundary layers, wave packets and turbulent spots, compressible flows, transition control, and transition modeling.
Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models.The wall geometry and flow conditions of the DNS are representative of the experimental data for a Mach 4.9 turbulent boundary layer that was tested on a two-dimensional planar convexwall model in the high-speed blow-down wind tunnel located at the National Aerothermochemistry Laboratory at Texas A&M University. The favorable pressure gradient is induced bythe streamwise surface curvature. The DNS results will be compared with experimental measurements and RANS computations for the same flow conditions and wall geometry. The DNS datasets will also be analyzed to provide an improved physical understanding of hypersonic boundary-layer turbulence subject to mechanical nonequilibrium, including the validity of Morkovin’s hypothesis and the various available compressibility transformations, as well as to obtain turbulence statistics that are relevant to RANS modeling, including the distribution ofeddy viscosity, the budgets of turbulent kinetic energy and Reynolds stresses, and the velocity-temperature correlations. The RANS predictions will be obtained by using both the commonly used one- and two-equation models that assume a constant turbulent Prandtl number and an algebraic, variable Prandtl number model for the turbulent energy flux.
Large-eddy simulation (LES) is conducted of a three-dimensional temporal mixing layer whose lower stream is initially laden with liquid drops which may evaporate during the simulation. The gas-phase equations are written in an Eulerian frame for two perfect gas species (carrier gas and vapour emanating from the drops), while the liquid-phase equations are written in a Lagrangian frame. The effect of drop evaporation on the gas phase is considered through mass, species, momentum and energy source terms. The drop evolution is modelled using physical drops, or using computational drops to represent the physical drops. Simulations are performed using various LES models previously assessed on a database obtained from direct numerical simulations (DNS). These LES models are for: (i) the subgrid-scale (SGS) fluxes and (ii) the filtered source terms (FSTs) based on computational drops. The LES, which are compared to filtered-and-coarsened (FC) DNS results at the coarser LES grid, are conducted with 64 times fewer grid points than the DNS, and up to 64 times fewer computational than physical drops. It is found that both constant-coefficient and dynamic Smagorinsky SGS-flux models, though numerically stable, are overly dissipative and damp generated small-resolved-scale (SRS) turbulent structures. Although the global growth and mixing predictions of LES using Smagorinsky models are in good agreement with the FC-DNS, the spatial distributions of the drops differ significantly. In contrast, the constant-coefficient scale-similarity model and the dynamic gradient model perform well in predicting most flow features, with the latter model having the advantage of not requiring a priori calibration of the model coefficient. The ability of the dynamic models to determine the model coefficient during LES is found to be essential since the constant-coefficient gradient model, although more accurate than the Smagorinsky model, is not consistently numerically stable despite using DNS-calibrated coefficients. With accurate SGS-flux models, namely scale-similarity and dynamic gradient, the FST model allows up to a 32-fold reduction in computational drops compared to the number of physical drops, without degradation of accuracy; a 64-fold reduction leads to a slight decrease in accuracy.
Direct numerical simulations (DNS) of adverse-pressure-gradient turbulent boundary layers over a planar concave wall are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The wall geometry and flow conditions of the DNS are representative of the experimental data for a Mach 4.9 turbulent boundary layer that was tested on a two-dimensional planar concave wall model in the high-speed blow-down wind tunnel located at the National Aerothermochemistry Laboratory at Texas A&M University (TAMU). The DNS was validated against the experimental results of TAMU for the same flow conditions and wall geometry. An analysis of the DNS datasets was also conducted to provide an assessment of the validity of Morkovin’s hypothesis and the strong Reynolds analog for turbulence subject to mechanical nonequilibrium. In addition to the DNS results, RANS predictions are obtained by using the Baldwin-Lomax (BL), Spalart-Allmaras (SA), and the k - w SST turbulence models. The comparisons between RANS and DNS showed little impact of an adverse pressure gradient on the accuracy of these models, at least up to an incompressible Clauser pressure gradient parameter of beta(sub inc) 1.22. While the Boussinesq assumption provided reasonable predictions for the Reynolds shear stress, it failed to adequately predict the normal components of the Reynolds stress.
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.
This second of three reports on a computational study of a mixing layer laden with evaporating liquid drops presents the evaluation of Large Eddy Simulation (LES) models. The LES models were evaluated on an existing database that had been generated using Direct Numerical Simulation (DNS). The DNS method and the database are described in the first report of this series, Part 1 of a Computational Study of a Drop-Laden Mixing Layer (NPO-30719), NASA Tech Briefs, Vol. 28, No.7 (July 2004), page 59. The LES equations, which are derived by applying a spatial filter to the DNS set, govern the evolution of the larger scales of the flow and can therefore be solved on a coarser grid. Consistent with the reduction in grid points, the DNS drops would be represented by fewer drops, called computational drops in the LES context. The LES equations contain terms that cannot be directly computed on the coarser grid and that must instead be modeled. Two types of models are necessary: (1) those for the filtered source terms representing the effects of drops on the filtered flow field and (2) those for the sub-grid scale (SGS) fluxes arising from filtering the convective terms in the DNS equations. All of the filtered-sourceterm models that were developed were found to overestimate the filtered source terms. For modeling the SGS fluxes, constant-coefficient Smagorinsky, gradient, and scale-similarity models were assessed and calibrated on the DNS database. The Smagorinsky model correlated poorly with the SGS fluxes, whereas the gradient and scale-similarity models were well correlated with the SGS quantities that they represented.
On the basis of direct numerical simulations of rapidly compressed turbulence, Zeman and Coleman have developed a model to represent rapid directional compression contribution to the pressure dilatation term in the turbulent kinetic energy equation. The model has been implemented in the CFD code for simulation of supersonic compression corner flow with an extended separated region. The computational results have shown a significant improvement with respect to the baseline solution given by the standard k- epsilon turbulence model which does not contain any compressibility corrections.
Most flows of aerodynamic interest are compressible and turbulent. However, our present knowledge on the structures and mechanisms of turbulence is mostly based on incompressible flows. In the present work, compressibility effects in turbulent, high-speed, boundary layer flows are systematically investigated using the Direct Numerical Simulation (DNS) approach. Three-dimensional, time-dependent, fully nonlinear, compressible Navier-Stokes equations were numerically integrated by high-order finite-difference methods; no modeling for turbulence is used during the solution because the available resolution is sufficient to capture the relevant scales. The boundary layer problem deals with fully-turbulent compressible flows over flat geometries. Apart from its practical relevance to technological flows, turbulent compressible boundary layer flow is the simplest experimentally realizable turbulent compressible flow. Still, measuring difficulties prohibit a detailed experimental description of the flow, especially in the near-wall region. DNS studies provide a viable means to probe the physics of compressible turbulence in this region. The focus of this work is to explore the paths of energy transfer through which compressible turbulence is sustained. The structural similarities and differences between the incompressible and compressible turbulence are also investigated. The energy flow patterns or energy cascades are found to be directly related to the evolution of vortical structures which are generated in the near-wall region. Near-wall structures, and mechanisms which are not readily accessible through physical experiments are analyzed and their critical role on the evolution and the behavior of the flow is documented extensively.
Pressure fluctuations in a turbulent channel flow are investigated by analyzing a database obtained from a direct numerical simulation. Detailed statistics associated with the pressure fluctuations are presented. Characteristics associated with the rapid (linear) and slow (nonlinear) pressure are discussed. It is found that the slow pressure fluctuations are larger than the rapid pressure fluctuations throughout the channel except very near the wall, where they are about the same magnitude. This is contrary to the common belief that the nonlinear source terms are negligible compared to the linear source terms. Probability density distributions, power spectra, and two-point correlations are examined to reveal the characteristics of the pressure fluctuations. The global dependence of the pressure fluctuations and pressure-strain correlations are also examined by evaluating the integral associated with Green's function representations of them. In the wall region where the pressure-strain terms are large, most contributions to the pressure-strain terms are from the wall region (i.e., local), whereas away from the wall where the pressure-strain terms are small, contributions are global. Structures of instantaneous pressure and pressure gradients at the wall and the corresponding vorticity field are examined.
Pressure fluctuations in a turbulent channel flow are investigated by analyzing a database obtained from a direct numerical simulation. Detailed statistics associated with the pressure fluctuations are presented. Characteristics associated with the rapid (linear) and slow (nonlinear) pressure are discussed. It is found that the slow pressure fluctuations are larger than the rapid pressure fluctuations throughout the channel except very near the wall, where they are about the same magnitude. This is contrary to the common belief that the nonlinear source terms are negligible compared to the linear source terms. Probability density distributions, power spectra, and two-point correlations are examined to reveal the characteristics of the pressure fluctuations. The global dependence of the pressure fluctuations and pressure-strain correlations are also examined by evaluating the integral associated with Green's function representations of them. In the wall region where the pressure-strain terms are large, most contributions to the pressure-strain terms are from the wall region (i.e., local), whereas away from the wall where the pressure-strain terms are small, contributions are global. Structures of instantaneous pressure and pressure gradients at the wall and the corresponding vorticity field are examined.
Despite decades of development of unstructured mesh methods, high-fidelity time-accurate simulations are still predominantly carried out on structured, or unstructured hexahedral meshes by using high-order finite-difference, weighted essentially non-oscillatory (WENO), or hybrid schemes formed by their combinations. In this work, the space-time conservation element solution element (CESE) method is used to simulate several flow problems including supersonic jet/shock interaction and its impact on launch vehicle acoustics, and direct numerical simulations of turbulent flows using tetrahedral meshes. This paper provides a status report for the continuing development of the space-time conservation element solution element (CESE) numerical and software framework under the Revolutionary Computational Aerosciences (RCA) project. Solution accuracy and large-scale parallel performance of the numerical framework is assessed with the goal of providing a viable paradigm for future high-fidelity flow physics simulations.
The use of simulation data bases for the examination of turbulent flows is an effective research tool. Studies of the structure of turbulence have been hampered by the limited number of probes and the impossibility of measuring all desired quantities. Also, flow visualization is confined to the observation of passive markers with limited field of view and contamination caused by time-history effects. Computer flow fields are a new resource for turbulence research, providing all the instantaneous flow variables in three-dimensional space. Simulation data bases also provide much-needed information for phenomenological turbulence modeling. Three dimensional velocity and pressure fields from direct simulations can be used to compute all the terms in the transport equations for the Reynolds stresses and the dissipation rate. However, only a few, geometrically simple flows have been computed by direct numerical simulation, and the inventory of simulation does not fully address the current modeling needs in complex turbulent flows. The availability of three-dimensional flow fields also poses challenges in developing new techniques for their analysis, techniques based on experimental methods, some of which are used here for the analysis of direct-simulation data bases in studies of the mechanics of turbulent flows.
New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.
New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 (exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.