Example h7: Simulation of aortic pressure in a simplified mesenteric artery network

In this example we simulate the aortic pressure cycle (Fig. 1) at the proximal end of the descending aorta in a simplified mesenteric artery network. The tree was created by Peng Du of BME by digitising the visible human data as explained in example_h4. The pressure at the distal ends of the tree is maintained constant at 9 kPa (7.2 mmHg) to maintain continuous blood flow within the network. For details of the governing flow equations and numerical methods etc. see example_h2 and example_h3 .

Fig. 1. Aortic Pressure variation

Fig. 2. Pressure Distribution at 1.09 s (Peak ejection) Fig. 3. Velocity Distribution at 1.09 s (Peak ejection)


The comfile run by this example is as follows:

# Author: Thusitha Mabotuwana
# Last modified: 15.06.06
# com file used to simulate blood flow in the mesenteric artery network.

fem def para;r;bifur;example
fem def coor;r;bifur;example
fem def base;r;bifur;example

fem def node;r;nodeFile;example
fem def elem;r;elementFile;example

fem def field;r;radiiValFile;example
fem def elem;r;radiiEleFile;example field

fem def grid;r;test;example coronary
fem update grid geometry;

fem def equat;r;bifur;example                         
fem defi mate;r;test;example
fem update material field material 4 no_field 1;

fem def time;r;periodic;example
fem def init;r;bifur;example
fem def solv;r;bifur;example

$time_initial=0.00;                       
$time_final=1;
$time_step = 0.00007;

$write_out_period= 0.087; 
# In order to generate Fig2 (at 1.09s), might have to make this much smaller than 0.087. 
# Also, if we want to reconstruct the pressure pulse which has a period of approx. 0.87s using our exnode value (only the 1st value however, ie. value at Node1), 
# we need to make this value sufficiently small so that we have enough values to capture all the characteristics of Fig1.

fem export points;$time_initial grid as vessel offset 0;         
fem export grid;vessel as vessel offset_node 0 elastic_tube; 

for ($time=$time_initial;$time<$time_final;$time=$time+$write_out_period)       
  {    
   fem solve restart to $write_out_period delta_t $time_step;        
   fem export points;$time grid as vessel offset 0           
  }                                                                           





Files used by this example are:

Name                 Modified     Size

example_h7.com 20-Jun-2006 1.5k bifur.ipbase 13-Jun-2006 2.0k bifur.ipcoor 13-Jun-2006 572 bifur.ipequa 13-Jun-2006 1.3k bifur.ipinit 13-Jun-2006 674 bifur.ippara 13-Jun-2006 5.9k bifur.ipsolv 13-Jun-2006 761 elementFile.ipelem 13-Jun-2006 62k nodeFile.ipnode 13-Jun-2006 77k periodic.iptime 13-Jun-2006 2.1k radiiEleFile.ipelfd 13-Jun-2006 23k radiiValFile.ipfiel 13-Jun-2006 15k test.ipgrid 13-Jun-2006 216 test.ipmate 13-Jun-2006 6.7k view.com 13-Jun-2006 1.2k

Download the entire example:

Name                  Modified     Size

examples_h_h7.tar.gz 19-Aug-2006 573k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:23:16 201646
cm-debugSuccessSat Mar 5 01:45:02 2016166
mips-irix
cmSuccessThu Aug 16 05:00:48 2007439
cm-debugSuccessWed Jul 25 05:11:28 20072023
cm-debug-clear-mallocSuccessSat Aug 18 05:12:41 20072034
cm-debug-clear-malloc7SuccessMon Aug 13 06:07:51 20072014
cm64SuccessThu Aug 16 04:53:23 2007430
cm64-debugSuccessTue Aug 7 06:30:13 20072055
rs6000-aix
cmSuccessWed Mar 4 01:24:58 200948
cm-debugSuccessMon Mar 2 02:18:39 2009522
cm64SuccessWed Mar 4 01:25:00 200947
cm64-debugSuccessTue Mar 3 02:11:46 2009517
x86_64-linux
cmSuccessSun Mar 6 00:36:16 20161657
cm-debugSuccessSat Mar 5 00:46:32 20162108

Testing status by file:


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

Input last modified: Tue Jun 20 11:22:41 2006


CMISS Help / Examples / h / h7