Marching cubes: A high resolution 3D surface construction algorithm
SIGGRAPH '87 Proceedings of the 14th annual conference on Computer graphics and interactive techniques
A multi-relaxation lattice kinetic method for passive scalar diffusion
Journal of Computational Physics
Advection-diffusion lattice Boltzmann scheme for hierarchical grids
Computers & Mathematics with Applications
Hi-index | 31.45 |
We describe an algorithm for simulating reactive flows in porous media, in which the pore space is mapped explicitly. Chemical reactions at the solid-fluid boundaries lead to dissolution (or precipitation), which makes it necessary to track the movement of the solid-fluid interface during the course of the simulation. We have developed a robust algorithm for constructing a piecewise continuous (C"1) surface, which enables a rapid remapping of the surface to the grid lines. The key components of the physics are the Navier-Stokes equations for fluid flow in the pore space, the convection-diffusion equation to describe the transport of chemical species, and rate equations to model the chemical kinetics at the solid surfaces. A lattice-Boltzmann model was used to simulate fluid flow in the pore space, with linear interpolation at the solid boundaries. A finite-difference scheme for the concentration field was developed, taking derivatives along the direction of the local fluid velocity. When the flow is not aligned with the grid this leads to much more accurate convective fluxes and surface concentrations than a standard Cartesian template. A robust algorithm for the surface reaction rates has been implemented, avoiding instabilities when the surface is close to a grid point. We report numerical tests of different aspects of the algorithm and assess the overall convergence of the method.