#Appendix A (Part 1 Input File)#

#-----------------------------------------------------------
~Simulation Title Card
#-----------------------------------------------------------
1,
STOMP Tutorial Problem 6a,
Mart Oostrom/Mark White,
PNNL,
June 03,
15:00,
1,
Simulation of NAPL spills in 2D domain,


#-----------------------------------------------------------
~Solution Control Card
#-----------------------------------------------------------
Normal,
Water-Oil,
2,
0,s,10,d,1,s,1,d,1.25,8,1.e-6,
10,d,200,d,1,s,10,d,1.25,8,1.e-6,
10000,
Variable Aqueous Diffusion,
,


#-----------------------------------------------------------
~Grid Card
#-----------------------------------------------------------
Uniform Cartesian,
25,1,20,
2,m,
1,m,
0.5,m,


#-----------------------------------------------------------
~Rock/Soil Zonation Card
#-----------------------------------------------------------
1,
Sand,1,25,1,1,1,20,


#-----------------------------------------------------------
~Mechanical Properties Card
#-----------------------------------------------------------
Sand,2650,kg/m^3,0.43,0.43,,,Millington and Quirk,


 
#-----------------------------------------------------------
~Hydraulic Properties Card
#-----------------------------------------------------------
Sand,10,hc m/day,,,10,hc m/day,


#-----------------------------------------------------------
~Saturation Function Card
#-----------------------------------------------------------
72.0,dynes/cm,,,32,dynes/cm,
Sand,Van Genuchten,0.1,1/cm,2.0,0.10,72.0,dynes/cm,,


#-----------------------------------------------------------
~Aqueous Relative Permeability Card
#-----------------------------------------------------------
Sand,Mualem,,


#-----------------------------------------------------------
~NAPL Relative Permeability Card
#-----------------------------------------------------------
Sand,Mualem,,


#-----------------------------------------------------------
~Oil Properties Card
#-----------------------------------------------------------
NAPL,
165.834,g/mol,251.,K,394.4,K,620.2,K,
47.6,bar,289.6,cm^3/mol,0.2758,0.2515,0.0,debyes,
-1.431e+1,5.506e-1,-4.513e-4,1.429e-7,
Equation 1,-7.36067,1.82732,-3.47735,-1.00033,
Constant,800.0,kg/m^3,
Constant,0.002,Pa s,
1.0e8,Pa,


#-----------------------------------------------------------
~Initial Conditions Card
#-----------------------------------------------------------
2,
Aqueous Pressure,140000,Pa,0.0,1/m,,,-9793.5192,1/m,1,25,1,1,1,20,
NAPL Pressure,-1.e9,Pa,,,,,,,1,25,1,1,1,20,


#-----------------------------------------------------------
~Boundary Conditions Card
#-----------------------------------------------------------
3,
top,zero flux,neumann,
13,13,1,1,20,20,2,
0,d,-1.e9,Pa,0.0,-0.25,m/day,
10.0,d,-1.e9,Pa,0.0,-0.25,m/day,
west,hydraulic gradient,dirichlet,
1,1,1,1,1,20,1,
0,d,140000,Pa,0.0,-1.e9,Pa,
east,hydraulic gradient,dirichlet,
25,25,1,1,1,20,1,
0,d,140000,Pa,0.0,-1.e9,Pa,


#-----------------------------------------------------------
~Output Options Card
#-----------------------------------------------------------
9,
13,1,20,
13,1,15,
13,1,10,
13,1,8,
13,1,6,
13,1,4,
1,1,7,
1,1,6,
1,1,4,
1,1,day,m,6,6,6,
2,
napl saturation,,
aqueous saturation,,
2,
10,d,
100,d,
3,
no restart,,
napl saturation,,
aqueous saturation,,


#-----------------------------------------------------------
~Surface Flux Card
#-----------------------------------------------------------
3,
NAPL volumetric flux,m^3/day,m^3,top,13,13,1,1,20,20,
NAPL volumetric flux,m^3/day,m^3,west,1,1,1,1,1,20,
NAPL volumetric flux,m^3/day,m^3,east,25,25,1,1,1,20,



