Meshless Local Petrov-Galerkin Research Articles

Disclosure of the intricacies in coupled groundwater flow and contaminant transport using mesh-less local Petrov–Galerkin method

Disclosure of the intricacies in coupled groundwater flow and contaminant transport using mesh-less local Petrov–Galerkin method

Fourth-order exponential time differencing Runge–Kutta scheme and local meshless method to investigate unsteady diffusion–convection problems of anisotropic functionally graded materials

The current article is devoted to proposing a local meshfree method for simulating the unsteady diffusion–convection problems of anisotropic functionally graded materials (AFGM). The first stage of the proposed numerical algorithm is based on discretizing the spatial variables by utilizing the direct meshless local Petrov–Galerkin (DMLPG) method. Then a system of ODEs is extracted which is associated to the time variable. Furthermore, the fourth-order exponential time differencing Runge–Kutta method is applied to solve the system of ODEs. The main aim of this work is to introduce a flexible and low-cost numerical procedure to solve unsteady diffusion–convection problems of AFGM on the complicated geometry. The new numerical technique has admissible resultants to simulate the studied model. This numerical formulation is applied to several examples with complex domains to check its ability.

Aquifer flow parameter estimation using coupled meshless methods and metaheuristic algorithms

A reliable analysis of various groundwater problems requires accurate input of aquifer parameters. However, field measurement of such parameters is tedious and expensive. Inverse modelling by Simulation-Optimization (SO) resolves this limitation. In this study, the unknown transmissivity of confined aquifers is estimated using SO models. Three simulation models of strong, weak and hybrid categories of meshless methods, namely, Radial Point Collocation Method (RPCM), Meshless Local Petrov Galerkin (MLPG) and Meshless Weak Strong (MWS) form, are coupled with metaheuristic algorithms of Differential Evolution (DE), Particle Swarm Optimization (PSO) and Whale Optimization Algorithm (WOA), resulting in nine SO models. Five of these nine models are novel SO models and first-time application of WOA for groundwater flow parameter estimation. The application of models to heterogeneous hypothetical and complex field-type aquifers prove that solutions are similar to true transmissivities. This study provides an insight into the selection of suitable SO model based on available resources and requirements.

OpenMP-based parallel MLPG solver for analysis of heat conduction

PurposeIn the present work, we focus on developing an in-house parallel meshless local Petrov-Galerkin (MLPG) code for the analysis of heat conduction in two-dimensional and three-dimensional regular as well as complex geometries.Design/methodology/approachThe parallel MLPG code has been implemented using open multi-processing (OpenMP) application programming interface (API) on the shared memory multicore CPU architecture. Numerical simulations have been performed to find the critical regions of the serial code, and an OpenMP-based parallel MLPG code is developed, considering the critical regions of the sequential code.FindingsBased on performance parameters such as speed-up and parallel efficiency, the credibility of the parallelization procedure has been established. Maximum speed-up and parallel efficiency are 10.94 and 0.92 for regular three-dimensional geometry (343,000 nodes). Results demonstrate the suitability of parallelization for larger nodes as parallel efficiency and speed-up are more for the larger nodes.Originality/valueFew attempts have been made in parallel implementation of the MLPG method for solving large-scale industrial problems. Although the literature suggests that message-passing interface (MPI) based parallel MLPG codes have been developed, the OpenMP model has rarely been touched. This work is an attempt at the development of OpenMP-based parallel MLPG code for the very first time.

Analysis and application of MLPG7 for diffusion equations with nonlinear reaction terms

