Radiation heat transfer
 

Dear all,

I have written a computer code for radiation transfer equation(RTE) in the cylindrical system.I used the discrete ordinate method(DOM)with the S4 direction. If you have worked in this case, we have a problem in the curved system so we should use a special difference in RTE and special algorithm for it. For solving, I start for solving RTE in zero weight and ETA=0, and MU=1-ZITA**2, this procedure is for finding I1/2. and after that I start solving from minus Mu, and so on. My resullt isn't correct and I have about 25% diffrence in some casese with correct results, If you know something about this method please send me.

Thank you, Mehdi