Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConvergence Analysis and Numerical Studies for Linearly Elastic Peridynamics with Dirichlet-Type Boundary Conditions
305
0
0.0
(
0
)
Added by
Yue Yu
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The nonlocal models of peridynamics have successfully predicted fractures and deformations for a variety of materials. In contrast to local mechanics, peridynamic boundary conditions must be defined on a finite volume region outside the body. Therefore, theoretical and numerical challenges arise in order to properly formulate Dirichlet-type nonlocal boundary conditions, while connecting them to the local counterparts. While a careless imposition of local boundary conditions leads to a smaller effective material stiffness close to the boundary and an artificial softening of the material, several strategies were proposed to avoid this unphysical surface effect. In this work, we study convergence of solutions to nonlocal state-based linear elastic model to their local counterparts as the interaction horizon vanishes, under different formulations and smoothness assumptions for nonlocal Dirichlet-type boundary conditions. Our results provide explicit rates of convergence that are sensitive to the compatibility of the nonlocal boundary data and the extension of the solution for the local model. In particular, under appropriate assumptions, constant extensions yield $frac{1}{2}$ order convergence rates and linear extensions yield $frac{3}{2}$ order convergence rates. With smooth extensions, these rates are improved to quadratic convergence. We illustrate the theory for any dimension $dgeq 2$ and numerically verify the convergence rates with a number of two dimensional benchmarks, including linear patch tests, manufactured solutions, and domains with curvilinear surfaces. Numerical results show a first order convergence for constant extensions and second order convergence for linear extensions, which suggests a possible room of improvement in the future convergence analysis.
rate research
Read More
We present an approach to handle Dirichlet type nonlocal boundary conditions for nonlocal diffusion models with a finite range of nonlocal interactions. Our approach utilizes a linear extrapolation of prescribed boundary data. A novelty is, instead of using local gradients of the boundary data that are not available a priori, we incorporate nonlocal gradient operators into the formulation to generalize the finite differences-based methods which are pervasive in literature; our particular choice of the nonlocal gradient operators is based on the interplay between a constant kernel function and the geometry of nonlocal interaction neighborhoods. Such an approach can be potentially useful to address similar issues in peridynamics, smoothed particle hydrodynamics and other nonlocal models. We first show the well-posedness of the newly formulated nonlocal problems and then analyze their asymptotic convergence to the local limit as the nonlocality parameter shrinks to zero. We justify the second order localization rate, which is the optimal order attainable in the absence of physical boundaries.
This paper is concerned with the multiplicity results to a class of $p$-Kirchhoff type elliptic equation with the homogeneous Neumann boundary conditions by an abstract linking lemma due to Br{e}zis and Nirenberg. We obtain the twofold results in subcritical and critical cases, which is a meaningful addition and completeness to the known results about Kirchhoff equation.
This is the first part of our study of inertial manifolds for the system of 1D reaction-diffusion-advection equations which is devoted to the case of Dirichlet or Neumann boundary conditions. Although this problem does not initially possess the spectral gap property, it is shown that this property is satisfied after the proper non-local change of the dependent variable. The case of periodic boundary conditions where the situation is principally different and the inertial manifold may not exist is considered in the second part of our study.
A hybrid surface integral equation partial differential equation (SIE-PDE) formulation without the boundary condition requirement is proposed to solve the electromagnetic problems. In the proposed formulation, the computational domain is decomposed into two emph{overlapping} domains: the SIE and PDE domains. In the SIE domain, complex structures with piecewise homogeneous media, e.g., highly conductive media, are included. An equivalent model for those structures is constructed through replacing them by the background medium and introducing a surface equivalent electric current density on an enclosed boundary to represent their electromagnetic effects. The remaining computational domain and homogeneous background medium replaced domain consist of the PDE domain, in which inhomogeneous or non-isotropic media are included. Through combining the surface equivalent electric current density and the inhomogeneous Helmholtz equation, a hybrid SIE-PDE formulation is derived. Unlike other hybrid formulations, where the transmission condition is usually used, no boundary conditions are required in the proposed SIE-PDE formulation, and it is mathematically equivalent to the original physical model. Through careful construction of basis functions to expand electric fields and the equivalent current density, the discretized formulation is compatible on the interface of the SIE and PDE domain. Finally, its accuracy and efficiency are validated through two numerical examples. Results show that the proposed SIE-PDE formulation can obtain accurate results including both near and far fields, and significant performance improvements in terms of CPU time and memory consumption compared with the FEM are achieved.
We introduce a technique to automatically convert local boundary conditions into nonlocal volume constraints for nonlocal Poissons and peridynamic models. The proposed strategy is based on the approximation of nonlocal Dirichlet or Neumann data with a local solution obtained by using available boundary, local data. The corresponding nonlocal solution converges quadratically to the local solution as the nonlocal horizon vanishes, making the proposed technique asymptotically compatible. The proposed conversion method does not have any geometry or dimensionality constraints and its computational cost is negligible, compared to the numerical solution of the nonlocal equation. The consistency of the method and its quadratic convergence with respect to the horizon is illustrated by several two-dimensional numerical experiments conducted by meshfree discretization for both the Poissons problem and the linear peridynamic solid model.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
