Graphs in Transport Problems (Zero Concentration)

[olivertheseira]

Hi,

I am a new user currently conducting geochemical modeling of CO2 sequestration in the cap rock of a certain formation. My graph of Calcite Concentration vs. Distance is showing zero concentration across the entire transport distance of 10 m. Can anyone tell me what my problem is based on my code below?

EQUILIBRIUM_PHASES 1 Carbonate assemblage
CO2(g)   -2.1   
Dolomite   0.0   32.5
Calcite   0.0   39.0   
Siderite   0.0   0.9

USER_GRAPH 1 Calcite Concentration vs. Distance from Aquifer Interface
-axis_titles      "Distance from aquifer interface (m)" "Calcite concentration (mol/L)"
-initial_solutions   true
-axis_scale x_axis
-axis_scale y_axis
-connect_simulations   true
-plot_concentration_vs   x
10 graph_x DIST
20 graph_y TOT("Calcite")

TRANSPORT
-cells                 10
-shifts         100
-time_step         100 year
-flow_direction      diffusion_only
-boundary_conditions   constant closed   
-lengths         10*1.0   
-dispersivities      0
-correct_disp      false
-diffusion_coefficient   3.33e-10
-stagnant         0
-initial_time      0
-print_cells      1-10
-print_frequency      100
-warnings         true
-multi_D         false
-interlayer_D      false

KINETICS 1
Calcite
-formula
-m      39
-m0      39
-parms   1.67e5   0.6   # cm^2/mol calcite, exp factor
-tol   1e-8

END

Thanks in advance,
Oliver Theseira

[dlparkhurst]

Check the definition of the TOT and EQUI functions in The Basic Interpreter section of the manual.