- Author:
- Soroush Safaei <ssaf006@aucklanduni.ac.nz>
- Date:
- 2016-11-26 15:27:38+13:00
- Desc:
- Adding the bond graph version of AE1
- Permanent Source URI:
- https://models.cellml.org/workspace/290/rawfile/449f8d5632de6c5745dd74d7b1fe8b703f5773f0/BG/AE1_BG.cellml
<?xml version='1.0'?>
<model name="AE1_BG" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
<import xlink:href="Units.cellml">
<units name="per_mol" units_ref="per_mol"/>
<units name="J_per_mol" units_ref="J_per_mol"/>
<units name="mol_per_s" units_ref="mol_per_s"/>
</import>
<import xlink:href="Parameters-AE1.cellml">
<component component_ref="param" name="param"/>
</import>
<component name="environment">
<variable name="time" public_interface="out" units="second"/>
</component>
<component name="main">
<variable name="t" public_interface="in" units="second"/>
<!-- Constants-->
<variable name="RT" public_interface="in" units="J_per_mol"/>
<!-- State variables-->
<variable initial_value="0.01" name="q_Cl_o" units="mole"/>
<variable initial_value="0.0" name="q_Cl_i" units="mole"/>
<variable initial_value="0.0" name="q_HCO3_i" units="mole"/>
<variable initial_value="0.01" name="q_HCO3_o" units="mole"/>
<variable initial_value="0.01" name="q_E_o" units="mole"/>
<variable initial_value="0.0" name="q_ECl_o" units="mole"/>
<variable initial_value="0.0" name="q_ECl_i" units="mole"/>
<variable initial_value="0.0" name="q_E_i" units="mole"/>
<variable initial_value="0.0" name="q_EHCO3_i" units="mole"/>
<variable initial_value="0.0" name="q_EHCO3_o" units="mole"/>
<variable name="nu_Re1" units="mol_per_s"/>
<variable name="nu_Re2" units="mol_per_s"/>
<variable name="nu_Re3" units="mol_per_s"/>
<variable name="nu_Re4" units="mol_per_s"/>
<variable name="nu_Re5" units="mol_per_s"/>
<variable name="nu_Re6" units="mol_per_s"/>
<variable name="mu_E_o" units="J_per_mol"/>
<variable name="mu_Cl_o" units="J_per_mol"/>
<variable name="mu_ExCl_o" units="J_per_mol"/>
<variable name="mu_ECl_o" units="J_per_mol"/>
<variable name="mu_ECl_i" units="J_per_mol"/>
<variable name="mu_ExCl_i" units="J_per_mol"/>
<variable name="mu_Cl_i" units="J_per_mol"/>
<variable name="mu_E_i" units="J_per_mol"/>
<variable name="mu_HCO3_i" units="J_per_mol"/>
<variable name="mu_ExHCO3_i" units="J_per_mol"/>
<variable name="mu_EHCO3_i" units="J_per_mol"/>
<variable name="mu_EHCO3_o" units="J_per_mol"/>
<variable name="mu_ExHCO3_o" units="J_per_mol"/>
<variable name="mu_HCO3_o" units="J_per_mol"/>
<!-- Constitutive parameters-->
<variable name="K_Cl_o" public_interface="in" units="per_mol"/>
<variable name="K_Cl_i" public_interface="in" units="per_mol"/>
<variable name="K_HCO3_i" public_interface="in" units="per_mol"/>
<variable name="K_HCO3_o" public_interface="in" units="per_mol"/>
<variable name="K_E_o" public_interface="in" units="per_mol"/>
<variable name="K_ECl_o" public_interface="in" units="per_mol"/>
<variable name="K_ECl_i" public_interface="in" units="per_mol"/>
<variable name="K_E_i" public_interface="in" units="per_mol"/>
<variable name="K_EHCO3_i" public_interface="in" units="per_mol"/>
<variable name="K_EHCO3_o" public_interface="in" units="per_mol"/>
<variable name="k_Re1" public_interface="in" units="mol_per_s"/>
<variable name="k_Re2" public_interface="in" units="mol_per_s"/>
<variable name="k_Re3" public_interface="in" units="mol_per_s"/>
<variable name="k_Re4" public_interface="in" units="mol_per_s"/>
<variable name="k_Re5" public_interface="in" units="mol_per_s"/>
<variable name="k_Re6" public_interface="in" units="mol_per_s"/>
<!-- Conservation laws-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>mu_ExCl_o</ci>
<apply>
<plus/>
<ci>mu_E_o</ci>
<ci>mu_Cl_o</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_ExCl_i</ci>
<apply>
<plus/>
<ci>mu_Cl_i</ci>
<ci>mu_E_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_ExHCO3_i</ci>
<apply>
<plus/>
<ci>mu_E_i</ci>
<ci>mu_HCO3_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_ExHCO3_o</ci>
<apply>
<plus/>
<ci>mu_HCO3_o</ci>
<ci>mu_E_o</ci>
</apply>
</apply>
<!-- Definitions-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Cl_o</ci>
</apply>
<apply>
<minus/>
<ci>nu_Re1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Cl_i</ci>
</apply>
<ci>nu_Re3</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_HCO3_i</ci>
</apply>
<apply>
<minus/>
<ci>nu_Re4</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_HCO3_o</ci>
</apply>
<ci>nu_Re6</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_E_o</ci>
</apply>
<apply>
<minus/>
<ci>nu_Re6</ci>
<ci>nu_Re1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_ECl_o</ci>
</apply>
<apply>
<minus/>
<ci>nu_Re1</ci>
<ci>nu_Re2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_ECl_i</ci>
</apply>
<apply>
<minus/>
<ci>nu_Re2</ci>
<ci>nu_Re3</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_E_i</ci>
</apply>
<apply>
<minus/>
<ci>nu_Re3</ci>
<ci>nu_Re4</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_EHCO3_i</ci>
</apply>
<apply>
<minus/>
<ci>nu_Re4</ci>
<ci>nu_Re5</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_EHCO3_o</ci>
</apply>
<apply>
<minus/>
<ci>nu_Re5</ci>
<ci>nu_Re6</ci>
</apply>
</apply>
<!-- Constitutive relations-->
<apply>
<eq/>
<ci>mu_E_o</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_E_o</ci>
<ci>q_E_o</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Cl_o</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Cl_o</ci>
<ci>q_Cl_o</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_ECl_o</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_ECl_o</ci>
<ci>q_ECl_o</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_ECl_i</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_ECl_i</ci>
<ci>q_ECl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Cl_i</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Cl_i</ci>
<ci>q_Cl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_E_i</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_E_i</ci>
<ci>q_E_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_HCO3_i</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_HCO3_i</ci>
<ci>q_HCO3_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_EHCO3_i</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_EHCO3_i</ci>
<ci>q_EHCO3_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_EHCO3_o</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_EHCO3_o</ci>
<ci>q_EHCO3_o</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_HCO3_o</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_HCO3_o</ci>
<ci>q_HCO3_o</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>nu_Re1</ci>
<apply>
<times/>
<ci>k_Re1</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_ExCl_o</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_ECl_o</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>nu_Re2</ci>
<apply>
<times/>
<ci>k_Re2</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_ECl_o</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_ECl_i</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>nu_Re3</ci>
<apply>
<times/>
<ci>k_Re3</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_ECl_i</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_ExCl_i</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>nu_Re4</ci>
<apply>
<times/>
<ci>k_Re4</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_ExHCO3_i</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_EHCO3_i</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>nu_Re5</ci>
<apply>
<times/>
<ci>k_Re5</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_EHCO3_i</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_EHCO3_o</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>nu_Re6</ci>
<apply>
<times/>
<ci>k_Re6</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_EHCO3_o</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_ExHCO3_o</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<connection>
<map_components component_1="main" component_2="param"/>
<map_variables variable_1="RT" variable_2="RT"/>
<map_variables variable_1="K_Cl_o" variable_2="K_Cl_o"/>
<map_variables variable_1="K_Cl_i" variable_2="K_Cl_i"/>
<map_variables variable_1="K_HCO3_i" variable_2="K_HCO3_i"/>
<map_variables variable_1="K_HCO3_o" variable_2="K_HCO3_o"/>
<map_variables variable_1="K_E_o" variable_2="K_E_o"/>
<map_variables variable_1="K_ECl_o" variable_2="K_ECl_o"/>
<map_variables variable_1="K_ECl_i" variable_2="K_ECl_i"/>
<map_variables variable_1="K_E_i" variable_2="K_E_i"/>
<map_variables variable_1="K_EHCO3_i" variable_2="K_EHCO3_i"/>
<map_variables variable_1="K_EHCO3_o" variable_2="K_EHCO3_o"/>
<map_variables variable_1="k_Re1" variable_2="k_Re1"/>
<map_variables variable_1="k_Re2" variable_2="k_Re2"/>
<map_variables variable_1="k_Re3" variable_2="k_Re3"/>
<map_variables variable_1="k_Re4" variable_2="k_Re4"/>
<map_variables variable_1="k_Re5" variable_2="k_Re5"/>
<map_variables variable_1="k_Re6" variable_2="k_Re6"/>
</connection>
<connection>
<map_components component_1="main" component_2="environment"/>
<map_variables variable_1="t" variable_2="time"/>
</connection>
</model>