Hi!
I needed to do a restart with ThermoRichardsMechanics model.
I took the output of the last time step of the first run as bulk mesh and used the corresponding MeshNode field as initial values for temperature, pressure and displacement. I used the time stamp of the last time step as t0. compensate_non_equilibrium_initial_residuum was set to false.
While the temperature looks as expected, the pressure makes a jump after restart. (It seems like a jump in the pressure derivative as well).
Is there anything I missed?
FE3d-test-press-pt1.pdf (11.7 KB)
FE3d-test-temp-pt1.pdf (9.4 KB)