OGS
th_decovalex.prj
<OpenGeoSysProject>
  <mesh>th_decovalex.vtu</mesh>
  <geometry>boundary.gml</geometry>
  <processes>
    <process>
      <name>th_decovalex</name>
      <type>HT</type>
      <coupling_scheme>staggered</coupling_scheme>
      <integration_order>2</integration_order>
      <process_variables>
        <temperature>T</temperature>
        <pressure>p</pressure>
      </process_variables>
      <specific_body_force>0 -9.81</specific_body_force>
      <secondary_variables>
        <secondary_variable internal_name="darcy_velocity" output_name="darcy_velocity" />
      </secondary_variables>
    </process>
  </processes>
  <media>
    <medium id="0">
      <phases>
        <phase>
          <type>AqueousLiquid</type>
          <properties>
            <property>
              <name>density</name>
              <type>Linear</type>
              <reference_value>1000.0</reference_value>
              <independent_variable>
                <variable_name>temperature</variable_name>
                <reference_condition>293.0</reference_condition>
                <slope>-4.3e-4</slope>
              </independent_variable>
            </property>
            <property>
              <name>viscosity</name>
              <type>Constant</type>
              <value>1.0e-3</value>
            </property>
            <property>
              <name>specific_heat_capacity</name>
              <type>Constant</type>
              <value>4000</value>
            </property>
            <property>
              <name>thermal_conductivity</name>
              <type>Constant</type>
              <value>18921600.0</value>
            </property>
          </properties>
        </phase>
        <phase>
          <type>Solid</type>
          <properties>
            <property>
              <name>storage</name>
              <type>Constant</type>
              <value>0.0</value>
            </property>
            <property>
              <name>density</name>
              <type>Parameter</type>
              <parameter_name>rho_solid</parameter_name>
            </property>
            <property>
              <name>thermal_conductivity</name>
              <type>Parameter</type>
              <parameter_name>lambda</parameter_name>
            </property>
            <property>
              <name>specific_heat_capacity</name>
              <type>Parameter</type>
              <parameter_name>c_p</parameter_name>
            </property>
          </properties>
        </phase>
      </phases>
      <properties>
        <property>
          <name>thermal_longitudinal_dispersivity</name>
          <type>Constant</type>
          <value>0.0</value>
        </property>
        <property>
          <name>thermal_transversal_dispersivity</name>
          <type>Constant</type>
          <value>0.0</value>
        </property>
        <property>
          <name>permeability</name>
          <type>Constant</type>
          <value>6.3072e-14 0. 0. 6.3072e-14</value>
        </property>
        <property>
          <name>porosity</name>
          <type>Constant</type>
          <value>0.01</value>
        </property>
        <property>
          <name>thermal_conductivity</name>
          <type>EffectiveThermalConductivityPorosityMixing</type>
        </property>
      </properties>
    </medium>
    <medium id="1">
      <phases>
        <phase>
          <type>AqueousLiquid</type>
          <properties>
            <property>
              <name>density</name>
              <type>Linear</type>
              <reference_value>1000.0</reference_value>
              <independent_variable>
                <variable_name>temperature</variable_name>
                <reference_condition>293.0</reference_condition>
                <slope>-4.3e-4</slope>
              </independent_variable>
            </property>
            <property>
              <name>viscosity</name>
              <type>Constant</type>
              <value>1.0e-3</value>
            </property>
            <property>
              <name>specific_heat_capacity</name>
              <type>Constant</type>
              <value>4000</value>
            </property>
            <property>
              <name>thermal_conductivity</name>
              <type>Constant</type>
              <value>18921600.0</value>
            </property>
          </properties>
        </phase>
        <phase>
          <type>Solid</type>
          <properties>
            <property>
              <name>density</name>
              <type>Parameter</type>
              <parameter_name>rho_solid</parameter_name>
            </property>
            <property>
              <name>storage</name>
              <type>Constant</type>
              <value>0.0</value>
            </property>
            <property>
              <name>thermal_conductivity</name>
              <type>Parameter</type>
              <parameter_name>lambda</parameter_name>
            </property>
            <property>
              <name>specific_heat_capacity</name>
              <type>Parameter</type>
              <parameter_name>c_p</parameter_name>
            </property>
          </properties>
        </phase>
      </phases>
      <properties>
        <property>
          <name>thermal_longitudinal_dispersivity</name>
          <type>Constant</type>
          <value>0.0</value>
        </property>
        <property>
          <name>thermal_transversal_dispersivity</name>
          <type>Constant</type>
          <value>0.0</value>
        </property>
        <property>
          <name>specific_heat_capacity</name>
          <type>Constant</type>
          <value>4000</value>
        </property>
        <property>
          <name>permeability</name>
          <type>Constant</type>
          <value>3.1536e-10 0. 0. 3.1536e-10</value>
        </property>
        <property>
          <name>porosity</name>
          <type>Constant</type>
          <value>0.41</value>
        </property>
        <property>
          <name>thermal_conductivity</name>
          <type>EffectiveThermalConductivityPorosityMixing</type>
        </property>
      </properties>
    </medium>
  </media>
  <time_loop>
    <global_process_coupling>
      <max_iter>6</max_iter>
      <convergence_criteria>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <reltol>1.e-14</reltol>
        </convergence_criterion>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <reltol>1.e-14</reltol>
        </convergence_criterion>
      </convergence_criteria>
    </global_process_coupling>
    <processes>
      <process ref="th_decovalex">
        <nonlinear_solver>basic_picard_T</nonlinear_solver>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <abstol>1.e-6</abstol>
        </convergence_criterion>
        <time_discretization>
          <type>BackwardEuler</type>
        </time_discretization>
        <time_stepping>
          <type>FixedTimeStepping</type>
          <t_initial>0.0</t_initial>
          <t_end>1e3</t_end>
          <timesteps>
            <pair>
              <repeat>4</repeat>
              <delta_t>0.001</delta_t>
            </pair>
            <pair>
              <repeat>1</repeat>
              <delta_t>0.006</delta_t>
            </pair>
            <pair>
              <repeat>9</repeat>
              <delta_t>0.01</delta_t>
            </pair>
            <pair>
              <repeat>9</repeat>
              <delta_t>0.1</delta_t>
            </pair>
            <pair>
              <repeat>39</repeat>
              <delta_t>1.</delta_t>
            </pair>
            <pair>
              <repeat>2</repeat>
              <delta_t>2.0</delta_t>
            </pair>
            <pair>
              <repeat>1</repeat>
              <delta_t>16.0</delta_t>
            </pair>
            <pair>
              <repeat>1</repeat>
              <delta_t>30.0</delta_t>
            </pair>
            <pair>
              <repeat>1</repeat>
              <delta_t>40.0</delta_t>
            </pair>
            <pair>
              <repeat>10</repeat>
              <delta_t>86.0</delta_t>
            </pair>
            <pair>
              <repeat>90</repeat>
              <delta_t>100.0</delta_t>
            </pair>
            <pair>
              <repeat>100</repeat>
              <delta_t>900</delta_t>
            </pair>
            <pair>
              <repeat>90</repeat>
              <delta_t>1.e4</delta_t>
            </pair>
          </timesteps>
        </time_stepping>
      </process>
      <process ref="th_decovalex">
        <nonlinear_solver>basic_picard_H</nonlinear_solver>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <abstol>1.e-6</abstol>
        </convergence_criterion>
        <time_discretization>
          <type>BackwardEuler</type>
        </time_discretization>
        <time_stepping>
          <type>FixedTimeStepping</type>
          <t_initial>0.0</t_initial>
          <t_end>1.e3</t_end>
          <timesteps>
            <pair>
              <repeat>4</repeat>
              <delta_t>0.001</delta_t>
            </pair>
            <pair>
              <repeat>1</repeat>
              <delta_t>0.006</delta_t>
            </pair>
            <pair>
              <repeat>9</repeat>
              <delta_t>0.01</delta_t>
            </pair>
            <pair>
              <repeat>9</repeat>
              <delta_t>0.1</delta_t>
            </pair>
            <pair>
              <repeat>39</repeat>
              <delta_t>1.</delta_t>
            </pair>
            <pair>
              <repeat>2</repeat>
              <delta_t>2.0</delta_t>
            </pair>
            <pair>
              <repeat>1</repeat>
              <delta_t>16.0</delta_t>
            </pair>
            <pair>
              <repeat>1</repeat>
              <delta_t>30.0</delta_t>
            </pair>
            <pair>
              <repeat>1</repeat>
              <delta_t>40.0</delta_t>
            </pair>
            <pair>
              <repeat>10</repeat>
              <delta_t>86.0</delta_t>
            </pair>
            <pair>
              <repeat>90</repeat>
              <delta_t>100.0</delta_t>
            </pair>
            <pair>
              <repeat>100</repeat>
              <delta_t>900</delta_t>
            </pair>
            <pair>
              <repeat>90</repeat>
              <delta_t>1.e4</delta_t>
            </pair>
          </timesteps>
        </time_stepping>
      </process>
    </processes>
    <output>
      <type>VTK</type>
      <prefix>th_decovalex</prefix>
      <suffix>ts{:timestep}_t_{:time}</suffix>
      <timesteps>
        <pair>
          <repeat>1</repeat>
          <each_steps>10</each_steps>
        </pair>
      </timesteps>
      <variables>
        <variable>T</variable>
        <variable>p</variable>
        <variable>darcy_velocity</variable>
      </variables>
    </output>
  </time_loop>
  <parameters>
    <parameter>
      <name>p0</name>
      <type>Constant</type>
      <value>4.7e6</value>
    </parameter>
    <parameter>
      <name>p_bc_top</name>
      <type>Constant</type>
      <value>4.3e+6</value>
    </parameter>
    <parameter>
      <name>p_bc_bottom</name>
      <type>Constant</type>
      <value>4.7e+6</value>
    </parameter>
    <parameter>
      <name>lambda</name>
      <type>MeshElement</type>
      <field_name>scaled_solid_thermal_conductivity</field_name>
    </parameter>
    <parameter>
      <name>lambda_fluid</name>
      <type>Constant</type>
      <value>18921600.0</value>
    </parameter>
    <parameter>
      <name>c_p</name>
      <type>MeshElement</type>
      <field_name>solid_heat_capacity</field_name>
    </parameter>
    <parameter>
      <name>rho_solid</name>
      <type>MeshElement</type>
      <field_name>solid_density</field_name>
    </parameter>
    <parameter>
      <name>n0</name>
      <type>Constant</type>
      <value>0.01</value>
    </parameter>
    <parameter>
      <name>n1</name>
      <type>Constant</type>
      <value>0.41</value>
    </parameter>
    <parameter>
      <name>T0</name>
      <type>Constant</type>
      <value>298</value>
    </parameter>
    <parameter>
      <name>TB_top</name>
      <type>Constant</type>
      <value>294.0</value>
    </parameter>
    <parameter>
      <name>TB_bottom</name>
      <type>Constant</type>
      <value>310</value>
    </parameter>
    <parameter>
      <name>mu0</name>
      <type>Constant</type>
      <values>6.3072e-14 0. 0. 6.3072e-14</values>
    </parameter>
    <parameter>
      <name>mu1</name>
      <type>Constant</type>
      <values>3.1536e-10 0. 0. 3.1536e-10</values>
    </parameter>
    <parameter>
      <name>hpower_on_arc</name>
      <type>Constant</type>
      <value>1.0</value>
    </parameter>
    <parameter>
      <name>heat_power</name>
      <type>CurveScaled</type>
      <curve>time_heat_power</curve>
      <parameter>hpower_on_arc</parameter>
    </parameter>
    <parameter>
      <name>heat_power_c</name>
      <type>Constant</type>
      <value>2477507279.46192</value>
    </parameter>
  </parameters>
  <curves>
    <curve>
      <name>time_heat_power</name>
      <coords>0 1.00E-03 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 2 3 4 5 6 7 8 9 10 20 30 40 50 60 70 80 90 100 200 300 400 500 600 700 800 900 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 20000 30000 40000 50000 60000 70000 80000 90000 100000 200000 300000 400000 500000 600000 700000 800000 900000 1000000 1.00E+10</coords>
      <values>2477507279.46192 2477459529.54288 2477411689.43088 2477363767.32528 2477315670.51024 2477267748.40464 2477219911.13088 2477171729.484 2477123892.21024 2477075710.56336 2477027873.2896 2472672907.16208 2467757091.56544 2462762401.59456 2457693385.056 2452552352.27712 2447343250.9344 2442069507.41424 2436734637.66528 2431341468.88992 2425894034.2272 2369163847.82544 2310711582.0096 2253930470.64432 2202214690.836 2158958419.37568 2112431166.82368 2078620739.1 2045550770.03376 2013204028.0392 1726032720.05136 1495088443.47744 1308713237.89344 1157694104.67264 1034749222.15824 934122524.20848 851262379.89264 782568686.43936 725193373.44672 448276352.8464 341492110.56144 277947936.256176 236509762.562928 209106613.956816 190948269.79944 178912492.01208 154295944.435872 143838731.939184 87693573.405456 68334737.77152 58916070.8136 52907269.428288 48548275.775136 45238544.066592 42688989.27792 41088819.808944 36699991.77528 21249719.02512 14044626.9792 9640538.832816 6944041.074528 5293033.113744 4282092.927936 3663114.947568 3284183.858832 3052109.993328 1972070.615664 1660867.028928 1472515.998336 1358407.99656 1289308.804704 1247386.360032 1222095.6864 1206749.827296 1197405.142848 1197405.142848</values>
    </curve>
  </curves>
  <process_variables>
    <process_variable>
      <name>T</name>
      <components>1</components>
      <order>1</order>
      <initial_condition>T0</initial_condition>
      <boundary_conditions>
        <boundary_condition>
          <geometrical_set>th_decovalex</geometrical_set>
          <geometry>top</geometry>
          <type>Dirichlet</type>
          <parameter>TB_top</parameter>
        </boundary_condition>
        <boundary_condition>
          <geometrical_set>th_decovalex</geometrical_set>
          <geometry>bottom</geometry>
          <type>Dirichlet</type>
          <parameter>TB_bottom</parameter>
        </boundary_condition>
        <boundary_condition>
          <geometrical_set>th_decovalex</geometrical_set>
          <geometry>arc</geometry>
          <type>Neumann</type>
          <parameter>heat_power</parameter>
        </boundary_condition>
      </boundary_conditions>
    </process_variable>
    <process_variable>
      <name>p</name>
      <components>1</components>
      <order>1</order>
      <initial_condition>p0</initial_condition>
      <boundary_conditions>
        <boundary_condition>
          <geometrical_set>th_decovalex</geometrical_set>
          <geometry>top</geometry>
          <type>Dirichlet</type>
          <parameter>p_bc_top</parameter>
        </boundary_condition>
        <boundary_condition>
          <geometrical_set>th_decovalex</geometrical_set>
          <geometry>bottom</geometry>
          <type>Dirichlet</type>
          <parameter>p_bc_bottom</parameter>
        </boundary_condition>
      </boundary_conditions>
    </process_variable>
  </process_variables>
  <nonlinear_solvers>
    <nonlinear_solver>
      <name>basic_picard_T</name>
      <type>Picard</type>
      <max_iter>10</max_iter>
      <linear_solver>linear_solver_T</linear_solver>
    </nonlinear_solver>
    <nonlinear_solver>
      <name>basic_picard_H</name>
      <type>Picard</type>
      <max_iter>10</max_iter>
      <linear_solver>linear_solver_H</linear_solver>
    </nonlinear_solver>
  </nonlinear_solvers>
  <linear_solvers>
    <linear_solver>
      <name>linear_solver_T</name>
      <lis>-i cg -p jacobi -tol 1e-20 -maxiter 10000</lis>
      <eigen>
        <solver_type>BiCGSTAB</solver_type>
        <precon_type>DIAGONAL</precon_type>
        <max_iteration_step>10000</max_iteration_step>
        <error_tolerance>1e-20</error_tolerance>
      </eigen>
      <petsc>
        <prefix>T</prefix>
        <parameters>-T_ksp_type bcgs -T_pc_type bjacobi -T_ksp_rtol 1e-16 -T_ksp_max_it 10000</parameters>
      </petsc>
    </linear_solver>
    <linear_solver>
      <name>linear_solver_H</name>
      <lis>-i cg -p jacobi -tol 1e-20 -maxiter 10000</lis>
      <eigen>
        <solver_type>CG</solver_type>
        <precon_type>DIAGONAL</precon_type>
        <max_iteration_step>10000</max_iteration_step>
        <error_tolerance>1e-20</error_tolerance>
      </eigen>
      <petsc>
        <prefix>H</prefix>
        <parameters>-H_ksp_type cg -H_pc_type bjacobi -H_ksp_rtol 1e-16 -H_ksp_max_it 10000</parameters>
      </petsc>
    </linear_solver>
  </linear_solvers>
</OpenGeoSysProject>