COMPARISON BETWEEN GLOBAL, CLASSICAL DOMAIN DECOMPOSITION AND LOCAL, SINGLE AND DOUBLE COLLOCATION METHODS BASED ON RBF INTERPOLATION FOR SOLVING CONVECTION-DIFFUSION EQUATION

2008 ◽  
Vol 19 (11) ◽  
pp. 1737-1751 ◽  
Author(s):  
GAIL GUTIERREZ ◽  
WHADY FLOREZ

This work presents a performance comparison of several meshless RBF formulations for convection-diffusion equation with moderate-to-high Peclet number regimes. For the solution of convection-diffusion problems, several comparisons between global (full-domain) meshless RBF methods and mesh-based methods have been presented in the literature. However, in depth studies between new local RBF collocation methods and full-domain symmetric RBF collocation methods are not reported yet. The RBF formulations included: global symmetric method, symmetric double boundary collocation method, additive Schwarz domain decomposition method (DDM) when it is incorporated into two anterior approaches, and local single and double collocation methods. It can be found that the accuracy of solutions deteriorates as Pe increases, if no special treatment is used. From the numerical tests, it seems that the local methods, especially the derived double collocation technique incorporating PDE operator, are more effective than full domain approaches even with iterative DDM in solving moderate-to-high Pe convection-diffusion problems subject to mixed boundary conditions.

Mathematics ◽  
2019 ◽  
Vol 8 (1) ◽  
pp. 13
Author(s):  
Rui Du ◽  
Jincheng Wang ◽  
Dongke Sun

We have tested the accuracy and stability of lattice-Boltzmann (LB) simulations of the convection-diffusion equation in a two-dimensional channel flow with reactive-flux boundary conditions. We compared several different implementations of a zero-concentration boundary condition using the Two-Relaxation-Time (TRT) LB model. We found that simulations using an interpolation of the equilibrium distribution were more stable than those based on Multi-Reflection (MR) boundary conditions. We have extended the interpolation method to include mixed boundary conditions, and tested the accuracy and stability of the simulations over a range of Damköhler and Péclet numbers.


Author(s):  
Sundar Namala ◽  
Rizwan Uddin

Abstract Nodal integral methods (NIM) are a class of efficient coarse mesh method that use transverse averaging to reduce the governing partial differential equation(s) (PDE) into a set of ordinary differential equations (ODE), and these ODEs or their approximations are analytically solved. Since this method depends on transverse averaging, the standard application of this approach gets restricted to domains that have boundaries that are parallel to one of the coordinate axes (2D) or coordinate planes (3D). The hybrid nodal-integral/finite-element method (NI-FEM) 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 solved using FEM. The crux of the hybrid NI-FEM is in developing interfacial conditions at the common interfaces between the regions solved by the NIM and the FEM. Since the discrete variables in the two numerical approaches are different, this requires special treatment of the discrete quantities on the interface between the two different types of discretized elements. We here report the development of hybrid NI-FEM in a parallel framework in Fortran using PETSc for the time-dependent convection-diffusion equation (CDE) in arbitrary domains. 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 efficient compared to standalone conventional numerical schemes like FEM.


2021 ◽  
pp. 399-408
Author(s):  
Zhenwei Zhu ◽  
Junjie Chen

The convection-diffusion equation is of primary importance in understanding transport phenomena within a physical system. However, the currently available methods for solving unsteady convection-diffusion problems are generally not able to offer excellent accuracy in both time and space variables. A procedure was given in detail to solve the unsteady one-dimensional convection-diffusion equation through a combination of Runge-Kutta methods and compact difference schemes. The combination method has fourth-order accuracy in both time and space variables. Numerical experiments were conducted and the results were compared with those obtained by the Crank-Nicolson method in order to check the accuracy of the combination method. The analysis results indicated that the combination method is numerically stable at low wave numbers and small CFL numbers. The combination method has higher accuracy than the Crank-Nicolson method.


Sign in / Sign up

Export Citation Format

Share Document