Example b2122: 2d circles coupled bidomain fem-bem-fem

Example b2122 - coupled bidomain circles fem-bem-fem This example solves the bidomain equation through to a point where an activation wave is well established on a fem inner layer with grid points and then iterates between the extracellular and the boundary element problem to get the epicardial fluxes to match. This is similar to example b8 but using concentric circles and a fem outer layer.


The comfile run by this example is as follows:


# Command updated by fixcom.sh on Wed Aug 23 17:50:04 NZT 2000
# Old command: assign variable BEM 1
$BEM = 1


# Command updated by fixcom.sh on Wed Aug 23 17:50:04 NZT 2000
# Old command: assign variable GRID 2
$GRID = 2


# Command updated by fixcom.sh on Wed Aug 23 17:50:04 NZT 2000
# Old command: assign variable FEMLAYER 3
$FEMLAYER = 3


# Command updated by fixcom.sh on Wed Aug 23 17:50:04 NZT 2000
# Old command: assign variable LAPLACE 1
$LAPLACE = 1


# Command updated by fixcom.sh on Wed Aug 23 17:50:04 NZT 2000
# Old command: assign variable TMEMBRANE 2
$TMEMBRANE = 2


# Command updated by fixcom.sh on Wed Aug 23 17:50:04 NZT 2000
# Old command: assign variable EXTCELLR 3
$EXTCELLR = 3


# Command updated by fixcom.sh on Wed Aug 23 17:50:04 NZT 2000
# Old command: assign variable EXTUPDATE 4
$EXTUPDATE = 4


fem define para;r;fembem;example

fem define regi;r;fembem;example
fem define node;r;test;example reg all
fem define base;r;fembem;example
fem define elem;r;test;example reg all
fem define fibr;d reg $GRID
fem define fibr;d reg $FEMLAYER
fem define elem;r;grid;example fibre reg $GRID
fem define elem;r;fem;example fibre reg $FEMLAYER
fem define grid;r;fem;example reg $GRID
fem update grid geometry reg 2
fem update grid metric reg 2
fem group grid external as boundary reg 2
fem define equa;r;bem;example reg $BEM class $LAPLACE lock
fem define equa;r;bem2;example reg $FEMLAYER class $LAPLACE
fem define equa;r;fem;example reg $GRID class $TMEMBRANE,$EXTCELLR,$EXTUPDATE lock
fem define mate;r;bem;example reg $BEM class $LAPLACE
fem define mate;r;bem2;example reg $FEMLAYER class $LAPLACE
fem define mate;r;fem;example reg $GRID class $TMEMBRANE,$EXTCELLR
fem define cell;r;fem;example reg $GRID class $TMEMBRANE,$EXTCELLR
fem update grid material reg $GRID class $TMEMBRANE
fem define init;r;fem;example reg $GRID class $TMEMBRANE,$EXTCELLR
fem define init;r;bem;example reg $BEM class $LAPLACE
fem define init;r;bem2;example reg $FEMLAYER class $LAPLACE
fem define coup;r;fembem;example
fem define solv;r;fem;example reg $GRID class $TMEMBRANE
fem define solv;r;fem2;example reg $GRID class $EXTCELLR
fem define solv;r;fem3;example reg $GRID class $EXTUPDATE
fem define solve;r;bem;example coupled class $LAPLACE reg $BEM,$FEMLAYER

fem solve to 10 reg $GRID class $TMEMBRANE,$EXTCELLR,$EXTUPDATE
fem update node potential reg $BEM,$GRID
fem solve coupled class $LAPLACE reg $BEM,$FEMLAYER
fem update node potential reg $BEM,$GRID check


# Command updated by fixcom.sh on Wed Aug 23 17:50:04 NZT 2000
# Old command: DO A=1..6
for $A ( 1..6 ) 
{

fem solve iterate 1 reg $GRID,$BEM class $TMEMBRANE,$EXTCELLR,$EXTUPDATE,$LAPLACE
fem update node potential reg $BEM,$GRID
fem solve coupled reg $BEM,$FEMLAYER class $LAPLACE
fem update node potential reg $BEM,$GRID check

# Command updated by fixcom.sh on Wed Aug 23 17:50:04 NZT 2000
# Old command: ENDDO
}



fem define window;c
fem dr field grid 1 zmin=-50 zmax=9 noscale reg $GRID class $EXTCELLR
fem dr field zmin=-50 zmax=9 noscale reg $FEMLAYER class $LAPLACE

Files used by this example are:

Name                  Modified     Size

example_b2122.com 20-Nov-2001 3.0k bem.ipcell 10-Apr-2000 590 bem.ipequa 26-May-2003 1.3k bem.ipinit 10-Apr-2000 380 bem.ipmate 10-Apr-2000 590 bem.ipmate_cell_old 10-Apr-2000 590 bem.irsolv 13-Apr-2007 1.5k bem2.ipcell 10-Apr-2000 590 bem2.ipequa 26-May-2003 1.4k bem2.ipinit 10-Apr-2000 480 bem2.ipmate 10-Apr-2000 590 bem2.ipmate_cell_old 10-Apr-2000 590 fem.ipcell 10-Apr-2000 2.3k fem.ipelfb 10-Apr-2000 1.2k fem.ipequa 26-May-2003 1.5k fem.ipgrid 06-Mar-2003 1.1k fem.ipinit 10-Apr-2000 3.7k fem.ipmate 10-Apr-2000 1.7k fem.ipmate_cell_old 10-Apr-2000 4.0k fem.ipsolv 13-Apr-2007 1.8k fem.ipsolv.old 10-Apr-2000 1.5k fem2.ipsolv 13-Apr-2007 1.8k fem2.ipsolv.old 10-Apr-2000 1.5k fem3.ipsolv 13-Apr-2007 1.8k fem3.ipsolv.old 10-Apr-2000 1.5k fembem.ipbase 10-Apr-2000 2.5k fembem.ipcoup 21-Aug-2002 437 fembem.ipcoup.old 10-Apr-2000 411 fembem.ippara 12-Nov-2002 5.9k fembem.ipregi 10-Apr-2000 93 grid.ipelfb 10-Apr-2000 1.2k output.com 10-Apr-2000 158 test.irelem 10-Apr-2000 8.5k test.irnode 10-Apr-2000 15k test_old.irelem 10-Apr-2000 8.5k

Download the entire example:

Name                                 Modified     Size

examples_b_b2_b21_b212_b2122.tar.gz 14-Apr-2007 15k

Html last generated: Sun Mar 6 05:51:10 2016

Input last modified: Fri Apr 13 10:38:23 2007


CMISS Help / Examples / b / b2 / b21 / b212 / b2122