TY - JOUR

T1 - SPH-ϵ simulation of 2D turbulence driven by a moving cylinder

AU - Monaghan, J. J.

PY - 2017/9/1

Y1 - 2017/9/1

N2 - The aim of this paper is to describe how effectively a turbulence model for the particle method SPH (Monaghan 2011 denoted by SPH-ϵ) simulates 2D turbulence driven by a moving cylinder inside a square box with no-slip walls. The model makes few assumptions, conserves linear and angular momentum, satisfies a discrete version of Kelvin's circulation theorem and is computationally efficient. Furthermore, the results from the model for decaying two dimensional turbulence are in good agreement with the experimental and computational results. The SPH turbulence model is based on a Lagrangian similar to that used for the Lagrangian averaged Navier–Stokes (LANS) turbulence model, but with a different smoothed velocity. The smoothed velocity preserves the shape of the spectrum of the unsmoothed velocity for short length scales, but reduces its magnitude by an amount which depends on a parameter ϵ, where 0≤ϵ≤1. When ϵ=0, there is no smoothing, and when ϵ=1 the smoothing is a maximum. We compare calculations with resolution given by the Reynolds length ℓR and those using SPH-ϵ with ϵ=0.8, and a coarser resolution 2ℓR. The second and fourth order velocity correlation functions, the energy spectrum and the mixing calculated using SPH-ϵ and the coarse resolution 2ℓR, are in good agreement with those using the resolution ℓR, but requires 8 times less work to integrate to the same time. The mixing induced by the turbulence is in good agreement with high resolution calculations.

AB - The aim of this paper is to describe how effectively a turbulence model for the particle method SPH (Monaghan 2011 denoted by SPH-ϵ) simulates 2D turbulence driven by a moving cylinder inside a square box with no-slip walls. The model makes few assumptions, conserves linear and angular momentum, satisfies a discrete version of Kelvin's circulation theorem and is computationally efficient. Furthermore, the results from the model for decaying two dimensional turbulence are in good agreement with the experimental and computational results. The SPH turbulence model is based on a Lagrangian similar to that used for the Lagrangian averaged Navier–Stokes (LANS) turbulence model, but with a different smoothed velocity. The smoothed velocity preserves the shape of the spectrum of the unsmoothed velocity for short length scales, but reduces its magnitude by an amount which depends on a parameter ϵ, where 0≤ϵ≤1. When ϵ=0, there is no smoothing, and when ϵ=1 the smoothing is a maximum. We compare calculations with resolution given by the Reynolds length ℓR and those using SPH-ϵ with ϵ=0.8, and a coarser resolution 2ℓR. The second and fourth order velocity correlation functions, the energy spectrum and the mixing calculated using SPH-ϵ and the coarse resolution 2ℓR, are in good agreement with those using the resolution ℓR, but requires 8 times less work to integrate to the same time. The mixing induced by the turbulence is in good agreement with high resolution calculations.

UR - http://www.scopus.com/inward/record.url?scp=85018692446&partnerID=8YFLogxK

U2 - 10.1016/j.euromechflu.2017.03.011

DO - 10.1016/j.euromechflu.2017.03.011

M3 - Article

AN - SCOPUS:85018692446

VL - 65

SP - 486

EP - 493

JO - European Journal of Mechanics, B/Fluids

JF - European Journal of Mechanics, B/Fluids

SN - 0997-7546

ER -