1,061 research outputs found
An advanced meshless method for time fractional diffusion equation
Recently, because of the new developments in sustainable engineering and renewable energy, which are usually governed by a series of fractional partial differential equations (FPDEs), the numerical modelling and simulation for fractional calculus are attracting more and more attention from researchers. The current dominant numerical method for modeling FPDE is Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings including difficulty in simulation of problems with the complex problem domain and in using irregularly distributed nodes. Because of its distinguished advantages, the meshless method has good potential in simulation of FPDEs. This paper aims to develop an implicit meshless collocation technique for FPDE. The discrete system of FPDEs is obtained by using the meshless shape functions and the meshless collocation formulation. The stability and convergence of this meshless approach are investigated theoretically and numerically. The numerical examples with regular and irregular nodal distributions are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling and simulation of fractional partial differential equations
Differential quadrature method for space-fractional diffusion equations on 2D irregular domains
In mathematical physics, the space-fractional diffusion equations are of
particular interest in the studies of physical phenomena modelled by L\'{e}vy
processes, which are sometimes called super-diffusion equations. In this
article, we develop the differential quadrature (DQ) methods for solving the 2D
space-fractional diffusion equations on irregular domains. The methods in
presence reduce the original equation into a set of ordinary differential
equations (ODEs) by introducing valid DQ formulations to fractional directional
derivatives based on the functional values at scattered nodal points on problem
domain. The required weighted coefficients are calculated by using radial basis
functions (RBFs) as trial functions, and the resultant ODEs are discretized by
the Crank-Nicolson scheme. The main advantages of our methods lie in their
flexibility and applicability to arbitrary domains. A series of illustrated
examples are finally provided to support these points.Comment: 25 pages, 25 figures, 7 table
Comparison of a Material Point Method and a Galerkin meshfree method for the simulation of cohesive-frictional materials
The simulation of large deformation problems, involving complex history-dependent constitutive laws, is of paramount importance in several engineering fields. Particular attention has to be paid to the choice of a suitable numerical technique such that reliable results can be obtained. In this paper, a Material Point Method (MPM) and a Galerkin Meshfree Method (GMM) are presented and verified against classical benchmarks in solid mechanics. The aim is to demonstrate the good behavior of the methods in the simulation of cohesive-frictional materials, both in static and dynamic regimes and in problems dealing with large deformations. The vast majority of MPM techniques in the literature are based on some sort of explicit time integration. The techniques proposed in the current work, on the contrary, are based on implicit approaches, which can also be easily adapted to the simulation of static cases. The two methods are presented so as to highlight the similarities to rather than the differences fromPeer ReviewedPostprint (published version
Addressing Integration Error for Polygonal Finite Elements Through Polynomial Projections: A Patch Test Connection
Polygonal finite elements generally do not pass the patch test as a result of
quadrature error in the evaluation of weak form integrals. In this work, we
examine the consequences of lack of polynomial consistency and show that it can
lead to a deterioration of convergence of the finite element solutions. We
propose a general remedy, inspired by techniques in the recent literature of
mimetic finite differences, for restoring consistency and thereby ensuring the
satisfaction of the patch test and recovering optimal rates of convergence. The
proposed approach, based on polynomial projections of the basis functions,
allows for the use of moderate number of integration points and brings the
computational cost of polygonal finite elements closer to that of the commonly
used linear triangles and bilinear quadrilaterals. Numerical studies of a
two-dimensional scalar diffusion problem accompany the theoretical
considerations
- …