pyrtid.forward.solve_geochem_implicit

pyrtid.forward.solve_geochem_implicit(grid: RectilinearGrid, tr_model: TransportModel, gch_params: GeochemicalParameters, time_params: TimeParameters, time_index: int) → None

Compute the mineral dissolution.

The equation reads:

\[\overline{c}_{i}^{n+1} = \overline{c}_{i}^{n} + \Delta t^{n} k_{v} A_{s} \overline{c}_{i}^{n} \left( 1 - \dfrac{c_{i}^{n+1}}{K_{s}}\right)\]