Numerical study of isothermal drying in porous media

This work is focus on the efficient numerical computation of solution of coupled quasilinear partial differential equations of parabolic type. These equations model the isothermal drying of porous media. The complete system consists of two primary variables along with 19 dependent variables which are solution dependent. The governing equations are two strongly coupled quasilinear partial differential equations and 19 nonlinear algebraic equations of state. First, the physical behavior of the variables through the numerical simulation in the one-dimensional case using finite volume method in space and implicit Euler in time is presented. A local change of drying states causes local steep gradients which restrict the time step size in those local regions. To increase the efficiency, next a numerical scheme with adaptive time stepping using linearly implicit methods is implemented. The positivity preservation and the stability of the solution is discussed under some assumptions. Numerical results in the two-dimensional case were established using dimensional splitting. © 2021 by Begell House, Inc. www.begellhouse.com