GeoEnergyLab-EPFL / PyFrac

a Python planar 3D hydraulic fracture growth solver
http://pyfrac.epfl.ch/
Other
34 stars 20 forks source link

Fracture closure #15

Closed Skonki closed 9 months ago

Skonki commented 11 months ago

Fracture closure in the case of re-meshing and/or the presence of fluid compressibility seems to have a bug.

Skonki commented 11 months ago

Edit: The issue is also there if we do not have re-meshing and no fluid compressibility.

This means that between the simulations run for the recession paper and this update there seems to be an issue (as we do not use iterative solve here).

Skonki commented 11 months ago

Figure_ClosureIssue2 Figure_ClosureIssue

See how the opening goes up shortly before closure. Oddly, the efficiency does not.

I will now investigate volume.

Note that this is without re-meshing..the result in case of re-meshing is identical.

Skonki commented 9 months ago

We observe similar issues with different (order one results) parameter sets.

However, this might be a known problem of leak-off continuing in cells with activated width constraint.

I add the task to check this.

See figure of the results

Bildschirmfoto 2024-01-12 um 17 24 25
Skonki commented 9 months ago

For the data of Mohsen (used in the recession paper), this was actually the bug, and the results now show the expected behavior of the opening.

The test on the example (with re-meshing and all, rather coarse) is now running; let's see.

Bildschirmfoto 2024-01-16 um 07 07 22

Skonki commented 9 months ago

The problem is solved, pull request follows.

See here the new graphics as a result.

Bildschirmfoto 2024-01-16 um 11 21 56