Numerical Performance of Stability Enhancing and Speed Increasing Steps in Radiative Transfer Solution Methods

• Methods for solving the radiative transfer problem, which are crucial for a number of sectors of industry, involve several numerical challenges. DORT2002 is a fast and numerically stable solution method for this problem, and has been implemented in MATLAB. This paper studies the numerical performance of DORT2002, in terms of stability, speed and accuracy. Focus is on the effects of the steps that are needed to make the method numerically efficient, and that differ from a naive implementation. Performance tests show that the steps that are included to improve stability and speed of DORT2002 are very successful. Together they give an unconditionally stable solution method to a problem previously considered numerically intractable, and decrease computation time compared to a naive implementation with a factor 1 000-10 000 in typical cases and with a factor up to and beyond 10 000 000 in extreme cases. It is also shown that the speed increasing steps are not introduced at the cost of reduced accuracy. Further studies and developments, that can have a positive impact on computation time, are suggested.

Ämnesord

NATURVETENSKAP  -- Matematik (hsv//swe)

NATURAL SCIENCES  -- Mathematics (hsv//eng)

Nyckelord

radiative transfer

solution method

numerical performance

stability

speed

accuracy

MATHEMATICS

MATEMATIK

Publikations- och innehållstyp

ref (ämneskategori)

art (ämneskategori)