K1c_func=K1c_func)
# injection parameters
Q0 = 0.001 # injection rate
Injection = InjectionProperties(Q0, Mesh)
# fluid properties
Fluid = FluidProperties(viscosity=1.1e-5)
# simulation properties
simulProp = SimulationProperties()
simulProp.finalTime = 500 # the time at which the simulation stops
simulProp.set_volumeControl(
True) # to set up the solver in volume control mode (inviscid fluid)
simulProp.tolFractFront = 4e-3 # increase tolerance for the anisotropic case
simulProp.remeshFactor = 1.5 # the factor by which the mesh will be compressed.
simulProp.set_outputFolder(
"./Data/ellipse") # the disk address where the files are saved
simulProp.set_simulation_name('anisotropic_toughness_benchmark')
simulProp.symmetric = True # solving with faster solver that assumes fracture is symmetric
# initializing fracture
gamma = (K1c_func(np.pi / 2) / K1c_func(0))**2 # gamma = (Kc1/Kc3)**2
Fr_geometry = Geometry('elliptical', minor_axis=2., gamma=gamma)
init_param = InitializationParameters(Fr_geometry, regime='E_K')
# creating fracture object
Fr = Fracture(Mesh, init_param, Solid, Fluid, Injection, simulProp)
# create a Controller
controller = Controller(Fr, Solid, Fluid, Injection, simulProp)
# simulation properties
simulProp = SimulationProperties()
simulProp.finalTime = 1000 # the time at which the simulation stops
simulProp.set_volumeControl(
True) # to set up the solver in volume control mode (inviscid fluid)
simulProp.tolFractFront = 4e-3 # increase tolerance for the anisotropic case
simulProp.frontAdvancing = "implicit" # set the front advancing scheme to implicit
simulProp.set_tipAsymptote('K') # set the tip asymptote to toughness dominated
simulProp.set_outputFolder(
"./data/TI_elasticity_ellipse") # setting the output folder
simulProp.set_simulation_name(
'TI_ellasticy_benchmark') # setting the simulation name
simulProp.TI_KernelExecPath = '../TI_Kernel/build/' # path to the executable that calculates TI kernel
simulProp.symmetric = True # solving with faster solver that assumes fracture is symmetric
simulProp.remeshFactor = 1.5 # the factor by which the domain is expanded
# initialization parameters
Fr_geometry = Geometry('elliptical', minor_axis=1, gamma=gamma)
init_param = InitializationParameters(Fr_geometry, regime='E_E')
# # 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()
####################
