phreeqc (version 3.8.7)

ex12: 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.

References

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

See Also

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

Examples

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