Location: Smith Chase 2004 @ 45733c42c7a9 / smith_chase_nokes_shaw_wake_2004.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2009-12-02 13:53:17+13:00
Desc:
Updated the model curation status as the model does not run in COR as we were initially claiming.
Permanent Source URI:
http://models.cellml.org/workspace/smith_chase_nokes_shaw_wake_2004/rawfile/45733c42c7a99b89ae541db553e37b9588c2f0c6/smith_chase_nokes_shaw_wake_2004.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 23/11/2009 at 9:23:44 at a.m. using:

COR (0.9.31.1333)
Copyright 2002-2009 Dr Alan Garny
http://cor.physiol.ox.ac.uk/ - cor@physiol.ox.ac.uk

CellML 1.0 was used to generate this model
http://www.cellml.org/
-->
<model name="smith_chase_nokes_shaw_wake_2004" cmeta:id="smith_chase_nokes_shaw_wake_2004" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Minimal haemodynamic system model including ventricular interaction and valve dynamics.</title>
				<author>
					<firstname>Geoffrey</firstname>
					<surname>Nunns</surname>
					<affiliation>
						<shortaffil>Bioengineering Institute, The University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para> This CellML model will run in PCEnv and OpenCell to recreate the published results. Note that the model will not run in COR because it contains differential algebraic equations (DAEs).  The CellML transaltion of the model is based on the original Matlab code obtained from the model authors. However, the CellML model will only simulate 14 heart beats, because if statements (piecewise equations) were used instead of a remainder operator to cycle the heart beat, but this number of beats can be increased if desired.               </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>Abstract: Characterising circulatory dysfunction and choosing a suitable treatment is often difficult and time consuming, and can result in a deterioration in patient condition, or unsuitable therapy choices. A stable minimal model of the human cardiovascular system (CVS) is developed with the ultimate specific aim of assisting medical staff for rapid, on site modelling to assist in diagnosis and treatment. Models found in the literature simulate specific areas of the CVS with limited direct usefulness to medical staff. Others model the full CVS as a closed loop system, but they were found to be very complex, difficult to solve, or unstable. This paper develops a model that uses a minimal number of governing equations with the primary goal of accurately capturing trends in the CVS dynamics in a simple, easily solved, robust model. The model is shown to have long term stability and consistency with non-specific initial conditions as a result. An "open on pressure close on flow" valve law is created to capture the effects of inertia and the resulting dynamics of blood flow through the cardiac valves. An accurate, stable solution is performed using a method that varies the number of states in the model depending on the specific phase of the cardiac cycle, better matching the real physiological conditions. Examples of results include a 9% drop in cardiac output when increasing the thoracic pressure from -4 to 0 mmHg, and an increase in blood pressure from 120/80 to 165/130 mmHg when the systemic resistance is doubled. These results show that the model adequately provides appropriate magnitudes and trends that are in agreement with existing data for a variety of physiologically verified test cases simulating human CVS function.   </para>
				<para>The complete original paper reference is cited below:</para>
				<para>Minimal haemodynamic system model including ventricular interaction and valve dynamics, Bram W Smith, J Geoffrey Chase, Roger I Nokes, Geoffrey M Shaw, Graeme Wake, 2004, <emphasis>Medical Engineering and Physics</emphasis>, 26, 131-139. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=15036180&amp;dopt=Abstract">PubMed ID: 15036180    </ulink>
				</para>
				<informalfigure float="0" id="fig_reaction_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>model diagram</title>
							</objectinfo>
							<imagedata fileref="shaw_2007.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the minimal closed loop model of the cardiovascular system.</caption>
				</informalfigure>
			</sect1>
		</article>
	</documentation>
	
	<units xmlns="http://www.cellml.org/cellml/1.0#" name="kPa">
      <unit units="pascal" prefix="kilo"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="mL">
      <unit units="liter" prefix="milli"/>
   </units>   
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_mL">
      <unit units="liter" prefix="milli" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="kPa_per_mL">
      <unit units="kPa"/>
      <unit units="per_mL"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="kPa_second_per_mL">
      <unit units="pascal" prefix="kilo"/>
      <unit units="second"/>
      <unit units="mL" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="mL_per_second">
      <unit units="mL"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="mL_per_second2">
      <unit units="mL"/>
      <unit units="second" exponent="-2"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_second2">
      <unit units="second" exponent="-2"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="kPa_second2_per_mL">
	  <unit units="pascal" prefix="kilo"/>
	  <unit units="second" exponent="2"/>
	  <unit units="mL" exponent="-1"/>
   </units>
   
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="environment">
      <variable name="time" units="second" public_interface="out"/>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="heart_parameters">
      <variable name="R_mt" units="kPa_second_per_mL" initial_value="0.0158" public_interface="out"/>
      <variable name="R_av" units="kPa_second_per_mL" initial_value="0.0180" public_interface="out"/>
      <variable name="R_tc" units="kPa_second_per_mL" initial_value="0.0237" public_interface="out"/>
      <variable name="R_pv" units="kPa_second_per_mL" initial_value="0.0055" public_interface="out"/>
      <variable name="R_pul" units="kPa_second_per_mL" initial_value="0.1552" public_interface="out"/>
      <variable name="R_sys" units="kPa_second_per_mL" initial_value="1.0889" public_interface="out"/>
      <variable name="L_tc" units="kPa_second2_per_mL" initial_value="8.0093e-5" public_interface="out"/>
      <variable name="L_pv" units="kPa_second2_per_mL" initial_value="1.4868e-4" public_interface="out"/>
      <variable name="L_mt" units="kPa_second2_per_mL" initial_value="7.6968e-5" public_interface="out"/>
      <variable name="L_av" units="kPa_second2_per_mL" initial_value="1.2189e-4" public_interface="out"/>
      <variable name="V_tot" units="mL" initial_value="5.5" public_interface="out"/>
      <variable name="P_th" units="kPa" initial_value="-4" public_interface="out"/>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="driver_function">
      <variable name="e_t" units="dimensionless" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="A" units="dimensionless" initial_value="1"/>
      <variable name="B" units="per_second2" initial_value="80"/>
      <variable name="C" units="second" initial_value="0.375"/>
      <variable name="tau" units="second"/>
      <variable name="period" units="second" initial_value="0.75"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau</ci>
            <piecewise>
               <piece>
                  <ci>time</ci>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <ci>period</ci>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <ci>period</ci>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">4</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">4</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">5</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">5</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">6</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">6</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">7</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">7</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">8</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">8</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">9</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">9</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">10</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">10</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">11</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">11</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">12</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">12</cn>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <apply>
                        <times/>
                        <ci>period</ci>
                        <cn cellml:units="dimensionless">13</cn>
                     </apply>
                  </apply>
               </piece>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>e_t</ci>
            <apply>
               <times/>
               <ci>A</ci>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <ci>B</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
                           <ci>tau</ci>
                           <ci>C</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="pericardium">
      <variable name="V_pcd" units="mL" public_interface="out"/>
      <variable name="P_pcd" units="kPa" public_interface="out"/>
      <variable name="P_peri" units="kPa" public_interface="out"/>
      <variable name="V_lv" units="mL" public_interface="in"/>
      <variable name="V_rv" units="mL" public_interface="in"/>
      <variable name="P_th" units="kPa" public_interface="in"/>
      <variable name="P_0_pcd" units="kPa" initial_value="0.5003"/>
      <variable name="V_0_pcd" units="mL" initial_value="200"/>
      <variable name="lambda_pcd" units="per_mL" initial_value="0.03"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_pcd</ci>
            <apply>
               <plus/>
               <ci>V_lv</ci>
               <ci>V_rv</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_pcd</ci>
            <apply>
               <times/>
               <ci>P_0_pcd</ci>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>lambda_pcd</ci>
                        <apply>
                           <minus/>
                           <ci>V_pcd</ci>
                           <ci>V_0_pcd</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_peri</ci>
            <apply>
               <plus/>
               <ci>P_pcd</ci>
               <ci>P_th</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="left_ventricle">
      <variable name="V_lvf" units="mL" public_interface="out" private_interface="out"/>
      <variable name="P_lvf" units="kPa" public_interface="out"/>
      <variable name="P_lv" units="kPa" public_interface="out"/>
      <variable name="V_lv" units="mL" initial_value="94.6812" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V_spt" units="mL" public_interface="in"/>
      <variable name="e_t" units="dimensionless" public_interface="in"/>
      <variable name="P_peri" units="kPa" public_interface="in"/>
      <variable name="P_es_lvf" units="kPa" private_interface="in"/>
      <variable name="P_ed_lvf" units="kPa" private_interface="in"/>
      <variable name="R_mt" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_av" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="P_pu" units="kPa" public_interface="in"/>
      <variable name="P_ao" units="kPa" public_interface="in"/>
      <variable name="E_es_lvf" units="kPa_per_mL" public_interface="out" private_interface="in"/>
      <variable name="lambda_lvf" units="per_mL" public_interface="out" private_interface="in"/>
      <variable name="P_0_lvf" units="kPa" public_interface="out" private_interface="in"/>
      <variable name="Q_mt" units="mL_per_second" public_interface="in"/>
      <variable name="Q_av" units="mL_per_second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_lvf</ci>
            <apply>
               <minus/>
               <ci>V_lv</ci>
               <ci>V_spt</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_lvf</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>e_t</ci>
                  <ci>P_es_lvf</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>e_t</ci>
                  </apply>
                  <ci>P_ed_lvf</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_lv</ci>
            <apply>
               <plus/>
               <ci>P_lvf</ci>
               <ci>P_peri</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V_lv</ci>
            </apply>
            <piecewise>
               <piece>
                  <cn cellml:units="mL_per_second">0</cn>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <ci>Q_mt</ci>
                        <cn cellml:units="mL_per_second">0</cn>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>Q_av</ci>
                        <cn cellml:units="mL_per_second">0</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>Q_av</ci>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>Q_mt</ci>
                     <cn cellml:units="mL_per_second">0</cn>
                  </apply>
               </piece>
               <piece>
                  <ci>Q_mt</ci>
                  <apply>
                     <lt/>
                     <ci>Q_av</ci>
                     <cn cellml:units="mL_per_second">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <minus/>
                     <ci>Q_mt</ci>
                     <ci>Q_av</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="lvf_calculator">
      <variable name="P_es_lvf" units="kPa" public_interface="out"/>
      <variable name="P_ed_lvf" units="kPa" public_interface="out"/>
      <variable name="E_es_lvf" units="kPa_per_mL" initial_value="2.8798" public_interface="out"/>
      <variable name="V_lvf" units="mL" public_interface="in"/>
      <variable name="V_d_lvf" units="mL" initial_value="0"/>
      <variable name="P_0_lvf" units="kPa" initial_value="0.1203" public_interface="out"/>
      <variable name="lambda_lvf" units="per_mL" initial_value="0.033" public_interface="out"/>
      <variable name="V_0_lvf" units="mL" initial_value="0"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>P_es_lvf</ci>
            <apply>
               <times/>
               <ci>E_es_lvf</ci>
               <apply>
                  <minus/>
                  <ci>V_lvf</ci>
                  <ci>V_d_lvf</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_ed_lvf</ci>
            <apply>
               <times/>
               <ci>P_0_lvf</ci>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>lambda_lvf</ci>
                        <apply>
                           <minus/>
                           <ci>V_lvf</ci>
                           <ci>V_0_lvf</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="right_ventricle">
      <variable name="V_rvf" units="mL" public_interface="out" private_interface="out"/>
      <variable name="P_rvf" units="kPa" public_interface="out"/>
      <variable name="P_rv" units="kPa" public_interface="out"/>
      <variable name="V_rv" units="mL" initial_value="90.7302" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V_spt" units="mL" public_interface="in"/>
      <variable name="e_t" units="dimensionless" public_interface="in"/>
      <variable name="P_peri" units="kPa" public_interface="in"/>
      <variable name="P_es_rvf" units="kPa" private_interface="in"/>
      <variable name="P_ed_rvf" units="kPa" private_interface="in"/>
      <variable name="R_tc" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_pv" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="P_pa" units="kPa" public_interface="in"/>
      <variable name="P_vc" units="kPa" public_interface="in"/>
      <variable name="E_es_rvf" units="kPa_per_mL" public_interface="out" private_interface="in"/>
      <variable name="lambda_rvf" units="per_mL" public_interface="out" private_interface="in"/>
      <variable name="P_0_rvf" units="kPa" public_interface="out" private_interface="in"/>
      <variable name="Q_tc" units="mL_per_second" public_interface="in"/>
      <variable name="Q_pv" units="mL_per_second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_rvf</ci>
            <apply>
               <plus/>
               <ci>V_rv</ci>
               <ci>V_spt</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_rvf</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>e_t</ci>
                  <ci>P_es_rvf</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>e_t</ci>
                  </apply>
                  <ci>P_ed_rvf</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_rv</ci>
            <apply>
               <plus/>
               <ci>P_rvf</ci>
               <ci>P_peri</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V_rv</ci>
            </apply>
            <piecewise>
               <piece>
                  <cn cellml:units="mL_per_second">0</cn>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <ci>Q_tc</ci>
                        <cn cellml:units="mL_per_second">0</cn>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>Q_pv</ci>
                        <cn cellml:units="mL_per_second">0</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>Q_pv</ci>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>Q_tc</ci>
                     <cn cellml:units="mL_per_second">0</cn>
                  </apply>
               </piece>
               <piece>
                  <ci>Q_tc</ci>
                  <apply>
                     <lt/>
                     <ci>Q_pv</ci>
                     <cn cellml:units="mL_per_second">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <minus/>
                     <ci>Q_tc</ci>
                     <ci>Q_pv</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="rvf_calculator">
      <variable name="P_es_rvf" units="kPa" public_interface="out"/>
      <variable name="P_ed_rvf" units="kPa" public_interface="out"/>
      <variable name="E_es_rvf" units="kPa_per_mL" initial_value="0.585" public_interface="out"/>
      <variable name="V_rvf" units="mL" public_interface="in"/>
      <variable name="V_d_rvf" units="mL" initial_value="0"/>
      <variable name="P_0_rvf" units="kPa" initial_value="0.2157" public_interface="out"/>
      <variable name="lambda_rvf" units="per_mL" initial_value="0.023" public_interface="out"/>
      <variable name="V_0_rvf" units="mL" initial_value="0"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>P_es_rvf</ci>
            <apply>
               <times/>
               <ci>E_es_rvf</ci>
               <apply>
                  <minus/>
                  <ci>V_rvf</ci>
                  <ci>V_d_rvf</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_ed_rvf</ci>
            <apply>
               <times/>
               <ci>P_0_rvf</ci>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>lambda_rvf</ci>
                        <apply>
                           <minus/>
                           <ci>V_rvf</ci>
                           <ci>V_0_rvf</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="septum">
      <variable name="V_spt" units="mL" public_interface="out"/>
      <variable name="P_lv" units="kPa" public_interface="in"/>
      <variable name="P_rv" units="kPa" public_interface="in"/>
      <variable name="e_t" units="dimensionless" public_interface="in"/>
      <variable name="E_es_spt" units="kPa_per_mL" initial_value="48.754"/>
      <variable name="V_d_spt" units="mL" initial_value="2"/>
      <variable name="P_0_spt" units="kPa" initial_value="1.1101"/>
      <variable name="lambda_spt" units="per_mL" initial_value="0.435"/>
      <variable name="V_0_spt" units="mL" initial_value="2"/>
      <variable name="one" units="dimensionless" initial_value="1"/>
      <variable name="E_es_lvf" units="kPa_per_mL" public_interface="in"/>
      <variable name="P_0_lvf" units="kPa" public_interface="in"/>
      <variable name="lambda_lvf" units="per_mL" public_interface="in"/>
      <variable name="V_lv" units="mL" public_interface="in"/>
      <variable name="E_es_rvf" units="kPa_per_mL" public_interface="in"/>
      <variable name="P_0_rvf" units="kPa" public_interface="in"/>
      <variable name="V_rv" units="mL" public_interface="in"/>
      <variable name="lambda_rvf" units="per_mL" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <cn cellml:units="kPa">0</cn>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>e_t</ci>
                           <ci>E_es_spt</ci>
                           <apply>
                              <minus/>
                              <ci>V_spt</ci>
                              <ci>V_d_spt</ci>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>one</ci>
                              <ci>e_t</ci>
                           </apply>
                           <ci>P_0_spt</ci>
                           <apply>
                              <minus/>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <times/>
                                    <ci>lambda_spt</ci>
                                    <apply>
                                       <minus/>
                                       <ci>V_spt</ci>
                                       <ci>V_0_spt</ci>
                                    </apply>
                                 </apply>
                              </apply>
                              <ci>one</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>e_t</ci>
                        <ci>E_es_lvf</ci>
                        <apply>
                           <minus/>
                           <ci>V_lv</ci>
                           <ci>V_spt</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>e_t</ci>
                     </apply>
                     <ci>P_0_lvf</ci>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>lambda_lvf</ci>
                              <apply>
                                 <minus/>
                                 <ci>V_lv</ci>
                                 <ci>V_spt</ci>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>e_t</ci>
                  <ci>E_es_rvf</ci>
                  <apply>
                     <plus/>
                     <ci>V_rv</ci>
                     <ci>V_spt</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>e_t</ci>
                  </apply>
                  <ci>P_0_rvf</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>lambda_rvf</ci>
                           <apply>
                              <plus/>
                              <ci>V_rv</ci>
                              <ci>V_spt</ci>
                           </apply>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="pulmonary_artery">
      <variable name="P_pa" units="kPa" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="E_es_pa" units="kPa_per_mL" initial_value="0.369"/>
      <variable name="V_pa" units="mL" initial_value="43.0123"/>
      <variable name="V_d_pa" units="mL" initial_value="0"/>
      <variable name="R_pv" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_pul" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="P_th" units="kPa" public_interface="in"/>
      <variable name="P_rv" units="kPa" public_interface="in"/>
      <variable name="P_pu" units="kPa" public_interface="in"/>
      <variable name="Q_pul" units="mL_per_second" public_interface="in"/>
      <variable name="Q_pv" units="mL_per_second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>P_pa</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>E_es_pa</ci>
                  <apply>
                     <minus/>
                     <ci>V_pa</ci>
                     <ci>V_d_pa</ci>
                  </apply>
               </apply>
               <ci>P_th</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V_pa</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <ci>Q_pul</ci>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>Q_pv</ci>
                     <cn cellml:units="mL_per_second">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <minus/>
                     <ci>Q_pv</ci>
                     <ci>Q_pul</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="pulmonary_vein">
      <variable name="P_pu" units="kPa" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="E_es_pu" units="kPa_per_mL" initial_value="0.0073"/>
      <variable name="V_pu" units="mL" initial_value="808.4579"/>
      <variable name="V_d_pu" units="mL" initial_value="0"/>
      <variable name="R_pul" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_mt" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="P_th" units="kPa" public_interface="in"/>
      <variable name="P_pa" units="kPa" public_interface="in"/>
      <variable name="P_lv" units="kPa" public_interface="in"/>
      <variable name="Q_mt" units="mL_per_second" public_interface="in"/>
      <variable name="Q_pul" units="mL_per_second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>P_pu</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>E_es_pu</ci>
                  <apply>
                     <minus/>
                     <ci>V_pu</ci>
                     <ci>V_d_pu</ci>
                  </apply>
               </apply>
               <ci>P_th</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V_pu</ci>
            </apply>
            <piecewise>
               <piece>
                  <ci>Q_pul</ci>
                  <apply>
                     <lt/>
                     <ci>Q_mt</ci>
                     <cn cellml:units="mL_per_second">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <minus/>
                     <ci>Q_pul</ci>
                     <ci>Q_mt</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="aorta">
      <variable name="P_ao" units="kPa" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="E_es_ao" units="kPa_per_mL" initial_value="0.6913"/>
      <variable name="V_ao" units="mL" initial_value="133.3381"/>
      <variable name="V_d_ao" units="mL" initial_value="0"/>
      <variable name="R_av" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_sys" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="P_lv" units="kPa" public_interface="in"/>
      <variable name="P_vc" units="kPa" public_interface="in"/>
      <variable name="Q_av" units="mL_per_second" public_interface="in"/>
      <variable name="Q_sys" units="mL_per_second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>P_ao</ci>
            <apply>
               <times/>
               <ci>E_es_ao</ci>
               <apply>
                  <minus/>
                  <ci>V_ao</ci>
                  <ci>V_d_ao</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V_ao</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <ci>Q_sys</ci>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>Q_av</ci>
                     <cn cellml:units="mL_per_second">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <minus/>
                     <ci>Q_av</ci>
                     <ci>Q_sys</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="vena_cava">
      <variable name="P_vc" units="kPa" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="E_es_vc" units="kPa_per_mL" initial_value="0.0059"/>
      <variable name="V_vc" units="mL" initial_value="329.7803"/>
      <variable name="V_d_vc" units="mL" initial_value="0"/>
      <variable name="R_sys" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_tc" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="P_ao" units="kPa" public_interface="in"/>
      <variable name="P_rv" units="kPa" public_interface="in"/>
      <variable name="Q_sys" units="mL_per_second" public_interface="in"/>
      <variable name="Q_tc" units="mL_per_second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>P_vc</ci>
            <apply>
               <times/>
               <ci>E_es_vc</ci>
               <apply>
                  <minus/>
                  <ci>V_vc</ci>
                  <ci>V_d_vc</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V_vc</ci>
            </apply>
            <piecewise>
               <piece>
                  <ci>Q_sys</ci>
                  <apply>
                     <lt/>
                     <ci>Q_tc</ci>
                     <cn cellml:units="mL_per_second">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <minus/>
                     <ci>Q_sys</ci>
                     <ci>Q_tc</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="flow">
      <variable name="Q_sys" units="mL_per_second" public_interface="out"/>
      <variable name="Q_pul" units="mL_per_second" public_interface="out"/>
      <variable name="Q_mt" units="mL_per_second" initial_value="245.5813" public_interface="out"/>
      <variable name="Q_av" units="mL_per_second" initial_value="0" public_interface="out"/>
      <variable name="Q_tc" units="mL_per_second" initial_value="190.0661" public_interface="out"/>
      <variable name="Q_pv" units="mL_per_second" initial_value="0" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="P_ao" units="kPa" public_interface="in"/>
      <variable name="P_vc" units="kPa" public_interface="in"/>
      <variable name="P_pa" units="kPa" public_interface="in"/>
      <variable name="P_pu" units="kPa" public_interface="in"/>
      <variable name="P_rv" units="kPa" public_interface="in"/>
      <variable name="P_lv" units="kPa" public_interface="in"/>
      <variable name="R_sys" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_pul" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_mt" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_av" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_tc" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="R_pv" units="kPa_second_per_mL" public_interface="in"/>
      <variable name="L_mt" units="kPa_second2_per_mL" public_interface="in"/>
      <variable name="L_av" units="kPa_second2_per_mL" public_interface="in"/>
      <variable name="L_tc" units="kPa_second2_per_mL" public_interface="in"/>
      <variable name="L_pv" units="kPa_second2_per_mL" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Q_sys</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>P_ao</ci>
                  <ci>P_vc</ci>
               </apply>
               <ci>R_sys</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Q_pul</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>P_pa</ci>
                  <ci>P_pu</ci>
               </apply>
               <ci>R_pul</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Q_mt</ci>
            </apply>
            <piecewise>
               <piece>
                  <cn cellml:units="mL_per_second2">0</cn>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <apply>
                           <minus/>
                           <ci>P_pu</ci>
                           <ci>P_lv</ci>
                        </apply>
                        <cn cellml:units="kPa">0</cn>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>Q_mt</ci>
                        <cn cellml:units="mL_per_second">0</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <ci>P_pu</ci>
                           <ci>P_lv</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Q_mt</ci>
                           <ci>R_mt</ci>
                        </apply>
                     </apply>
                     <ci>L_mt</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Q_av</ci>
            </apply>
            <piecewise>
               <piece>
                  <cn cellml:units="mL_per_second2">0</cn>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <apply>
                           <minus/>
                           <ci>P_lv</ci>
                           <ci>P_ao</ci>
                        </apply>
                        <cn cellml:units="kPa">0</cn>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>Q_av</ci>
                        <cn cellml:units="mL_per_second">0</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <ci>P_lv</ci>
                           <ci>P_ao</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Q_av</ci>
                           <ci>R_av</ci>
                        </apply>
                     </apply>
                     <ci>L_av</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Q_tc</ci>
            </apply>
            <piecewise>
               <piece>
                  <cn cellml:units="mL_per_second2">0</cn>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <apply>
                           <minus/>
                           <ci>P_vc</ci>
                           <ci>P_rv</ci>
                        </apply>
                        <cn cellml:units="kPa">0</cn>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>Q_tc</ci>
                        <cn cellml:units="mL_per_second">0</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <ci>P_vc</ci>
                           <ci>P_rv</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Q_tc</ci>
                           <ci>R_tc</ci>
                        </apply>
                     </apply>
                     <ci>L_tc</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Q_pv</ci>
            </apply>
            <piecewise>
               <piece>
                  <cn cellml:units="mL_per_second2">0</cn>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <apply>
                           <minus/>
                           <ci>P_rv</ci>
                           <ci>P_pa</ci>
                        </apply>
                        <cn cellml:units="kPa">0</cn>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>Q_pv</ci>
                        <cn cellml:units="mL_per_second">0</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <ci>P_rv</ci>
                           <ci>P_pa</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Q_pv</ci>
                           <ci>R_pv</ci>
                        </apply>
                     </apply>
                     <ci>L_pv</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <group xmlns="http://www.cellml.org/cellml/1.0#">
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="left_ventricle">
         <component_ref component="lvf_calculator"/>
      </component_ref>
      <component_ref component="right_ventricle">
         <component_ref component="rvf_calculator"/>
      </component_ref>
   </group>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="driver_function"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="left_ventricle"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="right_ventricle"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="pulmonary_artery"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="pulmonary_vein"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="aorta"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="vena_cava"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="flow"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="driver_function" component_2="left_ventricle"/>
      <map_variables variable_1="e_t" variable_2="e_t"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="driver_function" component_2="right_ventricle"/>
      <map_variables variable_1="e_t" variable_2="e_t"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="driver_function" component_2="septum"/>
      <map_variables variable_1="e_t" variable_2="e_t"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="septum" component_2="left_ventricle"/>
      <map_variables variable_1="P_lv" variable_2="P_lv"/>
      <map_variables variable_1="V_lv" variable_2="V_lv"/>
      <map_variables variable_1="V_spt" variable_2="V_spt"/>
      <map_variables variable_1="E_es_lvf" variable_2="E_es_lvf"/>
      <map_variables variable_1="lambda_lvf" variable_2="lambda_lvf"/>
      <map_variables variable_1="P_0_lvf" variable_2="P_0_lvf"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="septum" component_2="right_ventricle"/>
      <map_variables variable_1="P_rv" variable_2="P_rv"/>
      <map_variables variable_1="V_rv" variable_2="V_rv"/>
      <map_variables variable_1="V_spt" variable_2="V_spt"/>
      <map_variables variable_1="E_es_rvf" variable_2="E_es_rvf"/>
      <map_variables variable_1="lambda_rvf" variable_2="lambda_rvf"/>
      <map_variables variable_1="P_0_rvf" variable_2="P_0_rvf"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="pericardium" component_2="left_ventricle"/>
      <map_variables variable_1="V_lv" variable_2="V_lv"/>
      <map_variables variable_1="P_peri" variable_2="P_peri"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="pericardium" component_2="right_ventricle"/>
      <map_variables variable_1="V_rv" variable_2="V_rv"/>
      <map_variables variable_1="P_peri" variable_2="P_peri"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="left_ventricle" component_2="lvf_calculator"/>
      <map_variables variable_1="P_ed_lvf" variable_2="P_ed_lvf"/>
      <map_variables variable_1="P_es_lvf" variable_2="P_es_lvf"/>
      <map_variables variable_1="V_lvf" variable_2="V_lvf"/>
      <map_variables variable_1="E_es_lvf" variable_2="E_es_lvf"/>
      <map_variables variable_1="lambda_lvf" variable_2="lambda_lvf"/>
      <map_variables variable_1="P_0_lvf" variable_2="P_0_lvf"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="right_ventricle" component_2="rvf_calculator"/>
      <map_variables variable_1="P_ed_rvf" variable_2="P_ed_rvf"/>
      <map_variables variable_1="P_es_rvf" variable_2="P_es_rvf"/>
      <map_variables variable_1="V_rvf" variable_2="V_rvf"/>
      <map_variables variable_1="E_es_rvf" variable_2="E_es_rvf"/>
      <map_variables variable_1="lambda_rvf" variable_2="lambda_rvf"/>
      <map_variables variable_1="P_0_rvf" variable_2="P_0_rvf"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="heart_parameters" component_2="left_ventricle"/>
      <map_variables variable_1="R_mt" variable_2="R_mt"/>
      <map_variables variable_1="R_av" variable_2="R_av"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="heart_parameters" component_2="right_ventricle"/>
      <map_variables variable_1="R_tc" variable_2="R_tc"/>
      <map_variables variable_1="R_pv" variable_2="R_pv"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="heart_parameters" component_2="pulmonary_artery"/>
      <map_variables variable_1="R_pv" variable_2="R_pv"/>
      <map_variables variable_1="R_pul" variable_2="R_pul"/>
      <map_variables variable_1="P_th" variable_2="P_th"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="heart_parameters" component_2="pulmonary_vein"/>
      <map_variables variable_1="R_mt" variable_2="R_mt"/>
      <map_variables variable_1="R_pul" variable_2="R_pul"/>
      <map_variables variable_1="P_th" variable_2="P_th"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="heart_parameters" component_2="aorta"/>
      <map_variables variable_1="R_av" variable_2="R_av"/>
      <map_variables variable_1="R_sys" variable_2="R_sys"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="heart_parameters" component_2="vena_cava"/>
      <map_variables variable_1="R_sys" variable_2="R_sys"/>
      <map_variables variable_1="R_tc" variable_2="R_tc"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="heart_parameters" component_2="pericardium"/>
      <map_variables variable_1="P_th" variable_2="P_th"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="heart_parameters" component_2="flow"/>
      <map_variables variable_1="R_pul" variable_2="R_pul"/>
      <map_variables variable_1="R_pv" variable_2="R_pv"/>
      <map_variables variable_1="R_tc" variable_2="R_tc"/>
      <map_variables variable_1="R_mt" variable_2="R_mt"/>
      <map_variables variable_1="R_av" variable_2="R_av"/>
      <map_variables variable_1="R_sys" variable_2="R_sys"/>
      <map_variables variable_1="L_tc" variable_2="L_tc"/>
      <map_variables variable_1="L_pv" variable_2="L_pv"/>
      <map_variables variable_1="L_mt" variable_2="L_mt"/>
      <map_variables variable_1="L_av" variable_2="L_av"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="left_ventricle" component_2="pulmonary_vein"/>
      <map_variables variable_1="P_lv" variable_2="P_lv"/>
      <map_variables variable_1="P_pu" variable_2="P_pu"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="left_ventricle" component_2="aorta"/>
      <map_variables variable_1="P_lv" variable_2="P_lv"/>
      <map_variables variable_1="P_ao" variable_2="P_ao"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="right_ventricle" component_2="pulmonary_artery"/>
      <map_variables variable_1="P_rv" variable_2="P_rv"/>
      <map_variables variable_1="P_pa" variable_2="P_pa"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="right_ventricle" component_2="vena_cava"/>
      <map_variables variable_1="P_rv" variable_2="P_rv"/>
      <map_variables variable_1="P_vc" variable_2="P_vc"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="pulmonary_vein" component_2="pulmonary_artery"/>
      <map_variables variable_1="P_pu" variable_2="P_pu"/>
      <map_variables variable_1="P_pa" variable_2="P_pa"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="aorta" component_2="vena_cava"/>
      <map_variables variable_1="P_ao" variable_2="P_ao"/>
      <map_variables variable_1="P_vc" variable_2="P_vc"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="aorta" component_2="flow"/>
      <map_variables variable_1="P_ao" variable_2="P_ao"/>
      <map_variables variable_1="Q_sys" variable_2="Q_sys"/>
      <map_variables variable_1="Q_av" variable_2="Q_av"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="vena_cava" component_2="flow"/>
      <map_variables variable_1="P_vc" variable_2="P_vc"/>
      <map_variables variable_1="Q_sys" variable_2="Q_sys"/>
      <map_variables variable_1="Q_tc" variable_2="Q_tc"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="pulmonary_vein" component_2="flow"/>
      <map_variables variable_1="P_pu" variable_2="P_pu"/>
      <map_variables variable_1="Q_mt" variable_2="Q_mt"/>
      <map_variables variable_1="Q_pul" variable_2="Q_pul"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="pulmonary_artery" component_2="flow"/>
      <map_variables variable_1="P_pa" variable_2="P_pa"/>
      <map_variables variable_1="Q_pul" variable_2="Q_pul"/>
      <map_variables variable_1="Q_pv" variable_2="Q_pv"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="left_ventricle" component_2="flow"/>
      <map_variables variable_1="P_lv" variable_2="P_lv"/>
      <map_variables variable_1="Q_mt" variable_2="Q_mt"/>
      <map_variables variable_1="Q_av" variable_2="Q_av"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="right_ventricle" component_2="flow"/>
      <map_variables variable_1="P_rv" variable_2="P_rv"/>
      <map_variables variable_1="Q_tc" variable_2="Q_tc"/>
      <map_variables variable_1="Q_pv" variable_2="Q_pv"/>
   </connection>
   
   
   <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
  <foo xmlns="http://www.example.org/foo/bar/baz" /> 
 <!-- 
      The following RDF block contains metadata that applies to this document
      as a whole, as indicated by the empty about attribute on the 
      <rdf:Description> element.
  
  --> 
 <rdf:Description rdf:about="">
