Example #1
0
    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)
Example #2
0
 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)
Example #3
0
    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)
Example #4
0
 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)
Example #6
0
    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)
Example #7
0
 def __init__(self):
     mesh = Rectangle(0.0,0.0,meshlength,meshheight,nx,ny)
     FSI.__init__(self,mesh)
Example #8
0
 def __init__(self):
     mesh = Rectangle(0.0, 0.0, meshlength, meshheight, nx, ny)
     FSI.__init__(self, mesh)