Next Article in Journal
Robust Multi-Label Classification with Enhanced Global and Local Label Correlation
Next Article in Special Issue
Magnetorheological Fluid of High-Speed Unsteady Flow in a Narrow-Long Gap: An Unsteady Numerical Model and Analysis
Previous Article in Journal
Centrally Essential Rings and Semirings
Previous Article in Special Issue
Multiple Comparison Procedures for Exponential Mean Lifetimes Compared with Several Controls
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Spatio-Temporal Fully Meshless Method for Parabolic PDEs

by
Juan José Benito
1,†,
Ángel García
1,†,
Mihaela Negreanu
2,†,
Francisco Ureña
1,*,† and
Antonio M. Vargas
2,†
1
Escuela Técnica de Ingenieros Industriales, Universidad Nacional de Educación a Distancia, 28040 Madrid, Spain
2
Departamento de Análisis Matemático y Matemática Aplicada, Instituto de Matemática Interdisciplinar, Universidad Complutense de Madrid, 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2022, 10(11), 1870; https://doi.org/10.3390/math10111870
Submission received: 18 April 2022 / Revised: 25 May 2022 / Accepted: 27 May 2022 / Published: 30 May 2022
(This article belongs to the Special Issue Applications of Partial Differential Equations in Engineering)

Abstract

:
We introduce a meshless method derived by considering the time variable as a spatial variable without the need to extend further conditions to the solution of linear and non-linear parabolic PDEs. The method is based on a moving least squares method, more precisely, the generalized finite difference method (GFDM), which allows us to select well-conditioned stars. Several 2D and 3D examples, including the time variable, are shown for both regular and irregular node distributions. The results are compared with explicit GFDM both in terms of errors and execution time.

                Zum Raum wird hier die Zeit
                (Here, time becomes space)
                (Richard Wagner, Parsifal)

1. Introduction

A wide variety of physical phenomena are studied mathematically through evolution partial differential equations. The difference between parabolic evolution equations and elliptic equations is important, both from the point of view of mathematical analysis and numerical analysis. In this article, we are interested in the numerical treatment of parabolic PDEs. In particular, we focus on the numerical solution at a time T of the general parabolic PDE
U t   +   L ( x , t , U )   =   h ( x , t ) , x Γ , 0   <   t     T
with the boundary and initial condition
U ( x , t )   =   g ( x , t ) , x Γ ; U ( x , 0 )   =   U 0 ( x ) .
Here, Γ is a bounded domain of R or R 2 , and the non-linear operator
L ( x , t , U )   =   f 1 ( x , t ) U x   +   f 2 ( x , t ) 2 U x 2 , if Γ R
or, if Γ R 2 ,
L ( x , t , U )   =   f 1 ( x , y , t ) U x   +   f 2 ( x , y , t ) U y   +   f 3 ( x , y , t ) 2 U x 2 + f 4 ( x , y , t ) 2 U y 2   +   f 5 ( x , y , t ) 2 U x y .
Throughout the paper, we assume enough regularity of the functions h , g , U 0 and f i in order to provide classical solutions.
The time evolution of the modeled physical problem creates the need to study the propagation of the numerical error. This, in some cases, produces consistent numerical schemes that generate unstable solutions. The most common way to solve parabolic partial differential equations numerically is to use a different discretisation for the time and space derivatives. For instance,
U t U ( x , y , t   +   Δ t ) U ( x , y , t ) Δ t
is frequently used for the time derivative, which entails a study of the stability of the numerical method.
In this paper, we introduce a meshless method (see [1,2,3,4] for a detailed description of meshless methods) for resolving parabolic PDEs that consists of using the formulas given by the generalized finite difference method to discretise all the derivatives of the equation, including the time derivative. The researchers Jensen [5], Liszka and Orkisz [6], Orkisz [7] and Perrone and Kao [8] have contributed to the development of the GFDM, as well as Benito, Gavete and Ureña [9] who studied the influence of several factors and developed the h-adaptive method for the solution of the PDEs in 2D. During the last two decades, the method has attracted growing attention, not only because of its applications, but also, for instance, because some improvements to the one dimensional method have been made [10].
By considering the time variable as a spatial variable in a one-dimensional problem (or two-dimensional), we can find the 2D (or 3D) Taylor series of the solution, in the variables x and t (or x , y and t), and at any nodal point a desired discretization, without adding any extra conditions. In this way, we avoid the stability issues related to the time evolution of the problem. We transform problem (1) and (2) into an elliptic problem with boundary conditions (see Figure 1).
B . C . = g ( x , t ) , ( x , t ) Γ   ×   [ 0 , T ] , U 0 ( x ) , ( x , t ) Γ   ×   { t   =   0 } , U ( x , t ) , ( x , t ) Γ   ×   { t   =   T } .
Our procedure consists of defining a GFD star at each point of Γ   ×   [ 0 , T ) and Γ   ×   { t   =   T } and solving the elliptic problem without adding any extra condition or applying any other numerical method. A pictorial definition of the selection of the stars is given in Figure 2.
For this reason, we have called this approach the “space time cloud method (STCM)”. The previous discretization of the computational domain for parabolic PDEs is not new. Several meshless numerical methods have recently been applied, such as a kernel-based method [11], the space–time diffuse approximation meshless method [12], the Trefftz method [13], the localized radial basis function collocation method [14] and the multiquadratic method with the radial basis collocation method [15]. In addition, the GFDM has been applied to solve space-time problems in [16,17,18]. The novelty of the STCM method is to use only the generalized finite difference formulas to solve a time-dependent problem (parabolic in this case), which makes it unnecessary to study the stability and the calculation of the time step in the explicit case, or not to have to solve for each time step a system of equations in the implicit case. With the solution of a single system, the problem is solved. The advantage of STCM over the explicit GFDM is essentially not having to study the stability and, therefore, not having to calculate the time step and perform the calculations with that value. In addition, the advantage over FEM is that it solves the equation without the need for integration, it does not require mesh and the boundary conditions are simpler to implement, and, with respect to FDM, its main advantage is that irregular meshes can be used.
The structure of the paper is as follows. In Section 2, we introduce the meshless space time cloud method and derive the discretization of the Equations (1)–(3) for two and three dimensions. We provide several examples where we compare the present method with the standard GFDM applied to the original problem. Finally, some conclusions are drawn.

2. Space Time Cloud Method (STCM)

