~Simulation Title Card
1,
STOMP Example Problem CO2-1 Zero Salinity,
M.D. White,
Pacific Northwest Laboratory,
21 May 2002,
09:45 AM PST,
10,
Intercomparison of simulation models for CO2 disposal in
underground storage reservoirs.
Test Problem 3: Radial Flow from a CO2 Injection Well
This problem addresses two-phase flow of CO2 and water
for simplified flow geometry and medium properties.  The 
aquifer into which injection is made is assumed infinite-acting, 
homogenoeus, and isotropic. Gravity and inertial effects are 
neglected, injection is made at a constant mass rate, and flow 
is assumed 1-D radial (line source).  Under the conditions 
stated the problem has a similarity solution where dependence on
radial distance R and time t occurs only through the similarity 
variable x = R2/t (O'Sullivan 1981; Doughty and Pruess 1992).

~Solution Control Card
Normal,
STOMP-CO2,
1,
0,day,1.e+5,day,1.e-3,s,1.e+4,day,1.15,16,1.e-06,
10000,
Variable Aqueous Diffusion,
Variable Gas Diffusion,
0,

~Grid Card
Cylindrical,
100,1,1,
0.3,m,0.34068267,m,0.386882272,m,0.439346951,m,0.498926308,m,0.566585156,m,
0.643419145,m,0.730672508,m,0.829758203,m,0.9422808,m,1.070062462,m,
1.215172455,m,1.379960655,m,1.567095601,m,1.779607712,m,2.020938356,m,
2.294995583,m,2.606217409,m,2.959643684,m,3.360997708,m,3.81677891,m,
4.334368098,m,4.922146988,m,5.589633925,m,6.347638032,m,7.208434242,m,
8.185962079,m,9.296051391,m,10.55667869,m,11.98825828,m,13.61397279,m,
15.46014866,m,17.55668242,m,19.9375248,m,22.64123061,m,25.71158298,m,
29.19830246,m,33.15785213,m,37.65435198,m,42.76061723,m,48.55933748,m,
55.14441582,m,62.62248937,m,71.11465626,m,80.75843656,m,91.70999929,m,
104.1466914,m,118.2699096,m,134.308362,m,152.5217712,m,173.2050808,m,
196.6932312,m,223.3665839,m,253.6570806,m,288.0552382,m,327.1180922,m,
371.4782167,m,421.853969,m,479.0611216,m,544.0260733,m,617.8008506,m,
701.5801442,m,796.7206557,m,904.7630673,m,1027.456991,m,1166.789304,m,
1325.016317,m,1504.700323,m,1708.751078,m,1940.472932,m,2203.618331,m,
2502.448588,m,2841.802888,m,3227.176651,m,3664.810527,m,4161.79145,m,
4726.16741,m,5367.077773,m,6094.901285,m,6921.424143,m,7860.030856,m,
8925.920993,m,10136.35532,m,11510.93531,m,13071.92059,m,14844.58935,m,
16857.64778,m,19143.69485,m,21739.75025,m,24687.85387,m,28035.74657,m,
31837.64332,m,36155.1111,m,41058.06594,m,46625.90509,m,52948.79278,m,
60129.12031,m,68283.16417,m,77542.96894,m,88058.48564,m,100000,m,
0.0,deg,45.0,deg,
0.0,m,100.0,m,

~Rock/Soil Zonation Card
1,
Aquifer,1,100,1,1,1,1,

~Mechanical Properties Card
Aquifer,2650,kg/m^3,0.12,0.12,Compressibility,4.5e-10,1/Pa,100.0,bar,Millington and Quirk,

~Hydraulic Properties Card
Aquifer,1.e-13,m^2,,,,,0.8,0.8,

~Saturation Function Card
Aquifer,van Genuchten,0.5,1/m,1.84162,0.0,0.457,0.0,

~Aqueous Relative Permeability Card
Aquifer,Mualem Irreducible,0.457,0.30,

~Gas Relative Permeability Card
Aquifer,Corey,0.30,0.05,

~Salt Transport Card
Aquifer,0.0,m,0.0,m,

~Initial Conditions Card
Gas Pressure,Aqueous Pressure,
3,
Gas Pressure,120.0,Bar,,,,,,,1,100,1,1,1,1,
Aqueous Pressure,120.0,Bar,,,,,,,1,100,1,1,1,1,
Temperature,45.0,C,,,,,,,1,100,1,1,1,1,

~Source Card
1,
Gas Mass Rate,Water-Vapor Mass Fraction,1,1,1,1,1,1,1,
0,s,120.0,bar,12.5,kg/s,0.0,

~Boundary Conditions Card
1,
East,Aqu. Dirichlet,Gas Dirichlet,Aqu. Mass Frac.,
100,100,1,1,1,1,1,
0,s,120.0,bar,0.0,120.0,bar,1.0,0.0,,

~Output Options Card
2,
33,1,1,
63,1,1,
1,1,s,m,deg,6,6,6,
12,
Aqueous Pressure,Pa,
Gas Pressure,Pa,
Aqueous Saturation,,
Gas Saturation,,
Aqueous Density,kg/m^3,
Gas Density,kg/m^3,
Aqueous Viscosity,Pa s,
Gas Viscosity,Pa s,
Aqueous Relative Permeability,,
Gas Relative Permeability,,
CO2 Aqueous Mass Fraction,,
Water Gas Mass Fraction,,
4,
30,day,
100,day,
1000,day,
10000,day,
13,
Aqueous Pressure,Pa,
Gas Pressure,Pa,
Aqueous Saturation,,
Gas Saturation,,
Aqueous Density,kg/m^3,
Gas Density,kg/m^3,
Aqueous Viscosity,Pa s,
Gas Viscosity,Pa s,
Aqueous Relative Permeability,,
Gas Relative Permeability,,
CO2 Aqueous Mass Fraction,,
Water Gas Mass Fraction,,
X Node Centroid,m,