Sparse Discrete Ordinates Method in Radiative Transfer

Grella, Konstantin; Schwab, Christoph
July 2011
Computational Methods in Applied Mathematics;2011, Vol. 11 Issue 3, p305
Academic Journal
The stationary monochromatic radiative transfer equation (RTE) is a partial differential transport equation stated on a five-dimensional phase space, the Cartesian product of physical and angular domain. We solve the RTE with a Galerkin FEM in physical space and collocation in angle, corresponding to a discrete ordinates method (DOM). To reduce the complexity of the problem and to avoid the "curse of dimension", we adapt the sparse grid combination technique to the solution space of the RTE and show that we obtain a sparse DOM which uses essentially only as many degrees of freedom as required for a purely spatial transport problem. For smooth solutions, the convergence rates deteriorate only by a logarithmic factor. We compare the sparse DOM to the standard full DOM and a sparse tensor product approach developed earlier with Galerkin FEM in physical space and a spectral method in angle. Numerical experiments confirm our findings.


Related Articles

  • Adaptive Image Colorization with Total Variation and Partial Differential Equation. Ying Chen; Shensu Li; Heng Yao // International Journal of Digital Content Technology & its Applic;Dec2011, Vol. 5 Issue 12, p316 

    An adaptive image colorization algorithm based on total variation and partial differential equations (PDE) is proposed to overcome problems of color blurring near edges in some colorization algorithms. A partial differential equation is established through minimization of a functional. An...

  • Nag Fortran Library Now Has Over 1,500 Routines.  // Database & Network Journal;Feb2005, Vol. 35 Issue 1, p23 

    The version of the NAG Numerical Library includes over 300 functions, taking the total number of routines in Mark 21 to 1500, which covers the world's most widely used numerical and statistical algorithms. The functions include a complete chapter covering mesh generation that incorporates...

  • Finite elements on evolving surfaces. Dziuk, G.; Elliott, C. M. // IMA Journal of Numerical Analysis;Apr2007, Vol. 27 Issue 2, p262 

    In this article, we define a new evolving surface finite-element method for numerically approximating partial differential equations on hypersurfaces G(t) in n+1 which evolve with time. The key idea is based on approximating G(t) by an evolving interpolated polyhedral (polygonal if n = 1)...

  • Semidiscrete Schemes of the Finite Element Method for Degenerate Hyperbolic Equations. Lyashko, A.; Fedotov, E. // Differential Equations;Jul2005, Vol. 41 Issue 7, p997 

    Considers semidiscrete schemes of the finite element method for degenerate hyperbolic equations. Analysis of the Dirichlet boundary conditions; Recommendation of a method for singling out the singularity in the scheme.

  • A C�-Conforming Finite Element Method for Nonlinear Fourth-Order hyperbolic Equation. Yang Liu; Hong Li; Siriguleng He; Wei Gao; Zhichao Fang // International Journal of Engineering & Natural Sciences;2011, Vol. 5 Issue 4, p234 

    In this paper, the C�-conforming finite element method is analyzed for a class of nonlinear fourth-order hyperbolic partial differential equation. Some a priori bounds are derived using Lyapunov functional, and existence, uniqueness and regularity for the weak solutions are proved. Optimal...

  • Computable error bounds for nonconforming Fortin–Soulie finite element approximation of the Stokes problem. Ainsworth, Mark; Allendes, Alejandro; Barrenechea, Gabriel R.; Rankin, Richard // IMA Journal of Numerical Analysis;Apr2012, Vol. 32 Issue 2, p417 

    We propose computable a posteriori error estimates for a second-order nonconforming finite element approximation of the Stokes problem. The estimator is completely free of unknown constants and gives a guaranteed numerical upper bound on the error in terms of a lower bound for the inf-sup...

  • A new local stabilized nonconforming finite element method for the Stokes equations. Li, Jian; Chen, Zhangxin // Computing;Mar2008, Vol. 82 Issue 2/3, p157 

    In this paper, we propose and study a new local stabilized nonconforming finite method based on two local Gauss integrations for the two-dimensional Stokes equations. The nonconforming method uses the lowest equal-order pair of mixed finite elements (i.e., NCP 1– P 1). After a stability...

  • 2D & 3D Finite Element Method Packages of CEMTool for Engineering PDE Problems. Choon Ki Ahn; Soohee Han; Wook Hyun Kwon // Enformatika;2006, Vol. 12, p109 

    CEMTool is a command style design and analyzing package for scientific and technological algorithm and a matrix based computation language. In this paper, we present new 2D & 3D finite element method (FEM) packages for CEMTool. We discuss the detailed structures and the important features of...

  • STIFFNESS OF THE LINEAR DIFFUSION AND WAVE-TYPE PARTIAL DIFFERENTIAL EQUATIONS. Celaya, Elisabete Alberdi; Anza Aguirrezabala, Juan José // ROMAI Journal;2013, Vol. 9 Issue 2, p11 

    After applying the Finite Element Method (FEM) to the one-dimensional diffusion-type and wave-type Partial Differential Equations (PDEs) with boundary conditions and initial conditions, a first order and a second order ODE systems are obtained respectively. The latter can be reduced to a first...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics