Direct Meshless Local Petrov–Galerkin Method Research Articles

Numerical study of the multi-dimensional Galilei invariant fractional advection–diffusion equation using direct mesh-less local Petrov–Galerkin method

This article presents a local mesh-less procedure for simulating the Galilei invariant fractional advection–diffusion (GI-FAD) equations in one, two, and three-dimensional spaces. The proposed method combines a second-order Crank–Nicolson scheme for time discretization and the second-order weighted and shifted Grünwald difference (WSGD) formula. This time discretization scheme ensures unconditional stability and convergence with an order of O(τ2). In the spatial domain, a mesh-less weak form is employed based on the direct mesh-less local Petrov–Galerkin (DMLPG) method. The DMLPG method employs the generalized moving least-square (GMLS) approximation in conjunction with the local weak form of the equation. By utilizing simple polynomials as shape functions in the GMLS approximation, the necessity for complex shape function construction in the MLS approximation is eliminated. To validate and demonstrate the effectiveness of the proposed algorithm, a variety of problems in one, two, and three dimensions are investigated on both regular and irregular computational domains. The numerical results obtained from these investigations confirm the accuracy and reliability of the developed approach in simulating GI-FAD equations.

Multi-material topology optimization with a direct meshless local Petrov–Galerkin method applied to a two-dimensional boundary value problem

ABSTRACT This work presents a novel approach for the solution of multi-material topology optimization problems of elastic structures coupling the Direct Meshless Local Petrov–Galerkin (DMLPG) method with the Gradual Bi-direction Evolutionary Structural Optimization (gBESO) method. In recent years, meshless methods have been used in several engineering fields, showing some advantages compared to mesh-based methods. DMLPG is a truly meshless method, which achieves results with good accuracy and computational efficiency, avoiding the numerical integration of complicated shape functions by adopting low-degree polynomials. In the proposed algorithm, DMLPG is used to obtain smooth nodal displacements, strains and stresses, and gBESO updates the structural geometry based on design sensitivity values of the compliance while it gradually increases the Young's modulus of the material. The optimization problem minimizes the compliance objective function subject to volume constraints. Numerical examples demonstrate the applicability and validity of the technique.

Simulation of predator–prey system with two-species, two chemicals and an additional chemotactic influence via direct meshless local Petrov–Galerkin method

PurposeThis paper aims to introduce a new numerical approach based on the local weak form and the Petrov–Galerkin idea to numerically simulation of a predator–prey system with two-species, two chemicals and an additional chemotactic influence.Design/methodology/approachIn the first proceeding, the space derivatives are discretized by using the direct meshless local Petrov–Galerkin method. This generates a nonlinear algebraic system of equations. The mentioned system is solved by using the Broyden’s method which this technique is not related to compute the Jacobian matrix.FindingsThis current work tries to bring forward a trustworthy and flexible numerical algorithm to simulate the system of predator–prey on the nonrectangular geometries.Originality/valueThe proposed numerical results confirm that the numerical procedure has acceptable results for the system of partial differential equations.

Numerical solution of an inverse source problem for a time-fractional PDE via direct meshless local Petrov–Galerkin method

In this paper, we employ an improved Meshless Local Petrov–Galerkin (MLPG), namely Direct Meshless Local Petrov–Galerkin (DMLPG) method, for solving an inverse time-dependent source problem for two-dimensional fractal mobile/immobile solute transport equation on regular and irregular domains. In the weak form DMLPG method, numerical integrations apply over low-degree polynomial basis functions instead of complicated moving least squares (MLS) shape functions of MLPG method. Therefore, the computational costs often reduce in the DMLPG method. In space domain, we employ the DMLPG method and the time derivatives of problem are discretized based on finite difference schemes. Numerical results demonstrate efficiency and accuracy of the proposed algorithm for solving the inverse source problem.

Numerical study of the unsteady 2D coupled magneto-hydrodynamic equations on regular/irregular pipe using direct meshless local Petrov–Galerkin method

The study of the conduction of fluids in magnetic fields has attracted a lot of interest in recent years. In this paper, we use new truly meshless methods for numerical solving of the unsteady 2D coupled magneto-hydrodynamic (MHD) equations with uniform and scattered distributions of nodes on regular and polar domains. These methods are implemented based on a generalization of the moving least squares approximation (GMLS) method and local weak forms of MHD equations. By observing the numerical results obtained from these methods and comparing them with others, one can find the efficiency, accuracy, and speed of these methods.

