1(
2  NAFEMS  Case for conduction/convection test T4
3  0.6m X 1.0 m rectangle  one 1.0m side insulated.
4  one 0.6m size at 100C, the rest are convecting to 0C with
5  750 W/m^2K. Conductivity is 52 W/mK.
6  Temp at the point (0.6,0.2) needs to be  18.3C
7
8  Osman Buyukisik  osman@fuse.net
9)
10
11echo -no
12number_of_space_dimensions 2
13condif_temperature
14end_initia
15
16start_define
17  insulated  geometry_line 4
18end_define
19start_define
20  fixed_temp geometry_line 3
21end_define
22start_define
23  convected_sides geometry_set 1
24end_define
25start_define
26  eps 0.01
27end_define
28
29group_type 0  -condif
30group_condif_conductivity 0 52.0
31group_condif_density      0 1.0
32
33geometry_line   1  0.6 0.0  0.6 1.0  eps
34geometry_line   2  0.6 1.0  0.0 1.0  eps
35fixed_temp         0.0 0.0  0.6 0.0  eps
36insulated          0.0 0.0  0.0 1.0  eps
37convected_sides  -geometry_line 1 -geometry_line 2
38
39condif_convection 0 750.0  0.0
40condif_convection_geometry 0 -convected_sides
41
42bounda_unknown 0  -fixed_temp -temp
43bounda_time    0  0. 100. 100. 100.
44
45post_point 1  0.6 0.2
46
47control_mesh_macro            10  -rectangle   0  4 6
48control_mesh_macro_parameters 10  0.3 0.5 0.6 1.0
49
50control_mesh_refine_globally  20 -p_refinement
51control_mesh_refine_globally  21 -p_refinement
52
53control_timestep              30  1.0 1.0
54
55target_item    1  -post_point_dof 1 -temp
56target_value   1  18.3 0.10
57end_data
58