Example 334: Contamination of Lake Pacific

The comfile run by this example is as follows:

#Example_334 Radiation contamination of Lake Pacific   
                                                                      
fem                             #Set the FEM environment 
fem define nodes;r;muroroa;example  #Read in the nodes to define Lake Pacific
fem define window;c             #Calculate the window size
fem define bases;r;;example     #Defines a bilinear basis function
fem define elements;r;;example  #Read in the elements
fem draw line                   #Draw the lines
fem draw nodes                  #Draw the nodes. The river enters the lake at
#                               !  node 1 and muroroa island is at node 13.
fem define equation;p           #Define a TIME INTEGRATION analysis of 
#                               !  the ADVECTION-DIFFUSION equation
#                               !  using GALERKIN FINITE ELEMENTS with 
#                               !  LINEAR ANALYSIS. Default all other options.
fem define material;p           #Define material properties. Set the constant
#                               !  and linear source terms to zero, the
#                               !  storage coefficient and both diffusion
#                               !  coefficients to 1, the mass flow in the
#                               !  x direction to 0.5 and the mass flow in 
#                               !  y direction to 1.0.
fem define initial;p            #Define some initial conditions. The initial
#                               !  solution is zero. Set node 13 to 1.0 (to
#                               !  represent maximum contamination), and
#                               !  node 1 to 0.0 (to represent minimum 
#                               !  contamination where the river enters the
#                               !  lake). Set all boundary nodes (other than
#                               !  1) to have integrated fluxes of zero 
#                               !  (the water cannot leave the lake).
#                               !  The boundary nodes are: 2,3,4,5,10,
#                               !  15,20,25,24,23,22,21,16,11,6
fem define solve;p              #Defines solver options. Default all options
#                               !  except set the final time to 1 seconds and
#                               !  the time increment to 0.05 seconds.
fem solve                       #Solve the time integration problem
fem draw field animate          #Show the contamination field over time.
#                               !  Note this may take a little while.

Files used by this example are:

Name                Modified     Size

example_334.com     20-Feb-2004  2.5k
muroroa.ipbase      10-Apr-2000  1.1k
muroroa.ipelem      10-Apr-2000  4.2k
muroroa.ipequa      26-May-2003  1.3k
muroroa.ipequa.old  18-Jun-2002  1.2k
muroroa.ipinit      15-Feb-2001  1.0k
muroroa.ipmate      07-Apr-2003  2.5k
muroroa.ipmate.old  15-Feb-2001  2.0k
muroroa.ipnode      10-Apr-2000  3.5k
muroroa.ipsolv      13-Apr-2007  1.8k

Download the entire example:

Name                      Modified     Size

examples_3_33_334.tar.gz  14-Apr-2007  5.0k

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

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