Numerical solution strategy for natural convection problems in a triangular cavity using a direct meshless local Petrov-Galerkin method combined with an implicit artificial-compressibility model

Thermal analysis of natural convection in a cavity is solved by numerical approach. The direct meshless local Petrov–Galerkin method combined with an implicit artificial-compressibility model is proposed here to simulate laminar natural convective heat transfer in triangular cavities. The spatial terms of the governing equation are discretized using the direct meshless local Petrov–Galerkin method with the Dirac function as the test function. The method uses an artificial-compressibility scheme to couple the pressure with the velocity components directly. The semi-algebraic system is solved using an implicit backward Euler pseudo-time method. Kovasznay flow problem is presented to examine the accuracy of proposed method. The validation of the obtained results by using experimental data is also provided, and four test problems — natural convection in an equilateral triangular cavity, a right-angled triangular cavity, an isosceles triangular cavity, and an irregular triangular cavity — are solved by using the proposed method. The numerical results obtained using the proposed method showed excellent agreement with those obtained using the conventional methods available in the literature.

Investigation of generalized Couette hydromagnetic flow of two-step exothermic chemical reaction in a channel via the direct meshless local Petrov–Galerkin method

This paper concerns to propose a high-order local meshless method to simulate the Couette flow based on the incompressible Navier–Stokes equation and the generalized Couette hydromagnetic flow of two-step exothermic chemical reaction in a channel. The numerical solution is performed in the two main stages. First, time variable of the incompressible Navier–Stokes equation is approximated by a finite difference scheme as in this procedure, the velocities in x- and y-directions must be updated by using an explicit difference formula. Also, the pressure component is calculated by solving a Poisson’s equation in each time step. Here is the second stage e.g. a local meshless method based on the direct meshless local Petrov–Galerkin method is employed to solve the Poisson’s equation which corresponds to the pressure value. Some famous phenomena such as Couette flow are studied by using the proposed numerical algorithm. The results show that the new numerical procedure not only has suitable and acceptable accuracy but also it simply solves such interesting models. Up to the best of our knowledge, this is the first time that the direct meshless local Petrov–Galerkin method is applied for these problems and it can have more accuracy than other meshless methods.

An efficient meshfree method based on Pascal polynomials and multiple-scale approach for numerical solution of 2-D and 3-D second order elliptic interface problems

In this work, we propose an efficient meshfree method based on Pascal polynomials and multiple-scale approach for numerical solutions of two-dimensional (2-D) and three-dimensional (3-D) elliptic interface problems which may have discontinuous coefficients and curved interfaces with or without sharp corners. The proposed method uses Pascal polynomials as basis functions and utilizes multiple-scale approach for stabilizing numerical solutions. It is a well known fact that using polynomial basis without any modifications to obtain numerical solutions of partial differential equations may be peculiar owing to ill conditioned resultant coefficient matrix which is formed after process of discretization. Hence, as a remedy to the highly ill-conditioned coefficient matrix we employ multiple-scale approach. The main idea behind the multiple-scale approach is to make norm of all columns of resultant coefficient matrix equal to each other. The proposed method is a truly strong-form meshfree method since we do not need any mesh or integration process in problem domain, these features makes the implementation of the method very simple in computer environment. The efficiency of the proposed method is tested by some test problems which may have smooth interface or interface with sharp corners. Stability of the proposed method is investigated numerically in the presence of noise effect. Further, to show accuracy of the proposed method we present some comparisons with available numerical methods in literature, such as direct meshless local Petrov-Galerkin method, matched interface and boundary method, spectral element method and some meshless methods based on radial basis functions. The obtained numerical results and their comparisons confirm applicability of the proposed method for 2-D and 3-D steady state elliptic interface problems.

Direct meshless local Petrov-Galerkin method to investigate anisotropic potential and plane elastostatic equations of anisotropic functionally graded materials problems

The meshless local technique is employed to simulate potential problems involving material anisotropy and also plane elastostatic equations of anisotropic functionally graded materials (PEE-AFGM). We utilize direct meshless local Petrov Galerkin (DMLPG) technique to discrete the spatial directions. The DMLPG5 (based on weak form) is applied to solve potential problems involving material anisotropy and the DMLPG2 (based on strong form) is employed to simulate PEE-AFGM. The coefficients of the anisotropic conductivity require the gradient term to be modified accordingly. The grading function can be any general function that varies smoothly from point to point in the material. The developed procedures are tested on various examples on non-rectangular domains to check their accuracy.

DMLPG method for specifying a control function in two-dimensional parabolic inverse PDEs

