A stabilizer free spatial weak Galerkin finite element methods for time-dependent convection-diffusion equations

Author(s):  
Ahmed Al-Taweel ◽  
Saqib Hussain ◽  
Xiaoshen Wang ◽  
Mohammed Cheichan

In this paper, we propose a stabilizer free spatial weak Galerkin (SFSWG) finite element method for solving time-dependent convection diffusion equations based on weak form Eq. (4). SFSWG method in spatial direction and Euler difference operator Eq. (37) in temporal direction are used. The main reason for using the SFSWG method is because of its simple formulation that makes this algorithm more efficient and its implementation easier. The optimal rates of convergence of 𝒪⁢(hk) and 𝒪⁢(hk+1) in a discrete H1 and L2 norms, respectively, are obtained under certain conditions if polynomial spaces (Pk⁢(K),Pk⁢(e),[Pj⁢(K)]2) are used in the SFSWG finite element method. Numerical experiments are performed to verify the effectiveness and accuracy of the SFSWG method.

Author(s):  
Sundar Namala ◽  
Rizwan Uddin

Abstract Nodal integral methods (NIM) are a class of efficient coarse mesh methods that use transverse averaging to reduce the governing partial differential equation(s) (PDE) into a set of ordinary differential equations (ODE). The standard application of NIM is restricted to domains that have boundaries parallel to one of the coordinate axes/palnes (in 2D/3D). The hybrid nodal-integral/finite-element method (NI-FEM) reported here has been developed to extend the application of NIM to arbitrary domains. NI-FEM is based on the idea that the interior region and the regions with boundaries parallel to the coordinate axes (2D) or coordinate planes (3D) can be solved using NIM, and the rest of the domain can be discretized and solved using FEM. The crux of the hybrid NI-FEM is in developing interfacial conditions at the common interfaces between the NIM regions and FEM regions. We here report the development of hybrid NI-FEM for the time-dependent convection-diffusion equation (CDE) in arbitrary domains. Resulting hybrid numerical scheme is implemented in a parallel framework in Fortran and solved using PETSc. The preliminary approach to domain decomposition is also discussed. Numerical solutions are compared with exact solutions, and the scheme is shown to be second order accurate in both space and time. The order of approximations used for the development of the scheme are also shown to be second order. The hybrid method is more efficient compared to standalone conventional numerical schemes like FEM.


Sign in / Sign up

Export Citation Format

Share Document