I am struggling in understanding how to implement the additions to the matrix beside the additions to the rhs.
It says deltaTime * J^T * Ms^-1 * J is added to the matrix but I don't quite understand what J^T and J are. I thought it would be fluid and solid fractions accordingly (divided by voxeldelta) but I am quite confused.
Is anybody able to help or understand how this equation expresses itself codewise (for example for the right face of a cell)?
Thanks in advance