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 ----- #