When solving integro-differential equations by means of numerical methods one has to deal with large systems of linear equations, such as happens in transport theory [10]. Many iterative techniques are now used in Transport Theory in order to solve problems of 2D and 3D dimensions. In this paper, we choose two problems to solve the following transport equation,
μψx(x,μ)+h(x,μ)ψ(x,μ)==11k(x,μ,μ')ψ(x,μ')dμ'+q(x,μ)(1)
where x: represents the spatial variable, μ: the cosine of the angle, ψ: the angular flux, h(x, μ): is the collision frequency, k(x, μ, μ'): the scattering kernel, q(x, μ): the source. The aim of this work is the straightforward application of the Krylov spaces technique [2] to the governing equation or to its discretizations derived of the discrete ordinates method (choosing a finite number of directions and then approximating the integral term by means of a proper sum). The equation (1) can be written in functional form as
Aψ=q(2)
with ψ in the Hilbert space L2 ([0,a] × [-1,1])., and q is the source function. The operator A˜ derived from a discrete ordinates scheme that approximates the operator
Aψ=μψx(x,μ)++h(x,μ)ψ(x,μ)11k(x,μ,μ')ψ(x,μ')dμ'(3)
generates the following subspace
K'''(A˜,q)=Span{q,Aq˜,A˜2q,,A˜m1q}(4)
i.e. the subspace generated by the iterations of order 0, 1, 2, ... , m−1 of the source function q. Two methods are specially outstanding, the Lanczos method to solve the problem given by equation (2) with certain boundary conditions, and the conjugate gradient method to solve the same problem with identical boundary conditions. We discuss and accelerate the basic iterative method [8]. An important conclusion is the generation of these methods to solve linear systems in Hilbert spaces, if verify the convergence conditions, which are outlined in this work. The first problem is a cubic domain with two regions, one with a source near the vertex at the origin and the shield region. In this case, the Cartesian planes (specifically 0<x<L, 0<y<L, 0<z<L) are reflexive boundaries and the rest faces of the cube are vacuum boundaries. The second problem has a geometry with three regions, one of the source, other is the void region and the third one is a shield domain.
This content is only available via PDF.
You do not currently have access to this content.