#Appendix B (Part 2 Input File)#


#------------------------------------------------------------------
~Simulation Title Card
#------------------------------------------------------------------
1,
STOMP Tutorial Problem 6b (steady-state calculations),
Mart Oostrom/Mark White,
PNNL,
June 03,
15:00,
2,
Steady state calculations to set up flow field,
Restart file will be used for subsequential NAPL infiltration simulation,


#------------------------------------------------------------------
~Solution Control Card
#------------------------------------------------------------------
Normal,
Water-Oil,
1,
0,s,1000,yr,1,d,1000,yr,1.25,8,1.e-6,
10000,
Variable Aqueous Diffusion,
,


#------------------------------------------------------------------
~Grid Card
#------------------------------------------------------------------
Uniform Cartesian,
25,1,20,
2,m,
1,m,
0.5,m,


#------------------------------------------------------------------
~Rock/Soil Zonation Card
#------------------------------------------------------------------
1,
Sand,1,25,1,1,1,20,


#------------------------------------------------------------------
~Mechanical Properties Card
#------------------------------------------------------------------
Sand,2650,kg/m^3,0.43,0.43,,,Millington and Quirk,


 
#------------------------------------------------------------------
~Hydraulic Properties Card
#------------------------------------------------------------------
Sand,10,hc m/day,,,10,hc m/day,


#------------------------------------------------------------------
~Saturation Function Card
#------------------------------------------------------------------
72.0,dynes/cm,,,32,dynes/cm,
Sand,Van Genuchten,0.1,1/cm,2.0,0.10,72.0,dynes/cm,,


#------------------------------------------------------------------
~Aqueous Relative Permeability Card
#------------------------------------------------------------------
Sand,Mualem,,


#------------------------------------------------------------------
~NAPL Relative Permeability Card
#------------------------------------------------------------------
Sand,Mualem,,


#------------------------------------------------------------------
~Oil Properties Card
#------------------------------------------------------------------
NAPL,
165.834,g/mol,251.,K,394.4,K,620.2,K,
47.6,bar,289.6,cm^3/mol,0.2758,0.2515,0.0,debyes,
-1.431e+1,5.506e-1,-4.513e-4,1.429e-7,
Equation 1,-7.36067,1.82732,-3.47735,-1.00033,
Constant,800.0,kg/m^3,
Constant,0.002,Pa s,
1.0e8,Pa,


#------------------------------------------------------------------
~Initial Conditions Card
#------------------------------------------------------------------
2,
Aqueous Pressure,140000,Pa,-40.0,1/m,,,-9793.5192,1/m,1,25,1,1,1,20,
NAPL Pressure,-1.e9,Pa,,,,,,,1,25,1,1,1,20,


#------------------------------------------------------------------
~Boundary Conditions Card
#------------------------------------------------------------------
3,
top,neumann,zero flux,
1,25,1,1,20,20,1,
0,d,-100.0,mm/year,0.0,-1.e9,Pa,
west,hydraulic gradient,dirichlet,
1,1,1,1,1,10,1,
0,d,140040,Pa,0.0,-1.e9,Pa,
east,hydraulic gradient,dirichlet,
25,25,1,1,1,10,1,
0,d,138040,Pa,0.0,-1.e9,Pa,


#------------------------------------------------------------------
~Output Options Card
#------------------------------------------------------------------
9,
1,1,20,
25,1,20,
13,1,20,
13,1,18,
13,1,16,
13,1,14,
13,1,12,
13,1,10,
13,1,8,
1,1,yr,m,6,6,6,
3,
napl saturation,,
aqueous saturation,,
z aqueous volumetric flux,mm/year,
0,
2,
aqueous saturation,,
z aqueous volumetric flux,mm/year,




#Appendix C (Part 3 Input File)#


#----------------------------------------------------
~Simulation Title Card
#----------------------------------------------------
1,
STOMP Tutorial Problem 6b (NAPL Infiltration),
Mart Oostrom/Mark White,
PNNL,
June 03,
15:00,
2,
Simulation of NAPL spill in 2D domain,
Simulation starts with Restart file,


