About this Journal Submit a Manuscript Table of Contents
International Journal of Stochastic Analysis
Volume 2010 (2010), Article ID 519684, 16 pages
doi:10.1155/2010/519684
Research Article

The Rothe's Method to a Parabolic Integrodifferential Equation with a Nonclassical Boundary Conditions

Department of Mathematics, University of Larbi Ben M'Hidi, Oum El Bouagui 04000, Algeria

Received 26 February 2009; Revised 17 August 2009; Accepted 10 December 2009

Academic Editor: Alexander M. Krasnosel'skii

Copyright © 2010 Abdelfatah Bouziani and Rachid Mechri. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

This paper is devoted to prove, in a nonclassical function space, the weak solvability of parabolic integrodifferential equations with a nonclassical boundary conditions. The investigation is made by means of approximation by the Rothes method which is based on a semidiscretization of the given problem with respect to the time variable.

1. Introduction

The purpose of this paper is to study the solvability of the following equation: 𝜕 𝑣 ( 𝜕 𝜕 𝑡 𝑥 , 𝑡 ) 2 𝑣 𝜕 𝑥 2 ( 𝑥 , 𝑡 ) = 𝑡 0 𝑎 ( 𝑡 𝑠 ) 𝑘 ( [ ] 𝑠 , 𝑣 ( 𝑥 , 𝑠 ) ) 𝑑 𝑠 + 𝑔 ( 𝑥 , 𝑡 ) , ( 𝑥 , 𝑡 ) ( 0 , 1 ) × 0 , 𝑇 , ( 1 . 1 ) with the initial condition 𝑣 ( 𝑥 , 0 ) = 𝑉 0 ( 𝑥 ) , 𝑥 ( 0 , 1 ) , ( 1 . 2 ) and the integral conditions 1 0 [ ] , 𝑣 ( 𝑥 , 𝑡 ) 𝑑 𝑥 = 𝐸 ( 𝑡 ) , 𝑡 0 , 𝑇 1 0 [ ] , 𝑥 𝑣 ( 𝑥 , 𝑡 ) 𝑑 𝑥 = 𝐺 ( 𝑡 ) , 𝑡 0 , 𝑇 ( 1 . 3 ) where 𝑣 is an unknown function, 𝐸 , 𝐺 , and 𝑉 0 are given functions supposed to be sufficiently regular, while 𝑘 and 𝑎 are suitably defined functions satisfying certain conditions to be specified later and 𝑇 is a positive constant.

Since 1930, various classical types of initial boundary value problems have been investigated by many authors using Rothe time-discretization method; see, for instance, the monographs by Rektorys [1] and Kačur [2] and references cited therein. The linear case of our problem, that is, 𝑡 0 𝑎 ( 𝑡 𝑠 ) 𝑘 ( 𝑠 , 𝑣 ( 𝑥 , 𝑠 ) ) 𝑑 𝑠 = 0 , appears, for instance, in the modelling of the quasistatic flexure of a thermoelastic rod (see [3]) and has been studied, firstly, by the second author with a more general second-order parabolic equation or a 2m-parabolic equation in [35] by means of the energy-integrals method and, secondly, by the two authors via the Rothe method [68]. For other models, we refer the reader, for instance, to [912], and references therein.

The paper is organized as follows. In Section 2, we transform problem (1.1)–(1.3) to an equivalent one with homogeneous integral conditions, namely, problem (2.3). Then, we specify notations and assumptions on data before stating the precise sense of the desired solution. In Section 3, by the Rothe discretization in time method, we construct approximate solutions to problem (2.3). Some a priori estimates for the approximations are derived in Section 4, while Section 5 is devoted to establish the existence and uniqueness of the solution.

2. Preliminaries, Notation, and Main Result

It is convenient at the beginning to reduce problem (1.1)–(1.3) with inhomogeneous integral conditions to an equivalent one with homogeneous conditions. For this, we introduce a new unknown function 𝑢 by setting

[ ] 𝑢 ( 𝑥 , 𝑡 ) = 𝑣 ( 𝑥 , 𝑡 ) 𝑅 ( 𝑥 , 𝑡 ) , ( 𝑥 , 𝑡 ) ( 0 , 1 ) × 0 , 𝑇 , ( 2 . 1 ) where

𝑅 ( 𝑥 , 𝑡 ) = 6 ( 2 𝐺 ( 𝑡 ) 𝐸 ( 𝑡 ) ) 𝑥 2 ( 3 𝐺 ( 𝑡 ) 2 𝐸 ( 𝑡 ) ) . ( 2 . 2 ) Then, the function 𝑢 is seen to be the solution of the following problem:

𝜕 𝑢 ( 𝜕 𝜕 𝑡 𝑥 , 𝑡 ) 2 𝑢 𝜕 𝑥 2 ( 𝑥 , 𝑡 ) = 𝑡 0 [ ] , 𝑎 ( 𝑡 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑥 , 𝑠 ) ) 𝑑 𝑠 + 𝑓 ( 𝑥 , 𝑡 ) , ( 𝑥 , 𝑡 ) ( 0 , 1 ) × 0 , 𝑇 𝑢 ( 𝑥 , 0 ) = 𝑈 0 ( 𝑥 ) , 𝑥 ( 0 , 1 ) , 1 0 [ ] , 𝑢 ( 𝑥 , 𝑡 ) 𝑑 𝑥 = 0 , 𝑡 0 , 𝑇 1 0 [ ] , 𝑥 𝑢 ( 𝑥 , 𝑡 ) 𝑑 𝑥 = 0 , 𝑡 0 , 𝑇 ( 2 . 3 ) where 𝑓 ( 𝑥 , 𝑡 ) = 𝑔 ( 𝑥 , 𝑡 ) 𝜕 𝑅 ( 𝑥 , 𝑡 ) , 𝑈 𝜕 𝑡 0 ( 𝑥 ) = 𝑉 0 ( 𝑥 ) 𝑅 ( 𝑥 , 0 ) , 𝑘 ( 𝑠 , 𝑢 ( 𝑥 , 𝑠 ) ) = 𝑘 ( 𝑠 , 𝑢 ( 𝑥 , 𝑠 ) ) 𝑅 ( 𝑥 , 𝑡 ) . ( 2 . 4 ) Hence, instead of looking for the function 𝑣 , we search for the function 𝑢 . The solution of problem (1.1)–(1.3) will be simply given by the formula 𝑣 ( 𝑥 , 𝑡 ) = 𝑢 ( 𝑥 , 𝑡 ) + 𝑅 ( 𝑥 , 𝑡 ) .

We introduce the function spaces, which we need in our investigation. Let 𝐿 2 ( 0 , 1 ) and 𝐿 2 ( 0 , 𝑇 ; 𝐿 2 ( 0 , 1 ) ) be the standard function spaces. We denote by 𝐶 0 ( 0 , 1 ) the linear space of continuous functions with compact support in ( 0 , 1 ) . Since such functions are Lebesgue integrable, we can define on 𝐶 0 ( 0 , 1 ) the bilinear form given by ( ( 𝑢 , 𝑣 ) ) = 1 0 𝑥 𝑢 𝑥 𝑣 𝑑 𝑥 , ( 2 . 5 ) where

