Fixed issue #2853 with a flag

[liyangrock]

This PR addressed issue #2853 by adding a flag to determine if the lagging update should be applied to the fracture stencil weights.

This was done specifically for the hydrofracture solver, as it was observed that issue #2853 was currently causing performance degradation only for this particular solver. Therefore, the flag was exclusively introduced to the hydrofracture solver.

For further discussions related to issue #2853, please see #2922 .

[Guotong-Ren]

To understand, is the flow area in fracture from last time step? Do you have any thoughts on evaluation of the time step size impact on the solution?

[liyangrock]

To understand, is the flow area in fracture from last time step? Do you have any thoughts on evaluation of the time step size impact on the solution?

Indeed, the calculation of fracture stencil weights is performed using the hydraulicAperture from the previous time step. I have conducted a comparison between dt=0.25s and dt=5.0s for viscosity-dominated KGD cases. The resulting net pressure and aperture at the mouth are displayed below.

The difference between the numerical result (with dt=0.25s) and the analytical result can be attributed to the fact that the numerical scenario employs a low toughness value (i.e., 1e4), whereas the analytical solution does not incorporate any toughness.

The input files are included here.
kgdVis.zip