with Dbulk and δbulk stand for the bulk electrolyte diffusivity and the bulk ionic conductivity, respectively. The Dbulk was set at 1 m 2 s -1 and δbulk was provided by the supplier and it is equal to 1.25 mScm -1 . All of these input parameters were considered isotropic. Periodic boundary conditions were considered for the outer xz and yz planes.