def create_solver(self):
kernel = CubicSpline(dim=2)
integrator = EPECIntegrator(mg=DEMStep(), al=DEMStep())
dt = 1e-6
print("DT: %s" % dt)
tf = 0.0013
solver = Solver(kernel=kernel,
dim=3,
integrator=integrator,
dt=dt,
tf=tf,
adaptive_timestep=False)
return solver
def create_solver(self):
kernel = CubicSpline(dim=2)
integrator = EPECIntegrator(al=DEMStep(), nylon=DEMStep())
dt = 1e-7
print("DT: %s" % dt)
tf = 0.0033
solver = Solver(kernel=kernel,
dim=3,
integrator=integrator,
dt=dt,
tf=tf,
adaptive_timestep=False,
pfreq=1000)
return solver
def create_solver(self):
kernel = CubicSpline(dim=2)
integrator = EPECIntegrator(cylinders=DEMStep(), t_wall=DEMStep())
# dt = 1e-4
dt = 1e-6
print("DT: %s" % dt)
tf = 0.51
solver = Solver(kernel=kernel,
dim=2,
integrator=integrator,
dt=dt,
tf=tf,
adaptive_timestep=False)
return solver