Search NASA⌕ Search

SEARCH · Search NASA

Results for “Turbulent 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 145 records · Page 8

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

Bui, Trong T.↗

The incremental unknowns-a multilevel scheme for the simulation of turbulent channel flows

In numerical simulation of complex flows, it is important to identify different length scales of the flow and treat them differently. In this report, we introduce a new multilevel scheme for simulating turbulent channel flows. Two different versions of the scheme, namely the spectral and finite difference versions, are presented. The spectral version of the scheme is based on a spectral-Galerkin formulation which provides a natural decomposition of the flow into small and large wavelength parts, and which leads to linear systems that can be solved with quasi-optimal computational complexity. In the finite difference version, the Incremental Unknown (IU) is used to separate the length scales. Preliminary numerical results indicate that the scheme is well suited for turbulence computations and provides results which are comparable to that by Direct Numerical Simulation (DNS) but with significantly less CPU time.

Chen, M.↗

A Numerical Study of Self-Similarity in a Turbulent Plane Wake Using Large-Eddy Simulation

Turbulent wakes are known to develop self-similarly sufficiently far downstream from obstacles that generate them. It has long been assumed that the spreading rate of the wake in the self-similar regime is independent of the details of the body generating the wake, being dependent only on the total drag (or momentum deficit). This assumption seems to be in contradiction with some recent experiments. In this study we attempt to complement these experimental investigations through a numerical study of a time-developing wake. A numerical study has the advantage of eliminating many of the uncontrolled factors present in experiments and allowing precise control of initial conditions. Large-eddy simulations employing the recently developed dynamic localization model are used to extend previous results from direct numerical simulations. The large-eddy simulation results are compared to the direct numerical simulation database, wherever such comparisons are feasible, as a check of the method. Like the experiments, the large-eddy simulations suggest that non-unique self-similar states, characterized by different spreading rates and turbulent statistics, are possible and that they can be maintained for significant time periods. The study also demonstrates the predictive capability of the dynamic localization subgrid model.

Ghosal, Sandip↗

Study of Lagrangian characteristic times using direct numerical simulation of turbulence

Direct numerical simulations and large eddy simulations of homogeneous isotropic turbulence are used to compute Lagrangian statistics of turbulence and, in particular, its time scales. The computed time scales are compared with the spectral times scales that are frequently used in Eddy Damped Quasi-Normal Markovian calculations of the spectrum. The time scale models are rather good at high wavenumber and the results point to directions for improvement of the time scales at low wavenumber.

Lee, C. H.↗

A regularized deconvolution method for turbulent closure modeling in implicitly filtered large-eddy simulation

Turbulence closure models that are based on a regularized deconvolution method (RDM) are proposed for application to implicitly filtered Large-Eddy Simulation (LES). This method reconstructs sub-grid scale contributions by approximately recovering filtered quantities through an optimization approach. Physical principles on the boundedness and conservation of reconstructed quantities are preserved by constraining the optimization problem. Further, to account for effects of scales that are not resolved in implicitly filtered LES, a reconstruction method is proposed that combines sub-grid projection and energy-similarity to enrich unresolved scales. With this method, closure models are developed for turbulent scalar fluxes and turbulence–chemistry interaction that are intrinsically consistent with the LES formulation. These closure models are evaluated in application to a piloted turbulent jet flame with inhomogeneous inlet composition, in which combustion is represented by a two-scalar manifold model. The RDM-formulation is benchmarked against simulations obtained using a presumed PDF method and eddy diffusivity closure model. In partially-premixed flame regions, it is observed that RDM improves the prediction of turbulent mixing, and faster grid convergence is observed for results with RDM compared to simulations employing a presumed PDF method.

42 ENGINEERING↗

Advances in turbulent transport modeling based on direct simulations of turbulence

Recent work using numerical simulations to study the modeling of scalar transport is outlined. A model developed by reference to simulations for homogeneous shear flow is presented. This model predicts cross-gradient transports that can give rise to the illusion of counter-gradient transport. The model is found to work remarkably well when compared to simulations for channel flow, but fails to predict important features of the scalar fluxes in mixing layers. Reasons for this failure are discussed and suggestions are made for future work.

Reynolds, W. C.↗

Development and application of an atmospheric turbulence model for use in flight simulators in flight simulators