𝑥 𝑢 = 𝑥 0 𝑢 ( 𝜁 , ) 𝑑 𝜁 . ( 2 . 6 ) The bilinear form (2.5) is considered as a scalar product on 𝐶 0 ( 0 , 1 ) for which 𝐶 0 ( 0 , 1 ) is not complete.

Definition 2.1. We denote by 𝐵 1 2 ( 0 , 1 ) a completion of 𝐶 0 ( 0 , 1 ) for the scalar product (2.5), which is denoted by ( , ) 𝐵 1 2 ( 0 , 1 ) , called the Bouziani space or the space of square integrable primitive functions on ( 0 , 1 ) . By the norm of function 𝑢 from 𝐵 1 2 ( 0 , 1 ) , we understand the nonnegative number 𝑢 𝐵 1 2 ( 0 , 1 ) = ( 𝑢 , 𝑢 ) 𝐵 1 2 ( 0 , 1 ) = 𝑥 𝑢 , ( 2 . 7 ) where 𝑣 denotes the norm of 𝑣 in 𝐿 2 ( 0 , 1 ) .

For 𝑢 𝐿 2 ( 0 , 1 ) , we have the elementary inequality 𝑢 𝐵 1 2 ( 0 , 1 ) 1 2 𝑢 . ( 2 . 8 )

We denote by 𝐿 2 ( 0 , 𝑇 ; 𝐵 1 2 ( 0 , 1 ) ) the space of functions which are square integrable in the Bochner sense, with the scalar product

( 𝑢 , 𝑣 ) 𝐿 2 ( 0 , 𝑇 ; 𝐵 1 2 ( 0 , 1 ) ) = 𝑇 0 ( 𝑢 ( , 𝑡 ) , 𝑣 ( , 𝑡 ) ) 𝐵 1 2 ( 0 , 1 ) 𝑑 𝑡 . ( 2 . 9 )

Since the space 𝐵 1 2 ( 0 , 1 ) is a Hilbert space, it can be shown that 𝐿 2 ( 0 , 𝑇 ; 𝐵 1 2 ( 0 , 1 ) ) is a Hilbert space as well. The set of all continuous abstract functions in [ 0 , 𝑇 ] equipped with the norm

s u p 0 𝜏 𝑇 𝑢 ( , 𝜏 ) 𝐵 1 2 ( 0 , 1 ) ( 2 . 1 0 ) is denoted 𝐶 ( 0 , 𝑇 ; 𝐵 1 2 ( 0 , 1 ) ) . Let 𝑉 be the set which we define as follows:

𝑉 = 𝑣 𝐿 2 ( 0 , 1 ) ; 1 0 𝑣 ( 𝑥 ) 𝑑 𝑥 = 1 0 𝑥 𝑣 ( 𝑥 ) 𝑑 𝑥 = 0 . ( 2 . 1 1 )

Since 𝑉 is the null space of the continuous linear mapping 𝑙 : 𝐿 2 ( 0 , 1 ) 2 , 𝜑 𝑙 ( 𝜑 ) = ( 1 0 𝜑 ( 𝑥 ) 𝑑 𝑥 , 1 0 𝑥 𝜑 ( 𝑥 ) 𝑑 𝑥 ) , it is a closed linear subspace of 𝐿 2 ( 0 , 1 ) , consequently 𝑉 is a Hilbert space endowed with the inner product ( , ) . Strong or weak convergence is denoted by or , respectively. The letter 𝐶 will stand for a generic positive constant which may be different in the same discussion.

