Representation of Classical Solutions to a Linear Wave Equation with Pure Delay
Abstract
For a wave equation with pure delay, we study an inhomogeneous initialboundary value problem in a bounded 1D domain. Under smoothness assumptions, we prove unique existence of classical solutions for any given finite time horizon and give their explicit representation. Continuous dependence on the data in a weak extrapolated norm is also shown.
1 Introduction
The wave equation is a typical linear hyperbolic secondorder partial differential equation which naturally arises when modeling phenomena of continuum mechanics such as sound, light, water or other kind of waves in acoustics, (electro)magnetics, elasticity and fluid dynamics, etc. (cf. [6, 13]). Providing a rather adequate description of physical processes, partial differential equations, or equations with distributed parameters in general, have found numerous applications in mechanics, medicine, ecology, etc. Introducing aftereffects such as delay into these equations has gained a lot of attention over several past decades. See e.g., [2, 3, 7, 8]. Mathematical treatment of such systems requires additional carefulness since distributed systems with delay often turn out to be even illposed (cf. [4, 5, 12]).
In the present paper, we consider an initialboundary value problem for a general linear wave equation with pure delay and constant coefficients in a bounded interval subject to nonhomogeneous Dirichlet boundary conditions. To solve the equation, we employ Fourier’s separation method as well as the special functions referred to as delay sine and cosine functions which were introduced in [9, 10]. We prove the existence of a unique classical solution on any finite time interval, show its continuous dependence on the data in a very weak extrapolated norm, give its representation as a Fourier series and prove its absolute and uniform convergence under certain conditions on the data.
2 Equation with pure delay
For , , we consider the following linear wave equation in a bounded interval with a single delay being a second order partial differencedifferential equation for an unknown function
(2.1) 
subject to nonhomogeneous Dirichlet boundary conditions and initial conditions
(2.2) 
Since we are interested in studying classical solutions, the following compatibility conditions are required to assure for smoothness of the solution on the boundary of timespace cylinder
Definition 2.1.
Remark 2.2.
The previous definition does not impose any continuity of time derivatives in . If the continuity or even smoothness are desired, additional compatibility conditions on the data, including , are required.
Let , , denote the standard Sobolev norm (cf. [1]) and denote the norm of corresponding negative Sobolev space. We introduce the norm and define the Hilbert space as a completion of with respect to . Obviously, , i.e., can be continuously embedded into the space of distributions.
With this notation, we easily see that (with denoting the distributional derivative) is a bounded linear operator on since
Theorem 2.3.
There exists a constant , dependent only on , such that the estimate
Proof.
Let be the classical solution to Equations (2.1), (2.2). We define
Then satisfies homogeneous Dirichlet boundary conditions and solves the equation
(2.3) 
in the extrapolated space with
We multiply the equation with in the scalar product of and use Young’s inequality to get the estimate
(2.4) 
As in [11], we introduce the history variable
and obtain
Multiplying these identities with in and performing a partial integration, we find
(2.5) 
Adding equations (2.4) and (2.5) to the trivial identity
we obtain
Thus, we have shown
(2.6) 
where
From Equation (2.6) we conclude
Using now the integral form of Gronwall’s inequality, we obtain
(2.7) 
for certain . Taking into account
for some constants and exploiting the definition of , the proof is a direct consequence of Equation (2.7). ∎
Corollary 2.4.
Remark 2.5.
Next, we want to establish conditions on the data allowing for the existence of a classical solution. Performing the substitution
(2.8) 
with a new unknown function (cp. [11]), the initial boundary value problem (2.1), (2.2) can be written in the following simplified form with a selfadjoint operator on the righthand side
(2.9) 
with complemented by the following boundary and initial conditions
(2.10)  
(2.11) 
and
The solution will be determined in the form
Here, is an arbitrary function with satisfying the boundary conditions
Assuming , we let
(2.12) 