This paper extends the recently proposed variant of meshless local Petrov Galerkin (MLPG) method, i.e., MLPG7, for solving time dependent PDEs. As test function, the method uses a novel modification of fundamental solution of Laplace operator that not only the test function itself but also its derivative vanish on boundary of local subdomains. Therefore, more stable local integral equations are obtained by replacing a boundary integral of derivatives of field variables with a domain integral of field variables itself. MLPG7 formulation converts the nonlinear governing equation into a system of first order ordinary differential equations (ODEs). For this system of ODEs, a theoretical stability analysis is carried out for the case of regularly spaced nodal points. Stability and convergence conditions for fully discretized equations by Euler (forward and backward) and Crank–Nicolson methods are investigated. The nonlinear term is treated iteratively. A numerical study investigates the convergence and stability of the method. A comparison is also done with some well-known methods and the results reveal the excellence of MLPG7. The convergence of the iterative approach is numerically verified by another illustrative test example. As nonlinear test problems, Allen–Cahn equations are studied.

A coupled flow and transport model for simulation of multi-species reactive transport in unconfined aquifer using meshless local Petrov Galerkin (MLPG) method.

An understanding of natural degradation of multiple reactive contaminants in the aquifers is essential before designing the monitoring or remediation programs for polluted aquifers. Since such reactive contaminants are ubiquitous, a number of research works has been performed in the past three decades for the modelling of multi-species reactive transport (MSRT) phenomenon. The widely used finite difference method (FDM) and finite element method (FEM)-based models suffer a drawback of relying on a grid/mesh, which makes the solution unstable. Addressing such difficulties, the latest research on the MSRT models is directed towards the meshless methods. In this study, the meshless local Petrov Galerkin (MLPG) method-based multi-species reactive transport model (MLPG-MSRT) is presented, with an objective to create a robust simulation tool for the prediction of fate of multiple contaminants of the first-order reaction network. The developed model is validated for reversible as well as irreversible reaction networks with the available analytical solutions. Also, the MLPG model for unconfined aquifer flow (UF) is developed, validated, and coupled with the MLPG-MSRT model. The MLPG-UF-MSRT model results are further compared with the established FDM-based MODFLOW-RT3D model solutions for a rectangular and a real field type study. The results showed that the proposed model can simulate MSRT as accurately as the FDM-based models with an additional advantage of simplicity and stability, and thus, is more efficient for complex field problems.

A new meshless local integral equation method

This paper proposes a novel meshless local integral equation (LIE) method for numerical solutions of two and three-dimensional Poisson equations. The proposed method can be regarded as a new variant of the meshless local Petrov-Galerkin (MLPG) method which already has six variants, and so the proposed method can be called MLPG7. A variant of MLPG called local boundary integral equation (LBIE) method, uses a localized fundamental solution of Laplace equation as test function. This test function vanishes on local boundaries and therefore, the derived local equations do not involve the gradient of field variables. The motivation for new formulation is using an average of LBIE equation over radius of local sub-domains instead of a single sub-domain. A new kernel is introduced which is a new modification of fundamental solutions of Laplace equation. It is proved that using the new kernel as test function leads to the same local equations as averaging. Infact, the new formulation is a modification of the LBIE method for which the boundary integral is replaced with a domain integral. It is proved that new local weak formulation is equivalent with strong form. The convergence of the method for the case of regularly spaced nodal points is proved. For trial approximation, this paper uses Gaussian radial basis functions (Gaussian RBF), however, derivative free property of the method allows the use of any non-differentiable basis functions. Stability and convergence of the method are numerically studied by considering both 2D and 3D problems with regular and irregular nodal points. The method is compared with LBIE, finite differences (FD) and RBF-FD methods and the numerical results reveal the significance of the proposed method.

Petrov–Galerkin zonal free element method for 2D and 3D mechanical problems

AbstractIn this article, a novel weak‐form zonal Petrov–Galerkin free element method is proposed for two‐ and three‐dimensional linear mechanical problems. By absorbing the advantages of finite block method and strong‐form finite element method, the block mapping technique is used in the free element method. Combining the characteristics of the meshless local Petrov–Galerkin method, the local Petrov–Galerkin formulation based on the zonal free element method is formed at last. Besides, the local integral domain selected in the local collocation element is circular or spherical to simplify programming. The transformation of the local integral domain between the physical and normalized spaces is given for two‐ and three‐dimensional problems. The comparison of accuracy and convergence between the new proposed Petrov–Galerkin method and the conventional methods is carried out. Some challenging examples including fracture mechanics problems and a complex 3D problem are given to validate the convergence and accuracy of the proposed method.