#----------------------------------------------------
~Solution Control Card
#----------------------------------------------------
Restart,
Water-Oil,
2,
0,s,10,d,1,s,1,d,1.25,8,1.e-6,
10,d,200,d,1,s,10,d,1.25,8,1.e-6,
10000,
Variable Aqueous Diffusion,
,


#----------------------------------------------------
~Grid Card
#----------------------------------------------------
Uniform Cartesian,
25,1,20,
2,m,
1,m,
0.5,m,


#----------------------------------------------------
~Rock/Soil Zonation Card
#----------------------------------------------------
1,
Sand,1,25,1,1,1,20,


#----------------------------------------------------
~Mechanical Properties Card
#----------------------------------------------------
Sand,2650,kg/m^3,0.43,0.43,,,Millington and Quirk,


 
#----------------------------------------------------
~Hydraulic Properties Card
#----------------------------------------------------
Sand,10,hc m/day,,,10,hc m/day,


#----------------------------------------------------
~Saturation Function Card
#----------------------------------------------------
72.0,dynes/cm,,,32,dynes/cm,
Sand,Van Genuchten,0.1,1/cm,2.0,0.10,72.0,dynes/cm,,


#----------------------------------------------------
~Aqueous Relative Permeability Card
#----------------------------------------------------
Sand,Mualem,,


#----------------------------------------------------
~NAPL Relative Permeability Card
#----------------------------------------------------
Sand,Mualem,,


#----------------------------------------------------
~Oil Properties Card
#----------------------------------------------------
NAPL,
165.834,g/mol,251.,K,394.4,K,620.2,K,
47.6,bar,289.6,cm^3/mol,0.2758,0.2515,0.0,debyes,
-1.431e+1,5.506e-1,-4.513e-4,1.429e-7,
Equation 1,-7.36067,1.82732,-3.47735,-1.00033,
Constant,800.0,kg/m^3,
Constant,0.002,Pa s,
1.0e8,Pa,


#----------------------------------------------------
~Initial Conditions Card
#----------------------------------------------------
0,


#----------------------------------------------------
~Boundary Conditions Card
#----------------------------------------------------
5,
top,neumann,neumann,
13,13,1,1,20,20,4,
0,d,-100,mm/year,0.0,-0.25,m/day,
10.0,d,-100,mm/year,0.0,-0.25,m/day,
10.0,d,-100,mm/year,0.0,0.,m/day,
100.0,d,-100,mm/year,0.0,0.,m/day,
top,neumann,zero flux,
1,12,1,1,20,20,1,
0,d,-100.0,mm/year,0.0,-1.e9,Pa,
top,neumann,zero flux,
14,25,1,1,20,20,1,
0,d,-100.0,mm/year,0.0,-1.e9,Pa,
west,hydraulic gradient,dirichlet,
1,1,1,1,1,10,1,
0,d,140040,Pa,0.0,-1.e9,Pa,
east,hydraulic gradient,dirichlet,
25,25,1,1,1,10,1,
0,d,138040,Pa,0.0,-1.e9,Pa,


#----------------------------------------------------
~Output Options Card
#----------------------------------------------------
9,
13,1,20,
13,1,15,
13,1,10,
13,1,8,
13,1,6,
13,1,4,
1,1,20,
1,1,15,
1,1,10,
1,1,day,m,6,6,6,
3,
napl saturation,,
aqueous saturation,,
z aqueous volumetric flux,mm/year,
2,
10,d,
100,d,
3,
no restart,,
napl saturation,,
aqueous saturation,,


#----------------------------------------------------
~Surface Flux Card
#----------------------------------------------------
3,
NAPL volumetric flux,m^3/day,m^3,top,13,13,1,1,20,20,
NAPL volumetric flux,m^3/day,m^3,west,1,1,1,1,1,20,
NAPL volumetric flux,m^3/day,m^3,east,25,25,1,1,1,20,