An official website of the United States government
We study the regularity of the solution to an obstacle problem for a class of integro–differential operators. The differential part is a second order elliptic operator, whereas the nonlocal part is given by the integral fractional Laplacian. The obtained smoothness is then used to design and analyze a finite element scheme.
more »
« less
Estimates for Dirichlet-to-Neumann Maps as Integro-differential Operators
Some linear integro-differential operators have old and classical representations as the Dirichlet-to-Neumann operators for linear elliptic equations, such as the 1/2-Laplacian or the generator of the boundary process of a reflected diffusion. In this work, we make some extensions of this theory to the case of a nonlinear Dirichlet-to-Neumann mapping that is constructed using a solution to a fully nonlinear elliptic equation in a given domain, mapping Dirichlet data to its normal derivative of the resulting solution. Here we begin the process of giving detailed information about the Lévy measures that will result from the integro-differential representation of the Dirichlet-to-Neumann mapping. We provide new results about both linear and nonlinear Dirichlet-to-Neumann mappings. Information about the Lévy measures is important if one hopes to use recent advancements of the integro-differential theory to study problems involving Dirichlet-to-Neumann mappings.
more »
« less
- PAR ID:
- 10097955
- Date Published:
- Journal Name:
- Potential analysis
- ISSN:
- 1572-929X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Grid-free Monte Carlo methods based on thewalk on spheres (WoS)algorithm solve fundamental partial differential equations (PDEs) like the Poisson equation without discretizing the problem domain or approximating functions in a finite basis. Such methods hence avoid aliasing in the solution, and evade the many challenges of mesh generation. Yet for problems with complex geometry, practical grid-free methods have been largely limited to basic Dirichlet boundary conditions. We introduce thewalk on stars (WoSt)algorithm, which solves linear elliptic PDEs with arbitrary mixed Neumann and Dirichlet boundary conditions. The key insight is that one can efficiently simulate reflecting Brownian motion (which models Neumann conditions) by replacing the balls used by WoS withstar-shapeddomains. We identify such domains via the closest point on the visibility silhouette, by simply augmenting a standard bounding volume hierarchy with normal information. Overall, WoSt is an easy modification of WoS, and retains the many attractive features of grid-free Monte Carlo methods such as progressive and view-dependent evaluation, trivial parallelization, and sublinear scaling to increasing geometric detail.more » « less
-
Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.more » « less
-
We prove higher Sobolev regularity for bounded weak solutions to a class of nonlinear nonlocal integro-differential equations. The leading operator exhibits nonuniform growth, switching between two different fractional elliptic "phases" that are determined by the zero set of a modulating coefficient. Solutions are shown to improve both in integrability and differentiability. These results apply to operators with rough kernels and modulating coefficients. To obtain these results we adapt a particular fractional version of the Gehring lemma developed by Kuusi, Mingione, and Sire in their work "Nonlocal self-improving properties" Analysis & PDE, 8(1):57–114 for the specific nonlinear setting under investigation in this manuscript.more » « less
-
We develop a monotone, two-scale discretization for a class of integro-differential operators of order2s,s \in (0,1). We apply it to develop numerical schemes, and derive pointwise convergence rates for linear and obstacle problems governed by such operators. As applications of the monotonicity, we provide error estimates for free boundaries and a convergent numerical scheme for a concave fully nonlinear, nonlocal, problem.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1007/s11118-019-09776-w
