HM2a_VPF.prj

<OpenGeoSysProject>
  <meshes>
    <mesh>domain.vtu</mesh>
    <mesh>physical_group_DSS1.vtu</mesh>
    <mesh>physical_group_DSS2.vtu</mesh>
    <mesh>physical_group_DSS3.vtu</mesh>
    <mesh>physical_group_DSS4.vtu</mesh>
    <mesh>physical_group_DSS5.vtu</mesh>
    <mesh>physical_group_DSS6.vtu</mesh>
    <mesh>physical_group_DSS7.vtu</mesh>
    <mesh>physical_group_DSS8.vtu</mesh>
    <mesh>physical_group_DSS1a.vtu</mesh>
    <mesh>physical_group_DSS2a.vtu</mesh>
    <mesh>physical_group_DSS3a.vtu</mesh>
    <mesh>physical_group_DSS4a.vtu</mesh>
    <mesh>physical_group_DSS5a.vtu</mesh>
    <mesh>physical_group_DSS6a.vtu</mesh>
    <mesh>physical_group_DSS7a.vtu</mesh>
    <mesh>physical_group_DSS8a.vtu</mesh>
    <mesh>physical_group_PEE1.vtu</mesh>
    <mesh>physical_group_PEE2.vtu</mesh>
    <mesh>physical_group_PEE3.vtu</mesh>
    <mesh>physical_group_PEE4.vtu</mesh>
    <mesh>physical_group_PEE5.vtu</mesh>
    <mesh>physical_group_PEE6.vtu</mesh>
    <mesh>physical_group_PEE7.vtu</mesh>
    <mesh>physical_group_PEE8.vtu</mesh>
    <mesh>physical_group_PEE1a.vtu</mesh>
    <mesh>physical_group_PEE2a.vtu</mesh>
    <mesh>physical_group_PEE3a.vtu</mesh>
    <mesh>physical_group_PEE4a.vtu</mesh>
    <mesh>physical_group_PEE5a.vtu</mesh>
    <mesh>physical_group_PEE6a.vtu</mesh>
    <mesh>physical_group_PEE7a.vtu</mesh>
    <mesh>physical_group_PEE8a.vtu</mesh>
    <mesh>physical_group_p_left.vtu</mesh>
    <mesh>physical_group_p_right.vtu</mesh>
    <mesh>physical_group_p_top.vtu</mesh>
    <mesh>physical_group_p_bottom.vtu</mesh>
    <mesh>physical_group_Outlet_R_fullFracture.vtu</mesh>
    <mesh>physical_group_Outlet_L_fullFracture.vtu</mesh>
    <mesh>physical_group_Inlet.vtu</mesh>
  </meshes>
  <processes>
    <process>
      <name>HMPhaseField</name>
      <type>HM_PHASE_FIELD</type>
      <coupling_scheme>staggered</coupling_scheme>
      <spatial_stabilization_parameter>1.0</spatial_stabilization_parameter>
      <integration_order>2</integration_order>
      <phasefield_model>AT1</phasefield_model>
      <energy_split_model>VolumetricDeviatoric</energy_split_model>
      <irreversible_threshold>0.1</irreversible_threshold>
      <constitutive_relation id="0">
        <type>LinearElasticIsotropic</type>
        <youngs_modulus>E1</youngs_modulus>
        <poissons_ratio>nu1</poissons_ratio>
      </constitutive_relation>
      <constitutive_relation id="1">
        <type>LinearElasticIsotropic</type>
        <youngs_modulus>E1</youngs_modulus>
        <poissons_ratio>nu1</poissons_ratio>
      </constitutive_relation>
      <constitutive_relation id="2">
        <type>LinearElasticIsotropic</type>
        <youngs_modulus>E2</youngs_modulus>
        <poissons_ratio>nu2</poissons_ratio>
      </constitutive_relation>
      <phasefield_parameters>
        <residual_stiffness>k</residual_stiffness>
        <crack_resistance>gc</crack_resistance>
        <crack_length_scale>ls</crack_length_scale>
      </phasefield_parameters>
      <width_init>w_ini</width_init>
      <fracture_threshold>0.05</fracture_threshold>
      <fluid_compressibility>0.</fluid_compressibility>
      <specific_fracture_direction>0. 0.</specific_fracture_direction>
      <diffused_range_parameter>2</diffused_range_parameter>
      <fracture_permeability_parameter>10.</fracture_permeability_parameter>
      <process_variables>
        <phasefield>phasefield</phasefield>
        <pressure>pressure</pressure>
        <displacement>displacement</displacement>
      </process_variables>
      <secondary_variables>
        <secondary_variable internal_name="sigma" output_name="sigma" />
        <secondary_variable internal_name="epsilon" output_name="epsilon" />
        <secondary_variable internal_name="width_node" output_name="width_node" />
      </secondary_variables>
      <specific_body_force>0 0</specific_body_force>
    </process>
  </processes>
  <media>
    <medium id="0">
      <phases>
        <phase>
          <type>AqueousLiquid</type>
          <properties>
            <property>
              <name>viscosity</name>
              <type>Constant</type>
              <value>1e-3</value>
            </property>
            <property>
              <name>density</name>
              <type>Constant</type>
              <value>1e3</value>
            </property>
          </properties>
        </phase>
        <phase>
          <type>Solid</type>
          <properties>
            <property>
              <name>density</name>
              <type>Constant</type>
              <value>2.750e3</value>
            </property>
            <property>
              <name>biot_coefficient</name>
              <type>Constant</type>
              <value>0.8</value>
            </property>
          </properties>
        </phase>
      </phases>
      <properties>
        <property>
          <name>reference_temperature</name>
          <type>Constant</type>
          <value>293.15</value>
        </property>
        <property>
          <name>porosity</name>
          <type>Constant</type>
          <value>0.01</value>
        </property>
        <property>
          <name>permeability</name>
          <type>Constant</type>
          <value>1.e-19</value>
        </property>
      </properties>
    </medium>
    <medium id="1">
      <phases>
        <phase>
          <type>AqueousLiquid</type>
          <properties>
            <property>
              <name>viscosity</name>
              <type>Constant</type>
              <value>1e-3</value>
            </property>
            <property>
              <name>density</name>
              <type>Constant</type>
              <value>1e3</value>
            </property>
          </properties>
        </phase>
        <phase>
          <type>Solid</type>
          <properties>
            <property>
              <name>density</name>
              <type>Constant</type>
              <value>2.750e3</value>
            </property>
            <property>
              <name>biot_coefficient</name>
              <type>Constant</type>
              <value>0.8</value>
            </property>
          </properties>
        </phase>
      </phases>
      <properties>
        <property>
          <name>reference_temperature</name>
          <type>Constant</type>
          <value>293.15</value>
        </property>
        <property>
          <name>porosity</name>
          <type>Constant</type>
          <value>0.01</value>
        </property>
        <property>
          <name>permeability</name>
          <type>Constant</type>
          <value>1.e-19</value>
        </property>
      </properties>
    </medium>
    <medium id="2">
      <phases>
        <phase>
          <type>AqueousLiquid</type>
          <properties>
            <property>
              <name>viscosity</name>
              <type>Constant</type>
              <value>1e-3</value>
            </property>
            <property>
              <name>density</name>
              <type>Constant</type>
              <value>1e3</value>
            </property>
          </properties>
        </phase>
        <phase>
          <type>Solid</type>
          <properties>
            <property>
              <name>density</name>
              <type>Constant</type>
              <value>2.750e3</value>
            </property>
            <property>
              <name>biot_coefficient</name>
              <type>Constant</type>
              <value>0.01</value>
            </property>
          </properties>
        </phase>
      </phases>
      <properties>
        <property>
          <name>reference_temperature</name>
          <type>Constant</type>
          <value>293.15</value>
        </property>
        <property>
          <name>porosity</name>
          <type>Constant</type>
          <value>0.0001</value>
        </property>
        <property>
          <name>permeability</name>
          <type>Constant</type>
          <value>1.e-9</value>
        </property>
      </properties>
    </medium>
  </media>
  <time_loop>
    <global_process_coupling>
      <max_iter>100</max_iter>
      <convergence_criteria>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <abstol>1.e-12</abstol>
          <reltol>1.e-4</reltol>
        </convergence_criterion>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <abstol>1.e-12</abstol>
          <reltol>1.e-4</reltol>
        </convergence_criterion>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <abstol>1.e-12</abstol>
          <reltol>1.e-4</reltol>
        </convergence_criterion>
      </convergence_criteria>
      <local_coupling_processes>
        <max_iter>5</max_iter>
        <process_name>pressure</process_name>
        <process_name>displacement</process_name>
      </local_coupling_processes>
    </global_process_coupling>
    <processes>
      <process ref="HMPhaseField">
        <process_name>phasefield</process_name>
        <nonlinear_solver>petsc_snes</nonlinear_solver>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <abstol>1e-12</abstol>
          <reltol>1.e-6</reltol>
        </convergence_criterion>
        <time_discretization>
          <type>BackwardEuler</type>
        </time_discretization>
        <time_stepping>
          <type>IterationNumberBasedTimeStepping</type>
          <t_initial>0</t_initial>
          <t_end>3500</t_end>
          <initial_dt>2</initial_dt>
          <minimum_dt>1</minimum_dt>
          <maximum_dt>200</maximum_dt>
          <number_iterations>1 2 5 10 20</number_iterations>
          <multiplier>1.6 1.4 1.0 0.7 0.5</multiplier>
        </time_stepping>
      </process>
      <process ref="HMPhaseField">
        <process_name>pressure</process_name>
        <nonlinear_solver>basic_newton_p</nonlinear_solver>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <abstol>1e-12</abstol>
          <reltol>1.e-6</reltol>
        </convergence_criterion>
        <time_discretization>
          <type>BackwardEuler</type>
        </time_discretization>
        <time_stepping>
          <type>IterationNumberBasedTimeStepping</type>
          <t_initial>0</t_initial>
          <t_end>3500</t_end>
          <initial_dt>2</initial_dt>
          <minimum_dt>1</minimum_dt>
          <maximum_dt>200</maximum_dt>
          <number_iterations>1 2 5 10 20</number_iterations>
          <multiplier>1.6 1.4 1.0 0.7 0.5</multiplier>
        </time_stepping>
      </process>
      <process ref="HMPhaseField">
        <process_name>displacement</process_name>
        <nonlinear_solver>basic_newton_u</nonlinear_solver>
        <convergence_criterion>
          <type>DeltaX</type>
          <norm_type>NORM2</norm_type>
          <abstol>1e-10</abstol>
          <reltol>5.e-5</reltol>
        </convergence_criterion>
        <time_discretization>
          <type>BackwardEuler</type>
        </time_discretization>
        <time_stepping>
          <type>IterationNumberBasedTimeStepping</type>
          <t_initial>0</t_initial>
          <t_end>3500</t_end>
          <initial_dt>2</initial_dt>
          <minimum_dt>1</minimum_dt>
          <maximum_dt>200</maximum_dt>
          <number_iterations>1 2 5 10 20</number_iterations>
          <multiplier>1.6 1.4 1.0 0.7 0.5</multiplier>
        </time_stepping>
      </process>
    </processes>
    <output>
      <type>VTK</type>
      <prefix>GreatCell_2D</prefix>
      <variables />
      <output_iteration_results>false</output_iteration_results>
      <timesteps>
        <pair>
          <repeat>1</repeat>
          <each_steps>1</each_steps>
        </pair>
      </timesteps>
    </output>
  </time_loop>
  <parameters>
    <parameter>
      <name>E1</name>
      <type>Constant</type>
      <value>8.39e10</value>
    </parameter>
    <parameter>
      <name>nu1</name>
      <type>Constant</type>
      <value>0.21</value>
    </parameter>
    <parameter>
      <name>E2</name>
      <type>Constant</type>
      <value>1.e8</value>
    </parameter>
    <parameter>
      <name>nu2</name>
      <type>Constant</type>
      <value>0.4</value>
    </parameter>
    <parameter>
      <name>k</name>
      <type>Constant</type>
      <value>1e-10</value>
    </parameter>
    <parameter>
      <name>gc</name>
      <type>Constant</type>
      <value>10000000</value>
    </parameter>
    <parameter>
      <name>ls</name>
      <type>Constant</type>
      <value>0.001</value>
    </parameter>
    <parameter>
      <name>w_ini</name>
      <type>MeshElement</type>
      <field_name>width_ic</field_name>
    </parameter>
    <parameter>
      <name>injectionFlowRate_Inlet</name>
      <type>Constant</type>
      <value>8.334e-4</value>
    </parameter>
    <parameter>
      <name>injectionFlowRate_Inlet_ramp</name>
      <type>CurveScaled</type>
      <curve>timeRamp</curve>
      <parameter>injectionFlowRate_Inlet</parameter>
    </parameter>
    <parameter>
      <name>p_outlet</name>
      <type>Constant</type>
      <values>3.45e6</values>
    </parameter>
    <parameter>
      <name>pressure_ic</name>
      <type>Constant</type>
      <value>0.1e6</value>
    </parameter>
    <parameter>
      <name>displacement0</name>
      <type>Constant</type>
      <values>0 0</values>
    </parameter>
    <parameter>
      <name>phasefield_ic</name>
      <type>MeshNode</type>
      <field_name>pf-ic</field_name>
    </parameter>
    <parameter>
      <name>phasefield_bc</name>
      <type>Constant</type>
      <value>1.</value>
    </parameter>
    <parameter>
      <name>zero</name>
      <type>Constant</type>
      <value>0</value>
    </parameter>
    <parameter>
      <name>loading</name>
      <type>Constant</type>
      <value>1.0</value>
    </parameter>
    <parameter>
      <name>PEE1_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE1_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE2_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE2_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE3_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE3_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE5_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE5_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE4_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE4_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE6_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE6_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE7_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE7_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE8_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE8_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE1a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE1a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE2a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE2a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE3a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE3a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE4a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE4a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE5a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE5a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE6a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE6a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE7a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE7a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>PEE8a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>PEE8a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS1_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS1_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS2_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS2_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS3_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS3_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS5_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS5_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS4_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS4_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS6_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS6_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS7_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS7_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS8_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS8_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS1a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS1a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS2a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS2a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS3a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS3a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS4a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS4a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS5a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS5a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS6a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS6a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS7a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS7a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
    <parameter>
      <name>DSS8a_SURFACE_CURVE</name>
      <type>CurveScaled</type>
      <curve>DSS8a_SURFACE_CURVE</curve>
      <parameter>loading</parameter>
    </parameter>
  </parameters>
  <curves>
    <curve>
      <name>timeRamp</name>
      <coords>0 2999.9 3000 3500</coords>
      <values>0 0.0 1.0 1.0</values>
    </curve>
    <curve>
      <name>PEE1_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-7730000.0 -7730000.0</values>
    </curve>
    <curve>
      <name>PEE2_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-5700000.0 -5700000.0</values>
    </curve>
    <curve>
      <name>PEE3_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-4390000.0 -4390000.0</values>
    </curve>
    <curve>
      <name>PEE4_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-2400000.0 -2400000.0</values>
    </curve>
    <curve>
      <name>PEE5_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-2300000.0 -2300000.0</values>
    </curve>
    <curve>
      <name>PEE6_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-4000000.0 -4000000.0</values>
    </curve>
    <curve>
      <name>PEE7_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-6400000.0 -6400000.0</values>
    </curve>
    <curve>
      <name>PEE8_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-7700000.0 -7700000.0</values>
    </curve>
    <curve>
      <name>PEE1a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-7730000.0 -7730000.0</values>
    </curve>
    <curve>
      <name>PEE2a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-5700000.0 -5700000.0</values>
    </curve>
    <curve>
      <name>PEE3a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-4390000.0 -4390000.0</values>
    </curve>
    <curve>
      <name>PEE4a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-2400000.0 -2400000.0</values>
    </curve>
    <curve>
      <name>PEE5a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-2300000.0 -2300000.0</values>
    </curve>
    <curve>
      <name>PEE6a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-4000000.0 -4000000.0</values>
    </curve>
    <curve>
      <name>PEE7a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-6400000.0 -6400000.0</values>
    </curve>
    <curve>
      <name>PEE8a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-7700000.0 -7700000.0</values>
    </curve>
    <curve>
      <name>DSS1_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-6715000.0 -6715000.0</values>
    </curve>
    <curve>
      <name>DSS2_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-5045000.0 -5045000.0</values>
    </curve>
    <curve>
      <name>DSS3_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-3395000.0 -3395000.0</values>
    </curve>
    <curve>
      <name>DSS4_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-2350000.0 -2350000.0</values>
    </curve>
    <curve>
      <name>DSS5_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-3150000.0 -3150000.0</values>
    </curve>
    <curve>
      <name>DSS6_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-5200000.0 -5200000.0</values>
    </curve>
    <curve>
      <name>DSS7_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-7050000.0 -7050000.0</values>
    </curve>
    <curve>
      <name>DSS8_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-7715000.0 -7715000.0</values>
    </curve>
    <curve>
      <name>DSS1a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-6715000.0 -6715000.0</values>
    </curve>
    <curve>
      <name>DSS2a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-5045000.0 -5045000.0</values>
    </curve>
    <curve>
      <name>DSS3a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-3395000.0 -3395000.0</values>
    </curve>
    <curve>
      <name>DSS4a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-2350000.0 -2350000.0</values>
    </curve>
    <curve>
      <name>DSS5a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-3150000.0 -3150000.0</values>
    </curve>
    <curve>
      <name>DSS6a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-5200000.0 -5200000.0</values>
    </curve>
    <curve>
      <name>DSS7a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-7050000.0 -7050000.0</values>
    </curve>
    <curve>
      <name>DSS8a_SURFACE_CURVE</name>
      <coords>0 500</coords>
      <values>-7715000.0 -7715000.0</values>
    </curve>
  </curves>
  <process_variables>
    <process_variable>
      <name>phasefield</name>
      <components>1</components>
      <order>1</order>
      <initial_condition>phasefield_ic</initial_condition>
      <boundary_conditions>
        <boundary_condition>
          <mesh>physical_group_PEE1</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE2</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE3</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE4</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE5</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE6</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE7</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE8</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE1a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE2a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE3a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE4a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE5a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE6a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE7a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE8a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS1</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS2</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS3</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS4</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS5</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS6</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS7</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS8</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS1a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS2a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS3a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS4a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS5a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS6a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS7a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS8a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>phasefield_bc</parameter>
        </boundary_condition>
      </boundary_conditions>
    </process_variable>
    <process_variable>
      <name>pressure</name>
      <components>1</components>
      <order>1</order>
      <initial_condition>pressure_ic</initial_condition>
      <source_terms>
        <source_term>
          <mesh>physical_group_Inlet</mesh>
          <type>Line</type>
          <parameter>injectionFlowRate_Inlet_ramp</parameter>
        </source_term>
      </source_terms>
      <boundary_conditions>
        <boundary_condition>
          <mesh>physical_group_Outlet_R_fullFracture</mesh>
          <type>DirichletWithinTimeInterval</type>
          <parameter>p_outlet</parameter>
          <time_interval>
            <start>0</start>
            <end>3500</end>
          </time_interval>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_Outlet_L_fullFracture</mesh>
          <type>DirichletWithinTimeInterval</type>
          <parameter>p_outlet</parameter>
          <time_interval>
            <start>0</start>
            <end>3500</end>
          </time_interval>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE1</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE2</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE3</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE4</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE5</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE6</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE7</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE8</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE1a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE2a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE3a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE4a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE5a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE6a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE7a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE8a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS1</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS2</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS3</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS4</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS5</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS6</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS7</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS8</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS1a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS2a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS3a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS4a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS5a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS6a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS7a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS8a</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>pressure_ic</parameter>
        </boundary_condition>
      </boundary_conditions>
    </process_variable>
    <process_variable>
      <name>displacement</name>
      <components>2</components>
      <order>1</order>
      <initial_condition>displacement0</initial_condition>
      <boundary_conditions>
        <boundary_condition>
          <mesh>physical_group_p_bottom</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>zero</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_p_bottom</mesh>
          <type>Dirichlet</type>
          <component>1</component>
          <parameter>zero</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_p_left</mesh>
          <type>Dirichlet</type>
          <component>0</component>
          <parameter>zero</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE1</mesh>
          <type>NormalTraction</type>
          <parameter>PEE1_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE2</mesh>
          <type>NormalTraction</type>
          <parameter>PEE2_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE3</mesh>
          <type>NormalTraction</type>
          <parameter>PEE3_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE4</mesh>
          <type>NormalTraction</type>
          <parameter>PEE4_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE5</mesh>
          <type>NormalTraction</type>
          <parameter>PEE5_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE6</mesh>
          <type>NormalTraction</type>
          <parameter>PEE6_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE7</mesh>
          <type>NormalTraction</type>
          <parameter>PEE7_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE8</mesh>
          <type>NormalTraction</type>
          <parameter>PEE8_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE1a</mesh>
          <type>NormalTraction</type>
          <parameter>PEE1a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE2a</mesh>
          <type>NormalTraction</type>
          <parameter>PEE2a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE3a</mesh>
          <type>NormalTraction</type>
          <parameter>PEE3a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE4a</mesh>
          <type>NormalTraction</type>
          <parameter>PEE4a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE5a</mesh>
          <type>NormalTraction</type>
          <parameter>PEE5a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE6a</mesh>
          <type>NormalTraction</type>
          <parameter>PEE6a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE7a</mesh>
          <type>NormalTraction</type>
          <parameter>PEE7a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_PEE8a</mesh>
          <type>NormalTraction</type>
          <parameter>PEE8a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS1</mesh>
          <type>NormalTraction</type>
          <parameter>DSS1_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS2</mesh>
          <type>NormalTraction</type>
          <parameter>DSS2_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS3</mesh>
          <type>NormalTraction</type>
          <parameter>DSS3_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS4</mesh>
          <type>NormalTraction</type>
          <parameter>DSS4_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS5</mesh>
          <type>NormalTraction</type>
          <parameter>DSS5_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS6</mesh>
          <type>NormalTraction</type>
          <parameter>DSS6_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS7</mesh>
          <type>NormalTraction</type>
          <parameter>DSS7_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS8</mesh>
          <type>NormalTraction</type>
          <parameter>DSS8_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS1a</mesh>
          <type>NormalTraction</type>
          <parameter>DSS1a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS2a</mesh>
          <type>NormalTraction</type>
          <parameter>DSS2a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS3a</mesh>
          <type>NormalTraction</type>
          <parameter>DSS3a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS4a</mesh>
          <type>NormalTraction</type>
          <parameter>DSS4a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS5a</mesh>
          <type>NormalTraction</type>
          <parameter>DSS5a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS6a</mesh>
          <type>NormalTraction</type>
          <parameter>DSS6a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS7a</mesh>
          <type>NormalTraction</type>
          <parameter>DSS7a_SURFACE_CURVE</parameter>
        </boundary_condition>
        <boundary_condition>
          <mesh>physical_group_DSS8a</mesh>
          <type>NormalTraction</type>
          <parameter>DSS8a_SURFACE_CURVE</parameter>
        </boundary_condition>
      </boundary_conditions>
    </process_variable>
  </process_variables>
  <nonlinear_solvers>
    <nonlinear_solver>
      <name>petsc_snes</name>
      <type>PETScSNES</type>
      <max_iter>50</max_iter>
      <linear_solver>linear_solver_d</linear_solver>
    </nonlinear_solver>
    <nonlinear_solver>
      <name>basic_newton_d</name>
      <type>Newton</type>
      <max_iter>50</max_iter>
      <linear_solver>linear_solver_d</linear_solver>
    </nonlinear_solver>
    <nonlinear_solver>
      <name>basic_newton_u</name>
      <type>Newton</type>
      <max_iter>659</max_iter>
      <linear_solver>linear_solver_u</linear_solver>
    </nonlinear_solver>
    <nonlinear_solver>
      <name>basic_newton_p</name>
      <type>Newton</type>
      <max_iter>500</max_iter>
      <linear_solver>linear_solver_p</linear_solver>
    </nonlinear_solver>
  </nonlinear_solvers>
  <linear_solvers>
    <linear_solver>
      <name>linear_solver_d</name>
      <eigen>
        <solver_type>BiCGSTAB</solver_type>
        <precon_type>ILUT</precon_type>
        <max_iteration_step>100</max_iteration_step>
        <error_tolerance>1e-16</error_tolerance>
      </eigen>
      <petsc>
        <parameters>-ksp_type bcgs -pc_type bjacobi -ksp_atol 1e-10 -ksp_rtol 1e-10 -snes_type vinewtonrsls -snes_linesearch_type l2 -snes_atol 1.e-5 -snes_rtol 1.e-5 -snes_max_it 1000 -snes_monitor</parameters>
      </petsc>
    </linear_solver>
    <linear_solver>
      <name>linear_solver_u</name>
      <eigen>
        <solver_type>BiCGSTAB</solver_type>
        <precon_type>ILUT</precon_type>
        <max_iteration_step>10000</max_iteration_step>
        <error_tolerance>1e-16</error_tolerance>
      </eigen>
      <petsc>
        <prefix>u</prefix>
        <parameters>-u_ksp_type bcgs -u_pc_type hypre -u_pc_hypre_type boomeramg -u_pc_hypre_boomeramg_strong_threshold 0.7 -u_ksp_atol 1e-12 -u_ksp_rtol 1e-12</parameters>
      </petsc>
    </linear_solver>
    <linear_solver>
      <name>linear_solver_p</name>
      <eigen>
        <solver_type>BiCGSTAB</solver_type>
        <precon_type>ILUT</precon_type>
        <max_iteration_step>10000</max_iteration_step>
        <error_tolerance>1e-16</error_tolerance>
      </eigen>
      <petsc>
        <prefix>p</prefix>
        <parameters>-p_ksp_type gmres
                        -p_ksp_gmres_restart 100
                        -p_pc_type hypre
                        -p_pc_hypre_type boomeramg
                        -p_pc_hypre_boomeramg_strong_threshold 0.25
                        -p_pc_hypre_boomeramg_max_iter 100
                        -p_ksp_atol 1e-12
                        -p_ksp_rtol 1e-8
                        -p_ksp_max_it 10000</parameters>
      </petsc>
    </linear_solver>
  </linear_solvers>
</OpenGeoSysProject>