def __init__(self): # Create mesh mesh = Rectangle(0.0, 0.0, 2.0, 1.0, 10*2**ref, 5*2**ref) # Create analytic expressions for various forces self.f_F = Expression(cpp_f_F, degree=2) self.g_F = Expression(cpp_g_F, degree=1) self.F_S = Expression(cpp_F_S, degree=2) self.F_M = Expression(cpp_F_M, degree=3) self.G_0 = Expression(cpp_G_S0, degree=5) # Exact solutions self.u_F = Expression(cpp_u_F, degree=2) self.p_F = Expression(cpp_p_F, degree=1) self.U_S = Expression(cpp_U_S, degree=2) self.P_S = Expression(cpp_P_S, degree=2) self.U_M = Expression(cpp_U_M, degree=3) # Initialize expressions forces = [self.f_F, self.g_F, self.F_S, self.F_M, self.G_0] solutions = [self.u_F, self.p_F, self.U_S, self.P_S, self.U_M] for f in forces + solutions: f.C = C f.t = 0.0 # Initialize base class FSI.__init__(self, mesh, application_parameters)
def __init__(self): mesh = Rectangle(0.0, 0.0, vessel_length, vessel_height, nx, ny, "crossed") self.P_Fwave = Expression(cpp_P_Fwave) self.P_Fwave.C = C self.endtime = 70.0 self.rho_S = 4.0 FSI.__init__(self,mesh,application_parameters)
def __init__(self): # Create mesh mesh = Rectangle(0.0, 0.0, 2.0, 1.0, 10 * 2**ref, 5 * 2**ref) # Create analytic expressions for various forces self.f_F = Expression(cpp_f_F, degree=2) self.g_F = Expression(cpp_g_F, degree=1) self.F_S = Expression(cpp_F_S, degree=2) self.F_M = Expression(cpp_F_M, degree=3) self.G_0 = Expression(cpp_G_S0, degree=5) # Exact solutions self.u_F = Expression(cpp_u_F, degree=2) self.p_F = Expression(cpp_p_F, degree=1) self.U_S = Expression(cpp_U_S, degree=2) self.P_S = Expression(cpp_P_S, degree=2) self.U_M = Expression(cpp_U_M, degree=3) # Initialize expressions forces = [self.f_F, self.g_F, self.F_S, self.F_M, self.G_0] solutions = [self.u_F, self.p_F, self.U_S, self.P_S, self.U_M] for f in forces + solutions: f.C = C f.t = 0.0 # Initialize base class FSI.__init__(self, mesh, application_parameters)
def __init__(self): mesh = Rectangle(0.0, 0.0, vessel_length, vessel_height, nx, ny, "crossed") self.P_Fwave = Expression(cpp_P_Fwave) self.P_Fwave.C = C self.endtime = 70.0 self.rho_S = 4.0 FSI.__init__(self, mesh, application_parameters)
def __init__(self): # Create mesh ny = 5 nx = 20 if application_parameters["crossed_mesh"]: mesh = Rectangle(0.0, 0.0, channel_length, channel_height, nx, ny, "crossed") else: mesh = Rectangle(0.0, 0.0, channel_length, channel_height, nx, ny) for i in range(3): mesh = refine(mesh) # Set material parameters self.E = 100.0 self.nu = 0.3 # Initialize base class FSI.__init__(self, mesh, application_parameters)
def __init__(self): mesh = Rectangle(0.0,0.0,meshlength,meshheight,nx,ny) FSI.__init__(self,mesh)
def __init__(self): mesh = Rectangle(0.0, 0.0, meshlength, meshheight, nx, ny) FSI.__init__(self, mesh)