I made a simple case of an 0.8 m x 0.8 m x 0.01 m plate on a 5 cm grid. rho=1000, c=1, and k=1000. I put an 0.2 m x 0.2 m radiant panel 0.4 m away which give a peak flux of ~5 kW/m2 at the center.
I ran three variants of the case with CELL_SIZE=0.05, 0.01, and 0.005. In a full blown 3D model, I would expect that as I reduce cell size that the temperature at the center of the plate should assymptotically increase. Here at 0.05 the center is 56.3, at 0.01 it is 46.3, and at 0.005 it is 60.3.
I had devices to measure the net flux on the front and back of the plate and the plate enthalpy. I didn't measure on the side since this is a thin obstruction; however, the total side area side is 3 % the area of the back and cooler overall than the back. The back loss was ~3 % of the front incident so at most the sides are ~0.1 %.
For the 0.05 case the enthalpy increased 72.85 and the net wall was 72.82 for a difference of 0.03
For the 0.005 case the values were 72.85 and 72.76 for a difference of 0.09
For the 0.01 case the values were 72.90 and 72.57 for a difference of 0.33.
Not accounting for the sides can't be more than about 0.07 so the 0.01 case error of 0.33 is much more than execpted. The 0.05 and 0.005 cases have an error that appears to be just mostly not including the sides.
plate_100_p05p01.txt
plate_100_p05p005.txt
plate_100_p05p05.txt
I made a simple case of an 0.8 m x 0.8 m x 0.01 m plate on a 5 cm grid. rho=1000, c=1, and k=1000. I put an 0.2 m x 0.2 m radiant panel 0.4 m away which give a peak flux of ~5 kW/m2 at the center.
I ran three variants of the case with CELL_SIZE=0.05, 0.01, and 0.005. In a full blown 3D model, I would expect that as I reduce cell size that the temperature at the center of the plate should assymptotically increase. Here at 0.05 the center is 56.3, at 0.01 it is 46.3, and at 0.005 it is 60.3.
I had devices to measure the net flux on the front and back of the plate and the plate enthalpy. I didn't measure on the side since this is a thin obstruction; however, the total side area side is 3 % the area of the back and cooler overall than the back. The back loss was ~3 % of the front incident so at most the sides are ~0.1 %.
For the 0.05 case the enthalpy increased 72.85 and the net wall was 72.82 for a difference of 0.03
For the 0.005 case the values were 72.85 and 72.76 for a difference of 0.09
For the 0.01 case the values were 72.90 and 72.57 for a difference of 0.33.
Not accounting for the sides can't be more than about 0.07 so the 0.01 case error of 0.33 is much more than execpted. The 0.05 and 0.005 cases have an error that appears to be just mostly not including the sides.
plate_100_p05p01.txt
plate_100_p05p005.txt
plate_100_p05p05.txt