User Tools

Site Tools


documentation:course:module6

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
documentation:course:module6 [12:06 11.01.2013] – [Morpheus models] Fabian Rostdocumentation:course:module6 [10:06 17.11.2017] (current) – [Description] Fabian Rost
Line 16: Line 16:
     * first explain CPM and let them play with that a little (e.g. cell sorting)     * first explain CPM and let them play with that a little (e.g. cell sorting)
     * explain PDE system we are going to use     * explain PDE system we are going to use
-      * I had to explain about phase plane analysis quiet a bit+      * I had to explain about phase plane analysis quite a bit
   * then on day 3 start to couple PDE & CPM   * then on day 3 start to couple PDE & CPM
     * start with coupling in one direction only     * start with coupling in one direction only
Line 29: Line 29:
   * Stan Marée (2002), From pattern formation to morphogenesis: multicellular coordination in //Dictyostelium discoideum// [[http://www-binf.bio.uu.nl/stan/Thesis/|link]]   * Stan Marée (2002), From pattern formation to morphogenesis: multicellular coordination in //Dictyostelium discoideum// [[http://www-binf.bio.uu.nl/stan/Thesis/|link]]
   * Alberto Quintero, Mirko Myllykoski, Anna Igolkina, Alexandra Rohde O’Sullivan Freltoft, Nitya Dixit, Fabian Rost (2012), Morphogenesis and Dynamics of Multicellular Systems, [[http://www.mafy.lut.fi/EcmiNL/issues.php?action=downart&ID=280|link]]   * Alberto Quintero, Mirko Myllykoski, Anna Igolkina, Alexandra Rohde O’Sullivan Freltoft, Nitya Dixit, Fabian Rost (2012), Morphogenesis and Dynamics of Multicellular Systems, [[http://www.mafy.lut.fi/EcmiNL/issues.php?action=downart&ID=280|link]]
 +
  
  
 ===== Morpheus models ===== ===== Morpheus models =====
 +
 +<code xml |h ClusterFormation.xml |h>
 +<?xml version='1.0' encoding='UTF-8'?>
 +<MorpheusModel version="1">
 +    <Description>
 +        <Title>Cluster Formation</Title>
 +    </Description>
 +    <Space>
 +        <Lattice class="square">
 +            <Size symbol="size" value="200 200 0"/>
 +            <BoundaryConditions>
 +                <Condition boundary="x" type="constant"/>
 +                <Condition boundary="y" type="constant"/>
 +                <Condition boundary="-x" type="constant"/>
 +                <Condition boundary="-y" type="constant"/>
 +            </BoundaryConditions>
 +            <NodeLength value="0.37"/>
 +        </Lattice>
 +        <SpaceSymbol symbol="space"/>
 +    </Space>
 +    <Time>
 +        <StartTime value="0"/>
 +        <StopTime symbol="end" value="5000"/>
 +        <TimeSymbol symbol="time"/>
 +    </Time>
 +    <CellTypes>
 +        <CellType class="medium" name="medium">
 +            <Property symbol="b" value="0"/>
 +            <Property symbol="max_c" value="0"/>
 +            <Property symbol="phase" value="0"/>
 +        </CellType>
 +        <CellType class="biological" name="ameba">
 +            <VolumeConstraint>
 +                <Strength value="1.0"/>
 +                <Target value="20"/>
 +            </VolumeConstraint>
 +            <Chemotaxis>
 +                <Layer symbol-ref="c"/>
 +                <Strength symbol-ref="mu"/>
 +            </Chemotaxis>
 +            <Property symbol="b" value="0"/>
 +            <Property symbol="max_c" value="0"/>
 +            <PDEReporter>
 +                <Input symbol-ref="c"/>
 +                <Maximum symbol-ref="max_c"/>
 +            </PDEReporter>
 +            <Property symbol="phase" value="1"/>
 +            <Property symbol="phaseTime" value="0"/>
 +            <Equation symbol-ref="phaseTime">
 +                <Expression>if(phase == 2, phaseTime+MCStime, 0)</Expression>
 +            </Equation>
 +            <Equation symbol-ref="phase" name="Rule_2_3">
 +                <Expression>if(phase == 1 and max_c > 0.1, 2,
 +if(phase == 2 and phaseTime>phase2duration, 3,
 +if(phase == 3 and max_c &lt; 0.05, 1, 
 +phase
 +))) </Expression>
 +            </Equation>
 +            <Property symbol="phase2duration" value="2"/>
 +            <Function symbol="mu">
 +                <Expression>if(phase == 2, 300, 0)</Expression>
 +            </Function>
 +        </CellType>
 +        <CellType class="biological" name="autoAmeba">
 +            <VolumeConstraint>
 +                <Strength value="1.0"/>
 +                <Target value="20"/>
 +            </VolumeConstraint>
 +            <!--
 +  <Disabled>
 +    <Chemotaxis>
 +        <Layer symbol-ref="c"/>
 +        <Strength symbol-ref="mu"/>
 +    </Chemotaxis>
 +  </Disabled>  -->
 +            <Property symbol="b" value="0.5"/>
 +            <Property symbol="max_c" value="0"/>
 +            <Property symbol="phase" value="0"/>
 +        </CellType>
 +    </CellTypes>
 +    <CPM>
 +        <Interaction>
 +            <Contact type1="ameba" type2="ameba" value="4"/>
 +            <Contact type1="ameba" type2="medium" value="2"/>
 +            <Contact type1="ameba" type2="autoAmeba" value="4"/>
 +            <Contact type1="autoAmeba" type2="medium" value="3"/>
 +            <Contact type1="autoAmeba" type2="autoAmeba" value="-15"/>
 +        </Interaction>
 +        <MCSDuration symbol="MCStime" value="0.1"/>
 +        <MetropolisKinetics temperature="1.0" yield="0.1" stepper="edgelist">
 +            <Neighborhood>
 +                <Order>2</Order>
 +            </Neighborhood>
 +        </MetropolisKinetics>
 +    </CPM>
 +    <PDE>
 +        <Layer symbol="c">
 +            <Diffusion rate="1"/>
 +            <Initial>
 +                <InitPDEExpression>
 +                    <Expression>0</Expression>
 +                </InitPDEExpression>
 +            </Initial>
 +            <BoundaryConditions>
 +                <Condition boundary="x" value="0"/>
 +                <Condition boundary="-x" value="0"/>
 +                <Condition boundary="y" value="0"/>
 +                <Condition boundary="-y" value="0"/>
 +            </BoundaryConditions>
 +        </Layer>
 +        <Layer symbol="r">
 +            <Diffusion rate="0"/>
 +            <Initial>
 +                <InitPDEExpression>
 +                    <Expression>0</Expression>
 +                </InitPDEExpression>
 +            </Initial>
 +            <BoundaryConditions>
 +                <Condition boundary="x" value="0"/>
 +                <Condition boundary="-x" value="0"/>
 +                <Condition boundary="y" value="0"/>
 +                <Condition boundary="-y" value="0"/>
 +            </BoundaryConditions>
 +        </Layer>
 +        <System solver="euler" time-step="0.1">
 +            <DiffEqn symbol-ref="c">
 +                <Expression>alpha*(-f -r)</Expression>
 +            </DiffEqn>
 +            <DiffEqn symbol-ref="r">
 +                <Expression>alpha*epsilon*(3.5*c-b-r)</Expression>
 +            </DiffEqn>
 +            <Function symbol="epsilon">
 +                <Expression>if(c &lt; 0.0065, 0.5,
 +if(c &lt; 0.841, 0.0589,
 +0.5))</Expression>
 +            </Function>
 +            <Function symbol="f">
 +                <Expression>if(c &lt; 0.0065, 20*c, 
 +if(c &lt; 0.841, -3*c+0.15,
 +15*(c-1)))</Expression>
 +            </Function>
 +            <Function symbol="alpha">
 +                <Expression>if(cell.id != 0, 1, 0)</Expression>
 +            </Function>
 +        </System>
 +    </PDE>
 +    <CellPopulations>
 +        <Population size="1" type="autoAmeba">
 +            <InitRectangle cells="1" type="regular">
 +                <Dimensions size="10 10 0" center="100 100 0"/>
 +            </InitRectangle>
 +        </Population>
 +        <Population size="1" type="ameba">
 +            <InitRectangle cells="500" type="random">
 +                <Dimensions size="200 200 0" origin="0 0 0"/>
 +            </InitRectangle>
 +        </Population>
 +    </CellPopulations>
 +    <Analysis>
 +        <Gnuplotter interval="10">
 +            <Terminal opacity="0.7" name="png"/>
 +            <PDE symbol-ref="c" min="-0.5" max="1.0">
 +                <ColorMap>
 +                    <Color value="-0.5" color="dark-violet"/>
 +                    <Color value="0" color="white"/>
 +                    <Color value="1.0" color="red"/>
 +                </ColorMap>
 +            </PDE>
 +            <Cells symbol-ref="phase" min="0" max="3"/>
 +            <Cells symbol-ref="cell.volume" min="15" max="25"/>
 +            <PDE symbol-ref="c" min="-0.5" max="1.0">
 +                <ColorMap>
 +                    <Color value="-0.5" color="dark-violet"/>
 +                    <Color value="0" color="white"/>
 +                    <Color value="1.0" color="red"/>
 +                </ColorMap>
 +            </PDE>
 +        </Gnuplotter>
 +    </Analysis>
 +</MorpheusModel>
 +</code>
 +
 +
 +  * {{documentation:course:clusters_new_parameters.mp4|Movie: Cluster Formation & cAmp Spirals}}
 +  * for the movie the system size was 500x500 with 4000 cells randomly placed, everything else was like in the example above
  
  
Line 270: Line 456:
  
 </code> </code>
 +
 +  * {{documentation:course:slug_new_parameters.mp4|Movie: Crawling Slug}}
  
documentation/course/module6.1357902376.txt.gz · Last modified: 12:06 11.01.2013 by Fabian Rost

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki