refactor(THM): Optimize the computation of the contribution of strain coupling to the RHS
In this MR, the computation of the contribution of strain coupling to the RHS is changed from \text{K}_{up}(\hat p - \hat p_{\text{prev}})/\Delta t to directly evaluating the following integral:
\int \alpha_B (\epsilon_v - \epsilon_v^{\text{prev}})/\Delta t \text{N} \mathrm{d} \Omega
- Tests covering your feature were added? with existing tests.
- Any new feature or behaviour change was documented?
- Add this merge request to the changelog
Edited by wenqing