This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

On discrete boundary value problems with nonlocal conditions in a quarter-plane

Vladimir Vasilyev Chair of Applied Mathematics and Computer Modeling
Belgorod State National Research University
Pobedy street 85, Belgorod 308015, Russia
[email protected]
 and  Anastasia Mashinets Chair of Applied Mathematics and Computer Modeling
Belgorod State National Research University
Pobedy street 85, Belgorod 308015, Russia
[email protected]
Abstract.

We consider discrete analogue of model pseudo-differential equations in discrete plane sector using discrete variant of Sobolev–Slobodetskii spaces. Starting from the concept of wave factorization for elliptic periodic symbol we describe solvability conditions for the equations and corresponding discrete boundary value problems. We give also a comparison between discrete and continuous solutions in appropriate discrete normed spaces.

Key words and phrases:
elliptic symbol, invertibility, two-dimensional periodic Riemann problem, digital pseudo-differential operator, discrete equation, periodic wave factorization
2010 Mathematics Subject Classification:
Primary: 35S15; Secondary: 65T50

1. Introduction

A theory of pseudo-differential operators and equations [9, 16, 17, 6] not so long history than other mathematical subjects of analysis. Nevertheless, these operators and related boundary value problems widely arise in a lot of applied problems in physics and technique (see, for example [14] and references therein). Discrete aspects of the theory are reflected in mathematical papers more weak [11, 2] although these studies are closely related to the theory of Fourier series [5]. In our opinion the discrete theory is very important since it permits to use computer calculations to solve concrete applied problems.

We interested in studying discrete pseudo-differential equations and their solvability in appropriate discrete functional spaces. There are certain approaches to studying discrete boundary value problems for partial differential equations including finite difference method [13, 12]. But these approaches are not applicable to studying discrete boundary value problems for elliptic pseudo-differential equations. According to this statement the first author with colleagues has started to develop discrete theory for elliptic pseudo-differential equations [21, 15]. This is main motivation, and we have started from certain canonical domains. First considerations were related to discrete mm-dimensional space and half-space, and here we consider discrete quadrant.

We consider a special type of boundary conditions, namely integral conditions on a boundary. These conditions are nonlocal, and it seems, such conditions are artificial. But there are a lot of applied problems for partial differential equations with such boundary conditions [1, 4, 3], therefore it is natural way. Moreover, these conditions appear in a natural way to determine arbitrary functions in a general solution of an elliptic pseudo-differential equation.

2. Cones, periodic symbols, digital operators and equations

2.1. Discrete spaces and transforms

Let 2\mathbb{Z}^{2} be an integer lattice in a plane. Let K={x2:x=(x1,x2),x1>0,x2>0}K=\{x\in\mathbb{R}^{2}:x=(x_{1},x_{2}),x_{1}>0,x_{2}>0\} be a quadrant, Kd=h2K,h>0K_{d}=h\mathbb{Z}^{2}\cap K,h>0. We consider functions of discrete variable ud(x~),x~=(x~1,x~2)h2u_{d}(\tilde{x}),\tilde{x}=(\tilde{x}_{1},\tilde{x}_{2})\in h\mathbb{Z}^{2}.

Let us denote 𝕋2=[π,π]2,=h1\mathbb{T}^{2}=[-\pi,\pi]^{2},\hbar=h^{-1}. We consider functions defined in 𝕋2\hbar\mathbb{T}^{2} as periodic functions defined in 2\mathbb{R}^{2} with basic square of periods 𝕋2\hbar\mathbb{T}^{2}.

One can define the discrete Fourier transform for the function udu_{d}

(Fdud)(ξ)u~d(ξ)=x~h2eix~ξud(x~)h2,ξ𝕋2,(F_{d}u_{d})(\xi)\equiv\tilde{u}_{d}(\xi)=\sum\limits_{\tilde{x}\in h\mathbb{Z}^{2}}e^{-i\tilde{x}\cdot\xi}u_{d}(\tilde{x})h^{2},\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \xi\in\hbar\mathbb{T}^{2},

if the latter series converges, and the function u~d(ξ)\tilde{u}_{d}(\xi) is a periodic function in 2\mathbb{R}^{2} with basic square of periods 𝕋2\hbar\mathbb{T}^{2}. Such discrete Fourier transform preserves all properties of integral Fourier transform, and the inverse discrete Fourier transform looks as follows

(Fd1u~d)(x~)=1(2π)2𝕋2eix~ξu~d(ξ)𝑑ξ,x~h2.(F_{d}^{-1}\tilde{u}_{d})(\tilde{x})=\frac{1}{(2\pi)^{2}}\int\limits_{\hbar\mathbb{T}^{2}}e^{i\tilde{x}\cdot\xi}\tilde{u}_{d}(\xi)d\xi,\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{x}\in h\mathbb{Z}^{2}.

The discrete Fourier transform gives one-to-one correspondence between spaces L2(h2)L_{2}(h\mathbb{Z}^{2}) and L2(𝕋2)L_{2}(\hbar\mathbb{T}^{2}) with norms

ud2=(x~h2|ud(x~)|2h2)1/2,u~d2=(ξ𝕋2|u~d(ξ)|2𝑑ξ)1/2.||u_{d}||_{2}=\left(\sum\limits_{\tilde{x}\in h\mathbb{Z}^{2}}|u_{d}(\tilde{x})|^{2}h^{2}\right)^{1/2},\quad||\tilde{u}_{d}||_{2}=\left(\int\limits_{\xi\in\hbar\mathbb{T}^{2}}|\tilde{u}_{d}(\xi)|^{2}d\xi\right)^{1/2}.

We need more general discrete functional spaces and we introduce such spaces using divided differences [13].

The divided differences of first order look as follows

(Δ1(1)ud)(x~)=h1(ud(x~1+h,x~2)ud(x~1,x~2)),(\Delta_{1}^{(1)}u_{d})(\tilde{x})=h^{-1}(u_{d}(\tilde{x}_{1}+h,\tilde{x}_{2})-u_{d}(\tilde{x}_{1},\tilde{x}_{2})),
(Δ2(1)ud)(x~)=h1(ud(x~1,x~2+h)ud(x~1,x~2)),(\Delta_{2}^{(1)}u_{d})(\tilde{x})=h^{-1}(u_{d}(\tilde{x}_{1},\tilde{x}_{2}+h)-u_{d}(\tilde{x}_{1},\tilde{x}_{2})),

and their discrete Fourier transforms are given by formulas

(Δk(1)ud)~(ξ)=h1(eihξk1)u~d(ξ),k=1,2.\widetilde{(\Delta_{k}^{(1)}u_{d})}(\xi)=h^{-1}(e^{-ih\cdot\xi_{k}}-1)\tilde{u}_{d}(\xi),k=1,2.

The divided difference of second order is a divided difference of first order from divided difference of first order

(Δ1(2)ud)(x~)=h2(ud(x~1+2h,x~2)2ud(x~1+h,x~2)+ud(x~1,x~2)),(\Delta_{1}^{(2)}u_{d})(\tilde{x})=h^{-2}(u_{d}(\tilde{x}_{1}+2h,\tilde{x}_{2})-2u_{d}(\tilde{x}_{1}+h,\tilde{x}_{2})+u_{d}(\tilde{x}_{1},\tilde{x}_{2})),
(Δ2(2)ud)(x~)=h2(ud(x~1,x~2+2h)2ud(x~1,x~2+h)+ud(x~1,x~2)),(\Delta_{2}^{(2)}u_{d})(\tilde{x})=h^{-2}(u_{d}(\tilde{x}_{1},\tilde{x}_{2}+2h)-2u_{d}(\tilde{x}_{1},\tilde{x}_{2}+h)+u_{d}(\tilde{x}_{1},\tilde{x}_{2})),

with the Fourier transform

(Δk(2)ud)~(ξ)=h2(eihξk1)2u~d(ξ),k=1,2.\widetilde{(\Delta_{k}^{(2)}u_{d})}(\xi)=h^{-2}(e^{-ih\cdot\xi_{k}}-1)^{2}\tilde{u}_{d}(\xi),k=1,2.

Discrete analogue of the Laplacian is the following

