Research

Boundary value problems with eigenvalue-dependent boundary and transmission conditions

Author Affiliations

Department of Mathematics, Faculty of Arts and Science, Gaziosmanpaşa University, Tokat, 60250, Turkey

Boundary Value Problems 2014, 2014:131  doi:10.1186/1687-2770-2014-131

The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2014/1/131

 Received: 4 March 2014 Accepted: 8 May 2014 Published: 22 May 2014

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Abstract

In this paper the operator-theoretical method to investigate a new type boundary value problems consisting of a two-interval Sturm-Liouville equation together with boundary and transmission conditions dependent on eigenparameter is developed. By suggesting our own approach, we construct modified Hilbert spaces and a linear operator in them in such a way that the considered problem can be interpreted as a spectral problem for this operator. Then we introduce so-called left- and right-definite solutions and give a representation of solution of the corresponding nonhomogeneous problem in terms of these one-hand solutions. Finally, we construct Green’s vector-function and investigate some important properties of the resolvent operator by using this Green’s vector-function.

Keywords:
Sturm-Liouville problems; eigenparameter-dependent boundary and transmission conditions; Green’s function; resolvent operator

1 Introduction

Many important special equations which appear in physics, such as Airy’s equation, Bessel’s equation, wave equation, heat equation, Schrödinger’s equation, Heun’s equation, advection-dispersion equation, etc., are associated with Sturm-Liouville type operators. For instance, the one-dimensional form of the advection-dispersion equation for a nonreactive dissolved solute in a saturated, homogeneous, isotropic porous medium under steady, uniform flow is

c t + ν c x = D c x x , 0 < x < L , t > 0 ,

where c ( x , t ) is the concentration of the solute, ν is the average linear groundwater velocity, D is the coefficient of hydrodynamic dispersion, and L is the length of the aquifer. Using the method of separation of variables, the problem can be written in the simplest Sturm-Liouville form

[ p ( x ) X ] + λ r ( x ) X = 0 , X ( 0 ) = 0 , X ( L ) = 0 , 0 < x < L .

This example makes it clear that the Sturm-Liouville problems are of broad interest. There is a well-developed theory for classical Sturm-Liouville problems (see, e.g., [1-5] and the references therein). Details of the derivation of the theory and of related background results can be found in the cited references. Although the subject of Sturm-Liouville problems is over 160 years old, these problems are an intensely active field of research today. The main tool for solvability analysis of such problems is the concept of Green’s function. Green’s functions have played an important role as a theoretical tool in the field of physics, since the possibility of a transition from the problems in mathematical physics to integral equations is based on the fundamental concept of Green’s function. Therefore, the powerful and unifying formalism of Green’s functions finds applications not only in standard physics subjects such as perturbation and scattering theory, bound-state formation, etc., but also at the forefront of current and, most likely, future developments (see [6]). Green’s function transforms the differential equation into the integral equation, which, at times, is more informative. In terms of Green’s function, the BVP with arbitrary data can be solved in a form that shows clearly the dependence of the solution on the data. Namely, Green’s function approach would allow us to have an integral representation of the solution instead of an infinite series. Determination of Green’s functions is also possible using Sturm-Liouville theory. This leads to series representation of Green’s functions (see, e.g., the monograph [1] as well as the recent results in [7] and the references therein).

Sturm-Liouville type problems with transmission conditions have become an important area of research in recent years because of the needs of modern technology, engineering and physics. Many of the mathematical problems encountered in the study of boundary-value-transmission problem cannot be treated with the usual techniques within the standard framework of boundary value problem (see [8-12]). In this study we shall consider a new type of Sturm-Liouville problems consisting of the two-interval Sturm-Liouville equation

L ( u ) : = ( p ( x ) u ( x ) ) + q ( x ) u ( x ) = μ 2 u ( x ) , x Ω Ω + (1)

together with eigenparameter-dependent boundary conditions of the form

1 ( u ) : = δ 10 u ( a ) δ 11 u ( a ) μ 2 ( δ ˜ 10 u ( a ) δ ˜ 11 u ( a ) ) = 0 , (2)

2 ( u ) : = δ 20 u ( b ) δ 21 u ( b ) + μ 2 ( δ ˜ 20 u ( b ) δ ˜ 21 u ( b ) ) = 0 , (3)

and eigenparameter-dependent transmission conditions at one interaction point x = c of the form

3 ( u ) : = γ 10 + u ( c + ) + γ 10 u ( c ) = 0 , (4)

