Modeling transport of charged species in pore networks: Solution of the Nernst-Planck equations coupled with fluid flow and charge conservation equations.