Parallel characteristic finite element method for time-dependent convection-diffusion problem

2010 ◽  
Vol 18 (4) ◽  
pp. 695-705 ◽  
Author(s):  
Jiansong Zhang ◽  
Danping Yang ◽  
Hongfei Fu ◽  
Hui Guo
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.


2001 ◽  
Vol 11 (02) ◽  
pp. 301-337 ◽  
Author(s):  
K. GERDES ◽  
J. M. MELENK ◽  
C. SCHWAB ◽  
D. SCHÖTZAU

The Streamline Diffusion Finite Element Method (SDFEM) for a two-dimensional convection–diffusion problem is analyzed in the context of the hp-version of the Finite Element Method (FEM). It is proved that the appropriate choice of the SDFEM parameters leads to stable methods on the class of "boundary layer meshes", which may contain anisotropic needle elements of arbitrarily high aspect ratio. Consistency results show that the use of such meshes can resolve layer components present in the solutions at robust exponential rates of convergence. We confirm these theoretical results in a series of numerical examples.


Sign in / Sign up

Export Citation Format

Share Document