*************************************************************************
TITLE: T.01 - TESTING SOLUBILITY LIMITED SOURCE OPTION
************************************************************************
/
GRID 12 BY 12
COORDINATE X RANGE = 10
COORDINATE Y RANGE = 10
HYDRAULIC properties: ss=0.001, k = 3*1.
TRANSport kd=0., dm=0.
SET P LINEar function 0. -1. X
/
INITial C = 0.
SELEct (6,6)
SOURce for C: SOLUbility: TABLe of 5 sets (TIME,VAL) SELEcted
(0, 0) (20,0) (20.1, 1) (100,1) (150,1)
total = 100 kg, time = 0.
DECAY HALF LIFE FOR C = 10
DIAgnostics U, C at: (6,6) every 20 steps
DEBUG GEOMETRY OFF
SOLVE 100 secs in steps of 1
output c
SOLVE 20 secs in steps of 1
output c
SOLVE 20 secs in steps of 1
output c
END
/
-Back To Top
************************************************************************
TITLE: T.02 - TESTING
SOLUBILITY LIMITED SOURCE OPTION
************************************************************************
/
GRID 12 BY 12
COORDINATE X RANGE = 10
COORDINATE Y RANGE = 10
HYDRAULIC properties: ss=0.001, k = 3*1.
TRANSport kd=0., dm=0.
SET P LINEar function 0. -1. X
/
INITial C = 0.
SELEct (6,6)
SOURce for C: SOLUbility = 1, total = 100 kg, time = 21. SELEcted
DECAY HALF LIFE FOR C = 20
DIAgnostics TIME U, C at: (6,6) every 20 steps
DEBUG GEOMETRY OFF
SOLVE 100 secs in steps of 0.1
output c
/SOLVE 20 secs in steps of 1
/output c
/SOLVE 20 secs in steps of 1
/output c
END
/
-Back To Top
************************************************************************
TITLE: T.03 - TESTING
SOLUBILITY LIMITED SOURCE OPTION
************************************************************************
/
GRID 12 BY 12
COORDINATE X RANGE = 10
COORDINATE Y RANGE = 10
HYDRAULIC properties: ss=0.001, k = 3*1.
TRANSport kd=0., dm=1.
SET P LINEar function 0. -1.E-30 X
/
SELECT (2,2) TO (11,11)
SET C = 1 SELECTED
/SELEct (6,6)
/SOURce for C: SOLUbility = 1, total = 100 kg, time = 21. SELEcted
CONVERGENCE C =1.E-6, MAX ITERATIONS: 20
DECAY HALF LIFE FOR C = 20
DIAgnostics TIME U, C at: (6,6) every 20 steps
DEBUG GEOMETRY OFF
SOLVE 100 secs in steps of 1.
output c
END
/
-Back To Top
************************************************************************
TITLE: T.04 - TESTING
DECAY RELATED OPTIONS
************************************************************************
/
Grid 100 x
COORDINATE X RANGE 1
MATERIAL DENSITY 2000
MATERIAL POROSITY 1
TRANSPORT C 0 0 0 0
BOUNDARY C X- FLUX = 0.
BOUNDARY C X+ FLUX = 0.
/
SELECT (2,1) TO (99,3) ! Option A
SOURCE C SOLU 1 1e+10 SELE ! Option A
DECAY C 0.693147 ! Option A
/
/ INITial C 1 ! Option B
/ DECAY half life C 1 ! Option B
/
DIAG (50,2) every 20 steps
CONVERGENCE C 1.E-10 100
/
DEBUG GEOMETRY OFF
SOLVE C 17 0.01
/
END
/
-Back To Top
************************************************************************
TITLE: T.05 -
Precipitation Test - Precipitation from liquid to solid
************************************************************************
/ Test 1: Initial Clqd=5.; Csld=0. No incoming Flux; K=1000 (fast
chem)
////////////////////////////////////////////////////////////////////////
GRID IS 3 BY 42
/
coor y RANGE 1
/
FOR MATErial 1
MATErial particle DENSity: 1
MATErial POROsity: NE=0.50
HYDRaulic properties as=1.0E-07, Kx=1 m/yr, Ky=1
/
DENSITY 1 !Kg/m3
/
RENAME C = "Clqd - Liquid State Concentration of ANDRA Species"
RENAME C2 = "Csld - Solid State Concentration of ANDRA Species"
ALLOCATE "PREAC - Liquid-Solid Precipitation Reaction Rate"
/
TRANSPORT Clqd Kd=0.0, Dm=0.1 m-sq/yr
TRANSPORT Csld in SOLID phase
/
SET Clqd = 5.0
SET Csld = 0.0
/
BOUNdary FOR Clqd: at Y- VALUE = 0
BOUNdary FOR Clqd: at Y+ flux = 0
BOUNdary FOR Csld: at Y+ flux = 0
/
CONV Clqd 1.E-5 20
PRECipitation reaction PREAC: Clqd & Csld freq=1000. Solublity =
.02
/
DIAGnostic time Clqd Csld BClqd (2,11) every 100 steps
/
DEBUG GEOMETRY OFF
SOLVE 1 YEARS IN STEPS OF 0.01
SELECT (2,1) (2,999)
OUTPut Clqd Csld in NARRow NOW SELEcted
/
save off
/
SOLVE 1 YEARS IN STEPS OF 0.01 1.005 0.02
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 3 YEARS IN STEPS OF 0.02 1.005 0.05
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 5 YEARS IN STEPS OF 0.05 1.005 0.10
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 10 YEARS IN STEPS OF 0.10 1.005 0.20
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 30 YEARS IN STEPS OF 0.20 1.005 0.50
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 50 YEARS IN STEPS OF 0.50 1.005 1.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 100 YEARS IN STEPS OF 1.0 1.005 2.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
END
/
-Back To Top
************************************************************************
TITLE: T.06 -
Precipitation Test - Uptake from solid to liquid
************************************************************************
/ Test 1: Initial Clqd=5.; Csld=0. No incoming Flux; K=1000 (fast
chem)
/ Test 2: Initial Clqd=0.; Csld=5. No incoming Flux; K=1000 (fast
chem)
////////////////////////////////////////////////////////////////////////
GRID IS 3 BY 42
/
coor y range 1.
/
FOR MATErial 1
MATErial particle DENSity: 1
MATErial POROsity: NE=0.50
HYDRaulic properties as=1.0E-07, Kx=1 m/yr, Ky=1
/
DENSITY 1 !Kg/m3
/
RENAME C = "Clqd - Liquid State Concentration of ANDRA Species"
RENAME C2 = "Csld - Solid State Concentration of ANDRA Species"
ALLOCATE "PREAC - Liquid-Solid Precipitation Reaction Rate"
/
TRANSPORT Clqd Kd=0.0, Dm=0.1 m-sq/yr
TRANSPORT Csld in SOLID phase
/
SET Clqd = 0.0
SET Csld = 5.0
/
BOUNdary FOR Clqd: at Y- VALUE = 0
BOUNdary FOR Clqd: at Y+ flux = 0
BOUNdary FOR Csld: at Y+ flux = 0
/
CONV Clqd 1.E-5 20
PRECipitation reaction PREAC: Clqd & Csld freq=1000. Solublity =
.02
/
DIAGnostic time Clqd Csld BClqd (2,11) every 100 steps
/
DEBUG GEOMETRY OFF
SOLVE 1 YEARS IN STEPS OF 0.01
SELECT (2,1) (2,999)
OUTPut Clqd Csld in NARRow NOW SELEcted
/
save off
/
SOLVE 1 YEARS IN STEPS OF 0.01 1.005 0.02
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 3 YEARS IN STEPS OF 0.02 1.005 0.05
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 5 YEARS IN STEPS OF 0.05 1.005 0.10
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 10 YEARS IN STEPS OF 0.10 1.005 0.20
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 30 YEARS IN STEPS OF 0.20 1.005 0.50
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 50 YEARS IN STEPS OF 0.50 1.005 1.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 100 YEARS IN STEPS OF 1.0 1.005 2.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
END
/
-Back To Top
************************************************************************
TITLE: T.07 -
Precipitation Test - Uptake from solid to liquid
************************************************************************
/ Test 1: Initial Clqd=5.; Csld=0. No incoming Flux; K=1000 (fast
chem)
/ Test 2: Initial Clqd=0.; Csld=5. No incoming Flux; K=1000 (fast
chem)
/ Test 3: Same as 2 but K=1 (slow chemistry)
////////////////////////////////////////////////////////////////////////
GRID IS 3 BY 42
/
coor y range 1.
/
FOR MATErial 1
MATErial particle DENSity: 1
MATErial POROsity: NE=0.50
HYDRaulic properties as=1.0E-07, Kx=1 m/yr, Ky=1
/
DENSITY 1 !Kg/m3
/
RENAME C = "Clqd - Liquid State Concentration of ANDRA Species"
RENAME C2 = "Csld - Solid State Concentration of ANDRA Species"
ALLOCATE "PREAC - Liquid-Solid Precipitation Reaction Rate"
/
TRANSPORT Clqd Kd=0.0, Dm=0.1 m-sq/yr
TRANSPORT Csld in SOLID phase
/
SET Clqd = 0.0
SET Csld = 5.0
/
BOUNdary FOR Clqd: at Y- VALUE = 0
BOUNdary FOR Clqd: at Y+ flux = 0
BOUNdary FOR Csld: at Y+ flux = 0
/
CONV Clqd 1.E-5 20
PRECipitation reaction PREAC: Clqd & Csld freq=1. Solublity = .02
/
DIAGnostic OFF time Clqd Csld BClqd (2,11) every 100 steps
/
DEBUG GEOMETRY OFF
SOLVE 1 YEARS IN STEPS OF 0.01
SELECT (2,1) (2,999)
OUTPut Clqd Csld in NARRow NOW SELEcted
/
save off
/
DEBUG GEOMETRY OFF
SOLVE 1 YEARS IN STEPS OF 0.01 1.005 0.02
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 3 YEARS IN STEPS OF 0.02 1.005 0.05
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 5 YEARS IN STEPS OF 0.05 1.005 0.10
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 10 YEARS IN STEPS OF 0.10 1.005 0.20
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 30 YEARS IN STEPS OF 0.20 1.005 0.50
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 50 YEARS IN STEPS OF 0.50 1.005 1.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 100 YEARS IN STEPS OF 1.0 1.005 2.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
END
/
-Back To Top
************************************************************************
TITLE: T.08 -
Precipitation Test - Precipitation from liquid to solid
************************************************************************
/ Test 1: Initial Clqd=5.; Csld=0. No incoming Flux; K=1000 (fast
chem)
/ Test 2: Initial Clqd=0.; Csld=5. No incoming Flux; K=1000 (fast
chem)
/ Test 3: Same as 2 but K=1 (slow chemistry)
/ Test 4: Initial Clqd = 5.; Csld = 0.; kinetic_K=1. (slow
reaction)
/ Liquid maintained at 5 at X- boundary
/ Thus significant diffusive inflow at X-: qd ~ 0.1*(5-0.2)
////////////////////////////////////////////////////////////////////////
GRID IS 3 BY 42
/
coor y range 1.
/
FOR MATErial 1
MATErial particle DENSity: 1
MATErial POROsity: NE=0.50
HYDRaulic properties as=1.0E-07, Kx=1 m/yr, Ky=1
/
DENSITY 1 !Kg/m3
/
RENAME C = "Clqd - Liquid State Concentration of ANDRA Species"
RENAME C2 = "Csld - Solid State Concentration of ANDRA Species"
ALLOCATE "PREAC - Liquid-Solid Precipitation Reaction Rate"
/
TRANSPORT Clqd Kd=0.0, Dm=0.1 m-sq/yr
TRANSPORT Csld in SOLID phase
/
SET Clqd = 5.0
SET Csld = 0.0
/
BOUNdary FOR Clqd: at Y- VALUE = 0
BOUNdary FOR Clqd: at Y+ VALUE = 5
BOUNdary FOR Csld: at Y+ flux = 0
/
CONV Clqd 1.E-5 20
PRECipitation reaction PREAC: Clqd & Csld freq=1. Solublity = .02
/
DIAGnostic time Clqd Csld BClqd (2,11) every 100 steps
/
DEBUG GEOMETRY OFF
SOLVE 1 YEARS IN STEPS OF 0.01
SELECT (2,1) (2,999)
OUTPut Clqd Csld in NARRow NOW SELEcted
/
save off
/
SOLVE 1 YEARS IN STEPS OF 0.01 1.005 0.02
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 3 YEARS IN STEPS OF 0.02 1.005 0.05
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 5 YEARS IN STEPS OF 0.05 1.005 0.10
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 10 YEARS IN STEPS OF 0.10 1.005 0.20
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 30 YEARS IN STEPS OF 0.20 1.005 0.50
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 50 YEARS IN STEPS OF 0.50 1.005 1.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 100 YEARS IN STEPS OF 1.0 1.005 2.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
END
/
-Back To Top
************************************************************************
TITLE: T.09 -
Precipitation Test - Precipitation from liquid to solid
************************************************************************
/ Test 1: Initial Clqd=5.; Csld=0. No incoming Flux; K=1000 (fast
chem)
/ Test 2: Initial Clqd=0.; Csld=5. No incoming Flux; K=1000 (fast
chem)
/ Test 3: Same as 2 but K=1 (slow chemistry)
/ Test 4: Initial Clqd = 5.; Csld = 0.; kinetic_K=1. (slow
reaction)
/ Liquid maintained at 5 at X- boundary
/ Thus significant diffusive inflow at X-: qd ~ 0.1*(5-0.2)
/ Test 5: Similar to 5 but carried to steady state with:
/ Diffusion in solid phase
/ Boundary C_lqd = 0.2; C_sld=4.8
////////////////////////////////////////////////////////////////////////
GRID IS 3 BY 42
/
coor y range 1.
/
FOR MATErial 1
MATErial particle DENSity: 1
MATErial POROsity: NE=0.50
HYDRaulic properties as=1.0E-07, Kx=1 m/yr, Ky=1
/
DENSITY 1 !Kg/m3
/
RENAME C = "Clqd - Liquid State Concentration of ANDRA Species"
RENAME C2 = "Csld - Solid State Concentration of ANDRA Species"
ALLOCATE "PREAC - Liquid-Solid Precipitation Reaction Rate"
/
TRANSPORT Clqd Kd=0.0, Dm=0.1 m^2/yr
TRANSPORT Csld in SOLID phase
DIFFusivity for Csld = 0.0001 m^2/yr
/
SET Clqd = 5.0
SET Csld = 0.0
/
BOUNdary FOR Clqd: at Y- VALUE = 0
BOUNdary FOR Csld: at Y- value = 0
BOUNdary FOR Clqd: at Y+ VALUE = 0.2
BOUNdary FOR Csld: at Y+ value = 4.8
/
CONV Clqd 1.E-5 20
CONV Csld 1.E-3 50
PRECipitation reaction PREAC: Clqd & Csld freq=1. Solublity = .02
/
DIAGnostic time dtime Clqd Csld BClqd (2,11) every 100 iterations
/
DEBUG GEOMETRY OFF
SOLVE 1 YEARS IN STEPS OF 0.01
SELECT (2,1) (2,999)
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
save off
/
SOLVE 1 YEARS IN STEPS OF 0.01 1.005 0.02
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 3 YEARS IN STEPS OF 0.02 1.005 0.05
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 5 YEARS IN STEPS OF 0.05 1.005 0.10
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 10 YEARS IN STEPS OF 0.10 1.005 0.20
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 30 YEARS IN STEPS OF 0.20 1.005 0.50
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 50 YEARS IN STEPS OF 0.50 1.005 1.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 100 YEARS IN STEPS OF 1.0 1.005 1.00
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 800 YEARS IN STEPS OF 1.0 1.005 1.00
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 1000 YEARS IN STEPS OF 1.0 1.005 1.00
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 1000 YEARS IN STEPS OF 1.0 1.005 1.00
OUTPut Clqd Csld in NARRow NOW SELEcted
/
END
/
-Back To Top
************************************************************************
TITLE: T.10 -
Precipitation Test - Precipitation from liquid to solid
************************************************************************
/ Test 1: Initial Clqd=5.; Csld=0. No incoming Flux; K=1000 (fast
chem)
/ Test 2: Initial Clqd=0.; Csld=5. No incoming Flux; K=1000 (fast
chem)
/ Test 3: Same as 2 but K=1 (slow chemistry)
/ Test 4: Initial Clqd = 5.; Csld = 0.; kinetic_K=1. (slow
reaction)
/ Liquid maintained at 5 at X- boundary
/ Thus significant diffusive inflow at X-: qd ~ 0.1*(5-0.2)
/ Test 5: Similar to 5 but carried to steady state with:
/ Diffusion in solid phase
/ Boundary C_lqd = 0.2; C_sld=4.8
/ Test 6: Same to 5: set boundary to 0 at 3000 yrs and deplete
value
////////////////////////////////////////////////////////////////////////
GRID IS 3 BY 42
/
coor y range 1.
/
FOR MATErial 1
MATErial particle DENSity: 1
MATErial POROsity: NE=0.50
HYDRaulic properties as=1.0E-07, Kx=1 m/yr, Ky=1
/
DENSITY 1 !Kg/m3
/
RENAME C = "Clqd - Liquid State Concentration of ANDRA Species"
RENAME C2 = "Csld - Solid State Concentration of ANDRA Species"
ALLOCATE "PREAC - Liquid-Solid Precipitation Reaction Rate"
/
TRANSPORT Clqd Kd=0.0, Dm=0.1 m^2/yr
TRANSPORT Csld in SOLID phase
DIFFusivity for Csld = 0.0001 m^2/yr
/
SET Clqd = 5.0
SET Csld = 0.0
/
BOUNdary FOR Clqd: at Y- VALUE = 0
BOUNdary FOR Csld: at Y- value = 0
BOUNdary FOR Clqd: at Y+ table 4 SETS TIME VALUE
0. 0.2
3000. 0.2
3010. 0.0
99999. 0.0
BOUNdary FOR Csld: at Y+ table 4 SETS TIME VALUE
0. 4.8
3000. 4.8
3010. 0.0
99999. 0.0
/
CONV Clqd 1.E-5 20
CONV Csld 1.E-3 50
PRECipitation reaction PREAC: Clqd & Csld freq=1. Solublity = .02
/
DIAGnostic time dtime Clqd Csld BClqd (2,11) every 100 iterations
/
DEBUG GEOMETRY OFF
SOLVE 1 YEARS IN STEPS OF 0.01
SELECT (2,1) (2,999)
OUTPut Clqd Csld in NARRow NOW SELEcted
/
save off
/
SOLVE 1 YEARS IN STEPS OF 0.01 1.005 0.02
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 3 YEARS IN STEPS OF 0.02 1.005 0.05
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 5 YEARS IN STEPS OF 0.05 1.005 0.10
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 10 YEARS IN STEPS OF 0.10 1.005 0.20
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 30 YEARS IN STEPS OF 0.20 1.005 0.50
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 50 YEARS IN STEPS OF 0.50 1.005 1.0
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 100 YEARS IN STEPS OF 1.0 1.005 1.00
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 800 YEARS IN STEPS OF 1.0 1.005 1.00
OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 1000 YEARS IN STEPS OF 1.0 1.005 1.00
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 1000 YEARS IN STEPS OF 1.0 1.005 1.00
/OUTPut Clqd Csld in NARRow NOW SELEcted
/
SOLVE 1000 YEARS IN STEPS OF 1.0 1.005 1.00
OUTPut Clqd Csld in NARRow NOW SELEcted
/
END
/
-Back To Top
************************************************************************
TITLE: T.11 - 2D
Precipitation Test with Flow
************************************************************************
/
GRID IS 22 BY 22
/
COORdinate y RANGE 1
/
FOR MATErial 1
MATErial particle DENSity: 1
MATErial POROsity: NE=0.50
HYDRaulic properties as=1.0E-03, Kx=0.1 m/yr, Ky=0.1
/
DENSITY 1 !Kg/m3
/
RENAME C = "Clqd - Liquid State Concentration of ANDRA Species"
RENAME C2 = "Csld - Solid State Concentration of ANDRA Species"
ALLOCATE "PREAC - Liquid-Solid Precipitation Reaction Rate"
/
TRANSPORT Clqd Kd=0.0, Dm=0.001 m-sq/yr
TRANSPORT Csld in SOLID phase
/
SET P = 0.
SET Clqd = 1.0
SET Csld = 0.0
/
BOUNdary FOR P: at X- VALUE = 1
BOUNdary FOR P: at X+ FLUX = 0
BOUNdary FOR P: at Y- FLUX = 0
SELECT (1,999) TO (12,999)
BOUNdary FOR P: at Y+ FLUX = 0 SELEcted
SELECT (13,999) TO (999,999)
BOUNdary FOR P: at Y+ VALUE = 0 SELEcted
/
BOUNdary FOR Clqd: at X- VALUE = 1
BOUNdary FOR Clqd: at X+ FLUX = 0
BOUNdary FOR Clqd: at Y- FLUX = 0
SELECT (1,999) TO (12,999)
BOUNdary FOR Clqd: at Y+ FLUX = 0 SELEcted
SELECT (13,999) TO (999,999)
BOUNdary FOR Clqd: at Y+ grad = 0 SELEcted !!outflow
/
BOUNdary FOR Csld: at X- VALUE = 0
BOUNdary FOR Csld: at X+ FLUX = 0
BOUNdary FOR Csld: at Y- FLUX = 0
SELECT (1,999) TO (12,999)
BOUNdary FOR Csld: at Y+ FLUX = 0 SELEcted
SELECT (13,999) TO (999,999)
BOUNdary FOR Csld: at Y+ grad = 0 SELEcted !!outflow
/
CONV Clqd 1.E-6 20
PRECipitation reaction PREAC: Clqd & Csld freq=1. Solublity = 0.2
DIFFusivity for Csld = 0.0001 !!! m^2/yr
/
DIAGnostic time dtime Clqd Csld BClqd (11,11) every 100 steps
/
DEBUG GEOMETRY OFF
SOLVE 1 YEARS IN STEPS OF 0.01
OUTPut Clqd Csld in NARRow NOW
/SAVE NOW 'prec.sav'
/
SOLVE 1 YEARS IN STEPS OF 0.01 1.005 0.02
/OUTPut Clqd Csld in NARRow NOW
/SAVE NOW 'prec.sav'
/
SOLVE 3 YEARS IN STEPS OF 0.02 1.005 0.05
/OUTPut Clqd Csld in NARRow NOW
/SAVE NOW 'prec.sav'
/
SOLVE 5 YEARS IN STEPS OF 0.05 1.005 0.10
/OUTPut Clqd Csld in NARRow NOW
/SAVE NOW 'prec.sav'
/
SOLVE 10 YEARS IN STEPS OF 0.10 1.005 0.20
/OUTPut Clqd Csld in NARRow NOW
/SAVE NOW 'prec.sav'
/
SOLVE 30 YEARS IN STEPS OF 0.20 1.005 0.50
/OUTPut Clqd Csld in NARRow NOW
/SAVE NOW 'prec.sav'
/
SOLVE 50 YEARS IN STEPS OF 0.50 1.005 1.0
/OUTPut Clqd Csld in NARRow NOW
/SAVE NOW 'prec.sav'
/
SOLVE 100 YEARS IN STEPS OF 1.0 1.005 2.0
/OUTPut Clqd Csld in NARRow NOW
/SAVE NOW 'prec.sav'
/
END
/
-Back To Top
************************************************************************
TITLE: T.12 - 2D
Precipitation Test with Flow
************************************************************************
/
GRID IS 3 BY 200
/
coor y nodes min 0.0 max 31.5
/
ALLOCATE "PREAC - Liquid-Solid Precipitation Reaction Rate"
/
SET V=3.16e-2
MATE 1 (1,1) to (3,2)
MATE 2 (1,3) to (3,200)
/
FOR MATErial 1
MATErial particle DENSity: 1
MATErial POROsity: NE=1.0e-3
TRANSPORT C Kd=0.0, Dm=0.0 m-sq/yr
/
FOR MATErial 2
MATErial particle DENSity: 1
MATErial POROsity: NE=1.0e-3
TRANSPORT C Kd=0.0, Dm=0. m-sq/yr 10.0
/
SELECT MATERIAL 1 as ID=MAT1
SELECT MATERIAL 2 as ID=MAT2
PRECipitation PREAC: C & C2 freq=1000. Solublity = 1.0 ID=MAT1
PRECipitation PREAC: C & C2 freq=1000. Solublity = 1.0e+10 ID=MAT2
TRANSPORT C2 in SOLID phase
/
SET C = 0.0 ID=MAT2
SET C2 = 10.0 ID=MAT1
/
BOUNdary FOR C: at Y- FLUX= 0.0
BOUNdary FOR C: at Y+ GRAD = 0.0
BOUNdary FOR C2: at Y+ FLUX = 0.0
/
CONV C 1.E-5 20
/
DIAGnostic time dtime C C2 BC (2,11) every 100 steps
/
DEBUG GEOMETRY OFF
SOLVE 1 YEARS IN STEPS OF 0.01
SELECT (2,1) (2,999)
OUTPut C C2 in NARRow NOW SELEcted
save C C2 now on 'PREC.SAV'
/
SOLVE 1 YEARS IN STEPS OF 0.01 1.005 0.02
OUTPut C C2 in NARRow NOW SELEcted
save C C2 now
/
SOLVE 3 YEARS IN STEPS OF 0.02 1.005 0.05
OUTPut C C2 in NARRow NOW SELEcted
save C C2 now
/
SOLVE 5 YEARS IN STEPS OF 0.05 1.005 0.10
OUTPut C C2 in NARRow NOW SELEcted
save C C2 now
/
SOLVE 10 YEARS IN STEPS OF 0.10 1.005 0.20
OUTPut C C2 in NARRow NOW SELEcted
save C C2 now
/
SOLVE 30 YEARS IN STEPS OF 0.20 1.005 0.50
OUTPut C C2 in NARRow NOW SELEcted
save C C2 now
/
SOLVE 50 YEARS IN STEPS OF 0.50 1.005 1.0
OUTPut C C2 in NARRow NOW SELEcted
save C C2 now
/
SOLVE 100 YEARS IN STEPS OF 1.0 1.005 2.0
OUTPut C C2 in NARRow NOW SELEcted
save C C2 now
/
END
/
-Back To Top
************************************************************************
TITLE: T.13 -
Precipitation Test - Cs function of variable
************************************************************************
GRID IS 3 BY 3
coor y nodes min 0.0 max 1.0
/
ALLOcate VAR1
ALLOCATE "PREAC - Liquid-Solid Precipitation Reaction Rate"
/
MATErial particle DENSity: 1
MATErial POROsity: 0.5
TRANSPORT C Kd=0.0, Dm=0.0 m-sq/yr
/
Set VAR1=0.1
PRECipitation PREAC: C & C2 Solublity Linear VAR1 freq 1000.0
/PRECipitation PREAC: C & C2 Solublity freq 1000.0 cslimit=0.1
TRANSPORT C2 in SOLID phase
/
SET C = 1.0
SET C2 = 0.0
/
CONV C C2 1.E-20 20 20
/
DIAGnostic time dtime C C2 VAR1 (2,2) every 1 iterations
DEBUG GEOMETRY OFF
SOLVE 1.0 0.01
/
END
/
-Back To Top
************************************************************************
TITLE: T.14 - TESTING
HISTORY OPTIONS
************************************************************************
GRID 22 by 12
/
COORdinate X faces RANGE 400
COORdinate Y faces RANGE 100
/
HYDRaulic props S=0.01, kx=1, ky=1
ROCK POROsity 0.01
TRANsport kd=0, Dm=0.1 m*m/yr, al= 10 m, at= 2 m
/
LOCATE (01,07) TO (99,99)
SET P LINEAR 0.4 -0.001 * X in SELEcted zone
/
LOCATE (01,01) TO (01,06) ID=LOBL
LOCATE (99,01) TO (99,06) ID=LOBR
LOCATE (01,07) TO (01,99) ID=HIBL
LOCATE (99,07) TO (99,99) ID=HIBR
/
BOUNdary cond for P: index X-, FLUX =0.0 ID=LOBL
BOUNdary cond for P: index X+, FLUX =0.0 ID=LOBR
BOUNdary cond for P: index Y-, FLUX =0.0
BOUNdary cond for P: index Y+, FLUX =0.0
/
SET C = 35.0
BOUNdary cond for C: X- 0.0 ID=HIBL
BOUNdary cond for C: X+ GRAD=0.0 ID=HIBR
BOUNdary cond for C: index X-, FLUX =0.0 ID=LOBL
BOUNdary cond for C: index X+, FLUX =0.0 ID=LOBR
BOUNdary cond for C: index Y-, FLUX =0.0
BOUNdary cond for C: index Y+, FLUX =0.0
/
DENSity rho=1000
/
HISTORY at 'TEST.HIS' time=0.5 tables
HISTORY at (6,2) (6,10) (16,2)
/
DIAGnostic node at TIME U, V P C (11,6) every 10 steps
CONVergence P 1.0E-15 iteration 2500
CONVergence REFE C 1.0E-3 10
/
DEBUG GEOMETRY OFF
SOLVE 100 0.01 1.005 1
/
END
/
-Back To Top
************************************************************************
TITLE: T.15 - TESTING
HISTORY OPTION: INITIAL READ FROM STEADY
************************************************************************
/
GRID 12 BY 12
COORDINATE X RANGE = 10
COORDINATE Y RANGE = 10
/
HYDRAULIC properties: ss=0.001, k = 0.097 0.0555
TRANSport kd=0., dm=0.001
BOUN P at X- 1.
BOUN P at X+ FLUX = 0.
BOUN P at Y- FLUX = 0.
BOUN P at Y+ VALUE = 0.
/
BOUN C at X- 1.
BOUN C at X+ FLUX = 0.
BOUN C at Y- FLUX = 0.
BOUN C at Y+ GRAD = 0.
/
DIAgnostics U, C at: (6,6) every 20 steps
output transport for C
/
DEBUG GEOMETRY OFF
SOLVE 100 STEADY !! Active for 1st run
SAVE 'test.sav' !! Active for 1st run
/
select 4, 2
/SOURce FLOW at Q=0.1 C = 5 !! Active for 2nd run
/read 'test.sav' !! Active for 2nd run
history on 'test.his' TABL PLOTs !! Active for 2nd run
history at
(4,3) (5,2) (3,4) (4,3) (2,5) (3,5) (5,3) (3,3)
(2,2) (4,4) (5,5) (6,6) (7,7) (9,9) (10,10) (11,11)
(1,1) (4,2) (5,6) (6,7) (8,7) (9,10) (10,10) (11,11)
(12,12) (11,12)
!! Active for 2nd run
/
SOLVE 1 0.01 !! Active for 2nd run
output u v kx ky
/
END
/
-Back To Top
************************************************************************
TITLE: T.16 - TESTING
RADIATION TYPE OF SOURCE
************************************************************************
/
GRID 12 BY 12
COORDINATE X RANGE = 10
COORDINATE Y RANGE = 10
HYDRAULIC properties: ss=0.001, k = 3*1.
TRANSport kd=0., dm=0.
SET P LINEar function 0. -1. X
/
INITial C = 0.
SELEct (6,6)
SOURce for C: RADIATION heq=0.001 feq = 0.2 positive selected
/
ALLOCATE Hcoeff
SET Hcoeff LINEAR 0.1 +10. * C
SOURce for C: RADIATION VARIABLE COEFF hcoeff feq = 0.2 positive
selected
/
ALLOCATE Cbase
SET Cbase LINEAR 0.1 +1 * C
SOURce for C: RADIATION HEQ=0.001 Cbase positive selected
/
SOURce for C: RADIATION VARIABLE COEFF hcoeff Cbase positive
selected
/
DIAgnostics U, C at: (6,6) every 20 steps
DEBUG GEOMETRY OFF
SOLVE P AND C ONLY 100 secs in steps of 1
/
END
/
-Back To Top
************************************************************************
TITLE: T.17 - HISTORY
STORAGE AND HISTORY SOURCE OPTIONS
************************************************************************
/
GRID 12 BY 12
COORDINATE X RANGE = 10
COORDINATE Y RANGE = 10
HYDRAULIC properties: ss=0.001, k = 3*1.
TRANSport kd=0., dm=0.
SET P LINEar function 0. -1. X
/
INITial C = 0.
SELEct (6,6)
SOURce for C: RADIATION heq=0.001 feq = 0.2 positive selected
/
ALLOCATE Hcoeff
SET Hcoeff LINEAR 0.1 +10. * C
SOURce for C: RADIATION VARIABLE COEFF hcoeff feq = 0.2 positive
selected
/
ALLOCATE Cbase
SET Cbase LINEAR 0.1 +1 * C
SOURce for C: RADIATION HEQ=0.001 Cbase positive selected
/
SOURce for C: RADIATION VARIABLE COEFF hcoeff Cbase positive
selected
/
HISTory STORage for P EVERY TIME= 7 units
HISTory STORage for C EVERY TIME=11 units for SELEcted region
HISTory of SOURCE for P EVERY TIME=2.5 for SELEcted region
HISTory of SOURCE for c every step
/
DIAgnostics U, C at: (6,6) every 20 steps
DEBUG GEOMETRY OFF
/
SELECT FROM (1,1) TO (99,99) (2,2)
OUTPUT FOR SELECTED REGION IN NARROW MODE
/
SOLVE P AND C ONLY 100 secs in steps of 1
/
END
/
-Back To Top
************************************************************************
TITLe T.18: Contaminant
Transport in a 3-Layer Aquifer System
/ Demonstration Example for PORFLOW Workshop
************************************************************************
GRID 50 by 22
Coordinate X minimum 0 maximum 2400 meters
Coordinate Y minimum 0 maximum 200 meters
/
*==================================================================
// SOIL ZONES
*==================================================================
Zone 1 everywhere !!! Upper unconfined aquifer
Zone 2 from ( 1,10) to (999,14) !!! Middle aquitard layer
Zone 3 from ( 1, 1) to (999, 9) !!! Lower confined aquifer
*==================================================================
// SOIL PROPERTIES
*==================================================================
FOR zone 1
ROCK density = 1600
ROCK porosity = 0.40
HYDRaulic Properties: Ss=0.01, Kx = 315, Ky = 315 m/yr
TRANsport Kd = 0., Dm = 5., Alpha_L = 5, Alpha_T = 1
/
FOR zone 2
ROCK density = 2000
ROCK porosity = 0.20
HYDRaulic Properties: Ss = 1.E-6, Kx = 3.15E-2, Ky = 3.15E-3
TRANsport Kd = 0., Dm = 5 m/yr^2
/
FOR zone 3
ROCK density = 1800
ROCK porosity = 0.30
HYDRaulic Properties: Ss=0.001, Kx = 157.5, Ky = 52.5
TRANsport Kd = 0., Dm = 5, Alpha_L = 1.0, Alpha_T = 0.2
*==================================================================
// INITIAL & BOUNDARY CONDITIONS
*==================================================================
SET LINEar function: P = 0. -0.005 * X
BOUNdary for P: index X- (left): GRAD = 0.005 !Regional gradient
BOUNdary for P: index Y+ (top): FLUX = 0.2 !infiltration m/yr
BOUNdary for P: index Y- (bottom): FLUX = 0.
/
SET C = 0. everywhere
BOUNdary for C: index X- (left): Value = 0. !Clean inflow
BOUNdary for C: index X+ (right): GRADient = 0. !Unimpeded outflow
BOUNdary for C: index Y-, (bottom): FLUX = 0. !No flux
*==================================================================
// SOURCE TERM
*==================================================================
LOCAte source region at (17,20) with ID=SOR1
SOURce for C: at ID=SOR1 strength = 1 in relative units
*==================================================================
// DIAGNOSTIC, OUTPUT & OPERATIONAL CONTROL
*==================================================================
DIAGnostic node (21,15) output every 20 steps
HISTory at (21,19), (21,13) (21,7)
HISTory for C on file 'T18.HIS'
TRACk particle starting at: (1000,180) on file 'ANST18.TRK'
TRACk particle starting at: (1000,120)
TRACk particle starting at: (1000,60) TABLe also
/
SELEct from (1,1) to (99,99) (5,3)
OUTPUT U C for SELEcted subregion
CONVergence REFErence for P: epsilon=0.001 max iterations=50
SOLVe for 100 years AUTO !0.01 1.01 0.1
SAVE NOW on file 'ANST18.SAV'
SOLVe for 400 years
/
END
/
-Back To Top
*******************************************************************
TITLe T.19: Contaminant
Transport in Vadose Zone of an Aquifer
/ Demonstration Example for PORFLOW Workshop
*******************************************************************
GRID 32 by 33
COORdinate X range 150
COORdinate Y RANGE 35
/
GRAVity 0., -1 (normalized; absolute value not required for
single-phase)
/
ROCK density 1800.0
ROCK porosity = 0.1875, 0.1875, 0.1875
HYDRaulic properties: ss = 1.E-7, (kx,ky,kz) --> 3*750
TRANsport properties for C: kd=0., dm=0.001, alpha_l=1, alpha_t=0.1
PROPerty averages by GEOMetric mean option
/
MULTiphase VAN Genuchten GENEral: n=2, alph=0.2, sr=0.25, sg=0.
m=1
MULTiphase CONDuctivity POWEr law: n=4, A=1., B=0, C=0, sr=0.25
/
INITial P -10. EVERywhere
INITIAL C = 0.
/
SELEct subdomain (1,1) to (1,17)
BOUNdary P: index X- LINEar function: 21. -1. * Y for SELEcted
segment
SELEct subdomain (1,18) to (1,99)
BOUNdary P: index X- FLUX = 0. for SELEcted segment
/
SELEct subdomain from (32,1) to (32,13)
BOUNdary P: index X+ LINEar function: 19. -1. * Y for SELEcted
segment
SELEct subdomain from (32,14) to (32,99)
BOUNdary P: index X+ FLUX = 0. for SELEcted segment
/
BOUNdary P: index Y+ FLUX = 1. m/yr
BOUNdary P: index Y- FLUX = 0.
/
BOUNdary for C: index X-, GRAD = 0.
BOUNdary for C: index X+, GRAD = 0.
/BOUNdary for C: index Y-, FLUX = 0.
/BOUNdary for C: index Y+, FLUX = 0.
/
SELEct region (17,31) to (18,31)
SOURce for C strength = 1000. in SELEcted region.
/
TRACk particle starting at: ( 5.,33.) dx=1000; dt=0.1 on
'ANST19.TRK'
TRACk particle starting at: (18.,32.) print TABLes
/
CONVergence REFErence variable P: local mode; epsilon=1.0E-4 10
/
DIAGnostic output for node (9, 25) every 1 steps
/
SAVE U, V, P, C, S and H on file 'ANST19.SAV'
SOLVE in STEADy state mode for 1000 steps
/SOLVE in 100 0.001 1.01 0.01 !AUTO steps
/
SELEct from (1,1) to (99,99) (4,3)
OUTPut U,V,C, P S and MOISture content for SELECted zone
/
END
/
-Back To Top
************************************************************************
TITLE: T.20: Serie #3,
CAS #2 - Precipitation with Transport & Filiation
************************************************************************
***** DEVELOPPEMENTS DU CODE PORFLOW POUR l'ANDRA
***** Introduction du modele de precipitation par famille
d'isotopes
***** Cas extrait des s‚ries de test de validation du
developpement
***** ACRI, mars 2000
************************************************************************
GRID NODEs is 38
/
COORdinate X corners :
0.00 0.01 0.02 0.03 0.04
0.05 0.06 0.07 0.08 0.09
0.10 0.12 0.14 0.17 0.20
0.25 0.30 0.40 0.50 0.60
0.70 0.75 0.80 0.83 0.86
0.88 0.90 0.91 0.92 0.93
0.94 0.95 0.96 0.97 0.98
0.99 1.00
/
FOR MATErial 1
MATErial dry solid DENSity = 2500 !kg/m3
MATErial POROsite effective = 0.12
HYDRaulic properties as=1.0E-07, Kx=1 m/an, Ky=1 m/an
/
DENSity reference eau 1000 !kg/m3
/
RENAME C = "C1L - Concentration premier isotope en phase liquide"
RENAME C2 = "C1S - Concentration premier isotope en phase solide"
RENAME C3 = "C2L - Concentration deuxieme isotope en phase liquide"
RENAME C4 = "C2S - Concentration deuxieme isotope en phase solide"
ALLOCATE "C3L - Concentration troisieme isotope en phase liquide"
ALLOCATE "C3S - Concentration troisieme isotope en phase solide"
/
ALLOCATE CLIQUID
ALLOCATE CTOT
ALLOCATE CSOLID
SET CLIQUID SUM OF C1L + C2L + C3L ALWAYS
SET CSOLID SUM OF C1S + C2S + C3S ALWAYS
SET CTOT SUM OF CLIQUID + CSOLID ALWAYS
/
TRANSPORT C1L Kd=0.0, Dm=3.1536e-5 !m2/year
TRANSPORT C1S in SOLID phase
TRANSPORT C2L Kd=0.0, Dm=3.1536e-5 !m2/year
TRANSPORT C2S in SOLID phase
TRANSPORT C3L Kd=0.0, Dm=3.1536e-5 !m2/year
TRANSPORT C3S in SOLID phase
/
BOUNdary FOR C1L: at X- VALUE = 1
BOUNdary FOR C1L: at X+ VALUE = 1
BOUNdary FOR C1S: at X- FLUX = 0
BOUNdary FOR C1S: at X+ FLUX = 0
BOUNdary FOR C2L: at X- VALUE = 1
BOUNdary FOR C2L: at X+ VALUE = 1
BOUNdary FOR C2S: at X- FLUX = 0
BOUNdary FOR C2S: at X+ FLUX = 0
BOUNdary FOR C3L: at X- VALUE = 0
BOUNdary FOR C3L: at X+ VALUE = 0
BOUNdary FOR C3S: at X- FLUX = 0
BOUNdary FOR C3S: at X+ FLUX = 0
/
BOUNdary FOR C1L: at Y- FLUX = 0
BOUNdary FOR C1L: at Y+ FLUX = 0
BOUNdary FOR C1S: at Y- FLUX = 0
BOUNdary FOR C1S: at Y+ FLUX = 0
BOUNdary FOR C2L: at Y- FLUX = 0
BOUNdary FOR C2L: at Y+ FLUX = 0
BOUNdary FOR C2S: at Y- FLUX = 0
BOUNdary FOR C2S: at Y+ FLUX = 0
BOUNdary FOR C3L: at Y- FLUX = 0
BOUNdary FOR C3L: at Y+ FLUX = 0
BOUNdary FOR C3S: at Y- FLUX = 0
BOUNdary FOR C3S: at Y+ FLUX = 0
/
SET C1L = 0.0
SET C1S = 0.0
SET C2L = 0.0
SET C2S = 0.0
SET C3L = 0.0
SET C3S = 0.0
/
ALLOCATE "CSUM - Somme des concentrations isotopes pour limite de
solubilite"
SET CSUM = SUM MASS of C1L C1S C2L C2S C3L C3S ALWAys
/
ALLOCATE "R1 - Reaction de precipitation - dissolution premier
isotope"
ALLOCATE "R2 - Reaction de precipitation - dissolution deuxieme
isotope"
ALLOCATE "R3 - Reaction de precipitation - dissolution troisieme
isotope"
PRECipitation R1: C1L & C1S freq = 1.E8 Solubility = 1.00 ELEMent
CSUM
PRECipitation R2: C2L & C2S freq = 1.E8 Solubility = 1.00 ELEMent
CSUM
PRECipitation R3: C3L & C3S freq = 1.E8 Solubility = 1.00 ELEMent
CSUM
/
DECAy rate for C2L with half LIFE of 1.0
REGEneration of C3L from decay of C2L factor 1.0
DECAy rate for C2S with half LIFE of 1.0
REGEneration of C3S from decay of C2S factor 1.0
/
LIMIT C1L MINImum 0. MAXImum 1.e50
LIMIT C2L MINImum 0. MAXImum 1.e50
LIMIT C3L MINImum 0. MAXImum 1.e50
LIMIT C1S MINImum 0. MAXImum 1.e50
LIMIT C2S MINImum 0. MAXImum 1.e50
LIMIT C3S MINImum 0. MAXImum 1.e50
/
MATRIX C1L C2L C3L C1S C2S C3S ADI X ONLY
/
CONV for SPECies in COUPled mode with a maximum of 3 iterations
CONV C1L 1.E-16 5
CONV C1S 1.E-16 5
CONV C2L 1.E-16 5
CONV C2S 1.E-16 5
CONV C3L 1.E-16 5
CONV C3S 1.E-16 5
/
DIAGnostic at (15,2) TIME CTOT CLIQ C1L C2L C3L every 1000 steps
DEBUg GEOMETRY OFF
/
SELECT from (1,2) to (12,2)
OUTPut C1L C2L C3L C1S C2S C3S CLIQ CSOL CTOT narrow selected
/
/ Diffusion criteria: DT < (DX*DX)/(2*Dm)
/ here: DT < 15854*DX*DX = 1.58E+00 years
/
SOLVE 1 years step 0.0005 1.00025 0.00175
/SAVE C1L C2L C3L CLIQ C1S C2S C3S CSOL CTOT file 'PORT20.SAV'
every 500 steps
/SOLVE 5 years step 0.0005 1.00025 0.00175 ! for full simulation
/
/SAVE RESTART
SAVE OFF
/
END
-Back To Top
1
2
3
4
5
6
7
8
9
10
# Back to PORFLOW EXAMPLES DATABASE TOC
Related Links:
# PORFLOW Applications
# PORFLOW Express
# PORFLOW Publications
# PORFLOW Users
# PORFLOW Price List
# Request CFDStudio/PORFLOW Demo
# CFDStudio/PORFLOW Tutorials
# PORFLOW Manual
# PORFLOW Validation Report
