Explicitly enforce zero flux in absence of other boundary information.

This commit is contained in:
Bård Skaflestad
2010-11-25 14:33:04 +01:00
parent d9493e94ce
commit 110b78b84b

View File

@@ -510,6 +510,9 @@ compute_flux(grid_t *G,
} else { } else {
dp = cpress[c1] - bc->bcval[f]; dp = cpress[c1] - bc->bcval[f];
} }
} else {
/* No BC -> no-flow (== zero pressure drop) */
dp = 0.0;
} }
fflux[f] = t * dp; fflux[f] = t * dp;