Lemma 2.2 (Gronwall's lemma). ( a 1 ) Let 𝑥 ( 𝑡 ) 0 , ( 𝑡 ) , 𝑦 ( 𝑡 ) be real integrable functions on the interval [ 𝑎 , 𝑏 ] . If 𝑦 ( 𝑡 ) ( 𝑡 ) + 𝑡 𝑎 𝑥 ( 𝑠 ) 𝑦 ( 𝑠 ) 𝑑 𝑠 , 𝑡 ( 𝑎 , 𝑏 ) , ( 2 . 1 2 ) then 𝑦 ( 𝑡 ) ( 𝑡 ) + 𝑡 𝑎 ( 𝑠 ) 𝑥 ( 𝑠 ) e x p 𝑡 𝑎 𝑥 ( 𝜏 ) 𝑑 𝜏 𝑑 𝑠 , 𝑡 ( 0 , 𝑇 ) . ( 2 . 1 3 )
In particular, if 𝑥 ( 𝑡 ) 𝐶 is a constant and ( 𝑡 ) is nondecreasing, then 𝑦 ( 𝑡 ) ( 𝑡 ) 𝑒 𝑐 ( 𝑡 𝑎 ) , 𝑡 ( 0 , 𝑇 ) . ( 2 . 1 4 )
( a 2 ) Let { 𝑎 𝑖 } 𝑖 be a sequence of real nonnegative numbers satisfying 𝑎 𝑖 𝐴 + 𝐵 𝑖 1 𝑘 = 1 𝑎 𝑘 , 𝑖 = 1 , 2 , , ( 2 . 1 5 ) where 𝐴 , 𝐵 , and are positive constants, such that 𝐵 < 1 . Then 𝑎 𝑖 [ ] 𝐴 e x p 𝐵 ( 𝑖 1 ) , ( 2 . 1 6 ) takes place for all 𝑖 = 1 , 2 , .

In the sequel, we make the following assumptions.

( 𝐻 1 ) Functions 𝑓 [ 0 , 𝑇 ] 𝐿 2 ( 0 , 1 ) and 𝑎 [ 0 , 𝑇 ] are Lipschitz continuous, that is,

𝑙 1 + ; 𝑓 𝑡 ( 𝑡 ) 𝑓 𝑙 1 | | 𝑡 𝑡 | | [ ] , , 𝑡 0 , 𝑇 𝑙 2 + ; | | 𝑡 𝑎 ( 𝑡 ) 𝑎 | | 𝑙 2 | | 𝑡 𝑡 | | [ ] . , 𝑡 0 , 𝑇 ( 2 . 1 7 )

( 𝐻 2 ) The mapping 𝑘 [ 0 , 𝑇 ] × 𝑉 𝐿 2 ( 0 , 1 ) is Lipschitz continuous in both variables, that is,

𝑙 3 + ; 𝑘 𝑡 ( 𝑡 , 𝑢 ) 𝑘 , 𝑢 𝑙 3 | | 𝑡 𝑡 | | + 𝑢 𝑢 , ( 2 . 1 8 ) for all 𝑡 , 𝑡 𝐼 , 𝑢 , 𝑢 𝑉 , and satisfies

𝑙 4 , 𝑙 5 + ; 𝑘 ( 𝑡 , 𝑢 ) 𝐵 1 2 ( 0 , 1 ) 𝑙 4 𝑢 𝐵 1 2 ( 0 , 1 ) + 𝑙 5 , ( 2 . 1 9 ) for all 𝑡 𝐼 and all 𝑢 𝑉 , where 𝑙 4 and 𝑙 5 are positive constants.

( 𝐻 3 ) 𝑈 0 𝐻 2 ( 0 , 1 ) and

1 0 𝑈 0 ( 𝑥 ) 𝑑 𝑥 = 1 0 𝑥 𝑈 0 ( 𝑥 ) 𝑑 𝑥 = 0 . ( 2 . 2 0 )

We will be concerned with a weak solution in the following sense.

Definition 2.3. A function 𝑢 𝐼 𝐿 2 ( 0 , 1 ) is called a weak solution to problem (2.3) if the following conditions are satisfied:(i) 𝑢 𝐿 ( 𝐼 , 𝑉 ) 𝐶 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) ,(ii) 𝑢 is strongly differentiable a.e. in 𝐼 and 𝑑 𝑢 / 𝑑 𝑡 𝐿 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) ,(iii) 𝑢 ( 0 ) = 𝑈 0 in 𝑉 ,(iv)the identity 𝑑 𝑢 𝑑 𝑡 ( 𝑡 ) , 𝑣 𝐵 1 2 ( 0 , 1 ) = + ( 𝑢 ( 𝑡 ) , 𝑣 ) 𝑡 0 𝑎 ( 𝑡 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝑑 𝑠 , 𝑣 𝐵 1 2 ( 0 , 1 ) + ( 𝑓 ( 𝑡 ) , 𝑣 ) 𝐵 1 2 ( 0 , 1 ) ( 2 . 2 1 ) holds for all 𝑣 𝑉 and a.e. 𝑡 [ 0 , 𝑇 ] .

To close this section, we announce the main result of the paper.

Theorem 2.4. Under assumptions ( 𝐻 1 ) - ( 𝐻 3 ) , problem (2.3) admits a unique weak solution 𝑢 , in the sense of Definition (2.3).

3. Construction of an Approximate Solution

In order to solve problem (2.3) by the Rothe method, we proceed as follows. Let 𝑛 be a positive integer, we divide the time interval 𝐼 = [ 0 , 𝑇 ] into 𝑛 subintervals 𝐼 𝑛 𝑗 = [ 𝑡 𝑛 𝑗 1 , 𝑡 𝑛 𝑗 ] , 𝑗 = 1 , , 𝑛 , where 𝑡 𝑛 𝑗 = 𝑗 𝑛 and 𝑛 = 𝑇 / 𝑛 . Then, for each 𝑛 1 , problem (2.3) may be approximated by the following recurrent sequence of time-discretized problems. Successively, for 𝑗 = 1 , , 𝑛 , we look for functions 𝑢 𝑛 𝑗 𝑉 such that 𝑢 𝑛 𝑗 𝑢 𝑛 𝑗 1 𝑛 𝑑 2 𝑢 𝑛 𝑗 𝑑 𝑥 2 = 𝑛 𝑗 1 𝑖 = 0 𝑎 𝑡 𝑛 𝑗 𝑡 𝑛 𝑖 𝑘 𝑡 𝑛 𝑖 , 𝑢 𝑛 𝑖 + 𝑓 𝑛 𝑗 , ( 3 . 1 ) 1 0 𝑢 𝑛 𝑗 ( 𝑥 ) 𝑑 𝑥 = 0 , ( 3 . 2 ) 1 0 𝑥 𝑢 𝑛 𝑗 ( 𝑥 ) 𝑑 𝑥 = 0 , ( 3 . 3 ) starting from

𝑢 𝑛 0 = 𝑈 0 , 𝛿 𝑢 𝑛 0 = 𝑑 2 𝑑 𝑥 2 𝑈 0 + 𝑓 ( 0 ) , ( 3 . 4 ) where 𝑢 𝑛 𝑗 ( 𝑥 ) = 𝑢 ( 𝑥 , 𝑡 𝑛 𝑗 ) , 𝛿 𝑢 𝑛 𝑗 = ( 𝑢 𝑛 𝑗 𝑢 𝑛 𝑗 1 ) / 𝑛 , 𝑓 𝑛 𝑗 ( 𝑥 ) = 𝑓 ( 𝑥 , 𝑡 𝑛 𝑗 ) . For this, multiplying for all 𝑗 = 1 , , 𝑛 , (3.1) by 2 𝑥 𝑣 = 𝑥 0 ( 𝜉 0 𝑣 ( 𝜏 ) 𝑑 𝜏 ) 𝑑 𝜉 and integrating over ( 0 , 1 ) , we get

1 0 𝛿 𝑢 𝑛 𝑗 ( 𝑥 ) 2 𝑥 𝑣 𝑑 𝑥 1 0 𝑑 2 𝑢 𝑛 𝑗 𝑑 𝑥 2 ( 𝑥 ) 2 𝑥 𝑣 𝑑 𝑥 = 𝑛 1 0 𝑗 1 𝑖 = 0 𝑎 𝑡 𝑛 𝑗 𝑡 𝑛 𝑖 𝑘 𝑡 𝑛 𝑖 , 𝑢 𝑛 𝑖 2 𝑥 𝑣 𝑑 𝑥 + 1 0 𝑓 𝑛 𝑗 2 𝑥 𝑣 𝑑 𝑥 . ( 3 . 5 ) Noting that, using a standard integration by parts, we have

2 1 𝑣 = 1 0 ( 1 𝜉 ) 𝑣 ( 𝜉 ) 𝑑 𝜉 = 1 0 𝑣 ( 𝜉 ) 𝑑 𝜉 1 0 𝜉 𝑣 ( 𝜉 ) 𝑑 𝜉 = 0 , 𝑣 𝑉 . ( 3 . 6 ) Carrying out some integrations by parts and invoking (3.6), we obtain for each term in (3.5) 1 0 𝛿 𝑢 𝑛 𝑗 2 𝑥 𝑣 𝑑 𝑥 = 𝛿 𝑢 𝑛 𝑗 , 𝑣 𝐵 1 2 ( 0 , 1 ) , 1 0 𝑑 2 𝑢 𝑛 𝑗 𝑑 𝑥 2 ( 𝑥 ) 2 𝑥 𝑢 𝑣 𝑑 𝑥 = 𝑛 𝑗 , , 𝑣 𝑛 1 0 𝑗 1 𝑖 = 0 𝑎 𝑡 𝑛 𝑗 𝑡 𝑛 𝑖 𝑘 𝑡 𝑛 𝑖 , 𝑢 𝑛 𝑖 ( 𝑥 ) 2 𝑥 𝑣 𝑑 𝑥 = 𝑛 𝑗 1 𝑖 = 0 𝑎 𝑡 𝑛 𝑗 𝑡 𝑛 𝑖 𝑘 𝑡 𝑛 𝑖 , 𝑢 𝑛 𝑖 , 𝑣 𝐵 1 2 ( 0 , 1 ) , ( 3 . 7 ) and for the last one

1 0 𝑓 𝑛 𝑗 ( 𝑥 ) 2 𝑥 𝑓 𝑣 ( 𝑥 ) 𝑑 𝑥 = 𝑛 𝑗 , 𝑣 𝐵 1 2 ( 0 , 1 ) . ( 3 . 8 ) By virtue of (3.7) and (3.8), (3.5) becomes 𝛿 𝑢 𝑛 𝑗 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑛 𝑗 , 𝑣 = 𝑛 𝑗 1 𝑖 = 0 𝑎 𝑡 𝑛 𝑗 𝑡 𝑛 𝑖 𝑘 𝑡 𝑛 𝑖 , 𝑢 𝑛 𝑖 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑓 𝑛 𝑗 , 𝑣 𝐵 1 2 ( 0 , 1 ) , ( 3 . 9 ) or

𝑢 𝑛 𝑗 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑛 𝑢 𝑛 𝑗 , 𝑣 = 2 𝑛 𝑗 1 𝑖 = 0 𝑎 𝑡 𝑛 𝑗 𝑡 𝑛 𝑖 𝑘 𝑡 𝑛 𝑖 , 𝑢 𝑛 𝑖 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑛 𝑓 𝑛 𝑗 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑛 𝑗 1 , 𝑣 𝐵 1 2 ( 0 , 1 ) . ( 3 . 1 0 ) Let 𝜂 ( , ) 𝑉 × 𝑉 and 𝐿 𝑗 ( ) 𝑉 be two functions defined by

𝜂 ( 𝑢 , 𝑣 ) = ( 𝑢 , 𝑣 ) 𝐵 1 2 ( 0 , 1 ) + 𝑛 𝐿 ( 𝑢 , 𝑣 ) , 𝑗 ( 𝑣 ) = 2 𝑛 𝑗 1 𝑖 = 0 𝑎 𝑡 𝑛 𝑗 𝑡 𝑛 𝑖 𝑘 𝑡 𝑛 𝑖 , 𝑢 𝑛 𝑖 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑛 𝑓 𝑛 𝑗 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑛 𝑗 1 , 𝑣 𝐵 1 2 ( 0 , 1 ) . ( 3 . 1 1 )

It is easy to see that the bilinear form 𝜂 ( , ) is continuous on 𝑉 and V-elliptic, and the form 𝐿 𝑗 ( ) is continuous for each 𝑗 = 1 , , 𝑛 . Then, Lax-Milgram lemma guarantees the existence and uniqueness of 𝑢 𝑛 𝑗 , for all 𝑗 = 1 , , 𝑛 .

4. A Priori Estimates

Lemma 4.1. There exists 𝐶 > 0 such that, for all 𝑛 1 and all 𝑗 = 1 , , 𝑛 , the solution 𝑢 𝑗 of the discretized problem (3.1)–(3.4) satisfies the estimates 𝑢 𝑛 𝑗 𝐶 , ( 4 . 1 ) 𝛿 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) 𝐶 . ( 4 . 2 )

