Location: 12 L Platform 1 model codes @ ad134c275417 / Airway_FTU / branched_no_Ca / branched_algeb1.cellml

Author:
aram148 <a.rampadarath@auckland.ac.nz>
Date:
2021-09-03 09:16:17+12:00
Desc:
renamed VSM
Permanent Source URI:
https://models.physiomeproject.org/workspace/6b0/rawfile/ad134c275417b15f67caabb85d67fd244645312b/Airway_FTU/branched_no_Ca/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>
    <units name="mmHg">
        <unit multiplier="133.32" units="pascal"/>
    </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/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <apply>
                                    <plus/>
                                    <ci>Da21</ci>
                                    <ci>Da22</ci>
                                </apply>
                                <ci>qhat</ci>
                            </apply>
                            <ci>lambda</ci>
                        </apply>
                        <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>
            <apply>
                <eq/>
                <ci>q22</ci>
                <apply>
                    <minus/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <apply>
                                <power/>
                                <ci>r22</ci>
                                <cn cellml:units="dimensionless">4</cn>
                            </apply>
                            <ci>ar22</ci>
                        </apply>
                        <ci>p</ci>
                    </apply>
                    <ci>pbot</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>