Configuration for solver: fully coupled Newton solver, with direct linear (MUMPS for FEDM, and PARDISO and MUMPS for COMSOL Multiphysics®) and iterative (gmres with hypre preconditioner) solver; nonlinear solver tolerance: 0.0001; time stepping tolerance: varied for different case studies