def set_rock_and_fluid(self):
"""
Set rock and fluid properties to granite and water.
"""
self.rock = pp.Granite()
self.rock.FRICTION_COEFFICIENT = 0.5
self.fluid = pp.Water()
self.rock.PERMEABILITY = 1e-16
def set_rock_and_fluid(self):
"""
Set rock and fluid properties to those of granite and water.
We ignore all temperature effects.
Credits: PorePy paper
"""
super().set_rock()
# Fluid. Temperature at ISC is 11 degrees average.
self.fluid = pp.Water(theta_ref=11)
def __init__(self, pb):
pb["theta_ref"] = pb["temperature_at_depth"]
pb["water"] = pp.Water(pb["theta_ref"])
pb["rock"] = pp.Granite(pb["theta_ref"])
# create the computational grid
gb = create_grid(pb)
# Darcy problem
self.assign_data(gb, pb, FlowData, "problem")
pb["file_name"] = "hydraulic_head"
flow = pp.EllipticModel(gb, **pb)
# transport problem
self.assign_data(gb, pb, TransportData, "transport_data")
pb["file_name"] = "temperature"
transport = TransportSolver(gb, **pb)
darcy_and_transport.DarcyAndTransport.__init__(self, flow, transport)
np.savetxt(folder + "/outflow.csv", advective._solver.outflow)
if __name__ == "__main__":
file_geo = "Algeroyna.csv"
aperture = pp.MILLIMETER
data = {"aperture": aperture, "kf": aperture ** 2 / 12}
folder = "upscaling"
upscaling(file_geo, data, folder, dfn=True)
file_geo = "Algeroyna.csv"
folder = "solution"
theta_ref = 30. * pp.CELSIUS
fluid = pp.Water(theta_ref)
rock = pp.Granite(theta_ref)
aperture = 2 * pp.MILLIMETER
# select the permeability depending on the selected test case
data = {
"aperture": aperture,
"kf": aperture ** 2 / 12,
"km": 1e7 * rock.PERMEABILITY,
"porosity_f": 0.85,
"dt": 1e6 * pp.SECOND,
"t_max": 1e7 * pp.SECOND,
"fluid": fluid,
"rock": rock,
}