(Δdud)(x~)=(Δ1(2)ud)(x~)+(Δ2(2)ud)(x~),(\Delta_{d}u_{d})(\tilde{x})=(\Delta_{1}^{(2)}u_{d})(\tilde{x})+(\Delta_{2}^{(2)}u_{d})(\tilde{x}),

so that its Fourier transform is

(Δdud)~(ξ)=h2((eihξ11)2+(eihξ21)2)u~d(ξ).\widetilde{(\Delta_{d}u_{d})}(\xi)=h^{-2}((e^{-ih\cdot\xi_{1}}-1)^{2}+(e^{-ih\cdot\xi_{2}}-1)^{2})\tilde{u}_{d}(\xi).

We use such discrete objects for constructing discrete Sobolev–Slobodetskii spaces to study wide class of discrete equations.

First, we introduce discrete analogue of the Schwartz space S(h2)S(h\mathbb{Z}^{2}) as a set of discrete functions with finite semi-norms

|ud|=supx~h2(1+|x~|)l|Δ(𝐤)ud(x~)||u_{d}|=\sup\limits_{\tilde{x}\in h\mathbb{Z}^{2}}(1+|\tilde{x}|)^{l}|\Delta^{({\bf k})}u_{d}(\tilde{x})|

for arbitrary l,𝐤=(k1,k2),kr,r=1,2l\in\mathbb{N},{\bf k}=(k_{1},k_{2}),k_{r}\in\mathbb{N},r=1,2,

Δ(𝐤)ud(x~)=Δ1k1Δ2k2ud(x~).\Delta^{({\bf k})}u_{d}(\tilde{x})=\Delta^{k_{1}}_{1}\Delta^{k_{2}}_{2}u_{d}(\tilde{x}).

Definition 1. A discrete distribution is called a linear continuous functional defined on the space S(h2)S(h\mathbb{Z}^{2}).

A set of such distributions will be denoted by S(h2)S^{\prime}(h\mathbb{Z}^{2}), and a value of the discrete distribution fdf_{d} on the test discrete function udS(h2)u_{d}\in S(h\mathbb{Z}^{2}) will be denoted by (fd,ud)(f_{d},u_{d}).

One can introduce a concept of a support for a discrete distribution. Namely, a support of the discrete function udS(h2)u_{d}\in S(h\mathbb{Z}^{2}) is a subset of the set h2h\mathbb{Z}^{2} such that ud(x~)0u_{d}(\tilde{x})\neq 0 for all points x~\tilde{x} from this subset. For an arbitrary set M2M\subset\mathbb{R}^{2} we denote Md=Mh2M_{d}=M\cap h\mathbb{Z}^{2}, and then one says that fd=0f_{d}=0 in the discrete domain MdM_{d} if (fd,ud)=0,udS(Md),(f_{d},u_{d})=0,\forall u_{d}\in S(M_{d}), where S(Md)S(h2)S(M_{d})\subset S(h\mathbb{Z}^{2}) consists of discrete functions with supports in MdM_{d}. If M~d\widetilde{M}_{d} is a union of such MdM_{d} where fd=0f_{d}=0 then support of the discrete distribution fdf_{d} is the set h2M~dh\mathbb{Z}^{2}\setminus\widetilde{M}_{d}.

Similarly [26] we can define standard operations in the space S(h2)S^{\prime}(h\mathbb{Z}^{2}), but differentiation will be changed by divided difference of first order. These operations are described in [21] in details, a convergence is meant as a weak convergence in the space S(h2)S^{\prime}(h\mathbb{Z}^{2}).

Example 1. If the function fd(x~)f_{d}(\tilde{x}) is locally summable then it generates the discrete distribution

(1) (fd,ud)=x~h2fd(x~)ud(x~)h2,udS(h2).(f_{d},u_{d})=\sum\limits_{\tilde{x}\in h\mathbb{Z}^{2}}f_{d}(\tilde{x})u_{d}(\tilde{x})h^{2},\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \forall u_{d}\in S(h\mathbb{Z}^{2}).

But there are different possibilities, for example, analogue of the Dirac mass-function

(δd,ud)=ud(0),(\delta_{d},u_{d})=u_{d}(0),

which can not be represented by the formula (1).

Let ζ2=h2((eihξ11)2+(eihξ21)2)\zeta^{2}=h^{-2}((e^{-ih\cdot\xi_{1}}-1)^{2}+(e^{-ih\cdot\xi_{2}}-1)^{2}). We introduce the following definition.

Definition 2. The space Hs(h2)H^{s}(h\mathbb{Z}^{2}) consists of discrete distributions and it is a closure of the space S(h2)S(h\mathbb{Z}^{2}) with respect to the norm

(2) uds=(𝕋2(1+|ζ2|)s|u~d(ξ)|2𝑑ξ)1/2.||u_{d}||_{s}=\left(\int\limits_{\hbar\mathbb{T}^{2}}(1+|\zeta^{2}|)^{s}|\tilde{u}_{d}(\xi)|^{2}d\xi\right)^{1/2}.

Let us remind that a lot of properties of such discrete spaces were studied in [7], Varying the parameter hh in (2) we obtain different norms which are equivalent to the L2L_{2}-norm. But constants in this equivalence depend on hh. In our constructions all constants do not depend on hh.

Definition 3. T͡he space Hs(Kd)H^{s}(K_{d}) consists of discrete distributions from Hs(h2)H^{s}(h\mathbb{Z}^{2}) such that their supports belong to the set Kd¯\overline{K_{d}}. A norm in the space Hs(Kd)H^{s}(K_{d}) is induced by the norm of the space Hs(h2)H^{s}(h\mathbb{Z}^{2}). The space H0s(Kd)H^{s}_{0}(K_{d}) consists of discrete distributions fdS(h2)f_{d}\in S^{\prime}(h\mathbb{R}^{2}) with supports inside of KdK_{d}, and these discrete distributions must admit a continuation into the space Hs(h2)H^{s}(h\mathbb{Z}^{2}). A norm in the space H0s(Kd)H^{s}_{0}(K_{d}) is given by the formula

fds+=inffds,||f_{d}||^{+}_{s}=\inf||\ell f_{d}||_{s},

where infimum is taken for all continuations \ell.

The Fourier image of the space Hs(Kd)H^{s}(K_{d}) will be denoted by H~s(Kd)\widetilde{H}^{s}(K_{d}).

2.2. Symbols, operators and projectors

Let A~d(ξ)\widetilde{A}_{d}(\xi) be a measurable periodic function in 2\mathbb{R}^{2} with basic square of periods 𝕋2\hbar\mathbb{T}^{2}. Such functions we call symbols.

Definition 4. A digital pseudo-differential operator AdA_{d} with the symbol Ad(ξ)A_{d}(\xi) in the discrete quadrant KdK_{d} is called an operator of the following type

(3) (Adud)(x~)=y~h2h2𝕋2A~d(ξ)ei(x~y~)ξu~d(ξ)𝑑ξ,x~Kd,(A_{d}u_{d})(\tilde{x})=\sum\limits_{\tilde{y}\in h\mathbb{Z}^{2}}h^{2}\int\limits_{\hbar\mathbb{T}^{2}}\widetilde{A}_{d}(\xi)e^{i(\tilde{x}-\tilde{y})\cdot\xi}\tilde{u}_{d}(\xi)d\xi,\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{x}\in K_{d},

We say that the operator AdA_{d} is elliptic one if

essinfξ𝕋2|Ad(ξ)|>0.ess\leavevmode\nobreak\ \inf_{\xi\in\hbar\mathbb{T}^{2}}|A_{d}(\xi)|>0.

A more general digital pseudo-differential operator with the symbol A~d(x~,ξ)\widetilde{A}_{d}(\tilde{x},\xi) depending on a spatial variable x~\tilde{x}

(Adud)(x~)=y~h2h2𝕋2Ad(x~,ξ)ei(x~y~)ξu~d(ξ)𝑑ξ,x~Kd,(A_{d}u_{d})(\tilde{x})=\sum\limits_{\tilde{y}\in h\mathbb{Z}^{2}}h^{2}\int\limits_{\hbar\mathbb{T}^{2}}A_{d}(\tilde{x},\xi)e^{i(\tilde{x}-\tilde{y})\cdot\xi}\tilde{u}_{d}(\xi)d\xi,\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{x}\in K_{d},