4 ( u ) : = γ 20 + u ( c + ) + γ 21 + u ( c + ) + γ 20 u ( c ) + γ 21 u ( c ) 4 ( u ) : = μ 2 ( γ ˜ 20 + u ( c + ) + γ ˜ 21 + u ( c + ) + γ ˜ 20 u ( c ) + γ ˜ 21 u ( c ) ) = 0 , (5)

where p ( x ) is a real-valued piecewise constant function, p ( x ) = p > 0 for x Ω = [ a , c ) , p ( x ) = p + > 0 for x Ω + = ( c , b ] , the potential q ( x ) is a real-valued function continuous in each of the intervals Ω and Ω + , and has finite limits q ( c ± ) = lim x c ± 0 q ( x ) , μ is a complex spectral parameter, the coefficients δ i j , δ ˜ i j , γ 2 j ± , γ ˜ 2 j ± ( i = 1 , 2 and j = 0 , 1 ), γ 10 ± are real numbers. This Sturm-Liouville problem is a non-classical eigenvalue problem since the eigenvalue parameter μ appears not only in the differential equation, but also in the boundary and transmission conditions. Moreover, in the differential equation there appears a singularity at one interior point. Because of these reasons the spectral theory of this problem is more complicate. Naturally, eigenfunctions of this problem may have discontinuity at the singular interior point. Some special cases of this problem arise after an application of the Fourier method to a varied assortment of physical problems. For instance, some boundary value problems with transmission conditions arise in heat and mass transfer problems [13], in vibrating string problems when the string is loaded additionally with point masses [14], in diffraction problems [12], in quantum mechanics [15], in thermal conduction problems for a thin laminated plate [16]etc. Such properties as isomorphism, coerciveness with respect to the spectral parameter, completeness and Abel bases property of a system of root functions of some boundary value problems with transmission conditions and its applications to the corresponding initial boundary value problems for parabolic equations have been investigated in [16-19]. For the background and applications of boundary value transmission problems to different areas, we refer the reader to the monographs and some recent contribution [8-11,17,18,20-25].

2 Hilbert space formulation of the problem

In certain cases the boundary value problem can be characterized by means of a uniquely determined unbounded self-adjoint operator. In these cases the eigenvalues and eigenfunctions of the boundary value problem are determined by the eigenvalues and eigenvectors of the corresponding operator; these will be called a self-adjoint case of the boundary value problem. In some cases such a characterization is not possible and these will be referred to as ‘symmetric’ cases in general. In classical point of view, our problem cannot be characterized as ‘self-adjoint case’. For ‘self-adjoint characterization’ of the considered problem (1)-(5), we shall define a new Hilbert space as follows.

Denote the determinant of the ith and jth columns of the matrix

θ = [ γ 20 + γ 21 + γ 20 γ 21 γ ˜ 20 + γ ˜ 21 + γ ˜ 20 γ ˜ 21 ]

by | θ i j | ( 1 i < j 4 ). Throughout the paper we shall assume that the conditions

γ 10 ( γ 10 | θ 12 | γ 10 + | θ 23 | ) > 0 , γ 10 + ( γ 10 + | θ 34 | γ 10 | θ 14 | ) < 0 , | θ 24 | = 0 , ρ 1 = det ( δ 11 δ 10 δ ˜ 11 δ ˜ 10 ) > 0 , ρ 2 = det ( δ 21 δ 20 δ ˜ 21 δ ˜ 20 ) > 0

hold. Define a new inner-product space ℋ as a direct sum space ( L 2 ( Ω ) L 2 ( Ω + ) ) C 3 equipped with the modified inner-product

U , V H : = 1 p ( | θ 34 | γ 10 γ 10 + | θ 14 | ) a c u ( x ) v ( x ) ¯ d x + 1 p + ( | θ 12 | γ 10 + γ 10 | θ 23 | ) c + b u ( x ) v ( x ) ¯ d x 1 ρ 1 ( | θ 34 | γ 10 γ 10 + | θ 14 | ) u 1 v ¯ 1 + 1 ρ 2 ( | θ 12 | γ 10 + γ 10 | θ 23 | ) u 2 v ¯ 2 + u 3 v ¯ 3 (6)

for U = ( u ( x ) , u 1 , u 2 , u 3 ) , V = ( v ( x ) , v 1 , v 2 , v 3 ) ( L 2 ( Ω ) L 2 ( Ω + ) ) C 3 . It is easy to see that the relation (6) really defines a new inner product in the direct sum space ( L 2 ( Ω ) L 2 ( Ω + ) ) C 3 .

Lemma 1is a Hilbert space.

