How to determine the corresponding relation between phreeqc cells and mesh nodes

[Yongqiang]

Hi everyone,

When we couple PHREEQC with a multiphase simulator, we can have n cells to represent the component at n mesh nodes of the multiphase simulator. So how to determine the mapping or projection relation between cells and nodes?

Best Regards,
Yongqiang

[dlparkhurst]

Usually, there is a natural numbering system for your flow and transport nodes, usually from 0 or 1 to n. If so, you can simply assign cells in PhreeqcRM with the same numbers.

In any case, assuming n nodes in the flow and transport simulator, PhreeqcRM will have a set of cells numbered from 0 to m, where m is less than or equal to n. You have to figure out how to assign concentrations calculated by the transport simulator at its nodes to cells of PhreeqcRM. You must develop an array with n concentrations of e elements that is submitted to to PhreeqcRM, which is mapped to the m cells of PhreeqcRM as you define with the method CreateMapping.

[Yongqiang]

Thank You, David!