- Author:
- Randall Britten <r.britten@auckland.ac.nz>
- Date:
- 2013-09-13 14:51:52+12:00
- Desc:
- Updated Mathematica to explore what equilibrium the textbook value of pK=6.1 gives compared to using the initial values read off Peter's slides.
Change pK = pKbyInitial to pK = pKfromTextbook to see the difference: textbook value does not match initial conditions on slide being at equilibrium, and ODE runs to a pH of around 4.5!
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/178/rawfile/494e3dc94373433a0bf3fdc7c37aaf20b64d095c/PressureConversion.cellml
<?xml version="1.0"?>
<model xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" cmeta:id="PressureConversion" name="PressureConversion">
<component name="ConcToPPres">
<math xmlns="http://www.w3.org/1998/Math/MathML" id="ConcToPPres">
<apply id="ConcToPPres">
<eq/>
<ci>Px
</ci>
<apply>
<divide/>
<ci>Cx
</ci>
<ci>solubility
</ci>
</apply>
</apply>
</math>
<variable name="Px" public_interface="out" units="pascal"/>
<variable name="Cx" public_interface="in" units="M"/>
<variable name="solubility" public_interface="in" units="M_per_pascal"/>
</component>
<units name="M">
<unit units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<component name="PPresToConc">
<math xmlns="http://www.w3.org/1998/Math/MathML" id="PPresToConc">
<apply id="PPresToConc">
<eq/>
<ci>Cx
</ci>
<apply>
<times/>
<ci>Px
</ci>
<ci>solubility
</ci>
</apply>
</apply>
</math>
<variable name="Cx" public_interface="out" units="M"/>
<variable name="Px" public_interface="in" units="pascal"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="Auxilliary">
<apply id="Px_proxy">
<eq/>
<ci>Px_proxy
</ci>
<ci>Px
</ci>
</apply>
</math>
<variable name="Px_proxy" units="pascal"/>
<variable name="solubility" public_interface="in" units="M_per_pascal"/>
</component>
<units name="M_per_pascal">
<unit units="M"/>
<unit exponent="-1" units="pascal"/>
</units>
</model>