We introduce the STCM for the two dimensional case (that is, Γ R ). Let Γ *   =   Γ   ×   [ 0 , T ] be a domain and
P = { x 1 , , x n } Γ *
a discretization of Γ * with n points (see Figure 1). We designate each point of the discretization P as a point. Λ s   =   { x 0 ; x 0 1 , , x 0 s } P with the central node x 0 P and x 0 i , ( i = 1 , , s ) P is a set of points located near x 0 . The selection of the points for each star can be made according to criteria as four quadrants or distance (see [9]).
With the previous notation, x 0   =   ( x 0 , t 0 ) and x 0 i   =   ( x i , t i ) , define h i   =   x i     x 0 , k i   =   t i     t 0 . Let us put U 0   =   U ( x 0 ) and U i   =   U ( x 0 i ) , then, by the Taylor series expansion, we have
U i   =   U 0   +   h i U 0 x   +   k i U 0 t + 1 2 h i 2 2 U 0 x 2   +   k i 2 2 U 0 t 2   +   2 h i k i 2 U 0 x t   +   R 2 ,
for i   =   1 , . . . , s , where R 2 is the remainder of the second order. Let us use the notations
υ i T   =   { h i , k i , h i 2 2 , k i 2 2 , h i k i }
and
Π 5 T = { u 0 x , u 0 t , 2 u 0 x 2 , 2 u 0 t 2 , 2 u 0 x t } .
By truncation in (4) up to the second order, we can obtain a second-order approximation of U i , which we shall denote u i . Then, we define the weighted residual function:
T ( u ) = i   =   1 s [ ( u 0     u i )   +   h i u 0 x   +   k i u 0 t   +   + 1 2 ( h i 2 2 u 0 x 2   +   k i 2 2 u 0 t 2   +   2 h i k i 2 u 0 x t ) ] 2 η i 2 ,
where η i   =   η ( h i , k i ) are positive weighting functions as [19]. We minimize the norm given by (5) with respect to the partial derivatives by considering the following linear system
W Π 5   =   β
where
W = h 1 h 2 h s k 1 k 2 k s h 1 k 1 h 2 k 2 h s k s η 1 2 η 2 2 η s 2 h 1 k 1 h 1 k 1 h 2 k 2 h 2 k 2 h s k s h s k s ,
and
β T   =   i = 1 s ( u 0   +   u i ) h i η i 2 , i = 1 s ( u 0   +   u i ) k i η i 2 , i = 1 s ( u 0   +   u i ) h i 2 η i 2 2 , i   =   1 s ( u 0   +   u i ) k i 2 η i 2 2 , i   =   1 s ( u 0   +   u i ) h i k i η i 2 .
In [20], the authors proved that W is a positive definite matrix and that the local truncation error is O ( h i 2 , k i 2 ) . Then, W admits a unique Cholesty decomposition, and, if we define
W 1   =   K K T ,
we have
Π 5   =   K K T β .
Thus, Equation (6) can be rewritten as
Π 5   =   u 0 K K T i   =   1 s η i 2 υ i   +   K K T i   =   1 s u i η i 2 υ i .
We denote the derivatives as follows
u ( x 0 , t 0 ) x   =   φ 01 u 0   +   i   =   1 s φ i 1 u i   +   O ( h i 2 , k i 2 ) , w i t h φ 01   =   i   =   1 s φ i 1 , u ( x 0 , t 0 ) t   =   φ 02 u 0   +   i = 1 s φ i 2 u i + O ( h i 2 , k i 2 ) , w i t h φ 02   =   i = 1 s φ i 2 , 2 u ( x 0 , t 0 ) x 2   =   φ 03 u 0   +   i   =   1 s φ i 3 u i   +   O ( h i 2 , k i 2 ) , w i t h φ 03   =   i = 1 s φ i 3 , 2 u ( x 0 , t 0 ) t 2   =   φ 04 u 0   +   i   =   1 s φ i 4 u i   +   O ( h i 2 , k i 2 ) , w i t h φ 03   =   i   =   1 s φ i 4 , 2 u ( x 0 , t 0 ) x t   =   φ 05 u 0   +   i   =   1 s φ i 5 u i + O ( h i 2 , k i 2 ) ) , w i t h φ 05   =   i   =   1 s φ i 5 ,
or, written in vectorial form,
Π 5 u ( x 0 , t 0 )   =   φ 0 u 0   +   i   =   1 s φ i u i + O ( h i 2 , k i 2 )
where φ 0 and φ i stand for
φ 0   =   { φ 01 , φ 02 , φ 03 , φ 04 , φ 05 } T , φ i = { φ i 1 , φ i 2 , φ i 3 , φ i 4 , φ i 5 } T ,
fulfilling
φ 0 = i   =   1 s φ i .
Since the method allows us to choose different numbers of points per star, greater or equal to 6, and the selection of these nodes can be made by different criteria (distance, quadrant, octant), in order to obtain well-conditioned stars, the STCM allows us to solve the problem without the requirement of adding further conditions.
Similarly for the 3D case, the approximation derivatives are:
u ( x 0 , y 0 , t 0 ) x   =   φ 01 u 0   +   i   =   1 s φ i 1 u i + O ( h i 2 , k i 2 , τ i 2 ) , u ( x 0 , y 0 , t 0 ) y = φ 02 u 0   +   i   =   1 s φ i 2 u i + O ( h i 2 , k i 2 , τ i 2 ) , u ( x 0 , y 0 , t 0 ) t   =   φ 03 u 0 + i   =   1 s φ i 3 u i   +   O ( h i 2 , k i 2 , τ i 2 ) , 2 u ( x 0 , y 0 , t 0 ) x 2   =   φ 04 u 0 + i   =   1 s φ i 4 u i + O ( h i 2 , k i 2 , τ i 2 ) , 2 u ( x 0 , y 0 , t 0 ) y 2   =   φ 05 u 0 + i   =   1 s φ i 5 u i + O ( h i 2 , k i 2 , τ i 2 ) , 2 u ( x 0 , y 0 , t 0 ) t 2   =   φ 06 u 0   +   i = 1 s φ i 6 u i   +   O ( h i 2 , k i 2 , τ i 2 ) , 2 u ( x 0 , y 0 , t 0 ) x y   =   φ 07 u 0   +   i   =   1 s φ i 7 u i + O ( h i 2 , k i 2 , τ i 2 ) , 2 u ( x 0 , y 0 , t 0 ) x t   =   φ 08 u 0   +   i = 1 s φ i 8 u i   +   O ( h i 2 , k i 2 , τ i 2 ) , 2 u ( x 0 , y 0 , t 0 ) y t   =   φ 09 u 0 + i   =   1 s φ i 9 u i + O ( h i 2 , k i 2 , τ i 2 ) ,
or, written in the vectorial form,
Π 9 u ( x 0 , y 0 , t 0 )   =   φ 0 u 0 + i   =   1 s φ i u i + O ( h i 2 , k i 2 , τ i 2 ) ,
where φ 0 and φ i stand for
φ 0   =   { φ 01 , φ 02 , φ 03 , φ 04 , φ 05 , φ 06 , φ 07 , φ 08 , φ 09 } T , φ i   =   { φ i 1 , φ i 2 , φ i 3 , φ i 4 , φ i 5 φ i 6 , φ i 7 , φ i 8 , φ i 9 } T ,
fulfilling
φ 0 = i   =   1 s φ i .
For the 3D case, the number of nodes per star should be greater or equal to 10 in order to obtain well-conditioned stars.
Remark 1.
The consistency of the GFD formulae was proved in [20] for the 2D case (which in our case is Γ R and Γ *   =   Γ   ×   [ 0 , T ] ) and for the 3D case (for our case, Γ R 2 and Γ *   =   Γ   ×   [ 0 , T ] ) in [21]. Moreover, the rate of convergence is quadratic, which is clear from the cited papers.

