×

A hybrid transport-diffusion method for Monte Carlo radiative-transfer simulations. (English) Zbl 1111.82308

Summary: Discrete Diffusion Monte Carlo (DDMC) is a technique for increasing the efficiency of Monte Carlo particle-transport simulations in diffusive media. If standard Monte Carlo is used in such media, particle histories will consist of many small steps, resulting in a computationally expensive calculation. In DDMC, particles take discrete steps between spatial cells according to a discretized diffusion equation. Each discrete step replaces many small Monte Carlo steps, thus increasing the efficiency of the simulation. In addition, given that DDMC is based on a diffusion equation, it should produce accurate solutions if used judiciously. In practice, DDMC is combined with standard Monte Carlo to form a hybrid transport-diffusion method that can accurately simulate problems with both diffusive and non-diffusive regions. In this paper, we extend previously developed DDMC techniques in several ways that improve the accuracy and utility of DDMC for nonlinear, time-dependent, radiative-transfer calculations. The use of DDMC in these types of problems is advantageous since, due to the underlying linearizations, optically thick regions appear to be diffusive. First, we employ a diffusion equation that is discretized in space but is continuous in time. Not only is this methodology theoretically more accurate than temporally discretized DDMC techniques, but it also has the benefit that a particle’s time is always known. Thus, there is no ambiguity regarding what time to assign a particle that leaves an optically thick region (where DDMC is used) and begins transporting by standard Monte Carlo in an optically thin region. Also, we treat the interface between optically thick and optically thin regions with an improved method, based on the asymptotic diffusion-limit boundary condition, that can produce accurate results regardless of the angular distribution of the incident Monte Carlo particles. Finally, we develop a technique for estimating radiation momentum deposition during the DDMC simulation, a quantity that is required to calculate correct fluid motion in coupled radiation-hydrodynamics problems. With a set of numerical examples, we demonstrate that our improved DDMC method is accurate and can provide efficiency gains of several orders of magnitude over standard Monte Carlo.

MSC:

