A line-by-line hybrid unstructured finite volume/Monte Carlo method for radiation transfer in 3D non-gray medium. (January 2018) - Sun, Hai-Feng
Highlights: A hybrid method incorporating the LBL model is developed for inhomogeneous and spectral media. The FVM is adopted for RTEs on spectral locations determined by the MCM. The TVD scheme is used to calculate the interface flux between control volumes. An algebraic multi-grid solution approach has been modified and used to accelerate the calculation. The presented method agrees well with the MCM with acceptable amount of spectral locations and computing time. Abstract: A hybrid method combing the unstructured finite volume method and the Monte Carlo method and incorporating the line-by-line model has been developed to simulate the radiative transfer in highly spectral and inhomogeneous medium. In this method, the unstructured finite volume method is adopted to solve the spectral radiative transfer equation at wave numbers or spectral locations determined by the Monte Carlo method. The Monte Carlo method takes effects by firstly defining the monotonic random number relations corresponding to the spectral emitted power density of every discretized element of the concerning medium, and then by reversing the spectral location through comparison of these relations with predefined random numbers. Through this Monte Carlo method, the actual number of spectral locations on which the spectral radiative transfer equations are solved may be reduced: only the spectral locations that have higher spectral emissive powers would be more possibly selected. To increase the performance of the presented method, the total variation diminishing scheme on unstructured grids is adopted in treating the spectral radiative intensity at interface between control volumes. And, the discretized radiative transfer equation is implicitly and iteratively solved by an algebraic multi-grid solution approach to accelerate the convergence of the equation. The presented method was applied to 3D homogeneous and inhomogeneous cases for the validation and performance studies. Results show that for both cases, the presented method agree well with pure Monte Carlo benchmark solutions with acceptable number of spectral locations and computing time.