Method of numerical investigation of saline groundwater dynamics

A numerical method for solving unsteady problems on two-component fluid dynamics in a porous medium that model transport of salt, dissolved in ground water, is developed. A mathematical model consists of the continuity equation, the Darcy’s low with the gravity, the salt diffusivity equation and the linear equation of state relating the solution density to the salt density. The numerical code is based on the finite-difference method realized on a straggled nonuniform grid. The continuity and Darcy’s equations are integrated with the use of the SIMPLE-type algorithm, then the salt transport equation is solved. The designed code is tested by solving three benchmark problems including diffusion of salt, transport of salt by the upflow and concentration convection.