Jean Ragusa

2papers

2 Papers

COMP-PHJun 6, 2019
Acceleration of Radiation Transport Solves Using Artificial Neural Networks

Mauricio Tano, Jean Ragusa

Discontinuous Finite Element Methods (DFEM) have been widely used for solving $S_n$ radiation transport problems in participative and non-participative media. In the DFEM $S_n$ methodology, the transport equation is discretized into a set of algebraic equations that have to be solved for each spatial cell and angular direction, strictly preserving the following of radiation in the system. At the core of a DFEM solver a small matrix-vector system (of 8 independent equations for tri-linear DFEM in 3D hexehdral cells) has to be assembled and solved for each cell, angle, energy group, and time step. These systems are generally solved by direct Gaussian Elimination. The computational cost of the Gaussian Elimination, repeated for each phase-space cell, amounts to a large fraction to the total compute time. Here, we have designed a Machine Learning algorithm based in a shallow Artificial Neural Networks (ANNs) to replace that Gaussian Elimination step, enabling a sizeable speed up in the solution process. The key idea is to train an ANN with a large set of solutions of random one-cell transport problems and then to use the trained ANN to replace Gaussian Elimination large scale transport solvers. It has been observed that ANNs decrease the solution times by at least a factor of 4, while introducing mean absolute errors between 1-3 \% in large scale transport solutions.

NAMay 8, 2019
Positive asymptotic preserving approximation of the radiation transport equation

Jean-Luc Guermond, Bojan Popov, Jean Ragusa

We introduce a (linear) positive and asymptotic preserving method or solving the one-group radiation transport equation. The approximation in space is discretization agnostic: the space approximation can be done with continuous or discontinuous finite elements (or finite volumes, or finite differences). The method is first-order accurate in space. This type of accuracy is coherent with Godunov's theorem since the method is linear. The two key theoretical results of the paper are Theorem~4.4 and Theorem~4.8. The method is illustrated with continuous finite elements. It is observed to converge with the rate $\calO(h)$ in the $L^2$-norm on manufactured solutions, and it is $\calO(h^2)$ in the diffusion regime. Unlike other standard techniques, the proposed method does not suffer from overshoots at the interfaces of optically thin and optically thick regions.