Location: 12 L Platform 1 model codes @ 9ac5ec0fe3f8 / Airway_FTU / Airway_FTU / branched_algeb1.cellml
- Author:
- aram148 <a.rampadarath@auckland.ac.nz>
- Date:
- 2021-12-17 11:46:34+13:00
- Desc:
- Added the Matlab 7 airway code
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/6b0/rawfile/9ac5ec0fe3f8ae60e6020b99146573defeb5a5e9/Airway_FTU/Airway_FTU/branched_algeb1.cellml
<?xml version='1.0'?>
<model name="branched_algeb" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<units name="mm">
<unit prefix="milli" units="metre"/>
</units>
<units name="mm3">
<unit exponent="3" units="mm"/>
</units>
<component name="algeb">
<variable name="Da1" public_interface="out" units="mm3"/>
<variable name="Da21" public_interface="out" units="mm3"/>
<variable name="Da22" public_interface="out" units="mm3"/>
<variable initial_value="10.19716213" name="aa" units="dimensionless"/>
<variable initial_value="1.9008e-8" name="moo1" units="dimensionless"/>
<variable name="r1" public_interface="in" units="mm"/>
<variable name="r21" public_interface="in" units="mm"/>
<variable name="r22" public_interface="in" units="mm"/>
<variable initial_value="1.8778" name="Lr1" units="mm"/>
<variable initial_value="1.7" name="Lr21" units="mm"/>
<variable initial_value="1.7" name="Lr22" units="mm"/>
<variable initial_value="-5000" name="qhat" units="dimensionless"/>
<variable name="p" public_interface="out" units="pascal"/>
<variable name="pbot" public_interface="out" units="pascal"/>
<variable name="Pmin_r1" public_interface="in" units="pascal"/>
<variable name="lambda" public_interface="out" units="mm3"/>
<variable name="temp" public_interface="out" units="mm3"/>
<variable name="ar1" units="mm"/>
<variable name="ar21" units="mm"/>
<variable name="ar22" units="mm"/>
<variable name="W" public_interface="out" units="mm3"/>
<variable name="Lambda" public_interface="out" units="mm3"/>
<variable name="q1p" public_interface="out" units="dimensionless"/>
<variable name="q21" public_interface="out" units="dimensionless"/>
<variable name="q22" public_interface="out" units="dimensionless"/>
<variable name="del_pr1" units="pascal"/>
<variable name="del_pr21" units="pascal"/>
<variable name="del_pr22" units="pascal"/>
<variable name="mu_r21" units="pascal"/>
<variable name="mu_r22" units="pascal"/>
<variable name="mu_r1" units="pascal"/>
<variable initial_value="5E1" name="A" units="dimensionless"/>
<variable name="ptop" units="pascal"/>
<variable name="tau_r1" public_interface="out" units="dimensionless"/>
<variable name="tau_r21" public_interface="out" units="dimensionless"/>
<variable name="tau_r22" public_interface="out" units="dimensionless"/>
<variable initial_value="0.445" name="R_ref" units="mm"/>
<variable name="pmidr1" public_interface="out" units="pascal"/>
<variable name="pmidr21" public_interface="out" units="pascal"/>
<variable name="pmidr22" public_interface="out" units="pascal"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>ptop</ci>
<ci>Pmin_r1</ci>
</apply>
<apply>
<eq/>
<ci>ar1</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">8</cn>
<ci>aa</ci>
<ci>moo1</ci>
<ci>Lr1</ci>
</apply>
<pi/>
</apply>
</apply>
<apply>
<eq/>
<ci>ar21</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">8</cn>
<ci>aa</ci>
<ci>moo1</ci>
<ci>Lr21</ci>
</apply>
<pi/>
</apply>
</apply>
<apply>
<eq/>
<ci>ar22</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">8</cn>
<ci>aa</ci>
<ci>moo1</ci>
<ci>Lr22</ci>
</apply>
<pi/>
</apply>
</apply>
<apply>
<eq/>
<ci>lambda</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<power/>
<ci>r21</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<ci>ar21</ci>
</apply>
<apply>
<divide/>
<apply>
<power/>
<ci>r22</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<ci>ar22</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Da1</ci>
<apply>
<divide/>
<apply>
<power/>
<ci>r1</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<ci>ar1</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Da21</ci>
<apply>
<divide/>
<apply>
<power/>
<ci>r21</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<ci>ar21</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Da22</ci>
<apply>
<divide/>
<apply>
<power/>
<ci>r22</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<ci>ar22</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>W</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>Da1</ci>
</apply>
<ci>Da21</ci>
</apply>
<ci>Da22</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>temp</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<apply>
<power/>
<ci>Da21</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<apply>
<power/>
<ci>Da22</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci>W</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<ci>lambda</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Lambda</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>temp</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>p</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>qhat</ci>
<apply>
<times/>
<ci>ptop</ci>
<ci>Da1</ci>
</apply>
</apply>
<apply>
<times/>
<ci>W</ci>
<ci>Lambda</ci>
</apply>
</apply>
</apply>
<!--*pow(W, -1{dimensionless})*pow(Lambda, -1{dimensionless});-->
<apply>
<eq/>
<ci>pbot</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>lambda</ci>
<ci>p</ci>
</apply>
<ci>qhat</ci>
</apply>
<ci>lambda</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q1p</ci>
<apply>
<minus/>
<ci>ptop</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<ci>r1</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<ci>ar1</ci>
</apply>
<ci>p</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q21</ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<ci>r21</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<ci>ar21</ci>
</apply>
<ci>p</ci>
</apply>
<ci>pbot</ci>
</apply>
</apply>
<!-- q22 = pow(r22, 4{dimensionless})/ar22*p-pbot;-->
<apply>
<eq/>
<ci>q22</ci>
<apply>
<minus/>
<ci>q1p</ci>
<ci>q21</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>del_pr1</ci>
<apply>
<minus/>
<ci>ptop</ci>
<ci>p</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>del_pr21</ci>
<apply>
<minus/>
<ci>p</ci>
<ci>pbot</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>del_pr22</ci>
<apply>
<minus/>
<ci>p</ci>
<ci>pbot</ci>
</apply>
</apply>
<!-- For parenchymal tethering-->
<apply>
<eq/>
<ci>mu_r21</ci>
<apply>
<divide/>
<apply>
<abs/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<divide/>
<ci>A</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>del_pr21</ci>
<apply>
<power/>
<ci>r21</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
<apply>
<times/>
<ci>del_pr22</ci>
<apply>
<power/>
<ci>r22</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_r22</ci>
<apply>
<divide/>
<apply>
<abs/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<divide/>
<ci>A</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>del_pr22</ci>
<apply>
<power/>
<ci>r22</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
<apply>
<times/>
<ci>del_pr21</ci>
<apply>
<power/>
<ci>r21</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_r1</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>mu_r21</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>mu_r22</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_r1</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>mu_r1</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<ci>R_ref</ci>
<ci>r1</ci>
</apply>
<ci>R_ref</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1.5</cn>
<apply>
<power/>
<apply>
<divide/>
<apply>
<minus/>
<ci>R_ref</ci>
<ci>r1</ci>
</apply>
<ci>R_ref</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_r21</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>mu_r21</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<ci>R_ref</ci>
<ci>r21</ci>
</apply>
<ci>R_ref</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1.5</cn>
<apply>
<power/>
<apply>
<divide/>
<apply>
<minus/>
<ci>R_ref</ci>
<ci>r21</ci>
</apply>
<ci>R_ref</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_r22</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>mu_r22</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<ci>R_ref</ci>
<ci>r22</ci>
</apply>
<ci>R_ref</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1.5</cn>
<apply>
<power/>
<apply>
<divide/>
<apply>
<minus/>
<ci>R_ref</ci>
<ci>r22</ci>
</apply>
<ci>R_ref</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pmidr1</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<plus/>
<ci>ptop</ci>
<ci>p</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pmidr21</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<plus/>
<ci>pbot</ci>
<ci>p</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pmidr22</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<plus/>
<ci>pbot</ci>
<ci>p</ci>
</apply>
</apply>
</apply>
</math>
</component>
</model>