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