Search

Many nonlocal models have adopted Euclidean balls as the nonlocal interaction neighborhoods. When solving them numerically, it is sometimes convenient to adopt polygonal approximations of such balls. A crucial question is, to what extent such approximations affect the nonlocal operators and the corresponding solutions. While recent works have analyzed this issue for a fixed horizon parameter, the question remains open in the case of a small or vanishing horizon parameter, which happens often in many practical applications and has significant impact on the reliability and robustness of nonlocal modeling and simulations. In this work, we are interested in addressing this issue and establishing the convergence of the nonlocal solutions associated with polygonally approximated interaction neighborhoods to the local limit of the original nonlocal solutions. Our finding reveals that the new nonlocal solution does not converge to the correct local limit when the number of sides of polygons is uniformly bounded. On the other hand, if the number of sides tends to infinity, the desired convergence can be established. These results may be used to guide future computational studies of nonlocal models.

The piecewise quadratic polynomial collocation is used to approximate the nonlocal model, which generally obtain the {\em nonsymmetric indefinite system} [Chen et al., IMA J. Numer. Anal., (2021)]. In this case, the discrete maximum principle is not satisfied, which might be trickier for the stability analysis of the high-order numerical schemes [D'Elia et al., Acta Numer., (2020); Leng et al., SIAM J. Numer. Anal., (2021)]. Here, we present the modified (shifted-symmetric) piecewise quadratic polynomial collocation for solving the linear nonlocal diffusion model, which has the {\em symmetric positive definite system} and satisfies the discrete maximum principle. Using Faulhaber's formula and Riemann zeta function, the perturbation error for symmetric positive definite system and nonsymmetric indefinite systems are given. Then the detailed proof of the convergence analysis for the nonlocal models with the general horizon parameter $\delta=\mathcal{O}\left(h^\beta\right)$, $\beta\geq0$ are provided. More concretely, the global error is $\mathcal{O}\left(h^{\min\left\{2,1+\beta\right\}}\right)$ if $\delta$ is not set as a grid point, but it shall recover $\mathcal{O}\left(h^{\max\left\{2,4-2\beta\right\}}\right)$ when $\delta$ is set as a grid point. We also prove that the shifted-symmetric scheme is asymptotically compatible, which has the global error $\mathcal{O}\left(h^{\min\left\{2,2\beta\right\}}\right)$ as $\delta,h\rightarrow 0$. The numerical experiments (including two-dimensional case) are performed to verify the convergence., Comment: 23 pages

Numerical solution of nonlocal constrained value problems with integrable kernels are considered. These nonlocal problems arise in nonlocal mechanics and nonlocal diffusion. The structure of the true solution to the problem is analyzed first. The analysis leads naturally to a new kind of discontinuous Galerkin method that efficiently solve the problem numerically. This method is shown to be asymptotically compatible. Moreover, it has optimal convergence rate for one dimensional case under very weak assumptions, and almost optimal convergence rate for two dimensional case under mild assumptions.

We propose a numerical strategy to generate the anisotropic meshes and select the appropriate stabilized parameters simultaneously for two dimensional convection-dominated convection-diffusion equations by stabilized continuous linear finite elements. Since the discretized error in a suitable norm can be bounded by the sum of interpolation error and its variants in different norms, we replace them by some terms which contain the Hessian matrix of the true solution, convective fields, and the geometric properties such as directed edges and the area of the triangle. Based on this observation, the shape, size and equidistribution requirements are used to derive the corresponding metric tensor and the stabilized parameters. It is easily found from our derivation that the optimal stabilized parameter is coupled with the optimal metric tensor on each element. Some numerical results are also provided to validate the stability and efficiency of the proposed numerical strategy., Comment: 24 pages, 4 figures

A uniform strategy to derive metric tensors in two spatial dimension for interpolation errors and their gradients in $L^p$ norm is presented. It generates anisotropic adaptive meshes as quasi-uniform ones in corresponding metric space, with the metric tensor being computed based on a posteriori error estimates in different norms. Numerical results show that the corresponding convergence rates are always optimal., Comment: 19 pages, 24 figures

A new anisotropic mesh adaptation strategy for finite element solution of elliptic differential equations is presented. It generates anisotropic adaptive meshes as quasi-uniform ones in some metric space, with the metric tensor being computed based on a posteriori error estimates proposed in \cite{YinXie}. The new metric tensor explores more comprehensive information of anisotropy for the true solution than those existing ones. Numerical results show that this approach can be successfully applied to deal with poisson and steady convection-dominated problems. The superior accuracy and efficiency of the new metric tensor to others is illustrated on various numerical examples of complex two-dimensional simulations., Comment: 18 pages, 18 figures

In this paper, we give a new type of a posteriori error estimators suitable for moving finite element methods under anisotropic meshes for general second-order elliptic problems. The computation of estimators is simple once corresponding Hessian matrix is recovered. Wonderful efficiency indices are shown in numerical experiments., Comment: 16 pages, 16 figures