Q0 = 0.05  # injection rate
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
Example #2
0
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')


# creating fracture object
Fr = Fracture(Mesh,
              init_param,
              Solid,
              Fluid,
              Injection,