return abs(x) < 75 and (y >= -75. - tolerance and y <= -75. + tolerance)
# injection parameters
Q0 = 0.001 # injection rate
Injection = InjectionProperties(Q0, Mesh, source_loc_func=source_location)
# fluid properties
Fluid = FluidProperties(viscosity=1.1e-3, density=1000)
# simulation properties
simulProp = SimulationProperties()
simulProp.finalTime = 1.1e4 # the time at which the simulation stops
simulProp.set_outputFolder(
"./Data/buoyant_line_source") # the disk address where the files are saved
simulProp.gravity = True # take the effect of gravity into account
simulProp.set_mesh_extension_direction(['top'])
simulProp.plotVar = ['w', 'regime']
simulProp.toleranceEHL = 1e-3
# initializing fracture
surv_cells, _, inner_cells = get_eliptical_survey_cells(Mesh,
80,
20,
center=[0.0, -75.0])
surv_cells_dist = surv_cells * 0 + (Mesh.hx**2 + Mesh.hy**2)**(1 / 2)
Fr_geometry = Geometry(shape='level set',
survey_cells=surv_cells,
tip_distances=surv_cells_dist,
inner_cells=inner_cells)
confining_stress_func=sigmaO_func,
minimum_width=1e-5)
# injection parameters
Q0 = np.asarray([[0.0, 500], [2000, 0]]) # injection rate
Injection = InjectionProperties(Q0, Mesh)
# fluid properties
Fluid = FluidProperties(viscosity=30, density=2400)
# simulation properties
simulProp = SimulationProperties()
simulProp.finalTime = 560000 # the time at which the simulation stops
simulProp.set_outputFolder(
"./Data/neutral_buoyancy") # the disk address where the files are saved
simulProp.gravity = True # set up the gravity flag
simulProp.tolFractFront = 3e-3 # increase the tolerance for fracture
# front iteration
simulProp.plotTSJump = 4 # plot every fourth time step
simulProp.saveTSJump = 2 # save every second time step
simulProp.maxSolverItrs = 200 # increase the Anderson iteration limit for the
# elastohydrodynamic solver
simulProp.tmStpPrefactor = np.asarray([[0, 80000], [0.5, 0.1]
]) # set up the time step prefactor
simulProp.timeStepLimit = 5000 # time step limit
simulProp.plotVar = ['w',
'v'] # plot fracture width and fracture front velocity
simulProp.set_mesh_extension_direction(
['top', 'horizontal']) # allow the fracture to extend in positive y and x
simulProp.set_mesh_extension_factor(1.2) # set the extension factor to 1.4
simulProp.useBlockToeplizCompression = True # use the Toepliz elasticity matrix to save memory