can be defined in the same way, but here we consider only operators of type (3).

We consider symbols satisfying the condition

(4) c1(1+|ζ2|)α/2|Ad(ξ)|c2(1+|ζ2|)α/2c_{1}(1+|\zeta^{2}|)^{\alpha/2}\leq|A_{d}(\xi)|\leq c_{2}(1+|\zeta^{2}|)^{\alpha/2}

with constants c1,c2c_{1},c_{2} non-depending on hh. The number α\alpha\in\mathbb{R} is called an order of digital pseudo-differential operator AdA_{d}.

The following simple result can be proved easily.

Lemma 1. A digital pseudo-differential operator AdA_{d} with the symbol A~d(ξ)\widetilde{A}_{d}(\xi) is a linear bounded operator Hs(h2)Hsα(h2)H^{s}(h\mathbb{Z}^{2})\to H^{s-\alpha}(h\mathbb{Z}^{2}) with a norm non-depending on hh.

We study a solvability of the discrete equation

(5) (Adud)(x~)=vd(x~),x~Kd,(A_{d}u_{d})(\tilde{x})=v_{d}(\tilde{x}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{x}\in K_{d},

in the space Hs(Kd)H^{s}(K_{d}) assuming that vdH0sα(Kd)v_{d}\in H^{s-\alpha}_{0}(K_{d}).

We will use certain special domain in two-dimensional complex space 2.\mathbb{C}^{2}. A domain of the type 𝒯h(K)=𝕋2+iK{\mathcal{T}}_{h}(K)=\hbar\mathbb{T}^{2}+iK is called a tube domain over the quadrant KK, and we will consider analytical functions f(x+iτ)f(x+i\tau) in the domain 𝒯h(K)=𝕋2+iK{\mathcal{T}}_{h}(K)=\hbar\mathbb{T}^{2}+iK.

Let us introduce the periodic Bochner kernel similar [26]

Bh(z)=x~Kdeix~(ξ+iτ)h2,ξ𝕋2,τK,B_{h}(z)=\sum\limits_{\tilde{x}\in K_{d}}e^{i\tilde{x}\cdot(\xi+i\tau)}h^{2},\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \xi\in\hbar\mathbb{T}^{2},\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tau\in K,

and corresponding integral operator

(Bhu~d)(ξ)=limτ0,τK14π2𝕋2Bh(ξ+iτη)u~d(η)𝑑η.(B_{h}\tilde{u}_{d})(\xi)=\lim\limits_{\tau\to 0,\tau\in K}\frac{1}{4\pi^{2}}\int\limits_{\hbar\mathbb{T}^{2}}B_{h}(\xi+i\tau-\eta)\tilde{u}_{d}(\eta)d\eta.

Lemma 2. For the quadrant KK the operator BhB_{h} has the following form

(Bhu~d)(ξ)=h28π2𝕋2u~d(η)𝑑η+limτ0+ih8π2𝕋2coth(ξ1η1+iτ1)2u~d(η)𝑑η+(B_{h}\tilde{u}_{d})(\xi)=\frac{h^{2}}{8\pi^{2}}\int\limits_{\mathbb{T}^{2}}\tilde{u}_{d}(\eta)d\eta+\lim\limits_{\tau\to 0+}\frac{ih}{8\pi^{2}}\int\limits_{\mathbb{T}^{2}}\cot\frac{h(\xi_{1}-\eta_{1}+i\tau_{1})}{2}\tilde{u}_{d}(\eta)d\eta+
+limτ0+ih8π2𝕋2coth(ξ2η2+iτ2)2u~d(η)𝑑η+\lim\limits_{\tau\to 0+}\frac{ih}{8\pi^{2}}\int\limits_{\mathbb{T}^{2}}\cot\frac{h(\xi_{2}-\eta_{2}+i\tau_{2})}{2}\tilde{u}_{d}(\eta)d\eta-
limτ0+h28π2𝕋2coth(ξ1η1+iτ1)2coth(ξ2η2+iτ2)2u~d(η)𝑑η,-\lim\limits_{\tau\to 0+}\frac{h^{2}}{8\pi^{2}}\int\limits_{\mathbb{T}^{2}}\cot\frac{h(\xi_{1}-\eta_{1}+i\tau_{1})}{2}\cot\frac{h(\xi_{2}-\eta_{2}+i\tau_{2})}{2}\tilde{u}_{d}(\eta)d\eta,

and BhB_{h} is a linear bounded operator Hs(𝕋2)Hs(𝕋2)H^{s}(\hbar\mathbb{T}^{2})\rightarrow H^{s}(\hbar\mathbb{T}^{2}) for |s|<1/2|s|<1/2. Moreover, the operator BhB_{h} is a projector H~s(h2)H~s(Kd)\widetilde{H}^{s}(h\mathbb{Z}^{2})\rightarrow\widetilde{H}^{s}(K_{d}).

Proof. Corresponding calculations for one-dimensional discrete cone were done in [20]. We use these evaluations adapting to our two-dimensional case. Since

x~kh+eix~kzkh=h2ih2cothzk2,zk=ξk+iτk,k=1,2.\sum\limits_{\tilde{x}_{k}\in h\mathbb{Z}_{+}}e^{-i\tilde{x}_{k}z_{k}}h=\frac{h}{2}-\frac{ih}{2}\cot\frac{hz_{k}}{2},\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ z_{k}=\xi_{k}+i\tau_{k},\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ k=1,2.

then multiplying two factors and applying the Fourier property on correspondence between a product and convolution we obtain the assertion.

Boundedness of the one-dimensional operator with the kernel hcothz2h\cot\frac{hz}{2}for |s|<1/2|s|<1/2 was proved in [7], Theorem 6; two-dimensional case can be considered by the same method. \blacksquare

Remark 1. The operator BhB_{h} is so called periodic bi-singular operator. Using classical results for Cauchy type integral [8, 10] one can evaluate the boundary value, but it is not important this time. Since these formulas are very huge we can do some simplifications without lost of generality. For example, we can consider the space S1(h2)S1(h2)S_{1}(h\mathbb{Z}^{2})\subset S_{1}(h\mathbb{Z}^{2}) with zeroes in coordinate axes and consider the space Hs(h2)H^{s}(h\mathbb{Z}^{2}) as closure of the set S1(h2)S_{1}(h\mathbb{Z}^{2}) assuming that all functions of discrete variable vanish on coordinate axes. For this case the first three summands in BhB_{h} will be zero.

Lemma 3. If |s|<1/2|s|<1/2 then the space H~s(h2)\widetilde{H}^{s}(h\mathbb{Z}^{2}) is uniquely represented as the direct sum

H~s(h2)=H~s(Kd)H~s(h2Kd)\widetilde{H}^{s}(h\mathbb{Z}^{2})=\widetilde{H}^{s}(K_{d})\oplus\widetilde{H}^{s}(h\mathbb{Z}^{2}\setminus K_{d})

Proof. It is simple consequence of Lemma 2. Indeed, the unique representation of the function f~H~(h2)\tilde{f}\in\widetilde{H}(h\mathbb{Z}^{2}) is the following

f~=Bhf~+(IBh)f~.\tilde{f}=B_{h}\tilde{f}+(I-B_{h})\tilde{f}.

A uniqueness of the such representation is possible only for |s|<1/2|s|<1/2. \blacksquare

To describe a solvability picture for the discrete equation (5) we need some additional elements of multidimensional complex analysis. We give it in the next section.

3. Periodic wave factorization

This concept is a periodic analogue of the wave factorization [18]. Some first preliminary considerations and results were described in [25, 24, 23, 22].

Definition 5. A periodic wave factorization for the elliptic symbol Ad(ξ)EαA_{d}(\xi)\in E_{\alpha} is called its representation in the form

Ad(ξ)=Ad,(ξ)Ad,=(ξ),A_{d}(\xi)=A_{d,\neq}(\xi)A_{d,=}(\xi),

where the factors Ad,(ξ),Ad,=(ξ)A_{d,\neq}(\xi),A_{d,=}(\xi) admit analytical continuation into tube domains 𝒯h(K),𝒯h(K){\mathcal{T}}_{h}(K),{\mathcal{T}}_{h}(-K) respectively with estimates

c1(1+|ζ^2|)æ2|Ad,(ξ+iτ)|c1(1+|ζ^2|)æ2,c_{1}(1+|\hat{\zeta}^{2}|)^{\frac{\ae}{2}}\leq|A_{d,\neq}(\xi+i\tau)|\leq c^{\prime}_{1}(1+|\hat{\zeta}^{2}|)^{\frac{\ae}{2}},
c2(1+|ζ^2|)αæ2|Ad,=(ξiτ)|c2(1+|ζ^2|)αæ2,c_{2}(1+|\hat{\zeta}^{2}|)^{\frac{{\alpha-\ae}}{2}}\leq|A_{d,=}(\xi-i\tau)|\leq c^{\prime}_{2}(1+|\hat{\zeta}^{2}|)^{\frac{{\alpha-\ae}}{2}},

and constants c1,c1,c2,c2c_{1},c^{\prime}_{1},c_{2},c^{\prime}_{2} non-depending on hh, where

ζ^22((eih(ξ1+iτ1)1)2+(eih(ξ2+iτ2)1)2),\hat{\zeta}^{2}\equiv\hbar^{2}\left((e^{-ih(\xi_{1}+i\tau_{1})}-1)^{2}+(e^{-ih(\xi_{2}+i\tau_{2})}-1)^{2}\right),
ξ=(ξ1,ξ2)𝕋2,τ(τ1,τ2)K.\xi=(\xi_{1},\xi_{2})\in\hbar\mathbb{T}^{2},\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tau-(\tau_{1},\tau_{2})\in K.

The number æ\ae\in{\mathbb{R}} is called an index of periodic wave factorization.

Unfortunately, we have no an algorithm to construct the factorization. But the are certain examples of periodic symbols which admit such factorization. We give one of them.

If ff is an arbitrary function of a discrete variable, fS(h2)f\in S(h\mathbb{Z}^{2}), suppfKd(Kd)supp\leavevmode\nobreak\ f\subset K_{d}\cup(-K_{d}) then we have

f=χ+f+χf,f=\chi_{+}f+\chi_{-}f,

where χ±\chi_{\pm} are indicators of ±Kd\pm K_{d}. Applying the discrete Fourier transform we obtain the representation f~=f~++f~\tilde{f}=\tilde{f}_{+}+\tilde{f}_{-}, and f~±\tilde{f}_{\pm} admit an analytical continuation into 𝒯h(±K){\mathcal{T}}_{h}(\pm K) according to Lemma 2. Thus, we can write expf~=expf~+expf~\exp\tilde{f}=\exp\tilde{f}_{+}\cdot\exp\tilde{f}_{-}, therefore we obtain periodic wave factorization with index zero for the function expf~\exp\tilde{f}.

Everywhere below we assume existence of such periodic wave factorization for the symbol Ad(ξ)A_{d}(\xi) with index æ\ae.

3.1. A unique solvability

This section is devoted to most simple case when a solution of the equation (5) exists and it is unique.

Theorem 1. Let |æs|<1/2|\ae-s|<1/2. Then the equation (5) has a unique solution for arbitrary right hand side vdH0sα(Kd)v_{d}\in H^{s-\alpha}_{0}(K_{d}), and it is given by the formula

u~d(ξ)Ad,1(ξ)Bh(Ad,=1(ξ)(vd)~(ξ)),\tilde{u}_{d}(\xi)-A^{-1}_{d,\neq}(\xi)B_{h}(A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi)),