This article describes a fast meshless technique, the direct meshless local Petrov–Galerkin (DMLPG) method, for identifying a control function in two-dimensional parabolic inverse PDEs. The DMLPG method is a truly meshless technique that directly approximates the local weak forms. In DMLPG, the local integrations are done over the low-degree polynomial basis functions instead of the complicated shape functions. Thus, the mass and stiff matrices are constructed by integration against polynomials. This feature overcomes the main drawback of meshless methods in comparison with the finite element methods (FEM) and significantly increases the computational efficiency in comparison with the other meshless weak form methods. In this paper, we apply the DMLPG and MLPG methods for solving two-dimensional parabolic inverse PDEs on regular and irregular domains. These methods are compared with each other and the superiority of the DMLPG over the classical MLPG is demonstrated. The numerical results confirm the good efficiency of the DMLPG method for solving inverse problems especially coefficient inverse problems. Finally, we will conclude that the DMLPG method can be considered as an attractive alternative to existing meshless weak form methods in solving inverse problems.

Numerical simulation of mixed convection in a square cavity using Direct Meshless Local Petrov-Galerkin (DMLPG) method.

This paper present the extension of the Direct Meshless Petrov-Galerkin (DMLPG) for solving mixed convection heat transfer in a square cavity for the flow up to Reynolds number up to Re 400 and Grashof number Gr = 100. The semi algebraic forms of the governing equations are solved using fractional step method. The numerical results of the DMLPG method are in good agreement with the conventional numerical results in the literature. Based on the numerical results it is shown that The DMLPG is very promising in dealing with the fluid flow and heat transfer problem.

DMLPG method for numerical simulation of soliton collisions in multi-dimensional coupled damped nonlinear Schrödinger system which arises from Bose–Einstein condensates

In this paper, the direct meshless local Petrov–Galerkin (DMLPG) method is applied to find the numerical solution of coupled damped nonlinear Schrödinger system in one, two and three-dimensional spaces. The propagation properties of single soliton, double and triple solitons of coupled damped nonlinear Schrödinger system are simulated and the interactions between these solitons are studied numerically. The efficient time differencing Runge–Kutta method is utilized for the time discretization. DMLPG shifts the numerical integrations over low-degree polynomials rather than over complicated shape functions and this significantly increases the computational efficiency of DMLPG in comparison with the other meshless local weak form methods especially in two and three dimensions. The main aim of this paper is to show that the DMLPG method can be simply used for solving high-dimensional system of non-linear partial differential equations especially coupled damped nonlinear Schrödinger system. The numerical results confirm the good efficiency of the proposed method for solving our model.

Direct meshless local Petrov–Galerkin (DMLPG) method for 2D complex Ginzburg–Landau equation

In this paper, the direct meshless local Petrov–Galerkin (DMLPG) approximation is proposed to solve the two-dimensional complex Ginzburg–Landau equation. DMLPG uses a generalized MLS (GMLS) method which directly approximates test functionals by values at nodes and therefore avoids integration over MLS shape functions and replaces it by a much cheaper integration over polynomials. An Implicit-Explicit method used to deal with the time derivative. Several examples are performed to compare the accuracy and efficiency of the DMLPG method with the traditional MLPG method.

Numerical investigation based on direct meshless local Petrov Galerkin (direct MLPG) method for solving generalized Zakharov system in one and two dimensions and generalized Gross–Pitaevskii equation

In the current investigation, we propose two efficient meshless numerical techniques for solving two models in optic engineering, i.e. the generalized Gross–Pitaevskii equation and the generalized Zakharov system. Two local meshless methods have been employed for solving these models: local radial basis functions collocation method and direct meshless local Petrov–Galerkin method. In this paper, we discrete the space direction using the local radial basis functions collocation and direct meshless local Petrov–Galerkin techniques and to obtain high-order numerical results, we use the fourth-order exponential time differencing Runge–Kutta method for discretizing the temporal direction. The obtained numerical results are compared with some well-known numerical techniques. Moreover, several examples are given that show the acceptable accuracy and efficiency of the proposed scheme.

Direct meshless local Petrov–Galerkin method for the two-dimensional Klein–Gordon equation

In this paper we apply the direct meshless local Petrov–Galerkin (DMLPG) method to solve the two dimensional Klein–Gordon equations in both strong and weak forms. Low computational cost is the main property of this method compared with the original MLPG technique. The reason lies behind the approach of generalized moving least squares approximation where the discretized functionals, obtained from the PDE problem, are directly approximated from nodal values. This shifts the integration over polynomials rather than the MLS shape functions, leading to an extremely faster scheme. We will see that this method can successfully solve the problem with a reasonable accuracy.

