NUMERICAL METHOD FOR SPECIATION AND FORWARD MODELING
For the explicit calculation of the diffuse layer, a charge-balance equation is used for each surface, ; the values of the master unknowns for each surface of the surface assemblage, and , are adjusted to achieve mole balance and charge balance for each surface. If the diffuse-layer composition is not explicitly included in the calculation, then the charge-potential equation, , is used in place of the surface charge-balance equation. If the non-electrostatic model is used for the surface assemblage, then neither the surface charge balance nor the charge-potential equation is included in the set of equations to be solved.
All equations for initial surface calculations are included as equality constraints in the solver. No equations are optimized and no inequality constraints are included.
An initial surface calculation is performed only if the surface initially is defined to be in equilibrium with a specified solution. The distribution of species for this solution has already been calculated, either by an initial solution calculation or by a reaction or transport calculation. Thus, the initial estimates of all master unknowns related to the aqueous phase are set equal to the values from the previous distribution of species. The initial estimate of the activity of the master species for each surface is set equal to one tenth of the number of moles of surface sites for that surface. For explicit and implicit diffuse-layer calculations, the initial estimate of the potential unknown, , for each surface is zero, which implies that the surface potential is zero.
For data input to PHREEQC, definition of the initial surface calculation is made with the SURFACE keyword data block. (See Description of Data Input.)