Preparatory work for fractures + wells

[paveltomin]

For now, assuming wells are not directly connected to the fractures.

In this version, only SinglePhasePoromechanicsConformingFracturesReservoir solver can run, although I am having hard time to build a well-behaving example. Other versions of solvers with wells are drafted but commented out for now.

Added poromechanics initialization for contact solver - technically runs but needs refinements.

Side thing: add default template types to avoid too long nested definitions

[jhuang2601]

Thanks for enabling fault initialization

[Guotong-Ren]

I am trying to understand the logic here. what is proposed is always using BASE besides singlePhaseProppantBase specialization? And I guess this implementation is assumed not only being restricted to SinglePhaseBase.

[paveltomin]

If I remember correctly, the change should be just related to the default template argument (SinglePhaseBase), no changes in other logic. SinglePhaseFVM< SinglePhaseProppantBase >::assembleFluxTerms is still down below.

[Guotong-Ren]

this question is not quite relevant to this PR. I am curious how much work do you estimate to prototype when well is connected to fractures? need adding extra blocks in the Jacobin matrix?

[paveltomin]

this question is not quite relevant to this PR. I am curious how much work do you estimate to prototype when well is connected to fractures? need adding extra blocks in the Jacobin matrix?

I would redirect this to @CusiniM for better estimate.
The code pieces should be there, but I am not exactly sure how to connect them.

[paveltomin]

this question is not quite relevant to this PR. I am curious how much work do you estimate to prototype when well is connected to fractures? need adding extra blocks in the Jacobin matrix?

I would redirect this to @CusiniM for better estimate. The code pieces should be there, but I am not exactly sure how to connect them.

I mean, I am ready to work on that, just need some guidance where to start.

[CusiniM]

this question is not quite relevant to this PR. I am curious how much work do you estimate to prototype when well is connected to fractures? need adding extra blocks in the Jacobin matrix?

I would redirect this to @CusiniM for better estimate. The code pieces should be there, but I am not exactly sure how to connect them.

I mean, I am ready to work on that, just need some guidance where to start.

I don't think it's a lot of work. However, it would be nice to first refactor the way wells are handled by introducing a WellStencil (or something similar) and then use a dedicated stencil for those connections. Once the connections are well defined adding the terms to the jacobian is trivial IMO.

[CusiniM]

is this necessary? Isn't it already being instantiated in SinglePhasePoromechanicsConformingFractures.cpp ?

[paveltomin]

I remember that somehow it was important but I can double-check
Would it miss versions of SinglePhasePoromechanics functions when only present in SinglePhasePoromechanicsConformingFractures?

[CusiniM]

ah yes, sounds possible