## Abstract In this paper we present the asymptotic analysis of the linear Boltzmann equation for neutrons with a small positive parameter ϵ related to the mean free path, based upon the Chapman–Enskog procedure of the kinetic theory. We prove that if proper initial conditions derived by considerin
Parallel solvers for the transient multigroup neutron diffusion equations
✍ Scribed by R. Scheichl
- Publisher
- John Wiley and Sons
- Year
- 2000
- Tongue
- English
- Weight
- 198 KB
- Volume
- 47
- Category
- Article
- ISSN
- 0029-5981
No coin nor oath required. For personal study only.
✦ Synopsis
For the safety and the control of a nuclear power plant it is necessary to simulate the constituent processes on a computer system. The three-dimensional multigroup neutron di usion equations are commonly used to describe the nuclear ÿssion in the reactor core. They form a complicated system of coupled parabolic partial di erential equations (PDEs) whose solution can involve very intensive computing. In this paper this system of PDEs is discretized using a special cell-centred mixed ÿnite volume method (NEM-M0) in space, and a method that combines Crank-Nicholson and the BDF(2)-method in time. The linear equation systems which arise are solved with multi-grid as well as with preconditioned BiCGStab. The kernel of both solution methods is an e ective Block-SOR method that makes use of the particular structure of the linear equation systems. The parallelization strategy is based on a grid partitioning that distributes the data and the work homogeneously on the processors. Finally, the program was tested for three typical reactor simulation problems on grids with di ering coarseness. The speedup achieved by parallelizing multi-grid and preconditioned Bi-CGStab was outstanding for all examples; even superlinear in some cases. Moreover, the parallel execution times were better than the parallel execution times of other established reactor simulation codes.
📜 SIMILAR VOLUMES
A Petrov4alerkin formulation based on two different perturbations to the weighting functions is presented. These perturbations stabilize the oscillations that are normally exhibited by the numerical solution of the transient advective-diffusive equation in the vicinity of sharp gradients produced by
Elliptic PDEs with variable coefficients in a domain with complex geometry occur in many ocean models. The parallelization of the elliptic solver by the Shur complement method is presented for the ice-ocean model BRIOS. The Schur complement method is usually employed as an iterative solver, but for
We present a specific use of domain decomposition and decomposition in function space combined with asymptotic analytical qualitative results to obtain, on parallel computers, efficient and accurate solvers [3] for rapidly varying quasi-planar unsteady combustion fronts in liquids. In particular, we