Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 163 records · Page 9

Turbulent Flow Simulation at the Exascale: Opportunities and Challenges Workshop: August 4-5, 2015, Washington, D.C.

This report details the impact exascale will bring to turbulent-flow simulations in applied science and technology. The need for accurate simulation of turbulent flows is evident across the DOE 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↗

Molecular-level simulations of turbulence and its decay

Here, we provide the first demonstration that molecular-level methods based on gas kinetic theory and molecular chaos can simulate turbulence and its decay. The direct simulation Monte Carlo (DSMC) method, a molecular-level technique for simulating gas flows that resolves phenomena from molecular to hydrodynamic (continuum) length scales, is applied to simulate the Taylor-Green vortex flow. The DSMC simulations reproduce the Kolmogorov –5/3 law and agree well with the turbulent kinetic energy and energy dissipation rate obtained from direct numerical simulation of the Navier-Stokes equations using a spectral method. This agreement provides strong evidence that molecular-level methods for gases can be used to investigate turbulent flows quantitatively.

74 ATOMIC AND MOLECULAR PHYSICS↗

Simulation of turbulent shear flows at Stanford and NASA-Ames - What can we do and what have we learned?

The capabilities and limitations of large eddy simulation (LES) and full turbulence simulation (FTS) are outlined. It is pointed out that LES, although limited at the present time by the need for periodic boundary conditions, produces large-scale flow behavior in general agreement with experiments. What is more, FTS computations produce small-scale behavior that is consistent with available experiments. The importance of the development work being done on the National Aerodynamic Simulator is emphasized. Studies at present are limited to situations in which periodic boundary conditions can be applied on boundaries of the computational domain where the flow is turbulent.

Reynolds, W. C.↗

WHISTLER TURBULENCE FORWARD CASCADE VERSUS INVERSE CASCADE: THREE-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS

We present the results of the first fully three-dimensional particle-in-cell simulations of decaying whistler turbulence in a magnetized, homogeneous, collisionless plasma in which both forward cascades to shorter wavelengths, and inverse cascades to longer wavelengths are allowed to proceed. For the electron beta β {sub e} = 0.10 initial value considered here, the early-time rate of inverse cascade is very much smaller than the rate of forward cascade, so that at late times the fluctuation energy in the regime of the inverse cascade is much weaker than that in the forward cascade regime. Similarly, the wavevector anisotropy in the inverse cascade regime is much weaker than that in the forward cascade regime.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct simulation of compressible turbulence

Several direct simulations of 3-D homogeneous, compressible turbulence are presented with emphasis on the differences with incompressible turbulent simulations. A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order Runge-Kutta time discretization scheme is sufficient to produce well-resolved flows at Taylor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of velocity is useful to differentiate between the purely compressible effects and those effects solely due to vorticity production. In the context of homogeneous flows, this decomposition in unique. Time-dependent energy and dissipation spectra of the compressible and solenoidal velocity components indicate the presence of localized small scale structures. These structures are strongly a function of the initial conditions. Researchers concentrate on a regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) and density and temperature fluctuations much greater than sq Ma. This leads to a state in which more than 70 percent of the kinetic energy is contained in the so-called compressible component of the velocity. Furthermore, these conditions lead to the formation of curved weak shocks (or shocklets) which travel at approximately the sound speed across the physical domain. Various terms in the vorticity and divergence of velocity production equations are plotted versus time to gain some understanding of how small scales are actually formed. Possible links with Burger turbulence are examined. To visualize better the dynamics of the flow, new graphic visualization techniques have been developed. The 3-D structure of the shocks are visualized with the help of volume rendering algorithms developed in-house. A combination of stereographic projection and animation greatly increase the number of visual cues necessary to properly interpret the complex flow.

Zang, T. A.↗

Large-eddy simulations of turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines rotating at low tip speed ratios

Turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines (VAWTs) rotating at low tip speed ratios (TSRs) are studied numerically. The turbulent flows are simulated using the large-eddy simulation (LES) model, and the rotating turbine blades are modeled using the actuator line method. The helical VAWT has identical key parameters as the straight VAWT except for the 135° helical twist of the blades over the 0.3m vertical span. A set of LES runs are performed for two TSRs, 0.6 and 0.4, and the results are reported and analyzed. At these low TSRs, the wake behind the straight-bladed VAWT exhibits two-dimensional dominant flow motions (in the horizontal plane perpendicular to the straight blades) in the near-wake region that cause considerable spanwise expansion of the wake as it extends downstream. In contrast, the helical-bladed VAWT generates highly three-dimensional (3D) wake flow structures and upward/downward mean flow motions within the wake that cause the wake to expand mainly in the vertical direction. Turbulence statistical analyses also show that the 3D wake flow features induced by the helical blades accelerate the wake transition to turbulence and enhance the small-scale turbulent dissipation (as shown by the subgrid-scale turbulent dissipation in the LES), which leads to a more rapid decay of the wake turbulence intensity than that in the straight-bladed VAWT case at the same TSR. Finally, compared with the straight-bladed VAWT, the helical-bladed VAWT also exhibits much smaller temporal variations for the torque and power coefficients during the rotation cycle, which can be beneficial for wind power generation.

