Solid = MaterialProperties(Mesh,
Eprime,
K_Ic,
confining_stress_func=sigmaO_func)
# injection parameters
Q0 = 0.001 # injection rate
Injection = InjectionProperties(Q0, Mesh)
# fluid properties
Fluid = FluidProperties(viscosity=1.1e-3)
# simulation properties
simulProp = SimulationProperties()
simulProp.finalTime = 145. # the time at which the simulation stops
simulProp.bckColor = 'sigma0' # setting the parameter according to which the mesh is color coded
simulProp.set_outputFolder("./Data/height_contained")
simulProp.tmStpPrefactor = 1.0 # decreasing the size of time step
simulProp.plotVar = ['footprint'] # plotting footprint
# initializing fracture
Fr_geometry = Geometry(shape='radial', radius=1.)
init_param = InitializationParameters(Fr_geometry, regime='M')
# creating fracture object
Fr = Fracture(Mesh,
init_param,
Solid,
Fluid,
Injection,
simulProp)
Injection = InjectionProperties(Q0, Mesh, source_coordinates=[0.1, 0.])
# fluid properties
Fluid = FluidProperties(viscosity=1.1e-3)
# simulation properties
simulProp = SimulationProperties()
simulProp.finalTime = 0.0003 # the time at which the simulation stops
myfolder = "./Data/toughness_jump_3p6"
simulProp.set_outputFolder(
myfolder) # the disk address where the files are saved
simulProp.projMethod = 'LS_continousfront'
simulProp.frontAdvancing = 'implicit'
simulProp.useBlockToeplizCompression = True
simulProp.saveToDisk = False
simulProp.bckColor = 'K1c'
simulProp.saveFluidVelAsVector = True
# To decide what you will see when you print:
#simulProp.plotVar = ['ffvf','regime']
#simulProp.plotVar = ['footprint','regime']
#simulProp.plotVar = ['footprint']
# setting up mesh extension options
simulProp.meshExtensionAllDir = True
simulProp.maxElementIn = 10000
simulProp.set_mesh_extension_factor(1.1)
simulProp.set_mesh_extension_direction(['all'])
simulProp.meshReductionPossible = True
simulProp.maxCellSize = 0.1
Eprime,
K1c_func=My_KIc_func,
confining_stress_func=sigmaO_func,
minimum_width=1e-8)
# injection parameters
Q0 = 0.001 # injection rate
Injection = InjectionProperties(Q0, Mesh)
# fluid properties
Fluid = FluidProperties(viscosity=1.1e-3)
# simulation properties
simulProp = SimulationProperties()
simulProp.bckColor = 'sigma0'
simulProp.finalTime = 0.28 # the time at which the simulation stops
simulProp.outputTimePeriod = 1e-4 # to save after every time step
simulProp.tmStpPrefactor = 0.5 # decrease the pre-factor due to explicit front tracking
simulProp.set_outputFolder("./Data/stress_heterogeneities") # the disk address where the files are saved
simulProp.saveFluidFluxAsVector = True
simulProp.plotVar = ['ffvf']
simulProp.projMethod = 'LS_continousfront' # <--- mandatory use
simulProp.saveToDisk = True
simulProp.set_mesh_extension_factor(1.1)
simulProp.set_mesh_extension_direction(['all'])
simulProp.useBlockToeplizCompression = True
# initialization parameters
Fr_geometry = Geometry('radial', radius=0.12)
init_param = InitializationParameters(Fr_geometry, regime='M')
confining_stress_func=sigmaO_func,
Carters_coef=1e-6)
# injection parameters
Q0 = np.asarray([[0, 6000], [0.001, 0]])
Injection = InjectionProperties(Q0, Mesh, source_coordinates=[0, -20])
# fluid properties
Fluid = FluidProperties(viscosity=1e-3)
# simulation properties
simulProp = SimulationProperties()
simulProp.finalTime = 1.6e4 # the time at which the simulation stops
simulProp.set_outputFolder(
"./Data/fracture_closure") # the disk address where the files are saved
simulProp.bckColor = 'confining stress' # setting the parameter for the mesh color coding
simulProp.plotTSJump = 4 # set to plot every four time steps
simulProp.plotVar = ['w', 'lk', 'footprint'
] # setting the parameters that will be plotted
simulProp.tmStpPrefactor = np.asarray(
[[0, 6000], [0.8, 0.4]]) # decreasing the time step pre-factor after 6000s
simulProp.maxSolverItrs = 120 # increase maximum iterations for the elastohydrodynamic solver
# initialization parameters
Fr_geometry = Geometry('radial', radius=20)
init_param = InitializationParameters(Fr_geometry, regime='M')
# creating fracture object
Fr = Fracture(Mesh, init_param, Solid, Fluid, Injection, simulProp)
# create a Controller
Solid = MaterialProperties(Mesh,
Eprime,
confining_stress_func=sigmaO_func,
minimum_width=1e-8)
# injection parameters
Q0 = np.asarray([[0, 31, 151], [0.0009e-6, 0.0065e-6, 0.0023e-6]])
Injection = InjectionProperties(Q0, Mesh)
# fluid properties
Fluid = FluidProperties(viscosity=30)
# simulation properties
simulProp = SimulationProperties()
simulProp.bckColor = 'confining stress' # the parameter according to which the background is color coded
simulProp.frontAdvancing = 'explicit'
simulProp.set_outputFolder('./Data/Wu_et_al')
simulProp.set_solTimeSeries(np.asarray([22., 60., 144., 376., 665.]))
simulProp.plotVar = ['footprint']
# initializing fracture
Fr_geometry = Geometry('radial', radius=0.019)
init_param = InitializationParameters(Fr_geometry, regime='M')
# creating fracture object
Fr = Fracture(Mesh, init_param, Solid, Fluid, Injection, simulProp)
# create a Controller
controller = Controller(Fr, Solid, Fluid, Injection, simulProp)