where Lc is the characteristic length and j is the current density. Sensitivity analysis of an electrodeoxidation model for titanium production requires the automatic generation of a finite element mesh using the Delaunay triangulation algorithm. The finite element method uses a weak formulation of field variable constraints in coupled systems of partial differential equations to make discontinuities integrable. This is important for accurate numerical simulation at corner discontinuities in a model domain. In addition, increasing the number of finite elements using mesh refinement techniques improve the accuracy to an acceptable level (i.e. error ≤ 5 % for simulations in 2D axisymmetric model domains).