You are here: Home Repository Vleaf reaction-diffusion models Fick's laws of diffusion model

Fick's laws of diffusion model

by Michael A. Guravage last modified Sep 13, 2011 04:31 PM


Fick's laws of diffusion.

Detail Description:

To implements Fick's law of chemical diffusion, insert the following code in to the myDiffusionmodel:CelltoCellTransoprt class method:

void Mydiffusionmodel::CelltoCellTransport(Wall *w, double *dchem_c1, double *dchem_c2) {
  // Passive fluxes (Fick's law)
  for (int c=0;cC1()->BoundaryPolP() || w->C2()->BoundaryPolP()) return;
    double phi = w->Length() * ( par->D[c] ) * ( w->C2()->Chemical(c) - w->C1()->Chemical(c) );
    dchem_c1[c] += phi;
    dchem_c2[c] -= phi;


Here 'w' indicates a cell wal separating two cells, w->C1() and w->C2(). The cell walls' length is given by w->Length().  

To see the details describing how to implement Fick's laws of chemical diffusion,  please refer to the section entitled 'Reaction-diffusion and cell differentiation' in 'Building simulation models of developing plant organs using Virtualleaf.'

The protocols are included with the VirtualLeaf framework. Each protocol occupies a unique revision in the mercurial source code repository, and is marked with a  mnemonic tag. Use mercurial's update command to move to the appropriate revision:

> cd virtualleaf
> hg update -r 'pin1 recycling'


The protocols are not included in the framework build process. To build and install our 'divide over arbitrary axis' model type:

> cd src/protocols
> make





Language: source code

Source: : rd fick laws

There are currently no items in this folder.

Document Actions
« February 2020 »