OHara model - minor comments by Michael Clerx
Hi all,
Michael recently published some comments on the OHara model: https://docs.google.com/document/d/111fqNzQGvGAjB_PrkvejEhzqwROrR6czz_OBz7Ep-iM/edit#heading=h.q1i97y5bekil
He points out that the calculation of x1
is slightly wrong in the original publication.
openCARP's implementation has the same error:
https://git.opencarp.org/openCARP/openCARP/-/blob/master/physics/limpet/models/OHara.model?ref_type=heads#L598
He also suggests a possible optimization for Jnak, reducing the E_x computation a bit.
Additionally (depending on units), the H+ concentration might be wrong, too. https://git.opencarp.org/openCARP/openCARP/-/blob/master/physics/limpet/models/OHara.model?ref_type=heads#L580
Fortunately, these errors seem to be of little consequence – well, unfortunate for us, because we are still chasing the steady-state of the model in openCARP ;(