添加参数文件

[mapengfei-nwpu]

from fenics import *
# parameters for fluid.
fluid_parameters = {
    'time_step' : 1e-4,
    'time_start': 0,
    'time_end'  : 10,
    'mesh_size' : [32, 32, 32],
    'viscosity' : 0.04,
    'point_0'   : [-1.0, -2.0,-2.0],
    'point_1'   : [ 3.0,  2.0, 2.0]
}

# parameters for solid mechanics.
pa1 = { 
    'C' : 10, # kPa
    'bf' : 8.0, # kPa
    'bt' : 2.0,
    'bfs' : 4.0,
    'kappa' : 5e6,
    'Tactive' : 60 # kPa
}

if __name__ == "__main__":
    # from LeftVentricleParams import fluid_parameters
    from fenics import *
    mesh_size  = fluid_parameters["mesh_size"]
    dt         = fluid_parameters["time_step"]
    T          = fluid_parameters["time_end"] - fluid_parameters["time_start"]
    nu         = fluid_parameters["viscosity"]
    point_0    = Point(fluid_parameters["point_0"])
    point_1    = Point(fluid_parameters["point_1"])
    fluid_mesh = BoxMesh(point_0, point_1, mesh_size[0], mesh_size[1], mesh_size[2])