#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,