A Multi-Scale Continuum Study for Boron Nitride Nanotubes Based on the Meshless Local Petrov–Galerkin Method

A Multi-Scale Continuum Study for Boron Nitride Nanotubes Based on the Meshless Local Petrov–Galerkin Method

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.

The direct meshless local Petrov–Galerkin technique with its error estimate for distributed-order time fractional Cable equation

Distributed-order fractional calculus is a quickly growing concept of the more general area of fractional calculus that has significant and extensive usage for designing complex systems. This work is used the direct meshless local Petrov–Galerkin (DMLPG) technique for the numerical solution of the distributed-order time fractional Cable equation. DMLPG implements a generalized moving least square (GMLS) process to discretize the equation in space variables. By using this scheme, the test function is approximated via the values at nodes, directly. Thus, this algorithm passes integration with the MLS shape functions substituting with a more inexpensive integration than polynomials. Here, the distributed integral is discretized by the M-point Gauss–Legendre quadrature rule. Then, the finite difference scheme is applied to approximate the fractional derivative discretization. Also, the unconditionally stability and rate of convergence O(τ2−max{α,β}) of the time-discrete technique are demonstrated. Moreover, the current method converts the problem into a system of linear algebraic equations. To demonstrate the capability and flexibility of our scheme, some examples with different geometric domains are supposed in two-dimensional cases.

Numerical study of the variable-order time-fractional mobile/immobile advection-diffusion equation using direct meshless local Petrov-Galerkin methods

In this paper, we use direct meshless local Petrov-Galerkin (DMLPG) methods for solving the variable-order time-fractional mobile/immobile advection-diffusion equation in two dimensions. The basis of the DMLPG methods is on the generalized moving least-square (GMLS) approximation and local weak form of the equation. Since the GMLS approximation uses basic polynomials as shape functions, there is no need to construct complex shape functions in the moving least-square (MLS) approximation. In addition, the calculation of numerical integrals in the local weak form of the equation is simple. We also give proof of the convergence and stability of the time-discrete scheme. The analytical results show that the semi-discrete time scheme is unconditionally stable and the convergence order is O(Δt2). Moreover, to show the accuracy and efficiency of the methods, we use the regular and irregular domains with distributions of uniform and scattered nodes in numerical examples.

Inverse modeling application for aquifer parameters estimation using a precise simulation–optimization model

In this research, a simulation–optimization (S/O) model is used in order to estimate aquifer parameters on two aquifers. In this model, meshless local Petrov–Galerkin (MLPG) is used for simulation purpose and modified teaching–learning-based optimization (MTLBO) algorithm is engaged as optimization model. Linking these two powerful models generates a S/O model named MLPG-MTLBO. The proposed model is applied on two aquifers: a standard and a real field aquifer. In standard aquifer, parameters are only transmissivity coefficients in x and y direction for three zones. The acquired results by MLPG-MTLBO are really close to true values. This fact presents the power of MLPG-MTLBO inverse model. Therefore, it is applied on field aquifer. Unconfined aquifer of Birjand recognized as real case study. Parameters which are needed to be estimated are specific yields and hydraulic conductivity coefficients. These parameters are computed by MLPG-MTLBO and entered to the groundwater flow model. The achieved groundwater table compared with observation data and RMSE is calculated. RMSE value is 0.356 m; however, this error criterion for MLPG and FDM is 0.757 m and 1.197 m, respectively. This means that estimation is precise and makes the RMSE to reduce from 0.757 to 0.356 m, and also, MLPG-MTLBO is an accurate model for this aim.

Numerical investigation of wave interaction with two closely spaced floating boxes using particle method

