saltvcon

    Test the Salt Variable Conductivity Macro

    Tests the calculations of thermal conductivity of crushed and intact salt.

    Intact salt:

    kxi = k_{t-300}(300/T)^1.14 Munson et al. (1990) Overtest for Simulate Defense High-Level Waste (Laboratories, SAND89-2671)

    Thermal conductivity of crushed salt from Asse mine:

    kx_asse = -270*phi^4+370*phi^3-136*phi^2+1.5*phi+5 Bechtold et al. (2004) Backfilling and sealing of underground respositories for radioactive waste in salt (BAMBUS II project), EUR 20621, ISBN 92-894-7767-9

    kx = (k_{t-300}/kx_asse)*(300/T)^1.14 if kx is less then 1.e-6, set to 1.e-6. If porosity is greater than 0.4, it is truncated as 0.4 since the Kx relationship is only valid within this range.

    The excel spreadsheet /information/saltvcon.xlsx contains the associated calculations.

    Test Directory: FEHM/fehmpytests/saltvcon

    Example File intact.dat

    
    title: 1-d heat pipe calculation               11/28/12
    #   All nodes are crushed salt
    #************************************************************************75
    zone
    	1
    nnum
    	6 1 2 3 4 5 6
    
    salt
    saltvapr
    1
    saltnum
    permavg 1.0
    poravg  1.0
    pormin 1.d-5
    
    saltvcon
      3   26.85  5.4 1.14
      4   26.85  1.08 -270. 370. -136. 1.5 5  1.14
    
      1 0 0 1
      1 0 0 2
    
    saltppor 
    7
      1 0 0   4.866e-9 4.637 1.e-3  0.20
     
    saltadif
     333
    saltend    
    node
    6
    1 2 3 4 5 6
    perm 
       1 0 0 -14 -14 -14
       1 0 0 -13 -13 -13
    
    rlp
            3  0.05 1.0  4   1.56   -10.  0.06
            1  0.0  0.   1.0    1.0   0.15  1.0
    
            1      0       0       2
    
    rock  
     1 1 1   2165.  931.  0.01
     2 2 1   2165.  931.  0.1
     3 3 1   2165.  931.  0.3
     4 4 1   2165.  931.  0.5
     5 5 1   2165.  931.  0.7
     6 6 1   2165.  931.  0.9
    
    #vcon
    #  3   26.85  5.4   1.14
    #  4   26.85  1.08 -270. 370. -136. 1.5 5  1.14
    #
    #  1 0 0 2
    #
    #ppor
    #7
    #  1 0 0   4.866e-9 4.637 1.e-3  0.20
    #
    flxo 
    5
    1 2 
    2 3
    3 4
    4 5
    5 6
    # - - - - - - - - - - - - - - - - 
    time
    1.e-4  200.    1    01  1995  5  0.0 
    0.5 -2.0 1 1
    1.0 -2.0 1 1
    5.0 -2.0 1 1
    10.0 -2.0 1 1
    
    ctrl
      -15   1.e-04   24  100 gmre
    1   0 0 2
    0 
    1.0   0.0  1.
    15   1.5  1.e-09 1.  
    0 +1 
    iter
    1.e-5 1.e-5 1.e-3 -1.e-3 1.2
    00 0 0 10 1000.
    sol
    1    -1
    #- - - - - - - - - - - - - - -
    #vapl    #gaz debug comment out
    #adif
    #333
    #- - - - - - - - - - - - -
    pres       #initial pres sat
    1 0 0 0.1   0.10  2
    
    ngas  reset P
    3
    1  1  1 -20.
    2  2  1 -40.
    3  3  1 -60.
    4  4  1 -80.
    5  5  1 -100.
    6  6  1 -120.
    
    
    
    hflx
     1 1 1 20.   1.e6
     6 6 1 120.   1.e6
    
    #- - - - - - - - - - - - -
    #cont
    #avsx  100  10000.
    #temp
    #sat
    #porosity
    #perm
    #mat
    #conc
    #pres
    #vap  
    #density
    #liquid
    #end
    #cden now in moles 226/7.5 = *****************************************
    cden
    1
    30.1
    #***********************************************************
    trac
     0  1  1.e-7  1.0
     0. 3652.5  1.e6  1.e6
     50  1.6 1.e-1 1.  1
     2
     1
    0 0 0 1  1.e-9 .33333 .33333 .33333
    
      1 0 0 1
    
      1 0 0  6.16 
    
    
    0
      1 0 0 17.2414
    
    
    rxn
    ** NCPLX,NUMRXN
             0, 1
    ** Coupling of the aqueous components (dRi/dUj)
    1
    1
    ** IDCPNT(IC),CPNTNAM(IC),IFXCONC(IC),CPNTPRT(IC) (comp,name,cond.; NCPNT rows)
        1     1    0      0    1.e-9
    ** IDCPLX(IX), CPLXNAM(IX),CPLXPRT(IX) (ID # and name of complex, NCPLX rows)
    ** IDIMM(IM), IMMNAM(IM),IMMPRT(IM)(ID # and name of immobile spec, NIMM rows)
        1  2    0
    ** IDVAP(IV), VAPNAM(IM),VAPPRT(IV) (ID # and name of vapor species, NVAP rows)
    ** skip nodes? **
    0   -1
    ** RSDMAX
       1.0e-9
    ******  Chemical reaction information  ********
    ** LOGKEQ (=0 if stability constants are given as K, =1 if given as log(K))
    ** CKEQ(IX), HEQ(IX) (Stability constants and Enthalpys, NCPLX rows)
    ** STOIC(IX,IC) (Stoichiometric coeff: NCPLX rows, NCPNT columns)
    ** Precipiation/Dissolution REACTION (type 7) **
            8
    ** Where does the reaction take place ? ***
    1 0 0
    
    ** immobile species participating in reaction **
            1
    ** the number of total aqueous species in reaction **
            1
    ** total aqueous species in reaction **
            1
    ** stoichiometry of the immobilie species **
            1
    ** stoichiometry of the aqueous species **
            1
    ** solubility product **
    lookup 8
         10 40 100 150 200 250 300 350
        6.12 6.23  6.65 7.21  8.00  9.06  10.45  12.27
    porosity change
    ** molecular weight of mineral (kg/mol), density of mineral (kg/m^3) SALT Wikipedia**
    0.0558 , 2165.
    ** rate constant (moles/(m^2*sec)) **
            0.01
    ** surface area of the mineral (m^2) **
            1
    stop
    
    ** rate constant (moles/(m^2*sec)) **
            0.1