Bug in AdTpfaFlux

[keileg]

There seems to be a bug in the AdTpfaFlux mixin class related to Neumann conditions on internal domains.

A likely partial cause is that, in AdTpfaFlux.diffusive_flux, filtering on Neumann conditions is only done on external boundaries (contrast this with AdTpfaFlux.pressure_trace, where internal boundaries is also filtered). It could also be that both methods need filtering on fracture tips, it seems likely this is currently not done.

To get this right, and properly tested, we need to solve a flow problem with both intersections and fracture tips.

@IvarStefansson: Feel free to add to the problem description.

[IvarStefansson]

Agree with the description. This setup might be relevant:

porepy/tests/models/test_thermoporomechanics.py

Line 325 in e4d484f

class ThermoporomechanicsWell(

If you use it, could you be bothered to fix #1108?

[IvarStefansson]

Or this one, depending on whether you want to include Fourier or not:

porepy/tests/models/test_poromechanics.py

Line 513 in e4d484f

class PoromechanicsWell(

[keileg]

Or this one, depending on whether you want to include Fourier or not:

porepy/tests/models/test_poromechanics.py

Line 513 in e4d484f

class PoromechanicsWell(

Will do.