Poor fitting results in ‘scuff-neq’

Hi, I am using scuff-neq to reproduce the radiative heat transfer between two gold circular plates in ' Fluctuating-surface-current formulation of radiative heat transfer: Theory and applications '. Where L = 0.2 μm, T1 = 300 K, T2 = 0 K.

For case of 2 R / L = 0.5 μm, the results as follows on the left:

My results can not fit with the reference value even if the mesh fineness lc = 0.015 μm.
When lc < 0.015 μm, the program do not stop even if omega is larger than 100.

Please give me some comments.

Other information ares as follows: In the Figure on the left, the ratio of the calculated value(lc = 0.02 μm) to the reference value as follows:

5.097101546
1.460421075
0.672121506
0.456759251
0.388060778
0.427822662
0.401782079
0.445790497
0.429814165
0.498550455
0.509916803
0.627869996
0.395448914
0.188555675

Cylinder.geo

R=0.05;   // disc radius
T=0.2;   // disc thickness

lc = 0.015;

Point(100) = {       0,   0,    0, lc };
Point(101) = {       R,   0,    0, lc };
Point(102) = {       0,   R,    0, lc };
Point(103) = {      -R,   0,    0, lc };
Point(104) = {       0,  -R,    0, lc };

Circle(101)    = { 101, 100, 102 };
Circle(102)    = { 102, 100, 103 };
Circle(103)    = { 103, 100, 104 };
Circle(104)    = { 104, 100, 101 };

Extrude{0,0,-T} { Line{101, 102, 103, 104}; }

Line Loop(1) = { 101, 102, 103, 104 };
Line Loop(2) = { 105, 109, 113, 117 };

Ruled Surface(1) = { 1 };
Ruled Surface(2) = { 2 };

Cylinders.scuffgeo

MATERIAL Gold
    wp = 1.37e16;
    gamma = 5.32e13;
    Eps(w) = 1 - wp^2 / (w * (w + i*gamma));
ENDMATERIAL

OBJECT UpperCylinder
    MESHFILE Cylinder_1461.msh
    MATERIAL Gold
        ROTATED 180 ABOUT 0 1 0
    DISPLACED 0 0 0.0
ENDOBJECT

OBJECT LowerCylinder
    MESHFILE Cylinder_1461.msh
    MATERIAL Gold
ENDOBJECT

Cylinders.trans

TRANS 0.02 OBJECT UpperCylinder DISP 0 0 0.02
TRANS 0.08 OBJECT UpperCylinder DISP 0 0 0.08
TRANS 0.16 OBJECT UpperCylinder DISP 0 0 0.16
TRANS 0.20 OBJECT UpperCylinder DISP 0 0 0.20
TRANS 0.40 OBJECT UpperCylinder DISP 0 0 0.40
TRANS 0.60 OBJECT UpperCylinder DISP 0 0 0.60
TRANS 0.80 OBJECT UpperCylinder DISP 0 0 0.80
TRANS 1.60 OBJECT UpperCylinder DISP 0 0 1.60
TRANS 2.00 OBJECT UpperCylinder DISP 0 0 2.00
TRANS 4.00 OBJECT UpperCylinder DISP 0 0 4.00

RunScript

#!/bin/bash

ARGS=""
ARGS="${ARGS} --geometry Cylinders.scuffgeo"
ARGS="${ARGS} --TransFile Cylinders.trans"
ARGS="${ARGS} --Temperature UpperCylinder 300"
ARGS="${ARGS} --OPFT"
ARGS="${ARGS} --Power"

scuff-neq ${ARGS}

Best, Thank you in advance.

HomerReid / scuff-em

Poor fitting results in ‘scuff-neq’ #89