solves the homogeneous equation
(2.13) subject to homogeneous boundary and nonhomogeneous initial conditions
(2.14) In particular, with the function selected as in Equation (2.12), we obtain
(2.15) 
solves the nonhomogeneous equation
(2.16) with
(2.17) subject to homogeneous boundary and initial conditions. For from Equation (2.12), we have
(2.18)
3 Homogeneous equation
In this section, we obtain a formal solution to the initialboundary value problem (2.13) with initial and boundary conditions given in Equations (2.10), (2.11). We exploit Fourier’s separation method to determine in the product form . After plugging this ansatz into Equation (2.13), we find
Hence,
By formally separating the variables, we deduce
Thus, the equation can be decoupled as follows
(3.1) 
These are linear second order ordinary (delay) differential equations with constant coefficients.
Due to the zero boundary conditions for , the boundary conditions for the second equation in (3.1) will also be homogeneous, i.e.,
Therefore, we obtain a Sturm & Liouville problem admitting nontrivial solutions only for the eigennumbers
and the corresponding eigenfunctions
Assuming
we denote
and consider the first equation in (3.1), i.e.,
(3.2) 
The initial conditions for each of the equations in (3.2) can be obtained by expanding the initial data into a Fourier series with respect to the eigenfunction basis of the second equation in (3.1)
(3.3) 
for . Let us further determine the solution of the Cauchy problem associated with each of the equations in (3.2) subject to the initial conditions from Equation (3.3).
First, we briefly present some useful results from the theory of second order delay differential equations with pure delay obtained in [9]. The authors considered a linear homogeneous second order delay differential equation
(3.4) 
They introduced two special functions referred to as delay cosine and sine functions. Exploiting these functions, a unique solution to the initial value problem (3.4) was obtained.
Definition 3.1.
Delay cosine is the function given as
(3.5) 
with order polynomials on each of the intervals continuously adjusted at the nodes , .
Definition 3.2.
Delay sine is the function given as
(3.6) 
with order polynomials on each of the intervals continuously adjusted at the nodes , .
There has further been proved that delay cosine uniquely solves the linear homogeneous second order ordinary delay differential equation with pure delay subject to the unit initial conditions in and the delay sine in its turn solves Equation (3.4) subject to the initial conditions for .
Using the fact above, the solution of the Cauchy problem was represented in the integral form. In particular, the solution to the homogeneous delay differential equation (3.4) with the initial conditions in for an arbitrary was shown to be given as
(3.7) 
Turning back to the delay differential equation (3.2) with the initial conditions (2.4), we obtain their unique solution in the form
(3.8) 
Thus, assuming sufficient smoothness of the data to be specified later, the solution to the homogeneous equation (2.13) satisfying homogeneous boundary and nonhomogeneous initial conditions in reads as
(3.9) 
4 Nonhomogeneous equation
Next, we consider the nonhomogeneous equation (2.16) with the righthand side from Equation (2.18) subject to homogeneous initial and boundary conditions
where
The solution will be constructed as as Fourier series with respect to the eigenfunctions of the Sturm & Liouville problem from the previous section, i.e.,
(4.1) 
Plugging the ansatz from (4.1) into Equation (2.6) and comparing the timedependent Fourier coefficients, we obtain a system of countably many second order delay differential equations
(4.2) 
In [9], the initial value problem for the nonhomogeneous delay differential equation
with homogeneous initial conditions in was shown to be uniquely solved by
(4.3) 
5 General case solution
The solution in the general case can thus formally be represented as the following series
(5.1) 
5.1 Convergence of the Fourier series
Theorem 5.1.
Let , and . Further, let the data functions , , and be such that their Fourier coefficients and given in Equations (3.3) and (4.5) satisfy the conditions
(5.2) 
for an arbitrary, but fixed . Let
Then the classical solution to problem (2.9)–(2.11) can be represented as an absolutely and uniformly convergent Fourier series given in Equation (5.1). The latter series is a twice continuously differentiable function with respect to both variables. Its derivatives of order less or equal two with respect to and can be obtained by a termwise differentiation of the series and the resulting series are also absolutely and uniformly convergent in .
Proof.
We regroup the series from Equation (5.1) into the following sum
where
and