3. Numerical Results

In the present section, we give several numerical results obtained by solving the parabolic equations in 2D and 3D using STCM. In order to show the applicability of the present method, we compare the present method with the results obtained previously in the literature. In the following numerical examples, the distance criterion has been used, the number of nodes per star is eight plus the central node and the weighting function is the inverse of the distance cubed, η ( d i )   =   1 d i 3 with d i being the distance from node i to the central node of the star. More precisely, we compare the errors and execution times of the STCM and the GFDM. In all the cases under consideration, the accuracy of the STCM and the GFDM are computed using the expressions:
l 2   =   i   =   1 Z ( u i     U i ) 2 Z 1 2 , l   =   m a x | u i     U i | ,
where Z is the number of points in ( Γ Γ )   ×   ( 0 , T ] .

3.1. 2D Problems

For the two-dimensional examples (one space variable and one temporal variable), we use the clouds of points of Figure 3 and Figure 4.

3.1.1. Example 1

We consider the equation
π 2 u ( x , t ) t     2 u ( x , t ) x 2   =   e t [ π 2 x ( 1     x ) + 2 ]
with initial and boundary conditions:
u ( x , 0 )   =   sin ( π x )   +   x ( 1     x ) ; u ( 0 , t )   =   u ( 1 , t )   =   0
and whose exact solution is
u ( x , 0 )   =   e t [ sin ( π x )   +   x ( 1     x ) ] .
Table 1 shows the errors and execution times for the resolution of Example 1 in the clouds of Figure 3 and Figure 4, respectively, using STCM and GFDM explicitly (with step size Δ t   =   0.001 s for cloud 1 and Δ t   =   0.00001 for cloud 2). It can be seen that the error norms are very similar and the execution times of the STCM are slightly higher.

3.1.2. Example 2

