PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

OPTIMAL ERROR ESTIMATES OF THE PARTITION OF UNITY METHOD WITH LOCAL POLYNOMIAL APPROXIMATION SPACES

In this paper, we provide a theoretical method（PUFEM）, which belongs to the analysis of the partition of unity finite element family of meshfree methods. The usual error analysis only shows the order of error estimate to the same as the local approximations[12]. Using standard linear finite element base functions as partition of unity and polynomials as local approximation space, in l-d case, we derive optimal order error estimates for PUFEM interpolants. Our analysis show that the error estimate is of one order higher than the local approximations. The interpolation error estimates yield optimal error estimates for PUFEM solutions of elliptic boundary value problems.

Partition of unity interpolation on multivariate convex domains

In this paper,we present an algorithm for multivariate interpolation of scattered data sets lying in convex domainsΩ⊆R^(N),for any N≥2.To organize the points in a multidimensional space,we build a kd-tree space-partitioning data structure,which is used to efficiently apply a partition of unity interpolant.This global scheme is combined with local radial basis function(RBF)approximants and compactly supported weight functions.A detailed description of the algorithm for convex domains and a complexity analysis of the computational procedures are also considered.Several numerical experiments show the performances of the interpolation algorithm on various sets of Halton data points contained inΩ,whereΩcan be any convex domain,like a 2D polygon or a 3D polyhedron.Finally,an application to topographical data contained in a pentagonal domain is presented.

CONVERGENCE ANALYSIS OF SOME FINITE ELEMENT PARALLEL ALGORITHMS FOR THE STATIONARY INCOMPRESSIBLE MHD EQUATIONS

By combination of iteration methods with the partition of unity method(PUM),some finite element parallel algorithms for the stationary incompressible magnetohydrodynamics(MHD)with different physical parameters are presented and analyzed.These algorithms are highly efficient.At first,a global solution is obtained on a coarse grid for all approaches by one of the iteration methods.By parallelized residual schemes,local corrected solutions are calculated on finer meshes with overlapping sub-domains.The subdomains can be achieved flexibly by a class of PUM.The proposed algorithm is proved to be uniformly stable and convergent.Finally,one numerical example is presented to confirm the theoretical findings.

A high order numerical manifold method and its application to linear elastic continuous and fracture problems

The numerical manifold method(NMM) is a partition of unity(PU) based method. For the purpose of obtaining better accuracy with the same mesh, high order global approximation can be adopted by increasing the order of local approximations(LAs). This,however, will cause the "linear dependence"(LD) issue, where the global matrix is rank deficient even after sufficient constraints are enforced. In this paper, through quadrilateral mesh to form the mathematical cover, a high order numerical manifold method called Quad4-COLS(NMM) is developed, where the constrained and orthonormalized least-squares method(CO-LS) is used to construct the LAs. The developed Quad4-COLS(NMM) does not need extra nodes or DOFs to construct high order global approximations, while is free from the LD issue. Nine numerical tests including five tests for linear elastic continuous problems and four tests for linear elastic fracture problems are carried out to validate the accuracy and robustness of the proposed Quad4-COLS(NMM).