STABLE BOUNDARY CONDITIONS AND DISCRETIZATION FOR PN EQUATIONS

Jonas Bünger, Neeraj Sarna, Manuel Torrilhon

1. Applied and Computational Mathematics (ACoM), Department of Mathematics, RWTH Aachen University, Schinkelstr 2, 52062 Aachen, Germany
• Received:2019-10-17 Online:2022-11-15 Published:2022-11-18
• Contact: Jonas Bünger, Email: buenger@acom.rwth-aachen.de
• Supported by:
The authors acknowledge funding of the German Research Foundation (DFG) under grant TO 414/4-1.

Jonas Bünger, Neeraj Sarna, Manuel Torrilhon. STABLE BOUNDARY CONDITIONS AND DISCRETIZATION FOR PN EQUATIONS[J]. Journal of Computational Mathematics, 2022, 40(6): 977-1003.

A solution to the linear Boltzmann equation satisfies an energy bound, which reflects a natural fact: The energy of particles in a finite volume is bounded in time by the energy of particles initially occupying the volume augmented by the energy transported into the volume by particles entering the volume over time. In this paper, we present boundary conditions (BCs) for the spherical harmonic (PN) approximation, which ensure that this fundamental energy bound is satisfied by the PN approximation. Our BCs are compatible with the characteristic waves of PN equations and determine the incoming waves uniquely. Both, energy bound and compatibility, are shown on abstract formulations of PN equations and BCs to isolate the necessary structures and properties. The BCs are derived from a Marshak type formulation of BC and base on a non-classical even/odd-classification of spherical harmonic functions and a stabilization step, which is similar to the truncation of the series expansion in the PN method. We show that summation by parts (SBP) finite differences on staggered grids in space and the method of simultaneous approximation terms (SAT) allows to maintain the energy bound also on the semi-discrete level.

