Abstract
A strategy for a general electrochemical simulator is presented based on formulating the finite difference form of coupled mass transport and kinetic equations as a single sparse matrix. A treatment of a general mechanism consisting of first- and second-order homogeneous steps and quasi-reversible heterogeneous couples is developed for both steady-state and transient conditions. This is illustrated for two-dimensional simulations with potential applications including microdisk and channel microband electrodes, though the method may easily be applied to higher or lower dimensions and could include terms describing mass transport due to migration. Nonlinearities due to second-order kinetics are treated using Newton's method with an analytically derived Jacobian. ILU and MILU preconditioned Krylov subspace methods (CGS, BICGSTAB(𝓁 ), and RGMRES) are used to solve the resulting nonsymmetric sparse linear system. The efficiency of these solvers is compared to “brute force” LU factorization, SIP, and multigrid methods for a typical steady-state channel microband problem.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.