Proof Let U n = ( u n ( x ) , u 1 n , u 2 n , u 3 n ) , n = 1 , 2 ,  , be any Cauchy sequence in ℋ. Then by (6) the sequences ( u n ( x ) ) and ( u 1 n , u 2 n , u 3 n ) will be Cauchy sequences in the Hilbert spaces L 2 ( Ω ) L 2 ( Ω + ) and C 3 , respectively. Therefore they are convergent. Let u 0 ( x ) and ( u 10 , u 20 , u 30 ) be limits of these sequences, respectively. Defining U 0 = ( u 0 ( x ) , u 10 , u 20 , u 30 ) we have that U 0 H and U n U 0 in ℋ. The proof is complete. □

Let us now define the boundary and transmission functionals B a ( u ) : = δ 10 u ( a ) δ 11 u ( a ) , B ˜ a ( u ) : = δ ˜ 10 u ( a ) δ ˜ 11 u ( a ) , B b ( u ) : = δ 20 u ( b ) δ 21 u ( b ) , B ˜ b ( u ) : = δ ˜ 20 u ( b ) δ ˜ 21 u ( b ) , T c ( u ) : = T c ( u ) T c + ( u ) , T ˜ c ( u ) : = T ˜ c ( u ) T ˜ c + ( u ) , where T c ± ( u ) = γ 20 ± u ( c ± ) + γ 21 ± u ( c ± ) , T ˜ c ± ( u ) = γ ˜ 20 ± u ( c ± ) + γ ˜ 21 ± u ( c ± ) and the linear operator : H H with the domain

dom ( ) : = { U = ( u ( x ) , u 1 , u 2 , u 3 ) : u ( x ) , u ( x ) A C loc ( Ω ) A C loc ( Ω + ) , with a finite left- and right-hand limits  u ( c 0 )  and  u ( c 0 ) ; L ( u ) L 2 [ a , b ] ; 3 ( u ) = 0 ; u 1 = B ˜ a ( u ) ; u 2 = B ˜ b ( u ) ; u 3 = T ˜ c ( u ) }

and action low

( u ( x ) , u 1 , u 2 , u 3 ) = ( L f , B a ( u ) , B b ( u ) , T c ( u ) ) .

Then problem (1)-(5) can be written in the operator equation form as U = μ 2 U , U dom ( ) in the Hilbert space ℋ.

Theorem 1The linear operatoris symmetric in the Hilbert space ℋ.

Proof By applying the method of [22] it is not difficult to prove that the operator ℜ is densely defined in ℋ, i.e., dom ( ) ¯ = H . Now let U = ( u ( x ) , B ˜ a ( u ) , B ˜ b ( u ) , T ˜ c ( u ) T ˜ c + ( u ) ) , V = ( v ( x ) , B ˜ a ( v ) , B ˜ b ( v ) , T ˜ c ( v ) T ˜ c + ( v ) ) dom ( ) . By partial integration we have

U , V H U , V H = ( | θ 34 | γ 10 γ 10 + | θ 14 | ) W ( u , v ¯ ; a ) ( | θ 34 | γ 10 γ 10 + | θ 14 | ) W ( u , v ¯ ; c ) + ( | θ 12 | γ 10 + γ 10 | θ 23 | ) W ( u , v ¯ ; b ) ( | θ 12 | γ 10 + γ 10 | θ 23 | ) W ( u , v ¯ ; c + ) 1 ρ 1 ( | θ 34 | γ 10 γ 10 + | θ 14 | ) ( B a ( u ) B ˜ a ( v ) ¯ B ˜ a ( u ) B a ( v ) ¯ ) + 1 ρ 2 ( | θ 12 | γ 10 + γ 10 | θ 23 | ) ( B ˜ b ( u ) B b ( v ) ¯ B b ( u ) B ˜ b ( v ) ¯ ) + ( T c ( u ) T ˜ c ( v ¯ ) T ˜ c ( u ) T c ( v ¯ ) ) , (7)

where, as usual, W ( u , v ¯ ; x ) denotes the Wronskians of the functions u and v ¯ . From the definitions of boundary functionals we get that

B a ( u ) B ˜ a ( v ) ¯ B ˜ a ( u ) B a ( v ) ¯ = ρ 1 W ( u , v ¯ ; a ) , (8)

B ˜ b ( u ) B b ( v ) ¯ B b ( u ) B ˜ b ( v ) ¯ = ρ 2 W ( u , v ¯ ; b ) . (9)

Further, taking in view the definition of ℜ and initial conditions (14)-(19) we can derive that

