Coupling of transport and diffusion models in linear transport theory

Guillaume Bal; Yvon Maday

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique (2002)

  • Volume: 36, Issue: 1, page 69-86
  • ISSN: 0764-583X

Abstract

top
This paper is concerned with the coupling of two models for the propagation of particles in scattering media. The first model is a linear transport equation of Boltzmann type posed in the phase space (position and velocity). It accurately describes the physics but is very expensive to solve. The second model is a diffusion equation posed in the physical space. It is only valid in areas of high scattering, weak absorption, and smooth physical coefficients, but its numerical solution is much cheaper than that of transport. We are interested in the case when the domain is diffusive everywhere except in some small areas, for instance non-scattering or oscillatory inclusions. We present a natural coupling of the two models that accounts for both the diffusive and non-diffusive regions. The interface separating the models is chosen so that the diffusive regime holds in its vicinity to avoid the calculation of boundary or interface layers. The coupled problem is analyzed theoretically and numerically. To simplify the presentation, the transport equation is written in the even parity form. Applications include, for instance, the treatment of clear or spatially inhomogeneous regions in near-infra-red spectroscopy, which is increasingly being used in medical imaging for monitoring certain properties of human tissues.

How to cite

top

Bal, Guillaume, and Maday, Yvon. "Coupling of transport and diffusion models in linear transport theory." ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique 36.1 (2002): 69-86. <http://eudml.org/doc/244754>.

@article{Bal2002,
abstract = {This paper is concerned with the coupling of two models for the propagation of particles in scattering media. The first model is a linear transport equation of Boltzmann type posed in the phase space (position and velocity). It accurately describes the physics but is very expensive to solve. The second model is a diffusion equation posed in the physical space. It is only valid in areas of high scattering, weak absorption, and smooth physical coefficients, but its numerical solution is much cheaper than that of transport. We are interested in the case when the domain is diffusive everywhere except in some small areas, for instance non-scattering or oscillatory inclusions. We present a natural coupling of the two models that accounts for both the diffusive and non-diffusive regions. The interface separating the models is chosen so that the diffusive regime holds in its vicinity to avoid the calculation of boundary or interface layers. The coupled problem is analyzed theoretically and numerically. To simplify the presentation, the transport equation is written in the even parity form. Applications include, for instance, the treatment of clear or spatially inhomogeneous regions in near-infra-red spectroscopy, which is increasingly being used in medical imaging for monitoring certain properties of human tissues.},
author = {Bal, Guillaume, Maday, Yvon},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique},
keywords = {linear transport; even parity formulation; diffusion approximation; domain decomposition; diffuse tomography; Boltzmann equation; diffusion equation},
language = {eng},
number = {1},
pages = {69-86},
publisher = {EDP-Sciences},
title = {Coupling of transport and diffusion models in linear transport theory},
url = {http://eudml.org/doc/244754},
volume = {36},
year = {2002},
}

TY - JOUR
AU - Bal, Guillaume
AU - Maday, Yvon
TI - Coupling of transport and diffusion models in linear transport theory
JO - ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique
PY - 2002
PB - EDP-Sciences
VL - 36
IS - 1
SP - 69
EP - 86
AB - This paper is concerned with the coupling of two models for the propagation of particles in scattering media. The first model is a linear transport equation of Boltzmann type posed in the phase space (position and velocity). It accurately describes the physics but is very expensive to solve. The second model is a diffusion equation posed in the physical space. It is only valid in areas of high scattering, weak absorption, and smooth physical coefficients, but its numerical solution is much cheaper than that of transport. We are interested in the case when the domain is diffusive everywhere except in some small areas, for instance non-scattering or oscillatory inclusions. We present a natural coupling of the two models that accounts for both the diffusive and non-diffusive regions. The interface separating the models is chosen so that the diffusive regime holds in its vicinity to avoid the calculation of boundary or interface layers. The coupled problem is analyzed theoretically and numerically. To simplify the presentation, the transport equation is written in the even parity form. Applications include, for instance, the treatment of clear or spatially inhomogeneous regions in near-infra-red spectroscopy, which is increasingly being used in medical imaging for monitoring certain properties of human tissues.
LA - eng
KW - linear transport; even parity formulation; diffusion approximation; domain decomposition; diffuse tomography; Boltzmann equation; diffusion equation
UR - http://eudml.org/doc/244754
ER -

References

