Example 845: Lax Wendroff finite difference flow solution

One dimensional solution to the coronary blood flow equations independent of any host mesh mechanical deformation using the Lax Wendroff Finte Differences.


The comfile run by this example is as follows:

fem define para;r;flow;example
fem reallocate
fem define coord;r;flow;example
fem define nodes;r;flow;example
fem define field;r;flow;example
fem define base;r;flow;example
fem define elements;r;flow;example
fem define elements;r;flow;example field        
fem define grid;r;flow;example coronary          
# generates the coronary grid geometry
fem update grid geometry                        
fem define equation;r;flow;example               
# flow in elastic tubes
fem list para nq                         
fem define material;r;flow;example               
# defines the material and fluid properties for the problem
fem update material field material 4 no_field 1; 
# updates unstressed radius (mat 4) from field variable 1
fem define time;r;flow;example                   
# defines boundary conditions
fem define initial;r;flow;example                 
# defines the inital conditions for fluid flow in the network
fem define solve;r;flow;example                   

$time_initial=0.00;
$time_final=0.004;
$time_step=0.0001;   
$write_out_period=0.0004;

fem export points;$time_initial grid as flow offset 0;  # exports initial soln
fem export grid;grid as flow 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+$write_out_period grid as flow offset 0           
}                                                                           



Additional testing commands:

fem list nodes

Files used by this example are:

Name             Modified     Size

example_845.com 14-Aug-2001 1.5k 0.0004.exnode 04-May-2005 5.1k 0.0008.exnode 04-May-2005 5.1k 0.0012.exnode 04-May-2005 5.1k 0.0016.exnode 04-May-2005 5.1k 0.002.exnode 04-May-2005 5.1k 0.0024.exnode 04-May-2005 5.1k 0.0028.exnode 04-May-2005 5.1k 0.0032.exnode 04-May-2005 5.1k 0.0036.exnode 04-May-2005 5.1k 0.004.exnode 04-May-2005 5.1k 0.exnode 14-Aug-2001 5.1k flow.ipbase 14-Aug-2001 893 flow.ipcoor 10-Apr-2000 572 flow.ipelem 14-Aug-2001 1.0k flow.ipelfd 14-Aug-2001 818 flow.ipequa 22-Nov-2005 1.3k flow.ipfiel 14-Aug-2001 1.3k flow.ipgrid 27-Nov-2001 224 flow.ipinit 01-Dec-2005 929 flow.ipmate 10-Apr-2000 6.8k flow.ipnode 14-Aug-2001 1.1k flow.ippara 12-Nov-2002 5.9k flow.ipsolv 27-Aug-2002 760 flow.iptime 30-Jul-2001 1.4k grid.exelem 16-Aug-2001 7.5k test_output.com 14-Aug-2001 15

Download the entire example:

Name                      Modified     Size

examples_8_84_845.tar.gz 19-Aug-2006 16k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:02:05 20160
cm-debugSuccessSat Mar 5 00:02:24 20160
mips-irix
cmSuccessSun Aug 19 01:29:43 20074
cm-debugSuccessWed Aug 15 01:22:21 20077
cm-debug-clear-mallocSuccessSat Aug 18 01:41:52 200723
cm-debug-clear-malloc7SuccessMon Aug 20 01:36:07 200733
cm64SuccessSun Aug 19 01:29:12 20074
cm64-debugSuccessTue Aug 21 01:23:08 20078
cm64-debug-clear-mallocSuccessThu Apr 1 10:50:51 20047
rs6000-aix
cmSuccessWed Mar 4 01:08:05 20091
cm-debugSuccessMon Mar 2 01:07:00 20091
cm64SuccessWed Mar 4 01:07:30 20091
cm64-debugSuccessTue Mar 3 01:12:14 20091
x86_64-linux
cmSuccessSun Mar 6 00:01:04 20161
cm-debugSuccessSat Mar 5 00:01:15 20161

Testing status by file:


Html last generated: Sun Mar 6 05:50:33 2016

Input last modified: Thu Dec 1 15:38:46 2005


CMISS Help / Examples / 8 / 84 / 845