T c ( u ) T ˜ c ( v ¯ ) T ˜ c ( u ) T c ( v ¯ ) = ( | θ 12 | γ 10 + γ 10 | θ 23 | ) W ( u , v ¯ ; c + ) + ( | θ 34 | γ 10 γ 10 + | θ 14 | ) W ( u , v ¯ ; c ) . (10)

Finally, substituting (8), (9) and (10) in (7) we obtain that

U , V H = U , V H for all  U , V dom ( ) .

The proof is complete. □

Theorem 2The linear operatoris self-adjoint in ℋ.

Proof Since ℜ is symmetric and densely defined on ℋ, it is sufficient to show that if

U , V H = U , W H (11)

for all U = ( u ( x ) , B ˜ a ( u ) , B ˜ b ( u ) , T ˜ c ( u ) T ˜ c + ( u ) ) D ( ) , then V D ( ) and V = W , where V = ( v ( x ) , ħ 1 , ħ 2 , ħ 3 ) and W = ( w ( x ) , κ 1 , κ 2 , κ 3 ) . Writing equality (11) for all U ( C 0 ( Ω ) C 0 ( Ω + ) ) C 3 D ( ) by standard Sturm-Liouville theory, we find that u , v L 2 = u , w L 2 and w ( x ) = L v . Then from equality (11) it follows that

L u , v L 2 = u , L v L 2 1 ρ 1 ( | θ 34 | γ 10 γ 10 + | θ 14 | ) ( B ˜ a ( u ) κ ¯ 1 B a ( u ) ħ ¯ 1 ) 1 ρ 2 ( | θ 12 | γ 10 + γ 10 | θ 23 | ) ( B ˜ b ( u ) κ ¯ 2 + B b ( u ) ħ ¯ 2 ) + [ ( T ˜ c ( u ) T ˜ c + ( u ) ) κ ¯ 3 ( T c ( u ) T c + ( u ) ) ħ ¯ 3 ] .

On the other hand, by two partial integrations we get

L u , v L 2 = u , L v L 2 + ( | θ 34 | γ 10 γ 10 + | θ 14 | ) W ( u , v ¯ ; a ) ( | θ 34 | γ 10 γ 10 + | θ 14 | ) W ( u , v ¯ ; c ) + ( | θ 12 | γ 10 + γ 10 | θ 23 | ) W ( u , v ¯ ; b ) ( | θ 12 | γ 10 + γ 10 | θ 23 | ) W ( u , v ¯ ; c + ) . (12)

Thus,

1 ρ 1 ( | θ 34 | γ 10 γ 10 + | θ 14 | ) ( B ˜ a ( u ) κ ¯ 1 B a ( u ) ħ ¯ 1 ) 1 ρ 2 ( | θ 12 | γ 10 + γ 10 | θ 23 | ) ( B ˜ b ( u ) κ ¯ 2 + B b ( u ) ħ ¯ 2 ) + [ ( T ˜ c ( u ) T ˜ c + ( u ) ) κ ¯ 3 ( T c ( u ) T c + ( u ) ) ħ ¯ 3 ] = ( | θ 34 | γ 10 γ 10 + | θ 14 | ) W ( u , v ¯ ; a ) ( | θ 34 | γ 10 γ 10 + | θ 14 | ) W ( u , v ¯ ; c ) + ( | θ 12 | γ 10 + γ 10 | θ 23 | ) W ( u , v ¯ ; b ) ( | θ 12 | γ 10 + γ 10 | θ 23 | ) W ( u , v ¯ ; c + ) . (13)

From this equality, by applying the technique of Theorem 2.5 in our previous work [11], it can be derived easily that ħ 1 = B ˜ a ( v ) , ħ 2 = B ˜ b ( v ) , ħ 3 = T ˜ c ( v ) T ˜ c + ( v ) , 3 v = 0 and κ 1 = B a ( w ) , κ 2 = B b ( w ) , κ 3 = T c ( w ) T c + ( w ) . The proof is complete. □

Theorem 3The operatorhas only point spectrum, i.e., σ ( ) = σ p ( ) .

Proof It suffices to prove that if μ 2 is not an eigenvalue of ℜ, then μ 2 is a regular point of ℜ, i.e., μ 2 σ ( ) . Let μ 2 not be an eigenvalue of ℜ. The resolvent operator R ( μ 2 , ) = ( μ 2 I ) 1 exists and is defined on all of ℋ. By Theorem 2 and the closed graph theorem, we get that R ( μ 2 , ) is bounded. Thus, μ 2 σ ( ) . Hence σ ( ) = σ p ( ) . □

3 Left-definite and right-definite solutions