Proof. Testing the difference (3.9) 𝑗 1 -(3.9) 𝑗 with 𝑣 = 𝛿 𝑢 𝑛 𝑗 ( 𝑉 ) , taking into account assumptions ( 𝐻 1 ) - ( 𝐻 3 ) and the Cauchy-Schwarz inequality, we obtain 𝛿 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑛 𝑗 𝑢 𝑛 𝑗 1 𝐵 1 2 ( 0 , 1 ) 𝛿 𝑢 𝑛 𝑗 1 𝐵 1 2 ( 0 , 1 ) + 𝐶 1 3 2 𝑛 𝑗 2 𝑖 = 0 𝑢 𝑛 𝑖 𝐵 1 2 ( 0 , 1 ) + 𝐶 1 3 𝑛 + 𝐶 1 3 𝑛 𝑢 𝑛 𝑗 1 𝐵 1 2 ( 0 , 1 ) , ( 4 . 3 ) where 𝐶 1 𝑙 = 3 m a x 2 𝜁 , 𝑇 𝑙 2 𝜁 + 𝑀 1 𝜁 + 𝑙 1 , 𝑀 1 = m a x 𝑡 𝐼 | | 𝑎 | | 𝑙 ( 𝑡 ) , 𝜁 = m a x 4 , 𝑙 5 . ( 4 . 4 ) Multiplying the left-hand side of the last inequality with ( 1 ( 𝐶 1 / 3 ) 𝑛 ) ( < 1 ) and adding the terme 2 3 𝐶 1 𝑛 𝑢 𝑛 𝑗 𝑢 𝑛 𝑗 1 𝐵 1 2 ( 0 , 1 ) 𝛿 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) ( < 0 ) , ( 4 . 5 ) we get 1 𝐶 1 𝑛 𝛿 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) 𝑢 𝑛 𝑗 1 𝐵 1 2 ( 0 , 1 ) + 𝛿 𝑢 𝑛 𝑗 1 𝐵 1 2 ( 0 , 1 ) + 𝐶 1 2 𝑛 𝑗 2 𝑖 = 0 𝑢 𝑛 𝑖 𝐵 1 2 ( 0 , 1 ) + 𝐶 1 𝑛 . ( 4 . 6 ) Applying the last inequality recursively, it follows that 1 𝐶 1 𝑛 𝑗 𝛿 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) 𝑢 𝑛 0 𝐵 1 2 ( 0 , 1 ) + 𝛿 𝑢 𝑛 0 𝐵 1 2 ( 0 , 1 ) + 𝐶 1 𝑇 + 𝑇 𝐶 1 𝑛 𝑗 2 𝑖 = 0 𝑢 𝑛 𝑖 𝐵 1 2 ( 0 , 1 ) , ( 4 . 7 ) or, by virtue of Lemma 2.2 , there exists 𝑛 0 such that 𝛿 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) 𝐶 2 , 𝑛 𝑛 0 , ( 4 . 8 ) where 𝐶 2 = e x p 𝑇 𝐶 1 + 1 𝛿 𝑢 𝑛 0 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑛 0 𝐵 1 2 ( 0 , 1 ) + 𝑇 𝐶 1 × e x p e x p 𝑇 𝐶 1 + 1 𝑇 𝐶 1 , ( 4 . 9 ) and so our proof is complete.

We address now the question of convergence and existence.

5. Convergence and Existence

Now let us introduce the Rothe function 𝑢 𝑛 ( 𝑡 ) 𝐼 𝑉 obtained from the functions 𝑢 𝑗 by piecewise linear interpolation with respect to time

𝑢 𝑛 ( 𝑡 ) = 𝑢 𝑛 𝑗 1 + 𝛿 𝑢 𝑛 𝑗 𝑡 𝑡 𝑛 𝑗 1 , i n 𝐼 𝑛 𝑗 , ( 5 . 1 ) as well the step functions ̃ 𝑢 𝑛 ( 𝑡 ) , ̂ 𝑢 𝑛 ( 𝑡 ) , 𝑓 𝑛 ( 𝑡 ) , and ̃ 𝑘 ( 𝑡 , ̃ 𝑢 𝑛 ( 𝑡 ) ) defined as follows:

̃ 𝑢 𝑛 𝑢 ( 𝑡 ) = 𝑛 0 𝑢 , f o r 𝑡 = 0 , 𝑛 𝑗 𝐼 , i n 𝑛 𝑗 𝑡 = 𝑛 𝑗 1 , 𝑡 𝑛 𝑗 , ̂ 𝑢 𝑛 𝑢 ( 𝑡 ) = 𝑛 0 𝑢 , f o r 𝑡 = 0 , 𝑛 𝑗 1 𝐼 , i n 𝑛 𝑗 , 𝑓 ( 5 . 2 ) 𝑛 𝑓 ( 𝑡 ) = 𝑓 ( 0 ) , f o r 𝑡 = 0 , 𝑛 𝑗 𝐼 , i n 𝑛 𝑗 , ̃ 𝑘 ( 5 . 3 ) 𝑛 ( 𝑡 ) = 0 , f o r 𝑡 = 0 , 𝑛 𝑗 1 𝑖 = 0 𝑎 𝑡 𝑛 𝑗 𝑡 𝑛 𝑖 𝑘 𝑡 𝑛 𝑖 , 𝑢 𝑛 𝑖 𝐼 , i n 𝑛 𝑗 = 𝑡 𝑛 𝑗 1 , 𝑡 𝑛 𝑗 . ( 5 . 4 )