where vd\ell v_{d} is an arbitrary continuation of vdv_{d} into Hsα(h2)H^{s-\alpha}(h\mathbb{Z}^{2}).

Proof Let vd\ell v_{d} be an arbitrary continuation of vdH0sα(Kd)v_{d}\in H^{s-\alpha}_{0}(K_{d}) into Hsα(h2)H^{s-\alpha}(h\mathbb{Z}^{2}). Let us introduce the function

wd(x~)=(vd)(x~)(Adud)(x~),w_{d}(\tilde{x})=(\ell v_{d})(\tilde{x})-(A_{d}u_{d})(\tilde{x}),

so that w(x~)=0w(\tilde{x})=0 for x~Kd\tilde{x}\notin K_{d}.

Now we write the equation (5) in the form

(Adud)(x~)+wd(x~)=(vd)(x~),x~h2,(A_{d}u_{d})(\tilde{x})+w_{d}(\tilde{x})=(\ell v_{d})(\tilde{x}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{x}\in h\mathbb{Z}^{2},

and after applying the discrete Fourier transform and periodic wave factorization we obtain

(6) Ad,(ξ)u~d(ξ)+Ad,=1(ξ)w~d(ξ)=Ad,=1(ξ)(vd)~(ξ),ξ𝕋2,A_{d,\neq}(\xi)\tilde{u}_{d}(\xi)+A^{-1}_{d,=}(\xi)\tilde{w}_{d}(\xi)=A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \xi\in\hbar\mathbb{T}^{2},

We have the following inclusions according to Lemma 1 and Lemma 2

Ad,(ξ)u~d(ξ)H~sæ(Kd),Ad,=1(ξ)w~d(ξ)H~sæ(h2Kd),A_{d,\neq}(\xi)\tilde{u}_{d}(\xi)\in\widetilde{H}^{s-\ae}(K_{d}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ A^{-1}_{d,=}(\xi)\tilde{w}_{d}(\xi)\in\widetilde{H}^{s-\ae}(h\mathbb{Z}^{2}\setminus K_{d}),
Ad,=1(ξ)(vd)~(ξ)H~sæ(h2),A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi)\in\widetilde{H}^{s-\ae}(h\mathbb{Z}^{2}),

and then according to Lemma 3 the right hand side of the equality (6) is uniquely represented by the sum

Ad,=1(ξ)(vd)~(ξ)=fd+(ξ)+fd(ξ),A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi)=f^{+}_{d}(\xi)+f^{-}_{d}(\xi),

where

fd+(ξ)=Bh(Ad,=1(ξ)(vd)~(ξ)),fd(ξ)=(IBh)(Ad,=1(ξ)(vd)~(ξ))f^{+}_{d}(\xi)=B_{h}(A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi)),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ f^{-}_{d}(\xi)=(I-B_{h})(A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi))

.

Further, we rewrite the equality (6)

Ad,(ξ)u~d(ξ)fd+(ξ)=fd(ξ)Ad,=1(ξ)w~d(ξ)A_{d,\neq}(\xi)\tilde{u}_{d}(\xi)-f^{+}_{d}(\xi)=f^{-}_{d}(\xi)-A^{-1}_{d,=}(\xi)\tilde{w}_{d}(\xi)

and using the uniqueness of the representation as the direct sum H~sæ(Kd)H~sæ(h2Kd)\widetilde{H}^{s-\ae}(K_{d})\oplus\widetilde{H}^{s-\ae}(h\mathbb{Z}^{2}\setminus K_{d}) we conclude that both left hand side and right hand side should be zero. Thus,

u~d(ξ)=Ad,1(ξ)Bh(Ad,=1(ξ)(vd)~(ξ)),\tilde{u}_{d}(\xi)=A^{-1}_{d,\neq}(\xi)B_{h}(A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi)),

and Theorem 1 is proved. \blacksquare

4. Discrete boundary value problem

In this section we consider more interesting case when the equation (5) has a lot of solutions.

4.1. Form of a discrete solution

This section uses some results from [21] concerning a form of a discrete distribution supported at the origin.

Theorem 2. Let æs=n+δ,n,|δ|<1/2\ae-s=n+\delta,n\in\mathbb{N},|\delta|<1/2. Then a general solution of the equation (5) has the following form