In this section we shall define two basic solutions υ ( x , μ ) and ω ( x , μ ) on the left interval Ω (so-called left-definite solutions) and two basic solutions υ + ( x , μ ) and ω + ( x , μ ) on the right interval Ω + (so-called right-definite solutions) by a special procedure as follows. Let υ ( x , μ ) and ω + ( x , μ ) be solutions of equation (1) on Ω and Ω + satisfying the initial conditions

υ ( a , μ ) = δ 11 μ 2 δ ˜ 11 , υ ( a , μ ) x = δ 10 μ 2 δ ˜ 10 (14)

and

ω + ( b , μ ) = δ 21 + μ 2 δ ˜ 21 , ω + ( b , μ ) x = δ 20 + μ 2 δ ˜ 20 , (15)

respectively. By using these solutions we shall define the other solutions υ + ( x , μ ) and ω ( x , μ ) by the initial conditions

υ + ( c + , μ ) = γ 10 γ 10 + υ ( c , μ ) , (16)

υ + ( c + , μ ) x = 1 γ 10 + ( γ 21 + μ 2 γ ˜ 21 + ) [ γ 10 + ( μ 2 γ ˜ 21 γ 21 ) υ ( c , μ ) x υ + ( c + , μ ) x = + ( γ 10 + ( μ 2 γ ˜ 20 γ 20 ) γ 10 ( μ 2 γ ˜ 20 + γ 20 + ) ) υ ( c , μ ) ] (17)

and

ω ( c , μ ) = γ 10 + γ 10 ω + ( c + , μ ) , (18)

ω ( c , μ ) x = 1 γ 10 ( γ 21 μ 2 γ ˜ 21 ) [ γ 10 ( μ 2 γ ˜ 21 + γ 21 + ) ω + ( c + , μ ) x ω ( c , μ ) x = + ( γ 10 ( μ 2 γ ˜ 20 + γ 20 + ) γ 10 + ( μ 2 γ ˜ 20 γ 20 ) ) ω + ( c + , μ ) ] , (19)

respectively. The existence of these solutions follows from the well-known Cauchy-Picard theorem of ordinary differential equation theory. Moreover, by applying the method of [20], we can prove that each of these solutions are entire functions of the parameter μ C for each fixed x.

4 Construction of Green’s function

In this section we develop the idea of a resolvent operator to solve nonhomogeneous boundary-value transmission problems (BVTP) as follows. Consider the operator equation

( μ 2 I ) U = F (20)

for arbitrary F = ( f ( x ) , f 1 , f 2 , f 3 ) H . This operator equation is equivalent to the following nonhomogeneous BVTP:

( μ 2 L ) u ( x ) = f ( x ) , x Ω Ω + , (21)

1 ( u ) = f 1 , 2 ( u ) = f 2 , 3 ( u ) = 0 , 4 ( u ) = f 3 . (22)

Let us define the Wronskians Δ ( μ ) : = W [ υ , ω ] x for x Ω and Δ + ( μ ) : = W [ υ + , ω + ] x for x Ω + and suppose that ( γ 21 ± μ 2 γ ˜ 21 ± ) Δ ± ( μ ) 0 . We shall search the resolvent function of this BVTP in the form

