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)