u~d(ξ)=Ad,1(ξ)Qn(ξ)Bh(Qn1(ξ)Ad,=1(ξ)(vd)~(ξ))+\tilde{u}_{d}(\xi)=A^{-1}_{d,\neq}(\xi)Q_{n}(\xi)B_{h}(Q_{n}^{-1}(\xi)A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi))+
+Ad,1(ξ)(k=0n1c~k(ξ1)ζ^2k+d~k(ξ2)ζ^1k),+A^{-1}_{d,\neq}(\xi)\left(\sum\limits_{k=0}^{n-1}\tilde{c}_{k}(\xi_{1})\hat{\zeta}_{2}^{k}+\tilde{d}_{k}(\xi_{2})\hat{\zeta}_{1}^{k}\right),

where Qn(ξ)Q_{n}(\xi) is an arbitrary polynomial of order nn of variables ζk=(eihξk1),k=1,2,\zeta_{k}=\hbar(e^{-ih\xi_{k}}-1),k=1,2, satisfying the condition (4) for α=n\alpha=n, c~k(ξ1),d~k(ξ2),k=0,1,,n1,\tilde{c}_{k}(\xi_{1}),\tilde{d}_{k}(\xi_{2}),k=0,1,\cdots,n-1, – are arbitrary functions from Hsk(h𝕋),sk=sæ+k1/2H^{s_{k}}(h\mathbb{T}),s_{k}=s-\ae+k-1/2.

The a priori estimate

udsconst(fsα++k=0n1([ck]sk+[dk]sk)),||u_{d}||_{s}\leq const\left(||f||^{+}_{s-\alpha}+\sum\limits_{k=0}^{n-1}([c_{k}]_{s_{k}}+[d_{k}]_{s_{k}})\right),

holds, where []sk[\cdot]_{s_{k}} denotes a norm in Hsk(h𝕋)H^{s_{k}}(h\mathbb{T}), and constconst does not depend on hh.

Proof. We start from the equality (6). Let Qn(ξ)Q_{n}(\xi) be an arbitrary polynomial of order nn of variables ζk=(eihξk1),k=1,2,\zeta_{k}=\hbar(e^{-ih\xi_{k}}-1),k=1,2, satisfying the condition (4) for α=n\alpha=n. We multiply the equality (6) by Qn1(ξ)Q^{-1}_{n}(\xi)

(7) Qn1(ξ)Ad,(ξ)u~d(ξ)+Qn1(ξ)Ad,=1(ξ)w~d(ξ)=Qn1(ξ)Ad,=1(ξ)(vd)~(ξ),ξ𝕋2,Q^{-1}_{n}(\xi)A_{d,\neq}(\xi)\tilde{u}_{d}(\xi)+Q^{-1}_{n}(\xi)A^{-1}_{d,=}(\xi)\tilde{w}_{d}(\xi)=Q^{-1}_{n}(\xi)A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \xi\in\hbar\mathbb{T}^{2},

We have in view of Lemma 1

Qn1(ξ)Ad,=1(ξ)(vd)~(ξ)H~sæ+n(h2),Q^{-1}_{n}(\xi)A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi)\in\widetilde{H}^{s-\ae+n}(h\mathbb{Z}^{2}),

and since sæ+n=δs-\ae+n=-\delta then according to Lemma 3 we write the unique decomposition

Qn1(ξ)Ad,=1(ξ)(vd)~(ξ)=Fd+(ξ)+Fd(ξ),Q^{-1}_{n}(\xi)A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi)=F^{+}_{d}(\xi)+F^{-}_{d}(\xi),

where

Fd+(ξ)=Bh(Qn1(ξ)Ad,=1(ξ)(vd)~(ξ)),Fd(ξ)=(IBh)(Qn1(ξ)Ad,=1(ξ)(vd)~(ξ)).F^{+}_{d}(\xi)=B_{h}(Q^{-1}_{n}(\xi)A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi)),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ F^{-}_{d}(\xi)=(I-B_{h})(Q^{-1}_{n}(\xi)A^{-1}_{d,=}(\xi)\widetilde{(\ell v_{d})}(\xi)).

Taking into account this fact we rewrite the equality (7) in the form

Ad,(ξ)u~d(ξ)+Ad,=1(ξ)w~d(ξ)=Qn(ξ)Fd+(ξ)+Qn(ξ)Fd(ξ),A_{d,\neq}(\xi)\tilde{u}_{d}(\xi)+A^{-1}_{d,=}(\xi)\tilde{w}_{d}(\xi)=Q_{n}(\xi)F^{+}_{d}(\xi)+Q_{n}(\xi)F^{-}_{d}(\xi),

and further,

Ad,(ξ)u~d(ξ)Qn(ξ)Fd+(ξ)=Qn(ξ)Fd(ξ)Ad,=1(ξ)w~d(ξ),A_{d,\neq}(\xi)\tilde{u}_{d}(\xi)-Q_{n}(\xi)F^{+}_{d}(\xi)=Q_{n}(\xi)F^{-}_{d}(\xi)-A^{-1}_{d,=}(\xi)\tilde{w}_{d}(\xi),

Since Fd+(ξ)H~sæ+n(Kd),Fd(ξ)H~sæ+n(h2Kd)F^{+}_{d}(\xi)\in\widetilde{H}^{s-\ae+n}(K_{d}),F^{-}_{d}(\xi)\in\widetilde{H}^{s-\ae+n}(h\mathbb{Z}^{2}\setminus K_{d}) then according to Lemma 1 we conclude Qn(ξ)Fd+(ξ)H~sæ(Kd),Qn(ξ)Fd(ξ)H~sæ(h2Kd)Q_{n}(\xi)F^{+}_{d}(\xi)\in\widetilde{H}^{s-\ae}(K_{d}),Q_{n}(\xi)F^{-}_{d}(\xi)\in\widetilde{H}^{s-\ae}(h\mathbb{Z}^{2}\setminus K_{d}). Applying the inverse discrete Fourier transform we obtain an equality for two discrete distributions. The left hand side vanishes at least under one condition x~1<0\tilde{x}_{1}<0 or x~2<0\tilde{x}_{2}<0, and the right hand side vanishes under the condition .x~1>0,x~2>0\tilde{x}_{1}>0,\tilde{x}_{2}>0, Thus, it should be a discrete distribution supported on sides of the discrete quadrant {(x~1,x~2)h2:{x~1>0,x~2=0}{x~1=0,x~2>0}}\{(\tilde{x}_{1},\tilde{x}_{2})\in h\mathbb{Z}^{2}:\{\tilde{x}_{1}>0,\tilde{x}_{2}=0\}\cup\{\tilde{x}_{1}=0,\tilde{x}_{2}>0\}\}. Using corresponding result from [21] we obtain the following form for this discrete distribution

k=0n1(ck(x~1)(Δ2(k)δd)(x~2)+dk(x~2)(Δ1(k)δd)(x~1)),\sum\limits_{k=0}^{n-1}\left(c_{k}(\tilde{x}_{1})(\Delta_{2}^{(k)}\delta_{d})(\tilde{x}_{2})+d_{k}(\tilde{x}_{2})(\Delta_{1}^{(k)}\delta_{d})(\tilde{x}_{1})\right),

where all summands should be elements of the space Hsæ(h2)H^{s-\ae}(h\mathbb{Z}^{2}).

The left question is how much summands we need in the right-hand side. Counting principle is a very simple because every summand should belong to the space H~s(𝕋2)\widetilde{H}^{s}(\hbar\mathbb{T}^{2}).

Let us consider the summand ck(ξ1)ζ2kc_{k}(\xi_{1})\zeta_{2}^{k}. Taking into account that order of Ad,+1(ξ)A^{-1}_{d,+}(\xi) is æ-\ae we need to verify the finiteness of the HsæH^{s-\ae}-norm for ck(ξ1)ζ2kc_{k}(\xi_{1})\zeta_{2}^{k}. We have