u ( x , μ ) = { ħ 11 ( x , μ ) υ ( x , μ ) + ħ 12 ( x , μ ) ω ( x , μ ) for  x Ω , ħ 21 ( x , μ ) υ + ( x , μ ) + ħ 22 ( x , μ ) ω + ( x , μ ) for  x Ω + , (23)

where the functions ħ 11 ( x , μ ) and ħ 12 ( x , μ ) are the solutions of the system of equations

{ ħ 11 ( x , μ ) x υ ( x , μ ) + ħ 12 ( x , μ ) x ω ( x , μ ) = 0 , ħ 11 ( x , μ ) x υ ( x , μ ) x + ħ 12 ( x , μ ) x ω ( x , μ ) x = f ( x ) p (24)

and the functions ħ 21 ( x , μ ) , ħ 22 ( x , μ ) are the solutions of the system of equations

{ ħ 21 ( x , μ ) x υ + ( x , μ ) + d ħ 22 ( x , μ ) x ω + ( x , μ ) = 0 , ħ 21 ( x , μ ) x υ + ( x , μ ) x + ħ 22 ( x , μ ) x ω + ( x , μ ) x = f ( x ) p + (25)

for x Ω and x Ω + , respectively. Since Δ ( μ ) 0 and Δ + ( μ ) 0 , from (24) and (25) we have

ħ 11 ( x , μ ) = 1 p Δ ( μ ) x c f ( s ) ω ( s , μ ) d s + ħ ˜ 11 ( μ ) , x Ω , ħ 12 ( x , μ ) = 1 p Δ ( μ ) a x f ( s ) υ ( s , μ ) d s + ħ ˜ 12 ( μ ) , x Ω , ħ 21 ( x , μ ) = 1 p + Δ + ( μ ) x b f ( s ) ω + ( s , μ ) d s + ħ ˜ 21 ( μ ) , x Ω + , ħ 22 ( x , μ ) = 1 p + Δ + ( μ ) c + x f ( s ) υ + ( s , μ ) d s + ħ ˜ 22 ( μ ) , x Ω + ,

where ħ ˜ i j ( μ ) ( i , j = 1 , 2 ) are unknown functions depending only on the parameter μ. Substituting into (23) gives

u ( x , μ ) = { ω ( x , μ ) p Δ ( μ ) a x υ ( s , μ ) f ( s ) d s + υ ( x , μ ) p Δ ( μ ) x c ω ( s , μ ) f ( s ) d s + ħ ˜ 11 ( μ ) υ ( x , μ ) + ħ ˜ 12 ( μ ) ω ( x , μ ) for  x Ω , ω + ( x , μ ) p + Δ + ( μ ) c + x υ + ( s , μ ) f ( s ) d s + υ + ( x , μ ) p + Δ + ( μ ) x b ω + ( s , μ ) f ( s ) d s + ħ ˜ 21 ( μ ) υ + ( x , μ ) + ħ ˜ 22 ( μ ) ω + ( x , μ ) for  x Ω + . (26)

By differentiating we have

u ( x , μ ) x = { 1 p Δ ( μ ) ω ( x , μ ) x a x υ ( s , μ ) f ( s ) d s + 1 p Δ ( μ ) υ ( x , μ ) x x c ω ( s , μ ) f ( s ) d s + ħ ˜ 11 ( μ ) υ ( x , μ ) x + ħ ˜ 12 ( μ ) ω ( x , μ ) x for  x Ω , 1 p + Δ + ( μ ) ω + ( x , μ ) x c + x υ + ( s , μ ) f ( s ) d s + 1 p + Δ + ( μ ) υ + ( x , μ ) x x b ω + ( s , μ ) f ( s ) d s + ħ ˜ 21 ( μ ) υ + ( x , μ ) x + ħ ˜ 22 ( μ ) ω + ( x , μ ) x for  x Ω + . (27)

By using (26), (27) and conditions (22) we can derive that

ħ ˜ 11 ( μ ) = 1 p + Δ + ( μ ) c + b ω + ( s , μ ) f ( s ) d s + f 2 Δ + ( μ ) + f 3 γ 10 + ( γ 21 + μ 2 γ ˜ 21 + ) Δ + ( μ ) , ħ ˜ 12 ( μ ) = f 1 Δ ( μ ) , ħ ˜ 21 ( μ ) = f 2 Δ + ( μ )

and

ħ ˜ 22 ( μ ) = 1 p Δ ( μ ) a c υ ( s , μ ) f ( s ) d s + f 1 Δ ( μ ) + f 3 γ 10 ( γ 21 μ 2 γ ˜ 21 ) Δ ( μ ) .

Putting in (26) gives

u ( x , μ ) = { ω ( x , μ ) p Δ ( μ ) a x υ ( s , μ ) f ( s ) d s + υ ( x , μ ) p Δ ( μ ) x c ω ( s , μ ) f ( s ) d s + υ ( x , μ ) ( 1 p + Δ + ( μ ) c + b ω + ( s , μ ) f ( s ) d s + f 2 Δ + ( μ ) + f 3 γ 10 + ( γ 21 + μ 2 γ ˜ 21 + ) Δ + ( μ ) ) + ω ( x , μ ) f 1 Δ ( μ ) for  x Ω , ω + ( x , μ ) p + Δ + ( μ ) c + x υ 2 ( s , μ ) f ( s ) d s + υ + ( x , μ ) p + Δ + ( μ ) x b ω + ( s , μ ) f ( s ) d s + υ + ( x , μ ) f 2 Δ + ( μ ) + ω + ( x , μ ) ( 1 p Δ ( μ ) a c υ ( s , μ ) f ( s ) d s + f 1 Δ ( μ ) + f 3 γ 10 ( γ 21 μ 2 γ ˜ 21 ) Δ ( μ ) ) for  x Ω + . (28)

Thus we find the needed resolvent function u ( x , μ ) in terms of the left- and right-define solutions υ ± ( x , μ ) and ω ± ( x , μ ) . By introducing Green’s function as

G ( x , s ; μ ) = { υ ( x , μ ) ω ( s , μ ) p Δ ( μ ) if  x , s Ω , s x , ω ( x , μ ) υ ( s , μ ) p Δ ( μ ) if  x , s Ω , s < x , υ ( x , μ ) ω + ( s , μ ) p Δ ( μ ) if  x Ω , s Ω + , ω ( x , μ ) υ + ( s , μ ) p + Δ + ( μ ) if  x Ω + , s Ω , ω + ( x , μ ) υ + ( s , μ ) p + Δ + ( μ ) if  x Ω + , s x , υ + ( x , μ ) ω + ( s , μ ) p + Δ + ( μ ) if  x Ω + , s > x (29)

from (28) and (29) we have that the considered problem (21)-(22) has a unique solution given by

u ( x , μ ) = 1 p ( | θ 34 | γ 10 γ 10 + | θ 14 | ) a c G ( x , s ; μ ) f ( s ) d s + 1 p + ( | θ 12 | γ 10 + γ 10 | θ 23 | ) c + b G ( x , s ; μ ) f ( s ) d s + 1 ρ 1 ( | θ 34 | γ 10 γ 10 + | θ 14 | ) B ˜ a ( G ( x , ; μ ) ) f 1 1 ρ 2 ( | θ 12 | γ 10 + γ 10 | θ 23 | ) B ˜ b ( G ( x , ; μ ) ) f 2 + T ˜ c ( G ( x , ; μ ) ) f 3 . (30)

5 Representations of the resolvent operator in terms of Green’s vector-function

We now shall define Green’s vector-function as follows:

G ˜ x , μ : = ( G ( x , ; μ ) , B ˜ a ( G ( x , ; μ ) ) , B ˜ b ( G ( x , ; μ ) ) , T ˜ c ( G ( x , ; μ ) ) ) . (31)

Consequently, for the solution U ( F , μ ) of nonhomogeneous operator equation (21), we obtain the following formula:

U ( F , μ ) = ( G ˜ x , μ , F ¯ , B ˜ a ( G ˜ x , μ , F ¯ ) , B ˜ b ( G ˜ x , μ , F ¯ ) , T ˜ c ( G ˜ x , μ , F ¯ ) ) . (32)

Using this, the resolvent function (30) can be written in the form

u ( x , μ ) = G ˜ x , μ , F ¯ H , (33)

where F ¯ = ( f ¯ ( x ) , f ¯ 1 , f ¯ 2 , f ¯ 3 ) H . Consequently, we have the following theorem.

Theorem 4For the resolvent operator R ( μ 2 , ) = ( μ 2 I ) 1 , the formula

R ( μ 2 , ) F = ( G ˜ x , μ , F ¯ , B ˜ a ( G ˜ x , μ , F ¯ ) , B ˜ b ( G ˜ x , μ , F ¯ ) , T ˜ c ( G ˜ x , μ , F ¯ ) ) , (34)

where F ¯ = ( f ¯ ( x ) , f ¯ 1 , f ¯ 2 , f ¯ 3 ) H , holds.

Theorem 5The estimation

R ( μ 2 , ) F H | Im μ 2 | 1 F H , F H

holds for all regular value μ 2 such that Im μ 2 0 .

Proof Let F = ( f ( x ) , f 1 , f 2 , f 3 ) H . Denote U = R ( μ 2 , ) F . Since U = μ 2 U F , taking into account that the operator ℜ is symmetric, we have

μ 2 U , U H F , U H = U , U H = U , U H = μ ¯ 2 U , U H F , U ¯ H .

Using the well-known Cauchy-Schwarz inequality, we conclude that

| Im μ 2 | U H 2 = | Im F , U | H F H U H .

Consequently,

R ( μ 2 , ) F H | Im μ 2 | 1 F H .

The proof is complete. □

Theorem 6The resolvent operator R ( μ 2 , ) is compact in the Hilbert space ℋ.

Proof Let λ 1 λ 2 be eigenvalues of ℜ and let E λ 1 , E λ 2 , be orthogonal projections onto the corresponding eigenspace. Since ℜ is a self-adjoint operator with discrete spectrum, we can write the spectral resolution of the resolvent operator R ( μ 2 , ) by

R ( μ 2 , ) = n = 1 ( λ n μ 2 ) 1 E λ n . (35)

Similarly to [22] we can easily show that λ n = O ( n 2 ) . Thus ( λ n μ 2 ) 1 = O ( n 2 ) for n . Consequently, the series (35) is strongly convergent. It is obvious that the orthogonal projections E λ n , n = 1 , 2 ,  , are compact operators since each of them are of finite rank. Consequently, the sum of series (35) is also compact in ℋ. The proof is complete. □

Competing interests

The author declares that she has no competing interests.

Acknowledgements

The author is grateful to anonymous referees for their constructive comments and suggestions, which led to the improvement of the original manuscript.

References

1. Levitan, BM, Sargsyan, IS: Sturm-Liouville and Dirac Operators, Springer, New York (1991)

2. Pryce, JD: Numerical Solution of Sturm-Liouville Problems, Oxford University Press, New York (1993) Oxford Science Publications

3. Hinton, D, Schaefer, PW: Spectral Theory and Computational Methods for Sturm-Liouville Problems, Dekker, New York (1997)

4. Titchmarsh, EC: Eigenfunction Expansions Associated with Second-Order Differential Equations. Part I, Oxford University Press, London (1962)

5. Zettl, A: Sturm-Liouville Theory, Am. Math. Soc., Providence (2005)

6. Duffy, DG: Green’s Functions with Applications, Chapman & Hall/CRC, Boca Raton (2001).

7. Stakgold, I, Holst, MJ: Green’s Functions and Boundary Value Problems, Wiley, New York (2011).

8. Ao, J, Sun, J, Zhang, M: Matrix representations of Sturm-Liouville problems with transmission conditions. Comput. Math. Appl.. 63, 1335–1348 (2012). Publisher Full Text

9. Aydemir, K, Mukhtarov, OS: Green’s function method for self-adjoint realization of boundary-value problems with interior singularities. Abstr. Appl. Anal.. 2013, Article ID 503267 (2013)

Article ID 503267

Publisher Full Text

10. Kong, Q, Wang, Q: Using time scales to study multi-interval Sturm-Liouville problems with interface conditions. Results Math.. 63, 451–465 (2013). Publisher Full Text

11. Mukhtarov, OS, Aydemir, K: New type Sturm-Liouville problems in associated Hilbert spaces. J. Funct. Spaces Appl.. 2014, Article ID 606815 (2014)

12. Voitovich, NN, Katsenelbaum, BZ, Sivov, AN: Generalized Method of Eigen-Vibration in the Theory of Diffraction, Nauka, Moscow (1997) (Russian)

13. Likov, AV, Mikhailov, YA: The Theory of Heat and Mass Transfer, Qosenergaizdat, Moscow (1963) (Russian)

14. Tikhonov, AN, Samarskii, AA: Equations of Mathematical Physics, Pergamon, New York (1963)

15. Albeverio, S, Gesztesy, F, Hoegh Krohn, R, Holden, H: Solvable Models in Quantum Mechanics, AMS Chelsea Publishing, Providence (2005)

16. Titeux, I, Yakubov, Y: Completeness of root functions for thermal conduction in a strip with piecewise continuous coefficients. Math. Models Methods Appl. Sci.. 7, 1035–1050 (1997). Publisher Full Text

17. Mukhtarov, OS, Demir, H: Coerciveness of the discontinuous initial-boundary value problem for parabolic equations. Isr. J. Math.. 114, 239–252 (1999). Publisher Full Text

18. Mukhtarov, OS, Yakubov, S: Problems for ordinary differential equations with transmission conditions. Appl. Anal.. 81, 1033–1064 (2002). Publisher Full Text

19. Rasulov, ML: Methods of Contour Integration, North-Holland, Amsterdam (1967)

20. Akdoğan, Z, Demirci, M, Mukhtarov, OS: Green function of discontinuous boundary-value problem with transmission conditions. Math. Methods Appl. Sci.. 30, 1719–1738 (2007). Publisher Full Text

21. Bairamov, E: On the characteristic values of the real component of a dissipative boundary value transmission problem. Appl. Math. Comput.. 218, 9657–9663 (2012). Publisher Full Text

22. Mukhtarov, OS, Kadakal, M: Some spectral properties of one Sturm-Liouville type problem with discontinuous weight. Sib. Math. J.. 46, 681–694 (2005). Publisher Full Text

23. Altınışık, N, Mukhtarov, OS, Kadakal, M: Asymptotic formulas for eigenfunctions of the Sturm-Liouville problems with eigenvalue parameter in the boundary conditions. Kuwait J. Sci. Eng.. 39, 1–19 (2012)

24. Muhtarov, FS, Aydemir, K: Distributions of eigenvalues for Sturm-Liouville problem under jump conditions. J. New Results Sci.. 1, 81–89 (2012)

25. Uǧurlu, E, Bairamov, E: Dissipative operators with impulsive conditions. J. Math. Chem.. 51, 1670–1680 (2013). Publisher Full Text