Vladimirov method

One of the most accurate numerical methods for solving the kinetic equation of neutron transfer in nuclear reactors, based on integration along characteristics. Proposed in 1952 by V.S. Vladimirov for the solution of integro-differential kinetic equations in the case of spherically-symmetric reactors. The principle of the method may be illustrated by taking a subcritical reactor with a neutron source as an example. For a one-dimensional spherically-symmetric geometry and a single velocity, the equation of the neutron flux (where is the radius, , and is the cosine of the angle between the vector of the neutron velocity and the radius) has the form

(1)

with the boundary condition

(2)

meaning that neutrons do not impinge from outside on the external boundary of the system; , and are given piecewise-continuous functions of . The substitution

(3)

yields the equation

(4)

where . This equation may be readily solved as an ordinary first-order differential equation, and

(5)

For each characteristic of the differential part of the kinetic equation (1) one chooses a specific system of nodes , where is the grid radius chosen. The solution of equation (5) is conducted by the method of successive approximation (cf. Sequential approximation, method of), beginning from the initial approximation of the function:

(6)

Here

(where ) is readily found with the aid of (5) at all nodes of the grid, after having replaced the integrals in (5) by sums and having obtained an expression interconnecting the values of and at two neighbouring points on the characteristic. In order to obtain the value of in the following approximation, one must compute , which is done using the quadrature formula involving the points of the circle . The rate of convergence of the successive approximations is determined by the size and the physical characteristics of the reactor.

The eigenvalue problem (the determination of the critical parameters of the reactor) is solved in a similar manner.

Vladimirov's method can be generalized to include problems involving several velocities and multi-dimensional problems, and can be readily programmed in a computer. Unlike the Carlson method, Vinogradov's method employs -variable grids for different 's, thus enhancing the computational accuracy on the reactor-vacuum boundary (close to ), as compared with the regions in the vicinity of , where the neutron flux is near-isotropic.

References

[1]	G.I. Marchuk, "Methods of calculation of nuclear reactors" , Moscow (1961) (In Russian)