For our second 2D examples, we choose
u ( x , t ) t     u ( x , t ) x     x 2 u ( x , t ) x 2 = e t [ ( x 2 x 1 ) sin x 3 cos x ]
together with
u ( x , 0 )   =   x sin ( x ) ; u ( 0 , t )   =   0 ; u ( 1 , t )   =   e t sin ( 1 ) .
The exact solution of problem (13) and (14) is
u ( x , t ) = e t [ x sin ( x ) ] .
The error norms, run times, and a comparison between the STCM and the explicit GFDM (step sizes are the same as in the previous example as the meshes are the same), for this second example in the clouds of points of Figure 3 and Figure 4, are depicted in Table 2. We observe that the errors are rather similar and the execution times are smaller for the STCM for a great number of points.
It is clear from the examples that, as the number of nodes increases, the execution times are shorter in the STCM than in the explicit GFDM, which is logical, since, as the time step decreases (in order to guarantee the stability conditions), the explicit GFDM has to make more iterations.

3.2. 3D Problems

To continue our study, we provide examples in the three-dimensional setting. We use both regular and irregular clouds of points, seen in Figure 5, Figure 6, Figure 7 and Figure 8. The black points denote the boundary nodes, the blue ones denote the interior nodes and the green points denote the Γ   ×   { t   =   T } boundary nodes. In the following numerical examples, the distance criterion has been used, the number of nodes per star is 24 plus the central node and the weighting function is the inverse of the distance squared, η ( d i )   =   1 d i 2 with d i being the distance from node i to the central node of the star.

3.2.1. Example 3

Consider the following equation
2 π 2 u ( x , y , t ) t 2 u ( x , y , t ) x 2 2 u ( x , y , t ) y 2   =   2 x ( 1     x ) ( 1     y )
with the initial condition
u ( x , y , 0 )   =   sin ( π x ) sin ( π y )   +   x ( 1     x ) ( 1     y ) .
By a direct check, the exact solution is (the Dirichlet boundary conditions are chosen so that the equation is fulfilled)
u ( x , y , t )   =   e t [ sin ( π x ) sin ( π y ) + x ( 1     x ) ( 1     y ) ]
We depict in Table 3 the error norms and run times for this third example. We found that for similar errors, the STCM consumes slightly less time than the explicit GFDM (where we use Δ t   =   0.001 s for cloud 3, Δ t = 0.00001 s for cloud 4, Δ t   =   0.05 s for cloud 5 and Δ t   =   0.0025 s for cloud 6). For this third example, we have also compared the STCM with the implicit scheme given by the GFDM formulae (with step size Δ t = 0.1 s for cloud 5 and Δ t   =   0.05 s for cloud 6). In this case, the discretisation of the derivatives is
u t   =   u 0 j   +   1     u 0 j 2 Δ t ; 2 u x 2   =   φ 04 u 0 j + 1 + i   =   1 s φ i 4 u i j   +   1 ; 2 u x 2   =   φ 05 u 0 j   +   1 + i   =   1 s φ i 5 u i j   +   1

3.2.2. Example 4

Finally, we solve the equation
u ( x , y , t ) t     u ( x , y , t ) x     u ( x , y , t ) y     x 2 u ( x , y , t ) x 2 y 2 u ( x , y , t ) y 2 = = e t [ ( ( x 2 x 1 ) sin x 3 x cos x ) y sin y + ( ( y 2 y 1 ) sin y 3 y cos y ) x sin x ] ,
with initial data
u ( x , y , 0 )   =   x y sin ( x ) sin ( y ) ,
and exact solution
u ( x , y , t ) = e t [ x y sin ( x ) sin ( y ) ] .
Again, in Table 4, we can observe that, even though the errors are rather similar for both methods, the STCM achieves better execution times.

4. Conclusions

The STCM allows the solution of the numerically parabolic PDEs in an efficient way without the need for any additional condition to the well-defined problem (initial and boundary conditions). In the various examples treated, both in 2D and 3D, linear and non-linear, in regular and irregular domains, the numerical results show that the errors are similar to those obtained with the explicit GFDM. The run times are also similar, but the numerical examples show that the time consumption of the STCM is smaller if the number of points in the computational domain is high.
The potential of the STCM consists of the elimination of the time step constraint forced by the stability of the explicit method, which is a real advance in the resolution of linear and non-linear parabolic PDEs.

