Conceptual Models > Design of conceptual models

EC Values in Column Experiment

(1/1)

mek001:
Dear all,
I am trying to create a 1-D flow model with PHREEQC and print the EC value at the output. I divided my column into 10 cells. I want Solution 0 to have only Na and Cl elements and minerals found in the core placed in Equilibrium Phases 1-10. I assumed the sample was saturated, so I just filled the sample with water in Solution 1-10. I assumed CO2 could infiltrate Solution 0, so I also created the Gas Phase. I tried to punch the Electrical Conductivity value by PUNCH with SC, but I failed. It only shows zero values. Meanwhile, I observe the dissolution and precipitation of minerals. Therefore, considering NaCl input and mineral dissolution I expect an EC value higher than 0 at the end. I will be very happy if you can point out my mistake. I use the wateq4f database.

--- Code: ---SOLUTION 0
    temp      25
    pH        7
    pe        4
    redox     pe
    units     mmol/kgw
    density   1
    Na        1000
    Cl        1000
    -water    0.0268 # kg
EQUILIBRIUM_PHASES 0
    CO2(g)    -3.5 10
SOLUTION 1-10 What is inside at the beginning
    temp      25
    pH        7
    pe        4
    redox     pe
    units     mmol/kgw
    density   1
    -water    0.0268 # kg
EQUILIBRIUM_PHASES 1-10
    Aragonite 0 0
    Calcite   0 0.203608752
    Chlorite7A 0 0.001749392
    Gypsum    0 0
    Halite    0 0
    Halloysite 0 0.000360242
    Illite    0 0.000282529
    Kaolinite 0 0.000374961
    Magnesite 0 0.007353813
    Quartz    0 0.001764314
    Gibbsite  0 0.001205128
    Goethite  0 0.001710748

SELECTED_OUTPUT 1
    -file                 ex12.sel
    -high_precision       true
    -reset                false
    -user_punch           true
USER_PUNCH 1
    -headings EC
    -start
30  PUNCH SC
    -end
TRANSPORT
    -cells                 10
    -shifts                55
    -time_step             550.45 # seconds
    -lengths               10*0.0093
    -dispersivities        10*0.00093
    -correct_disp          true
    -punch_cells           10


END


--- End code ---
Regards

dlparkhurst:
Wateq4f.dat does not have the definitions needed to calculate specific conductance.

Only phreeqc.dat, Amm.dat, and pitzer.dat have the necessary definitions in the SOLUTION_SPECIES data block.

Navigation

[0] Message Index