def init_ogs_project(rad, task_root=None): # model setup model = OGS(task_root="pump_3d_steady", task_id="model", output_dir="out") model.pcs.add_block( # set the process type PCS_TYPE="GROUNDWATER_FLOW", NUM_TYPE="NEW", TIM_TYPE="STEADY") # generate a radial 3D mesh and conductivity field model.msh.generate("radial", dim=3, angles=angles, rad=rad, z_arr=-np.arange(11)) return model
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='Thies_quad_2d_root', task_id='Thies_quad_2d', output_dir='out', ) model.msh.read_file('Thies_quad_2d.msh') model.gli.read_file('Thies_quad_2d.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='GROUNDWATER_FLOW', NUM_TYPE='NEW', PRIMARY_VARIABLE='HEAD', ) model.rfd.read_file('Thies_quad_2d.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='GROUNDWATER_FLOW', PRIMARY_VARIABLE='HEAD', GEO_TYPE=['POLYLINE', 'left_bc'], DIS_TYPE=['CONSTANT', 20.0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='GROUNDWATER_FLOW', PRIMARY_VARIABLE='HEAD', GEO_TYPE=['POLYLINE', 'right_bc'], DIS_TYPE=['CONSTANT', 20.0], )
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='creep3d_root', task_id='creep3d', output_dir='out', ) model.msh.read_file('creep3d.msh') model.gli.read_file('creep3d.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='HEAT_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='DEFORMATION', ) model.rfd.read_file('creep3d.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='HEAT_TRANSPORT', PRIMARY_VARIABLE='TEMPERATURE1', GEO_TYPE=['SURFACE', 'TOP'], DIS_TYPE=['CONSTANT', 27.0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='HEAT_TRANSPORT', PRIMARY_VARIABLE='TEMPERATURE1',
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='2units2faults_FCT_root', task_id='2units2faults_FCT', output_dir='out', ) model.msh.read_file('2units2faults_FCT.msh') model.gli.read_file('2units2faults_FCT.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='LIQUID_FLOW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='HEAT_TRANSPORT', ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='HEAT_TRANSPORT', PRIMARY_VARIABLE='TEMPERATURE1', GEO_TYPE=['SURFACE', 'South'], DIS_TYPE=[ ['LINEAR', 4], [0, 40], [1, 80], [2, 80], [3, 40], ], )
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='elder_root', task_id='elder', output_dir='out', ) model.msh.read_file('elder.msh') model.gli.read_file('elder.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='LIQUID_FLOW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MASS_TRANSPORT', ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='LIQUID_FLOW', PRIMARY_VARIABLE='PRESSURE1', GEO_TYPE=['POLYLINE', 'LEFT_TOP_CORNER'], DIS_TYPE=['CONSTANT', 0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='MASS_TRANSPORT', PRIMARY_VARIABLE='CONCENTRATION1', GEO_TYPE=['SURFACE', 'TOP_HALF_RIGHT'], DIS_TYPE=['CONSTANT', 1.0],
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='3D_HGW_root', task_id='3D_HGW', output_dir='out', ) model.msh.read_file('3D_HGW.msh') model.gli.read_file('3D_HGW.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='GROUNDWATER_FLOW', ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='GROUNDWATER_FLOW', PRIMARY_VARIABLE='HEAD', GEO_TYPE=['SURFACE', 'L_SF'], DIS_TYPE=['CONSTANT', 10], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='GROUNDWATER_FLOW', PRIMARY_VARIABLE='HEAD', GEO_TYPE=['SURFACE', 'R_SF'], DIS_TYPE=['CONSTANT', 9], ) model.ic.add_block( main_key='INITIAL_CONDITION', PCS_TYPE='GROUNDWATER_FLOW',
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='1d_pri_root', task_id='1d_pri', output_dir='out', ) model.msh.read_file('1d_pri.msh') model.gli.read_file('1d_pri.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='GROUNDWATER_FLOW', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='FLUID_MOMENTUM', ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='GROUNDWATER_FLOW', PRIMARY_VARIABLE='HEAD', GEO_TYPE=['POINT', 'POINT2'], DIS_TYPE=['CONSTANT', 1], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='GROUNDWATER_FLOW', PRIMARY_VARIABLE='HEAD', GEO_TYPE=['POINT', 'POINT0'],
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='BGRa_root', task_id='BGRa', output_dir='out', ) model.msh.read_file('BGRa.msh') model.gli.read_file('BGRa.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='HEAT_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='DEFORMATION', ELEMENT_MATRIX_OUTPUT=0, ) model.rfd.read_file('BGRa.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='HEAT_TRANSPORT', PRIMARY_VARIABLE='TEMPERATURE1', GEO_TYPE=['POLYLINE', 'TOP'], DIS_TYPE=['CONSTANT', 27.0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='HEAT_TRANSPORT',
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='HBr_10C_Diff_new_root', task_id='HBr_10C_Diff_new', output_dir='out', ) model.msh.read_file('HBr_10C_Diff_new.msh') model.gli.read_file('HBr_10C_Diff_new.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='GROUNDWATER_FLOW', NUM_TYPE='NEW', BOUNDARY_CONDITION_OUTPUT=[], ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MASS_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='HEAT_TRANSPORT', TEMPERATURE_UNIT='KELVIN', ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='GROUNDWATER_FLOW', PRIMARY_VARIABLE='HEAD', GEO_TYPE=['POINT', 'POINT0'],
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='h_gas_ele_root', task_id='h_gas_ele', output_dir='out', ) model.msh.read_file('h_gas_ele.msh') model.gli.read_file('h_gas_ele.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='AIR_FLOW', NUM_TYPE='NEW', ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='AIR_FLOW', PRIMARY_VARIABLE='PRESSURE1', GEO_TYPE=['POINT', 'POINT0'], DIS_TYPE=['CONSTANT', 95500.0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='AIR_FLOW', PRIMARY_VARIABLE='PRESSURE1', GEO_TYPE=['POINT', 'POINT2'], DIS_TYPE=['CONSTANT', 95500.0], ) model.ic.add_block( main_key='INITIAL_CONDITION',
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='CO2phase_gen_E300_root', task_id='CO2phase_gen_E300', output_dir='out', ) model.msh.read_file('CO2phase_gen_E300.msh') model.gli.read_file('CO2phase_gen_E300.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MULTI_PHASE_FLOW', NUM_TYPE='NEW', USE_PRECALCULATED_FILES=[], SAVE_ECLIPSE_DATA_FILES=[], SIMULATOR='ECLIPSE', SIMULATOR_PATH='C:\ecl\2012.1\bin\pc_x86_64\e300.exe', SIMULATOR_MODEL_PATH='./eclipse/ECL.DATA', ELEMENT_MATRIX_OUTPUT=0, ST_RHS=1, BOUNDARY_CONDITION_OUTPUT=[], DISSOLVED_CO2_PCS_NAME='C(4)', DISSOLVED_CO2_INGAS_PCS_NAME='C(4)_gas', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MASS_TRANSPORT', NUM_TYPE='NEW', BOUNDARY_CONDITION_OUTPUT=[], )
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='m_e_transiso_2D_root', task_id='m_e_transiso_2D', output_dir='out', ) model.msh.read_file('m_e_transiso_2D.msh') model.gli.read_file('m_e_transiso_2D.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='DEFORMATION', ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_X1', GEO_TYPE=['POLYLINE', 'LEFT'], DIS_TYPE=['CONSTANT', 0.0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_Y1', GEO_TYPE=['POINT', 'POINT0'], DIS_TYPE=['CONSTANT', 0.0], ) model.st.add_block( main_key='SOURCE_TERM', PCS_TYPE='DEFORMATION',
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='h2_Liako_root', task_id='h2_Liako', output_dir='out', ) model.msh.read_file('h2_Liako.msh') model.gli.read_file('h2_Liako.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MULTI_PHASE_FLOW', NUM_TYPE='NEW', ELEMENT_MATRIX_OUTPUT=0, ) model.rfd.read_file('h2_Liako.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='MULTI_PHASE_FLOW', PRIMARY_VARIABLE='PRESSURE2', GEO_TYPE=['POINT', 'POINT1'], DIS_TYPE=['CONSTANT', 101325], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='MULTI_PHASE_FLOW', PRIMARY_VARIABLE='PRESSURE2', GEO_TYPE=['POINT', 'POINT0'], DIS_TYPE=['CONSTANT', 101325], )
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='m1_3Dload_root', task_id='m1_3Dload', output_dir='out', ) model.msh.read_file('m1_3Dload.msh') model.gli.read_file('m1_3Dload.gli') model.ddc.read_file('m1_3Dload.ddc') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='DEFORMATION', NUM_TYPE='NEW', ) model.rfd.read_file('m1_3Dload.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_X1', GEO_TYPE=['SURFACE', 'SURFACE1'], DIS_TYPE=[ ['LINEAR', 4], [0, 0.0], [1, 0.0], [2, 0.0], [3, 0.0], ], TIM_TYPE=['CURVE', 1], )
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='Ogata-Banks_root', task_id='Ogata-Banks', output_dir='out', ) model.msh.read_file('Ogata-Banks.msh') model.gli.read_file('Ogata-Banks.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='LIQUID_FLOW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='HEAT_TRANSPORT', TEMPERATURE_UNIT='KELVIN', ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='LIQUID_FLOW', PRIMARY_VARIABLE='PRESSURE1', GEO_TYPE=['POINT', 'POINT1'], DIS_TYPE=['CONSTANT', 100000], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='HEAT_TRANSPORT', PRIMARY_VARIABLE='TEMPERATURE1', GEO_TYPE=['POINT', 'POINT0'],
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='3D_UniAxiLoad_root', task_id='3D_UniAxiLoad', output_dir='out', ) model.msh.read_file('3D_UniAxiLoad.msh') model.gli.read_file('3D_UniAxiLoad.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='DEFORMATION', BOUNDARY_CONDITION_OUTPUT=[], ) model.rfd.read_file('3D_UniAxiLoad.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_Z1', GEO_TYPE=['SURFACE', 'BOTTOM'], DIS_TYPE=['CONSTANT', 0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_Z1', GEO_TYPE=['SURFACE', 'TOP'], DIS_TYPE=['CONSTANT', -5e-07], TIM_TYPE=['CURVE', 2], )
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='govin_line_root', task_id='govin_line', output_dir='out', ) model.msh.read_file('govin_line.msh') model.gli.read_file('govin_line.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='OVERLAND_FLOW', NUM_TYPE='NEW', ) model.ic.add_block( main_key='INITIAL_CONDITION', PCS_TYPE='OVERLAND_FLOW', PRIMARY_VARIABLE='HEAD', GEO_TYPE='DOMAIN', DIS_TYPE=['CONSTANT', 1e-06], ) model.st.add_block( main_key='SOURCE_TERM', PCS_TYPE='OVERLAND_FLOW', PRIMARY_VARIABLE='HEAD', GEO_TYPE=['POLYLINE', 'CHANNEL'], DIS_TYPE=['CONSTANT_NEUMANN', 0.004], ) model.st.add_block( main_key='SOURCE_TERM',
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='heat_transfer_root', task_id='heat_transfer', output_dir='out', ) model.msh.read_file('heat_transfer.msh') model.gli.read_file('heat_transfer.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='TNEQ', TEMPERATURE_UNIT='KELVIN', ELEMENT_MATRIX_OUTPUT=0, ) model.rfd.read_file('heat_transfer.rfd') model.ic.add_block( main_key='INITIAL_CONDITION', PCS_TYPE='TNEQ', PRIMARY_VARIABLE='PRESSURE1', GEO_TYPE='DOMAIN', DIS_TYPE=['CONSTANT', 100000.0], ) model.ic.add_block( main_key='INITIAL_CONDITION', PCS_TYPE='TNEQ', PRIMARY_VARIABLE='TEMPERATURE1', GEO_TYPE='DOMAIN', DIS_TYPE=['CONSTANT', 773.15], )
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='h_frac_line2D_root', task_id='h_frac_line2D', output_dir='out', ) model.msh.read_file('h_frac_line2D.msh') model.gli.read_file('h_frac_line2D.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='LIQUID_FLOW', NUM_TYPE='NEW', ELEMENT_MATRIX_OUTPUT=0, ) model.rfd.read_file('h_frac_line2D.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='LIQUID_FLOW', PRIMARY_VARIABLE='PRESSURE1', GEO_TYPE=['POLYLINE', 'BOTTOM'], DIS_TYPE=[ ['LINEAR', 2], [1, 10000.0], [2, 100000.0], ], ) model.ic.add_block( main_key='INITIAL_CONDITION', PCS_TYPE='LIQUID_FLOW',
name=str(overall_count) + '_border_' + str(border) + '_stor_' + str(storage) + '_rech_' + str(rech_abv) dim_no = 2 parent_dir = CWD + '/setup' # name of directory (entire path) of one single ogs run dire = parent_dir + '/' + name # make folders if not os.path.exists(parent_dir): os.mkdir(parent_dir) if not os.path.exists(dire): os.mkdir(dire) pcs_type_flow = 'GROUNDWATER_FLOW' var_name_flow = 'HEAD' t_id = 'transect' # ------------------------generate ogs base class----------------------------- # ogs = OGS(task_root=dire+"/", task_id=t_id, output_dir=dire+"/"+"steady"+"/") # ------------------------ MSH -------------------------------- # # generate a rectangular mesh in x-z-plane # first part of mesh, configuration based on an element size of 1x1 m, because variable "length" is used for number of nodes ogs.msh.generate("rectangular", dim=2, mesh_origin=(0., 0.), element_no=(border, n_cellsz), element_size=(1, 1)) ogs.msh.MATERIAL_ID = 0 # second part of mesh second_mesh = MSH() second_mesh.generate("rectangular", dim=2, mesh_origin=(border, 0.), element_no=(length - border, n_cellsz), element_size=(1, 1))
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='m_mises_root', task_id='m_mises', output_dir='out', ) model.msh.read_file('m_mises.msh') model.gli.read_file('m_mises.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='DEFORMATION', ) model.rfd.read_file('m_mises.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_X1', GEO_TYPE=['POLYLINE', 'LEFT'], DIS_TYPE=['CONSTANT', 0.0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_Y1', GEO_TYPE=['POLYLINE', 'BOTTOM'], DIS_TYPE=['CONSTANT', 0.0], ) model.st.add_block( main_key='SOURCE_TERM',
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='1D_TPF_resS_trans_root', task_id='1D_TPF_resS_trans', output_dir='out', ) model.msh.read_file('1D_TPF_resS_trans.msh') model.gli.read_file('1D_TPF_resS_trans.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='PS_GLOBAL', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MASS_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MASS_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MASS_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block(
output = "output" dim_no = 2 parent_dir = CWD + '/setup' # name of directory (entire path) of one single ogs run dire = parent_dir + '/' + name # make folders if not os.path.exists(parent_dir): os.mkdir(parent_dir) if not os.path.exists(dire): os.mkdir(dire) pcs_type_flow = 'GROUNDWATER_FLOW' var_name_flow = 'HEAD' t_id = 'transect' # ------------------------generate ogs base class----------------------------- # ogs = OGS(task_root=dire + "/", task_id=t_id, output_dir=dire + "/" + output + "/") # ------------------------ MSH -------------------------------- # # generate a rectangular mesh in x-y-plane ogs.msh.generate("rectangular", dim=2, mesh_origin=(0., 0.), element_no=(n_cellsx, n_cellsz), element_size=(1, 1)) # rotate mesh to obtain a cross section in x-z-plane ogs.msh.rotate(angle=np.pi / 2.0, rotation_axis=(1., 0., 0.)) # round nodes ogs.msh.NODES[:, 1] = np.around(ogs.msh.NODES[:, 1], 0) ogs.msh.NODES[:, 0] = np.around(ogs.msh.NODES[:, 0], 4) ogs.msh.NODES[:, 2] = np.around(ogs.msh.NODES[:, 2], 4)
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='hm_foot_tet_root', task_id='hm_foot_tet', output_dir='out', ) model.msh.read_file('hm_foot_tet.msh') model.gli.read_file('hm_foot_tet.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE=[ ['LIQUID_FLOW'], ['DEFORMATION'], ], MEMORY_TYPE=0, ST_RHS=2, ) model.rfd.read_file('hm_foot_tet.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION_FLOW', PRIMARY_VARIABLE='PRESSURE1', GEO_TYPE=['SURFACE', 'TOP'], DIS_TYPE=['CONSTANT', 0.0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION_FLOW', PRIMARY_VARIABLE='DISPLACEMENT_X1',
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='1d_isofrac_AS_root', task_id='1d_isofrac_AS', output_dir='out', ) model.msh.read_file('1d_isofrac_AS.msh') model.gli.read_file('1d_isofrac_AS.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='GROUNDWATER_FLOW', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MASS_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MASS_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='MASS_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block(
# Use the path to the recharge.txt as top comment in the .rfd-file rfd_top_com = recharge_path # Name the folder name = (str(overall_count) + "_b_" + str("{:0>4d}".format(thickness)) + "_Lc_" + str(c_length) + "_stor_" + str(storage) + "_rech_" + str(rech_abv)) # Name of directory (entire path) of one single ogs run dire = parent_dir + "/" + name # Make folders if not os.path.exists(dire): os.mkdir(dire) # -------------------- generate ogs base class ogs = OGS( task_root=dire + "/", task_id=t_id, output_dir=dire + "/" + "steady" + "/", ) # -------------------- MSH # Generate a rectangular mesh in x-z-plane. The first and the # second part of the mesh will be generated independently. # Afterwards, the two parts will be combined to one single # mesh. First part on the left side. This generates a mesh in # x-y-plane which will be rotated in a few steps. if thickness >= 100: # The size of cells in x and z direction (any float). s_cellsx = 10 s_cellsz = 10 # Number of cells is different for each model: n_cellsz = int(thickness / s_cellsx)
test_cnt += 1 relative_root, task_id = os.path.split(name) root = os.path.join(rootdir, relative_root) if use_old_dir: tmp_root = relative_root else: tmp_root = ("test_{:03}".format(test_cnt) + "_model") tmp_log = (str(test_cnt) + " " + tmp_root + " - load: " + os.path.join(root, task_id)) log_str.append(tmp_log) if print_progress: print(tmp_log) model = OGS(task_root=os.path.join(out_dir, tmp_root), output_dir="out") load_success = False try: model.load_model( verbose=False, task_root=root, task_id=task_id, encoding="ISO-8859-15", # lots of windows files (non utf8) use_task_id=True, # use the given task_id for some reference ) except ValueError as valerr: print("FAIL..." + str(valerr)) log_str.append("FAIL..." + str(valerr)) else: load_success = True print("OK...loading successful")
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='TDiff-Wall_root', task_id='TDiff-Wall', output_dir='out', ) model.msh.read_file('TDiff-Wall.msh') model.gli.read_file('TDiff-Wall.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='HEAT_TRANSPORT', TEMPERATURE_UNIT='KELVIN', ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='HEAT_TRANSPORT', PRIMARY_VARIABLE='TEMPERATURE1', GEO_TYPE=['POINT', 'POINT1'], DIS_TYPE=['CONSTANT', 25.0], ) model.ic.add_block( main_key='INITIAL_CONDITION', PCS_TYPE='HEAT_TRANSPORT', PRIMARY_VARIABLE='TEMPERATURE1', GEO_TYPE='DOMAIN', DIS_TYPE=['CONSTANT', 25], ) model.st.add_block( main_key='SOURCE_TERM',
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='tm1_3Dorigin_root', task_id='tm1_3Dorigin', output_dir='out', ) model.msh.read_file('tm1_3Dorigin.msh') model.gli.read_file('tm1_3Dorigin.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='HEAT_TRANSPORT', NUM_TYPE='NEW', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='DEFORMATION', NUM_TYPE='NEW', ) model.rfd.read_file('tm1_3Dorigin.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_X1', GEO_TYPE=['SURFACE', 'SURFACE1'], DIS_TYPE=[ ['LINEAR', 4], [0, 0.0], [1, 0.0], [2, 0.0],
# -*- coding: utf-8 -*- from ogs5py import OGS model = OGS( task_root='uc_creep01_root', task_id='uc_creep01', output_dir='out', ) model.msh.read_file('uc_creep01.msh') model.gli.read_file('uc_creep01.gli') model.pcs.add_block( main_key='PROCESS', PCS_TYPE='HEAT_TRANSPORT', ) model.pcs.add_block( main_key='PROCESS', PCS_TYPE='DEFORMATION', ) model.rfd.read_file('uc_creep01.rfd') model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_Y1', GEO_TYPE=['POLYLINE', 'PLY_2'], DIS_TYPE=['CONSTANT', 0], ) model.bc.add_block( main_key='BOUNDARY_CONDITION', PCS_TYPE='DEFORMATION', PRIMARY_VARIABLE='DISPLACEMENT_X1', GEO_TYPE=['POLYLINE', 'PLY_0'],