Numerical studies towards practical large-eddy simulation

Large-eddy simulation developments and validations are presented for an improved simulation of turbulent internal flows. Numerical methods are proposed according to two competing criteria: numerical qualities (precision and spectral characteristics), and adaptability to complex configurations. First, methods are tested on academic test-cases, in order to abridge with fundamental studies. Consistent results are obtained using adaptable finite volume method, with higher order advection fluxes, implicit grid filtering and "low-cost" shear-improved Smagorinsky model. This analysis particularly focuses on mean flow, fluctuations, two-point correlations and spectra. Moreover, it is shown that exponential averaging is a promising tool for LES implementation in complex geometry with deterministic unsteadiness. Finally, adaptability of the method is demonstrated by application to a configuration representative of blade-tip clearance flow in a turbomachine

Anomalous Scaling of Structure Functions and Dynamic Constraints on Turbulence Simulations

The connection between anomalous scaling of structure functions (intermittency) and numerical methods for turbulence simulations is discussed. It is argued that the computational work for direct numerical simulations (DNS) of fully developed turbulence increases as $Re^{4}$, and not as $Re^{3}$ expected from Kolmogorov's theory, where $Re$ is a large-scale Reynolds number. Various relations for the moments of acceleration and velocity derivatives are derived. An infinite set of exact constraints on dynamically consistent subgrid models for Large Eddy Simulations (LES) is derived from the Navier-Stokes equations, and some problems of principle associated with existing LES models are highlighted.Comment: 18 page

GPU-Accelerated Large-Eddy Simulation of Turbulent Channel Flows

High performance computing clusters that are augmented with cost and power efficient graphics processing unit (GPU) provide new opportunities to broaden the use of large-eddy simulation technique to study high Reynolds number turbulent flows in fluids engineering applications. In this paper, we extend our earlier work on multi-GPU acceleration of an incompressible Navier-Stokes solver to include a large-eddy simulation (LES) capability. In particular, we implement the Lagrangian dynamic subgrid scale model and compare our results against existing direct numerical simulation (DNS) data of a turbulent channel flow at Reτ = 180. Overall, our LES results match fairly well with the DNS data. Our results show that the Reτ = 180 case can be entirely simulated on a single GPU, whereas higher Reynolds cases can benefit from a GPU cluster

Turbulent Compressible Convection with Rotation - Penetration above a Convection Zone

We perform Large eddy simulations of turbulent compressible convection in stellar-type convection zones by solving the Navi\'{e}r-Stokes equations in three dimensions. We estimate the extent of penetration into the stable layer above a stellar-type convection zone by varying the rotation rate ({\boldmath$\rm\Omega$}), the inclination of the rotation vector ($\theta$) and the relative stability ($S$) of the upper stable layer. The computational domain is a rectangular box in an f-plane configuration and is divided into two regions of unstable and stable stratification with the stable layer placed above the convectively unstable layer. Several models have been computed and the penetration distance into the stable layer above the convection zone is estimated by determining the position where time averaged kinetic energy flux has the first zero in the upper stable layer. The vertical grid spacing in all the model is non-uniform, and is less in the upper region so that the flows are better resolved in the region of interest. We find that the penetration distance increases as the rotation rate increases for the case when the rotation vector is aligned with the vertical axis. However, with the increase in the stability of the upper stable layer, the upward penetration distance decreases. Since we are not able to afford computations with finer resolution for all the models, we compute a number of models to see the effect of increased resolution on the upward penetration. In addition, we estimate the upper limit on the upward convective penetration from stellar convective cores.Comment: Accepted for Publication in Asttrophysics & Space Scienc

Dust Devil Sediment Transport: From Lab to Field to Global Impact

The impact of dust aerosols on the climate and environment of Earth and Mars is complex and forms a major area of research. A difficulty arises in estimating the contribution of small-scale dust devils to the total dust aerosol. This difficulty is due to uncertainties in the amount of dust lifted by individual dust devils, the frequency of dust devil occurrence, and the lack of statistical generality of individual experiments and observations. In this paper, we review results of observational, laboratory, and modeling studies and provide an overview of dust devil dust transport on various spatio-temporal scales as obtained with the different research approaches. Methods used for the investigation of dust devils on Earth and Mars vary. For example, while the use of imagery for the investigation of dust devil occurrence frequency is common practice for Mars, this is less so the case for Earth. Modeling approaches for Earth and Mars are similar in that they are based on the same underlying theory, but they are applied in different ways. Insights into the benefits and limitations of each approach suggest potential future research focuses, which can further reduce the uncertainty associated with dust devil dust entrainment. The potential impacts of dust devils on the climates of Earth and Mars are discussed on the basis of the presented research results

Orbital Observations of Dust Lofted by Daytime Convective Turbulence

Over the past several decades, orbital observations of lofted dust have revealed the importance of mineral aerosols as a climate forcing mechanism on both Earth and Mars. Increasingly detailed and diverse data sets have provided an ever-improving understanding of dust sources, transport pathways, and sinks on both planets, but the role of dust in modulating atmospheric processes is complex and not always well understood. We present a review of orbital observations of entrained dust on Earth and Mars, particularly that produced by the dust-laden structures produced by daytime convective turbulence called “dust devils”. On Earth, dust devils are thought to contribute only a small fraction of the atmospheric dust budget; accordingly, there are not yet any published accounts of their occurrence from orbit. In contrast, dust devils on Mars are thought to account for several tens of percent of the planet’s atmospheric dust budget; the literature regarding martian dust devils is quite rich. Because terrestrial dust devils may temporarily contribute significantly to local dust loading and lowered air quality, we suggest that martian dust devil studies may inform future studies of convectively-lofted dust on Earth