This paper applies Improved Meshless Local Petrov-Galerkin method based on Rankine source function (IMLPG_R) to simulate two-dimensional wave interaction with two closely spaced floating boxes. The present study considers four different box configurations based on the degree of freedom, namely FF (both boxes fixed), HH (only heave allowed), RR (only roll allowed), and MM (both heave and roll allowed). The influences of the different degrees of freedom of the boxes on the wave elevation in the gap, and forces acting on the boxes, are investigated. Regular waves with a wave period ranging from 0.6s to 1.6s and two focused wave groups based on the JONSWAP spectrum with a peak frequency of 0.869Hz and 1.081Hz, respectively, are employed for the analysis. The piston mode resonance frequencies of different configurations were observed to vary slightly depending on the motion restrained. Compared to the FF configuration, the MM configuration showed a 28.3% increase in wave amplification in the gap. Some unique characteristics observed in the response amplitude operator (RAO) and flow patterns around the boxes are explained using instantaneous velocity and vorticity diagrams. The response of various configurations under a nonlinear transient-focused wave is also examined. The first-order and second-order time series of wave elevation in the gap is obtained by using the four-phase combination method.

Simulation of reactive transport in porous media using Meshless Local Petrov Galerkin (MLPG) and combination of Meshless Weak and Strong (MWS) form methods

In this paper, two meshless methods, namely, a weak form Meshless Local Petrov Galerkin (MLPG) method, and Meshless Weak Strong (MWS) form method, obtained by combining MLPG with a strong form Radial Point Collocation Method (RPCM), are presented for simulation of advection-dispersion-reaction phenomena of the contaminants in the porous media. The first-order decay and sorption reactions are considered in this study. The Crank Nicolson scheme is applied for the time discretization. The weak form MLPG is a truly meshless and robust numerical technique, that can be applied to complex aquifer systems with derivative boundaries. However, in this method, the computational time is increased due to the integration, which is not essential for simple problems. Thus, the MLPG method is further coupled with a strong form RPCM with an aim of decreasing the background integration, by modelling only the nodes around the derivative boundaries using MLPG method and the other nodes by a direct RPCM which do not require integration. The proposed MWS model automatically converts into a complete RPCM model if there are no derivative boundaries. Thus, this model being both accurate and computationally efficient is suitable for simple and moderately complex aquifer systems and MLPG is the most stable and reliable method for modelling the most complex aquifer problems. Both the developed models are tested with available analytical solutions and applied for hypothetical case studies. The results prove the efficiency of the models and the applicability of each model is described in detail.

Meshless local Petrov-Galerkin method for 2D fractional Fokker-Planck equation involved with the ABC fractional derivative

This paper examines a time fractional version of the 2D Fokker-Planck equation involved with the Atangana-Baleanu-Caputo fractional derivative, under the Dirichlet boundary conditions. A fully discretization approach based on the meshless local Petrov-Galerkin method and ( 3 − β ) -order approximation is proposed for this equation. More precisely, we apply the meshless local Petrov-Galerkin method based on the Moving Kriging interpolation to discretize the space domain, and utilize the ( 3 − β ) -order approximation together with the θ -weighted finite difference method to discretize the temporal domain. By implementing this method, we get a solution for the problem by solving a system of algebraic equations. The validity of the method is investigated by solving four examples.

A strong-form meshfree computational method for plane elastostatic equations of anisotropic functionally graded materials via multiple-scale Pascal polynomials

A strong-form meshfree method is proposed for solving plane elastostatic equations of anisotropic functionally graded materials. Any general function may be the grading function and it is changing smoothly from location to location in the material. The proposed method is based on Pascal polynomial basis and multiple-scale technique and it is a genuinely meshfree method since no numerical integrations over domains and meshing processes are required for considered problems. Implementation of the proposed method is straightforward and the method gives very accurate results. Stability of the solutions are examined numerically in occurrence of random noise. Some certain test problems with known exact solutions are solved both on regular and irregular geometries. Acquired solutions by the suggested method are compared with the exact solutions as well as with solutions of some existing numerical techniques in literature, such as boundary element, meshless local Petrov–Galerkin and radial basis function based meshless methods, to show accuracy of the proposed method.

