Example 312: Steady state temperature distribution inside an annulus using FEM

This example calculates the steady state temperature distribution inside an annulus using FEM.

The comfile run by this example is as follows:

#Example_312 Steady State Temperature distribution inside an annulus
#            using FEM
 
fem                              #Set the FEM environment
fem define node;r;femannulus;example #Read in the fem annulus nodal coordinates
fem define base;r;;example           #Define a bilinear basis function.
fem list bases                   #List the basis function.
fem define element;r;;example        #Define the fem annulus elements
fem define window;c              #calculates the size of the window to fit
#                                !  the mesh
fem draw axes                    #Draw axes
fem draw node                    #Draw the nodes
fem draw line                    #Draw the lines
fem define equation;r;;example       #Defines the equation to be solved.  We are
#                                !  undertaking a static analysis involving
#                                !  Laplaces equation and we are going to
#                                !  solve it using Galerkin FEM.
fem list equation                #List the equation
fem define material;r;;example       #Define the materials. Since we are dealing
#                                !  with Laplaces equation there are no
#                                !  material properties to be defined.
fem define initial;r;;example        #Define the boundary conditions to be used.
#                                !  The conditions consist of a specified
#                                !  temperature distribution around the inside
#                                !  of the annulus and zero flux on the outside
#                                !  of the annulus with the exception of zero
#                                !  temperature at node 17.
fem draw node                    #Draw the nodes
fem list initial                 #List the initial conditions.
fem define solve;r;;example          #Define the solution options.
fem solve                        #Solves the problem.
fem list node solution           #List the node solution. Compare these with
#                                !  the equivalent BEM example in example_324.
#                                !  Note: the analytic answers for selected
#                                !  (the problem is symmetric) nodes is:
#                                !  Node 17 :  0.00000
#                                !  Node 18 : -0.06091
#                                !  Node 19 : -0.11254
#                                !  Node 20 : -0.14704
#                                !  Node 21 : -0.15915
fem list node flux               #List the node fluxes. Compare these with
#                                !  the equivalent BEM example in example_324.
fem draw field                   #Draw the temperature field

Additional testing commands:

FEM list node solution
FEM list node flux

Files used by this example are:

Name               Modified     Size

example_312.com    10-Apr-2000  2.7k
femannulus.ipanal  10-Apr-2000  416 
femannulus.ipbase  10-Apr-2000  1.1k
femannulus.ipelem  10-Apr-2000  4.2k
femannulus.ipequa  26-May-2003  1.3k
femannulus.ipinit  10-Apr-2000  8.6k
femannulus.ipmate  10-Apr-2000   54 
femannulus.ipmesh  10-Apr-2000  717 
femannulus.ipnode  10-Apr-2000  5.3k
femannulus.ipsolv  13-Apr-2007  1.1k
test_output.com    10-Apr-2000   42 

Download the entire example:

Name                      Modified     Size

examples_3_31_312.tar.gz  14-Apr-2007  6.8k

Testing status by version:

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.

  
  

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

Input last modified: Fri Apr 13 10:31:34 2007