- Author:
- Randall Britten <r.britten@auckland.ac.nz>
- Date:
- 2013-09-10 19:27:27+12:00
- Desc:
- Including pH in model, setting initial pH to 7.4, which I think is physiologically realistic. Other parameters still seem far from physiological though.
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/178/rawfile/e1501746c35cafc2cfc25f935e022ef65fe8a492/pHcontrol.cellml
<?xml version="1.0"?>
<model xmlns="http://www.cellml.org/cellml/1.1#" name="pHcontrol">
<component name="Environment">
<variable name="Time" public_interface="out" units="second"/>
</component>
<component name="Flux">
<variable initial_value="2" name="CO2" units="mM"/>
<variable initial_value="0" name="HCO3" units="mM"/>
<variable initial_value="initialH" name="H" units="mM"/>
<variable initial_value="0.15" name="kf" units="per_second"/>
<variable initial_value="0.5" name="kb" units="per_mM_per_second"/>
<variable name="J" units="mM_per_second"/>
<variable name="Time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J
</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>kf
</ci>
<ci>CO2
</ci>
</apply>
<apply>
<times/>
<ci>kb
</ci>
<ci>HCO3
</ci>
<ci>H
</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>Time
</ci>
</bvar>
<ci>CO2
</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>J
</ci>
</apply>
<ci>CO2source
</ci>
</apply>
<ci>CO2sink
</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>Time
</ci>
</bvar>
<ci>HCO3
</ci>
</apply>
<ci>J
</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>Time
</ci>
</bvar>
<ci>H
</ci>
</apply>
<apply>
<plus/>
<ci>J
</ci>
<ci>protonSource
</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>pH
</ci>
<apply>
<minus/>
<apply>
<log/>
<apply>
<divide/>
<ci>H
</ci>
<ci>c1
</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>initialH
</ci>
<apply>
<times/>
<apply>
<power/>
<ci>ten
</ci>
<apply>
<minus/>
<ci>initialPh
</ci>
</apply>
</apply>
<ci>c1
</ci>
</apply>
</apply>
</math>
<variable initial_value="0" name="CO2source" units="mM_per_second"/>
<variable initial_value="0" name="CO2sink" units="mM_per_second"/>
<variable initial_value="0" name="protonSource" units="mM_per_second"/>
<variable name="pH" units="dimensionless"/>
<variable initial_value="1000" name="c1" units="mM"/>
<variable name="initialH" units="mM"/>
<variable initial_value="10" name="ten" units="dimensionless"/>
<variable initial_value="7.4" name="initialPh" units="dimensionless"/>
</component>
<units name="mM">
<unit prefix="milli" units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="per_second">
<unit exponent="-1" units="second"/>
</units>
<units name="per_mM_per_second">
<unit exponent="-1" units="mM"/>
<unit exponent="-1" units="second"/>
</units>
<units name="mM_per_second">
<unit units="mM"/>
<unit exponent="-1" units="second"/>
</units>
<connection>
<map_components component_1="Environment" component_2="Flux"/>
<map_variables variable_1="Time" variable_2="Time"/>
</connection>
</model>
