def create_scheme(self):
aha = AdamiHuAdamsScheme(['fluid'], ['dam'],
dim=2,
rho0=ro,
c0=co,
alpha=alp,
gy=-9.81,
nu=0.0,
h0=0.03,
gamma=1.0)
wcsph = WCSPHScheme(['fluid'], ['dam'],
dim=2,
rho0=ro,
c0=co,
h0=0.03,
hdx=1.3,
hg_correction=True,
gy=-9.81,
alpha=alp,
gamma=gamma,
update_h=True)
edac = EDACScheme(['fluid'], ['dam'],
dim=2,
rho0=ro,
c0=co,
gy=-9.81,
alpha=0.0,
nu=0.0,
h=0.03,
clamp_p=True)
return SchemeChooser(default='wcsph', wcsph=wcsph, aha=aha, edac=edac)
def create_scheme(self):
wcsph = WCSPHScheme(['fluid'], ['wall', 'paddle', 'obstacle'],
dim=2,
rho0=ro,
c0=co,
h0=0.01,
hdx=1.3,
gy=-9.81,
hg_correction=True,
alpha=alp,
gamma=gamma,
update_h=True)
edac = EDACScheme(['fluid'], ['wall', 'paddle', 'obstacle'],
dim=2,
rho0=ro,
c0=co,
gy=-9.81,
alpha=alp,
nu=0.0,
h=0.01,
clamp_p=True)
aha = AdamiHuAdamsScheme(['fluid'], ['wall', 'paddle', 'obstacle'],
dim=2,
rho0=ro,
h0=0.01,
gamma=1.0,
alpha=alp,
gy=-9.81,
nu=0.0,
c0=co)
return SchemeChooser(default='wcsph', aha=aha, wcsph=wcsph, edac=edac)
def create_scheme(self):
wcsph = WCSPHScheme(
['fluid'], ['boundary'], dim=2, rho0=ro, c0=co,
h0=h, hdx=1.3, gy=-9.81, alpha=alpha, beta=beta,
gamma=gamma, hg_correction=True, update_h=True
)
aha = AdamiHuAdamsScheme(
fluids=['fluid'], solids=['boundary'], dim=2, c0=co, nu=nu,
rho0=ro, h0=h, p0=0.0, gy=-g, gamma=1.0, tdamp=0.0, alpha=alpha
)
edac = EDACScheme(
fluids=['fluid'], solids=['boundary'], dim=2, c0=co, nu=nu,
rho0=ro, h=h, pb=0.0, gy=-g, eps=0.0, clamp_p=True
)
iisph = IISPHScheme(
fluids=['fluid'], solids=['boundary'], dim=2, nu=nu,
rho0=ro, gy=-g
)
gtvf = GTVFScheme(
fluids=['fluid'], solids=['boundary'], dim=2, nu=nu,
rho0=ro, gy=-g, h0=None, c0=co, pref=None
)
sisph = SISPHScheme(
fluids=['fluid'], solids=['boundary'], dim=2, nu=nu,
c0=co, rho0=ro, alpha=0.05, gy=-g, pref=ro*co**2,
internal_flow=False, hg_correction=True, gtvf=True, symmetric=True
)
s = SchemeChooser(default='wcsph', wcsph=wcsph, aha=aha, edac=edac,
iisph=iisph, gtvf=gtvf, sisph=sisph)
return s
def create_scheme(self):
wcsph = WCSPHScheme(['fluid'], ['boundary'],
dim=2,
rho0=ro,
c0=co,
h0=h,
hdx=1.3,
gy=-9.81,
alpha=alpha,
beta=beta,
gamma=gamma,
hg_correction=True,
update_h=True)
aha = AdamiHuAdamsScheme(fluids=['fluid'],
solids=['boundary'],
dim=2,
c0=co,
nu=nu,
rho0=ro,
h0=h,
p0=0.0,
gy=-g,
gamma=1.0,
tdamp=0.0,
alpha=alpha)
edac = EDACScheme(fluids=['fluid'],
solids=['boundary'],
dim=2,
c0=co,
nu=nu,
rho0=ro,
h=h,
pb=0.0,
gy=-g,
eps=0.0,
clamp_p=True)
iisph = IISPHScheme(fluids=['fluid'],
solids=['boundary'],
dim=2,
nu=nu,
rho0=ro,
gy=-g)
s = SchemeChooser(default='wcsph',
wcsph=wcsph,
aha=aha,
edac=edac,
iisph=iisph)
return s