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:38 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 =====
-{{documentation:course:clusters_new_parameters.mp4|Movie: Cluster Formation & cAmp Spirals}}+ 
 +<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 
  
 <code xml |h CrawlingSlug.xml |h> <code xml |h CrawlingSlug.xml |h>
Line 271: Line 457:
 </code> </code>
  
-{{documentation:course:slug_new_parameters.mp4|Movie: Crawling Slug}}+  * {{documentation:course:slug_new_parameters.mp4|Movie: Crawling Slug}}
  
documentation/course/module6.1357904294.txt.gz · Last modified: 12:38 11.01.2013 by Fabian Rost

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki