def initialize(self):
self.geom = DamBreak3DGeometry(
container_height=5, container_width=20, container_length=20,
fluid_column_height=4, fluid_column_width=3, fluid_column_length=3,
obstacle_center_x=0, obstacle_center_y=0, # parameters object .stl
#obstacle_length= , obstacle_height= , obstacle_width= ,
dx=dx, nboundary_layers=nboundary_layers, hdx=hdx, rho0=rho0,
with_obstacle=False)
def consume_user_options(self):
dx = self.options.dx
self.dx = dx
self.hdx = self.options.hdx
self.geom = DamBreak3DGeometry(dx=dx,
nboundary_layers=nboundary_layers,
hdx=self.hdx,
rho0=ro)
self.co = 10.0 * self.geom.get_max_speed(g=9.81)
def consume_user_options(self):
dx = self.options.dx
self.dx = dx
self.hdx = self.options.hdx
self.h0 = self.hdx * self.dx
self.geom = DamBreak3DGeometry(container_height=0.4,
container_width=0.61,
container_length=1.6,
fluid_column_height=0.3,
fluid_column_width=0.61,
fluid_column_length=0.4,
obstacle_center_x=0.96,
obstacle_center_y=0,
obstacle_length=0.12,
obstacle_height=0.75,
obstacle_width=0.12,
nboundary_layers=nboundary_layers,
with_obstacle=True,
dx=dx,
hdx=hdx,
rho0=ro)
self.co = c0
def create_particles(self):
geom = DamBreak3DGeometry(container_height=height,
container_width=width,
container_length=length,
fluid_column_height=height * 0.75,
fluid_column_width=width,
fluid_column_length=length,
nboundary_layers=n_layers,
with_obstacle=False,
dx=self.dx,
hdx=self.hdx,
rho0=rho)
[fluid, boundary] = geom.create_particles()
fluid.x = fluid.x - length * 0.5
boundary.x = boundary.x - length * 0.5
# Setting up intital velocity of the tank
omega0 = theta_0 * omega_r
boundary.u, boundary.w = boundary.z * omega0, -boundary.x * omega0
self.scheme.setup_properties([fluid, boundary])
particles = [fluid, boundary]
return particles
def initialize(self):
self.geom = DamBreak3DGeometry(dx=dx,
nboundary_layers=nboundary_layers,
hdx=hdx,
rho0=rho0,
with_obstacle=False)
def initialize(self):
self.geom = DamBreak3DGeometry(dx=dx,
nboundary_layers=nboundary_layers,
hdx=hdx,
rho0=ro)
self.co = 10.0 * self.geom.get_max_speed(g=9.81)