Lax Wendroff finite difference flow solution

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:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:02:05 2016	0
cm-debug	Success	Sat Mar 5 00:02:24 2016	0
mips-irix
cm	Success	Sun Aug 19 01:29:43 2007	4
cm-debug	Success	Wed Aug 15 01:22:21 2007	7
cm-debug-clear-malloc	Success	Sat Aug 18 01:41:52 2007	23
cm-debug-clear-malloc7	Success	Mon Aug 20 01:36:07 2007	33
cm64	Success	Sun Aug 19 01:29:12 2007	4
cm64-debug	Success	Tue Aug 21 01:23:08 2007	8
cm64-debug-clear-malloc	Success	Thu Apr 1 10:50:51 2004	7
rs6000-aix
cm	Success	Wed Mar 4 01:08:05 2009	1
cm-debug	Success	Mon Mar 2 01:07:00 2009	1
cm64	Success	Wed Mar 4 01:07:30 2009	1
cm64-debug	Success	Tue Mar 3 01:12:14 2009	1
x86_64-linux
cm	Success	Sun Mar 6 00:01:04 2016	1
cm-debug	Success	Sat Mar 5 00:01:15 2016	1

Testing status by file:

cmiss_test.log

i686-linux
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
mips-irix
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
Success	cm-debug-clear-malloc:	cmiss_test.log.retain.
Success	cm-debug-clear-malloc7:	cmiss_test.log.retain.
Success	cm64:	cmiss_test.log.retain.
Success	cm64-debug:	cmiss_test.log.retain.
Success	cm64-debug-clear-malloc:	cmiss_test.log.retain.
rs6000-aix
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
Success	cm64:	cmiss_test.log.retain.
Success	cm64-debug:	cmiss_test.log.retain.
x86_64-linux
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.

cmiss_test.out

i686-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
mips-irix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
rs6000-aix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
x86_64-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.

0.004.exnode

i686-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
mips-irix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
rs6000-aix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
x86_64-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.

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