<!-- 
        The Model Builder Metadata.  The Dublin Core "creator" element is used  
        to indicate the person who translated the model into CellML.
      

  --> 
 <dc:creator rdf:parseType="Resource">
 <vCard:N rdf:parseType="Resource">
  <vCard:Family>Nunns</vCard:Family> 
  <vCard:Given>Geoffrey</vCard:Given> 
  <vCard:Other>Rogan</vCard:Other> 
  </vCard:N>
<vCard:EMAIL rdf:parseType="Resource">
  <rdf:value>gnunns1@jhu.edu</rdf:value> 
  <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet" /> 
  </vCard:EMAIL>
 <vCard:ORG rdf:parseType="Resource">
  <vCard:Orgname>The University of Auckland</vCard:Orgname> 
  <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit> 
  </vCard:ORG>
  </dc:creator>
 <!-- 
        The Creation Date metadata. This is the date on which the model
        was translated into CellML.
      

  --> 
 <dcterms:created rdf:parseType="Resource">
  <dcterms:W3CDTF>2009-12-01</dcterms:W3CDTF> 
  </dcterms:created>
 <!-- 
        The Modification History metadata. This lists the changes that have been
        made to the document, who made the changes, and when they were made.
      

  --> 
  </rdf:Description>
 <!-- 
      The following metadata refers to the model itself, as indicated by the
      reference to the ID "smith_chase_nokes_shaw_wake_2004", which is    
      declared on the <model> element.
   
  --> 
 <rdf:Description rdf:about="#smith_chase_nokes_shaw_wake_2004">
 <!--  A human readable name for the model. 
  --> 
  <dc:title>The Smith et al. 2004 minimal haemodynamic system model including valve mechanics and ventricular interactions</dc:title> 
 <!--  A comment regarding the model. 
  --> 
 <cmeta:comment rdf:parseType="Resource">
  <rdf:value>This is the CellML description of Smith et al.'s 2004 minimal haemodynamic system model including valve mechanics and ventricular interactions</rdf:value> 
 <!--  The creator of the comment. 
  --> 
 <dc:creator rdf:parseType="Resource">
  <vCard:FN>Geoffrey Nunns</vCard:FN> 
  </dc:creator>
  </cmeta:comment>
  <cmeta:species>Human</cmeta:species> 
  <cmeta:bio_entity>heart</cmeta:bio_entity> 
 <!--   Keyword(s) 
  --> 
 <bqs:reference rdf:parseType="Resource">
 <dc:subject rdf:parseType="Resource">
  <bqs:subject_type>keyword</bqs:subject_type> 
 <rdf:value>
 <rdf:Bag>
  <rdf:li>cardiovascular circulation</rdf:li> 
  </rdf:Bag>
  </rdf:value>
  </dc:subject>
  </bqs:reference>
 <bqs:reference rdf:parseType="Resource">
  <bqs:Pubmed_id>15036180</bqs:Pubmed_id> 
 <bqs:JournalArticle rdf:parseType="Resource">
 <dc:creator>
 <rdf:Seq>
 <rdf:li rdf:parseType="Resource">
 <bqs:Person rdf:parseType="Resource">
 <vCard:N rdf:parseType="Resource">
  <vCard:Family>Smith</vCard:Family> 
  <vCard:Given>Bram</vCard:Given> 
  <vCard:Other>W</vCard:Other> 
  </vCard:N>
  </bqs:Person>
  </rdf:li>
  <rdf:li rdf:parseType="Resource">
  <bqs:Person rdf:parseType="Resource">
  <vCard:N rdf:parseType="Resource">
  <vCard:Family>Chase</vCard:Family> 
  <vCard:Given>J</vCard:Given> 
  <vCard:Other>Geoffrey</vCard:Other> 
  </vCard:N>
  </bqs:Person>
  </rdf:li>
   <rdf:li rdf:parseType="Resource">
  <bqs:Person rdf:parseType="Resource">
  <vCard:N rdf:parseType="Resource">
  <vCard:Family>Nokes</vCard:Family> 
  <vCard:Given>Roger</vCard:Given> 
  <vCard:Other>I</vCard:Other> 
  </vCard:N>
  </bqs:Person>
  </rdf:li>
   <rdf:li rdf:parseType="Resource">
  <bqs:Person rdf:parseType="Resource">
  <vCard:N rdf:parseType="Resource">
  <vCard:Family>Shaw</vCard:Family> 
  <vCard:Given>Geoffrey</vCard:Given> 
  <vCard:Other>M</vCard:Other> 
  </vCard:N>
  </bqs:Person>
  </rdf:li>
    <rdf:li rdf:parseType="Resource">
  <bqs:Person rdf:parseType="Resource">
  <vCard:N rdf:parseType="Resource">
  <vCard:Family>Wake</vCard:Family> 
  <vCard:Given>Graeme</vCard:Given> 
  </vCard:N>
  </bqs:Person>
  </rdf:li>

 
 
  </rdf:Seq>
  </dc:creator>
  <dc:title>Minimal haemodynamic system model including ventricular interaction and valve dynamics.</dc:title> 
 <dcterms:issued rdf:parseType="Resource">
  <dcterms:W3CDTF>2004-03-26</dcterms:W3CDTF> 
  </dcterms:issued>
 <bqs:Journal rdf:parseType="Resource">
  <dc:title>Medical Engineering and Physics</dc:title> 
  </bqs:Journal>
  <bqs:volume>26</bqs:volume> 
  <bqs:first_page>131</bqs:first_page> 
  <bqs:last_page>139</bqs:last_page> 
  </bqs:JournalArticle>
  </bqs:reference>
  </rdf:Description>
  </rdf:RDF>
</model>