A multiple-resolution strategy for Direct Numerical Simulation of scalar turbulence

In this paper a numerical procedure to simulate low diffusivity scalar turbulence is presented. The method consists of using a grid for the advected scalar with a higher spatial resolution than that of the momentum. The latter usually requires a less refined mesh and integrating both fields on a single grid tailored to the most demanding variable produces an unnecessary computational overhead. A multiple resolution approach is used also in the time integration in order to maintain the stability of the scalars on the finer grid. The method is the more advantageous the less diffusive the scalar is with respect to momentum, therefore it is particularly well suited for large Prandtl or Schmidt number flows. However, even in the case of equal diffusivities the present procedure gives CPU time and memory occupation savings, due to the increased gradients and more intermittent behaviour of the scalars when compared to momentum.

[1]  Katepalli R. Sreenivasan,et al.  Under consideration for publication in J. Fluid Mech. 1 Very fine structures in scalar mixing , 2008 .

[2]  P. Linden.  THE FLUID MECHANICS OF NATURAL VENTILATION , 1999 .

[3]  J. S. Andersen,et al.  Passive scalar intermittency in low temperature helium flows. , 2000, Physical review letters.

[4]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[5]  Georges-Henri Cottet,et al.  Hybrid spectral-particle method for the turbulent transport of a passive scalar , 2014, J. Comput. Phys..

[6]  Siddhartha Verma,et al.  On filtering in the viscous-convective subrange for turbulent mixing of high Schmidt number passive scalars , 2013 .

[7]  Vergassola,et al.  Universality and saturation of intermittency in passive scalar turbulence , 2000, Physical review letters.

[8]  J. Burgers.  A mathematical model illustrating the theory of turbulence , 1948 .

[9]  Diego Rossinelli,et al.  GPU and APU computations of Finite Time Lyapunov Exponent fields , 2012, J. Comput. Phys..

[10]  Diego Donzis,et al.  Resolution effects and scaling in numerical simulations of passive scalar mixing in turbulence , 2010 .

[11]  J. Andrzej Domaradzki,et al.  Direct numerical simulations of passive scalars with Pr>1 advected by turbulent flow , 1997, Journal of Fluid Mechanics.

[12]  Detlef Lohse,et al.  Heat transfer and large scale dynamics in turbulent Rayleigh-Bénard convection , 2008, 0811.0471.

[13]  Toshiyuki Gotoh,et al.  Spectral compact difference hybrid computation of passive scalar in isotropic turbulence , 2012, J. Comput. Phys..

[14]  Detlef Lohse,et al.  Prandtl and Rayleigh number dependence of heat transport in high Rayleigh number thermal convection , 2011, Journal of Fluid Mechanics.

[15]  F Toschi,et al.  Optimal interpolation schemes for particle tracking in turbulence. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  Shraiman,et al.  Heat transport in high-Rayleigh-number convection. , 1990, Physical review. A, Atomic, molecular, and optical physics.

[17]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[18]  Detlef Lohse,et al.  Plume emission statistics in turbulent Rayleigh–Bénard convection , 2014, Journal of Fluid Mechanics.

[19]  D. Lohse,et al.  Small-Scale Properties of Turbulent Rayleigh-Bénard Convection , 2010 .

[20]  Roberto Verzicco,et al.  A pencil distributed finite difference code for strongly turbulent wall-bounded flows , 2015, 1501.01247.

[21]  Katepalli R. Sreenivasan,et al.  On local isotropy of passive scalars in turbulent shear flows , 1991, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[22]  Eric D. Siggia,et al.  Scalar turbulence , 2000, Nature.

[23]  Alain Pumir,et al.  Structure of the three-point correlation function of a passive scalar in the presence of a mean gradient , 1998 .

[24]  S. Arya.  Air Pollution Meteorology and Dispersion , 1998 .

[25]  Alain Pumir,et al.  A numerical study of the mixing of a passive scalar in three dimensions in the presence of a mean gradient , 1994 .

[26]  W. Peltier,et al.  MIXING EFFICIENCY IN STRATIFIED SHEAR FLOWS , 2003 .

[27]  Chao Sun,et al.  Salinity transfer in bounded double diffusive convection , 2015, Journal of Fluid Mechanics.

[28]  R. Verzicco,et al.  A Finite-Difference Scheme for Three-Dimensional Incompressible Flows in Cylindrical Coordinates , 1996 .

[29]  Charles R. Doering,et al.  On upper bounds for infinite Prandtl number convection with or without rotation , 2001 .

[30]  Richard J. A. M. Stevens,et al.  Boundary layer structure in turbulent thermal convection and its consequences for the required numerical resolution , 2010, 1109.6870.

[31]  Z. Warhaft.  Passive Scalars in Turbulent Flows , 2000 .

[32]  Eric D. Siggia,et al.  High Rayleigh Number Convection , 1994 .

[33]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[34]  Frank P. Incropera,et al.  Convection heat transfer in electronic equipment cooling , 1988 .

[35]  Raymond W. Schmitt,et al.  DOUBLE DIFFUSION IN OCEANOGRAPHY , 1994 .

[36]  Detlef Lohse,et al.  Radial boundary layer structure and Nusselt number in Rayleigh–Bénard convection , 2009, Journal of Fluid Mechanics.