- Author:
- aram148 <42922407+aram148@users.noreply.github.com>
- Date:
- 2022-07-19 16:08:46+12:00
- Desc:
- Added CO2 exchange and pH bg models
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/8c0/rawfile/afeeca156d18d9435303b544ea071f7cf868595d/bg_ph_co2.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="bg_ph_co2" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<units name="per_mol">
<unit exponent="-1" units="mole"/>
</units>
<units name="per_s">
<unit exponent="-1" units="second"/>
</units>
<units name="mol_per_s">
<unit units="mole"/>
<unit exponent="-1" units="second"/>
</units>
<units name="J_per_C">
<unit units="joule"/>
<unit exponent="-1" units="coulomb"/>
</units>
<units name="per_m">
<unit exponent="-1" units="meter"/>
</units>
<component name="main">
<variable initial_value="0" name="t" units="second"/>
<!-- State variables-->
<variable initial_value="0.0012" name="q1" units="mole"/>
<variable initial_value="0.024" name="q2" units="mole"/>
<variable initial_value="3.98e-8" name="q3" units="mole"/>
<variable name="v1" units="mol_per_s"/>
<!-- Constitutive parameters-->
<variable initial_value="0.00363384" name="K_q1" units="per_mol"/>
<variable initial_value="0.0602813" name="K_q2" units="per_mol"/>
<variable initial_value="0.0602813" name="K_q3" units="per_mol"/>
<variable initial_value="0.00021905" name="K_v1" units="per_s"/>
<variable name="pH" units="dimensionless"/>
<variable name="pHp" units="dimensionless"/>
<variable name="pKHCO3" units="dimensionless"/>
<variable name="pKHCO3_p" units="dimensionless"/>
<variable initial_value="0.225" name="alpha_p" units="dimensionless"/>
<variable initial_value="0.191" name="alpha_e" units="dimensionless"/>
<variable name="CO2_p" units="mole"/>
<variable name="CO2_T" units="mole"/>
<variable initial_value="9.3" name="Hb_p" units="mole"/>
<variable initial_value="21" name="Hb_e" units="mole"/>
<variable name="f_e" units="dimensionless"/>
<variable name="f_p" units="dimensionless"/>
<variable name="HCO3_p" units="mole"/>
<!-- Membrane parameters-->
<!-- Conservation laws-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q1</ci>
</apply>
<apply>
<minus/>
<ci>v1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q2</ci>
</apply>
<ci>v1</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q3</ci>
</apply>
<ci>v1</ci>
</apply>
<!-- Constitutive relations-->
<apply>
<eq/>
<ci>v1</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<ci>K_v1</ci>
</apply>
<ci>K_q1</ci>
<ci>q1</ci>
</apply>
<apply>
<times/>
<ci>K_v1</ci>
<ci>K_q2</ci>
<ci>K_q3</ci>
<ci>q2</ci>
<ci>q3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pH</ci>
<apply>
<minus/>
<apply>
<log/>
<ci>q3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pHp</ci>
<apply>
<minus/>
<ci>pH</ci>
<apply>
<log/>
<apply>
<minus/>
<cn cellml:units="dimensionless">3.094</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.335</cn>
<ci>pHp</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pKHCO3</ci>
<apply>
<minus/>
<ci>pH</ci>
<apply>
<log/>
<apply>
<divide/>
<ci>q2</ci>
<ci>q1</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pKHCO3_p</ci>
<ci>pKHCO3</ci>
</apply>
<apply>
<eq/>
<ci>CO2_p</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>alpha_e</ci>
<ci>alpha_p</ci>
</apply>
<ci>q1</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>HCO3_p</ci>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pHp</ci>
<ci>pKHCO3_p</ci>
</apply>
</apply>
<ci>CO2_p</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>f_p</ci>
<apply>
<divide/>
<ci>Hb_p</ci>
<ci>Hb_e</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>f_e</ci>
<apply>
<minus/>
<cn cellml:units="dim">1</cn>
<ci>f_p</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>CO2_T</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<ci>q1</ci>
<ci>q2</ci>
</apply>
<ci>f_p</ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci>CO2_p</ci>
<ci>HCO3_p</ci>
</apply>
<ci>f_e</ci>
</apply>
</apply>
</apply>
</math>
</component>
</model>