Example b213b: Two BEM one GRID Laplace equation test

This example demonstrates how to couple multiple boundary element regions by direct assembly and tests this against an analytic solution to Laplace's equation.

The comfile run by this example is as follows:

#Example b213b
# This example demonstrates how to couple multiple ventriclar
#  regions by direct assembly.

# Problem regions

# Problem classes

# Define problem geometry
fem define para;r;minimal;example
fem define regi;r;torsoslice;example
fem define node;r;LV3;example reg $LVENTL
fem define node;r;RV3;example reg $RVENTL
fem define node;r;bemhrt3;example reg $BEMHRT
fem define node;r;gridhtr;example reg $GRIDHT
fem define base;r;linear;example
fem define elem;r;LV3;example reg $LVENTL
fem define elem;r;RV3;example reg $RVENTL
fem define elem;r;bemhrt3;example reg $BEMHRT
fem define elem;r;gridhtr;example reg $GRIDHT
fem define fibr;r;gridhtr;example reg $GRIDHT
fem define elem;r;gridhtr;example reg $GRIDHT fibre

# Define finite difference grid
fem define grid;r;gridhtr;example reg $GRIDHT
fem update grid geomet reg $GRIDHT
fem update grid metric reg $GRIDHT
fem group grid xi2 0 elem 11..20,27,28 as LVendo reg $GRIDHT
fem group grid xi2 1 elem 11,13,19,20 as RVendo1 reg $GRIDHT
fem group grid xi2 0 elem 1,3,9,10,21,22,25,26 as RVendo2 reg $GRIDHT
fem group grid grid RVendo1,RVendo2 as RVendo reg $GRIDHT
fem group grid grid LVendo,RVendo as endocardium reg $GRIDHT
fem group grid xi2 1 elem 1..10,21..26 as epicardium  reg $GRIDHT
fem group grid external as boundary reg $GRIDHT
fem group grid line 1845 xidirn 2 as stimulus reg $GRIDHT

# Create node to grid mapping
fem update node grid grid_reg $GRIDHT bem_reg $BEMHRT exclude node 77,78
fem update node geomet from grid reg $BEMHRT except node 77,78

# Define equations
fem define equa;r;solve8;example reg $GRIDHT class $GRDCLS lock
fem define equa;r;LV;example reg $LVENTL class $LVECLS lock
fem define equa;r;RV;example reg $RVENTL class $RVECLS lock

# Define material parameters
fem define mate;r;LV;example reg $LVENTL class $LVECLS
fem define mate;d reg $GRIDHT class $GRDCLS
fem update grid materi reg $GRIDHT class $GRDCLS
fem define mate;r;RV;example reg $RVENTL class $RVECLS

# Define initial conditions
fem define init;d reg $GRIDHT class $GRDCLS
fem define init;r;LV;example reg $LVENTL class $LVECLS
fem define init;r;RV;example reg $RVENTL class $RVECLS

# Define right ventricle normal reversals
fem define norm;r;ventcl;example

# Define solution parameters
fem define solv;r;solve8;example reg $GRIDHT class $GRDCLS
fem define solv;r;LV;example reg $LVENTL class $LVECLS
fem define solv;r;RV;example reg $RVENTL class $RVECLS

# Solve the passive torso problems
fem solve reg $LVENTL class $LVECLS
fem solve reg $RVENTL class $RVECLS

# Solve the bidomain problem for 5ms
fem solve activ reg $GRIDHT,$LVENTL,$RVENTL class $GRDCLS,$LVECLS,$RVECLS grid endocardium cptype 4

Additional testing commands:

#Testing example b213b
fem list node sol reg 2 cl 2
fem list node sol reg 3 cl 3

Files used by this example are:

Name               Modified     Size

example_b213b.com 06-Jun-2001 2.7k LV.ipequa 26-May-2003 1.2k LV.ipinit 06-Jun-2001 1.2k LV.ipmate 06-Jun-2001 606 LV.ipsolv 13-Apr-2007 1.2k LV3.ipelem 06-Jun-2001 24k LV3.ipnode 06-Jun-2001 28k RV.ipequa 26-May-2003 1.2k RV.ipinit 06-Jun-2001 1.2k RV.ipmate 06-Jun-2001 606 RV.ipsolv 13-Apr-2007 1.2k RV3.ipelem 06-Jun-2001 24k RV3.ipnode 06-Jun-2001 28k bemhrt3.ipelem 06-Jun-2001 80k bemhrt3.ipnode 06-Jun-2001 94k gridhtr.ipelem 06-Jun-2001 7.5k gridhtr.ipelfb 06-Jun-2001 4.3k gridhtr.ipfibr 06-Jun-2001 7.0k gridhtr.ipgrid 06-Mar-2003 561 gridhtr.ipnode 06-Jun-2001 14k linear.ipbase 06-Jun-2001 2.2k minimal.ippara 12-Nov-2002 5.9k solve8.ipequa 26-May-2003 1.2k solve8.ipsolv 13-Apr-2007 1.2k test_output.com 06-Jun-2001 82 torsoslice.ipregi 06-Jun-2001 93 ventcl.ipnorm 06-Jun-2001 418

Download the entire example:

Name                                 Modified     Size

examples_b_b2_b21_b213_b213b.tar.gz 14-Apr-2007 64k

Testing status by version:

StatusTestedReal time (s)
cmSuccessSun Mar 6 00:02:26 20161
cm-debugSuccessSat Mar 5 00:06:22 20162
cmSuccessSun Aug 19 01:40:45 200711
cm-debugSuccessWed Aug 15 01:36:49 200726
cm-debug-clear-mallocSuccessSat Aug 18 01:43:32 200731
cm-debug-clear-malloc7SuccessMon Aug 20 01:40:32 200732
cm64SuccessSun Aug 19 01:41:46 200711
cm64-debugSuccessTue Aug 21 01:35:05 200728
cm64-debug-clear-mallocSuccessThu Apr 1 11:19:50 200413
cmSuccessWed Mar 4 01:09:01 20092
cm-debugSuccessMon Mar 2 01:08:58 20097
cm64SuccessWed Mar 4 01:08:36 20092
cm64-debugSuccessTue Mar 3 01:14:13 20096
cmSuccessSun Mar 6 00:01:06 20161
cm-debugSuccessSat Mar 5 00:01:43 20161

Testing status by file:

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

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

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