top
  1. [1] M.L. Adams and E.W. Larsen, Fast iterative methods for deterministic particle transport computations. Preprint (2001). 
  2. [2] R.E. Alcouffe, Diffusion synthetic acceleration methods for the diamond-differenced discrete-ordinates equations. Nucl. Sci. Eng. 64 (1977) 344. 
  3. [3] G. Allaire and G. Bal, Homogenization of the criticality spectral equation in neutron transport. ESAIM: M2AN 33 (1999) 721–746. Zbl0931.35010
  4. [4] S.R. Arridge, H. Dehghani, M. Schweiger, and E. Okada, The finite element model for the propagation of light in scattering media: A direct method for domains with nonscattering regions. Med. Phys. 27 (2000) 252–264. 
  5. [5] G. Bal, Couplage d’équations et homogénéisation en transport neutronique. Thèse de Doctorat de l’Université Paris 6 (1997). In French. 
  6. [6] G. Bal, First-order corrector for the homogenization of the criticality eigenvalue problem in the even parity formulation of the neutron transport. SIAM J. Math. Anal. 30 (1999) 1208–1240. Zbl0937.35007
  7. [7] G. Bal,Spatially varying discrete ordinates methods in X Y - geometry. M 3 AS (Math. Models Methods Appl. Sci.) 10 (2000) 1277–1303. Zbl0969.65126
  8. [8] G. Bal, Transport through diffusive and non-diffusive regions, embedded objects, and clear layers. To appear in SIAM J. Appl. Math. Zbl1020.45004MR1918572
  9. [9] G. Bal, V. Freilikher, G. Papanicolaou, and L. Ryzhik, Wave transport along surfaces with random impedance. Phys. Rev. B 6 (2000) 6228–6240. 
  10. [10] G. Bal and L. Ryzhik, Diffusion approximation of radiative transfer problems with interfaces. SIAM J. Appl. Math. 60 (2000) 1887–1912. Zbl0976.45008
  11. [11] A. Bensoussan, J.-L. Lions, and G. Papanicolaou, Boundary layers and homogenization of transport processes. Res. Inst. Math. Sci. Kyoto Univ. 15 (1979)53–157. Zbl0408.60100
  12. [12] J.-F. Bourgat, P. Le Tallec, B. Perthame, and Y. Qiu, Coupling Boltzmann and Euler equations without overlapping, in Domain Decomposition Methods in Science and Engineering, The Sixth International Conference on Domain Decomposition, Como, Italy, June 15-19, 1992, Contemp. Math. 157, American Mathematical Society, Providence, RI (1994) 377–398. Zbl0796.76063
  13. [13] M. Cessenat, Théorèmes de trace L p pour des espaces de fonctions de la neutronique. C. R. Acad. Sci. Paris Sér. I Math. 299 (1984) 831–834. Zbl0568.46030
  14. [14] S. Chandrasekhar, Radiative Transfer. Dover Publications, New York (1960). MR111583
  15. [15] P.G. Ciarlet, The Finite Element Method for Elliptic Problems. North-Holland Publishing Company, Amsterdam, New York, Oxford (1978). Zbl0383.65058MR520174
  16. [16] R. Dautray and J.-L. Lions, Mathematical Analysis and Numerical Methods for Science and Technology 6. Springer-Verlag, Berlin (1993). Zbl0755.35001MR1295030
  17. [17] J.J. Duderstadt and W.R. Martin, Transport Theory. Wiley-Interscience, New York (1979). Zbl0407.76001MR551868
  18. [18] M. Firbank, S.A. Arridge, M. Schweiger, and D.T. Delpy, An investigation of light transport through scattering bodies with non-scattering regions. Phys. Med. Biol. 41 (1996) 767–783. 
  19. [19] F. Gastaldi, A. Quarteroni, and G. Sacchi Landriani, On the coupling of two-dimensional hyperbolic and elliptic equations: analytical and numerical approach, in Domain Decomposition Methods for Partial Differential Equations, Third International Symposium on Domain Decomposition Methods for Partial Differential Equations, Houston, TX, 1989, SIAM, Philadelphia, PA (1990) 22–63. Zbl0709.65092
  20. [20] F. Golse, P.-L. Lions, B. Perthame, and R. Sentis, Regularity of the moments of the solution of a transport equation. J. Funct. Anal. 76 (1988) 110–125. Zbl0652.47031
  21. [21] A.H. Hielscher, R.E. Alcouffe, and R.L. Barbour, Comparison of finite-difference transport and diffusion calculations for photon migration in homogeneous and heterogeneous tissues. Phys. Med. Biol. 43 (1998) 1285–1302. 
  22. [22] A. Ishimaru, Wave Propagation and Scattering in Random Media. Academics, New York (1978). Zbl0873.65115
  23. [23] B. Lapeyre, E. Pardoux and R. Sentis, Méthodes de Monte-Carlo pour les équations de transport et de diffusion, in Mathématiques & Applications 29, Springer-Verlag, Berlin (1998). Zbl0886.65124MR1621249
  24. [24] C.D. Levermore, W.J. Morokoff and B.T. Nadiga, Moment realizability and the validity of the Navier-Stokes equations for rarefied gas dynamics. Phys. Fluids 10 (1998) 3214–3226. Zbl1185.76858
  25. [25] E.E. Lewis and W.F. Miller Jr., Computational Methods of Neutron Transport. John Wiley & Sons, New York (1984). Zbl0594.65096
  26. [26] L.D. Marini and A. Quarteroni, A Relaxation procedure for domain decomposition methods using finite elements. Numer. Math. 55 (1989) 575–598. Zbl0661.65111
  27. [27] J. Planchard, Méthodes mathématiques en neutronique, in Collection de la Direction des Études et Recherches d’EDF, Eyrolles (1995). In French. 
  28. [28] L. Ryzhik, G. Papanicolaou, and J.B. Keller, Transport equations for elastic and other waves in random media. Wave Motion 24 (1996) 327–370. Zbl0954.74533
  29. [29] H. Sato and M.C. Fehler, Seismic Wave Propagation and Scattering in the Heterogeneous Earth, in AIP Series in Modern Acoustics and Signal Processing, AIP Press, Springer, New York (1998). Zbl0894.73001MR1488700
  30. [30] M. Tidriri, Asymptotic analysis of a coupled system of kinetic equations. C. R. Acad. Sci. Paris Sér. I 328 (1999) 637–642. Zbl0945.76074
  31. [31] S. Tiwari, Application of moment realizability criteria for the coupling of the Boltzmann and Euler equations. Transport Theory Statist. Phys. 29 (2000) 759–783. Zbl1174.82330

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.