Statistical properties of an ideal subgrid-scale correction for Lagrangian particle tracking in turbulent channel flow

One issue associated with the use of Large-Eddy Simulation (LES) to investigate the dispersion of small inertial particles in turbulent flows is the accuracy with which particle statistics and concentration can be reproduced. The motion of particles in LES fields may differ significantly from that observed in experiments or direct numerical simulation (DNS) because the force acting on the particles is not accurately estimated, due to the availability of the only filtered fluid velocity, and because errors accumulate in time leading to a progressive divergence of the trajectories. This may lead to different degrees of inaccuracy in the prediction of statistics and concentration. We identify herein an ideal subgrid correction of the a-priori LES fluid velocity seen by the particles in turbulent channel flow. This correction is computed by imposing that the trajectories of individual particles moving in filtered DNS fields exactly coincide with the particle trajectories in a DNS. In this way the errors introduced by filtering into the particle motion equations can be singled out and analyzed separately from those due to the progressive divergence of the trajectories. The subgrid correction term, and therefore the filtering error, is characterized in the present paper in terms of statistical moments. The effects of the particle inertia and of the filter type and width on the properties of the correction term are investigated.Comment: 15 pages,24 figures. Submitted to Journal of Physics: Conference Serie

Large-eddy simulation and wall modelling of turbulent channel flow

We report large-eddy simulation (LES) of turbulent channel flow. This LES neither resolves nor partially resolves the near-wall region. Instead, we develop a special near-wall subgrid-scale (SGS) model based on wall-parallel filtering and wall-normal averaging of the streamwise momentum equation, with an assumption of local inner scaling used to reduce the unsteady term. This gives an ordinary differential equation (ODE) for the wall shear stress at every wall location that is coupled with the LES. An extended form of the stretched-vortex SGS model, which incorporates the production of near-wall Reynolds shear stress due to the winding of streamwise momentum by near-wall attached SGS vortices, then provides a log relation for the streamwise velocity at the top boundary of the near-wall averaged domain. This allows calculation of an instantaneous slip velocity that is then used as a ‘virtual-wall’ boundary condition for the LES. A Kármán-like constant is calculated dynamically as part of the LES. With this closure we perform LES of turbulent channel flow for Reynolds numbers Re_τ based on the friction velocity u_τ and the channel half-width δ in the range 2 × 10^3 to 2 × 10^7. Results, including SGS-extended longitudinal spectra, compare favourably with the direct numerical simulation (DNS) data of Hoyas & Jiménez (2006) at Re_τ = 2003 and maintain an O(1) grid dependence on Re_τ

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

An efficient Two-Layer wall model for accurate numerical simulations of aeronautical applications

Two-Layer wall models have been widely studied since they allow wall modeled Large Eddy Simulationsof general non-equilibrium flows. However, they are plagued by two persistent problems, the "log-layermismatch" and the resolved Reynolds stresses inflow. Several methodologies have been proposed so far todeal with these problems separately. In this work, a time-filtering methodology is used to tackle both issuesat once with a single and low-computational-cost step, easily applicable to complex three-dimensionalgeometries. Additionally, it is shown that the techniques intended to suppress the Reynolds stresses inflowproposed so far, were not sufficient to completely mitigate their detrimental effects.Peer ReviewedPostprint (published version

Some issues concerning Large-Eddy Simulation of inertial particle dispersion in turbulent bounded flows

The problem of an accurate Eulerian-Lagrangian modeling of inertial particle dispersion in Large Eddy Simulation (LES) of turbulent wall-bounded flows is addressed. We run Direct Numerical Simulation (DNS) for turbulent channel flow at shear Reynolds numbers equal to 150 and 300 and corresponding a-priori and a-posteriori LES on differently coarse grids. We then tracked swarms of different inertia particles and we examined the influence of filtering and of Sub-Grid Scale (SGS) modeling for the fluid phase on particle velocity and concentration statistics. We also focused on how particle preferential segregation is predicted by LES. Results show that even ``well-resolved'' LES is unable to reproduce the physics as demonstrated by DNS, both for particle accumulation at the wall and for particle preferential segregation. Inaccurate prediction is observed for the entire range of particles considered in this study, even when the particle response time is much larger than the flow timescales not resolved in LES. Both a-priori and a-posteriori tests indicate that recovering the level of fluid and particle velocity fluctuations is not enough to have accurate prediction of near-wall accumulation and local segregation. This may suggest that reintroducing the correct amount of higher-order moments of the velocity fluctuations is also a key point for SGS closure models for the particle equation. Another important issue is the presence of possible flow Reynolds number effects on particle dispersion. Our results show that, in small Reynolds number turbulence and in the case of heavy particles, the shear fluid velocity is a suitable scaling parameter to quantify these effects