<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 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>