Bug in calculation resistance factor resfac in LeakyLineDoublet
#71
Labels
Comments
|
Other way to put it: The resistance of the leaky layer is really res / kaq. That is not desired, of course and should be fixed. |
mbakker7
pushed a commit
that referenced
this issue
Jul 5, 2023
* Add LeakyBuildingPit element - add LeakyIntHeadDiffEquation mixin class - add LeakyIntHeadDiffLineSink element for leaky pit boundaries with leaky walls - add LeakyBuildingPit element, inherits from BuildingPit but adds resistance term - add notebook testing LeakyBuildingPit allow BuildingPit to be placed in models with 1 layer isort __init__ * use self.disveclayers for RHS * fix #71 - resfac fixed in LineDoublet(1d) - resfac fixed in LineSink(1d) * update leaky building pit nb * allow passing kwargs to contour func * add integration of normal flux along line - add disvecnorm: calculate disvec normal to angle theta - intdisvecnorm: quad/legendre integration of disvecnorm - add quad integrand function for integrating from 0->r * rename *disvecnorm to *normflux - improve code and docstring * fix overwrite of label kwarg * add intnormflux notebook * add intnormflux notebook * add quad_vec import * improve normal flux nb * optimize code a bit and reorder coordinates * update normalflux nb * docstrings normflux functions * add circular leaky building pit notebook * rename to IntLeakyWallEquation * rename to IntLeakyWallEquation * delete nb * remove py3.7 testing, add 3.10 and 3.11 * fix buildingpit nb * add shapely to ci reqs
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
The resistance factor in LeakyLineDoublet is calculated as T / res, when this should be aquifer thickness / res. The value of the flux is therefore off by factor
kaq.The text was updated successfully, but these errors were encountered: