arxiv: v1 [math.ca] 27 Mar 2013

Transcription

1 Boundary value problem with transmission conditions arxiv: v1 [math.ca] 27 Mar 213 K. Aydemir and O. Sh. Mukhtarov Department of Mathematics, Faculty of Science, Gaziosmanpaşa University, 625 Tokat, Turkey omukhtarov@yahoo.com, kadriye.aydemir@gop.edu.tr Abstract : One important innovation here is that for the Sturm-Liouville considered equation together with eigenparameter dependent boundary conditions and two supplementary transmission conditions at one interior point. We develop Green s function method for spectral analysis of the considered problem in modified Hilbert space. Keywords : Sturm-Liouville problems, Green s function. 1 Introduction In physics many problems arise in the form of boundary value problems involving second order ordinary differential equations. This derivation based on that in [6], however, is more thorough than that in most elementary physics texts; while most parameters such as density and other thermal properties are treated as constant in such treatments, the following allows fundamental properties of the bar to vary as a function of the bars length, which will lead to a Sturm-Liouville problem of a more general nature. In this study we shall investigate one discontinuous eigenvalue problem which consists of Sturm-Liouville equation, Γ(y) := y (x,λ)+q(x)y(x,λ) = λy(x,λ) (1.1) to hold in finite interval (,π) except at one inner point (,π), where discontinuity in u and u are prescribed by transmission conditions Γ 1 (y) := a 1 y (,λ)+a 2 y(,λ)+a 3 y (+,λ)+a 4 y(+,λ) =, (1.2) Γ 2 (y) := b 1 y (,λ)+b 2 y(,λ)+b 3 y (+,λ)+b 4 y(+,λ) =, (1.3) together with the boundary conditions Γ 3 (y) := cosαy(,λ)+sinαy (,λ) =, (1.4)

2 Γ 4 (y) := cosβy(π,λ)+sinβy (π,λ) =, (1.5) where the potential q(x) is real-valued, continuous in each interval [ 1, ) and (,1] and has a finite limits q( c) ; α,β,β 1 are real numbers; λ is a complex eigenparameter. In this study by using an own technique we introduce a new equivalent inner product in the Hilbert space L 2 ( 1,) L 2 (,1) and a linear operator in it such a way that the considered problem can be interpreted as eigenvalue problem for this operator. 2 Preliminary Results [ ] a1 a Let T= 2 a 3 a 4. Denote the determinant of the i-th and j-th b 1 b 2 b 3 b 4 columns of the matrix T by ρ ij. Note that throughout this study we shall assume that ρ 12 > and ρ 34 >. In this section we shall define two basic solutions φ(x,λ) and χ(x,λ) by own technique as follows. At first, let us consider solutions of the equation (1.1) on the left-hand [,) of the considered interval [π,) (,π] satisfying the initial conditions y(,λ) = sinα, y(, λ) x = cosα (2.1) By virtue of well-known existence and uniqueness theorem of ordinary differential equation theory this initial-value problem for each λ has a unique solution φ 1 (x,λ). Moreover [[12], Teorem 7] this solution is an entire function of λ for each fixed x [,). Using this solutions we can prove that the equation (1.1) on the right-hand interval (,π] of the considered interval [,) (,π] has the solution u = φ 2 (x,λ) satisfying the initial conditions y(,λ) = 1 φ 1 (,λ) (ρ 23 φ 1 (,λ)+ρ 24 ) ρ 12 x (2.2) y (,λ) = 1 φ 1 (,λ) (ρ 13 φ 1 (,λ)+ρ 14 ). ρ 12 x (2.3) By applying the method of [5] we can prove that the equation (1.1) on (,π] has an unique solution φ 2 (x,λ) satisfying the conditions (2.2)-(2.3) which also is an entire function of the parameter λ for each fixed x (,π]. Consequently, the function φ(x, λ) defined by φ(x,λ) = { φ 1(x,λ) for x [,) φ 2 (x,λ) for x (,π]. 2 (2.4)

3 satisfies equation (1.1), the first boundary condition (1.4) and the both transmission conditions (1.2) and (1.3). Similarly, χ 2 (x,λ) be solutions of equation (1.1) on the left-right interval (, π] subject to initial conditions y(π,λ) = sinβ, y(π, λ) x = cosβ. (2.5) By virtue of [[12], Teorem 7] each of these solutions are entire functions of λ for fixed x. By applying the same technique we can prove there is an unique solution χ 1 (x,λ) of equation (1.1) the left-hand interval [,) satisfying the initial condition y(,λ) = 1 χ 2 (,λ) (ρ 14 χ 2 (,λ)+ρ 24 ), ρ 34 x (2.6) y (,λ) = 1 χ 2 (,λ) (ρ 13 χ 2 +ρ 23 ). ρ 34 x (2.7) By applying the similar technique as in [5] we can prove that the solutions χ 1 (x,λ) are also an entire functions of parameter λ for each fixed x. Consequently, the function χ(x, λ) defined as { χ1 (x,λ), x [,) χ(x,λ) = χ 2 (x,λ), x (,π] satisfies the equation (1.1) on whole [,) (,π], the other boundary condition (1.5) and the both transmission conditions (1.2) and (1.3). In the Hilbert Space H = L 2 [ 1,) L 2 (,1] of two-component vectors we define an inner product by < y,z > H := ρ 12 y(x)z(x)dx+ρ 34 y(x)z(x)dx for y = y(x), z = z(x) H. We introduce the linear operator A : H H with domain of definition satisfying the following conditions i) y and y are absolutely continuous in each of intervals [,) and (,π] and has a finite limits y(c ) and y (c ) ii) Γy(x) H, Γ 1 y(x) = Γ 2 y(x) = Γ 3 y(x) = Γ 4 y(x) =, Obviously D(A) is a linear subset dense in H. We put (Ay)(x) = Γy(x), x H for y D(A). Then the problem (1.1) (1.5) is equivalent to the equation Ay = λy 3