||ck(Δ2(k)δd)||sæ2=𝕋2(1+|ζ2|)sæ||ck(ξ1)ζ2k|2dξ||c_{k}(\Delta^{(k)}_{2}\delta_{d})||^{2}_{s-\ae}=\int\limits_{\hbar\mathbb{T}^{2}}(1+|\zeta^{2}|)^{s-\ae}||c_{k}(\xi_{1})\zeta_{2}^{k}|^{2}d\xi
=𝕋2(1+|ζ2|)sæ||ck(ξ1)|2|ζ2k|2dξa12(sæ+k+1/2)𝕋|ck(ξ1)|2dξ1=\int\limits_{\hbar\mathbb{T}^{2}}(1+|\zeta^{2}|)^{s-\ae}||c_{k}(\xi_{1})|^{2}|\zeta_{2}^{k}|^{2}d\xi\leq a_{1}\hbar^{2(s-\ae+k+1/2)}\int\limits_{\hbar\mathbb{T}}|c_{k}(\xi_{1})|^{2}d\xi_{1}
a2𝕋(1+|ζ12|)sæ+k+1/2|ck(ξ1)|2𝑑ξ1,\leq a_{2}\int\limits_{\hbar\mathbb{T}}(1+|{\zeta_{1}}^{2}|)^{s-\ae+k+1/2}|c_{k}(\xi_{1})|^{2}d\xi_{1},

and the constants a1,a2a_{1},a_{2} do not depend on hh.

The last summand should be (n1)(n-1)th because for nnth summand we obtain a positive growth: for k=nk=n we have sn=sæn+1/2=nδ+n+1/2=δ+1/2>0s_{n}=s-\ae-n+1/2=-n-\delta+n+1/2=-\delta+1/2>0.

A priori estimates can be obtained in the same way described in [21]. \blacksquare

4.2. The Dirichlet discrete boundary condition

We consider here first simple case with discrete Dirichlet boundary conditions. We suppose in this section that æs=1+δ,|δ|<1/2,vd0.\ae-s=1+\delta,|\delta|<1/2,v_{d}\equiv 0. It follows from Theorem 2 that we have the following general solution of the equation (5)

(8) u~d(ξ)=Ad,1(ξ)(c~0(ξ1)+d~0(ξ2)),\tilde{u}_{d}(\xi)=A^{-1}_{d,\neq}(\xi)(\tilde{c}_{0}(\xi_{1})+\tilde{d}_{0}(\xi_{2})),

where c0,d0Hsæ1/2()c_{0},d_{0}\in H^{s-\ae-1/2}(\hbar\mathbb{Z}) are arbitrary functions. To determine uniquely these functions we add the discrete Dirichlet conditions on angle sides

