Example b213a: One BEM one grid coupled Laplace

This example directly couples a circular boundary element region to an annular bidomain grid region. A Laplace equation is solved to test the coupling against an analytic solution.


The comfile run by this example is as follows:

#Example b213a
#
# This example directly couples a circular boundary element
#  region to an annular bidomain grid region. 
#

# Set up problem regions
#  (1 grid and 1 bem)
$GRIDHT=1 
$VENTCL=2 

# Set up problem classes
#  (1 laplace and 2 bidomain)
$TMEM=1 
$BLPL=2 

# Define problem geometry
fem define para;r;torsoslice;example
fem define regi;r;torsoslice;example
fem define node;r;gridht;example  reg $GRIDHT
fem define node;r;ventcl3;example reg $VENTCL
fem define base;r;linear;example
fem define elem;r;gridht;example  reg $GRIDHT
fem define elem;r;ventcl3;example reg $VENTCL

# Generate the grid mesh
fem define grid;r;gridht;example reg $GRIDHT
fem update grid geom reg $GRIDHT
fem update grid metric reg $GRIDHT
fem group grid xi2=0 as endocardium reg $GRIDHT
fem group grid external as boundary reg $GRIDHT
fem group grid line 1 xidirn 2 as stimulus reg $GRIDHT

# Calculate the mapping from grid points to
#  boundary element nodes
fem update node grid grid_reg $GRIDHT bem_reg $VENTCL
fem update node geomet from grid reg $VENTCL

# Define the problem setup
fem define equa;r;gridht;example reg $GRIDHT class $TMEM lock
fem define equa;r;ventcl;example reg $VENTCL class $BLPL lock
fem define mate;r;ventcl;example reg $VENTCL class $BLPL
fem define mate;d reg $GRIDHT class $TMEM
fem update grid material reg $GRIDHT class $TMEM
fem define init;r;ventcl;example reg $VENTCL class $BLPL
fem define init;d reg $GRIDHT class $TMEM
fem define solv;r;gridht;example reg $GRIDHT class $TMEM
fem define solv;r;ventcl;example reg $VENTCL class $BLPL

# Solve the laplace problem
fem solve reg $VENTCL class $BLPL

# Solve the coupled bidomain problem
fem solve activ reg $GRIDHT,$VENTCL class $TMEM,$BLPL grid endocardium cptype 2

Additional testing commands:

#Testing example b213a
fem list node sol reg 2 class 2


Files used by this example are:

Name               Modified     Size

example_b213a.com 06-Jun-2001 1.7k gridht.ipelem 06-Jun-2001 2.2k gridht.ipequa 26-May-2003 1.2k gridht.ipgrid 06-Mar-2003 539 gridht.ipnode 06-Jun-2001 5.1k gridht.ipsolv 13-Apr-2007 1.1k linear.ipbase 06-Jun-2001 2.2k test_output.com 06-Jun-2001 56 torsoslice.ippara 12-Nov-2002 5.9k torsoslice.ipregi 06-Jun-2001 93 ventcl.ipequa 26-May-2003 1.2k ventcl.ipinit 06-Jun-2001 572 ventcl.ipmate 06-Jun-2001 606 ventcl.ipsolv 13-Apr-2007 1.2k ventcl3.ipelem 06-Jun-2001 8.2k ventcl3.ipnode 06-Jun-2001 9.8k

Download the entire example:

Name                                 Modified     Size

examples_b_b2_b21_b213_b213a.tar.gz 14-Apr-2007 9.6k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:02:25 20160
cm-debugSuccessSat Mar 5 00:04:23 20160
mips-irix
cmSuccessSun Aug 19 01:29:49 20074
cm-debugSuccessWed Aug 15 01:24:20 20077
cm-debug-clear-mallocSuccessSat Aug 18 01:35:03 200717
cm-debug-clear-malloc7SuccessMon Aug 20 01:32:51 200720
cm64SuccessSun Aug 19 01:30:24 20074
cm64-debugSuccessTue Aug 21 01:23:36 20077
cm64-debug-clear-mallocSuccessThu Apr 1 11:19:37 20047
rs6000-aix
cmSuccessWed Mar 4 01:07:48 20091
cm-debugSuccessMon Mar 2 01:06:27 20092
cm64SuccessWed Mar 4 01:07:48 20091
cm64-debugSuccessTue Mar 3 01:12:17 20091
x86_64-linux
cmSuccessSun Mar 6 00:01:05 20160
cm-debugSuccessSat Mar 5 00:01:37 20160

Testing status by file:


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

Input last modified: Fri Apr 13 10:39:14 2007


CMISS Help / Examples / b / b2 / b21 / b213 / b213a