|  | OGS
    | 
 <OpenGeoSysProject>
   <meshes>
     <mesh>mixed_dimensional_grid_2_updated.vtu</mesh>
     <mesh>x_inlet.vtu</mesh>
     <mesh>x_outlet.vtu</mesh>
     <mesh>borehole_left_Layer11.vtu</mesh>
     <mesh>borehole_right_Layer10.vtu</mesh>
   </meshes>
   <processes>
     <process>
       <name>hc</name>
       <type>ComponentTransport</type>
       <integration_order>2</integration_order>
       <process_variables>
         <concentration>Si</concentration>
         <pressure>pressure</pressure>
       </process_variables>
       <specific_body_force>0 0 -9.81</specific_body_force>
       <secondary_variables>
         <secondary_variable name="darcy_velocity" />
       </secondary_variables>
       <aperture_size>
         <parameter>fracture_thickness_const</parameter>
       </aperture_size>
     </process>
   </processes>
   <media>
     <medium id="*">
       <phases>
         <phase>
           <type>AqueousLiquid</type>
           <components>
             <component>
               <name>Si</name>
               <properties>
                 <property>
                   <name>pore_diffusion</name>
                   <type>Constant</type>
                   <value>1e-9</value>
                 </property>
                 <property>
                   <name>retardation_factor</name>
                   <type>Constant</type>
                   <value>1</value>
                 </property>
                 <property>
                   <name>decay_rate</name>
                   <type>Parameter</type>
                   <parameter_name>decay</parameter_name>
                 </property>
               </properties>
             </component>
           </components>
           <properties>
             <property>
               <name>density</name>
               <type>Constant</type>
               <value>1e3</value>
             </property>
             <property>
               <name>viscosity</name>
               <type>Constant</type>
               <value>1e-3</value>
             </property>
           </properties>
         </phase>
       </phases>
       <properties>
         <property>
           <name>permeability</name>
           <type>Parameter</type>
           <parameter_name>kappa1_frac</parameter_name>
         </property>
         <property>
           <name>porosity</name>
           <type>Parameter</type>
           <parameter_name>constant_porosity_parameter</parameter_name>
         </property>
         <property>
           <name>longitudinal_dispersivity</name>
           <type>Constant</type>
           <value>10</value>
         </property>
         <property>
           <name>transversal_dispersivity</name>
           <type>Constant</type>
           <value>1</value>
         </property>
       </properties>
     </medium>
   </media>
   <time_loop>
     <processes>
       <process ref="hc">
         <nonlinear_solver>basic_picard</nonlinear_solver>
         <convergence_criterion>
           <type>PerComponentDeltaX</type>
           <norm_type>NORM2</norm_type>
           <reltols>1e-12 1e-12</reltols>
         </convergence_criterion>
         <time_discretization>
           <type>BackwardEuler</type>
         </time_discretization>
         <time_stepping>
           <type>IterationNumberBasedTimeStepping</type>
           <t_initial>0.0</t_initial>
           <t_end>864000000</t_end>
           <initial_dt>86400</initial_dt>
           <minimum_dt>10</minimum_dt>
           <maximum_dt>8.64e10</maximum_dt>
           <number_iterations>1 7 11 15</number_iterations>
           <multiplier>1.4 1.2 0.9 0.5</multiplier>
         </time_stepping>
       </process>
     </processes>
     <output>
       <type>VTK</type>
       <prefix>DFN_HC</prefix>
       <timesteps>
         <pair>
           <repeat>1000</repeat>
           <each_steps>1</each_steps>
         </pair>
       </timesteps>
       <variables />
       <suffix>ts{:timestep}_t_{:time}</suffix>
     </output>
   </time_loop>
   <parameters>
     <parameter>
       <name>fracture_thickness_const</name>
       <type>MeshElement</type>
       <field_name>width_ic</field_name>
     </parameter>
     <parameter>
       <name>kappa1_frac</name>
       <type>MeshElement</type>
       <field_name>permeability_ic</field_name>
     </parameter>
     <parameter>
       <name>constant_porosity_parameter</name>
       <type>Constant</type>
       <value>1</value>
     </parameter>
     <parameter>
       <name>decay</name>
       <type>Constant</type>
       <value>1e-10</value>
     </parameter>
     <parameter>
       <name>c0</name>
       <type>Constant</type>
       <value>0.0</value>
     </parameter>
     <parameter>
       <name>p0</name>
       <type>Function</type>
       <expression>1000*9.810*(500-z)</expression>
     </parameter>
     <parameter>
       <name>pinlet</name>
       <type>Function</type>
       <expression>1000*9.810*(500-z)+294300</expression>
     </parameter>
     <parameter>
       <name>poutlet</name>
       <type>Function</type>
       <expression>1000*9.810*(500-z)</expression>
     </parameter>
     <parameter>
       <name>c_bottom</name>
       <type>Constant</type>
       <value>1</value>
     </parameter>
   </parameters>
   <process_variables>
     <process_variable>
       <name>Si</name>
       <components>1</components>
       <order>1</order>
       <initial_condition>c0</initial_condition>
       <boundary_conditions>
         <boundary_condition>
           <type>Dirichlet</type>
           <mesh>borehole_right_Layer10</mesh>
           <parameter>c_bottom</parameter>
         </boundary_condition>
       </boundary_conditions>
     </process_variable>
     <process_variable>
       <name>pressure</name>
       <components>1</components>
       <order>1</order>
       <initial_condition>p0</initial_condition>
       <boundary_conditions>
         <boundary_condition>
           <type>Dirichlet</type>
           <mesh>x_outlet</mesh>
           <parameter>poutlet</parameter>
         </boundary_condition>
         <boundary_condition>
           <type>Dirichlet</type>
           <mesh>x_inlet</mesh>
           <parameter>pinlet</parameter>
         </boundary_condition>
       </boundary_conditions>
     </process_variable>
   </process_variables>
   <nonlinear_solvers>
     <nonlinear_solver>
       <name>basic_picard</name>
       <type>Picard</type>
       <max_iter>30</max_iter>
       <linear_solver>general_linear_solver</linear_solver>
     </nonlinear_solver>
   </nonlinear_solvers>
   <linear_solvers>
     <linear_solver>
       <name>general_linear_solver</name>
       <eigen>
         <solver_type>SparseLU</solver_type>
       </eigen>
     </linear_solver>
   </linear_solvers>
 </OpenGeoSysProject>