Do we need Neumann boundary condition for the pore network model?
The pore network model does not support pore throats at the boundary and consequently does not support Neumann boundary condition. Instead, we use source terms at the boundary pore bodies, when we want to inject to network. Although having a source term at the boundary pore body should give the same results as a Neumann-type boundary condition (when we think about the computing of the residual...), sometimes it causes confusion for the users when they want to interpret the simulation results (e.g, when a fluid is injected to a boundary pore body using a source term, but the flow direction in the connected throat is also into the pore body). That leads to these questions:
- Do we need Neumann boundary condition for the pore network model? An implementation which gives the user an option to assign the desirable flux to the throats connected to the boundary pore bodies?
- Or even do we need to have an option to have pore throats at the boundary instead of pore bodies?
Edited by Maziar Veyskarami