17 WIND ENERGY↗

Numerical simulation of turbulent jet noise, part 2

Results on the numerical simulation of jet flow fields were used to study the radiated sound field, and in addition, to extend and test the capabilities of the turbulent jet simulation codes. The principal result of the investigation was the computation of the radiated sound field from a turbulent jet. In addition, the computer codes were extended to account for the effects of compressibility and eddy viscosity, and the treatment of the nonlinear terms of the Navier-Stokes equations was modified so that they can be computed in a semi-implicit way. A summary of the flow model and a description of the numerical methods used for its solution are presented. Calculations of the radiated sound field are reported. In addition, the extensions that were made to the fundamental dynamical codes are described. Finally, the current state-of-the-art for computer simulation of turbulent jet noise is summarized.

Metcalfe, R. W.↗

Direct Numerical Simulations of Turbulent Flow over Misaligned Traveling Waves

Direct numerical simulations of turbulent flow over prescribed traveling waves in a half-channel flow setup subject to a combined streamwise and spanwise pressure gradient are performed at a friction Reynolds number of 180. The simulations undertaken in this study consider both flow-aligned (including opposing and following) and misaligned waves as we attempt to quantify the effect of wave misalignment on turbulence statistics, the pressure drag force, and the wave attenuation rate. For the simulations, we consider three characteristic wave age values corresponding to slow-, intermediate-, and fast-moving waves. Wave misalignment is taken into account by applying a spanwise pressure gradient vertical to the traveling waves, which results in a three-dimensional turbulent flow field above the moving waves. Key flow quantities such as the mean velocity, velocity variances, and momentum fluxes are found to vary with the wave parameters, confirming the findings of previous studies.

direct numerical simulation↗

Numerical simulation of turbulent trailing edge flows

Numerical simulations of the time-dependent, Reynolds-averaged, Navier-Stokes equations, employing a two-equation turbulence model, are presented and compared with measurements from a series of trailing edge experiments at transonic Mach numbers. The test flows include an asymmetric flow with no separation, an asymmetric flow with a small region of separation and a symmetric flow with a large shock-wave induced separated zone. Comparisons are made for mean surface quantities as well as for mean and fluctuating flow-field quantities. For the trailing-edge flows with little or no separation, the solutions correctly predict all the major features of the flow field. Treatment of the viscous-inviscid interaction was found to be important for predicting these test cases. Two-equation eddy-viscosity turbulence models were found to be adequate for these flows. However, for the shock-wave induced separation case, these turbulence models were inadequate to predict this flow field. Modifications of the turbulence model to correct these deficiencies are discussed.

Horstman, C. C.↗

On direct numerical simulations of turbulent reacting flows

A description of the emerging field of direct numerical simulations of turbulent, chemically reacting flows is presented. The types of direct numerical simulations, physical issues related to implementing the simulations, as well as the various numerical methods used are described. Examples are presented of recent applications of direct numerical simulations to a variety of problems, displaying both the potential of the method and also some of its limitations. Finally, our view of the potential role of direct numerical simulations in future research on turbulent, chemically reacting flows is presented.

Jou, W.-H.↗

Grid resolution requirement for resolving rare and high intensity wall-shear stress events in direct numerical simulations

Turbulent signals are intermittent with large instantaneous fluctuations. Such large fluctuations lead to small Kolmogorov scales that are hard to resolve in numerical simulations [P. K. Yeung, K. R. Sreenivasan, and S. B. Pope, Effects of finite spatial and temporal resolution in direct numerical simulations of incompressible isotropic turbulence, Phys. Rev. Fluids 3, 064603 (2018)]. Here, the present paper follows the above basic logic, but instead of dissipation events in isotropic turbulence, we study wall-shear stress events in plane channel flow. Wall-shear stress fluctuations are increasingly more intermittent as the Reynolds number increases. Hence, one has to employ higher grid resolutions as the Reynolds number increases in order to resolve a given percentage of wall-shear stress events. The objective of this paper is to quantify effects of the grid resolutions on the rare and high intensity wall-shear stress events. We find that the standard grid resolution resolves about 99% of the wall-shear stress events at Re τ = 180. A slightly higher grid resolution has to be employed in order to resolve 99% of the wall-shear stress events at higher Reynolds numbers, and if the standard grid resolution is used for, e.g., a Re τ = 10000 channel flow, one resolves about 90%–95% wall-shear stress events.

