An asymptotically compatible formulation for local-to-nonlocal coupling problems without overlapping regions

Abstract

In this paper we design and analyze an explicit partitioned procedure for a 2D dynamic local-to-nonlocal (LtN) coupling problem, based on a new nonlocal Robin-type transmission condition. The nonlocal subproblem is modeled by the nonlocal heat equation with a finite horizon parameter δ characterizing the range of nonlocal interactions, and the local subproblem is described by the classical heat equation. We consider a heterogeneous system where the local and nonlocal subproblems present different physical properties, and employ no overlapping region between the two subdomains. We first propose a new generalization of classical local Neumann-type condition by converting the local flux to a correction term in the nonlocal model, and show that the proposed Neumann-type boundary formulation recovers the local case as O(δ 2 ) in the L ∞ norm. We then extend the nonlocal Neumann-type boundary condition to a Robin-type boundary condition, and develop a local-to-nonlocal coupling formulation with Robin-Dirichlet transmission conditions. To stabilize the explicit coupling procedure and to achieve asymptotic compatibility, the choice of the coefficient in the Robin condition is obtained via amplification factor analysis for the discretized system with coarse grids. Employing a high-order meshfree discretization method in the nonlocal solver and a linear finite element method in the local solver, the selection of optimal Robin coefficients are verified with numerical experiments on heterogeneous and complicated domains. With the developed optimal coupling strategy, we numerically demonstrate the coupling framework's asymptotic convergence to the local limit with an O(δ) = O(h) rate, when there is a fixed ratio between the horizon size δ and the spatial discretization size h.

Numerical methods for nonlocal and fractional models

“…We note that there have been other works concerning AC schemes for nonlocal models based on the strong form. For example, Du, Zhang and Zheng (2018 d ), You, Yu and Kamensky (2019) and Tao et al (2019) have studied AC schemes for coupled local and nonlocal models. Du, Han, Zheng and Zhang (2018 a ), Zhang, Yang, Zhang and Du (2017) and Du et al (2018 d ) presented implementation techniques and numerical experiments of AC schemes for problems defined in infinite domains via the development of nonlocal artificial boundary conditions.…”

A data-driven peridynamic continuum model for upscaling molecular dynamics

“…In this work, without loss of generality, we consider the Dirichlet condition B I = I, where I is the identity operator. Other types of conditions, e.g., Neumann [40,42,46,47], Robin [43,48] or periodic [32], are also compatible with our learning algorithm.…”

“…The LPS model has known well-posedness properties under certain assumptions [24]. This section summarizes the mathematical formulation for the LPS model and illustrates a meshfree discretization [39][40][41][42][43][44].…”

Hybrid analysis and modeling, eclecticism, and multifidelity computing toward digital twin revolution

“…This can be effectively achieved using a domain decomposition approach [49,63‐65,85,100,108,125,143,190,251,261,322,331,332], a powerful concept often used to solve partial differential equations (PDEs) with complex geometries. Here, we can interpret the term domain decomposition in a broad sense such that it also includes the simultaneous advances in the multiscale and multiphysics problems [18,61,91,103,120,148,149,211,252,311,349,357]. Many approaches in the field of domain decomposition are well‐suited for heterogeneous problems, which have dissimilar governing equations or models (a.k.a.…”