Author:
Calloo Atyab A.,Hébert Alain
Abstract
This paper presents an implementation and a comparison of two spatial discretisation schemes over a hexagonal geometry for the two-dimensional discrete ordinates transport equation. The methods are a high-order Discontinuous Galerkin (DG) finite element scheme and a high-order Diamond Differencing (DD) scheme. The DG method has been, and is being, studied on the hexagonal geometry, also called a honeycomb mesh – but not the DD method. In this research effort, it was chosen to divide the hexagons into (at least) three lozenges. An affine transformation is then applied onto said lozenges to cast them into the reference quadrilaterals usually studied in finite elements. In practice, this effectively means that the equations used in Cartesian geometry have their terms and operators altered using the Jacobian matrix of the transformation. This was implemented in the discrete ordinates solver of the code DRAGON5. Two 2D benchmark problems were then used for the verification and validation, including one based on the Monju 3D reactor benchmark. It was found that the diamond-differencing scheme seemed better. It converged much faster towards the solution at comparable mesh refinements for first-order expansion of the flux. Even if this difference was not present for second-order, DG was slower, about two to four times slower.
Reference7 articles.
1. Calloo A. and Hébert A.. “Development and Evaluation of a Discontinuous Galerkin Method in the DRAGON5 Code.” In 38th Annual Conference of the Canadian Nuclear Society and 42nd Annual CNS/CNA Student Conference. Canadian Nuclear Society, Saskatoon, SK, Canada (2018).
2. Reed W. H. and Hill T. R.. “Triangular mesh methods for the neutron transport for a scalar hyperbolic equation.” Technical Report LA-UR-73-479, Los Alamos Scientific Lab., N. Mex. (1973).
3. High order diamond differencing schemes
Cited by
2 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献