4 in the Hilbert space H. Taking in view that the Wronskians W(φ i,χ i ;x) := φ i (x,λ)χ i (x,λ) φ i (x,λ)χ i(x,λ) are independent of variable x we shall denote w i (λ) = W(φ i,χ i ;x) (i = 1,2). By using (1.2) and (1.3) we have ρ 12 w 1 (λ) = ρ 34 w 2 (λ) for each λ C. It is convenient to introduce the notation w(λ) := ρ 34 w 1 (λ) = ρ 12 w 2 (λ). (2.8) Theorem 2.1. For all y,z D(A), the equality holds. < Ay,z >=< y,az > (2.9) Proof. Integrating by parts we have for all y,z D(A), < Ay,z > = ρ 12 = ρ 12 Γy(x)z(x)dx+ρ 34 Γy(x)z(x)dx y(x)γz(x)dx+ρ 34 y(x)γz(x)dx+ρ 12 W[y,z; ] ρ 12 W[y,z;]+ρ 34 W[y,z;π] ρ 34 W[y,z;] = < y,az > +ρ 12 W[y,z;] ρ 12 W[y,y ;] + ρ 34 W[y,y ;π] ρ 34 W[y,y ;] (2.1) From the boundary conditions (1.2)-(1.3) it is follows obviously that The direct calculation gives W(y,z;) = and W(y,z;π) = (2.11) ρ 12 W(y,z;) = ρ 34 W(y,z;). (2.12) Substituting (2.11) and (2.12) in (2.1) we obtain the equality (2.9). Relation (2.9) shows that the operator A is symmetric. Therefore all eigenvalues of the operator A are real and two eigenfunctions corresponding to the distinct eigenvalues are orthogonal in the sense of the following equality ρ 12 y(x)z(x)dx+ρ 34 y(x)z(x)dx =. (2.13) Theorem 2.2. D(A) is densely in the Hilbert space H. Proof. Theorem 2.3. The operator A is self-adjoint. Proof. 4

