#------------------------------------------------------------------------
~Simulation Title Card
#------------------------------------------------------------------------
1,
STOMP Tutorial Problem 8,
Mart Oostrom/Mark White,
PNNL,
June 03,
15:00,
4,
Simulation of 1D partitioning tracer experiment,
Partitioning tracer test for detection, estimation, and,
assessment of subsurface nonaqueous phase liquids,
Water Resources Research, Vol. 31, No. 5, Pages 1201-1211, May 1995,


#------------------------------------------------------------------------
~Solution Control Card
#------------------------------------------------------------------------
Normal,
Water-Oil w/ Transport,
1,
0,s,2400,min,0.1,s,5,min,1.25,8,1.e-6,
10000,
Variable Aqueous Diffusion,
,


#------------------------------------------------------------------------
~Grid Card
#------------------------------------------------------------------------
Uniform Cartesian,
1,1,50,
1.958562,cm,
1.958562,cm,
0.61,cm,


#------------------------------------------------------------------------
~Rock/Soil Zonation Card
#------------------------------------------------------------------------
1,
Ottawa sand,1,1,1,1,1,50,


#------------------------------------------------------------------------
~Mechanical Properties Card
#------------------------------------------------------------------------
Ottawa sand,2650,kg/m^3,0.362,0.362,,,Millington and Quirk,


#------------------------------------------------------------------------
~Hydraulic Properties Card
#------------------------------------------------------------------------
Ottawa sand,,,,,1.5e-4,hc m/s,


#------------------------------------------------------------------------
~Saturation Function Card
#------------------------------------------------------------------------
72.0,dynes/cm,,,35.43,dynes/cm,
Ottawa sand,Entrapment Van Genuchten,2.7,1/m,1.23,.26,72.0,dynes/cm,,0.3,,,,


#------------------------------------------------------------------------
~Aqueous Relative Permeability Card
#------------------------------------------------------------------------
#Rel perm estimated from Lenhard and Parker 1987 Eq. (5)
Ottawa sand,Constant,0.8,


#------------------------------------------------------------------------
~NAPL Relative Permeability Card
#------------------------------------------------------------------------
Ottawa sand,Mualem,0.5,


#------------------------------------------------------------------------
~Solute/Fluid Interactions Card
#------------------------------------------------------------------------
3,
Tritium,1.0e-10,m^2/s,1.0e-10,m^2/s,
Linear Isotherm,0.0,
1.e+8,yr,
IPA,1.0e-10,m^2/s,1.0e-10,m^2/s,
Linear Isotherm,0.04,
1.e+8,yr,
DMB,1.0e-10,m^2/s,1.0e-10,m^2/s,
Linear Isotherm,2.76,
1.e+8,yr,
0,


#------------------------------------------------------------------------
~Dissolved Oil Transport Card
#------------------------------------------------------------------------
Ottawa sand,0.17,cm,,cm,linear kd,0.0,m^3/kg,


#------------------------------------------------------------------------
~Solute/Porous Media Interactions Card
#------------------------------------------------------------------------
Ottawa sand,0.17,cm,,cm,
tritium,Linear Isotherm,0.0,
IPA,Linear Isotherm,0.0,
DMB,Linear Isotherm,0.0,


#------------------------------------------------------------------------
~Oil Properties Card
#------------------------------------------------------------------------
PCE,
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,
Modified Rackett,0.2758,0.2515,
Constant,0.97e-3,Pa s,
9.463e7,Pa,


#------------------------------------------------------------------------
~Initial Conditions Card
#------------------------------------------------------------------------
3,
Aqueous Pressure,104282.15,Pa,,,,,-9793.5192,1/m,1,1,1,1,1,50,
NAPL Pressure,-1.e9,Pa,,,,,,,1,1,1,1,1,50,
Trapped NAPL Saturation,0.20,,,,,,,,1,1,1,1,1,50,


#------------------------------------------------------------------------
~Boundary Conditions Card
#------------------------------------------------------------------------
2,
top,neumann,zero flux,inflow aqueous,inflow aqueous,inflow aqueous,
1,1,1,1,50,50,4,
0,d,-0.013035,cm/min,0.0,-1.e9,Pa,1.0,1/cm^3,1.0,1/cm^3,1.0,1/cm^3,
84.705,min,-0.013035,cm/min,0.0,-1.e9,Pa,1.0,1/cm^3,1.0,1/cm^3,1.0,1/cm^3,
84.705,min,-0.013035,cm/min,0.0,-1.e9,Pa,0.0,1/cm^3,0.0,1/cm^3,0.0,1/cm^3,
2400.,min,-0.013035,cm/min,0.0,-1.e9,Pa,0.0,1/cm^3,0.0,1/cm^3,0.0,1/cm^3,
bottom,dirichlet,zero flux,outflow,outflow,outflow,
1,1,1,1,1,1,1,
0,d,104312.02,Pa,0.0,-1.e9,Pa,,,,,,,


#------------------------------------------------------------------------
~Output Options Card
#------------------------------------------------------------------------
3,
1,1,1,
1,1,25,
1,1,50,
1,1,min,cm,6,6,6,
6,
napl saturation,,
aqueous courant number,,
z aqueous volumetric flux,cm/min,
solute aqueous concentration,tritium,1/cm^3,
solute aqueous concentration,IPA,1/cm^3,
solute aqueous concentration,DMB,1/cm^3,
4,
84.705,min,
6,hr,
12,hr,
24,hr,
5,
no restart,,
napl saturation,,
solute aqueous concentration,tritium,1/cm^3,
solute aqueous concentration,IPA,1/cm^3,
solute aqueous concentration,DMB,1/cm^3,


#------------------------------------------------------------------------
~Surface Flux Card
#------------------------------------------------------------------------
5,
solute flux,tritium,1/min,,top,1,1,1,1,50,50,
solute flux,tritium,1/min,,bottom,1,1,1,1,1,1,
solute flux,ipa,1/min,,bottom,1,1,1,1,1,1,
solute flux,dmb,1/min,,bottom,1,1,1,1,1,1,
aqueous volumetric flux,ml/min,ml,top,1,1,1,1,50,50,