[cszssj]

Here is a simple two dimensional case, the DO model is uesd, the inlet and outlet boundary is semitransparent. A source term is added to the solid phase (3e7 W/m3) .
When i use the default isotropic scattering phase function, the sum of radiation flux on the boundaries is equal to the integral of source term（3.6e5 W）,which is right. But when i use the UDF( copied from the user manual), the sum of radiation flux changes , which is wrong. Even I use the "ScatIso" UDF, which basically is isotropic scattering phase function, the answer still changes.
So please help me to solve this problem, thanks advance.



https://drive.google.com/file/d/0B-G...ew?usp=sharing