cannot get the power transfer between two SiO2Spheres

[ZhihaoJia16]

Hi, Homer

The results are all 0.000000e+00 when I run the example about the power transfer between two SiO2Spheres, http://homerreid.github.io/buff-em-documentation/examples/SiO2Spheres/
The results are same for methods of --OPFT, --EMTPFT, and --MomentPFT .
RunScript

#!/bin/bash

ARGS=""
ARGS="${ARGS} --geometry SiO2Spheres_677.buffgeo"
ARGS="${ARGS} --OmegaFile OmegaFile"
ARGS="${ARGS} --EMTPFT"
ARGS="${ARGS} --PAbs"

buff-neq ${ARGS}

results

# buff-neq run on jia-Precision-Tower-7810 (12/09/16::21:48:36)
# data file columns: 
# 1 omega 
# 2 transform tag
# 3 (sourceObject,destObject) 
# 4 absorbed power spectral density
# 5 radiated power spectral density
# 6 x-force flux spectral density
# 7 y-force flux spectral density
# 8 z-force flux spectral density
# 9 x-torque flux spectral density
# 10 y-torque flux spectral density
# 11 z-torque flux spectral density
DEFAULT 1.000000e-02 11 -1.670797e-09 1.670797e-09 2.965621e-11 -1.909764e-11 9.200753e-08 -1.000000e-02 2.561126e-10 -1.000000e-04 
DEFAULT 1.000000e-02 12 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 
DEFAULT 1.000000e-02 21 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 
DEFAULT 1.000000e-02 22 -1.670758e-09 1.670758e-09 1.992007e-14 -5.163125e-14 2.739386e-14 -2.334779e-12 5.712766e-13 1.269466e-12 
DEFAULT 1.258925e-02 11 -4.201814e-09 4.201814e-09 2.948099e-11 -1.836088e-11 9.224337e-08 -1.258925e-02 2.572150e-10 -1.584892e-04 
DEFAULT 1.258925e-02 12 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 

I hve looked through the codes, but i still do not know the reasons.
Is it a problem about SCUFF_PFT_EMT (https://github.com/HomerReid/buff-em/blob/master/src/applications/buff-neq/CreateBNEQData.cc#L177)? Or anything else?

Thank you.

Zhihao

[linxiaozhu]

Hi Homer

I met similar issue as discussed by Zhihao. When I tried to do an example run for calculating radiative heat transfer between two SiO2 spheres, the item corresponding to the heat transfer from one sphere to another is shown as zero.

Do you know if there is something I am missing that cause this problem?

Thanks,
--Linxiao