Example b2146: One ventricle and one torso region directly coupled

This example demonstrates how to couple one ventricular region and one torso cavity region by direct assembly.


The comfile run by this example is as follows:

#Example b2146
#
# This example demonstrates how to couple one ventricular
#  region and one torso cavity region by direct assembly.
#

# Problem classes
$CAVITY=1 
$VENTCL=2 
$BEMHRT=3 
$GRIDHT=4 

# Problem regions
$TLPL=1 
$BLPL=2 
$TMEM=3 
$EXTC=4 

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

# Define finite difference grid
fem define grid;r;gridht;example reg $GRIDHT
fem update grid geom reg $GRIDHT
fem update grid metric reg $GRIDHT
fem group grid external as boundary reg $GRIDHT
fem group grid line 1 xidirn 2 as stimulus reg $GRIDHT

# Create node to grid mapping
fem update node grid grid_reg $GRIDHT bem_reg $BEMHRT
fem update node geomet from grid reg $BEMHRT

# Define equations
fem define equa;r;cavity;example reg $CAVITY class $TLPL lock
fem define equa;r;ventcl;example reg $VENTCL class $BLPL lock
fem define equa;r;gridht;example reg $GRIDHT class $TMEM,$EXTC lock

# Define material parameters
fem define mate;r;cavity;example reg $CAVITY class $TLPL
fem define mate;r;ventcl;example reg $VENTCL class $BLPL
fem define mate;r;gridht;example reg $GRIDHT class $TMEM,$EXTC
fem define cell;r;gridht;example reg $GRIDHT class $TMEM,$EXTC
fem update grid material reg $GRIDHT class $TMEM,$EXTC

# Define initial conditions
fem define init;r;cavity;example reg $CAVITY class $TLPL
fem define init;r;ventcl;example reg $VENTCL class $BLPL
fem define init;r;gridht;example reg $GRIDHT class $TMEM,$EXTC

# Define solution parameters
fem define solv;r;gridht;example reg $GRIDHT class $TMEM
fem define solv;r;gridht;example reg $GRIDHT class $EXTC
fem define solv;r;cavity;example reg $CAVITY class $TLPL
fem define solv;r;ventcl;example reg $VENTCL class $BLPL

# Solve the passive torso problems
fem solve reg $CAVITY class $TLPL
fem solve reg $VENTCL class $BLPL

# Solve the bidomain problem for 10ms
fem solve activ coupled to 10 reg $GRIDHT,$VENTCL,$CAVITY class $TMEM,$EXTC,$BLPL,$TLPL grid boundary cptype 6 fixed_node 17

Additional testing commands:

#Testing example b2146
fem list node sol reg 1 cl 1


Files used by this example are:

Name               Modified     Size

example_b2146.com 06-Jun-2001 2.3k bemhrt3.ipelem 06-Jun-2001 16k bemhrt3.ipnode 06-Jun-2001 19k cavity.ipequa 26-May-2003 1.2k cavity.ipinit 06-Jun-2001 695 cavity.ipmate 06-Jun-2001 606 cavity.ipsolv 13-Apr-2007 1.2k cavity3.ipelem 06-Jun-2001 16k cavity3.ipnode 06-Jun-2001 19k gridht.ipcell 06-Jun-2001 1.7k gridht.ipelem 06-Jun-2001 2.2k gridht.ipequa 26-May-2003 1.4k gridht.ipgrid 06-Mar-2003 539 gridht.ipinit 06-Jun-2001 841 gridht.ipmate 06-Jun-2001 1.8k gridht.ipnode 06-Jun-2001 5.1k gridht.ipsolv 13-Apr-2007 1.7k gridht.ipsolv.old 06-Jun-2001 1.4k linear.ipbase 06-Jun-2001 2.2k test_output.com 06-Jun-2001 53 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_b214_b2146.tar.gz 14-Apr-2007 65k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:02:28 20161
cm-debugSuccessSat Mar 5 00:08:20 20162
mips-irix
cmSuccessSun Aug 19 01:41:55 200711
cm-debugSuccessWed Aug 15 01:33:32 200727
cm-debug-clear-mallocSuccessSat Aug 18 01:44:39 200732
cm-debug-clear-malloc7SuccessMon Aug 20 01:41:08 200732
cm64SuccessSun Aug 19 01:42:35 200711
cm64-debugSuccessTue Aug 21 01:35:49 200727
cm64-debug-clear-mallocSuccessThu Apr 1 11:22:02 200414
rs6000-aix
cmSuccessWed Mar 4 01:09:17 20092
cm-debugSuccessMon Mar 2 01:08:34 20095
cm64SuccessWed Mar 4 01:09:17 20092
cm64-debugSuccessTue Mar 3 01:14:12 20096
x86_64-linux
cmSuccessSun Mar 6 00:01:07 20161
cm-debugSuccessSat Mar 5 00:01:44 20161

Testing status by file:


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

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


CMISS Help / Examples / b / b2 / b21 / b214 / b2146