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


CMISS Help / Examples / 3 / 33 / 334