ex12 {phreeqc} | R Documentation

## Example 12–Advective and Diffusive Flux of Heat and Solutes

### Description

This example demonstrates the capability of PHREEQC to calculate transient transport of heat and solutes in a column or along a 1D flowline. A column is initially filled with a dilute KCl solution at 25 degrees C in equilibrium with a cation exchanger. A KNO3 solution then advects into the column and establishes a new temperature of 0 degrees C. Subsequently, a sodium chloride solution at 24 degrees C is allowed to diffuse from both ends of the column, assuming no heat is lost through the column walls. At one end, a constant boundary condition is imposed, and at the other end, the final cell is filled with the sodium chloride solution and a closed boundary condition is prescribed. For the column end with a constant boundary condition, an analytical solution is compared with PHREEQC results, for unretarded Cl- (R = 1.0) and retarded Na+ and temperature (R = 3.0). Finally, the second-order accuracy of the numerical method is verified by increasing the number of cells by a factor of three and demonstrating a decrease in the error of the numerical solution by approximately one order of magnitude relative to the analytical solution. The example can be run using the `phrRunString` routine.

### Source

http://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc

### References

http://pubs.usgs.gov/tm/06/a43/pdf/tm6-A43.pdf

### See Also

Other Examples: `ex10`, `ex11`, `ex13a`, `ex14`, `ex15`, `ex16`, `ex17`, `ex18`, `ex19`, `ex1`, `ex20a`, `ex21`, `ex22`, `ex2`, `ex3`, `ex4`, `ex5`, `ex6`, `ex7`, `ex8`, `ex9`

### Examples

    phrLoadDatabaseString(phreeqc.dat)
    phrSetOutputStringsOn(TRUE)
    phrRunString(ex12)
    phrGetOutputStrings()

[Package phreeqc version 3.4.10 Index]