approximate porosity change by forward Euler method
use central differences for porosity change