Multiscale Numerical Modeling of Solute Transport with TwoPhase Flow in a Porous Cavity

This paper introduces dimensional and numerical investigation of the problem of solute transport within the two-phase flow in a porous cavity. The model consists of momentum equations (Darcy’s law), mass (saturation) equation, and solute transport equation. The cavity boundaries are constituted by mixed Dirichlet-Neumann boundary conditions. The governing equations have been converted into a dimensionless form such that a group of dimensionless physical numbers appear including Lewis, Reynolds, Bond, capillary, and Darcy numbers. A time-splitting multiscale scheme has been developed to treat the time derivative discretization. Also, we use the Courant-Friedrichs-Lewy (CFL) stability condition to adapt the time step size. The pressure is calculated implicitly by coupling Darcy’s law and the continuity equation, then, the concentration equation is solved implicitly. Numerical experiments have been conducted and the effects of the dimensionless numbers have been on the saturation, concentration, pressure, velocity, and Sherwood number have been investigated.