Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 84
PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY Edited by: B.H.V. Topping, G. Montero and R. Montenegro
Paper 216
Sequential Numerical Study of some Variations of a SecondDegree Iterative Method: Application to the Neutron Diffusion Equation O. Flores^{1}^{2}, V. Vidal^{2}, G. Verdú^{3}, J. Garayoa^{2} and P. Flores^{4}
^{1}Department of Computer Science, Tuxtepec Institute of Technology, Oaxaca, Mexico
O. Flores, V. Vidal, G. Verdú, J. Garayoa, P. Flores, "Sequential Numerical Study of some Variations of a SecondDegree Iterative Method: Application to the Neutron Diffusion Equation", in B.H.V. Topping, G. Montero, R. Montenegro, (Editors), "Proceedings of the Fifth International Conference on Engineering Computational Technology", CivilComp Press, Stirlingshire, UK, Paper 216, 2006. doi:10.4203/ccp.84.216
Keywords: seconddegree iterative methods, neutron diffusion equation, largesparse linear systems, Krylov subspace methods, transient analysis, energy groups.
Summary
For design and safety reasons, nuclear power plants need fast and accurate plant simulators. The centre point of concern in the simulation of a nuclear power plant is the reactor core. Since it is the source of the energy that is produced in the reactor, a very accurate model of its processes is needed.
The neutron population into the reactor core is modeled using the Boltzmann transport equation. This threedimensional problem is modeled by a complicated system of coupled partial differential equations, the multigroup neutron diffusion equation[1,2], that have been discretised using a nodal collocation method in space and onestep backwarddifference method in time. The solution of these equations can involve very intensive computing. Therefore, it is necessary to find effective algorithms for the solution of the threedimensional model.
Bru et al. [3] applied two SecondDegree methods (A and B methods) [4] to solve the linear system of equations related to a twodimensional neutrondiffusion equation case. In our work, we test these methods to solve the linear system of equations related to a threedimensional neutrondiffusion equation for the case of the commercial reactor of Leibstadt [5]. Thus, when the neutrondiffusion equation is discretised in space and time, it is necessary to solve a large and sparse system of linear equations with the following block structure for each time step: In addition to the A and B methods to solve (116), we propose the following operations In scheme (117), named method C, we have added two different parameters and for each system to be computed, in order to accelerate its convergence. Besides to method C, we have carried out experiments with an `adaptable' technique, achieving some improvements in the process efficiency. This technique computes the solution of the systems and with a cheap precision at initial stages of the method. Then, this precision is 'adapted' or 'improved' towards a more demanding one in successive iterations. The application of this technique to method C gives the following algorithm, which we have named method D in present work, and we have obtained an acceptable performance in comparison to the others methods presented in this work.
Algorithm (Adaptable version) (1) Set ; { } (2) Set where ; (2) Solve (3) Solve (4) For (4a) Solve for with tolerance (4b) Solve for with tolerance (4c) if
end if until and References
purchase the fulltext of this paper (price £20)
go to the previous paper 