The influence of simulated turbulence on aircraft handling qualities was investigated. Pilot opinion of the handling qualities of a light general aviation aircraft were evaluated in a motion-base simulator using a simulated turbulence environment. A realistic representation of turbulence disturbances is described in terms of rms intensity and scale length and their random variations with time. The time histories generated by the proposed turbulence models showed characteristics which appear to be more similar to real turbulence than the frequently-used Gaussian turbulence model. In addition, the proposed turbulence models can flexibly accommodate changes in atmospheric conditions and be easily implemented in flight simulator studies. Six turbulence time histories, including the conventional Gaussian model, were used in an IFR-tracking task. The realism of each of the turbulence models and the handling qualities of the simulated airplane were evaluated. Analysis of pilot opinions shows that at approximately the same rms intensities of turbulence, the handling quality ratings transit from the satisfactory level, for the simple Gaussian model, to an unacceptable level for more realistic and compositely structured turbulence models.

Jacobson, I. D.↗

A new method for simulating atmospheric turbulence for rotorcraft applications

Simulation of atmospheric turbulence as seen by a rotating blade element involves treatment of cyclostationary processes. Conventional filtering techniques do not lend themselves well to the generation of such turbulence sample functions as are required in rotorcraft flight dynamics simulation codes. A method to generate sample functions containing second order statistics of mean and covariance is presented. Compared to ensemble averaging involving excessive computer time, the novelty is to exploit cycloergodicity and thereby, replace ensemble averaging by averaging over a single path sample function of long duration. The method is validated by comparing its covariance results with the analytical and ensemble averaged results for a widely used 1-D turbulence approximation.

J. Riaz↗

A new method for simulating atmospheric turbulence for rotorcraft applications

Simulation of atmospheric turbulence as seen by a rotating blade element involves treatment of cyclostationary processes. Conventional filtering techniques do not lend themselves well to the generation of such turbulence sample functions as are required in rotorcraft flight dynamics simulation codes. A method to generate sample functions containing second-order statistics of mean and covariance is presented. Compared to ensemble averaging involving excessive computer time, the novelty is to exploit cycloergodicity and thereby, replace ensemble averaging by averaging over a single-path sample function of long duration. The method is validated by comparing its covariance results with the analytical and ensemble-averaged results for a widely used one-dimensional turbulence approximation.

Prasad, J. V. R.↗

Direct numerical simulation and Reynolds-averaged Navier-Stokes modeling of the sudden viscous dissipation for multicomponent turbulence

Simulations of a turbulent multicomponent fluid mixture undergoing isotropic deformations are carried out to investigate the sudden viscous dissipation. This dissipative mechanism was originally demonstrated using simulations of an incompressible single-component fluid [S. Davidovits and N. J. Fisch, Phys. Rev. Lett. 116, 105004 (2016)]. By accounting for the convective and diffusive transfer of various species, the current work aims to increase the physical fidelity of previous simulations and their relevance to inertial confinement fusion applications. Direct numerical simulations of the compressed fluid show that the sudden viscous dissipation of turbulent kinetic energy is unchanged from the single-component scenario. More importantly, the simulations demonstrate that the mass fraction variance and covariance for the various species also exhibit a sudden viscous decay. Reynolds-averaged Navier-Stokes simulations were carried out using the $k - l$ model to assess its ability to reproduce the sudden viscous dissipation. Results show that the standard $k - l$ formulation does not capture the sudden decay of turbulent kinetic energy, mass-fraction variance, and mass-fraction covariance for simulations with various compression and expansion rates, or different exponents for the power-law model of viscosity. A new formulation of the $k - l$ model that is based on previous improvements to the $k - ε$ family of models is proposed, which leads to consistently good agreement with the direct numerical simulations for all the isotropic deformations under consideration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulations of turbulent premixed H2/O2/N2 flames with complex chemistry

Premixed stoichiometric H2/O2/N2 flames propagating in two-dimensional turbulence were studied using direct numerical simulation (simulations in which all fluid and thermochemical scales are fully resolved) including realistic chemical kinetics and molecular transport. Results are compared with earlier zero-chemistry (flame sheet) and one-step chemistry simulations. Consistent with the simpler models, the turbulent flame with realistic chemistry aligns preferentially with extensive strain rates in the tangent plane and flame curvature probability density functions are close to symmetric with near-zero means. By contrast to simple-chemistry results with non-unity Lewis numbers (ratio of thermal to species diffusivity), local flame structure does not correlate with curvature but rather with tangential strain rate. Turbulent straining results in substantial thinning of the flame relative to the steady unstrained laminar case. Heat release and H2O2 contours remain thin and connected ('flamelet-like') while species including H-atom and OH are more diffuse. Peak OH concentration occurs well behind the peak heat-release zone. The feasibility of incorporating realistic chemistry into full turbulence simulations to address issues such as pollutant formation in hydrocarbon-air flames is suggested.

