Pollution Dispersion in a Harbour

Example 336: Pollution Dispersion in a Harbour

This example investigates the dispersion of a pollutant in a harbour by diffusion alone. We consider the two dimensional case only, ie. a depth average model is considered. The pollutant is introduced to the harbour as a continuous supply from a point source. Pollution entering the harbour at node number 4 (bottom right corner). An obstacle in the harbour (ie an island) is represented at node 12 with a zero BC.

The comfile run by this example is as follows:

#Example_336 Pollution Dispersion in an Open Harbour
#
fem define node;r;harbour;example  #Reads in grid of nodes
fem define base;r;;example         #Defines bilinear basis functions
fem define elem;r;;example         #Defines elements
fem define win;c                   #Calculates window size
fem draw line                      #Draws the lines of the elements
fem define equa;r;;example         #Define TIME INTEGRATION analysis of the ADVECTION DIFFUSION
#				   !equation. The solution is defined by GALERKIN FE's and a 
#				   !LINEAR analysis is used. Default all other options
fem define mate;r;;example         #Defines the material properties. Set constant and linear 
#				   !source terms to 0. The storage coefficient is set to 1.
#	 		 	   !Diffusion is both directions is set to 1 and the mass flow
#				   !is 0 in both directions.
fem de init;r;;example             #
fem de solve;r;;example            #Defines the solver with no solution output
fem solve                          #
fem dr field                       #Draws a field of the pollutant
fem define init;w;harbour-sol      #Save solution as it gets overwritten by
fem display history;harbour        #Displays the history of the nodes

# export geometry and field to cmgui
fem define init;r;harbour-sol      #Recover saved solution
fem export no;harbour as harbour
fem export elem;harbour as harbour
fem export no;pollution field as pollution
fem export elem;pollution field as pollution

Additional testing commands:

fem list node solution
#

Files used by this example are:

Name                     Modified     Size

example_336.com          02-Sep-2002  1.4k
harbour.cgi              10-Apr-2000  1.6k
harbour.exelem           10-Apr-2000  5.6k
harbour.exnode           10-Apr-2000  1.9k
harbour.form             10-Apr-2000  1.1k
harbour.html             10-Apr-2000  1.1k
harbour.ipbase           10-Apr-2000  1.1k
harbour.ipelem           10-Apr-2000  4.3k
harbour.ipequa           26-May-2003  1.3k
harbour.ipequa.old       18-Jun-2002  1.2k
harbour.iphist           10-Apr-2000  162k
harbour.ipinit           10-Apr-2000  1.3k
harbour.ipinit_template  10-Apr-2000  1.2k
harbour.ipmate           07-Apr-2003  2.6k
harbour.ipmate.old       10-Apr-2000  2.0k
harbour.ipnode           10-Apr-2000  3.4k
harbour.ipsolv           13-Apr-2007  1.7k
mesh.gif                 10-Apr-2000  3.7k
test_output.com          10-Apr-2000   30 
view.com                 03-Sep-2002  422

Download the entire example:

Name                      Modified     Size

examples_3_33_336.tar.gz  14-Apr-2007   42k

Testing status by version:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:01:50 2016	1
cm-debug	Success	Sat Mar 5 00:02:07 2016	0
mips-irix
cm	Success	Sun Aug 19 01:36:57 2007	8
cm-debug	Success	Wed Aug 15 01:28:52 2007	12
cm-debug-clear-malloc	Success	Sat Aug 18 01:33:49 2007	15
cm-debug-clear-malloc7	Success	Mon Aug 20 01:31:15 2007	15
cm64	Success	Sun Aug 19 01:37:01 2007	9
cm64-debug	Success	Tue Aug 21 01:27:47 2007	14
cm64-debug-clear-malloc	Success	Thu Apr 1 10:13:26 2004	6
rs6000-aix
cm	Success	Wed Mar 4 01:08:51 2009	1
cm-debug	Success	Mon Mar 2 01:07:12 2009	3
cm64	Success	Wed Mar 4 01:08:51 2009	2
cm64-debug	Success	Tue Mar 3 01:12:55 2009	2
x86_64-linux
cm	Success	Sun Mar 6 00:01:03 2016	1
cm-debug	Success	Sat Mar 5 00:01:13 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.

harbour.exelem

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.

harbour.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.

pollution.exelem

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.

pollution.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.

Graphical output from this problem is given here.

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

Input last modified: Fri Apr 13 10:32:37 2007

CMISS Help / Examples / 3 / 33 / 336