Corollary 5.1. There exist 𝐶 > 0 such that the estimates 𝑢 𝑛 ( 𝑡 ) 𝐶 , ̃ 𝑢 𝑛 ( 𝑡 ) 𝐶 , 𝑡 𝐼 , ( 5 . 5 ) 𝑑 𝑢 𝑛 𝑑 𝑡 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝐶 , f o r a . e . 𝑡 𝐼 , ( 5 . 6 ) ̃ 𝑢 𝑛 ( 𝑡 ) 𝑢 𝑛 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 , ̂ 𝑢 𝑛 ( 𝑡 ) 𝑢 𝑛 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 ̃ 𝑘 , 𝑡 𝐼 , ( 5 . 7 ) 𝑛 ( 𝑡 ) 𝐶 , 𝑡 𝐼 , ( 5 . 8 ) hold for all 𝑛 .

Proof. For the inequalities (5.5), (5.6), and (5.7) see [6, Corollary  4.2.], whereas for the last inequality, assumption ( 𝐻 2 ) and estimate (4.1) guarantee the desired result.

Proposition 5.2. The sequence ( 𝑢 𝑛 ) 𝑛 converges in the norm of the space 𝐶 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) to some function 𝑢 𝐶 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) and the error estimate 𝑢 𝑛 𝑢 𝐶 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) 𝐶 𝑛 ( 5 . 9 ) takes place for all 𝑛 𝑛 0 .