(9) ud|x~1=0=fd(x~2),ud|x~2=0=gd(x~1).{u_{d}}_{|_{\tilde{x}_{1}=0}}=f_{d}(\tilde{x}_{2}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ {u_{d}}_{|_{\tilde{x}_{2}=0}}=g_{d}(\tilde{x}_{1}).

Thus, we have the discrete Dirichlet problem (5),(9).

First, we apply the discrete Fourier transform to discrete conditions (9) and obtain the following form

(10) ππu~d(ξ1,ξ2)𝑑ξ1=f~d(ξ2),ππu~d(ξ1,ξ2)𝑑ξ2=g~d(ξ1).\int\limits_{-\hbar\pi}^{\hbar\pi}\tilde{u}_{d}(\xi_{1},\xi_{2})d\xi_{1}=\tilde{f}_{d}(\xi_{2}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \int\limits_{-\hbar\pi}^{\hbar\pi}\tilde{u}_{d}(\xi_{1},\xi_{2})d\xi_{2}=\tilde{g}_{d}(\xi_{1}).

Substituting (10) into (8) we obtain the following relations

ππu~d(ξ)𝑑ξ1=ππAd,1(ξ)c~0(ξ1)𝑑ξ1+d~0(ξ2)ππAd,1(ξ)𝑑ξ1\int\limits_{-\hbar\pi}^{\hbar\pi}\tilde{u}_{d}(\xi)d\xi_{1}=\int\limits_{-\hbar\pi}^{\hbar\pi}A^{-1}_{d,\neq}(\xi)\tilde{c}_{0}(\xi_{1})d\xi_{1}+\tilde{d}_{0}(\xi_{2})\int\limits_{-\hbar\pi}^{\hbar\pi}A^{-1}_{d,\neq}(\xi)d\xi_{1}
ππu~d(ξ)𝑑ξ2=c~0(ξ1)ππAd,1(ξ)𝑑ξ2+ππAd,1(ξ)d~0(ξ2)𝑑ξ2\int\limits_{-\hbar\pi}^{\hbar\pi}\tilde{u}_{d}(\xi)d\xi_{2}=\tilde{c}_{0}(\xi_{1})\int\limits_{-\hbar\pi}^{\hbar\pi}A^{-1}_{d,\neq}(\xi)d\xi_{2}+\int\limits_{-\hbar\pi}^{\hbar\pi}A^{-1}_{d,\neq}(\xi)\tilde{d}_{0}(\xi_{2})d\xi_{2}

Let us denote

ππAd,1(ξ)𝑑ξ1a~0(ξ2),ππAd,1(ξ)𝑑ξ2b~0(ξ1)\int\limits_{-\hbar\pi}^{\hbar\pi}A^{-1}_{d,\neq}(\xi)d\xi_{1}\equiv\tilde{a}_{0}(\xi_{2}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \int\limits_{-\hbar\pi}^{\hbar\pi}A^{-1}_{d,\neq}(\xi)d\xi_{2}\equiv\tilde{b}_{0}(\xi_{1})

and suppose that a~0(ξ2),b~0(ξ1)0,ξ10,ξ20\tilde{a}_{0}(\xi_{2}),\tilde{b}_{0}(\xi_{1})\neq 0,\forall\xi_{1}\neq 0,\xi_{2}\neq 0.

Therefore, we have the following system of two linear integral equations with respect to two unknown functions c~0(ξ1),d~0(ξ2)\tilde{c}_{0}(\xi_{1}),\tilde{d}_{0}(\xi_{2})

(11) {ππM1(ξ)c~0(ξ1)𝑑ξ1+d~0(ξ2)=F~d(ξ2)c~0(ξ1)+ππM2(ξ)d~0(ξ2)𝑑ξ2=G~d(ξ1),\left\{\begin{array}[]{rcl}\int\limits_{-\hbar\pi}^{\hbar\pi}M_{1}(\xi)\tilde{c}_{0}(\xi_{1})d\xi_{1}+\tilde{d}_{0}(\xi_{2})=\tilde{F}_{d}(\xi_{2})\\ \tilde{c}_{0}(\xi_{1})+\int\limits_{-\hbar\pi}^{\hbar\pi}M_{2}(\xi)\tilde{d}_{0}(\xi_{2})d\xi_{2}=\tilde{G}_{d}(\xi_{1}),\end{array}\right.

where we have used the following notations

F~d(ξ2)=f~d(ξ2)a~01(ξ2),G~d(ξ1)=g~d(ξ1)b~01(ξ1),\tilde{F}_{d}(\xi_{2})=\tilde{f}_{d}(\xi_{2})\tilde{a}_{0}^{-1}(\xi_{2}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{G}_{d}(\xi_{1})=\tilde{g}_{d}(\xi_{1})\tilde{b}_{0}^{-1}(\xi_{1}),
M1(ξ)=Ad,1(ξ)a~01(ξ2),M2(ξ)=Ad,1(ξ)b~01(ξ1).M_{1}(\xi)=A^{-1}_{d,\neq}(\xi)\tilde{a}_{0}^{-1}(\xi_{2}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ M_{2}(\xi)=A^{-1}_{d,\neq}(\xi)\tilde{b}_{0}^{-1}(\xi_{1}).

Unique solvability conditions for the system (11) will be equivalent to unique solvability for the discrete Dirichlet problem (5),(9).

Thus, we obtain the following result.

Proposition 1. Let fd,gdHs1/2(+),s>1/2,vd0f_{d},g_{d}\in H^{s-1/2}(\mathbb{R}_{+}),s>1/2,v_{d}\equiv 0. Then the discrete Dirichlet problem (5),(9) is reduced to the equivalent system of linear integral equations (11).

4.3. Non-local discrete boundary condition

We consider here the æs=1+δ,|δ|<1/2\ae-s=1+\delta,|\delta|<1/2 for the equation (5) with different boundary conditions, namely

(12) x~1h+ud(x~1,x~2)h=fd(x~2),x~2h+ud(x~1,x~2)h=gd(x~1),x~h++ud(x~1,x~2)h2=0.\begin{array}[]{rcl}\sum\limits_{\tilde{x}_{1}\in h\mathbb{Z}_{+}}u_{d}(\tilde{x}_{1},\tilde{x}_{2})h=f_{d}(\tilde{x}_{2}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \sum\limits_{\tilde{x}_{2}\in h\mathbb{Z}_{+}}u_{d}(\tilde{x}_{1},\tilde{x}_{2})h=g_{d}(\tilde{x}_{1}),\\ \sum\limits_{\tilde{x}\in h\mathbb{Z}_{++}}u_{d}(\tilde{x}_{1},\tilde{x}_{2})h^{2}=0.\end{array}

These additional conditions will help us to determine uniquely the unknown functions c0,d0c_{0},d_{0} in the solution (8).

Indeed, using the discrete Fourier transform we rewrite the conditions (12) as follows

(13) u~d(0,ξ2)=f~d(ξ2),u~d(ξ1,0)=g~d(ξ1),u~d(0,0)=0.\tilde{u}_{d}(0,\xi_{2})=\tilde{f}_{d}(\xi_{2}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{u}_{d}(\xi_{1},0)=\tilde{g}_{d}(\xi_{1}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{u}_{d}(0,0)=0.

Now we substitute the formulas (13) into (8). The first two equality are

u~d(0,ξ2)=Ad,1(0,ξ2)(c~0(0)+d~0(ξ2))=f~d(ξ2),\tilde{u}_{d}(0,\xi_{2})=A^{-1}_{d,\neq}(0,\xi_{2})(\tilde{c}_{0}(0)+\tilde{d}_{0}(\xi_{2}))=\tilde{f}_{d}(\xi_{2}),
u~d(ξ1,0)=Ad,1(ξ1,0)(c~0(ξ1)+d~0(0))=g~d(ξ1).\tilde{u}_{d}(\xi_{1},0)=A^{-1}_{d,\neq}(\xi_{1},0)(\tilde{c}_{0}(\xi_{1})+\tilde{d}_{0}(0))=\tilde{g}_{d}(\xi_{1}).

It implies the following relations according to the third condition

f~d(0)=g~d(0),and from whichc~0(0)+d~0(0)=0,andc~0(0)=d~0(0)=0.\tilde{f}_{d}(0)=\tilde{g}_{d}(0),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \text{and from which}\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{c}_{0}(0)+\tilde{d}_{0}(0)=0,\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \text{and}\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \tilde{c}_{0}(0)=\tilde{d}_{0}(0)=0.

Then we have at least formally

(14) u~d(ξ)=Ad,1(ξ)(Ad,(ξ1,0)g~d(ξ1)+Ad,(0,ξ2)f~d(ξ2))\tilde{u}_{d}(\xi)=A^{-1}_{d,\neq}(\xi)\left(A_{d,\neq}(\xi_{1},0)\tilde{g}_{d}(\xi_{1})+A_{d,\neq}(0,\xi_{2})\tilde{f}_{d}(\xi_{2})\right)

It is left to formulate and to prove exactly the obtained result.

Theorem 3. Let fd,gdHs+1/2(h),vd0f_{d},g_{d}\in H^{s+1/2}(h\mathbb{Z}),v_{d}\equiv 0. Then the discrete problem (5),(12) has unique solution which is given by the formula (14).

The a priori estimate

udsconst(fds+1/2+gds+1/2)||u_{d}||_{s}\leq const(||f_{d}||_{s+1/2}+||g_{d}||_{s+1/2})

holds with a const non-depending on hh

Proof. We need to prove the a priori estimate only. Let us consider the first summand

Ad,1(ξ)Ad,(ξ1,0)g~d(ξ1)s2=||A^{-1}_{d,\neq}(\xi)A_{d,\neq}(\xi_{1},0)\tilde{g}_{d}(\xi_{1})||^{2}_{s}=
=𝕋2|Ad,1(ξ1,ξ2)Ad,(ξ1,0)g~d(ξ1)|2(1+|ζ2|)s𝑑ξ1𝑑ξ2=\int\limits_{\hbar\mathbb{T}^{2}}|A^{-1}_{d,\neq}(\xi_{1},\xi_{2})A_{d,\neq}(\xi_{1},0)\tilde{g}_{d}(\xi_{1})|^{2}(1+|\zeta^{2}|)^{s}d\xi_{1}d\xi_{2}\leq
C2s𝕋2|gd(ξ1)|2𝑑ξC12s+1ππ|gd(ξ1)|2𝑑ξ1\leq C\hbar^{2s}\int\limits_{\hbar\mathbb{T}^{2}}|g_{d}(\xi_{1})|^{2}d\xi\leq C_{1}\hbar^{2s+1}\int\limits_{-\hbar\pi}^{\hbar\pi}|g_{d}(\xi_{1})|^{2}d\xi_{1}\leq
C2ππ|gd(ξ1)|2(1+|ζ12|)s+1/2dξ1=||gd|s+1/22.\leq C_{2}\int\limits_{-\hbar\pi}^{\hbar\pi}|g_{d}(\xi_{1})|^{2}(1+|\zeta_{1}^{2}|)^{s+1/2}d\xi_{1}=||g_{d}|^{2}_{s+1/2}.

The second summand has the same estimate. \blacksquare

5. A comparison between discrete and continuous solutions

The continuous analogue of the discrete boundary value problem is the following [19].

Let AA be a pseudo-differential operator with the symbol A(ξ),ξ=(ξ1,ξ2)A(\xi),\xi=(\xi_{1},\xi_{2}) satisfying the condition

c1(1+|ξ)α|A(ξ)|c2(1+|ξ)α.c_{1}(1+|\xi)^{\alpha}\leq|A(\xi)|\leq c_{2}(1+|\xi)^{\alpha}.

and admitting the wave factorization with respect to the quadrant KK with index æ\ae.

We consider the equation

(15) (Au)(x)=0,xK,(Au)(x)=0,\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ x\in K,

with the following additional conditions

(16) 0+u(x1,x2)𝑑x1=f(x2),0+u(x1,x2)𝑑x2=g(x1),Ku(x)𝑑x=0.\int\limits_{0}^{+\infty}u(x_{1},x_{2})dx_{1}=f(x_{2}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \int\limits_{0}^{+\infty}u(x_{1},x_{2})dx_{2}=g(x_{1}),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \int\limits_{-K}u(x)dx=0.

A solution of the problem (15),(16) is sought in the space Hs(K)H^{s}(K) [18] and boundary functions are taken from the space Hs+1/2(+)H^{s+1/2}(\mathbb{R}_{+}). Such problem was considered in [19] and it has the solution

(17) u~(ξ)=A1(ξ)(A(ξ1,0)g~(ξ1)+A(0,ξ2)f~(ξ2))\tilde{u}(\xi)=A^{-1}_{\neq}(\xi)\left(A_{\neq}(\xi_{1},0)\tilde{g}(\xi_{1})+A_{\neq}(0,\xi_{2})\tilde{f}(\xi_{2})\right)

under condition that the symbol A(ξ)A(\xi) admits the wave factorization with respect to the quadrant KK

A(ξ)=A(ξ)A=(ξ)A(\xi)=A_{\neq}(\xi)A_{=}(\xi)

with index æ\ae such that æs=1+δ,|δ|<1/2\ae-s=1+\delta,|\delta|<1/2.

To construct a discrete boundary value problem which is good approximation to (15),(16) we need to choose Ad(ξ)A_{d}(\xi) and fd,gdf_{d},g_{d} in a special way. First, we introduce the operator lhl_{h} which acts as follows. For a function uu defined in \mathbb{R} we take its Fourier transform f~\tilde{f} then we take its restriction on T\hbar T and periodically extend it to \mathbb{R}. Finally, we take its inverse discrete Fourier transform and obtain the function of discrete variable (lhu)(x~),x~h(l_{h}u)(\tilde{x}),\tilde{x}\in h\mathbb{R}. Thus, we put

fd=lhf,gd=lhg.f_{d}=l_{h}f,\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ g_{d}=l_{h}g.

Second, the symbol of digital operator AdA_{d} we construct in the same way. If we have the wave factorization for the symbol A(ξ)A(\xi) then we take restrictions of factors on 𝕋2\hbar\mathbb{T}^{2} and the periodic symbol Ad(ξ)A_{d}(\xi) is a product of these restrictions. For such fd,gdf_{d},g_{d} and the symbol Ad(ξ)A_{d}(\xi) we obtain the following result.

Theorem 4. Let f,gS(),æ>1f,g\in S(\mathbb{R}),\ae>1. Then we have the following estimate for solutions uu andudu_{d} of the continuous problem (15),(16) and the discrete one (5),(12)

|u(x~)ud(x~)|C(f,g)hβ,|u(\tilde{x})-u_{d}(\tilde{x})|\leq C(f,g)h^{\beta},

where the const C(f,g)C(f,g) depends on functions f,gf,g, β>0\beta>0 can be an arbitrary number.

Proof. We need to compare two functions (14) and (17), more exactly their inverse discrete Fourier transform and inverse Fourier transform at points x~Kd\tilde{x}\in K_{d}. We have

ud(x~)u(x~)=14π2(𝕋2eix~ξu~d(ξ)𝑑ξ2eix~ξu~(ξ)𝑑ξ)=u_{d}(\tilde{x})-u(\tilde{x})=\frac{1}{4\pi^{2}}\left(\int\limits_{\hbar\mathbb{T}^{2}}e^{i\tilde{x}\cdot\xi}\tilde{u}_{d}(\xi)d\xi-\int\limits_{\mathbb{R}^{2}}e^{i\tilde{x}\cdot\xi}\tilde{u}(\xi)d\xi\right)=
=14π22𝕋2eix~ξA1(ξ)(A(ξ1,0)g~(ξ1)+A(0,ξ2)f~(ξ2))𝑑ξ,=\frac{1}{4\pi^{2}}\int\limits_{\mathbb{R}^{2}\setminus\hbar\mathbb{T}^{2}}e^{i\tilde{x}\cdot\xi}A^{-1}_{\neq}(\xi)\left(A_{\neq}(\xi_{1},0)\tilde{g}(\xi_{1})+A_{\neq}(0,\xi_{2})\tilde{f}(\xi_{2})\right)d\xi,

since according to our choice for Ad,fd,gdA_{d},f_{d},g_{d} the functions u~d\tilde{u}_{d} and u~\tilde{u} coincide in points ξ𝕋2\xi\in\hbar\mathbb{T}^{2}.

We will estimate one summand.

|14π22𝕋2eix~ξA1(ξ)A(ξ1,0)g~(ξ1)𝑑ξ|\left|\frac{1}{4\pi^{2}}\int\limits_{\mathbb{R}^{2}\setminus\hbar\mathbb{T}^{2}}e^{i\tilde{x}\cdot\xi}A^{-1}_{\neq}(\xi)A_{\neq}(\xi_{1},0)\tilde{g}(\xi_{1})d\xi\right|\leq
Cπ+dξ2(1+|ξ1|+|ξ2|)æπ+|ξ1|γ𝑑ξ1,\leq\leavevmode\nobreak\ C\int\limits_{\hbar\pi}^{+\infty}\frac{d\xi_{2}}{(1+|\xi_{1}|+|\xi_{2}|)^{\ae}}\int\limits_{\hbar\pi}^{+\infty}|\xi_{1}|^{-\gamma}d\xi_{1},

since g~S()\tilde{g}\in S(\mathbb{R}). It implies the required estimate. \blacksquare

Conclusion

In this paper we have considered two-dimensional cone only, but the authors continue to work in multidimensional situations and we hope to obtain results similar to a discrete half-space [21, 15].

As first practical applications the authors plan to study discrete variant of a quarter-plane problem in diffraction theory [14] and elasticity theory [18]. We hope it will useful application of the developed theory.

References

  • [1] B Alexeev, Nonlocal astrophysics, Elsevier, Amsterdam, 2017.
  • [2] L. Botchway, R. Kibiti, and M Ruzhansky, Difference equations and pseudo-differential equations on 𝐳n{\bf z}^{n}, arXiv:1705.07564v1 [math.FA] (2017), 1–29.
  • [3] F Camilli, R.D. Maio, and A Tosin, Measure-valued solutions to nonlocal transport equations on networks, J. Differ. Equations 264 (2018), no. 12, 7213–7241.
  • [4] W. Chen and S Kim, Bubble-tower solutions to asymptotically critical nonlocal elliptic equations on conformal infinities, J. Differ. Equations 265 (2018), no. 19, 5115–5155.
  • [5] R Edwards, Fourier series. a modern introduction, Springer-Verlag, Heidelberg, 1982.
  • [6] G Eskin, Boundary value problems for elliptic pseudodifferential equations, AMS, Providence, 1981.
  • [7] L Frank, Spaces of network functions, Math. USSR Sb. 15 (1971), no. 2, 182–226.
  • [8] F.D Gakhov, Boundary value problems, Dover Publications, Mineola, 1981.
  • [9] L Hoermander, Analysis of partial differential operators, vol. i–iv, Springer-Verlag, Berlin, 1983.
  • [10] N.I Muskhelishvili, Singular integral equations, North Holland, Amsterdam, 1976.
  • [11] V Rabinovich, Wiener algebra of operators on the lattice μ𝐳n\mu{\bf z}^{n} depending on the small parameter μ>0\mu>0, Complex Var. Ell. Equ. 58 (2013), no. 6, 751–766.
  • [12] V Ryaben’kii, Method of difference potentials and its applications, Springer-Verlag, Berlin–Heidelberg, 2002.
  • [13] A Samarskii, The theory of difference schemes, CRC Press, Boca Raton, 2001.
  • [14] F.-O Speck, From sommerfeld diffraction problems to operator factorisation, Constr. Math. Anal.. 2 (2019), no. 4, 183–216.
  • [15] O. A. Tarasova and V. B Vasilyev, To the theory of discrete boundary value problems, 4open 2 (2019), no. 17, 1–7.
  • [16] M Taylor, Pseudodifferential operators, Princeton University Press, Princeton, 1981.
  • [17] F Treves, Introduction to pseudodifferential operators and fourier integral operators, Springer, New York, 1980.
  • [18] V. B Vasil’ev, Wave factorization of elliptic symbols: Theory and applications, Kluwer Academic Publishers, Dordrecht–Boston–London, 2000.
  • [19] by same author, On some new boundary-value problems in nonsmooth domains, J. Math. Sci. 173 (2011), no. 2, 225–230.
  • [20] A. Vasilyev and V Vasilyev, Discrete singular operators and equations in a half-space, Azerb. J. Math. 3 (2013), no. 1, 84–93.
  • [21] A. V. Vasilyev and V. B Vasilyev, Pseudo-differential operators and equations in a discrete half-space, Math. Model. Anal. 23 (2018), no. 3, 492–506.
  • [22] V Vasilyev, Discrete equations and periodic wave factorization, Proceedings of the International Conference on Analysis and Applied Mathematics (ICAAM-2016) (Melville) (A. Ashyralyev and A. Lukashov, eds.), AIP Conf. Proc., vol. 1759, AIP Publishing, 2016, p. 020126.
  • [23] by same author, Discreteness, periodicity, holomorphy, and factorization, Integral Methods in Science and Engineering (New York) (C. Constanda, M. Dalla Riva, P.D. Lamberti, and P. Musolino, eds.), Theoretical Technique, vol. 1, Birkhauser, 2017, pp. 315–324.
  • [24] by same author, On discrete boundary value problems, Proceedings of the International Conference “Functional Analysis in Interdisciplinary Applications” (FAIA2017) (Melville) (T. Kal’menov and M. Sadybekov, eds.), AIP Conf. Proc., vol. 1880, AIP Publishing, 2017, p. 050010.
  • [25] by same author, The periodic cauchy kernel, the periodic bochner kernel, discrete pseudo-differential operators, Proceedings of the International Conference on Numerical Analysis and Applications (ICNAAM-2016) (Melville) (T. Simos and C. Tsitouras, eds.), AIP Conf. Proc., vol. 1863, AIP Publishing, 2017, p. 140014.
  • [26] V.S Vladimirov, Generalized functions in mathematical physics, Mir, Moscow, 1979.