coefficients = RANS3PF.Coefficients(
epsFact=epsFact_viscosity,
sigma=0.0,
rho_0=rho_0,
nu_0=nu_0,
rho_1=rho_1,
nu_1=nu_1,
g=g,
nd=nd,
ME_model=V_model,
PRESSURE_model=PRESSURE_model,
SED_model=SED_model,
VOS_model=VOS_model,
VOF_model=VOF_model,
LS_model=LS_model,
Closure_0_model=Closure_0_model,
Closure_1_model=Closure_1_model,
epsFact_density=epsFact_density,
stokes=False,
useVF=useVF,
useRBLES=useRBLES,
useMetrics=useMetrics,
eb_adjoint_sigma=1.0,
eb_penalty_constant=weak_bc_penalty_constant,
forceStrongDirichlet=ns_forceStrongDirichlet,
turbulenceClosureModel=ns_closure,
movingDomain=movingDomain,
dragAlpha=dragAlpha,
PSTAB=1.0,
cE=cE,
cMax=cMax,
CORRECT_VELOCITY=CORRECT_VELOCITY)
coefficients = RANS3PF.Coefficients(
epsFact=epsFact_viscosity,
sigma=0.0,
rho_0=rho_0,
nu_0=nu_0,
rho_1=rho_1,
nu_1=nu_1,
g=g,
nd=nd,
ME_model=V_model,
PRESSURE_model=PRESSURE_model,
SED_model=SED_model,
VOF_model=VOF_model,
VOS_model=VOS_model,
LS_model=LS_model,
Closure_0_model=Closure_0_model,
Closure_1_model=Closure_1_model,
epsFact_density=epsFact_density,
stokes=False,
useVF=useVF,
useRBLES=useRBLES,
useMetrics=useMetrics,
eb_adjoint_sigma=1.0,
eb_penalty_constant=weak_bc_penalty_constant,
forceStrongDirichlet=ns_forceStrongDirichlet,
turbulenceClosureModel=ns_closure,
movingDomain=movingDomain,
dragAlpha=dragAlpha,
PSTAB=ct.opts.PSTAB,
aDarcy=sedClosure.aDarcy,
betaForch=sedClosure.betaForch,
grain=sedClosure.grain,
packFraction=sedClosure.packFraction,
maxFraction=sedClosure.maxFraction,
frFraction=sedClosure.frFraction,
sigmaC=sedClosure.sigmaC,
C3e=sedClosure.C3e,
C4e=sedClosure.C4e,
eR=sedClosure.eR,
fContact=sedClosure.fContact,
mContact=sedClosure.mContact,
nContact=sedClosure.nContact,
angFriction=sedClosure.angFriction,
vos_function=ct.vos_function,
vos_limiter=ct.sedClosure.vos_limiter,
mu_fr_limiter=ct.sedClosure.mu_fr_limiter,
)
coefficients = RANS3PF.Coefficients(
epsFact=epsFact_viscosity,
sigma=0.0,
rho_0=rho_0,
nu_0=nu_0,
rho_1=rho_1,
nu_1=nu_1,
g=g,
nd=nd,
ME_model=V_model,
PRESSURE_model=PRESSURE_model,
SED_model=SED_model,
VOF_model=VOF_model,
VOS_model=VOS_model,
LS_model=LS_model,
Closure_0_model=Closure_0_model,
Closure_1_model=Closure_1_model,
epsFact_density=epsFact_density,
stokes=False,
useVF=useVF,
useRBLES=useRBLES,
useMetrics=useMetrics,
eb_adjoint_sigma=1.0,
eb_penalty_constant=weak_bc_penalty_constant,
forceStrongDirichlet=ns_forceStrongDirichlet,
turbulenceClosureModel=ns_closure,
movingDomain=movingDomain,
dragAlpha=dragAlpha,
PSTAB=0.0,
nParticles=1,
particle_epsFact=1.5,
particle_alpha=1e6,
particle_beta=1e6,
particle_penalty_constant=100.0,
particle_sdfList=[particle_sdf],
particle_velocityList=[particle_vel],
use_ball_as_particle=1,
ball_center=numpy.array([[0.2, 0.2, 0.0]]),
ball_radius=numpy.array([0.05]),
ball_velocity=numpy.array([[0.0, 0.0, 0.0]]),
ball_angular_velocity=numpy.array([[0.0, 0.0, 0.0]]),
use_sbm=USE_SBM)