Proof. By virtue of (5.2), (5.3), and (5.4) the variational equation (3.9) may be rewritten in the form 𝑑 𝑢 𝑛 𝑑 𝑡 ( 𝑡 ) , 𝑣 𝐵 1 2 ( 0 , 1 ) + ̃ 𝑢 𝑛 = ̃ 𝑘 ( 𝑡 ) , 𝑣 𝑛 ( 𝑡 ) , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑓 𝑛 ( 𝑡 ) , 𝑣 𝐵 1 2 ( 0 , 1 ) , ( 5 . 1 0 ) for a.e. 𝑡 [ 0 , 𝑇 ] . In view of (5.10), using (5.6) and (5.8) with the fact that 𝑓 𝑛 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝑀 2 = m a x 𝑡 𝐼 ( 𝑓 𝑡 ) 𝐵 1 2 ( 0 , 1 ) < , ( 5 . 1 1 ) we obtain | | ̃ 𝑢 𝑛 | | ̃ 𝑘 ( 𝑡 ) , 𝑣 𝑛 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) + 𝑓 𝑛 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) + 𝑑 𝑢 𝑛 𝑑 𝑡 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝑣 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑣 𝐵 1 2 ( 0 , 1 ) [ ] . , a . e . 𝑡 0 , 𝑇 ( 5 . 1 2 ) Now, for 𝑛 , 𝑚 being two positive integers, testing the difference (5.10) 𝑛 -(5.10) 𝑚 with 𝑣 = 𝑢 𝑛 ( 𝑡 ) 𝑢 𝑚 ( 𝑡 ) which is in 𝑉 , with the help of the Cauchy-Schwarz inequality and taking into account that 2 𝑑 𝑑 𝑡 𝑢 ( 𝑡 ) , 𝑢 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) = 𝑑 𝑑 𝑡 𝑢 ( 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) [ ] , a . e . 𝑡 0 , 𝑇 , ( 5 . 1 3 ) and, by virtue of (5.12) we obtain after some rearrangements 1 2 𝑑 𝑑 𝑡 𝑢 𝑛 ( 𝑡 ) 𝑢 𝑚 ( 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) + ̃ 𝑢 𝑛 ( 𝑡 ) ̃ 𝑢 𝑚 ( 𝑡 ) 2 𝑢 𝐶 𝑚 ( 𝑡 ) ̃ 𝑢 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) + 𝐶 ̃ 𝑢 𝑛 ( 𝑡 ) 𝑢 𝑛 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) + ̃ 𝑘 𝑛 ̃ 𝑘 ( 𝑡 ) 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝑢 𝑛 ( 𝑡 ) 𝑢 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) + 𝑓 𝑛 𝑓 ( 𝑡 ) 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝑢 𝑛 ( 𝑡 ) 𝑢 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) [ ] . , a . e . 𝑡 0 , 𝑇 ( 5 . 1 4 ) To derive the required result, we need to estimate the third and the last terms in the right-hand side, for this, let 𝑡 be arbitrary but fixed in ( 0 , 𝑇 ] , without loss of generality we can suppose that there exist three positive integers 𝑝 , 𝑞 and 𝛽 , such that 𝑡 𝑡 𝑛 𝑝 1 , 𝑡 𝑛 𝑝 𝑡 𝑚 𝑞 1 , 𝑡 𝑚 𝑞 , 𝑛 = 𝛽 𝑚 , 𝑡 𝑛 𝑝 = 𝑡 𝑚 𝑞 . ( 5 . 1 5 ) Hence, using (5.4) we can write ̃ 𝑘 𝑛 ̃ 𝑘 ( 𝑡 ) 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) = 𝑚 𝑝 1 𝑗 = 0 𝛽 ( 𝑗 + 1 ) 1 𝑖 = 𝑗 𝛽 𝑎 𝑡 𝑛 𝑝 𝑡 𝑛 𝑗 𝑘 𝑡 𝑛 𝑗 , 𝑢 𝑛 𝑗 𝑡 𝑎 𝑚 𝑞 𝑡 𝑚 𝑖 𝑘 𝑡 𝑚 𝑖 , 𝑢 𝑚 𝑖 𝐵 1 2 ( 0 , 1 ) . ( 5 . 1 6 ) By virtue of assumption ( 𝐻 1 ) and the fact that | 𝑎 ( 𝑡 𝑛 𝑝 𝑡 𝑛 𝑗 ) 𝑎 ( 𝑡 𝑚 𝑞 𝑡 𝑚 𝑖 ) | 𝐶 𝑛 , there exist 𝜀 𝑛 [ 0 , 𝐶 𝑛 ] such that ̃ 𝑘 𝑛 ̃ 𝑘 ( 𝑡 ) 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝑚 𝑝 1 𝑗 = 0 𝛽 ( 𝑗 + 1 ) 1 𝑖 = 𝑗 𝛽 𝐶 𝑛 𝜀 𝑛 𝑘 𝑡 𝑛 𝑗 , 𝑢 𝑛 𝑗 𝐵 1 2 ( 0 , 1 ) + | | 𝑎 𝑡 𝑚 𝑞 𝑡 𝑚 𝑖 | | 𝑘 𝑡 𝑛 𝑗 , 𝑢 𝑛 𝑗 𝑡 𝑘 𝑚 𝑖 , 𝑢 𝑚 𝑖 𝐵 1 2 ( 0 , 1 ) . ( 5 . 1 7 ) Therefore, recalling assumptions ( 𝐻 1 ) , ( 𝐻 2 ) and having in mind that 𝜀 𝑛 [ 0 , 𝐶 𝑛 ] , we estimate ̃ 𝑘 𝑛 ̃ 𝑘 ( 𝑡 ) 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝑚 𝑝 1 𝑗 = 0 𝛽 ( 𝑗 + 1 ) 1 𝑖 = 𝑗 𝛽 𝐶 𝑛 + 𝐶 𝑛 + 𝑢 𝑛 𝑗 𝑢 𝑚 𝑖 𝐵 1 2 ( 0 , 1 ) , ( 5 . 1 8 ) from where, we derive for all 𝑠 ( 𝑡 𝑚 𝑖 , 𝑡 𝑚 𝑖 + 1 ] ̃ 𝑘 𝑛 ̃ 𝑘 ( 𝑡 ) 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝑚 𝑝 1 𝑗 = 0 𝛽 ( 𝑗 + 1 ) 1 𝑖 = 𝑗 𝛽 𝐶 𝑛 + 𝐶 𝑛 + ̃ 𝑢 𝑛 ( 𝑠 ) 𝑢 𝑛 ( 𝑠 ) 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑛 ( 𝑠 ) 𝑢 𝑚 ( 𝑠 ) 𝐵 1 2 ( 0 , 1 ) + 𝑢 𝑚 ( 𝑠 ) ̃ 𝑢 𝑚 ( 𝑠 ) 𝐵 1 2 ( 0 , 1 ) . ( 5 . 1 9 ) Taking the supremum with respect to 𝑠 from 0 to 𝑡 in the right-hand side, invoking the fact that 𝑠 ( 𝑡 𝑚 𝑖 , 𝑡 𝑚 𝑖 + 1 ] ( 𝑡 𝑛 𝑗 1 , 𝑡 𝑛 𝑗 ] and estimate (5.7), we obtain ̃ 𝑘 𝑛 ̃ 𝑘 ( 𝑡 ) 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝑚 𝑞 1 𝑖 = 0 𝐶 𝑛 + 𝐶 s u p 0 𝑠 𝑡 𝑢 𝑛 ( 𝑠 ) 𝑢 𝑚 ( 𝑠 ) 𝐵 1 2 ( 0 , 1 ) , ( 5 . 2 0 ) so that ̃ 𝑘 𝑛 ̃ 𝑘 ( 𝑡 ) 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 + 𝐶 s u p 0 𝑠 𝑡 𝑢 𝑛 ( 𝑠 ) 𝑢 𝑚 ( 𝑠 ) 𝐵 1 2 ( 0 , 1 ) . ( 5 . 2 1 ) Let 𝑡 ( 𝑡 𝑛 𝑝 1 , 𝑡 𝑛 𝑝 ] ( 𝑡 𝑚 𝑞 1 , 𝑡 𝑚 𝑞 ] , from assumption ( 𝐻 1 ) it follows that 𝑓 𝑛 ( 𝑓 𝑡 ) 𝑚 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) = 𝑓 𝑡 𝑛 𝑝 𝑡 𝑓 𝑚 𝑞 𝐵 1 2 ( 0 , 1 ) 𝑙 1 | | 𝑡 𝑛 𝑝 𝑡 𝑚 𝑞 | | 𝑙 1 𝑛 . ( 5 . 2 2 ) Ignoring the second term in the left-hand side of (5.14) which is clearly positive and using estimates (5.5), (5.7), (5.21), and (5.22) yield 𝑑 𝑑 𝑡 𝑢 𝑛 ( 𝑡 ) 𝑢 𝑚 ( 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 + 𝑚 + 𝐶 s u p 0 𝑠 𝑡 𝑢 𝑛 ( 𝑠 ) 𝑢 𝑚 ( 𝑠 ) 2 𝐵 1 2 ( 0 , 1 ) [ ] , a . e . 𝑡 0 , 𝑇 . ( 5 . 2 3 ) Integrating this inequality with respect to time from 0 to 𝑡 and invoking the fact that 𝑢 𝑛 ( 0 ) = 𝑢 𝑚 ( 0 ) = 𝑈 0 , we get 𝑢 𝑛 ( 𝑡 ) 𝑢 𝑚 ( 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 + 𝑚 + 𝐶 𝑡 0 s u p 0 𝜉 𝑡 𝑢 𝑛 ( 𝜉 ) 𝑢 𝑚 ( 𝜉 ) 2 𝐵 1 2 ( 0 , 1 ) 𝑑 𝜉 , ( 5 . 2 4 ) whence s u p 0 𝑠 𝑡 𝑢 𝑛 ( 𝑠 ) 𝑢 𝑚 ( 𝑠 ) 2 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 + 𝑚 + 𝐶 𝑡 0 s u p 0 𝜉 𝑡 𝑢 𝑛 ( 𝜉 ) 𝑢 𝑚 ( 𝜉 ) 2 𝐵 1 2 ( 0 , 1 ) 𝑑 𝜉 . ( 5 . 2 5 ) Accordingly, by Gronwall's lemma we obtain s u p 0 𝑠 𝑡 𝑢 𝑛 ( 𝑠 ) 𝑢 𝑚 ( 𝑠 ) 2 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 + 𝑚 [ ] e x p ( 𝐶 𝑡 ) , 𝑡 0 , 𝑇 , ( 5 . 2 6 ) consequently s u p 0 𝑠 𝑡 𝑢 𝑛 ( 𝑠 ) 𝑢 𝑚 ( 𝑠 ) 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 + 𝑚 ( 5 . 2 7 ) takes place for all 𝑛 , 𝑚 . This implies that ( 𝑢 𝑛 ( 𝑡 ) ) 𝑛 is a Cauchy sequence in the Banach space 𝐶 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) , and hence it converges in the norm of this latter to some function 𝑢 𝐶 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) . Besides, passing to the limit 𝑚 in (5.27), we obtain the desired error estimate, which finishes the proof.

Now, we present some properties of the obtained solution.

The limit-function 𝑢 from Proposition  5.2, possesses the following properties:

(i) 𝑢 𝐶 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) 𝐿 ( 𝐼 , 𝑉 ) ) ,(ii) 𝑢 is strongly differentiable a.e. in 𝐼 and 𝑑 𝑢 / 𝑑 𝑡 𝐿 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) ,(iii) ̃ 𝑢 𝑛 ( 𝑡 ) 𝑢 ( 𝑡 ) in 𝐵 1 2 ( 0 , 1 ) for all 𝑡 𝐼 ,(iv) 𝑢 𝑛 ( 𝑡 ) , ̃ 𝑢 𝑛 ( 𝑡 ) 𝑢 ( 𝑡 ) in 𝑉 for all 𝑡 𝐼 ,(v) ( 𝑑 𝑢 𝑛 / 𝑑 𝑡 ) ( 𝑡 ) ( 𝑑 𝑢 / 𝑑 𝑡 ) ( 𝑡 ) in 𝐿 2 ( 𝐼 , 𝐵 1 2 ( 0 , 1 ) ) .

Proof. On the basis of estimates (5.5) and (5.6), uniform convergence statement from Proposition  5.2, and the continuous embedding 𝑉 𝐵 1 2 ( 0 , 1 ) , the assertions of the present theorem are a direct consequence of [2, Lemma  1.3.15].

