def create_solver(self):
kernel = QuinticSpline(dim=2)
integrator = PECIntegrator(liquid=TransportVelocityStep(),
gas=TransportVelocityStep())
solver = Solver(
kernel=kernel, dim=dim, integrator=integrator,
dt=dt, tf=tf, adaptive_timestep=False)
return solver
def create_solver(self):
"""Set up the default integrator for fiber particles."""
kernel = QuinticSpline(dim=3)
integrator = EPECIntegrator(
fiber1=TransportVelocityStep(),
fiber2=TransportVelocityStep())
solver = Solver(
kernel=kernel, dim=3, integrator=integrator, dt=self.dt,
tf=self.tf, N=200)
return solver
def create_solver(self):
kernel = WendlandQuintic(dim=2)
integrator = PECIntegrator(fluid=TransportVelocityStep())
solver = Solver(
kernel=kernel, dim=dim, integrator=integrator,
dt=dt, tf=tf, adaptive_timestep=False)
return solver
def create_solver(self):
kernel = QuinticSpline(dim=2)
integrator = EPECIntegrator(fluid=TransportVelocityStep())
solver = Solver(
kernel=kernel, dim=dim, integrator=integrator,
dt=dt, tf=tf, adaptive_timestep=False,
output_at_times=[0., 0.08, 0.16, 0.26])
return solver
def create_solver(self):
kernel = QuinticSpline(dim=2)
integrator = Integrator(fluid=TransportVelocityStep(),
obstacle=TwoStageRigidBodyStep())
solver = Solver(kernel=kernel,
dim=2,
integrator=integrator,
tf=tf,
dt=dt,
adaptive_timestep=False,
output_at_times=[1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0])
return solver