# Test goal: Sorption type comparison with simple 2D transport.
# Author: Pavel Exner
# Source for reference:
#   - qualitative, check in Paraview
#   - simple convective transport from right to left
#   - compare different sorption types as they slow down the convection of liquid phase
flow123d_version: 4.0.0a01
problem: !Coupling_Sequential
  description: Sorption type comparison with simple 2D transport.
  mesh:
    mesh_file: ../00_mesh/rectangle_5x1.msh
  flow_equation: !Flow_Darcy_LMH
    nonlinear_solver:
      linear_solver: !Petsc
        a_tol: 1.0e-07
    input_fields:
      - region: rectangle
        conductivity: 0.125
      - region: .BOUNDARY
        bc_type: dirichlet
        bc_pressure: !FieldFormula
          value: X[0]
    output:
      fields:
        - pressure_p0
        - velocity_p0
    output_stream:
      file: test_20.pvd
      format: !vtk
        variant: ascii
  solute_equation: !Coupling_OperatorSplitting
    transport: !Solute_Advection_FV
      input_fields:
        - region: ALL
          init_conc: 0
          porosity: 0.25
        - region: .BOUNDARY
          bc_conc: !FieldFormula
            value: 0.05*t
    substances:
      - A
      - As-lin
      - B
      - Bs-lang
      - C
      - Cs-freund
    output_stream:
      file: transport_test_20.pvd
      format: !vtk
        variant: ascii
      times:
        - step: 1.5
    time:
      end_time: 15.0
    reaction_term: !Sorption
      solvent_density: 1.0
      substances:
        - A
        - As-lin
        - Bs-lang
        - Cs-freund
      solubility:
        - 1.0
        - 1.0
        - 1.0
        - 1.0
      input_fields:
        - region: ALL
          rock_density: 1.0
          sorption_type:
            - none
            - linear
            - langmuir
            - freundlich
          distribution_coefficient:
            - 0
            - 0.6
            - 0.4
            - 0.6
          isotherm_other:
            - 0
            - 0
            - 0.6
            - 0.4
      output:
        fields: []
    balance:
      cumulative: true