Author Contributions

Formal analysis, A.M.V.; Methodology, F.U.; Project administration, J.J.B.; Resources, J.J.B.; Software, Á.G.; Writing—original draft, M.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge the support of the Escuela Técnica Superior de Ingenieros Industriales (UNED) of Spain, project 2021-IFC02. This work is also partially supported by the Project MTM2017-42907-P from MICINN (Spain).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, S.; Liu, W.K. Meshfree Particle Methods; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  2. Nguyen, V.P.; Rabczuk, T.; Bordas, S.; Duflot, M. Meshless methods: A review and computer implementation aspects. Math. Comput. Simul. 2008, 79, 763–813. [Google Scholar] [CrossRef] [Green Version]
  3. Huang, T.H.; Wei, H.; Chen, J.S.; Hillman, M.C. RKPM2D: An open-source implementation of nodally integrated reproducing kernel particle method for solving partial differential equations. Comput. Part. Mech. 2020, 7, 393–433. [Google Scholar] [CrossRef]
  4. Zhang, C.; Rezavand, M.; Zhu, Y.; Yu, Y.; Wu, D.; Zhang, W.; Wang, J.; Hu, X. SPHinXsys: An open-source multi-physics and multi-resolution library based on smoothed particle hydrodynamics. Comput. Phys. Commun. 2021, 267, 108066. [Google Scholar] [CrossRef]
  5. Jensen, P.S. Finite difference technique for variable grids. Comput. Struct. 1972, 2, 17–29. [Google Scholar] [CrossRef]
  6. Liszka, T.; Orkisz, J. The finite difference method at arbitrary irregular grids and its application in applied mechanics. Comput. Struct. 1980, 11, 83–95. [Google Scholar] [CrossRef]
  7. Orkisz, J. Finite difference method (Part, III). In Handbook of Computational Solid Mechanics; Kleiber, M., Ed.; Spriger: Berlin/Heidelberg, Germany, 1998. [Google Scholar]
  8. Perrone, N.; Kao, R. A general finite difference method for arbitrary meshes. Comput. Struct. 1975, 5, 45–58. [Google Scholar] [CrossRef]
  9. Benito, J.J.; Ureña, F.; Gavete, L. Influence of several factors in the generalized finite difference method. Appl. Math. Model. 2001, 25, 1039–1053. [Google Scholar] [CrossRef]
  10. Albuquerque-Ferreira, A.C.; Ureña, M.; Ramos, H. The generalized finite difference method with third-and fourth-order approximations and treatment of ill-conditioned stars. Eng. Anal. Bound. Elem. 2021, 127, 29–39. [Google Scholar] [CrossRef]
  11. Uddin, M.; Ali, H. The space-time kernel-based numerical method for Burgers’ equations. Mathematics 2018, 6, 212. [Google Scholar] [CrossRef] [Green Version]
  12. Sophy, T.; Silva, A.D.; Kribèche, A.; Sadat, H. An alternative space-time meshless method for solving transient heat transfer problems with high discontinuous moving sources. Numer. Heat Transf. Part B Fundam. 2016, 69, 377–388. [Google Scholar] [CrossRef] [Green Version]
  13. Ku, C.-Y.; Liu, C.-Y.; Yeih, W.; Liu, C.-S.; Fan, C.-M. A novel space-time meshless method for solving the backward heat conduction problem. Int. J. Heat Mass Transf. 2019, 130, 109–122. [Google Scholar] [CrossRef]
  14. Hamaidi, M.; Naji, A.; Charafi, A. Space-time localized radial basis function collocation method for solving parabolic and hyperbolic equations. Eng. Anal. Bound. Elem. 2016, 67, 152–163. [Google Scholar] [CrossRef]
  15. Li, Z.; Mao, X.Z. Global space-time multiquadric method for inverse heat conduction problem. Int. J. Numer. Methods Eng. 2011, 85, 355–379. [Google Scholar] [CrossRef]
  16. Lei, J.; Wei, X.; Wang, Q.; Gu, Y.; Fan, C.-M. A novel space-time generalized FDM for dynamic coupled thermoelasticity problems in heterogeneous plates. Arch. Appl. Mech. 2022, 92, 287–307. [Google Scholar] [CrossRef]
  17. Li, P.-W. Space-time generalized finite difference nonlinear model for solving unsteady Burgers’ equations. Appl. Math. Lett. 2021, 114, 106896. [Google Scholar] [CrossRef]
  18. Qu, W.; He, H. A spatial-temporal GFDM with an additional condition for transient heat conduction analysis of FGMs. Appl. Math. Lett. 2020, 110, 106579. [Google Scholar] [CrossRef]
  19. Ureña, F.; Gavete, L.; Garcia, A.; Benito, J.J.; Vargas, A.M. Solving second order non-linear parabolic PDEs using generalized finite difference method (GFDM). J. Comput. Appl. Math. 2019, 354, 221–241. [Google Scholar] [CrossRef]
  20. Gavete, L.; Ureña, F.; Benito, J.J.; Garcia, A.; Ureña, M.; Salete, E. Solving second order non-linear elliptic partial differential equations using generalized finite difference method. J. Comput. Appl. Math. 2017, 318, 378–387. [Google Scholar] [CrossRef]
  21. Ureña, F.; Gavete, L.; Garcia, A.; Benito, J.J.; Vargas, A.M. Non-linear Fokker-Planck equation solved with generalized finite differences in 2D and 3D. Appl. Math. Comput. 2020, 368, 124801. [Google Scholar] [CrossRef]