CLC Number:

 [1] C. Cercignani, The Boltzmann Equation and its Applications. Applied Mathematical Sciences, Springer New York, 2012.[2] K. Case and P. Zweifel, Linear Transport Theory. Addison-Wesley Series in Nuclear Engineering, Addison-Wesley Pub. Co., 1967.[3] E.W. Larsen, M.M. Miften, B.A. Fraass and I.A.D. Bruinvis, Electron dose calculations using the Method of Moments, Med. Phys., 24:1(1997), 111-125.[4] C. Ringhofer, C. Schmeiser and A. Zwirchmayer, Moment methods for the semiconductor boltzmann equation in bounded Position domains, SIAM Journal on Numerical Analysis, 39:3(2001), 1078-1095.[5] N. Sarna and M. Torrilhon, On Stable Wall Boundary Conditions for the Hermite discretization of the linearised Boltzmann equation, Journal of Statistical Physics, 170(2018), 101-126.[6] A. S. Rana and H. Struchtrup, Thermodynamically admissible boundary conditions for the regularized 13 moment equations, Physics of Fluids, 28:2(2016), 027105.[7] H. Struchtrup and M. Torrilhon, H Theorem, Regularization, and boundary conditions for linearized 13 moment equations, Physical Review Letters, 99(2007), 014502.[8] B. Davison, Neutron Transport Theory. International Series of Monographs on Physics, Clarendon Press, 1958.[9] K. Atkinson and W. Han, Spherical Harmonics and Approximations on the Unit Sphere: An Introduction, vol. 2014. Springer, 012012.[10] C. Müller, Spherical Harmonics. Springer-Verlag Berlin; New York, 1966.[11] J.C. Mark, The spherical harmonics method, part I, Tech. Report MT 92, National Research Council of Canada, 1944.[12] J.C. Mark, The spherical harmonics method, part II, Tech. Report MT 97, National Research Council of Canada, 1945.[13] R.E. Marshak, Note on the spherical harmonic method as applied to the Milne problem for a sphere, Phys. Rev., 71(1947), 443-446.[14] M. Williams, Mathematical Methods in Particle Transport Theory, Wiley-Interscience, 1971.[15] H. Egger and M. Schlottbom, A mixed variational framework for the radiative transfer equation, Mathematical Models and Methods in Applied Sciences, 22:3(2012), 1150014.[16] M. Schlottbom and H. Egger, A class of Galerkin schemes for time-dependent radiative transfer, SIAM Journal on Numerical Analysis, 54:6(2016), 3577-3599.[17] H.O. Kreiss and G. Scherer, Finite element and finite difference methods for hyperbolic partial differential equations, in Mathematical Aspects of Finite Elements in Partial Differential Equations (C. de Boor, ed.), 195-212, Academic Press, 1974.[18] B. Strand, Summation by parts for finite difference approximations for d=dx, Journal of Computational Physics, 110:1(1994), 47-67.[19] P. Olsson, Summation by parts, projections, and stability. I, Math. Comput., 64(1995), 1035-1065.[20] D.C.D.R. Fernández, J.E. Hicken and D.W. Zingg, Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations, Computers & Fluids, 95(2014), 171-196.[21] M.H. Carpenter, D. Gottlieb and S. Abarbanel, Time-stable boundary conditions for finitedifference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes, Journal of Computational Physics, 111:2(1994), 220-236.[22] B. Seibold and M. Frank, StaRMAP — A second order staggered grid method for spherical harmonics moment equations of radiative transfer, ACM Trans. Math. Softw., 41(2014), 4:1-4:28. Stable Boundary Conditions and Discretization for PN Equations 1003[23] O. O’Reilly, T. Lundquist, E.M. Dunham and J. Nordström, Energy stable and high-order-accurate finite difference methods on staggered grids, Journal of Computational Physics, 346(2017), 572-589.[24] L. Reimer, Scanning Electron Microscopy, Berlin: Springer, 2nd ed., 1998.[25] N.W.M. Ritchie, A new Monte Carlo application for complex sample geometries, Surface and Interface Analysis, 37:11(2005), 1006-1011.[26] N.W.M. Ritchie, Spectrum simulation in DTSA-II, Microscopy and Microanalysis, 15:5(2009), 454-522.[27] J. Howlet, Neutron Transport Theory. By B. Davison and J.B. Sykes. 75s. 1957. (Oxford University Press), The Mathematical Gazette, 43:345(1959), 221-222.[28] F. Salvat, M.J. Berger, A. Jablonski, I.K. Bronic, J. Mitroy, C.J. Powell and L. Sanche, Elastic scattering of electrons and positrons, vol. 7 of ICRU Report. Bethesda, MD: ICRU, 2007.[29] F. Salvat, A. Jablonski, and C.J. Powell, elsepa-Dirac partial-wave calculation of elastic scattering of electrons and positrons by atoms, positive ions and molecules, Computer Physics Communications, 165:2(2005), 157-190.[30] A. Jablonski, F. Salvat and C. Powell, NIST Electron Elastic-Scattering Cross-Section DatabaseVersion 3:1, National Institute of Standards and Technology, 2003.
 [1] Jie Chen, Zhengkang He, Shuyu Sun, Shimin Guo, Zhangxin Chen. EFFICIENT LINEAR SCHEMES WITH UNCONDITIONAL ENERGY STABILITY FOR THE PHASE FIELD MODEL OF SOLID-STATE DEWETTING PROBLEMS [J]. Journal of Computational Mathematics, 2020, 38(3): 452-468. [2] Kun Li, Youngju Lee, Christina Starkey. A NEW BOUNDARY CONDITION FOR RATE-TYPE NON-NEWTONIAN DIFFUSIVE MODELS AND THE STABLE MAC SCHEME [J]. Journal of Computational Mathematics, 2018, 36(4): 605-626. [3] Jie Yang, Guannan Zhang, Weidong Zhao. A FIRST-ORDER NUMERICAL SCHEME FOR FORWARD-BACKWARD STOCHASTIC DIFFERENTIAL EQUATIONS IN BOUNDED DOMAINS [J]. Journal of Computational Mathematics, 2018, 36(2): 237-258. [4] Xiao Li, Zhonghua Qiao, Hui Zhang. A SECOND-ORDER CONVEX SPLITTING SCHEME FOR A CAHN-HILLIARD EQUATION WITH VARIABLE INTERFACIAL PARAMETERS [J]. Journal of Computational Mathematics, 2017, 35(6): 693-710. [5] Qiujin Peng, Zhonghua Qiao, Shuyu Sun. STABILITY AND CONVERGENCE ANALYSIS OF SECOND-ORDER SCHEMES FOR A DIFFUSE INTERFACE MODEL WITH PENG-ROBINSON EQUATION OF STATE [J]. Journal of Computational Mathematics, 2017, 35(6): 737-765. [6] Jiwei Zhang, Zhenli Xu, Xiaonan Wu, Desheng Wang. STABILITY ANALYSIS FOR NONLINEAR SCHRÖDINGER EQUATIONS WITH NONLINEAR ABSORBING BOUNDARY CONDITIONS [J]. Journal of Computational Mathematics, 2017, 35(1): 1-18. [7] Qian Zhang, Ran Zhang. A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR SECOND-ORDER ELLIPTIC EQUATIONS WITH ROBIN BOUNDARY CONDITIONS [J]. Journal of Computational Mathematics, 2016, 34(5): 532-548. [8] Tao Tang, Jiang Yang. IMPLICIT-EXPLICIT SCHEME FOR THE ALLEN-CAHN EQUATION PRESERVES THE MAXIMUM PRINCIPLE [J]. Journal of Computational Mathematics, 2016, 34(5): 451-461. [9] Ruihan Guo, Liangyue Ji, Yan Xu. HIGH ORDER LOCAL DISCONTINUOUS GALERKIN METHODS FOR THE ALLEN-CAHN EQUATION: ANALYSIS AND SIMULATION [J]. Journal of Computational Mathematics, 2016, 34(2): 135-158. [10] Kwangil Kim, Yonghai Li. CONVERGENCE OF FINITE VOLUME SCHEMES FOR HAMILTON-JACOBI EQUATIONS WITH DIRICHLET BOUNDARY CONDITIONS [J]. Journal of Computational Mathematics, 2015, 33(3): 227-247. [11] V. Thomée, A.S. Vasudeva Murthy. FINITE DIFFERENCE METHODS FOR THE HEAT EQUATION WITH A NONLOCAL BOUNDARY CONDITION [J]. Journal of Computational Mathematics, 2015, 33(1): 17-32. [12] Francisco Guillén-González, Giordano Tier. SPLITTING SCHEMES FOR A NAVIER-STOKES-CAHN-HILLIARD MODEL FOR TWO FLUIDS WITH DIFFERENT DENSITIES [J]. Journal of Computational Mathematics, 2014, 32(6): 643-664. [13] Malte Braack, Piotr Boguslaw Mucha. DIRECTIONAL DO-NOTHING CONDITION FOR THE NAVIER-STOKES EQUATIONS [J]. Journal of Computational Mathematics, 2014, 32(5): 507-521. [14] Benyu Guo, Tao Sun, Chao Zhang. SPECTRAL AND SPECTRAL ELEMENT METHODS FOR HIGH ORDER PROBLEMS WITH MIXED BOUNDARY CONDITIONS [J]. Journal of Computational Mathematics, 2014, 32(4): 392-411. [15] Wensheng Zhang, Eric T. Chung, Chaowei Wang. STABILITY FOR IMPOSING ABSORBING BOUNDARY CONDITIONS IN THE FINITE ELEMENT SIMULATION OF ACOUSTIC WAVE PROPAGATION [J]. Journal of Computational Mathematics, 2014, 32(1): 1-20.
Viewed
Full text

Abstract