A composite nodal finite element for hexagons

J. P. Hennart, E. H. Mund, E. Del Valle

Research output: Contribution to journalArticlepeer-review

9 Scopus citations

Abstract

A nodal algorithm for the solution of the multigroup diffusion equations in hexagonal arrays is analyzed. Basically, the method consists of dividing each hexagon into four quarters and mapping the hexagon quarters onto squares. The resulting boundary value problem on a quadrangular domain is solved in primal weak formulation. Nodal finite element methods like the Raviart-Thomas RTk schemes provide accurate analytical expansions of the solution in the hexagons. Transverse integration cannot be performed on the equations in the quadrangular domain as simply as it is usually done on squares because these equations have essentially variable coefficients. However, by considering an auxiliary problem with constant coefficients (on the same quadrangular domain) and by using a "preconditioning" approach, transverse integration can be performed as for rectangular geometry. A description of the algorithm is given for a one-group diffusion equation. Numerical results are presented for a simple model problem with a known analytical solution and for keff evaluations of some benchmark problems proposed in the literature. For the analytical problem, the results indicate that the theoretical convergence orders of RTk schemes (k = 0,1) are obtained, yielding accurate solutions at the expense of a few preconditioning iterations.

Original languageEnglish
Pages (from-to)139-153
Number of pages15
JournalNuclear Science and Engineering
Volume127
Issue number2
DOIs
StatePublished - Oct 1997

Fingerprint

Dive into the research topics of 'A composite nodal finite element for hexagons'. Together they form a unique fingerprint.

Cite this