First, we consider the coefficient functions . For an arbitrary with , we find
If and , , are such that the condition
holds true, the series as well as its derivatives of order less or equal 2 converge absolutely and uniformly. Note that a single differentiation with respect to corresponds, roughly speaking, to a multiplication with .

Next, we consider the coefficients . For an arbitrary with , we perform the substitution and exploit the mean value theorem to estimate
Applying the theorem on differentiation under the integral sign to and taking into account that is twice weakly differentiable in , namely: , its derivatives are polynomials of order lower than those of and their convolution with is continuous, analogous estimates can be obtained for and which, in their turn, also follow to be continuous functions.
Now, if the condition
is satisfied, the series as well as its derivatives of order less or equal 2 converge absolutely and uniformly.

Finally, we look at the Fourier coefficients . Again, for an arbitrary period of time with , , we substitute . Once again, using the mean value theorem, we estimate
As before, can be shown to be twice continuously differentiable. If now
(5.3) is satisfied, then both and its derivatives of order less or equal 2 converge absolutely and uniformly.
Since all three conditions are guaranteed by the assumptions of the Theorem due to the fact , the proof is finished. ∎
References
 [1] Adams, R. A. Sobolev spaces, Pure and Applied Mathematics, Vol. 65, New YorkLondon: Academic Press, 1975
 [2] Azizbayov, E. I., Khusainov, D. Ya. Solution to a heat equation with delay (in Russian), Bulletin of Taras Shevchenko National University of Kyiv, Series: Cybernetics, 12, 2012, pp. 4–14
 [3] Bátkai, A., Piazzera, S. Semigroups for Delay Equations, Resarch Notes in Mathematics, 10, A.K. Peters: Wellesley MA, 2005
 [4] Datko, R. Two examples of illposedness with respect to time delays revisited. IEEE Trans. Automatic Control, 50, 1997, pp. 1374 – 1379
 [5] Dreher, M., Quintanilla, R., Racke, R. Illposed problems in thermomechanics. Appl. Math. Lett., 22(9), 2009, pp. 1374–1379
 [6] Eck, Ch., Garcke, H., Knabber, P. Mathematische Modellierung. SpringerVerlag Berlin Heidelberg, 2008
 [7] Els’gol’ts, L. E., Norkin, S. B. Introduction to the Theory and Application of Differential Equations with Deviating Arguments. Mathematics in Science and Engineering, 105, Academic Press, 1973, pp. 1–357
 [8] Hale, J. K. Theory of Functional Differential Equations, Applied Mathematical Sciences Series, 3, 1977, pp. 1–365
 [9] Khusainov, D. Ya., Diblík, J., Růžičkova, M., Lukáčová, J. Representation of a solution of the Cauchy problem for an oscillating system with pure delay. Nonlinear Oscillations, 11(2), 2008, pp. 276–285
 [10] Khusainov, D. Ya., Ivanov, A. F., Kovarzh, I. V. The solution of wave equation with delay (in Ukrainian with English summary), Bulletin of Taras Shevchenko National University of Kyiv. Series: Physics and Mathematics, 4, 2006, pp. 243–248
 [11] Khusainov, D. Ya., Pokojovy, M., Azizbayov, E, On Classical Solvability for a Linear 1D Heat Equation with Constant Delay, to appear in: Journal of Computational and Applied Mathematics, 2013
 [12] Racke, R. Instability of coupled systems with delay, Commun. Pure. Appl. Anal., 11(5), 2012, pp. 1753–1773
 [13] Tikhonov, A. N., Samarskii, A. A. Equations of Mathematical Physics. Dover Publications, 1990, pp. 1–765