Search for: All records

We use the polyharmonic extension approach to develop a numerical technique for discretizing higher-order powers of the spectral fractional Laplacian $$(-\Delta)^s$$ with $$s \in (1,2)$$. 

We present a technique for approximating solutions to the spectral fractional Lapla- cian, which is based on the Caffarelli-Silvestre extension and diagonalization. Our scheme uses the analytic solution to the associated eigenvalue problem in the extended dimension. We show its relation to a quadrature scheme. Numerical examples demonstrate the performance of the method 

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. 

We study optimal design problems where the design corresponds to a coefficient in the principal part of the state equation. The state equation, in addition, is parameter dependent, and we allow it to change type in the limit of this (modeling) parameter. We develop a framework that guarantees asymptotic compatibility, that is unconditional convergence with respect to modeling and discretization parameters to the solution of the corresponding limiting problems. This framework is then applied to two distinct classes of problems where the modeling parameter represents the degree of nonlocality. Specifically, we show unconditional convergence of optimal design problems when the state equation is either a scalar-valued fractional equation, or a strongly coupled system of nonlocal equations derived from the bond-based model of peridynamics. 

Abstract We consider the interaction between a poroelastic structure, described using the Biot model in primal form, and a free-flowing fluid, modelled with the time-dependent incompressible Stokes equations. We propose a diffuse interface model in which a phase field function is used to write each integral in the weak formulation of the coupled problem on the entire domain containing both the Stokes and Biot regions. The phase field function continuously transitions from one to zero over a diffuse region of width $$\mathcal{O}(\varepsilon)$$ around the interface; this allows the weak forms to be integrated uniformly across the domain, and obviates tracking the subdomains or the interface between them. We prove convergence in weighted norms of a finite element discretization of the diffuse interface model to the continuous diffuse model; here the weight is a power of the distance to the diffuse interface. We, in turn, prove convergence of the continuous diffuse model to the standard, sharp interface, model. Numerical examples verify the proven error estimates, and illustrate application of the method to fluid flow through a complex network, describing blood circulation in the circle of Willis. 

We prove the well posedness in weighted Sobolev spaces of certain linear and nonlinear elliptic boundary value problems posed on convex domains and under singular forcing. It is assumed that the weights belong to the Muckenhoupt class with ). We also propose and analyze a convergent finite element discretization for the nonlinear elliptic boundary value problems mentioned above. As an instrumental result, we prove that the discretization of certain linear problems are well posed in weighted spaces. 

Abstract We propose a monotone and consistent numerical scheme for the approximation of the Dirichlet problem for the normalized infinity Laplacian, which could be related to the family of the so-called two-scale methods. We show that this method is convergent and prove rates of convergence. These rates depend not only on the regularity of the solution, but also on whether or not the right-hand side vanishes. Some extensions to this approach, like obstacle problems and symmetric Finsler norms, are also considered. 

We study a non-local optimal control problem involving a linear, bond-based peridynamics model. In addition to existence and uniqueness of solutions to our problem, we investigate their behavior as the horizon parameter 𝛿, which controls the degree of nonlocality, approaches zero. We then study a finite element-based discretization of this problem, its convergence, and the so-called asymptotic compatibility as the discretization parameter h and the horizon parameter 𝛿 tend to zero simultaneously. 

We consider the interaction between a free flowing fluid and a porous medium flow, where the free flowing fluid is described using the time dependent Stokes equations, and the porous medium flow is described using Darcy’s law in the primal formulation. To solve this problem numerically, we use a diffuse interface approach, where the weak form of the coupled problem is written on an extended domain which contains both Stokes and Darcy regions. This is achieved using a phase-field function which equals one in the Stokes region and zero in the Darcy region, and smoothly transitions between these two values on a diffuse region of width (ϵ) around the Stokes-Darcy interface. We prove convergence of the diffuse interface formulation to the standard, sharp interface formulation, and derive rates of convergence. This is performed by deriving a priori error estimates for discretizations of the diffuse interface method, and by analyzing the modeling error of the diffuse interface approach at the continuous level. The convergence rates are also shown computationally in a numerical example.