feat: Separate Mobile and Immobile Parts of Aggregate Species Concentrations

[frankfeifan]

This PR aims to separate the mobile and immobile parts of each aggregate species in HPCReact. It is motivated by the fact that in some systems, the secondary species can be either mobile or immobile, so when doing reaction transport on the aggregate species concentrations, we should only transport the mobile part in the flux assembly while use the entire aggregate concentration in the accumulation assembly.

The main changes include:

• A new flag array called m_mobileSecondarySpeciesFlag in both EquilibriumReactionsParameters and MixedReactionParameters to indicate the mobility of the secondary species: 1 if mobile and 0 if immobile.
• A new function in massActions that separates the update for the total and mobile aggregate species concentration.
• Call the new function in massActions in the mixed system update.