Figure 1. Computational domain of the problem, where the green nodes belongs to Γ × [ 0 , T ] and Γ , and the blue nodes correspond to Γ × [ 0 , T ] .The solution of the original problem corresponds to the solution at Γ × { t = T } (green nodes).
Figure 1. Computational domain of the problem, where the green nodes belongs to Γ × [ 0 , T ] and Γ , and the blue nodes correspond to Γ × [ 0 , T ] .The solution of the original problem corresponds to the solution at Γ × { t = T } (green nodes).
Mathematics 10 01870 g001
Figure 2. On the left, the distance criterion to choose the nodes of a star in the interior of Γ × [ 0 , T ) where we choose the eight nearest nodes. On the right, the star is centered at a node of Γ × { t = T } (the star is centered at the black point, the yellow points belong to Γ × [ 0 , T ) and the red ones to Γ × [ 0 , T ] ).
Figure 2. On the left, the distance criterion to choose the nodes of a star in the interior of Γ × [ 0 , T ) where we choose the eight nearest nodes. On the right, the star is centered at a node of Γ × { t = T } (the star is centered at the black point, the yellow points belong to Γ × [ 0 , T ) and the red ones to Γ × [ 0 , T ] ).
Mathematics 10 01870 g002
Figure 3. First cloud of points for the 2D examples (Cloud 1).
Figure 3. First cloud of points for the 2D examples (Cloud 1).
Mathematics 10 01870 g003
Figure 4. Second cloud of points for the 2D examples (Cloud 2).
Figure 4. Second cloud of points for the 2D examples (Cloud 2).
Mathematics 10 01870 g004
Figure 5. First regular cloud of points for the 3D examples (Cloud 3).
Figure 5. First regular cloud of points for the 3D examples (Cloud 3).
Mathematics 10 01870 g005
Figure 6. Second regular cloud of points for the 3D examples (Cloud 4).
Figure 6. Second regular cloud of points for the 3D examples (Cloud 4).
Mathematics 10 01870 g006
Figure 7. First irregular cloud of points for the 3D examples (Cloud 5).
Figure 7. First irregular cloud of points for the 3D examples (Cloud 5).
Mathematics 10 01870 g007
Figure 8. Second irregular cloud of points for the 3D examples (Cloud 6).
Figure 8. Second irregular cloud of points for the 3D examples (Cloud 6).
Mathematics 10 01870 g008
Table 1. Errors and execution time of the Example 1 in the clouds of Figure 3 and Figure 4.
Table 1. Errors and execution time of the Example 1 in the clouds of Figure 3 and Figure 4.
Cloud 1 l 2 l t ( s )
STCM0.0009020.0016370.183955
Explicit-GFDM0.0008130.0013980.029527 s
Cloud 2 l 2 l t ( s )
STCM0.0002480.00057320.2965 s
Explicit-GFDM0.0002770.00040111.723 s
Table 2. Errors and execution time of the Example 2 in the clouds of Figure 3 and Figure 4.
Table 2. Errors and execution time of the Example 2 in the clouds of Figure 3 and Figure 4.
Cloud 1 l 2 l t ( s )
STCM0.0633480.0840770.10168 s
Explicit-GFDM0.0590030.0803380.096604 s
Cloud 2 l 2 l t ( s )
STCM0.0096180.01320426.145s
Explicit-GFDM0.0097870.01751044.409 s
Table 3. Errors and execution time of the Example 3 in the clouds of Figure 5, Figure 6, Figure 7 and Figure 8.
Table 3. Errors and execution time of the Example 3 in the clouds of Figure 5, Figure 6, Figure 7 and Figure 8.
Cloud 3 l 2 l t ( s )
STCM0.0165090.0325670.047252 s
Explicit-GFDM0.0152630.0252340.04283 s
Cloud 4 l 2 l t ( s )
STCM0.0008060.0062752.1954 s
Explicit-GFDM0.0008220.0015412.3099 s
Cloud 5 l 2 l t ( s )
STCM0.0036720.0055810.010077 s
Explicit-GFDM0.0039230.0058180.007449 s
Implicit-GFDM0.0038040.0057420.01121 s
Cloud 6 l 2 l t ( s )
STCM0.0008160.0016230.18481s
Explicit-GFDM0.0007330.0012830.26016 s
Implicit-GFDM0.0008310.0024320.19145 s
Table 4. Errors and execution time of Example 4 in the clouds of Figure 5, Figure 6, Figure 7 and Figure 8.
Table 4. Errors and execution time of Example 4 in the clouds of Figure 5, Figure 6, Figure 7 and Figure 8.
Cloud 3 l 2 l t ( s )
STCM0.0043380.0083770.87247 s
Explicit-GFDM0.0040980.0078601.1949 s
Cloud 4 l 2 l t ( s )
STCM0.0041360.0063943.2172 s
Explicit-GFDM0.0041660.0065673.4186 s
Cloud 5 l 2 l t ( s )
STCM0.0033830.0051830.35261 s
Explicit-GFDM0.0029230.0049870.5059 s
Cloud 6 l 2 l t ( s )
STCM0.0008710.0012761.94437 s
Explicit-GFDM0.0009310.0019052.2854 s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Benito, J.J.; García, Á.; Negreanu, M.; Ureña, F.; Vargas, A.M. A Novel Spatio-Temporal Fully Meshless Method for Parabolic PDEs. Mathematics 2022, 10, 1870. https://doi.org/10.3390/math10111870

AMA Style

Benito JJ, García Á, Negreanu M, Ureña F, Vargas AM. A Novel Spatio-Temporal Fully Meshless Method for Parabolic PDEs. Mathematics. 2022; 10(11):1870. https://doi.org/10.3390/math10111870

Chicago/Turabian Style

Benito, Juan José, Ángel García, Mihaela Negreanu, Francisco Ureña, and Antonio M. Vargas. 2022. "A Novel Spatio-Temporal Fully Meshless Method for Parabolic PDEs" Mathematics 10, no. 11: 1870. https://doi.org/10.3390/math10111870

APA Style

Benito, J. J., García, Á., Negreanu, M., Ureña, F., & Vargas, A. M. (2022). A Novel Spatio-Temporal Fully Meshless Method for Parabolic PDEs. Mathematics, 10(11), 1870. https://doi.org/10.3390/math10111870

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop