Approximation of 3D convection diffusion equation using DQM based on modified cubic trigonometric B-splines

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.

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.


2020 ◽  
Vol 28 (1) ◽  
pp. 71-92
Author(s):  
Mourad Bellassoued ◽  
Imen Rassas

AbstractWe consider the inverse boundary value problem for the dynamical steady-state convection-diffusion equation. We prove that the first-order coefficient and the scalar potential are uniquely determined by the Dirichlet-to-Neumann map. More precisely, we show in dimension {n\geq 3} a log-type stability estimate for the inverse problem under consideration. The method is based on reducing our problem to an auxiliary inverse problem and the construction of complex geometrical optics solutions of this problem.


2014 ◽  
Vol 125 (2) ◽  
pp. 548-550 ◽  
Author(s):  
H. Caglar ◽  
N. Caglar ◽  
M. Ozer

2021 ◽  
Vol 0 (0) ◽  
pp. 0
Author(s):  
Soumen Senapati ◽  
Manmohan Vashisth

<p style='text-indent:20px;'>In this article, we study the stability in the inverse problem of determining the time-dependent convection term and density coefficient appearing in the convection-diffusion equation, from partial boundary measurements. For dimension <inline-formula><tex-math id="M1">\begin{document}$ n\ge 2 $\end{document}</tex-math></inline-formula>, we show the convection term (modulo the gauge term) admits log-log stability, whereas log-log-log stability estimate is obtained for the density coefficient.</p>


2021 ◽  
Vol 6 (5) ◽  
pp. 4370-4393
Author(s):  
Asifa Tassaddiq ◽  
◽  
Muhammad Yaseen ◽  
Aatika Yousaf ◽  
Rekha Srivastava ◽  
...  

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.


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.


Sign in / Sign up

Export Citation Format

Share Document