5 3 The Green s function Let us consider the inhomogeneous differential equation y +(λ q(x)y = f(x), x [,) (,π] (3.1) together with the boundary conditions (1.4)-(1.5) and the transmission conditions (1.2)-(1.3). Making use of the definitions of the functions φ i,χ i (i = 1, 2) we see that the general solution of the differential equation (3.1) can be represented in the form y(x,λ) = φ 1(ξ,λ)f(ξ)dξ + φ 1(x,λ) ω 1 (λ) +c 1 φ 1 (x,λ)+d 1 χ 1 (x,λ), for x [,) χ 1 (x,λ) ω 1 (λ) χ 2 (x,λ) ω 2 (λ) x χ 1(ξ,λ)f(ξ)dξ φ 2(ξ,λ)f(ξ)dξ + φ 2(x,λ) π ω 2 (λ) x χ 2(ξ,λ)f(ξ)dξ +c 2 φ 2 (x,λ)+d 2 χ 2 (x,λ), for x (,π] (3.2) where c i,d i (i = 1,2) are arbitrary constants. By substitution into the boundary conditions (1.4)-(1.5) we have at once that d 1 =, c 2 =. (3.3) Further, substitution (3.2) into transmission conditions (1.2)-(1.3) we have the inhomogeneous linear system of equations for c 1 and d 1, the determinant of which is equal to and therefore is not vanish by assumption. Solving that system we find c 1 = 1 ω 2 (λ) y f (x,λ) = χ 2 (ξ,λ)f(ξ)dξ, d 2 = 1 ω 1 (λ) φ 1 (ξz,λ)f(ξ)dξ (3.4) Putting this equations in (3.2) we deduce that the problem(3.1), (1.2)-(1.4) has an unique solution y := y f (x,λ) in the form ρ 34 χ 1 (x,λ) 1 φ 1(ξ,λ)f(ξ)dξ + ρ 34φ 1 (x,λ) x χ 1(ξ,λ)f(ξ)dξ + ρ 12φ 1 (x,λ) χ 2(ξ,λ)f(ξ)dξ for x [,) ρ 12 χ 2 (x,λ) + ρ 34χ 2 (x,λ) (3.5) φ 2(ξ,λ)f(ξ)dξ + ρ 12φ 2 (x,λ) 1 x χ 2(ξ,λ)f(ξ)dξ φ 1(ξ,λ)f(ξ)dξ for x (,π] From this formula wederive that thegreen sfunction of theproblem(3.1), (1.2)-(1.4) can be represented as { φ(ξ,λ)χ(x,λ) for ξ x π, x,ξ G(x,ξ;λ) = φ(x,λ)χ(ξ,λ) (3.6) for x ξ π, x,ξ 5

6 and the formula (3.5) can be rewritten in terms of this Green s function as y f (x,λ) = ρ 12 G(x,ξ;λ)f(ξ)dξ +ρ 34 G(x,ξ;λ)f(ξ)dξ (3.7) Theorem 3.1. The resolvent operator R(λ, A) is compact. Proof. References [1] F.V. Atkinson., SOn bounds for the Titchmarsh-Weyl m-coecients and for spectral functions for second-order differential equations,proc. Royal Soc. Edinburgh. (A) 97 (1984), 1-7. [2] N.I. Akhiezer and I.M. Glazman, Theory of linear operators in Hilbert space: I and II, Pitman, London and Scottish Academic Press, Edinburgh: [3] E.A. Coddington and N. Levinson. Theory of ordinary dierential equations, McGraw-Hill, New York: [4] I. G. Petrovsky, Lectures on Partial Differential Equations (translated from Russian by A. Shenitzer), 1st English ed., New York: Interscience Publishers, [5] O. Sh. Mukhtarov, M. Kadakal and F.S. Muhtarov On discontinuous Sturm-Liouville Problems with transmission conditions. J. Math. Kyoto Univ., Vol. 44, No. 4(24), [6] E.Kreysig, Advanced Engineering Mathematics, 8th Edition. John Wiley Sons, New York, [7] B.M. Levitan, SThe Eigenfunction Expansion for the Second Order Differential Operator, M.L [8] B.M. Levitan and I.S. Sargsjan, Sturm-Liouville and Dirac operators, Mathematics and its Applications (Soviet Series),59. Kluwer Academic Publishers Group, Dordrecht:1991. [9] M.A. Naimark, Linear dierential operators: II. A(Ungar, New York: Translated from the second Russian edition). 6

7 [1] M. Rotenberg, Theory and application of Sturmian functions. Adv. Atomic. and Mol. Phys. 6 (197). [11] A. I. Sherstyuk, Problems of Theoretical Physics (Leningrad., Gos. Univ., Leningrad, 1988), Vol. 3. [12] E. C. Titchmarsh, Eigenfunctions Expansion Associated with Second Order Differential Equations I, second edn. Oxford Univ. Press, London, (1962). 7