[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]


Your setup runs KCl into the column for 2 years, which should move the
front 1 meter. Then you introduce 0 Cl water for 1 year, which should move
the interface between 0 Cl and 50 Cl a distance of .5 m. The velocity you
are using is .05/.1 = .5 m/yr.  Dispersivity is 0.1 m.  If all of these are
correct, PHREEQC should be solving the transport correctly.  Choosing a
constant or flux boundary for the inflow boundary condition will make a

In the first transport you are printing results at approximately 0.5 m. The
second transport the punch_frequency is greater than the number of shifts.
You have not defined SELECTED_OUTPUT, so nothing is written to the
selected-output file anyway.

Here is a plot for dispersivity of 0.1 m at the end of each transport
simulation. I think it looks ok, the large dispersivity has reduced the
peak concentration to less than 40 mmol/kgw at the end of the second
transport simulation. The flux boundary condition at the lower end may be
part of the difference with CXTFIT. It may be simulating a different
condition. The peak at the end is asymmetric because the length of time for
dispersion for the leading edge is longer than for the trailing edge.

(Embedded image moved to file: pic29755.pcx)

If you decrease the dispersivity to 0.1, the peak concentration is not
(Embedded image moved to file: pic05445.pcx)

I think PHREEQC is doing the problem you are requesting. Check that all of
the parameters and the velocity is correct. Perhaps the bottom boundary
condition is not consistent with CXTFIT and you should lengthen the column
to avoid it. Discretization could also cause differences. Otherwise, I
don't see any problem. We have tested PHREEQC against analytical solutions,
and the calculation approach the analytical solution as the cell size
decreases (see manual).

Let me know if you still think that the PHREEQC calculations are incorrect.


David Parkhurst (dlpark@xxxxxxxx)
U.S. Geological Survey
Box 25046, MS 413
Denver Federal Center
Denver, CO 80225

Project web page: https://wwwbrr.cr.usgs.gov/projects/GWC_coupled
TITLE Flow in column
      -file disp1.sel
      -reset false
      -total Cl
SOLUTION 0    # 50 mmol/l KCl
        units   mmol/l
        K       50.0
        Cl      50.0
SOLUTION 1-40  # Column without Cl
        units   mmol/l
        K       0
        Cl      0
        -cells  40
        -shifts 20
        -flow_d   forward
        -timest 3.15e6
        -bcon   constant  flux
        -length 0.05
        -disp   0.1
#           -punch_cells 10
      -punch_frequency 20
SOLUTION 0  # Original solution reenters
        units   mmol/l
        K       0
        Cl      0
        -shifts 10
        -punch_frequency    10
#        -punch  1-40

Attachment: pic29755.pcx
Description: Binary data

Attachment: pic05445.pcx
Description: Binary data

Project Home Page
Complete Water Resources Division Software
USGS Home Page
Water Resources Division Home Page
NRP Home Page
Help Page
USGS Privacy Statement       

Please note that some U.S. Geological Survey (USGS) information accessed through this page may be preliminary in nature and presented prior to final review and approval by the Director of the USGS. This information is provided with the understanding that it is not guaranteed to be correct or complete and conclusions drawn from such information are the sole responsibility of the user.

Any use of trade, product, or firm names in this publication is for descriptive purposes only and does not imply endorsement by the U.S. Government. 

The URL of this page is: https://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc/mail/msg00421.html
Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $
Visitor number 4646 since Jan 22, 1998.