( 
  NAFEMS  Case for conduction/convection test T4
  0.6m X 1.0 m rectangle  one 1.0m side insulated.
  one 0.6m size at 100C, the rest are convecting to 0C with
  750 W/m^2K. Conductivity is 52 W/mK. 
  Temp at the point (0.6,0.2) needs to be  18.3C 

  Osman Buyukisik  osman@fuse.net
)

echo -no
number_of_space_dimensions 2
condif_temperature
end_initia

start_define
  insulated  geometry_line 4
end_define
start_define
  fixed_temp geometry_line 3
end_define
start_define
  convected_sides geometry_set 1
end_define
start_define
  eps 0.01
end_define

group_type 0  -condif
group_condif_conductivity 0 52.0
group_condif_density      0 1.0

geometry_line   1  0.6 0.0  0.6 1.0  eps
geometry_line   2  0.6 1.0  0.0 1.0  eps
fixed_temp         0.0 0.0  0.6 0.0  eps
insulated          0.0 0.0  0.0 1.0  eps
convected_sides  -geometry_line 1 -geometry_line 2

condif_convection 0 750.0  0.0
condif_convection_geometry 0 -convected_sides

bounda_unknown 0  -fixed_temp -temp
bounda_time    0  0. 100. 100. 100.

post_point 1  0.6 0.2

control_mesh_macro            10  -rectangle   0  4 6
control_mesh_macro_parameters 10  0.3 0.5 0.6 1.0 

control_mesh_refine_globally  20 -p_refinement
control_mesh_refine_globally  21 -p_refinement

control_timestep              30  1.0 1.0

target_item    1  -post_point_dof 1 -temp
target_value   1  18.3 0.10
end_data