Baum, M.↗

Simulations of turbulent mixing and reacting flows and their applications to turbulence modeling

The method of full simulation is applied to reacting turbulent flows. Full simulation has proven of great value as a complement to experiments for the study of nonreacting turbulent flows. It provides insight into the physics of turbulent flows and their modeling. It is natural to try to extend these methods to the simulation of reacting turbulent flows. Because this is one of the first attempts at this type of simulation, a subsidiary goal of this work is to demonstrate the feasibility of using simulation to study turbulent reacting flows. In addition, it is shown that such simulations can be used to provide physical insight into the nature of turbulent combustion and to provide data that will help to construct models that can be used in engineering simulations of turbulent reacting flows.

Ferziger, J. H.↗

The effect of artificial bulk viscosity in simulations of forced compressible turbulence

The use of an artificial bulk viscosity for shock stabilization is a common approach employed in turbulence simulations with high-order numerics. The effect of the artificial bulk viscosity is analyzed in the context of large eddy simulations by using as a test case simulations of linearly-forced compressible homogeneous turbulence (Petersen and Livescu, 2010 [12]). This case is unique in that it allows for the specification of a priori target values for total dissipation and ratio of solenoidal to dilatational dissipation. A comparison between these target values and the true predicted levels of dissipation is thus used to investigate the performance of the artificial bulk viscosity. Results show that the artificial bulk viscosity is effective at achieving stable solutions, but also leads to large values of artificial dissipation that outweigh the physical dissipation caused by fluid viscosity. An alternate approach, which employs the artificial thermal conductivity only, shows that the dissipation of dilatational modes is entirely due to the fluid viscosity. However, this method leads to unwanted Gibbs oscillations around the shocklets. The use of shock sensors that further localize the artificial bulk viscosity did not reduce the amount of artificial dissipation introduced by the artificial bulk viscosity. Finally, an improved forcing function that explicitly accounts for the role of the artificial bulk viscosity in the budget of turbulent kinetic energy was explored.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulations of transition and turbulence on the Navier-Stokes computer

The Navier-Stokes Computer (NSC) consists of multiple local memory parallel processors interconnected in a hypercube network. Efficient implementation of algorithms on the NSC thus requires the effective utilization of both the coarse and fine grain paralelism inherent in the architectural design. The basic approach to implementing an algorithm on the NSC is presented herein. The particular finite-difference algorithm considered was developed for performing transition and turbulence simulations by direct solution of the time-dependent incompressible Navier-Stokes equations. The suitability of this algorithm for performing simulations of the isotropic turbulence problem is verified from computations performed on a Cray 2. Projected timing results for the algorithm on the NSC itself are presented for both the isotropic turbulence and laminar turbulent transition problems.

Krist, S. E.↗

Outcomes from the DOE Workshop on Turbulent Flow Simulation at the Exascale

This paper summarizes the outcomes from the Turbulent Flow Simulation at the Exascale: Opportunities and Challenges Workshop, which was held 4-5 August 2015, and was sponsored by the U.S. Department of Energy Office of Advanced Scientific Computing Research. The workshop objective was to define and describe the challenges and opportunities that computing at the exascale will bring to turbulent-flow simulations in applied science and technology. The need for accurate simulation of turbulent flows is evident across the U.S. Department of Energy applied-science and engineering portfolios, including combustion, plasma physics, nuclear-reactor physics, wind energy, and atmospheric science. The workshop brought together experts in turbulent-flow simulation, computational mathematics, and high-performance computing. Building upon previous ASCR workshops on exascale computing, participants defined a research agenda and path forward that will enable scientists and engineers to continually leverage, engage, and direct advances in computational systems on the path to exascale computing.

97 MATHEMATICS AND COMPUTING↗

Cross-separatrix simulations of turbulent transport in the field-reversed configuration

Recent local simulations of the field-reversed configuration (FRC) have reported drift-wave stability in the core and instability in the scrape-off layer (SOL). However, experimental measurements indicate the existence of fluctuations in both FRC core and SOL, with much lower amplitude fluctuations measured in the core. With the updated cross-separatrix capabilities of the simulation code used in this paper, nonlinear turbulence simulations find that linear instabilities grow in the SOL, generating fluctuations which spread from SOL to core. After saturation of the linear instabilities, a balance of the inward spread and local damping in the core is achieved. The steady state toroidal wavenumber spectrum shows lower amplitude core fluctuations and larger SOL fluctuations with amplitude decreasing towards shorter wavelengths, which are consistent with experimental measurements.

Lau, C. K.↗