There seems to be a problem with well boundary conditions when simulating reservoir production processes by the Darcy’s law module. I assigned Neumanna boundary conditions (constant mass flow) to the injection and production wells, and then the calculation was perfomed smoothly. In post-processing, I used the Darcy’s velocity to solve the mass flux of the production well, but I found that the mass flux calculated by the integral was inconsistent with the pre-prescribed mass flux. So, when I simulate the production process using the well boundary conditions with the specified production pressure, is it suitble to calculate the productivity by the the Darcy's velocity? Can someone give me some advice?
Thanks!