Solid = MaterialProperties(Mesh, Eprime, anisotropic_K1c=True, 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)
# injection parameters Q0 = 0.001 # injection rate Injection = InjectionProperties(Q0, Mesh) # fluid properties Fluid = FluidProperties(viscosity=0) # simulation properties simulProp = SimulationProperties() simulProp.finalTime = 1e2 # the time at which the simulation stops simulProp.set_tipAsymptote( 'K' ) # the tip asymptote is evaluated with the toughness dominated assumption simulProp.set_volumeControl( True ) # use the inviscid fluid solver(toughness dominated), imposing volume balance simulProp.set_outputFolder( "./Data/K_radial_symmetric") # the disk address where the files are saved simulProp.symmetric = True # assume fracture geometry to be symmetric (only available for volume control) # initializing fracture Fr_geometry = Geometry('radial', radius=0.15) init_param = InitializationParameters(Fr_geometry, regime='K') # creating fracture object Fr = Fracture(Mesh, init_param, Solid, Fluid, Injection, simulProp) # create a Controller controller = Controller(Fr, Solid, Fluid, Injection, simulProp)