def run_radial_simul(self, my_front_reconstruction, my_front_advancement,
my_vertex_or_path, my_param):
# setting up the verbosity level of the log at console
# setup_logging_to_console(verbosity_level='error')
outputfolder = "./Temp_Data/" + my_vertex_or_path + "_radial_" + my_front_advancement + "_" + my_front_reconstruction
self.remove(outputfolder)
# creating mesh
Mesh = CartesianMesh(my_param['Lx'], my_param['Ly'], my_param['Nx'],
my_param['Ny'])
# solid properties
nu = my_param['nu'] # Poisson's ratio
youngs_mod = my_param['youngs_mod'] # Young's modulus
Eprime = youngs_mod / (1 - nu**2) # plain strain modulus
K_Ic = my_param['K_Ic'] # fracture toughness
Cl = my_param['Cl'] # Carter's leak off coefficient
# material properties
Solid = MaterialProperties(Mesh, Eprime, K_Ic, Carters_coef=Cl)
# injection parameters
Q0 = my_param['Q0'] # injection rate
Injection = InjectionProperties(Q0, Mesh)
# fluid properties
Fluid = FluidProperties(viscosity=my_param['viscosity'])
# simulation properties
simulProp = SimulationProperties()
simulProp.finalTime = my_param[
'finalTime'] # the time at which the simulation stops
simulProp.set_tipAsymptote(
my_vertex_or_path
) # tip asymptote is evaluated with the viscosity dominated assumption
simulProp.frontAdvancing = my_front_advancement # to set explicit front tracking
simulProp.plotFigure = False
simulProp.set_solTimeSeries(np.asarray([2, 200, 5000, 30000, 100000]))
simulProp.saveTSJump, simulProp.plotTSJump = 5, 5 # save and plot after every five time steps
simulProp.set_outputFolder(outputfolder)
simulProp.projMethod = my_front_reconstruction
simulProp.log2file = False
# initialization parameters
Fr_geometry = Geometry('radial', radius=my_param['initialR'])
init_param = InitializationParameters(Fr_geometry,
regime=my_vertex_or_path)
# creating fracture object
Fr = Fracture(Mesh, init_param, Solid, Fluid, Injection, simulProp)
# create a Controller
controller = Controller(Fr, Solid, Fluid, Injection, simulProp)
# run the simulation
exitcode = controller.run()
return exitcode, outputfolder
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)
# run the simulation
controller.run()
# fluid properties
Fluid = FluidProperties(viscosity=1e-3)
# simulation properties
simulProp = SimulationProperties()
simulProp.finalTime = 1.8e4 # 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 = 3 # set to plot every four time steps
simulProp.plotVar = ['regime', '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
simulProp.set_solTimeSeries = np.asarray(
[7672, 9660, 12435, 14693, 15342, 15835])
# 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
controller = Controller(Fr, Solid, Fluid, Injection, simulProp)
# run the simulation
controller.run()
####################