def test_add_entry(self):
prjfile = "tunnel_ogs6py_add_entry.prj"
model = ogs6py.OGS(INPUT_FILE="tests/tunnel_ogs6py.prj", PROJECT_FILE=prjfile)
model.add_entry(tag="geometry", parent_xpath=".", text="geometry.gml")
model.write_input()
root = ET.parse(prjfile)
find = root.findall("./geometry")
self.assertEqual("geometry.gml", find[0].text)
def test_remove_element(self):
prjfile = "tunnel_ogs6py_remove_element.prj"
model = ogs6py.OGS(INPUT_FILE="tests/tunnel_ogs6py.prj", PROJECT_FILE=prjfile)
model.remove_element(xpath="./parameters/parameter[name='E']")
model.write_input()
root = ET.parse(prjfile)
find = root.findall("./parameters/parameter[name='E']/value")
self.assertEqual(0, len(find))
def test_replace_parameter(self):
prjfile = "tunnel_ogs6py_replace.prj"
model = ogs6py.OGS(INPUT_FILE="tests/tunnel_ogs6py.prj", PROJECT_FILE=prjfile)
model.replace_parameter(name="E", value=32)
model.write_input()
root = ET.parse(prjfile)
find = root.findall("./parameters/parameter[name='E']/value")
self.assertEqual("32", find[0].text)
def test_replace_medium_property(self):
prjfile = "tunnel_ogs6py_replace.prj"
model = ogs6py.OGS(INPUT_FILE="tests/tunnel_ogs6py.prj", PROJECT_FILE=prjfile)
model.replace_medium_property(mediumid=0, name="porosity", value=42)
model.write_input()
root = ET.parse(prjfile)
find = root.findall("./media/medium/properties/property[name='porosity']/value")
self.assertEqual("42", find[0].text)
def test_replace_phase_property(self):
prjfile = "tunnel_ogs6py_replace.prj"
model = ogs6py.OGS(INPUT_FILE="tests/tunnel_ogs6py.prj", PROJECT_FILE=prjfile)
model.replace_phase_property(mediumid=0, phase="Solid", name="thermal_expansivity", value=5)
model.write_input()
root = ET.parse(prjfile)
find = root.findall("./media/medium/phases/phase[type='Solid']/properties/property[name='thermal_expansivity']/value")
self.assertEqual("5", find[0].text)
def test_replace_text(self):
prjfile = "tunnel_ogs6py_replace.prj"
model = ogs6py.OGS(INPUT_FILE="tests/tunnel_ogs6py.prj", PROJECT_FILE=prjfile)
model.replace_text("tunnel_replace", xpath="./time_loop/output/prefix")
model.write_input()
root = ET.parse(prjfile)
find = root.findall("./time_loop/output/prefix")
self.assertEqual("tunnel_replace", find[0].text)
def test_add_block(self):
prjfile = "tunnel_ogs6py_add_block.prj"
model = ogs6py.OGS(INPUT_FILE="tests/tunnel_ogs6py.prj", PROJECT_FILE=prjfile)
model.add_block("parameter", parent_xpath="./parameters", taglist=["name", "type", "value"],
textlist=["mu","Constant","0.001"])
model.write_input()
root = ET.parse(prjfile)
find = root.findall("./parameters/parameter[name='mu']/value")
self.assertEqual("0.001", find[0].text)
def test_buildfromscratch(self):
model = ogs6py.OGS(PROJECT_FILE="tunnel_ogs6py.prj", MKL=True, OMP_NUM_THREADS=4)
model.mesh.add_mesh(filename="tunnel.vtu")
model.mesh.add_mesh(filename="tunnel_left.vtu")
model.mesh.add_mesh(filename="tunnel_right.vtu")
model.mesh.add_mesh(filename="tunnel_bottom.vtu")
model.mesh.add_mesh(filename="tunnel_top.vtu")
model.mesh.add_mesh(filename="tunnel_inner.vtu")
model.processes.set_process(
name="THERMO_RICHARDS_MECHANICS",
type="THERMO_RICHARDS_MECHANICS",
integration_order="3",
specific_body_force="0 0",
initial_stress="Initial_stress")
model.processes.set_constitutive_relation(type="LinearElasticIsotropic",
youngs_modulus="E",
poissons_ratio="nu")
model.processes.add_process_variable(process_variable="displacement",
process_variable_name="displacement")
model.processes.add_process_variable(process_variable="pressure",
process_variable_name="pressure")
model.processes.add_process_variable(process_variable="temperature",
process_variable_name="temperature")
model.processes.add_process_variable(secondary_variable="sigma",
output_name="sigma")
model.processes.add_process_variable(secondary_variable="epsilon",
output_name="epsilon")
model.processes.add_process_variable(secondary_variable="velocity",
output_name="velocity")
model.processes.add_process_variable(secondary_variable="saturation",
output_name="saturation")
model.media.add_property(medium_id="0",
phase_type="AqueousLiquid",
name="specific_heat_capacity",
type="Constant",
value="4280.0")
model.media.add_property(medium_id="0",
phase_type="AqueousLiquid",
name="thermal_conductivity",
type="Constant",
value="0.6")
model.media.add_property(medium_id="0",
phase_type="AqueousLiquid",
name="density",
type="Linear",
reference_value="999.1",
variable_name="phase_pressure",
reference_condition="1e5",
slope="4.5999999999999996e-10")
model.media.add_property(medium_id="0",
phase_type="AqueousLiquid",
name="thermal_expansivity",
type="Constant",
value="3.98e-4")
model.media.add_property(medium_id="0",
phase_type="AqueousLiquid",
name="viscosity",
type="Constant",
value="1.e-3")
model.media.add_property(medium_id="0",
name="permeability",
type="Constant",
value="1e-17")
model.media.add_property(medium_id="0",
name="porosity",
type="Constant",
value="0.15")
model.media.add_property(medium_id="0",
phase_type="Solid",
name="density",
type="Constant",
value="2300")
model.media.add_property(medium_id="0",
phase_type="Solid",
name="thermal_conductivity",
type="Constant",
value="1.9")
model.media.add_property(medium_id="0",
phase_type="Solid",
name="specific_heat_capacity",
type="Constant",
value="800")
model.media.add_property(medium_id="0",
name="biot_coefficient",
type="Constant",
value="0.6")
model.media.add_property(medium_id="0",
phase_type="Solid",
name="thermal_expansivity",
type="Constant",
value="1.7e-5")
model.media.add_property(medium_id="0",
name="thermal_conductivity",
type="EffectiveThermalConductivityPorosityMixing")
model.media.add_property(medium_id="0",
name="saturation",
type="Constant",
value="1")
model.media.add_property(medium_id="0",
name="relative_permeability",
type="Constant",
value="1")
model.media.add_property(medium_id="0",
name="bishops_effective_stress",
type="BishopsPowerLaw",
exponent="1")
model.timeloop.add_process(process="THERMO_RICHARDS_MECHANICS",
nonlinear_solver_name="nonlinear_solver",
convergence_type="PerComponentDeltaX",
norm_type="NORM2",
abstols="1e-4 1e-4 1e-10 1e-10",
time_discretization="BackwardEuler")
model.timeloop.set_stepping(process="THERMO_RICHARDS_MECHANICS",
type="IterationNumberBasedTimeStepping",
t_initial=0,
t_end=8,
initial_dt=0.1,
minimum_dt=1e-7,
maximum_dt=0.1,
number_iterations=[1, 4, 10, 20],
multiplier=[1.2, 1.0, 0.9, 0.8])
model.timeloop.add_output(
type="VTK",
prefix="tunnel",
repeat="10000",
each_steps="1",
variables=["displacement", "pressure", "temperature",
"sigma", "epsilon", "velocity", "saturation"],
fixed_output_times=[1, 2, 3],
suffix="_ts_{:timestep}_t_{:time}")
model.parameters.add_parameter(name="Initial_stress", type="Function",
mesh="tunnel", expression=["-5e6", "-5e6", "-5e6", "0"])
model.parameters.add_parameter(name="E", type="Constant", value="2e9")
model.parameters.add_parameter(name="nu", type="Constant", value="0.3")
model.parameters.add_parameter(name="T0", type="Constant", value="273.15")
model.parameters.add_parameter(name="displacement0",
type="Constant",
values="0 0")
model.parameters.add_parameter(name="pressure_ic", type="Constant", value="1e6")
model.parameters.add_parameter(name="dirichlet0", type="Constant", value="0")
model.parameters.add_parameter(name="temperature_ic",
type="Constant",
value="293.15")
model.parameters.add_parameter(name="pressure_bc",
type="CurveScaled",
curve="excavation_curve",
parameter="pressure_ic")
model.parameters.add_parameter(name="PressureLoad",
type="CurveScaled",
curve="excavation_curve",
parameter="PressureLoadValue")
model.parameters.add_parameter(
name="PressureLoadValue", type="Constant", value="-5.6e6")
model.parameters.add_parameter(name="heat_bv",
type="Constant",
value="25")
model.curves.add_curve(name="excavation_curve", coords=[0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0],
values=[1.0, 1.0, 0.8, 0.6, 0.4, 0.2, 0.1, 0.0])
model.processvars.set_ic(process_variable_name="temperature",
components="1",
order="1",
initial_condition="temperature_ic")
model.processvars.add_bc(process_variable_name="temperature",
mesh="tunnel_right",
type="Dirichlet",
component="0",
parameter="temperature_ic")
model.processvars.add_bc(process_variable_name="temperature",
mesh="tunnel_left",
type="Dirichlet",
component="0",
parameter="temperature_ic")
model.processvars.add_bc(process_variable_name="temperature",
mesh="tunnel_bottom",
type="Dirichlet",
component="0",
parameter="temperature_ic")
model.processvars.add_bc(process_variable_name="temperature",
mesh="tunnel_top",
type="Dirichlet",
component="0",
parameter="temperature_ic")
model.processvars.set_ic(process_variable_name="displacement",
components="2",
order="1",
initial_condition="displacement0")
model.processvars.add_bc(process_variable_name="displacement",
mesh="tunnel_right",
type="Dirichlet",
component="0",
parameter="dirichlet0")
model.processvars.add_bc(process_variable_name="displacement",
mesh="tunnel_left",
type="Dirichlet",
component="0",
parameter="dirichlet0")
model.processvars.add_bc(process_variable_name="displacement",
mesh="tunnel_top",
type="Dirichlet",
component="1",
parameter="dirichlet0")
model.processvars.add_bc(process_variable_name="displacement",
mesh="tunnel_bottom",
type="Dirichlet",
component="1",
parameter="dirichlet0")
model.processvars.add_bc(process_variable_name="displacement",
mesh="tunnel_inner",
type="NormalTraction",
parameter="PressureLoad")
model.processvars.set_ic(process_variable_name="pressure",
components="1",
order="1",
initial_condition="pressure_ic")
model.processvars.add_bc(process_variable_name="pressure",
mesh="tunnel_right",
type="Dirichlet",
component="0",
parameter="pressure_ic")
model.processvars.add_bc(process_variable_name="pressure",
mesh="tunnel_left",
type="Dirichlet",
component="0",
parameter="pressure_ic")
model.processvars.add_bc(process_variable_name="pressure",
mesh="tunnel_top",
type="Dirichlet",
component="0",
parameter="pressure_ic")
model.processvars.add_bc(process_variable_name="pressure",
mesh="tunnel_bottom",
type="Dirichlet",
component="0",
parameter="pressure_ic")
model.processvars.add_bc(process_variable_name="pressure",
mesh="tunnel_inner",
type="Dirichlet",
component="0",
parameter="pressure_bc")
model.nonlinsolvers.add_non_lin_solver(name="nonlinear_solver",
type="Newton",
max_iter="50",
linear_solver="general_linear_solver")
model.linsolvers.add_lin_solver(name="general_linear_solver",
kind="eigen",
solver_type="PardisoLU")
model.write_input()
with open("tunnel_ogs6py.prj", "rb") as f:
file_hash = hashlib.md5()
while chunk := f.read(8192):
file_hash.update(chunk)