Example b21528: Broken activation on knife collapse

This example solves an activation problem using the Beeler-Reuter model on a "knife" type element. There are two ways of ordering the collapse in this type of element, which reverse the order of xi collapsing. See b21527 to compare. The order used in this example does not work for activation and there is a floating point error due to numerical instability at time 3.7

The comfile run by this example is as follows:

# This is the example of a knife type element that gives incorrect
# result.

# Set up for bidomain.

fem define parameters;r;knife
fem define coordinates;r;knife

fem define node;r;knife

fem define base;r;knife

fem define elem;r;knife

fem define grid;r;knife
fem update grid geometry
fem update grid metric

fem ex data;normals normals as normals
fem group grid external as boundary
fem group grid xi2 high as stim_point oneoff

fem def equa;r;knife class $VM,$PHI_E
fem def mate;r;knife class $VM,$PHI_E
fem def cell;r;knife class $VM,$PHI_E

fem update grid material class $VM,$PHI_E

fem def init;r;knife class $VM,$PHI_E

fem def solve;r;knife class $VM
fem def solve;r;knife class $PHI_E
fem solve to 0 class $VM,$PHI_E

# Export the grid elements and nodes
fem export node;Domain as Domain reg 1
fem export elem;Domain grid_numbers as Domain reg 1
fem export elem;"MemPot_000" field as Solution class $VM

  for $M ( 0..10 )
   for $L ( 0..9 )
   #Solve for new time
   print "Time step: $M.$L\n"; 
   fem solve restart to "${M}.$L" class $VM,$PHI_E

   # Export solution
   fem export elem;"MemPot_${M}.$L" field as Solution class $VM


Files used by this example are:

Name                Modified     Size

example_b21528.com 29-Apr-2003 1.2k Domain.exelem 29-Apr-2003 74k Domain.exnode 29-Apr-2003 1.4k MemPot_0.0.exelem 29-Apr-2003 95k MemPot_0.1.exelem 29-Apr-2003 95k MemPot_0.2.exelem 29-Apr-2003 95k MemPot_0.3.exelem 29-Apr-2003 95k MemPot_0.4.exelem 29-Apr-2003 95k MemPot_0.5.exelem 29-Apr-2003 95k MemPot_0.6.exelem 29-Apr-2003 95k MemPot_0.7.exelem 29-Apr-2003 95k MemPot_0.8.exelem 29-Apr-2003 95k MemPot_0.9.exelem 29-Apr-2003 95k MemPot_000.exelem 29-Apr-2003 95k MemPot_1.0.exelem 29-Apr-2003 95k MemPot_1.1.exelem 29-Apr-2003 95k MemPot_1.2.exelem 29-Apr-2003 95k MemPot_1.3.exelem 29-Apr-2003 95k MemPot_1.4.exelem 29-Apr-2003 95k MemPot_1.5.exelem 29-Apr-2003 95k MemPot_1.6.exelem 29-Apr-2003 95k MemPot_1.7.exelem 29-Apr-2003 95k MemPot_1.8.exelem 29-Apr-2003 95k MemPot_1.9.exelem 29-Apr-2003 95k MemPot_2.0.exelem 29-Apr-2003 95k MemPot_2.1.exelem 29-Apr-2003 95k MemPot_2.2.exelem 29-Apr-2003 95k MemPot_2.3.exelem 29-Apr-2003 95k MemPot_2.4.exelem 29-Apr-2003 95k MemPot_2.5.exelem 29-Apr-2003 95k MemPot_2.6.exelem 29-Apr-2003 95k MemPot_2.7.exelem 29-Apr-2003 95k MemPot_2.8.exelem 29-Apr-2003 95k MemPot_2.9.exelem 29-Apr-2003 95k MemPot_3.0.exelem 29-Apr-2003 95k MemPot_3.1.exelem 29-Apr-2003 95k MemPot_3.2.exelem 29-Apr-2003 95k MemPot_3.3.exelem 29-Apr-2003 95k MemPot_3.4.exelem 29-Apr-2003 95k MemPot_3.5.exelem 29-Apr-2003 95k MemPot_3.6.exelem 29-Apr-2003 95k knife.ipbase 29-Apr-2003 2.4k knife.ipcell 29-Apr-2003 2.0k knife.ipcoor 29-Apr-2003 570 knife.ipelem 29-Apr-2003 721 knife.ipequa 26-May-2003 1.6k knife.ipgrid 29-Apr-2003 606 knife.ipinit 29-Apr-2003 903 knife.ipmate 29-Apr-2003 1.8k knife.ipnode 29-Apr-2003 2.4k knife.ippara 29-Apr-2003 5.9k knife.ipsolv 13-Apr-2007 1.8k normals.exdata 29-Apr-2003 190k

Download the entire example:

Name                                        Modified     Size

examples_b_b2_b21_b215_b2152_b21528.tar.gz 14-Apr-2007 401k

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

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

CMISS Help / Examples / b / b2 / b21 / b215 / b2152 / b21528