Theorem 5.3. Under Assumptions ( 𝐻 1 ) - ( 𝐻 3 ) , (2.3) admits a unique weak solution, namely, the limit function 𝑢 from Proposition  5.2, in the sense of Definition  2.3.

Proof. We have to show that the limit function 𝑢 satisfies all the conditions (i), (ii), (iii), and (iv) of Definition  2.3. Obviously, in light of the properties of the function 𝑢 listed in Theorem 5.3, the first two conditions of Definition  2.3 are already seen. On the other hand, since 𝑢 𝑛 𝑢 in 𝐶 ( 𝐼 , 𝑉 ) as 𝑛 and, by construction, 𝑢 𝑛 ( 0 ) = 𝑈 0 , it follows that 𝑢 ( 0 ) = 𝑈 0 , so the initial condition is also fulfilled, that is, Definition  2.3(iii) takes place. It remains to see that the integral identity (2.21) is obeyed by 𝑢 . For this, integrating (5.10) over ( 0 , 𝑡 ) and using the fact that 𝑢 𝑛 ( 0 ) = 𝑈 0 , we get 𝑢 𝑛 ( 𝑡 ) 𝑈 0 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑡 0 ̃ 𝑢 𝑛 ( 𝜏 ) , 𝑣 𝑑 𝜏 = 𝑡 0 ̃ 𝑘 𝜏 , ̃ 𝑢 𝑛 ( 𝜏 ) , 𝑣 𝐵 1 2 ( 0 , 1 ) 𝑑 𝜏 + 𝑡 0 𝑓 𝑛 ( 𝜏 ) , 𝑣 𝐵 1 2 ( 0 , 1 ) 𝑑 𝜏 , ( 5 . 2 8 ) consequently, after some rearrangements 𝑢 𝑛 ( 𝑡 ) 𝑈 0 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑡 0 ̃ 𝑢 𝑛 ( = 𝜏 ) , 𝑣 𝑑 𝜏 𝑡 0 𝜏 0 𝑎 ( 𝜏 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝑑 𝑠 , 𝑣 𝐵 1 2 ( 0 , 1 ) 𝑑 𝜏 + 𝑡 0 ( 𝑓 ( 𝜏 ) , 𝑣 ) 𝐵 1 2 ( 0 , 1 ) + 𝑑 𝜏 𝑡 0 ̃ 𝑘 𝜏 , ̃ 𝑢 𝑛 ( 𝜏 ) 𝜏 0 𝑎 ( 𝜏 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝑑 𝑠 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑑 𝜏 𝑡 0 𝑓 𝑛 ( 𝜏 ) 𝑓 ( 𝜏 ) , 𝑣 𝐵 1 2 ( 0 , 1 ) 𝑑 𝜏 . ( 5 . 2 9 ) Let ̂ 𝑠 𝑛 𝐼 𝐼 and ̂ 𝑠 𝑛 𝐼 𝐼 denote the functions ̂ 𝑠 𝑛 𝑡 ( 𝑡 ) = 0 , f o r 𝑡 = 0 , 𝑛 𝑗 1 𝐼 , i n 𝑛 𝑗 , ̃ 𝑠 𝑛 𝑡 ( 𝑡 ) = 0 , f o r 𝑡 = 0 , 𝑛 𝑗 𝐼 , i n 𝑛 𝑗 . ( 5 . 3 0 ) To investigate the desired result, we prove some convergence statements. Using (5.2), (5.4), and (5.30) we have for all 𝑡 ( 𝑡 𝑛 𝑗 1 , 𝑡 𝑛 𝑗 ] ̃ 𝑘 𝑡 , ̃ 𝑢 𝑛 ( 𝑡 ) 𝑡 0 = 𝑎 ( 𝑡 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝑑 𝑠 𝑡 𝑛 𝑗 0 𝑎 𝑡 𝑛 𝑗 ̂ 𝑠 𝑛 𝑘 ( 𝑠 ) ̂ 𝑠 𝑛 ( 𝑠 ) , ̂ 𝑢 𝑛 ( 𝑠 ) 𝑎 ( 𝑡 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝑑 𝑠 + 𝑡 𝑛 𝑗 𝑡 𝑎 ( 𝑡 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝑑 𝑠 . ( 5 . 3 1 ) Taking into account (5.5), (5.9), and assumptions ( 𝐻 1 ) , ( 𝐻 2 ) it follows that 𝑎 𝑡 𝑛 𝑗 ̂ 𝑠 𝑛 𝑘 ( 𝑠 ) ̂ 𝑠 𝑛 ( 𝑠 ) , ̂ 𝑢 𝑛 ( 𝑠 ) 𝑎 ( 𝑡 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 . ( 5 . 3 2 ) Thanks to (5.31) and (5.32) we obtain ̃ 𝑘 𝑡 , ̃ 𝑢 𝑛 ( 𝑡 ) 𝑡 0 𝑎 ( 𝑡 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝑑 𝑠 𝐵 1 2 ( 0 , 1 ) 𝐶 𝑛 . ( 5 . 3 3 ) On the other hand, in view of the assumed Lipschitz continuity of 𝑓 , we have 𝑓 𝑛 ( 𝜏 ) 𝑓 ( 𝜏 ) 𝐵 1 2 ( 0 , 1 ) 𝑓 ̃ 𝑠 𝑛 ( 𝜏 ) 𝑓 ( 𝜏 ) 𝐵 1 2 ( 0 , 1 ) 𝑙 1 𝑛 . ( 5 . 3 4 ) Now, the sequences { ( ̃ 𝑢 𝑛 ( 𝜏 ) , 𝑣 ) } , 𝑓 { ( 𝑛 ( 𝜏 ) , 𝑣 ) 𝐵 1 2 ( 0 , 1 ) } , and ̃ { ( 𝑘 ( 𝜏 , ̃ 𝑢 𝑛 ( 𝜏 ) ) , 𝑣 ) 𝐵 1 2 ( 0 , 1 ) } are uniformly bounded with respect to both 𝜏 and 𝑛 , so the Lebesgue theorem of majorized convergence is applicable to (5.29). Thus, having in mind (5.7), (5.9), (5.33), and (5.34), we derive that 𝑢 ( 𝑡 ) 𝑈 0 , 𝑣 𝐵 1 2 ( 0 , 1 ) + 𝑡 0 ( = 𝑢 ( 𝜏 ) , 𝑣 ) 𝑑 𝜏 𝑡 0 𝜏 0 𝑎 ( 𝜏 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝑑 𝑠 , 𝑣 𝐵 1 2 ( 0 , 1 ) 𝑑 𝜏 + 𝑡 0 ( 𝑓 ( 𝜏 ) , 𝑣 ) 𝐵 1 2 ( 0 , 1 ) 𝑑 𝜏 ( 5 . 3 5 ) takes place for all 𝑣 𝑉 and 𝑡 [ 0 , 𝑇 ] . Finally, differentiating (5.35) with respect to 𝑡 , we get 𝑑 𝑑 𝑡 𝑢 ( 𝑡 ) , 𝑣 𝐵 1 2 ( 0 , 1 ) = + ( 𝑢 ( 𝑡 ) , 𝑣 ) 𝑡 0 𝑎 ( 𝑡 𝑠 ) 𝑘 ( 𝑠 , 𝑢 ( 𝑠 ) ) 𝑑 𝑠 , 𝑣 𝐵 1 2 ( 0 , 1 ) + ( 𝑓 ( 𝑡 ) , 𝑣 ) 𝐵 1 2 ( 0 , 1 ) [ ] . , a . e . 𝑡 0 , 𝑇 ( 5 . 3 6 ) The uniqueness may be argued in the usual manner. Indeed, exploiting an idea in [11], consider 𝑢 1 and 𝑢 2 two different solutions of (2.3), and define 𝑤 = 𝑢 1 𝑢 2 then, we have 𝑑 𝑑 𝑡 𝑤 ( 𝑡 ) , 𝑣 𝐵 1 2 ( 0 , 1 ) + ( 𝑤 ( 𝑡 ) , 𝑣 ) = 𝑡 0 𝑘 𝑎 ( 𝑡 𝑠 ) 𝑠 , 𝑢 1 ( 𝑠 ) 𝑘 𝑠 , 𝑢 2 ( 𝑠 ) 𝑑 𝑠 , 𝑣 𝐵 1 2 ( 0 , 1 ) . ( 5 . 3 7 ) Choosing 𝑣 = 𝑤 ( 𝑡 ) as a test function, with the aid of Cauchy-Schwarz inequality and assumption ( 𝐻 1 ) , we obtain 1 2 𝑑 ( 𝑑 𝑡 𝑤 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) + 𝑤 ( 𝑡 ) 2 𝐶 𝑡 0 𝑘 𝑠 , 𝑢 1 ( 𝑠 ) 𝑘 𝑠 , 𝑢 2 ( 𝑠 ) 𝐵 1 2 ( 0 , 1 ) 𝑑 𝑠 𝑤 ( 𝑡 ) 𝐵 1 2 ( 0 , 1 ) . ( 5 . 3 8 ) Let 𝜉 [ 0 , 𝑝 ] such that 𝑤 ( 𝜉 ) 𝐵 1 2 ( 0 , 1 ) = m a x [ ] 𝑠 0 , 𝑝 𝑤 ( 𝑠 ) 𝐵 1 2 ( 0 , 1 ) , ( 5 . 3 9 ) integrating (5.38) over ( 0 , 𝑝 ) , 0 𝑝 𝑇 , using (5.39), and invoking assumption ( 𝐻 2 ) , we get 𝑝 0 1 2 𝑑 𝑑 𝑡 𝑤 ( 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) + 𝑤 ( 𝑡 ) 2 𝑑 𝑡 𝐶 𝑝 2 𝑤 ( 𝜉 ) 2 𝐵 1 2 ( 0 , 1 ) , ( 5 . 4 0 ) consequently, with the fact that 𝑤 ( 0 ) = 0 𝑝 0 1 2 𝑑 𝑑 𝑡 𝑤 ( 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) + 𝑤 ( 𝑡 ) 2 𝑑 𝑡 𝐶 𝑝 2 𝜉 0 𝑑 𝑑 𝑡 𝑤 ( 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) 𝑑 𝑡 . ( 5 . 4 1 ) Choosing 𝑝 as a constant verifying the condition 𝛼 , 𝑇 = 𝛼 𝑝 , 𝐶 𝑝 2 1 2 , ( 5 . 4 2 ) we have, by virtue of (5.41) 𝑝 0 1 2 𝑑 𝑑 𝑡 𝑤 ( 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) 𝑑 𝑡 + 𝑝 0 𝑤 ( 𝑡 ) 2 𝑑 𝑡 𝜉 0 1 2 𝑑 𝑑 𝑡 𝑤 ( 𝑡 ) 2 𝐵 1 2 ( 0 , 1 ) 𝑑 𝑡 , ( 5 . 4 3 ) taking into account that 𝜉 𝑝 , we obtain [ ] 𝑤 ( 𝑡 ) = 0 , o n 0 , 𝑝 . ( 5 . 4 4 ) Following the same lines as for [ 0 , 𝑝 ] , we deduce that [ ] 𝑤 ( 𝑡 ) = 0 , o n 𝑖 𝑝 , ( 𝑖 + 1 ) 𝑝 , 𝑖 = 1 , 2 , 3 , , ( 5 . 4 5 ) therefore, we derive 𝑤 ( 𝑡 ) 0 , on [ 0 , 𝑇 ] , then 𝑢 1 𝑢 2 . This achives the proof.

References

  1. K. Rektorys, The Method of Discretization in Time and Partial Differential Equations, vol. 4 of Mathematics and Its Applications (East European Series), D. Reidel, Dordrecht, The Netherlands, 1982. View at MathSciNet
  2. J. Kačur, Method of Rothe in Evolution Equations, vol. 80 of Teubner-Texte zur Mathematik, BSB B. G. Teubner Verlagsgesellschaft, Leipzig, Germany, 1985. View at MathSciNet
  3. A. Bouziani, “On the quasi static flexure of a thermoelastic rod,” Communications in Applied Analysis for Theory and Applications, vol. 6, no. 4, pp. 549–568, 2002. View at MathSciNet
  4. A. Bouziani, “Mixed problem with boundary integral conditions for a certain parabolic equation,” Journal of Applied Mathematics and Stochastic Analysis, vol. 9, no. 3, pp. 323–330, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  5. A. Bouziani and N.-E. Benouar, “Sur un problème mixte avec uniquement des conditions aux limites intégrales pour une classe d'équations paraboliques,” Maghreb Mathematical Review, vol. 9, no. 1-2, pp. 55–70, 2000. View at MathSciNet
  6. N. Merazga and A. Bouziani, “Rothe method for a mixed problem with an integral condition for the two-dimensional diffusion equation,” Abstract and Applied Analysis, vol. 2003, no. 16, pp. 899–922, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  7. N. Merazga and A. Bouziani, “Rothe time-discretization method for a nonlocal problem arising in thermoelasticity,” Journal of Applied Mathematics and Stochastic Analysis, vol. 2005, no. 1, pp. 13–28, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  8. N. Merazga and A. Bouziani, “On a time-discretization method for a semilinear heat equation with purely integral conditions in a nonclassical function space,” Nonlinear Analysis, vol. 66, no. 3, pp. 604–623, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. D. Bahaguna, A. K. Pani, and V. Raghavendra, “Rothe's method to semilinear hyperbolic integrodifferential equations,” Journal of Applied Mathematics and Stochastic Analysis, vol. 3, no. 4, pp. 245–252, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  10. D. Bahuguna and V. Raghavendra, “Rothe's method to parabolic integrodifferential equations via abstract integrodifferential equations,” Applicable Analysis, vol. 33, no. 3-4, pp. 153–167, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  11. D. Bahuguna and R. Shukla, “Method of semidiscretization in time for quasilinear integrodifferential equations,” International Journal of Mathematics and Mathematical Sciences, vol. 2004, no. 9, pp. 469–478, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  12. D. Bahuguna and S. K. Srivastava, “Approximation of solutions to evolution integrodifferential equations,” Journal of Applied Mathematics and Stochastic Analysis, vol. 9, no. 3, pp. 315–322, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet