- Department of Computational Mathematics and Programming, Institute of Applied Mathematics and Fundamental Sciences, Lviv Polytechnic National University, Lviv, Ukraine
The initial-boundary and the inverse coefficient problems for the semilinear hyperbolic equation with strong damping are considered in this study. The conditions for the existence and uniqueness of solutions in Sobolev spaces to these problems have been established. The inverse problem involves determining the unknown time-dependent parameter in the right-hand side function of the equation using an additional integral type overdetermination condition.
1 Introduction
Propagation of sound in a viscous gas and other similar processes of the same nature can be described by the model hyperbolic equation of the third order, which includes a mixed derivative with respect to spatial and time variables
where η is a positive constant, and ηΔxut represents low viscosity.
Many important physical phenomena can be modeled with the use of Equation 1 and its generalizations. These are, in particular, processes that occur in viscous media (propagation of disturbances in viscoelastic and viscous-plastic rods, movement of a viscous compressible fluid, sound propagation in a viscous gas), wave processes in different media, acoustic waves in environments where wave propagation disrupts the state of thermodynamic and mechanical equilibrium, liquid filtration processes in porous media, heat transfer in a heterogeneous environment, moisture transfer in soils, and longitudinal vibrations in a homogenous bar with viscosity. The term Δxut indicates that the level of stress is proportional to the level of strains and to the strain rate [1–5].
Due to its wide range of applications, different problems for Equation 1 were investigated by many authors. For example, the unique solvability of the direct initial-boundary value problems for Equation 1 and its nonlinear generalizations with power nonlinearities have been studied in other research [1, 2, 4–11].
The inverse problems, with the integral overdetermination conditions, of identifying of the coefficients in the right-hand side function of hyperbolic equations without damping or for other types of equations have been investigated in many studies [12–18]. Their unique solvability has been solved with the use of the methods such as integral equations, the Green function, regularization, and the Shauder principle [14] and successive approximations [18]. The unique solvability of a two-dimensional inverse problem for the linear third-order hyperbolic equation with constant coefficients and with the unknown time-dependent lower coefficient has been proved in Mehraliyev et al. [19].
The main objective of this study is to determine the sufficient conditions for the existence and uniqueness of the solution to the inverse problem for the third-order semilinear hyperbolic equation with an unknown time- dependent function on its right-hand side. The unknown function is determined from the equation, subject to initial, boundary, and integral type overdetermination conditions. To prove the main results of the study, we use the properties of the solution for the corresponding initial-boundary value problem and the method of successive approximations. These results are new for semilinear n-dimensional third-order hyperbolic equations with non-constant coefficients and an unknown function on their right-hand side. The unique solvability of the initial-boundary value problem has been proved using of the method of Galerkin approximations and the methods of monotonicity and compactness.
2 Problem setting
Let Ω ⊂ ℝn, n ∈ ℕ, be a bounded domain with the smooth boundary ∂Ω ∈ C1 and 0 < T < ∞. Denote Qτ = Ω × (0, τ), τ ∈ (0, T]; Qt1, t2 = Ω × (t1, t2), t1, t2 ∈ (0, T]. In this study, we consider the following inverse problem: find the sufficient conditions for the existence of a pair of functions (u(x, t), g(t)) that satisfies the equation with strong damping (in the sense of Definition 3.1).
and the initial, boundary, and overdetermination conditions
We shall use Lebesgue and Sobolev spaces L∞(·), L2(·), H1(·): = W1, 2(·), Ck(·), C([0, T];L2(G)), (see, e.g., Gajewski et al. [20]).
Suppose that the data of the problem (2–5) satisfy the following conditions.
(H1): aij(x, t) = aji(x, t), bij(x, t) = bji(x, t), and
for all ξ ∈ ℝn, almost all x ∈ Ω, all t ∈ [0, T], and i, j = 1, ..., n, where α0, α1 and β0, β1 are positive constants.
(H2): functions φ1(x, ξ), φ2(x, ξ) are measurable with respect to x ∈ Ω for all ξ ∈ ℝ1 and continuously differentiable concerning ξ ∈ ℝ. Moreover,
for almost all x ∈ Ω and ξ, η ∈ ℝ, where Li,0, Li,1 are positive constants.
(H3): , ,
(H4): E ∈ C2([0, T]), , .
(H5): .
Denote
Let γ0 = γ0(Ω) be a coefficient in Friedrich's inequality.
3 Initial-boundary value problem
Definition 3.1. A function u(x, t) is considered to be a solution of problem 2–4 if , u satisfies (3), and
for all functions and τ ∈ (0, T].
Theorem 3.2. Under the assumptions (H1)–(H3) and g ∈ L2(0, T), aijt ≤ 0 for all i, j = 1, 2, …, n, the problem (2–4) has a unique solution.
Proof. First, using Galerkin method, we prove the existence of a solution for the problem. Let be a basis in , orthonormal in L2(Ω). We will consider the sequence of functions
where the set is a solution of the initial value problem
Here and
The solution of system (8) exists on some interval [0, τ0] (Carathéodory's Theorem [21, p. 43]). The estimation (13) from below implies that this solution could be extended on [0, T]. Multiplying each equation of (8) on function respectively, summing up for k from 1 to N and integrating to t on interval [0, τ], τ ≤ τ0, we obtain
After transformations of terms from (9), we get
Note that
then from (10) we obtain
We rewrite the last inequality in the form
where
Then by Grönwall's lemma, from (12), we get
Therefore, from (11) we also get
Multiplying each equation of (8) on function respectively, summing up with respect to k from 1 to N and integrating on interval [0, τ], τ ≤ T, we obtain
After transformations in all terms from (15), we get
where Taking into account (13), (14), from (16) we obtain
where
The right-hand sides of the estimates (13), (14), and (17) are positive constants, independent of N. Therefore, there exists a subsequence of (which will be denoted by the same notation), such that as N → ∞
It follows from (18) that and therefore, as N → ∞. Besides, , and weakly in
Equations 8 and 18 imply the equality
for all functions and τ ∈ (0, T].
Let us prove that χ = φ2(x, ut).
Note that for all k ∈ ℕ. Due to (18), is fundamental in . So, for any ε > 0, there exists such a number N0 that for all N, k ∈ ℕ, N > N0 the inequality holds; thus, is also fundamental in and, therefore,
Consider the sequence
where . From (9), it follows that
After substitution (22) in (21), passing to the limit as N → ∞, taking into account (18), (20), and the assumptions of Theorem 3.2, we obtain
Choosing here , dividing the result on ϰ and then tending ϰ → 0 we obtain χ = φ2(x, ut). Hence, from (19), it follows (7).
Now we prove the uniqueness of the solution for the problems (2–4). On the contrary, suppose that there exist two solutions u(1)(x, t) and u(2)(x, t) of problems (2–4). Then ũ: = ũ(x, t) = u(1)(x, t)−u(2)(x, t) satisfies the conditions ũ(x, 0) ≡ 0, ũt(x, 0) ≡ 0, and the equality
holds for all τ ∈ (0, T].
After choosing v = ũt in (23) we get
From (24) by the same way as from (11) we got (12), we find the following estimate
Then from Grönwall's lemma and (25) we obtain and hence, ũ ≡ 0, and, therefore, u(1) = u(2) in QT.
4 Inverse problem
Definition 4.1. A pair of functions (u(x, t), g(t)) is a solution to the problem (2–5), if , and g ∈ C([0, T]), and it satisfies (5) and
holds for all functions and τ ∈ (0, T].
4.1 The equivalent problem
In this section, we shall find the equivalent problem for the problem (2–5).
Lemma 4.2. Let the assumptions of Theorem 3.2, (H4), and (H5) hold. A pair of functions (u(x, t), g(t)), where , g ∈ C([0, T]), is a solution to the problem (2–5) if and only if it satisfies (26) for all functions and for τ ∈ (0, T], the equality
holds for t ∈ [0, T].
Proof. Necessity: Let (u(x, t), g(t)) be a solution to the problem (2–5). From (26) and (5), it follows that
By integrating by parts in (28) and using the condition (H4), we get the equality
From (29), we can obtain (27).
Sufficiency: Let g* ∈ C([0, T]), , and they satisfy (4), (26), and (27). Then u* is a solution to the problem (2–4) with g* instead of g in (2).
We set t ≥ 0. In exactly the same way as in the proof of necessity, we obtain
On the other hand g*(t) and u*(x, t) satisfy (27)
It follows from (30), (31) that
Integrating (32) with the use of the equalities , implies E*(t) = E(t), t ≥ 0. Hence, u*(x, t) satisfies the overdetermination condition (5).
4.2 Main results
Let Denote
Theorem 4.3. Let aijt ≤ 0 for all i, j = 1, 2, …, n, and the assumptions (H1) – (H5) hold. Then there exists a unique solution to the problem (2–5).
Proof. I. In the first step, we shall prove the theorem for T ≤ T0.
We construct an approximation (um(x, t), gm(t)) of the solution of problem (2–5), where g1(t): = 0, the functions gm(t), m ≥ 2, satisfy the equality
and um satisfies the equality
for all , and the conditions
It follows from Theorem 3.2 that for each m ∈ ℕ there exists a unique function , that satisfies (34), (35). Now we show that converges to the solution of the problem (2–5). Denote
Equation 33 for t ∈ (0, T0] and m ≥ 3, implies the equality
We square both sides of equality (36) and integrate the result with respect to t, taking into account the hypotheses (H5), then we obtain
Then (37) implies the estimate
It follows from (35) that Hence, from (34) with , τ ∈ (0, T0], we get
The last term in (39)
Besides,
and
Then, taking into account (H1) – (H5), from (39) we get inequality
According to Grönwall's Lemma, we obtain
Interating (41) with respect to τ, we get the estimate
Besides, (40) and (42) for m ≥ 2 imply the estimates
and
where
Note that r2(t) = g2(t). Then, taking into account (33), we have
where M6 is a positive constant. It follows from (42) and (38) that for m ≥ 3
By using (45) and the assumption M4T0 < 1, we can show the estimate
Due to (42)
Besides, (43) and (44) for m ≥ 2 imply the estimates
and
And, therefore,
and
and for k ∈ ℕ, m ≥ 2
Besides, we square both sides of equality (36), taking into account the hypotheses (H5) and obtain
From (53), we can conclude that:
Therefore,
It follows from (46), (50), (51), (52), and (55) that for any ε > 0, there exists m0 such that for all k, m ∈ ℕ, m > m0, the following inequalities hold:
are true. Hence, the sequence is fundamental in and in C([0, T0]), is fundamental in and in , and is fundamental in and in . Therefore, as m → ∞
Theorem 3.2 implies the following estimate
and, by virtue of (56) , where M7 is independent on m, and therefore the right-hand side of (57) is bounded with the constant, independent on m. Hence, we can select a subsequence of sequence (we preserve the same notation for this subsequence), such that
Taking into account (56) and (58), from (33), (34) we get that the pair of functions (u(x, t), g(t)) satisfies (27) and (26). By virtue of Lemma 4.2 (u(x, t), g(t)) is a solution of the problem (2–5) in QT0.
II. Uniqueness of solution of the problem (2–5), with T ≤ T0.
Assume that (u(1)(x, t), g(1)(t)) and (u(2)(x, t), g(2)(t)) be two solutions of problem (2–5). Then the pair of functions , where ũ(x, t) = u(1)(x, t)−u(2)(x, t), satisfies the conditions ũ(x, 0) ≡ 0, ũt(x, 0) ≡ 0, the equality
for all and the equality
holds.
After choosing v = ũt in (59) we get
It is easy to get from (60) and (H5) inequalities
From (61) by the same way as from (39) we got (42), we find the following estimate
and taking into account (62) from (63), we obtain Since M4T0 < 1, we conclude that hence, u(1) = u(2) in QT0. Then (62) implies and, therefore, g(1)(t) ≡ g(2)(t) on [0, T0].
III. Let now T > 0 be arbitrary number.
Let us divide the interval [0, T] into a finite number of intervals [0, T1], [T1, 2T1], …, [(N − 1)T1, NT1], where NT1 = T, and T1 ≤ T0. According to I and II, there exists a unique solution (u1(x, t), g0,1(t)) to the problem (2–5) in the domain QT1.
Now, we will prove that there exists a unique solution in the domain QT1, 2T1: = Ω × (T1; 2T1) for the problem for the Equation 2 with conditions (4) and (5) as t ∈ [T1; 2T1], and with the initial condition u(x, T1) = u1(x, T1), ut(x, T1) = u1t(x, T1), and x ∈ Ω.
Let us change the variables t = τ+T1, τ ∈ [0;T1] in this problem. We will denote G0(τ) = g(τ+T1), U(x, τ) = u(x, τ+T1), and For the pair (U(x, τ), G0(τ)) we obtain the problem:
It is obvious that all coefficients of the Equation 64 and the functions E(1)(τ) satisfy the same conditions as the functions from (2) and (5). According to I and II, there exists a unique solution to the problem (64–67) in QT1, and, thus for the problems for the Equation 2 with conditions (4) and (5) as t ∈ [T1; 2T1] and with the initial condition u(x, T1) = u1(x, T1), ut(x, T1) = u1t(x, T1), and x ∈ Ω, in the domain QT1, 2T1. Denote it by (u2(x, t), g0,2(t)). By following similar reasoning on the intervals [2T1; 3T1], …, [(N − 1)T1; NT1], we can prove the existence and uniqueness of weak solutions (uk(x, kt), g0,k(t)), k = 3, …, N, in the domain Q(k−1)T1,kT1: = Ω × ((k − 1)T1, kT1) of the inverse problem for the Equation 2 with conditions (4) and (5) as t ∈ [(k − 1)T1; kT1] and the initial condition u(x, (k − 1)T1) = uk−1(x, (k − 1)T1), x ∈ Ω. It is clear that a pair of functions (u(x, t), g0(t)), where
is a solution for the problem (2–5) in the domain QT.
IV. The uniqueness of solution is proved similar as in II, III: Assume that (u(1)(x, t), g(1)(t)) and (u(2)(x, t), g(2)(t)) be two solutions of problem (2–5). Then the pair of functions , where ũ(x, t) = u(1)(x, t)−u(2)(x, t), satisfies the conditions ũ(x, 0) ≡ 0, ũt(x, 0) ≡ 0, the equality
for all , τ ∈ (0, T], and the equality
holds.
Let us divide the interval [0, T] into a finite number of intervals [0, T1], [T1, 2T1], …, [(N − 1)T1, NT1], where NT1 = T, and T1 ≤ T0.
Let us choose τ ∈ [0, T1] in (68). After choosing here v = ũt, we get
It is easy to get from (69) and (H5) inequalities
From (70), by the same way as from (39), we got (42). We find the following estimate:
and taking into account (71) from (72), we obtain Since M4T1 < 1, we conclude that hence, u(1) = u(2) in QT1. Then (71) implies and, therefore, g(1)(t) ≡ g(2)(t) on [0, T1].
Let us choose τ ∈ [0, 2T1] in (68). After choosing here v = ũt, we get
Note that ũ ≡ 0 in QT1 and on [0;T1], therefore, from (73) it follows that
It is easy to get from (69) and (H5) inequalities
From (74), by the same way as from (39), we got (42). We find the following estimate
and taking into account (75) from (76), we obtain Since M4T1 < 1, we conclude that hence, u(1) = u(2) in QT1, 2T1. Then, (75) implies and, therefore, g(1)(t) ≡ g(2)(t) on [T1, 2T1]. Therefore, u(1) = u(2) in Q2T1, g(1)(t) ≡ g(2)(t) on [0, 2T1].
Considering τ ∈ [0, 3T1], … , τ ∈ [0, NT1] in (68), by the same arguments, we find that u(1) = u(2) in Q(k−1)T1,kT1, g(1)(t) ≡ g(2)(t) on [(k − 1)T1, kT1], k = 1, 2, …, N. Therefore, u(1) = u(2) in QT, g(1)(t) ≡ g(2)(t) on [0, T].
5 Conclusions
In this study, we have derived the necessary conditions for the existence and the uniqueness of the solution for the initial-boundary value problem, as well as the inverse problem, for semilinear hyperbolic equation of the third order with an unknown parameter in its right-hand side function. To determine the unknown function, an additional integral-type overdetermination condition have been introduced. These results were obtained by utilizing the properties of the solutions to the initial-boundary value problem and the method of successive approximations.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found at: no.
Author contributions
NP: Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Hai DD. On a strongly damped quasilinear wave equation. Demonstr Math. (1986) 19:327–40. doi: 10.1515/dema-1986-0208
2. Angt DD, Dinh APN. On the strongly damped wave equation: utt − Δu − Δut + f(u) = 0. Siam J Math Anal. (1988) 19:1409–18. doi: 10.1137/0519103
3. Dvornyk AV, Tkachenko VI. Almost periodic solutions of the wave equation with damping and impulsive action. Ukr Math J. (2023) 75:68–79. doi: 10.1007/s11253-023-02186-7
4. Greenberg JM, McCamy RC, Mizel V. On the existence, uniqueness and stability of solutions of the equation . J Math Mech. (1968) 17:707–28.
5. Graham A. On the existence of solutions to the equation utt = uxxt + σ(ux). J Differ Equ. (1980) 35:200–31. doi: 10.1016/0022-0396(80)90040-6
6. Carvalho AN, Cholewa JW. Strongly damped wave equations in W1, p(Ω) × Lp(Ω). Discrete Contin Dyn Syst. (2007) 2007:230–9. doi: 10.3934/proc.2007.2007.230
7. Dafermos CM. The mixed initial-boundary value problem for the equations of nonlinear one-dimensional viscoelasticity. J Differ Eq. (1969) 6:71–86. doi: 10.1016/0022-0396(69)90118-1
8. Domanska H, Lavrenyuk S, Protsakh N. Problem for nonlinear third order hyperbolic equation. Nauk Bull Chernivtsi Univ Math. (2005) 269:34–42.
9. Ferreira J, Pereira DC. On a nonlinear degenerate evolution equation with strong damping. Int J Math Math Sci. (1992) 15:543–52. doi: 10.1155/S016117129200070X
10. Gao H, Ma TF. Global solutions for a nonlinear wave equation with the p-Laplacian operator. EJQTDE. (1999) 11:1–13. doi: 10.14232/ejqtde.1999.1.11
11. Protsakh N. Mixed problem for anisotropic third order equation. Math Methods Fiz-Mech Fields. (2007) 50:40–53.
12. Akhundov AY, Habibova AS. On the solvability of an inverse problem for a hyperbolic heat equation. Azerb J Math. (2023) 13:205–15.
13. Oussaeif T-E, Bouziani A. Inverse problem for a hyperbolic equation with an integral overdetermination condition. Electron J Differ Eq. (2016) 2016:1–7.
14. Lopushanska H, Lopushansky A. Inverse problem with a time-integral condition for a fractional diffusion equation. Math Methods Appl Sci. (2019) 42:3327–40. doi: 10.1002/mma.5587
15. Protsakh N. Determining of three unknown functions of a semilinear ultraparabolic equation. Math Methods Appl Sci. (2021) 44:617–33. doi: 10.1002/mma.6768
16. Protsakh N. Inverse problem for higher order ultraparabolic equation with unknown minor coefficient and right-hand side function. Miskolc Mathem Notes. (2020) 21:335–50. doi: 10.18514/MMN.2020.2879
17. Protsakh NP. Inverse problem for a semilinear ultraparabolic equation with unknown right-hand side. Ukr Math J. 66:371–90. doi: 10.1007/s11253-014-0937-4
18. Protsakh NP, Ptashnyk BY. Nonlinear Ultraparabolic Equations and Variational Inequalities. Kyiv: Naukova dumka. (2017).
19. Mehraliyev YT, Sadikhzade SR, Ramazanova AT. Two-dimensional inverse boundary value problem for a third-order pseudo-hyperbolic equation with an additional integral condition. Eur J Pure Appl Math. (2023) 16:670–86. doi: 10.29020/nybg.ejpam.v16i2.4743
20. Gajewski H, Gröger K, Zacharias K. Nichtlineare Operatorgleichungen und Operatordifferentialgleichungen. Berlin: Akademie-Verlag, Mathematische Lehrbücher und Monographien (1974). doi: 10.1515/9783112717899
Keywords: semilinear wave equation, inverse coefficient problem, existence of solutions, initial-boundary value problem, strong damping
Citation: Protsakh N (2024) Inverse problem for semilinear wave equation with strong damping. Front. Appl. Math. Stat. 10:1467441. doi: 10.3389/fams.2024.1467441
Received: 19 July 2024; Accepted: 27 August 2024;
Published: 24 September 2024.
Edited by:
Tamara Fastovska, V. N. Karazin Kharkiv National University, UkraineReviewed by:
Aysel Ramazanova, University of Duisburg-Essen, GermanyMansur I. Ismailov, Gebze Technical University, Türkiye
Copyright © 2024 Protsakh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nataliya Protsakh, protsakh@ukr.net