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

Patricia Cargo; Gérald Samba

ESAIM: Mathematical Modelling and Numerical Analysis (2010)

  • Volume: 44, Issue: 6, page 1193-1224
  • ISSN: 0764-583X

Abstract

top
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.

How to cite

top

Cargo, Patricia, and Samba, Gérald. "Resolution of the time dependent Pn equations by a Godunov type scheme having the diffusion limit." ESAIM: Mathematical Modelling and Numerical Analysis 44.6 (2010): 1193-1224. <http://eudml.org/doc/250800>.

@article{Cargo2010,
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. },
author = {Cargo, Patricia, Samba, Gérald},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis},
keywords = {Time-dependent transport; Pn model; diffusion limit; finite volume method; Riemann solver; time-dependent transport; -model; neutron transport equations},
language = {eng},
month = {10},
number = {6},
pages = {1193-1224},
publisher = {EDP Sciences},
title = {Resolution of the time dependent Pn equations by a Godunov type scheme having the diffusion limit},
url = {http://eudml.org/doc/250800},
volume = {44},
year = {2010},
}

TY - JOUR
AU - Cargo, Patricia
AU - Samba, Gérald
TI - Resolution of the time dependent Pn equations by a Godunov type scheme having the diffusion limit
JO - ESAIM: Mathematical Modelling and Numerical Analysis
DA - 2010/10//
PB - EDP Sciences
VL - 44
IS - 6
SP - 1193
EP - 1224
AB - 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.
LA - eng
KW - Time-dependent transport; Pn model; diffusion limit; finite volume method; Riemann solver; time-dependent transport; -model; neutron transport equations
UR - http://eudml.org/doc/250800
ER -

References

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

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

                
Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.