multi_solute

    Test Multi-Solute Transport with Chemical Reaction

    The test case is constructed from the VV Test Suite multi_solute model problem. It uses the trac and trxn macros and is a comparison of FEHM and PDREACT Solution for Concentration vs Time. Test compares the generated tracer files to old tracer files known to be correct. All concentraction values are tested.

    Test Cases for Multi-Solute Transport with Chemical Reaction - trac macro

    Species Co_aq
    Species Fe_aq
    Species EDTA_aq
    Species CoEDTA_aq
    Species FeEDTA_aq
    Species CoEDTA_s
    Species Co_s
    Species FeEDTA_s

    Multi-Solute Transport with Chemical Reaction - trxn macro

    Species Co_aq
    Species Fe_aq
    Species EDTA_aq
    Species CoEDTA_aq
    Species FeEDTA_aq
    Species CoEDTA_s
    Species Co_s
    Species FeEDTA_s

    Comparison files:

    multi_solute_trac_Cobalt[aq].trc	       multi_solute_trxn_Cobalt_s.trc
    multi_solute_trac_Cobalt[s].trc            multi_solute_trxn_Cobalt.trc
    multi_solute_trac_Co-EDTA[aq].trc	       multi_solute_trxn_Co-EDTA_a.trc
    multi_solute_trac_Co-EDTA[s].trc	       multi_solute_trxn_Co-EDTA_s.trc
    multi_solute_trac_EDTA[aq].trc             multi_solute_trxn_EDTA.trc
    multi_solute_trac_Fe-EDTA[aq].trc          multi_solute_trxn_Fe-EDTA_a.trc
    multi_solute_trac_Fe-EDTA[s].trc	       multi_solute_trxn_Fe-EDTA_s.trc
    multi_solute_trac_FreeIon_Cobalt[aq].trc   multi_solute_trxn_FreeIon_Cobalt.trc
    multi_solute_trac_FreeIon_EDTA[aq].trc     multi_solute_trxn_FreeIon_EDTA.trc
    multi_solute_trac_FreeIon_Iron[aq].trc     multi_solute_trxn_FreeIon_Iron.trc
    multi_solute_trac_Iron[aq].trc             multi_solute_trxn_Iron.trc        
    

    Test Directory: FEHM/fehmpytests/multi_solute

    Example File multi_solute.in

    
    COMPARE FEHMN and PDREACT: Linear Sorption w/ Surface Exchange  
    cond
    1 202 1 2.7 2.7 2.7
    
    ctrl
    50 1e-6 8
    1 202 1 2
    
    1 0 0.5
    25 2. 1.e-10 1.e-1
    1 0
    flow
    1 202 101 -0.05556 -25 0
    101 202 101 1. -25 -1
    
    init
    1. 25 25 0 1000 25 0 0
    node
    1
    202
    hist 
    concentration
    end
    perm
    1 202 1 5.0e-13 5.0e-30 5.0e-30
    
    rock
    1 202 1 1500 1000 0.4
    
    sol
    1 -1
    time
    1.e-6 7.25 1000 10 92 11
    
    # solute 1: Total Cobalt Concentration
    # solute 2: Total Iron Concentration
    # solute 3: Total EDTA Concentration
    # solute 4: CoEDTA adsorbed concentration
    # solute 5: Co adsorbed concentration
    # solute 6: FeEDTA adsorbed concentration
    trac
    1.d-80 1.0 1.e-6 0.5
    1. 2000 1.0 2000
    5 5.0 1.e-6 2.8333e-3
    6
    1
    1 0. 0. 1. 1.e-9 .05 1.e-34 1.e-34
    
    1 202 1 1
    
    
    1 202 101 3.1623e-5 1.0 4.16667
    
    1
    1 0. 0. 1. 1.e-9 .05 1.e-34 1.e-34
    
    1 202 1 1
    
    
    1 202 101 1.e-13 1.0 4.16667
    
    1
    1 0. 0. 1. 1.e-9 .05 1.e-34 1.e-34
    
    1 202 1 1
    
    
    1 202 101 3.1623e-5 1.0 4.16667
    
    0
    
    
    0
    
    
    0
    
    
    rxn
    ** NCPLX, NUMRXN
    2,4
    ** Coupling of the aqueous components (dRi/dUj)
    2
    1 0 1
    0 1 0
    ** IDCPNT(IC),CPNTNAM(IC),IFXCONC(IC),CPNTPRT(IC) (comp,name,cond.; NCPNT rows)
    1    Cobalt[aq]     0     0     1.e-9
    2      Iron[aq]     0     0     1.e-9
    3      EDTA[aq]     0     0     1.e-9
    ** IDCPLX(IX), CPLXNAM(IX),CPLXPRT(IX) (ID # and name of complex, NCPLX rows)
    101   Co-EDTA[aq]   0
    102   Fe-EDTA[aq]   0
    ** IDIMM(IM), IMMNAM(IM),IMMPRT(IM)(ID # and name of immoblie spec, NIMM rows)
    1   Co-EDTA[s]     0
    2   Fe-EDTA[s]     0
    3    Cobalt[s]     0
    ** IDVAP(IV), VAPNAM(IM), VAPPRT(IV) (ID # and name of vapor spec, NVAP rows)
    ** Skip nodes
    0
    ** RSDMAX
    1.0e-10
    **** Chemical reaction information for equilibrium reactions ******
    ** LOGKEQ (=0 if stability constants are given as K, =1 if given as log(K))
             0
    ** CKEQ(IX) ,HEQ(IX) (Stability constants and Enthaplys, NCPLX rows)
    1.0e+18    0
    6.31e+27   0
    ** STOIC(IX,IC) (Stoichiometric coeff: NCPLX rows, NCPNT columns)
    1.0       0.0       1.0
    0.0       1.0       1.0
    ** LINEAR KINETIC REACTION (type 1) **
            1
    ** Where does the reaction take place? **
     1 0 0
    
    ** Aqueous Component/Complex #, Solid Component #
          101      1
    ** Distribution coeffienct (kg water/ kg rock) **  
        0.533
    ** Mass transfer coefficient (1/hr) **
          1.0
    ** LINEAR KINETIC REACTION (type 1) **
            1
    ** Where does the reaction take place? **
     1 0 0
    
    ** Aqueous Component/Complex #, Solid Component #
            1      3
    ** Distribution coeffienct (kg rock/ kg water) **  
         5.07
    ** Mass transfer coefficient (1/hr) **
          1.0
    ** LINEAR KINETIC REACTION (type 1) **
            1
    ** Where does the reaction take place? **
     1 0 0
    
    ** Aqueous Component/Complex #, Solid Component #
          102      2
    ** Distribution coeffienct (kg rock/ kg water) **  
        0.427
    ** Mass transfer coefficient (1/hr) **
          1.0
    ** GENERAL EXCHANGE REACTION (type 3) **
            3
    ** Where does the reaction take place? **
     1 0 0
    
    ** # of solid, liquid and vapor species **
            3   0   0
    ** forward and reverse rate constants (1/hr) **
      1.26e-2 0
    ** Solid Species in reaction **
      1      2      3
    ** Stoichiometry **
      1.0   -1.0  -1.0
    coor   n/a     
       202
             1        0.00000        1.00000        0.00000
             2        0.10000        1.00000        0.00000
             3        0.20000        1.00000        0.00000
             4        0.30000        1.00000        0.00000
             5        0.40000        1.00000        0.00000
             6        0.50000        1.00000        0.00000
             7        0.60000        1.00000        0.00000
             8        0.70000        1.00000        0.00000
             9        0.80000        1.00000        0.00000
            10        0.90000        1.00000        0.00000
            11        1.00000        1.00000        0.00000
            12        1.10000        1.00000        0.00000
            13        1.20000        1.00000        0.00000
            14        1.30000        1.00000        0.00000
            15        1.40000        1.00000        0.00000
            16        1.50000        1.00000        0.00000
            17        1.60000        1.00000        0.00000
            18        1.70000        1.00000        0.00000
            19        1.80000        1.00000        0.00000
            20        1.90000        1.00000        0.00000
            21        2.00000        1.00000        0.00000
            22        2.10000        1.00000        0.00000
            23        2.20000        1.00000        0.00000
            24        2.30000        1.00000        0.00000
            25        2.40000        1.00000        0.00000
            26        2.50000        1.00000        0.00000
            27        2.60000        1.00000        0.00000
            28        2.70000        1.00000        0.00000
            29        2.80000        1.00000        0.00000
            30        2.90000        1.00000        0.00000
            31        3.00000        1.00000        0.00000
            32        3.10000        1.00000        0.00000
            33        3.20000        1.00000        0.00000
            34        3.30000        1.00000        0.00000
            35        3.40000        1.00000        0.00000
            36        3.50000        1.00000        0.00000
            37        3.60000        1.00000        0.00000
            38        3.70000        1.00000        0.00000
            39        3.80000        1.00000        0.00000
            40        3.90000        1.00000        0.00000
            41        4.00000        1.00000        0.00000
            42        4.10000        1.00000        0.00000
            43        4.20000        1.00000        0.00000
            44        4.30000        1.00000        0.00000
            45        4.40000        1.00000        0.00000
            46        4.50000        1.00000        0.00000
            47        4.60000        1.00000        0.00000
            48        4.70000        1.00000        0.00000
            49        4.80000        1.00000        0.00000
            50        4.90000        1.00000        0.00000
            51        5.00000        1.00000        0.00000
            52        5.10000        1.00000        0.00000
            53        5.20000        1.00000        0.00000
            54        5.30000        1.00000        0.00000
            55        5.40000        1.00000        0.00000
            56        5.50000        1.00000        0.00000
            57        5.60000        1.00000        0.00000
            58        5.70000        1.00000        0.00000
            59        5.80000        1.00000        0.00000
            60        5.90000        1.00000        0.00000
            61        6.00000        1.00000        0.00000
            62        6.10000        1.00000        0.00000
            63        6.20000        1.00000        0.00000
            64        6.30000        1.00000        0.00000
            65        6.40000        1.00000        0.00000
            66        6.50000        1.00000        0.00000
            67        6.60000        1.00000        0.00000
            68        6.70000        1.00000        0.00000
            69        6.80000        1.00000        0.00000
            70        6.90000        1.00000        0.00000
            71        7.00000        1.00000        0.00000
            72        7.10000        1.00000        0.00000
            73        7.20000        1.00000        0.00000
            74        7.30000        1.00000        0.00000
            75        7.40000        1.00000        0.00000
            76        7.50000        1.00000        0.00000
            77        7.60000        1.00000        0.00000
            78        7.70000        1.00000        0.00000
            79        7.80000        1.00000        0.00000
            80        7.90000        1.00000        0.00000
            81        8.00000        1.00000        0.00000
            82        8.10000        1.00000        0.00000
            83        8.20000        1.00000        0.00000
            84        8.30000        1.00000        0.00000
            85        8.40000        1.00000        0.00000
            86        8.50000        1.00000        0.00000
            87        8.60000        1.00000        0.00000
            88        8.70000        1.00000        0.00000
            89        8.80000        1.00000        0.00000
            90        8.90000        1.00000        0.00000
            91        9.00000        1.00000        0.00000
            92        9.10000        1.00000        0.00000
            93        9.20000        1.00000        0.00000
            94        9.30000        1.00000        0.00000
            95        9.40000        1.00000        0.00000
            96        9.50000        1.00000        0.00000
            97        9.60000        1.00000        0.00000
            98        9.70000        1.00000        0.00000
            99        9.80000        1.00000        0.00000
           100        9.90000        1.00000        0.00000
           101       10.00000        1.00000        0.00000
           102        0.00000        0.00000        0.00000
           103        0.10000        0.00000        0.00000
           104        0.20000        0.00000        0.00000
           105        0.30000        0.00000        0.00000
           106        0.40000        0.00000        0.00000
           107        0.50000        0.00000        0.00000
           108        0.60000        0.00000        0.00000
           109        0.70000        0.00000        0.00000
           110        0.80000        0.00000        0.00000
           111        0.90000        0.00000        0.00000
           112        1.00000        0.00000        0.00000
           113        1.10000        0.00000        0.00000
           114        1.20000        0.00000        0.00000
           115        1.30000        0.00000        0.00000
           116        1.40000        0.00000        0.00000
           117        1.50000        0.00000        0.00000
           118        1.60000        0.00000        0.00000
           119        1.70000        0.00000        0.00000
           120        1.80000        0.00000        0.00000
           121        1.90000        0.00000        0.00000
           122        2.00000        0.00000        0.00000
           123        2.10000        0.00000        0.00000
           124        2.20000        0.00000        0.00000
           125        2.30000        0.00000        0.00000
           126        2.40000        0.00000        0.00000
           127        2.50000        0.00000        0.00000
           128        2.60000        0.00000        0.00000
           129        2.70000        0.00000        0.00000
           130        2.80000        0.00000        0.00000
           131        2.90000        0.00000        0.00000
           132        3.00000        0.00000        0.00000
           133        3.10000        0.00000        0.00000
           134        3.20000        0.00000        0.00000
           135        3.30000        0.00000        0.00000
           136        3.40000        0.00000        0.00000
           137        3.50000        0.00000        0.00000
           138        3.60000        0.00000        0.00000
           139        3.70000        0.00000        0.00000
           140        3.80000        0.00000        0.00000
           141        3.90000        0.00000        0.00000
           142        4.00000        0.00000        0.00000
           143        4.10000        0.00000        0.00000
           144        4.20000        0.00000        0.00000
           145        4.30000        0.00000        0.00000
           146        4.40000        0.00000        0.00000
           147        4.50000        0.00000        0.00000
           148        4.60000        0.00000        0.00000
           149        4.70000        0.00000        0.00000
           150        4.80000        0.00000        0.00000
           151        4.90000        0.00000        0.00000
           152        5.00000        0.00000        0.00000
           153        5.10000        0.00000        0.00000
           154        5.20000        0.00000        0.00000
           155        5.30000        0.00000        0.00000
           156        5.40000        0.00000        0.00000
           157        5.50000        0.00000        0.00000
           158        5.60000        0.00000        0.00000
           159        5.70000        0.00000        0.00000
           160        5.80000        0.00000        0.00000
           161        5.90000        0.00000        0.00000
           162        6.00000        0.00000        0.00000
           163        6.10000        0.00000        0.00000
           164        6.20000        0.00000        0.00000
           165        6.30000        0.00000        0.00000
           166        6.40000        0.00000        0.00000
           167        6.50000        0.00000        0.00000
           168        6.60000        0.00000        0.00000
           169        6.70000        0.00000        0.00000
           170        6.80000        0.00000        0.00000
           171        6.90000        0.00000        0.00000
           172        7.00000        0.00000        0.00000
           173        7.10000        0.00000        0.00000
           174        7.20000        0.00000        0.00000
           175        7.30000        0.00000        0.00000
           176        7.40000        0.00000        0.00000
           177        7.50000        0.00000        0.00000
           178        7.60000        0.00000        0.00000
           179        7.70000        0.00000        0.00000
           180        7.80000        0.00000        0.00000
           181        7.90000        0.00000        0.00000
           182        8.00000        0.00000        0.00000
           183        8.10000        0.00000        0.00000
           184        8.20000        0.00000        0.00000
           185        8.30000        0.00000        0.00000
           186        8.40000        0.00000        0.00000
           187        8.50000        0.00000        0.00000
           188        8.60000        0.00000        0.00000
           189        8.70000        0.00000        0.00000
           190        8.80000        0.00000        0.00000
           191        8.90000        0.00000        0.00000
           192        9.00000        0.00000        0.00000
           193        9.10000        0.00000        0.00000
           194        9.20000        0.00000        0.00000
           195        9.30000        0.00000        0.00000
           196        9.40000        0.00000        0.00000
           197        9.50000        0.00000        0.00000
           198        9.60000        0.00000        0.00000
           199        9.70000        0.00000        0.00000
           200        9.80000        0.00000        0.00000
           201        9.90000        0.00000        0.00000
           202       10.00000        0.00000        0.00000
             0        0.00000        0.00000        0.00000
    elem
       4  100  0
           1     102     103       2       1
           2     103     104       3       2
           3     104     105       4       3
           4     105     106       5       4
           5     106     107       6       5
           6     107     108       7       6
           7     108     109       8       7
           8     109     110       9       8
           9     110     111      10       9
          10     111     112      11      10
          11     112     113      12      11
          12     113     114      13      12
          13     114     115      14      13
          14     115     116      15      14
          15     116     117      16      15
          16     117     118      17      16
          17     118     119      18      17
          18     119     120      19      18
          19     120     121      20      19
          20     121     122      21      20
          21     122     123      22      21
          22     123     124      23      22
          23     124     125      24      23
          24     125     126      25      24
          25     126     127      26      25
          26     127     128      27      26
          27     128     129      28      27
          28     129     130      29      28
          29     130     131      30      29
          30     131     132      31      30
          31     132     133      32      31
          32     133     134      33      32
          33     134     135      34      33
          34     135     136      35      34
          35     136     137      36      35
          36     137     138      37      36
          37     138     139      38      37
          38     139     140      39      38
          39     140     141      40      39
          40     141     142      41      40
          41     142     143      42      41
          42     143     144      43      42
          43     144     145      44      43
          44     145     146      45      44
          45     146     147      46      45
          46     147     148      47      46
          47     148     149      48      47
          48     149     150      49      48
          49     150     151      50      49
          50     151     152      51      50
          51     152     153      52      51
          52     153     154      53      52
          53     154     155      54      53
          54     155     156      55      54
          55     156     157      56      55
          56     157     158      57      56
          57     158     159      58      57
          58     159     160      59      58
          59     160     161      60      59
          60     161     162      61      60
          61     162     163      62      61
          62     163     164      63      62
          63     164     165      64      63
          64     165     166      65      64
          65     166     167      66      65
          66     167     168      67      66
          67     168     169      68      67
          68     169     170      69      68
          69     170     171      70      69
          70     171     172      71      70
          71     172     173      72      71
          72     173     174      73      72
          73     174     175      74      73
          74     175     176      75      74
          75     176     177      76      75
          76     177     178      77      76
          77     178     179      78      77
          78     179     180      79      78
          79     180     181      80      79
          80     181     182      81      80
          81     182     183      82      81
          82     183     184      83      82
          83     184     185      84      83
          84     185     186      85      84
          85     186     187      86      85
          86     187     188      87      86
          87     188     189      88      87
          88     189     190      89      88
          89     190     191      90      89
          90     191     192      91      90
          91     192     193      92      91
          92     193     194      93      92
          93     194     195      94      93
          94     195     196      95      94
          95     196     197      96      95
          96     197     198      97      96
          97     198     199      98      97
          98     199     200      99      98
          99     200     201     100      99
         100     201     202     101     100
           0       0       0       0       0
    stop