def create_scheme(self): Re = 420 nu = 0.125 / Re c0 = 1.25 hdx = 1.2 wcsph = WCSPHScheme( ['fluid'], ['boundary'], dim=2, c0=c0, h0=hdx, hdx=hdx, rho0=1, alpha=0.1, beta=0, gamma=7, ) dt = 5e-3 tf = 15 wcsph.configure_solver( tf=tf, dt=dt, ) return wcsph
def create_scheme(self): s = WCSPHScheme(['fluid'], ['top_wall', 'bottom_wall'], dim=3, rho0=rho0, c0=co, h0=dx * hdx, hdx=hdx, gamma=7.0, alpha=0.5, beta=0.0, nu=mu / rho0, tensile_correction=True) kernel = WendlandQuintic(dim=3) dt_cfl = 0.25 * (hdx * dx) / co dt_viscous = 0.125 * (hdx * dx)**2 * rho0 / mu print("CFL based time step: %.3E" % dt_cfl) print("Viscous time step: %.3E" % dt_viscous) dt = min(dt_cfl, dt_viscous) s.configure_solver(kernel=kernel, integrator_cls=EPECIntegrator, tf=tf, dt=dt, adaptive_timestep=True, n_damp=50, pfreq=200) return s
def create_scheme(self): s = WCSPHScheme( ['fluid'], [], dim=2, rho0=self.ro, c0=self.co, h0=self.dx*self.hdx, hdx=self.hdx, gamma=7.0, alpha=0.1, beta=0.0 ) dt = 5e-6 tf = 0.0076 s.configure_solver(dt=dt, tf=tf) return s
def create_scheme(self): s = WCSPHScheme( ['fluid'], ['solid'], dim=3, rho0=1000, c0=c0, h0=hdx*dx, hdx=hdx, gz=-9.81, alpha=0.25, beta=0.0, gamma=7, hg_correction=True, tensile_correction=False ) extra_steppers = dict(solid=OneStageRigidBodyStep()) s.configure_solver( extra_steppers=extra_steppers, tf=10.0, dt=dt ) return s
def create_scheme(self): s = WCSPHScheme( ['fluid'], ['boundary'], dim=dim, rho0=ro, c0=self.co, h0=h0, hdx=hdx, gz=-9.81, alpha=alpha, beta=beta, gamma=gamma, hg_correction=True, tensile_correction=True ) kernel = WendlandQuintic(dim=dim) s.configure_solver( kernel=kernel, integrator_cls=EPECIntegrator, tf=tf, dt=dt, adaptive_timestep=True, n_damp=50 ) return s
def create_scheme(self): s = WCSPHScheme( ['fluid'], [], dim=2, rho0=self.ro, c0=self.co, h0=self.dx*self.hdx, hdx=self.hdx, gamma=7.0, alpha=0.1, beta=0.0 ) kernel = Gaussian(dim=2) dt = 5e-6; tf = 0.0076 s.configure_solver( kernel=kernel, integrator_cls=EPECIntegrator, dt=dt, tf=tf, adaptive_timestep=True, cfl=0.3, n_damp=50, output_at_times=[0.0008, 0.0038] ) return s
def create_scheme(self): s = WCSPHScheme(['fluid'], [], dim=2, rho0=1.0, c0=1400, h0=1.3 * 0.025, hdx=1.3, gamma=7.0, alpha=0.1, beta=0.0) dt = 5e-6 tf = 0.0076 s.configure_solver( dt=dt, tf=tf, ) return s