Three-dimensional coupling between Boussinesq (FEM) and Navier–Stokes (particle based) models for wave structure interaction

The paper presents coupling between a mesh-based finite-element model for Boussinesq equations (FEBOUSS Agarwal et al., 2022) with a meshless local Petrov–Galerkin model for the Navier–Stokes equations (MLPG_R Agarwal et al., 2021) in 3D. Boussinesq equation models are widely used for simulating wave-propagation over large domains with uneven topography using a 2D surface mesh. Mesh-less models inherently capture large free-surface deformations and have shown promise in simulating wave-structure interaction, run-up and breaking phenomenon. The hybrid approach in this paper assumes a 3D MLPG_R sub-domain surrounded by the 2D mesh of FEBOUSS. The coupling interface in MLPG_R consists of relaxation zones that can be placed along multiple boundaries of the sub-domain for exchanging particle velocity from FEBOUSS. This hybrid model is therefore capable of simulating directional waves, that has not been reported previously.The paper first presents the procedure for calculating the depth-resolved velocities in 3D from the Boussinesq model. The resultant velocities are compared against theory, experiments and other models. The following sections present the coupling algorithm along a single and multiple coupling interfaces in MLPG_R. Validation results for this hybrid model are provided using surface elevation and velocity measurements for regular waves, including directional cases. In general, the results from the hybrid model are reported to have marginal over-prediction of peaks compared to purely MLPG_R simulation. Finally, the interaction of a vertical cylinder with direction regular wave is simulated using the 3D hybrid model.

On the dynamics and stability of size-dependent symmetric FGM plates with electro-elastic coupling using meshless local Petrov-Galerkin method

In present study, the equations of motion of sandwich symmetric functionally graded size-dependent plates integrated with piezoelectric layers are derived using Hamilton’s variational principle (with modified Lagrangian) based on a higher-order nonlocal strain gradient theory in conjunction with Reddy’s third-order shear deformation plate kinematics. For the first time, a set of sixth-order partial differential equations are derived and solved by meshless local Petrov-Galerkin method to investigate the bifurcation buckling and free vibration of plates under electromechanical in-plane forces with considered nanostructure and various boundary conditions. The shifted and scaled moving least-squares approach is utilized for approximation of unknown variables, and the Gauss weight function is used as test function for deriving local weak form of governing equations. The Gauss-Legendre quadrature is employed for numerical evaluation of weak forms. A power-law distribution and a half cosine variation are applied to model the variation of materials properties and electric potentials, respectively. The obtained results are verified with well-known solutions in the literature. The effect of nonlocal parameter, material length scale parameter, power-law index, predefined electric field, axial compressive and tensile forces, volume ratio of FG core and piezoelectric layers, and stiffness of Pasternak foundation on the critical buckling loads and natural frequencies are presented and discussed comprehensively. The results obtained can be used in the analysis, design and control of composite nanostructures that are widely used in MEMS and NEMS devices.

Flexible GMRES solver for interpolating MLPG analysis of heat conduction

This article evaluates the performance of flexible generalized minimum residual (FGMRES) solver in meshless local Petrov–Galerkin (MLPG) analysis of heat conduction. Interpolating MLPG formulation has been used for the meshless discretization, resulting in an asymmetric algebraic equation system. FGMRES is a variant of GMRES, which allows the preconditioner to vary from step to step. Also, the GMRES solver can be used as a preconditioner. The FGMRES method has been used with Jacobi, SOR(k), and GMRES preconditioners for MLPG analysis of heat conduction in two-dimensional and three-dimensional geometries. The performance of the FGMRES solver is compared with the restarted GMRES solver (and BiCGSTAB solver in 3-D complex geometry). Jacobi and SOR(k) preconditioners are used in the restarted GMRES and BICGSTAB solvers. Results demonstrate that the FGMRES solver performs better than the GMRES and BiCGSTAB solvers, except for two-dimensional complex geometry.