def mixedBC(self, boundary, userData):
"""Apply mixed boundary conditions."""
if boundary.marker() != pg.MARKER_BOUND_MIXED:
return 0
sourcePos = pg.center(userData['sourcePos'])
k = userData['k']
r1 = boundary.center() - sourcePos
# Mirror on surface at depth=0
r2 = boundary.center() - pg.RVector3(1.0, -1.0, 1.0) * sourcePos
r1A = r1.abs()
r2A = r2.abs()
rho = 1.
if self.resistivity is not None:
rho = self.resistivity[boundary.leftCell().id()]
n = boundary.norm()
if r1A > 1e-12 and r2A > 1e-12:
if (pg.besselK0(r1A * k) + pg.besselK0(r2A * k)) > 1e-12:
return 1./rho * k * (r1.dot(n) / r1A * pg.besselK1(r1A * k) +
r2.dot(n) / r2A * pg.besselK1(r2A * k)) /\
(pg.besselK0(r1A * k) + pg.besselK0(r2A * k))
else:
return 0.
else:
return 0.
def mixedBC(self, boundary, userData):
"""Apply mixed boundary conditions."""
if boundary.marker() != pg.MARKER_BOUND_MIXED:
return 0
sourcePos = pg.center(userData['sourcePos'])
k = userData['k']
r1 = boundary.center() - sourcePos
# Mirror on surface at depth=0
r2 = boundary.center() - pg.RVector3(1.0, -1.0, 1.0) * sourcePos
r1A = r1.abs()
r2A = r2.abs()
rho = 1.
if self.resistivity is not None:
rho = self.resistivity[boundary.leftCell().id()]
n = boundary.norm()
if r1A > 1e-12 and r2A > 1e-12:
if (pg.besselK0(r1A * k) + pg.besselK0(r2A * k)) > 1e-12:
return 1./rho * k * (r1.dot(n) / r1A * pg.besselK1(r1A * k) +
r2.dot(n) / r2A * pg.besselK1(r2A * k)) /\
(pg.besselK0(r1A * k) + pg.besselK0(r2A * k))
else:
return 0.
else:
return 0.
def mixedBC(boundary, userData):
sourcePos = userData['sourcePos']
k = userData['k']
r1 = boundary.center() - sourcePos
# Mirror on surface at depth=0
r2 = boundary.center() - pg.RVector3(1.0, -1.0, 1.0) * sourcePos
r1A = r1.abs()
r2A = r2.abs()
n = boundary.norm()
# need rho here !!!!!!!!!!!!!!!!!!!!!!!!!!!1
if r1A > 1e-12 and r2A > 1e-12:
return k * ((r1.dot(n)) / r1A * pg.besselK1(r1A * k) +
(r2.dot(n)) / r2A * pg.besselK1(r2A * k)) / \
(pg.besselK0(r1A * k) + pg.besselK0(r2A * k))
else:
return 0.
def mixedBC(boundary, userData):
sourcePos = userData['sourcePos']
k = userData['k']
r1 = boundary.center() - sourcePos
# Mirror on surface at depth=0
r2 = boundary.center() - pg.RVector3(1.0, -1.0, 1.0) * sourcePos
r1A = r1.abs()
r2A = r2.abs()
n = boundary.norm()
# need rho here !!!!!!!!!!!!!!!!!!!!!!!!!!!1
if r1A > 1e-12 and r2A > 1e-12:
return k * ((r1.dot(n)) / r1A * pg.besselK1(r1A * k) +
(r2.dot(n)) / r2A * pg.besselK1(r2A * k)) / \
(pg.besselK0(r1A * k) + pg.besselK0(r2A * k))
else:
return 0.
def mixedBC(boundary, userData):
"""Mixed boundary conditions.
Define the derivative of the analytical solution regarding the outer normal
direction :math:`\vec{n}`. So we can define the value for the Neumann type
Boundary conditions for the boundaries in the subsurface.
"""
sourcePos = userData['sourcePos']
k = userData['k']
r1 = boundary.center() - sourcePos
# Mirror on surface at depth=0
r2 = boundary.center() - pg.RVector3(1.0, -1.0, 1.0) * sourcePos
r1A = r1.abs()
r2A = r2.abs()
n = boundary.norm()
if r1A > 1e-12 and r2A > 1e-12:
return k * ((r1.dot(n)) / r1A * pg.besselK1(r1A * k) +
(r2.dot(n)) / r2A * pg.besselK1(r2A * k)) / \
(pg.besselK0(r1A * k) + pg.besselK0(r2A * k))
else:
return 0.
def mixedBC(boundary, userData):
"""Mixed boundary conditions.
Define the derivative of the analytical solution regarding the outer normal
direction :math:`\vec{n}`. So we can define the value for the Neumann type
Boundary conditions for the boundaries in the subsurface.
"""
sourcePos = userData['sourcePos']
k = userData['k']
r1 = boundary.center() - sourcePos
# Mirror on surface at depth=0
r2 = boundary.center() - pg.RVector3(1.0, -1.0, 1.0) * sourcePos
r1A = r1.abs()
r2A = r2.abs()
n = boundary.norm()
if r1A > 1e-12 and r2A > 1e-12:
return k * ((r1.dot(n)) / r1A * pg.besselK1(r1A * k) +
(r2.dot(n)) / r2A * pg.besselK1(r2A * k)) / \
(pg.besselK0(r1A * k) + pg.besselK0(r2A * k))
else:
return 0.