Hello!

I have a question with regards to a simple surface integral on, for example, a boxmesh geometry. If I solve a Poisson equation and the calculate the flow through the surface the standard technique is as:

flow_side_surface = dot(-K*grad(up), n)*ds(1)

where up is the solution of the Poisson equation, K a coefficient and ds(1) is the surface on the boundary I wish to integrate over.

Now, when I have a source on that boundary, how does the integral deal with it? I am asking this question because, on a simple test boxmesh I test for mass conservation by dividing the flow provided from the source by the sum of flow through the surfaces of the mesh. The ratio should be approximately 1. However, if the source is on a surface the ratio becomes 1.333. If the source is at the intersection of 2 surfaces the ratio is 2. And if the source is at the intersection of 3 surfaces (corner node) the ratio is 4. There’s a clear pattern but I can’t quite figure out why this is. Does anyone have any experience with this?

UPDATE: Clarification on what source means here. I have a numerical value for a source and I know where the source is in geometric space, so I create a function space for f and change the vector value at the boundary point to some numerical value. When this is done on the boundaries I get the problems given above, but not when I change an interior node value.