salt_perm_poro
Test the Salt Permeability and Porosity Macro
The porosity-permeability function for compacted salt from Cinar et at. (2006) is tested using a six node problem with porosities from 0.01 to 0.2. The excel spreadsheet in information/salt-perm-poro.xlsx contains calculations of the perm-poro function.
Cinar, Y, G Pusch and V Reitenbach (2006) Petrophysical and capillary properties of compacted salt. Transport in Porous Media. 64, p. 199-228, doi: 10.1007/s11242-005-2848-1
Test Directory: FEHM/fehmpytests/salt_perm_poro
Example File run.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.05
3 3 1 2165. 931. 0.1
4 4 1 2165. 931. 0.15
5 5 1 2165. 931. 0.155
6 6 1 2165. 931. 0.2
#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
surf 100 10000.
porosity
perm
end
#mat
#conc
#pres
#vap
#density
#liquid
#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 A[aq] 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 A[s] 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