Interior penalty discontinuous Galerkin method for a homogenized diffusion equation in reactor simulations
Paper i proceeding, 2015

Full core reactor simulations are generally based on a (at least) two-scales process, the first one consisting of a fine-mesh lattice calculation, and after an homogenization process the homogenized data are used for coarse-mesh full core reactor calculation. The discontinuity factors can be considered as part of these homogenized data, widely used to minimize the error due to spatial homogenization of the cross sections. Thus, the implementation of the discontinuity factors in Finite Element Methods is necessary in order to use these methods for homogenized core calculations. Here we proposed a variation of an Interior Penalty Discontinuous Galerkin Finite Element method to allow forcing the discontinuity of the neutron flux determined by the discontinuity factors. The proposed method is tested solving different one-dimensional benchmark problems, showing that the discontinuity factors technique can be successfully introduced in the Interior Penalty Discontinuous Galerkin Finite Element Method.

Discontinuity Factors

Finite Element Method

Neutron diffusion equation


Sebastian Gonzalez-Pintor

Chalmers, Teknisk fysik, Nukleär teknik

Antoni Vidal-Ferràndiz

Universitat Politecnica de Valencia (UPV)

Damian Ginestar

Universitat Politecnica de Valencia (UPV)

Christophe Demaziere

Chalmers, Teknisk fysik, Nukleär teknik

Mohammad Asadzadeh

Göteborgs universitet

Chalmers, Matematiska vetenskaper, Matematik

Gumersindo Verdú

Universitat Politecnica de Valencia (UPV)

Proc. Joint Int. Conf. Mathematics and Computation (M&C), Supercomputing in Nuclear Applications (SNA) and the Monte Carlo (MC) Method (MC2015)



Annan teknik

Annan fysik