Karthik R

Couple of thoughts on how to debug:

  • Are you printing out the values of the source term and are they consistent with your hand calculations (I'm assuming that they are not significantly different). I'm also curious about these values just before the solution diverges. Do some of these values become really small or larger?

  • What happens if you slowly ramp-up the source term over time?

  • Just before the solution diverges, could you please try and plot the values of all physical parameters that feature in your source term and identify if there is an anomaly somewhere?