Hello,

I am trying to write out the f0p residual kernel for poroelasticity in a manner that is usable for both quasi-static and dynamic formulations. As it stands, the QS formulation solves for three variables displacement, pressure, and volumetric strain (u, p, \epsilon_{v}), while the dynamic solves for (u, p, v), where v stands for velocity.

The f0p residual kernel is the one for time variation in fluid content, and as such contains a term that includes \frac{\partial \epsilon_{v}}{\partial t}. I would prefer to solve this by dividing the trace of the strain tensor by the timestep length (dt), but am unsure as how to access the value for dt within a kernel callback function, or even if it is possible.

Please advise,

Robert