Is it possible to use PhreeqcRM transport capabilities?

[omarchaabi]

Dear all,

I am Masters Student, and I am still new to Phreeqc in general and specially PhreeqcRM. I have simple Matlab flow code (Advection) that doesn't account for component concentration calculations. I am trying to include geochemical reactions through the use of PhreeqcRM. The updated concentrations of grid cells will be used to calculate a weight factor that will affect my flow equations in each time step. So is it possible that I use phreeqc transport capabilities for component concentration shifting between cells?.  Some other point that I am wondering about is that for my flow code I am using a very small time step, do I have to use the same time step in phreeqc?

Any insight would be highly appreciated.
I would also like to thank everyone contributing to this forum as it has been of a great help to me so far.

Regards,
Omar

[dlparkhurst]

PhreeqcRM is intended to be coupled with a multicomponent flow and transport code. The idea is that you solve the flow and multicomponent transport equations and then use PhreeqcRM to solve the chemical reactions.

PHREEQC has 1D reactive transport capabilities. This would be the simplest way to do 1D calculations, and you would not need PhreeqcRM or IPhreeqc. The 1D capabilities would not be useful for 2D and 3D calculations. PHAST and VS2DRT are available for 3D and 2D calculations with limitations. PHAST is constant density saturated flow, transport, and reactions in 3D, and VS2DRT accounts for simplified flow, heat, solute transport, and chemical reactions in 2D.