Numerical simulation of an electrochemical system and semi-analytical method

Electrochemical systems are often simulated by using numerical methods based on finite element solution of differential equations. We developed a method to decrease the needed computational resources, based on the analytical solution of a part of the system: the obtained analytical equations are applied as boundary conditions to the finite element calculation. The part of the system that is analytically solved is the region of the diffuse double layer. We provide two COMSOL models: i) a 1d fully numerical calculation and ii) the same simulation performed with the semi-analytical method. For each model, we provide the .mph file (Comsol) and the .java file (Java+Comsol). We also provide examples of the resulting impedances and frequency responses of various parameters.