Fully resolved scalar transport for high Prandtl number flows using adaptive mesh refinement

Abstract In multiphase systems, boundary layers occur at fluid-fluid or fluid-solid interfaces. In a direct numerical simulation, the grid requirements are often dictated by the thickness of these boundary layers. Systems that are characterized by high Prandtl (or Schmidt number) exhibit temperature (or mass) boundary layers that are much thinner than the momentum boundary layers. In this paper, a hybrid computational approach is presented that uses a fixed Cartesian grid for the Navier-Stokes and continuity equations and an adaptive mesh for scalar transport, thus reducing the memory and CPU requirements tremendously while resolving all boundary layers. We describe the key aspects that need to be addressed in this hybrid approach, related to discretization, grid mapping, velocity interpolation along with detailed verification tests. Finally, the robustness and accuracy of our hybrid methodology is demonstrated for forced-convection heat transfer over stationary spherical particles at high Prandtl numbers.

[1]  Ng Niels Deen,et al.  A DNS study of flow and heat transfer through slender fixed-bed reactors randomly packed with spherical particles , 2017 .

[2]  Gretar Tryggvason,et al.  Multiscale computations of mass transfer from buoyant bubbles , 2012 .

[3]  Xiaoming Cai,et al.  Coupling dynamics and chemistry in the air pollution modelling of street canyons: A review. , 2016, Environmental pollution.

[4]  Ng Niels Deen,et al.  Immersed boundary method (IBM) based direct numerical simulation of open-cell solid foams , 2017 .

[5]  A. W. Vreman.  Particle-resolved direct numerical simulation of homogeneous isotropic turbulence modified by small fixed spheres , 2016, Journal of Fluid Mechanics.

[6]  F. Recasens,et al.  CFD study on particle-to-fluid heat transfer in fixed bed reactors: Convective heat transfer at low and high pressure , 2006 .

[7]  J. Kuipers,et al.  Direct numerical simulation of particulate flow with heat transfer , 2013 .

[8]  Ng Niels Deen,et al.  Direct numerical simulation of flow and heat transfer in dense fluid-particle systems , 2012 .

[9]  W. E. Ranz,et al.  Evaporation from drops , 1952 .

[10]  Ng Niels Deen,et al.  Immersed boundary method (IBM) based direct numerical simulation of open-cell solid foams:hydrodynamics , 2017 .

[11]  Georg Stadler,et al.  The Dynamics of Plate Tectonics and Mantle Flow: From Local to Global Scales , 2010, Science.

[12]  Rodolfo Ostilla Monico,et al.  A multiple-resolution strategy for Direct Numerical Simulation of scalar turbulence , 2015, J. Comput. Phys..

[13]  Carsten Burstedde,et al.  p4est: Scalable Algorithms for Parallel Adaptive Mesh Refinement on Forests of Octrees , 2011, SIAM J. Sci. Comput..

[14]  Huijin Xu,et al.  Review on heat conduction, heat convection, thermal radiation and phase change heat transfer of nanofluids in porous media: Fundamentals and applications , 2019, Chemical Engineering Science.

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

[16]  S. Popinet.  Gerris: a tree-based adaptive solver for the incompressible Euler equations in complex geometries , 2003 .

[17]  Li-Tien Cheng,et al.  A second-order-accurate symmetric discretization of the Poisson equation on irregular domains , 2002 .

[18]  Carsten Burstedde,et al.  Recursive Algorithms for Distributed Forests of Octrees , 2014, SIAM J. Sci. Comput..

[19]  R. Chhabra,et al.  Effect of blockage on drag and heat transfer from a single sphere and an in-line array of three spheres , 2006 .

[20]  J. Centrella,et al.  Planar numerical cosmology. II. The difference equations and numerical tests , 1984 .

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

[22]  R. Andricevic,et al.  Contaminant dilution measure for the solute transport in an estuary , 2018, Advances in Water Resources.

[23]  Vinesh H. Gada,et al.  On a Novel Dual-Grid Level-Set Method for Two-Phase Flow Simulation , 2011 .

[24]  S. Whitaker.  Forced convection heat transfer correlations for flow in pipes, past flat plates, single cylinders, single spheres, and for flow in packed beds and tube bundles , 1972 .

[25]  A. W. Vreman,et al.  Turbulent channel flow past a moving array of spheres , 2018, Journal of Fluid Mechanics.

[26]  Anthony G. Dixon,et al.  Packed Tubular Reactor Modeling and Catalyst Design using Computational Fluid Dynamics , 2006 .

[27]  Mikhail Shashkov,et al.  Mimetic finite difference methods for diffusion equations on non-orthogonal non-conformal meshes , 2004 .

[28]  Kai Leong Chong,et al.  Multiple-resolution scheme in finite-volume code for active or passive scalar turbulence , 2018, J. Comput. Phys..

[29]  Dieter Bothe,et al.  Advanced subgrid-scale modeling for convection-dominated species transport at fluid interfaces with application to mass transfer from rising bubbles , 2017, J. Comput. Phys..

[30]  Robert D. Falgout,et al.  hypre: A Library of High Performance Preconditioners , 2002, International Conference on Computational Science.

[31]  S. Osher,et al.  Spatially adaptive techniques for level set methods and incompressible flow , 2006 .