scholarly journals Space–Time Radial Basis Function–Based Meshless Approach for Solving Convection–Diffusion Equations

Mathematics ◽  
2020 ◽  
Vol 8 (10) ◽  
pp. 1735
Author(s):  
Cheng-Yu Ku ◽  
Jing-En Xiao ◽  
Chih-Yu Liu

This article proposes a space–time meshless approach based on the transient radial polynomial series function (TRPSF) for solving convection–diffusion equations. We adopted the TRPSF as the basis function for the spatial and temporal discretization of the convection–diffusion equation. The TRPSF is constructed in the space–time domain, which is a combination of n–dimensional Euclidean space and time into an n + 1–dimensional manifold. Because the initial and boundary conditions were applied on the space–time domain boundaries, we converted the transient problem into an inverse boundary value problem. Additionally, all partial derivatives of the proposed TRPSF are a series of continuous functions, which are nonsingular and smooth. Solutions were approximated by solving the system of simultaneous equations formulated from the boundary, source, and internal collocation points. Numerical examples including stationary and nonstationary convection–diffusion problems were employed. The numerical solutions revealed that the proposed space–time meshless approach may achieve more accurate numerical solutions than those obtained by using the conventional radial basis function (RBF) with the time-marching scheme. Furthermore, the numerical examples indicated that the TRPSF is more stable and accurate than other RBFs for solving the convection–diffusion equation.

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.


Author(s):  
Mohammad Tamsir ◽  
Neeraj Dhiman ◽  
F.S. Gill ◽  
Robin

This paper presents an approximate solution of 3D convection diffusion equation (CDE) using DQM based on modified cubic trigonometric B-spline (CTB) basis functions. The DQM based on CTB basis functions are used to integrate the derivatives of space variables which transformed the CDE into the system of first order ODEs. The resultant system of ODEs is solved using SSPRK (5,4) method. The solutions are approximated numerically and also presented graphically. The accuracy and efficiency of the method is validated by comparing the solutions with existing numerical solutions. The stability analysis of the method is also carried out.


2013 ◽  
Vol 11 (01) ◽  
pp. 1350053 ◽  
Author(s):  
NACHIKETA MISHRA ◽  
Y. V. S. S. SANYASIRAJU

Exponential compact higher-order schemes have been developed for unsteady convection-diffusion equation (CDE). One of the developed scheme is sixth-order accurate which is conditionally stable for the Péclet number 0 ≤ Pe ≤ 2.8 and the other is fourth-order accurate which is unconditionally stable. Schemes for two-dimensional (2D) problems are made to use alternate direction implicit (ADI) algorithm. Example problems are solved and the numerical solutions are compared with the analytical solutions for each case.


2011 ◽  
Vol 21 (10) ◽  
pp. 2049-2068 ◽  
Author(s):  
ERIK BURMAN ◽  
GILES SMITH

We consider space semi-discretization of the transient convection–diffusion equation using a Streamline Upwind Petrov–Galerkin (SUPG) finite element method. We show stability and convergence both for the standard SUPG formulation with elementwise evaluated Laplacian and for a weakly consistent formulation where the elementwise Laplacian is replaced by a discrete, reconstructed, Laplacian. Some numerical examples are presented to illustrate the theory.


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.


Sign in / Sign up

Export Citation Format

Share Document