A robust model and numerical approach for solving solid oxide fuel cell (SOFC) problems