# American Institute of Mathematical Sciences

March  2009, 4(1): 35-65. doi: 10.3934/nhm.2009.4.35

## A uniformly second order numerical method for the one-dimensional discrete-ordinate transport equation and its diffusion limit with interface

 1 Department of Mathematics, University of Wisconsin, Madison, WI 53706 2 Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China, China

Received  August 2008 Revised  October 2008 Published  February 2009

In this paper, we propose a uniformly second order numerical method for the discete-ordinate transport equation in the slab geometry in the diffusive regimes with interfaces. At the interfaces, the scattering coefficients have discontinuities, so suitable interface conditions are needed to define the unique solution. We first approximate the scattering coefficients by piecewise constants determined by their cell averages, and then obtain the analytic solution at each cell, using which to piece together the numerical solution with the neighboring cells by the interface conditions. We show that this method is asymptotic-preserving, which preserves the discrete diffusion limit with the correct interface condition. Moreover, we show that our method is quadratically convergent uniformly in the diffusive regime, even with the boundary layers. This is 1) the first sharp uniform convergence result for linear transport equations in the diffusive regime, a problem that involves both transport and diffusive scales; and 2) the first uniform convergence valid up to the boundary even if the boundary layers exist, so the boundary layer does not need to be resolved numerically. Numerical examples are presented to justify the uniform convergence.
Citation: Shi Jin, Min Tang, Houde Han. A uniformly second order numerical method for the one-dimensional discrete-ordinate transport equation and its diffusion limit with interface. Networks & Heterogeneous Media, 2009, 4 (1) : 35-65. doi: 10.3934/nhm.2009.4.35
 [1] Giada Basile, Tomasz Komorowski, Stefano Olla. Diffusion limit for a kinetic equation with a thermostatted interface. Kinetic & Related Models, 2019, 12 (5) : 1185-1196. doi: 10.3934/krm.2019045 [2] Liu Liu. Uniform spectral convergence of the stochastic Galerkin method for the linear semiconductor Boltzmann equation with random inputs and diffusive scaling. Kinetic & Related Models, 2018, 11 (5) : 1139-1156. doi: 10.3934/krm.2018044 [3] Xu Yang, François Golse, Zhongyi Huang, Shi Jin. Numerical study of a domain decomposition method for a two-scale linear transport equation. Networks & Heterogeneous Media, 2006, 1 (1) : 143-166. doi: 10.3934/nhm.2006.1.143 [4] Stefan Possanner, Claudia Negulescu. Diffusion limit of a generalized matrix Boltzmann equation for spin-polarized transport. Kinetic & Related Models, 2011, 4 (4) : 1159-1191. doi: 10.3934/krm.2011.4.1159 [5] Shi Jin, Xu Yang, Guangwei Yuan. A domain decomposition method for a two-scale transport equation with energy flux conserved at the interface. Kinetic & Related Models, 2008, 1 (1) : 65-84. doi: 10.3934/krm.2008.1.65 [6] Alina Chertock, Changhui Tan, Bokai Yan. An asymptotic preserving scheme for kinetic models with singular limit. Kinetic & Related Models, 2018, 11 (4) : 735-756. doi: 10.3934/krm.2018030 [7] Hélène Hivert. Numerical schemes for kinetic equation with diffusion limit and anomalous time scale. Kinetic & Related Models, 2018, 11 (2) : 409-439. doi: 10.3934/krm.2018019 [8] Matthieu Alfaro, Arnaud Ducrot. Sharp interface limit of the Fisher-KPP equation. Communications on Pure & Applied Analysis, 2012, 11 (1) : 1-18. doi: 10.3934/cpaa.2012.11.1 [9] Pedro Aceves-Sánchez, Christian Schmeiser. Fractional diffusion limit of a linear kinetic equation in a bounded domain. Kinetic & Related Models, 2017, 10 (3) : 541-551. doi: 10.3934/krm.2017021 [10] Karl Kunisch, Markus Müller. Uniform convergence of the POD method and applications to optimal control. Discrete & Continuous Dynamical Systems - A, 2015, 35 (9) : 4477-4501. doi: 10.3934/dcds.2015.35.4477 [11] Qiang Du, Manlin Li. On the stochastic immersed boundary method with an implicit interface formulation. Discrete & Continuous Dynamical Systems - B, 2011, 15 (2) : 373-389. doi: 10.3934/dcdsb.2011.15.373 [12] Lizhi Ruan, Changjiang Zhu. Boundary layer for nonlinear evolution equations with damping and diffusion. Discrete & Continuous Dynamical Systems - A, 2012, 32 (1) : 331-352. doi: 10.3934/dcds.2012.32.331 [13] Champike Attanayake, So-Hsiang Chou. An immersed interface method for Pennes bioheat transfer equation. Discrete & Continuous Dynamical Systems - B, 2015, 20 (2) : 323-337. doi: 10.3934/dcdsb.2015.20.323 [14] Bhargav Kumar Kakumani, Suman Kumar Tumuluri. Asymptotic behavior of the solution of a diffusion equation with nonlocal boundary conditions. Discrete & Continuous Dynamical Systems - B, 2017, 22 (2) : 407-419. doi: 10.3934/dcdsb.2017019 [15] Shu Wang, Chundi Liu. Boundary Layer Problem and Quasineutral Limit of Compressible Euler-Poisson System. Communications on Pure & Applied Analysis, 2017, 16 (6) : 2177-2199. doi: 10.3934/cpaa.2017108 [16] So-Hsiang Chou. An immersed linear finite element method with interface flux capturing recovery. Discrete & Continuous Dynamical Systems - B, 2012, 17 (7) : 2343-2357. doi: 10.3934/dcdsb.2012.17.2343 [17] Masahiro Suzuki. Asymptotic stability of a boundary layer to the Euler--Poisson equations for a multicomponent plasma. Kinetic & Related Models, 2016, 9 (3) : 587-603. doi: 10.3934/krm.2016008 [18] Yangang Chen, Justin W. L. Wan. Numerical method for image registration model based on optimal mass transport. Inverse Problems & Imaging, 2018, 12 (2) : 401-432. doi: 10.3934/ipi.2018018 [19] Raffaele D’Ambrosio, Giuseppe De Martino, Beatrice Paternoster. A symmetric nearly preserving general linear method for Hamiltonian problems. Conference Publications, 2015, 2015 (special) : 330-339. doi: 10.3934/proc.2015.0330 [20] Yves Bourgault, Damien Broizat, Pierre-Emmanuel Jabin. Convergence rate for the method of moments with linear closure relations. Kinetic & Related Models, 2015, 8 (1) : 1-27. doi: 10.3934/krm.2015.8.1

2018 Impact Factor: 0.871