Numerical Solutions of Steady Tensorial Transport Equations Using Discontinuous Galerkin Method Implemented in FreeFem++

Abstract

The main purpose of this paper is to approximate the solution of the steady tensorial transport equations using discontinuous Galerkin finite element method implemented with the finite element solver FreeFem++. After introducing the formulations of the tensorial transport equations, the analysis of its componentwise equations, i.e., advection-reaction equations have been discussed. Discretizing the transport problem using discontinuous Galerkin finite element method, the iterative fixed-point method is used to obtain the solutions. We present the numerical simulations of two-dimensional benchmark problem and observe the instability of elasticity. All the simulations are done using the script developed in FreeFem++.