82C70 Transport processes in time-dependent statistical mechanics
80A20 Heat and mass transfer, heat flow (MSC2010)
78A40 Waves and radiation in optics and electromagnetic theory
82C80 Numerical methods of time-dependent statistical mechanics (MSC2010)
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Fleck, J. A.; Cummings, J. D., An Implicit Monte Carlo scheme for calculating time and frequency dependent nonlinear radiation transport, J. Comput. Phys., 8, 313 (1971) · Zbl 0229.65087
[2] Larsen, E. W.; Mercier, B., Analysis of a Monte Carlo method for nonlinear radiative transfer, J. Comput. Phys., 71, 50 (1987) · Zbl 0628.65128
[3] Martin, W. R.; Brown, F. B., Error modes in Implicit Monte Carlo, Trans. Amer. Nucl. Soc., 85, 508 (2002)
[4] Mosher, S. W.; Densmore, J. D., Stability and monotonicity conditions for linear, grey, 0-D Implicit Monte Carlo calculations, Trans. Amer. Nucl. Soc., 93, 520 (2005)
[5] Fleck, J. A.; Canfield, E. H., A random walk procedure for improving the computational efficiency of the Implicit Monte Carlo method for nonlinear radiation transport, J. Comput. Phys., 54, 508 (1984) · Zbl 0558.65096
[6] Giorla, J.; Sentis, R., A random walk method for solving radiative transfer equations, J. Comput. Phys., 70, 145 (1987) · Zbl 0616.65145
[7] T.J. Urbatsch, J.E. Morel, J. Gulick, Monte Carlo solution of spatially-discrete transport equations, Part II: Diffusion and transport-diffusion, in: Proceedings of the International Conference on Mathematics and Computation, Reactor Physics, and Environmental Analysis in Nuclear Applications, Madrid, Spain, 27-30 September 1999, vol. 1, p. 262.; T.J. Urbatsch, J.E. Morel, J. Gulick, Monte Carlo solution of spatially-discrete transport equations, Part II: Diffusion and transport-diffusion, in: Proceedings of the International Conference on Mathematics and Computation, Reactor Physics, and Environmental Analysis in Nuclear Applications, Madrid, Spain, 27-30 September 1999, vol. 1, p. 262.
[8] T.M Evans, T.J. Urbatsch, H. Lichtenstein, 1-D equilibrium Discrete Diffusion Monte Carlo, in: Proceedings of MC2000—International Conference, Lisbon, Portugal, July, 2000.; T.M Evans, T.J. Urbatsch, H. Lichtenstein, 1-D equilibrium Discrete Diffusion Monte Carlo, in: Proceedings of MC2000—International Conference, Lisbon, Portugal, July, 2000.
[9] Gentile, N. A., Implicit Monte Carlo Diffusion—an acceleration method for Monte Carlo time-dependent radiative transfer simulations, J. Comput. Phys., 172, 543 (2001) · Zbl 0986.65148
[10] J.D. Densmore, T.J. Urbatsch, T.M. Evans, M.W. Buksas, Discrete Diffusion Monte Carlo for grey Implicit Monte Carlo simulations, in: Proceedings of the International Topical Meeting on Mathematics and Computation, Supercomputing, Reactor Physics and Nuclear and Biological Applications, Avignon, France, 12-15 September 2005.; J.D. Densmore, T.J. Urbatsch, T.M. Evans, M.W. Buksas, Discrete Diffusion Monte Carlo for grey Implicit Monte Carlo simulations, in: Proceedings of the International Topical Meeting on Mathematics and Computation, Supercomputing, Reactor Physics and Nuclear and Biological Applications, Avignon, France, 12-15 September 2005.
[11] Larsen, E. W.; Pomraning, G. C.; Badham, V. C., Asymptotic analysis of radiative transfer problems, J. Quant. Spectrosc. Radiat. Transfer, 29, 285 (1983)
[12] Densmore, J. D., An improved method for treating Monte Carlo-diffusion interfaces, Trans. Amer. Nucl. Soc., 91, 139 (2004)
[13] Densmore, J. D., Interface methods for hybrid Monte Carlo-diffusion radiation transport simulations, Ann. Nucl. Energy, 33, 343 (2006)
[14] Habetler, G. J.; Matkowsky, B. J., Uniform asymptotic expansions in transport theory with small mean free paths, and the diffusion approximation, J. Math. Phys., 16, 846 (1975) · Zbl 0337.45015
[15] Bell, G. I.; Glasstone, S., Nuclear Reactor Theory (1985), Krieger Publishing: Krieger Publishing Malabar, FL
[16] Mihalas, D.; Weibel-Mihalas, B., Foundations of Radiation Hydrodynamics (1999), Dover Publications: Dover Publications Mineola, NY · Zbl 0651.76005
[17] Castor, J. I., Radiation Hydrodynamics (2004), Cambridge University Press: Cambridge University Press Cambridge, UK · Zbl 1323.78006
[18] Pomraning, G. C., The Equations of Radiation Hydrodynamics (1973), Pergamon Press: Pergamon Press Oxford, UK
[19] Lewis, E. E.; Miller, W. F., Computational Methods of Neutron Transport (1993), American Nuclear Society: American Nuclear Society La Grange Park, IL · Zbl 0594.65096
[20] Szilard, R. H.; Pomraning, G. C., Numerical transport and diffusion methods in radiative transfer, Nucl. Sci. Eng., 112, 256 (1992)
[21] Pomraning, G. C.; Foglesong, G. M., Transport-diffusion interfaces in radiative transfer, J. Comput. Phys., 32, 420 (1979) · Zbl 0439.45013
[22] C.J. Gesh, M.L. Adams, Finite element solutions of second-order forms of the transport equation at the interface between diffusive and non-diffusive regions, in: Proceedings of the International Meeting on Mathematical Methods for Nuclear Applications, Salt Lake City, UT, 9-13 September 2001.; C.J. Gesh, M.L. Adams, Finite element solutions of second-order forms of the transport equation at the interface between diffusive and non-diffusive regions, in: Proceedings of the International Meeting on Mathematical Methods for Nuclear Applications, Salt Lake City, UT, 9-13 September 2001.
[23] Davidson, G.; Densmore, J. D.; Prinja, A. K.; Morel, J. E., Asymptotically correct angular distributions for Monte Carlo-diffusion interfaces, Trans. Amer. Nucl. Soc., 94, 517 (2006)
[24] Densmore, J. D.; Davidson, G.; Carrington, D. B., Emissivity of discretized diffusion problems, Ann. Nucl. Energy, 33, 583 (2006)
[25] Andreev, E. S.; Kozmanov, M. Yu.; Rachilov, E. B., The maximum principle for a system of equations of energy and non-stationary radiation transfer, USSR Comput. Math. Math. Phys., 23, 104 (1983)
[26] Mercier, B., Application of accretive operators theory to the radiative transfer equations, SIAM J. Math. Anal., 18, 393 (1987) · Zbl 0654.47036
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.