42 ENGINEERING↗

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗

Gyrokinetic Particle Simulation of Turbulent Transport in Burning Plasmas

The SciDAC project at the IFS advanced the state of high performance computing for turbulent structures and turbulent transport. The team project with Prof Zhihong Lin [PI] at Univ California Irvine produced new understanding of the turbulent electron transport. The simulations were performed at the Texas Advanced Computer Center TACC and the NERSC facility by Wendell Horton, Lee Leonard and the IFS Graduate Students working in that group. The research included a Validation of the electron turbulent transport code using the data from a steady state university experiment at the University of Columbia in which detailed probe measurements of the turbulence in steady state were used for wide range of temperature gradients to compare with the simulation data. These results were published in a joint paper with Texas graduate student Dr. Xiangrong Fu using the work in his PhD dissertation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large-Eddy Simulation of Turbulent Wall-Pressure Fluctuations

Large-eddy simulations of a turbulent boundary layer with Reynolds number based on displacement thickness equal to 3500 were performed with two grid resolutions. The computations were continued for sufficient time to obtain frequency spectra with resolved frequencies that correspond to the most important structural frequencies on an aircraft fuselage. The turbulent stresses were adequately resolved with both resolutions. Detailed quantitative analysis of a variety of statistical quantities associated with the wall-pressure fluctuations revealed similar behavior for both simulations. The primary differences were associated with the lack of resolution of the high-frequency data in the coarse-grid calculation and the increased jitter (due to the lack of multiple realizations for averaging purposes) in the fine-grid calculation. A new curve fit was introduced to represent the spanwise coherence of the cross-spectral density.

Singer, Bart A.↗

Global gyrokinetic simulations of kinetic-ballooning-mode turbulence in Wendelstein 7-X

We apply the global gyrokinetic code GENE-3D to perform kinetic-ballooning-mode (KBM) turbulence simulations in the Wendelstein 7-X (W7-X) stellarator using experimentally relevant parameters. In contrast to the flux-tube approach, we observe saturation of the turbulent fluctuations at finite amplitude. Analyzing the turbulent characteristics, we identify the presence of weakly driven KBM turbulence hybridizing with ion-temperature-gradient and trapped-electron-mode turbulence in the scenario under consideration, showing that it is possible that KBM activity can be driven in W7-X, opening up the path for studying this type of turbulence, which is highly relevant for future fusion power plants.

Electrostatics↗

Outer midplane scrape-off layer profiles and turbulence in simulations of Alcator C-Mod inner-wall limited discharges

A region of steep plasma gradients, the so-called ”narrow-feature”, has been found in the near scrape-off layer (SOL) of inner-wall limited (IWL) discharges. Dedicated IWL discharges were carried out in Alcator C-Mod [E.S. Marmar et al., Nucl. Fusion 55, (2015)] to study this phenomenon, allowing detailed observations of the plasma profiles and fluctuations. Langmuir probe (LP) measurements show a clear two decay length n e and T e profile structure at the outer midplane. The Gas-Puff Imaging (GPI) diagnostic shows large turbulent fluctuations across the last closed flux-surface, hence supporting the hypothesis that turbulent phenomena play a role in setting the profile steepness. We have carried out the flux-driven non-linear turbulence simulations of two C-Mod discharges which allows a three-way comparison between LP, GPI, and simulation data. Observations and simulations correlate the steep gradient region characterizing the narrow feature with sheared poloidal flows and a deviation of the plasma potential from its floating value. Furthermore, the E x B shear rate exceeds the linear ballooning growth rate, indicating that the narrow feature could result from the effects of sheared flows, although causality could not be established. The fluctuation level in the narrow feature remains of order unity across the entire SOL, indicating that the transport reduction in the near-SOL cannot result from a simple quench rule.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The use of direct numerical simulation in the study of turbulent, chemically-reacting flows

Full turbulence simulations are used here to study the effects of chemical heat release on the large-scale structures in turbulent mixing layers. In agreement with laboratory results, it is found that the heat release lowers the rate at which the mixing layer grows and reduces the rate at which chemical products are formed. The baroclinic torque and thermal expansion in the mixing layer produce changes in the flame vortex structure that act to produce more diffuse vortices than in the constant density case, resulting in lower rotation rates of the large-scale structures. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layer are shown to result from vorticity generation by baroclinic torques. Calculations of the energy in the various wavenumbers shows that the heat release has a stabilizing effect on the growth rate of individual modes. This methodology can be applied for Reynolds numbers less than several hundred and for Damkoehler numbers less than about ten.

Riley, J. J.↗