def __init__(self):
from proteus.mprans.SedClosure import HsuSedStress
self.C4e = 1.
self.C3e = 1.2
self.eR = 0.8
self.aDarcy = 1.
self.bForch = 1.
self.grain = 0.1
self.packFraction = 0.2
self.packMargin = 0.01
self.sigmaC = 1.1
self.maxFraction = 0.635
self.frFraction = 0.57
self.fContact = 0.02
self.mContact = 2.
self.nContact = 5.
self.angFriction = np.pi / 6.
self.vos_limiter = 0.6
self.mu_fr_limiter = 0.1
self.sedSt = HsuSedStress(self.aDarcy, self.bForch, self.grain,
self.packFraction, self.packMargin,
self.maxFraction, self.frFraction,
self.sigmaC, self.C3e, self.C4e, self.eR,
self.fContact, self.mContact, self.nContact,
self.angFriction, self.vos_limiter,
self.mu_fr_limiter)
self.sedSt_nl = HsuSedStress(
self.aDarcy, self.bForch, self.grain, self.packFraction,
self.packMargin, self.maxFraction, self.frFraction, self.sigmaC,
self.C3e, self.C4e, self.eR, self.fContact, self.mContact,
self.nContact, self.angFriction, self.vos_limiter, 1e100)
Ut = opts.inflow_vel * sqrt(cf / 2.)
Yplus = Y_ * Ut / opts.nu_0
kappaP = (Ut**2) / sqrt(opts.Cmu)
dissipationP = (Ut**3) / (opts.K * Y_)
# ----- Sediment stress ----- #
sedClosure = HsuSedStress(
aDarcy=opts.alphaSed,
betaForch=opts.betaSed,
grain=opts.grain,
packFraction=opts.packFraction,
packMargin=opts.packMargin,
maxFraction=opts.maxFraction,
frFraction=opts.frFraction,
sigmaC=opts.sigmaC,
C3e=opts.C3e,
C4e=opts.C4e,
eR=opts.eR,
fContact=opts.fContact,
mContact=opts.mContact,
nContact=opts.nContact,
angFriction=opts.angFriction,
vos_limiter=opts.vos_limiter,
mu_fr_limiter=opts.mu_fr_limiter,
)
# ----- DOMAIN ----- #
domain = Domain.PlanarStraightLineGraphDomain()
# ----- Phisical constants ----- #
steady_current = wt.SteadyCurrent(U=[opts.inflow_vel, 0, 0],
mwl=opts.waterLevel,
rampTime=0.8)
# ----- Sediment stress ----- #
sedClosure = HsuSedStress(
aDarcy=150.0,
betaForch=0.0,
grain=0.0025,
packFraction=0.2,
packMargin=0.01,
maxFraction=0.635,
frFraction=0.57,
sigmaC=1.1,
C3e=1.2,
C4e=1.0,
eR=0.8,
fContact=0.05,
mContact=3.0,
nContact=5.0,
angFriction=pi / 6.,
vos_limiter=0.62,
mu_fr_limiter=1e-3,
)
# ----- DOMAIN ----- #
domain = Domain.PlanarStraightLineGraphDomain()
# ----- Phisical constants ----- #