Example 324: Steady state temperature distribution inside an annulus using BEM

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


The comfile run by this example is as follows:

#Example_324 Steady State Temperature distribution inside an annulus
#            using BEM
 
fem                              #Set the FEM environment
fem define node;r;bemannulus;example #Read in the bem annulus nodal coordinates
fem define base;r;;example           #Define a linear BE basis function.
fem list bases                   #List the basis function.
fem list bases family 1          #List the all the BE quadrature bases
fem define element;r;;example        #Define the bem 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 direct BEM.
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 FEM example in example_312.
#                                !  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 FEM example in example_312.
fem draw field                   #Draw the solution (temperature) field

Additional testing commands:

FEM list bases family 1
FEM list node solution
FEM list node flux

Files used by this example are:

Name               Modified     Size

example_324.com 10-Apr-2000 2.8k bemannulus.ipanal 10-Apr-2000 492 bemannulus.ipbase 10-Apr-2000 855 bemannulus.ipelem 10-Apr-2000 8.0k bemannulus.ipequa 26-May-2003 1.2k bemannulus.ipinit 10-Apr-2000 8.3k bemannulus.ipmate 10-Apr-2000 54 bemannulus.ipmesh 10-Apr-2000 796 bemannulus.ipnode 10-Apr-2000 5.3k bemannulus.ipsolv 13-Apr-2007 1.2k test_output.com 10-Apr-2000 66

Download the entire example:

Name                      Modified     Size

examples_3_32_324.tar.gz 14-Apr-2007 7.3k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:47 20161
cm-debugSuccessSat Mar 5 00:02:03 20161
mips-irix
cmSuccessSun Aug 19 01:25:52 20073
cm-debugSuccessWed Aug 15 01:20:51 20076
cm-debug-clear-mallocSuccessSat Aug 18 01:27:12 200712
cm-debug-clear-malloc7SuccessMon Aug 20 01:24:59 200711
cm64SuccessSun Aug 19 01:25:28 20075
cm64-debugSuccessTue Aug 21 01:21:46 20076
cm64-debug-clear-mallocSuccessThu Apr 1 10:13:03 20048
rs6000-aix
cmSuccessWed Mar 4 01:07:03 20091
cm-debugSuccessMon Mar 2 01:06:17 20091
cm64SuccessWed Mar 4 01:07:03 20090
cm64-debugSuccessTue Mar 3 01:11:52 20091
x86_64-linux
cmSuccessSun Mar 6 00:01:01 20160
cm-debugSuccessSat Mar 5 00:01:12 20160

Testing status by file:


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

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


CMISS Help / Examples / 3 / 32 / 324