Low-diffusivity scalar transport using a WENO scheme and dual meshing

JOURNAL OF COMPUTATIONAL PHYSICS

卷:240

Article

Kubrak, B.¹  Herlina, H.²  Greve, F.^3,4  Wissink, J. G.¹ 

[1] Brunel Univ London, Uxbridge UB8 3PH, Middx, England

[2] Karlsruhe Inst Technol, D-76021 Karlsruhe, Germany

[3] Tech Univ Dresden, D-01069 Dresden, Germany

[4] KIT Karlsruhe, D-76131 Karlsruhe, Germany

关键词: Air-water interface;    DNS;    Gas transfer;    WENO scheme;    Scalar transport;    High Schmidt number;   

DOI  :  10.1016/j.jcp.2012.12.039

来源: Elsevier

【 摘 要 】

Interfacial mass transfer of low-diffusive substances in an unsteady flow environment is marked by a very thin boundary layer at the interface and other regions with steep concentration gradients. A numerical scheme capable of resolving accurately most details of this process is presented. In this scheme, the fifth-order accurate WENO method developed by [13] was implemented on a non-uniform staggered mesh to discretize the scalar convection while for the scalar diffusion a fourth-order accurate central discretization was employed. The discretization of the scalar convection-diffusion equation was combined with a fourth-order Navier-Stokes solver which solves the incompressible flow. A dual meshing strategy was employed, in which the scalar was solved on a finer mesh than the incompressible flow. The order of accuracy of the solver for one-dimensional scalar transport was tested on both stretched and uniform grids. Compared to the fifth-order WENO implementation of [10], the [13] method was found to be superior on very coarse meshes. The solver was further tested by performing a number of two-dimensional simulations. At first a grid refinement test was performed at zero viscosity with shear acting on an initially axisymmetric scalar distribution. A second refinement test was conducted for an unstably stratified flow with low diffusivity scalar transport. The unstable stratification led to buoyant convection which was modelled using a Boussinesq approximation with a linear relationship between flow temperature and density. The results show that for the method presented a relatively coarse mesh is sufficient to accurately describe the fluid flow, while the use of a refined dual mesh for the low-diffusive scalars is found to be beneficial in order to obtain a highly accurate resolution with negligible numerical diffusion. (C) 2013 Elsevier Inc. All rights reserved.

【 授权许可】

Free