Hostname: page-component-8448b6f56d-c4f8m Total loading time: 0 Render date: 2024-04-24T12:22:41.949Z Has data issue: false hasContentIssue false
  • English
  • Français

Resolution of the time dependent Pn equations by a Godunov type scheme having the diffusion limit

Published online by Cambridge University Press:  15 April 2010

Patricia Cargo  and
Gérald Samba
Patricia Cargo
Affiliation:
CEA, DAM, DIF, 91297 Arpajon, France. gerald.samba@cea.fr
Gérald Samba
Affiliation:
CEA, DAM, DIF, 91297 Arpajon, France. gerald.samba@cea.fr
Get access

Abstract

We consider the Pn model to approximate the time dependent transport equation in one dimension of space. In a diffusive regime, the solution of this system is solution of a diffusion equation. We are looking for a numerical scheme having the diffusion limit property: in a diffusive regime, it has to give the solution of the limiting diffusion equation on a mesh at the diffusion scale. The numerical scheme proposed is an extension of the Godunov type scheme proposed by Gosse to solve the P1 model without absorption term. It requires the computation of the solution of the steady state Pn equations. This is made by one Monte-Carlo simulation performed outside the time loop. Using formal expansions with respect to a small parameter representing the inverse of the number of mean free path in each cell, the resulting scheme is proved to have the diffusion limit. In order to avoid the CFL constraint on the time step, we give an implicit version of the scheme which preserves the positivity of the zeroth moment.

Keywords

Time-dependent transportPn modeldiffusion limitfinite volume methodRiemann solver
Type
Research Article
Information
ESAIM: Mathematical Modelling and Numerical Analysis , Volume 44 , Issue 6 , November 2010 , pp. 1193 - 1224
Copyright
© EDP Sciences, SMAI, 2010

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

G. Bell and S. Glasstone, Nuclear Reactor Theory. Van Nostrand, Princeton (1970).
Berthon, C., Charrier, P. and Dubroca, B., HLLC, An scheme to solve the M 1 model of radiative transfer in two space dimensions. J. Sci. Comp. 31 (2007) 347–389. CrossRef
T.A. Brunner, Riemann solvers for time-dependent transport based on the maximum entropy and spherical harmonics closures. Ph.D. Thesis, University of Michigan (2000).
Buet, C. and Despres, B., Asymptotic preserving and positive schemes for radiation hydrodynamics. J. Comput. Phys. 215 (2006) 717–740. CrossRef
C. Buet and S. Cordier, Asymptotic preserving scheme and numerical methods for radiative hydrodynamic models. C. R. Acad. Sci., Sér. 1 Math. 338 (2004) 951–956.
K.M. Case and P.F. Zweifel, Linear Transport Theory. Addison-Wesley Publishing Co., Inc. Reading (1967).
S. Chandrasekhar, Radiative transfer. Dover, New York (1960).
R. Dautray and J.L. Lions, Analyse mathématique et calcul numérique. Chap. 21, Masson, Paris (1988).
J.J. Duderstadt and W.R. Martin, Transport theory. Wiley-Interscience, New York (1979).
E.M. Gelbard, Simplified spherical harmonics equations and their use in shielding problems. Technical report WAPD-T-1182, Bettis Atomic Power Laboratory, USA (1961).
L. Gosse and G. Toscani, An asymptotic-preserving well-balanced scheme for the hyperbolic heat equations. C. R. Acad. Sci., Sér. 1 Math. 334 (2002) 337–342.
Gosse, L. and Toscani, G., Space localization and well-balanced schemes for discrete kinetic models in diffusive regimes. SIAM J. Numer. Anal. 41 (2003) 641–658. CrossRef
Greenberg, J.M. and Leroux, A.Y., A well-balanced scheme for the numerical processing of source terms in hyperbolic equations. SIAM J. Numer. Anal. 33 (1996) 1–16. CrossRef
Keller, H.B., Approximate solutions of transport problems. II. Convergence and applications of the discrete-ordinate method. J. Soc. Indust. Appl. Math. 8 (1960) 43–73. CrossRef
Larsen, E.W., On numerical solutions of transport problems in the diffusion limit. Nucl. Sci. Eng. 83 (1983) 90. CrossRef
Larsen, E.W. and Keller, J.B., Asymptotic solution of neutron transport problems for small mean free paths. J. Math. Phys. 15 (1974) 75. CrossRef
Larsen, E.W., Pomraning, G.C. and Badham, V.C., Asymptotic analysis of radiative transfer problems. J. Quant. Spectrosc. Radiat. Transfer 29 (1983) 285. CrossRef
Lathrop, K.D., Ray effects in discrete ordinates equations. Nucl. Sci. Eng. 32 (1968) 357. CrossRef
Levermore, C.D., Relating Eddington factors to flux limiters. J. Quant. Spec. Rad. Transfer. 31 (1984) 149–160. CrossRef
R. McClarren, J.P. Holloway, T.A. Brunner and T. Melhorn, An implicit Riemann solver for the time-dependent P n equations, in International Topical Meeting on Mathematics and Computation, Supercomputing, Reactor Physics and Nuclear and Biological Applications, American Nuclear Society, Avignon, France (2005).
R. McClarren, J.P. Holloway and T.A. Brunner, Establishing an asymptotic diffusion limit for Riemann solvers on the time-dependent P n equations, in International Topical Meeting on Mathematics and Computation, Supercomputing, Reactor Physics and Nuclear and Biological Applications, American Nuclear Society, Avignon, France (2005).
McClarren, R., Holloway, J.P. and Brunner, T.A., On solutions to the P n equations for thermal radiative transfer. J. Comput. Phys. 227 (2008) 2864–2885. CrossRef
Pomraning, G.C., Diffusive limits for linear transport equations. Nucl. Sci. Eng. 112 (1992) 239–255. CrossRef