Remediation of contaminated groundwater by meshless local weak forms

In this paper, four meshless local weak form methods such as direct meshless local Petrov–Galerkin method (DMLPG), meshless local Petrov–Galerkin method (MLPG), local weak radial point interpolation method (LWRPIM) and local weak moving kriging method (LWMK) are applied to find the numerical solution of coupled non-linear advection–diffusion–reaction system arising in the prevention of groundwater contamination problem. A comparison between these methods is done from the perspective of accuracy and computational efficiency. An efficient fourth-order exponential time differencing Runge–Kutta method is utilized for the time discretization. The computational efficiency is the most significant advantage of the DMLPG method in comparison with the other meshless local weak form methods, because DMLPG reduces the computational costs, significantly. This is due to the fact that DMLPG shifts the numerical integrations over low-degree polynomials rather than over complicated shape functions. The main aim of this paper is to show that the meshless local weak form methods can be used for solving the system of non-linear partial differential equations especially coupled non-linear advection–diffusion–reaction system. The numerical results confirm the good efficiency of the proposed methods for solving our model.

Application of direct meshless local Petrov–Galerkin (DMLPG) method for some Turing-type models

Mathematical modeling of pattern formation in developmental biology leads to non-linear reaction---diffusion systems which are usually highly stiff in both diffusion and reaction terms. In this paper, the Direct Meshless Local Petrov---Galerkin (DMLPG) procedure is applied to find the numerical solution of some non-linear time-dependent reaction---diffusion systems such as Schnakenberg model, Gierer---Meinhardt model, FitzHugh---Nagumo model and Gray---Scott model. As far as we know, it is the first time that DMLPG method is applied for solving non-linear partial differential equations (PDEs) and systems of PDEs. Computational efficiency is the most significant advantage of the DMLPG method in comparison with the classic Meshless Local Petrov---Galerkin (MLPG) method. This is due to the fact that DMLPG shifts the numerical integrations over low-degree polynomials instead of over complicated moving least squares (MLS) shape functions and this reduces the computational costs, significantly. The main aim of this paper is to show that the DMLPG method is also suitable for solving the non-linear time-dependent systems, especially reaction---diffusion systems. Numerical results support the good efficiency of the proposed method for solving non-linear reaction---diffusion systems. Also it is shown that DMLPG provides considerable savings in computational time in comparison with the classical MLPG method.

Direct meshless local Petrov–Galerkin method for elliptic interface problems with applications in electrostatic and elastostatic

In recent years, there have been extensive efforts to find the numerical methods for solving problems with interface. The main idea of this work is to introduce an efficient truly meshless method based on the weak form for interface problems. The proposed method combines the direct meshless local Petrov–Galerkin method with the visibility criterion technique to solve the interface problems. It is well-known in the classical meshless local Petrov–Galerkin method, the numerical integration of local weak form based on the moving least squares shape function is computationally expensive. The direct meshless local Petrov–Galerkin method is a newly developed modification of the meshless local Petrov–Galerkin method that any linear functional of moving least squares approximation will be only done on its basis functions. It is done by using a generalized moving least squares approximation, when approximating a test functional in terms of nodes without employing shape functions. The direct meshless local Petrov–Galerkin method can be a very attractive scheme for computer modeling and simulation of problems in engineering and sciences, as it significantly uses less computational time in comparison with the classical meshless local Petrov–Galerkin method. To create the appropriate generalized moving least squares approximation in the vicinity of an interface, we choose the visibility criterion technique that modifies the support of the weight (or kernel) function. This technique, by truncating the support of the weight function, ignores the nodes on the other side of the interface and leads to a simple and efficient computational procedure for the solution of closed interface problems. In the proposed method, the essential boundary conditions and the jump conditions are directly imposed by substituting the corresponding terms in the system of equations. Also, the Heaviside step function is applied as the test function in the weak form on the local subdomains. Some numerical tests are given including weak and strong discontinuities in the Poisson interface problem. To demonstrate the application of these problems, linearized Poisson–Boltzmann and linear elasticity problems with two phases are studied. The proposed method is compared with analytical solution and the meshless local Petrov–Galerkin method on several test problems taken from the literature. The numerical results confirm the effectiveness of the proposed method for the interface problems and also provide significant savings in computational time rather than the classical meshless local Petrov–Galerkin method.