consThermERel :: Relation
consThermERel = negate (sy gradient) $. sy thFluxVect + sy volHtGen $=
  sy density * sy heatCapSpec * pderiv (sy temp) time
consThermERel :: Relation
consThermERel = negate (sy gradient) $. sy thFluxVect + sy volHtGen $=
  sy density * sy heatCapSpec * pderiv (sy temp) time
consThermERel = negate (sy gradient) $. sy thFluxVect + sy volHtGen $=
  sy density * sy heatCapSpec * pderiv (sy temp) time