On local approximations of the pressure-strain term in turbulence models

The results of numerical simulations of turbulent channel flows were used to examine the validity of the local approximation of the pressure-strain term in the Reynolds stress transport equation. Outside of the viscous sublayer the local approximation compares very well with the exact pressure strain. This agreement is due, at least in part, to the high correlation between the rapid pressure and its Laplacian, which suggests that only the near parts of the flow contribute to the rapid pressure at a point. In the viscous sublayer the distance over which the mean shear can be considered constant is comparable to the length scale in the normal direction of the correlations of velocity gradients, leading to failure of the local approximation

Dynamic subfilter-scale stress model for large-eddy simulations

We present a modification of the integral length-scale approximation (ILSA) model originally proposed by Piomelli et al. [Piomelli et al., J. Fluid Mech. 766, 499 (2015)] and apply it to plane channel flow and a backward-facing step. In the ILSA models the length scale is expressed in terms of the integral length scale of turbulence and is determined by the flow characteristics, decoupled from the simulation grid. In the original formulation the model coefficient was constant, determined by requiring a desired global contribution of the unresolved subfilter scales (SFSs) to the dissipation rate, known as SFS activity; its value was found by a set of coarse-grid calculations. Here we develop two modifications. We de-fine a measure of SFS activity (based on turbulent stresses), which adds to the robustness of the model, particularly at high Reynolds numbers, and removes the need for the prior coarse-grid calculations: The model coefficient can be computed dynamically and adapt to large-scale unsteadiness. Furthermore, the desired level of SFS activity is now enforced locally (and not integrated over the entire volume, as in the original model), providing better control over model activity and also improving the near-wall behavior of the model. Application of the local ILSA to channel flow and a backward-facing step and comparison with the original ILSA and with the dynamic model of Germano et al. [Germano et al., Phys. Fluids A 3, 1760 (1991)] show better control over the model contribution in the local ILSA, while the positive properties of the original formulation (including its higher accuracy compared to the dynamic model on coarse grids) are maintained. The backward-facing step also highlights the advantage of the decoupling of the model length scale from the mesh

Accuracy control for large-eddy simulation of turbulent mixing - Integral length-scale approach

Turbulent flow at high Reynolds numbers is currently not accessible on the basis of direct numerical simulation (DNS) of the Navier-Stokes equations - the computational complexity is too high to allow DNS in most realistic flow conditions. Instead, Large-Eddy Simulation (LES) offers an alternative in which the focus is on capturing the larger dynamic scales of a problem. However, the fundamental closure problem in LES induced by spatial filtering of nonlinear terms, and the role of discretization errors in the numerical treatment of the LES equations, induce a principal uncertainty in any LES prediction. This uncertainty requires quantification and control. We investigate error control capabilities of the Integral Length-Scale Approximation (ILSA) and apply this modeling to transitional and turbulent mixing, focussing on the achieved reliability of LES as function of the grid resolution and ‘sub-filter activity’.</p

Nested Dissection Method on Transputer

Nested dissection method is an elimination method for a set of linear algebraic equations with minimum fillins. Physically it divides a domain into four subdomains, and each subdomain is again divided into four. This procedure is repeated till all nodes are included in some subdomains. Using this characteristic, the authors examine the efficiency of the method on the transputer

Immersed boundary method for generalised finite volume and finite difference Navier-Stokes solvers

In Immersed Boundary Methods (IBM) the effect of complex geometries is introduced through the forces added in the Navier-Stokes solver at the grid points in the vicinity of the immersed boundaries. Most of the methods in the literature have been used with Cartesian grids. Moreover many of the methods developed in the literature do not satisfy some basic conservation properties (the conservation of torque, for instance) on non-uniform meshes. In this paper we will follow the RKPM method originated by Liu et al. [1] to build locally regularized functions that verify a number of integral conditions. These local approximants will be used both for interpolating the velocity field and for spreading the singular force field in the framework of a pressure correction scheme for the incompressible Navier-Stokes equations. We will also demonstrate the robustness and effectiveness of the scheme through various examples

Experimental assessment of a new form of scaling law for near-wall turbulence

Scaling laws and intermittency in the wall region of a turbulent flow are addressed by analyzing moderate Reynolds number data obtained by single component hot wire anemometry in the boundary layer of a flat plate. The paper aims in particular at the experimental validation of a new form of refined similarity recently proposed for the shear dominated range of turbulence, where the classical Kolmogorov-Oboukhov inertial range theory is inappropriate. An approach inspired to the extended self-similarity allows for the extraction of the different power laws for the longitudinal structure functions at several wall normal distances. A double scaling regime is found in the logarithmic region, confirming previous experimental results. Approaching the wall, the scaling range corresponding to the classical cascade-dominated range tends to disappear and, in the buffer layer, a single power law is found to describe the available range of scales. The double scaling is shown to be associated with two different forms of refined similarity. The classical form holds below the shear scale L s . The other, originally introduced on the basis of DNS data for a turbulent channel, is experimentally confirmed to set up above L s . Given the experimental diffulties in the evaluation of the instantaneous dissipation rate, some care is devoted to check that its one-dimensional surrogate does not bias the results. The increased intermittency as the wall is approached is experimentally found entirely consistent with the failure of the refined Kolmogorov-Oboukhov similarity and the establishment of its new form near the wall.Comment: 27 pages, 9 figure

Large-Eddy Simulation: Current Capabilities, Recommended Practices, and Future Research

This paper presents the results of an activity by the Large Eddy Simulation (LES) Working Group of the AIAA Fluid Dynamics Technical Committee to (1) address the current capabilities of LES, (2) outline recommended practices and key considerations for using LES, and (3) identify future research needs to advance the capabilities and reliability of LES for analysis of turbulent flows. To address the current capabilities and future needs, a survey comprised of eleven questions was posed to LES Working Group members to assemble a broad range of perspectives on important topics related to LES. The responses to these survey questions are summarized with the intent not to be a comprehensive dictate on LES, but rather the perspective of one group on some important issues. A list of recommended practices is also provided, which does not treat all aspects of a LES, but provides guidance on some of the key areas that should be considered

Experimental study of Taylor's hypothesis in a turbulent soap film

An experimental study of Taylor's hypothesis in a quasi-two-dimensional turbulent soap film is presented. A two probe laser Doppler velocimeter enables a non-intrusive simultaneous measurement of the velocity at spatially separated points. The breakdown of Taylor's hypothesis is quantified using the cross correlation between two points displaced in both space and time; correlation is better than 90% for scales less than the integral scale. A quantitative study of the decorrelation beyond the integral scale is presented, including an analysis of the failure of Taylor's hypothesis using techniques from predictability studies of turbulent flows. Our results are compared with similar studies of 3D turbulence.Comment: 27 pages, + 19 figure