( 
Solid 2D Axisymm Thermal model with heat generation and convection.
10mm radiusX10mm thick slab [ thickness does not matter: infinite cylinder ]. 
1.0 W/mm^3 heat gen in the element 
Tc = 0  in the surrounding fluid, with heat transfer
coeff = 1.0 W/mm^2K
top & bottom are adiabatic
k=1, c=0, density=0
The exact solution is 
T = 30 - 0.25 * R^2   where R is the distance from the center.
At R=0, T=30 
At R=10, T=5
This axisymm model includes the center. One element model
does not give an accurate center temp (peak). Hence we need to
use the mesh refinement command. Also the BCs are defined on the
geometry not on the initial mesh.

Osman Buyukisik        osman@fuse.net

)

echo -no
number_of_space_dimensions 2
condif_temperature
end_initia

node 1  0 0
node 2  10 0
node 3  0  10
node 4  10 10

element 1  -quad4 1 2 3 4

group_type 0  -condif
group_condif_conductivity 0 1
group_axisymmetric 0 -yes

geometry_line          0  0.0 0.0 10.0 0.0 0.1
geometry_line          1  10.0 0.0 10.0 10.0 0.1
geometry_line          2  0.0 10.0 10.0 10.0 0.1

geometry_quadrilateral 0  0 0 10 0 0 10 10 10 0.1

(only one side  convects)
condif_convection  0   1.0 0.0
condif_convection_geometry 0  -geometry_line 1

(heat gen in element 1)
force_element_volume          0  1.0
force_element_volume_geometry 0  -geometry_quadrilateral 0

control_mesh_refine_globally  1  -h_refinement
control_mesh_refine_globally  2  -h_refinement
control_mesh_refine_globally  3  -h_refinement

control_timestep             10  1. 1.

post_point                    1  0 0
post_point                    2  10 0
post_point                    3  0 10
post_point                    4  10 10

(
control_print                10  -time_current -node_dof
)

target_item    1  -post_point_dof 1 -temp
target_value   1  30 1.
target_item    2  -post_point_dof 2 -temp
target_value   2  5 .1
target_item    3  -post_point_dof 3 -temp
target_value   3  30 1.
target_item    4  -post_point_dof 4 -temp
target_value   4  5 .1
end_data
