scholarly journals A finite difference method for fractional diffusion equations with Neumann boundary conditions

2015 ◽  
Vol 13 (1) ◽  
Author(s):  
Béla J. Szekeres ◽  
Ferenc Izsák

AbstractA finite difference numerical method is investigated for fractional order diffusion problems in one space dimension. The basis of the mathematical model and the numerical approximation is an appropriate extension of the initial values, which incorporates homogeneous Dirichlet or Neumann type boundary conditions. The wellposedness of the obtained initial value problem is proved and it is pointed out that each extension is compatible with the original boundary conditions. Accordingly, a finite difference scheme is constructed for the Neumann problem using the shifted Grünwald–Letnikov approximation of the fractional order derivatives, which is based on infinite many basis points. The corresponding matrix is expressed in a closed form and the convergence of an appropriate implicit Euler scheme is proved.

2020 ◽  
Vol 0 (0) ◽  
Author(s):  
Mohammed Alaoui ◽  
Abdelkarim Hajjaj ◽  
Lahcen Maniar ◽  
Jawad Salhi

AbstractIn this paper, we study an inverse source problem for a degenerate and singular parabolic system where the boundary conditions are of Neumann type. We consider a problem with degenerate diffusion coefficients and singular lower-order terms, both vanishing at an interior point of the space domain. In particular, we address the question of well-posedness of the problem, and then we prove a stability estimate of Lipschitz type in determining the source term by data of only one component. Our method is based on Carleman estimates, cut-off procedures and a reflection technique.


2020 ◽  
Vol 54 (4) ◽  
pp. 1373-1413 ◽  
Author(s):  
Huaiqian You ◽  
XinYang Lu ◽  
Nathaniel Task ◽  
Yue Yu

In this paper we consider 2D nonlocal diffusion models with a finite nonlocal horizon parameter δ characterizing the range of nonlocal interactions, and consider the treatment of Neumann-like boundary conditions that have proven challenging for discretizations of nonlocal models. We propose a new generalization of classical local Neumann conditions by converting the local flux to a correction term in the nonlocal model, which provides an estimate for the nonlocal interactions of each point with points outside the domain. While existing 2D nonlocal flux boundary conditions have been shown to exhibit at most first order convergence to the local counter part as δ → 0, the proposed Neumann-type boundary formulation recovers the local case as O(δ2) in the L∞ (Ω) norm, which is optimal considering the O(δ2) convergence of the nonlocal equation to its local limit away from the boundary. We analyze the application of this new boundary treatment to the nonlocal diffusion problem, and present conditions under which the solution of the nonlocal boundary value problem converges to the solution of the corresponding local Neumann problem as the horizon is reduced. To demonstrate the applicability of this nonlocal flux boundary condition to more complicated scenarios, we extend the approach to less regular domains, numerically verifying that we preserve second-order convergence for non-convex domains with corners. Based on the new formulation for nonlocal boundary condition, we develop an asymptotically compatible meshfree discretization, obtaining a solution to the nonlocal diffusion equation with mixed boundary conditions that converges with O(δ2) convergence.


2021 ◽  
Vol 10 (8) ◽  
pp. 3013-3022
Author(s):  
C.A. Gomez ◽  
J.A. Caicedo

In this work, we consider the rescaled nonlocal diffusion problem with Neumann Boundary Conditions \[ \begin{cases} u_t^{\epsilon}(x,t)=\displaystyle\frac{1}{\epsilon^2} \int_{\Omega}J_{\epsilon}(x-y)(u^\epsilon(y,t)-u^\epsilon(x,t))dy\\ \qquad \qquad+\displaystyle\frac{1}{\epsilon}\int_{\partial \Omega}G_\epsilon(x-y)g(y,t)dS_y,\\ u^\epsilon(x,0)=u_0(x), \end{cases} \] where $\Omega\subset\mathbb{R}^{N}$ is a bounded, connected and smooth domain, $g$ a positive continuous function, $J_\epsilon(z)=C_1\frac{1}{\epsilon^N}J(\frac{z}{\epsilon}), G_\epsilon(x)=C_1\frac{1}{\epsilon^N}G(\frac{x}{\epsilon}),$ $J$ and $G$ well defined kernels, $C_1$ a normalization constant. The solutions of this model have been used without prove to approximate the solutions of a family of nonlocal diffusion problems to solutions of the respective analogous local problem. We prove existence and uniqueness of the solutions for this problem by using the Banach Fixed Point Theorem. Finally, some conclusions are given.


Geophysics ◽  
1978 ◽  
Vol 43 (6) ◽  
pp. 1099-1110 ◽  
Author(s):  
Albert C. Reynolds

Many finite difference models in use for generating synthetic seismograms produce unwanted reflections from the edges of the model due to the use of Dirichlet or Neumann boundary conditions. In this paper we develop boundary conditions which greatly reduce this edge reflection. A reflection coefficient analysis is given which indicates that, for the specified boundary conditions, smaller reflection coefficients than those obtained for Dirichlet or Neumann boundary conditions are obtained. Numerical calculations support this conclusion.


2004 ◽  
Vol 2004 (9) ◽  
pp. 777-792 ◽  
Author(s):  
Jiří Benedikt

We are interested in a nonlinear boundary value problem for(|u″|p−2u″)′​′=λ|u|p−2uin[0,1],p>1, with Dirichlet and Neumann boundary conditions. We prove that eigenvalues of the Dirichlet problem are positive, simple, and isolated, and form an increasing unbounded sequence. An eigenfunction, corresponding to thenth eigenvalue, has preciselyn−1zero points in(0,1). Eigenvalues of the Neumann problem are nonnegative and isolated,0is an eigenvalue which is not simple, and the positive eigenvalues are simple and they form an increasing unbounded sequence. An eigenfunction, corresponding to thenth positive eigenvalue, has preciselyn+1zero points in(0,1).


Sign in / Sign up

Export Citation Format

Share Document