Variational Methods in Semilinear Wave Equation with Nonlinear Boundary Conditions and Stability Questions

We present variational approach to the semilinear equation of the vibrating string xtt(t,y)-Δx(t,y)+l(t,y,x(t,y))=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_{tt}(t,y)-{\Delta } x(t,y)+l(t,y,x(t,y))=0$$\end{document} in bounded domain and certain type of nonlinearity on the boundary. To this effect we derive new dual variational method. Next the question of stability of solutions with respect to initial conditions is discussed.


Introduction
We assume Ω to be an open bounded domain in R n , n ≥ 2, with sufficiently smooth boundary Γ . In the paper we study, over a finite interval [0, T ], T < ∞ (T -arbitrary but fixed) the second order semilinear wave equation with nonlinearity containing source terms in the interior of the domain and dissipation and source terms on its boundary: x tt (t, y) − x(t, y) + l(t, y, x(t, y)) = 0, in (0, T ) × Ω, ∂ ν x(t, y) + x(t, y) + g(t, y, x t (t, y)) = h(t, y, x(t, y)) on Σ = (0, T ) × Γ, Studying the above problems have long history see e.g. [2][3][4]6,[8][9][10]. In most of these papers nonlinearities do not depend on (t, y) and are of special type and they concern well-posedness of the system given by (1) on the finite energy space, i.e. H 1 (Ω)×L2 (Ω). The well-posedness included existence and uniqueness of both local and global solutions. However, as it is shown in [9] if we want to admit more general nonlinearities then we cannot expect well-posedness in Hadamard sense. The main difficulty and the novelty of the problems considered here relate to the presence of the boundary nonlinear term h (t, y, x with the fact that Lopatinski condition does not hold for the Neumann problem. The above translates into the fact that in the absence of the damping, the linear map h → (x(t), x t (t)), is not bounded from L 2 (Σ) → H 1 (Ω) × L2 (Ω), unless the dimension of Ω is equal to one. More details on this problem are discussed in [4,9]. In this paper we study a global existence question over a finite interval [0, T ], T < ∞ for different types of nonlinearities l, g, h, than in [2,4,9,20] and depending on (t, y) and so we are not interested in well-posedness in Hadamard sense. But, we will prove stability of the system with respect to initial conditions i.e. continuous dependence (in some new sense) with respect to initial conditions.
It is well known that in general the nonlinearities l, h may drive the solution of (1) to blow up in finite time [2,12,14,23,[25][26][27][28] (see also discussion in [9] in context of appearance of the dissipation g). In order to exclude such a case we impose some interaction between structure of l, quantity of initial condition x 0 (·) ∈ H 1 (Ω) and length of time T (see the assumption As). The importance of the problem with nonlinearity on the boundary appears in optimal control theory see e.g. [9,10,[15][16][17] and the references therein. The problems like (1) were studied mostly by topological method, semigroup theory or monotone operator theory (see [3,9,10] for discussion about that). As it is known to the author there are no papers studying (1) (with nonlinearity on the boundary) by variational method except the author's papers [19][20][21][22] where a different type of nonlinearities were investigated (linear combination of monotone function or difference of two functions with strong assumptions). We would like to stress that just using the variational approach the dissipation term g on the boundary has no influence on the existence of solution to (1). This is why we do not discuss the interaction between dissipation and source on the boundary as it is usual done while approaching other methods (see e.g. [2,4,6,9,24]). What is essential in the method used here is that, we consider first the equation ∂ ν x(t, y) + x(t, y) + g(t, y, x t (t, y)) = h(t, y, x(t, y)) on Σ = (0, T ) × Γ and put H = 1 2 x 2 , which is a convex function. Convexity is exploited in the paper in different forms and just convexity eliminates other assumptions. The general nonlinearity of the interior source has the form: where F, G are C 1 and F is additionally a convex functions with respect to third variable. It turns out that in this case a size of initial conditions is not essential, however some restrictions on initial conditions are hidden in (HΓ ) and (As) below. We shall study (1) by variational method. However, we shall not consider (1) as the Euler-Lagrange equations of any functional. Instead, we shall investigate two functionals one in interior domain and the second on the boundary y, x t (t, y)) + h(t, y, x(t, y))))dydt (4) defined on some subspaces of the space C([0, T ]; H 1 (Ω)) and C([0, T ]; H 1 (Γ )), respectively, discussed below. H * is the Fenchel conjugate to H . Just studying functionals (3) and (4) we do not need to use energy functionals for system (1), which are standard tools to study relations between nonlinearity, length of time T and quantity of initial condition x 0 . Our purpose is to investigate (1) by studying critical points of functional (3). To this effect we apply a new duality approach. As it is easy to see, functional (3) is unbounded in C([0, T ]; H 1 (Ω)) and this is a reason for which we are looking for critical points of J F different than extremum or minmax type. Our aim is to find a nonlinear subspace X F of C([0, T ]; H 1 (Ω)) and study (3) just only on X F . In last section we prove stability of solutions of (1) when x 0 (·) and x 1 (·) are changing in a suitable way.
First we will study the equation with the help of the functional J Γ (x) (see (4))-Sect. 3. Next we solve where x Γ is a solution of (5), with the help of the functional J F (x) (see (3))-Sect. 4 for the case l(t, y,

Main Results
We will focus on the case when n ≥ 3 ( n = 2 being the least interesting, since the concept of criticality of Sobolev's embedding is much less pronounced). First we formulate results concerning problem (5), its relation (see [1]) to the functional J Γ (x) and problem (6). To this effect we need to explain what we mean by x 0 (·) ∈ H 1 (Γ ) in (5), x 0 (·) ∈ H 1 (Ω) in (6) and the normal derivative ∂ ν x(t, y) in (5) and in (4). To this effect define H 1 (Γ ) = x ∈ H 1 (Ω) : x | Γ ∈ H 1 (Γ ) and we shall consider the space C([0, T ]; H 1 (Γ )). Therefore by the normal derivative ∂ ν x(t, y) in (5) and also in (4) we mean normal derivative of a function x(·, We use a set and the functional y, x t (t, y)) + h(t, y, x(t, y))))dy considered on U t , t ∈ [0, T ] and H (t, y, x) = 1 2 x 2 . That means we are looking for solutions to (5) in U t , t ∈ [0, T ]. Therefore, if such a solution belonging to U t , t ∈ [0, T ], exists it is a strong solution to (5).

Remark 1
The unique restriction on nonlinearity is that only the growth on g and h are imposed. Moreover we would like to stress that x 0 ∈ H 1 (Γ ), thus we impose a little bit more regularity on initial condition x 0 ∈ H 1 (Ω).

Remark 2
It is worth to note that the damping term g may be zero i.e. it has not influence on existence of solutions to (5) in U t , t ∈ [0, T ]. Usually (see [4,5,9]), damping term g is monotonic and source term h is smooth but not necessarily monotonic, in our case we have, the source term of the form x − h and dissipation term g is only continuous and satisfying some growth condition. The reason is that we use quite different approach to solve (1). We would like to stress, once more, that just convexity of x → H (t, y, x) = 1 2 x 2 is very important here and plays a crucial role in the proof-it eliminates a role of damping term g. Just different method used here allows significantly relaxed assumptions concerning problem (1). We assume only that the term h is only continuous in x and satisfies the growth condition (8). Moreover we admit (this is done first time -to the best knowledge of the author) dependence of h on (t, y) (compare [5,20] and literature therein). In our case solutions of (5) are independent of solutions of (6). It is necessary to stress that solution to (5) exists on each finite interval [0, T ], as long as assumption (HΓ ) is satisfied on this interval, i.e. our assumptions depend on the interval [0, T ].
We introduce the definition of a weak solution to (6): Definition 1 (weak solution) By a weak solution to (6), defined on a given interval [0, T ], we mean a function x ∈ U , where We shall consider U with a topology induced by the norm

Remark 3
The weak solution considered in this paper is stronger than e.g. in [5] where The main contribution of this paper is a relaxation of assumptions on nonlinearities l and h. Moreover, in spite of lack of uniqueness results we have global (for given finite interval [0, T ]) existence theorem to (1) and continuous dependence of solutions with respect to initial data in the following sense: respectively, there is a subsequence of {x n }-solutions to (1) corresponding to x 0 n , x 1 n , which we denote again by {x n } weakly convergent in H 1 ((0, T )×Ω) and strongly in L2 ((0, T ) × Ω) to an element x ∈ U being a solution to (1) corresponding to x 0 , x 1 .
To formulate assumptions concerning equation (6) we need to recall some theorems from [13]-linear case of (6) which we use as a starting point to study (6).
Then there exists x being a unique weak solution to Shortly the solution x from the above theorem may be estimated by where ) and in the paper we will just use the last estimations. Everywhere below constants C and D will always denote those occurring in (10), (11). Note that then Assumptions concerning Eq. (1).
, a > 0, (v) assume that our original nonlinearity l (see (6)) has the form For every x ∈ X F consider a weak solution x ∈ U of the problem Define the map H assigning to x ∈ X F a weak solution x ∈ U of the above linear problem. Let us put X F = H(X F ).

Theorem 2 (Main theorem). Under (As) and (HΓ ) there exists x ∈ X F such that
and x is a weak solution to (6).
Remark 4 We would like to stress that the result obtained in the paper is of global type-for each given, but fixed T < ∞. However assumption (As) shows us that there is a certain relation between the type of nonlinearity and the quantity of , have to be bounded in some ball of U ∩ L 1 (0, T ; Lq (Ω)), q > 1. Moreover, we do not admit the case T = ∞. The known results concerning existence of global solutions to (1) (without uniqueness and damping term) assume that nonlinearity l is locally Lipschitz and sublinear at infinity or of polynomial growth |s| p−1 s ( p such that H 1 (Ω) ⊂ L p+1 (Ω)) with additional boundedness for initial data (see e.g. [9]). We do not assume that l is locally Lipschitz from H 1 (Ω) → L p+1 (Ω), as x → l(x) is only continuous in R. Instead, we assume the special structure of l (see (12))-a combination of monotone function and continuous ones. Moreover the assumptions: F x (x), G x (x) are bounded in some ball imply growth conditions for l. Some type of boundedness for initial data are hidden in the definitions of the set X F . We would like to stress that l is only continuous in x, usually l is at least of C 1 (see e.g. [4,5,9]). Moreover we are looking for solutions in the bounded set X F , so all solutions are bounded in the norm of C([0, T ]; H 1 (Ω)). This is a case which is often assumed to get global solution (see [4,9] and references therein).
We ilustrate the above theorem (its assumptions) by an example.
We have no uniqueness of solutions to (1), but a kind of continuous dependence on initial data is still possible.
Let F * be the Fenchel conjugate of F. Define a dual to J F functional where ·, · is the inner product in R n and ν = (ν 1 , . . . , ν n ) is the unit outward normal to Γ and x is the minimizer of J (see Theorem 2). Now we can formulate theorem which gives us additional informations on solutions to (1) important in classical mechanics. This theorem is absolutely new for problem (1).
Theorem 4 (Variational principle and duality result) Assume (As). Let x ∈ X F be such that q(t, y) = ∇x(t, y), and The proofs of theorems are given in Sects. 3 and 4. They consist of several steps. First we prove Propositions 1, 2 and Corollary 1, i.e. we solve problem (5). Next we prove Theorem 2 (Main theorem).

Proof of Existence for Problem (5)
As mentioned before, the main difficulty of the problem under study is the fact that the Neumann problem does not satisfy Lopatinski condition and therefore, the map from the boundary data in L 2 (Σ) into finite energy space is not bounded (unless dimension of Ω is equal to one). In order to cope with this problem in [4,9] a regularizing term-strongly monotone dissipation is introduced, whose effect is to 'force' the Lopatinski condition. We follow, in quite, different way. We convert the problem (5) into variational one and this allow us to omit the meaning of the dissipation term. However, the price we pay for that is the form of nonlinearity for the boundary source term equals h − x. We start with the proof of Proposition 1 Proof of Proposition 1 It is a simple consequence of the equivalence of the following two relations in (0, T ) × Γ : H (t, y, x(t, y)), y, x(t, y)) + H * (t, y, −∂ ν x(t, y) − g(t, y, x t (t, y)) + h(t, y, x(t, y))) = x(t, y) (−∂ ν x(t, y) − g(t, y, x t (t, y)) + h(t, y, x(t, y))) (19) and the inequality H (t, y, x(t, y)) + H * (t, y, −∂ ν x(t, y) − g(t, y, x t (t, y)) + h(t, y, x(t, y))) ≥ x(t, y) (−∂ ν x(t, y) − g(t, y, x t (t, y)) + h(t, y, x(t, y))) (20) which holds for all x(t, ·) ∈ U t , t ∈ [0, T ].
Let us notice that the sets U t are nonempty. We will study the functional J t Γ (x) on U t , t ∈ [0, T ] and we prove under hypothesis (HΓ ) that there exists a minimum to the functional J t Γ , i.e. the proof of Proposition 2. (t, ·)) ≥ 0 and weakly lower semicontinuous in U t , t ∈ [0, T ]. Moreover, let us observe, by (HΓ ), J t Γ (x(t, ·)) → ∞ when x(t, ·) H 1 (Γ ) → ∞. Really, it is enough to notice that for t ∈ [0, T ] and x ∈ U t J t Γ (x(t, ·)) = y, x(t, y))) 2 dy
It is clear that J t Γ (x(t, ·)) = J t Γ D (−x(t, ·)) for all x(t, ·) ∈ U t and so By the duality theory for convex functionals (see [11,18] Proof of Corollary 1 This is a direct consequence of two above Propositions 1 and 2.

Proof of Existence for Problem (6)
Consider equation (with x Γ ∈ H 1 (Σ) being a solution to (5)) x tt (t, y) − x(t, y) + F x (t, y, x(t, y)) − G x (t, y, x(t, y)) = 0, in (0, T ) × Ω, x(t, y) = x Γ (t, y), (t, y) ∈ Σ, and corresponding to (23) the functional y, x(t, y)))dydt defined on U. Notice that (23) is not the Euler-Lagrange equation for the action functional J F . The dual functional to (24), for any given w ∈ U , now reads y, w(t, y)))) dydt We prove the following where v is the weak solution of the problem y, x(t, y)) + G x (t, y, x(t, y) Proof Fix arbitrary x ∈ X F . Since x ∈ U and by the assumptions on F and G, see (As), it follows that F x (·, ·, x(·, ·)), G x (·, ·, x(·, ·)) ∈ L 1 (0, T ; L2 (Ω)). Hence by Theorem 1 and (10) there exists a unique solution v ∈ U of our problem for the Eq. (26) satisfying Taking into account the definition of the set X F we get the following estimation, independent of and since for some Hence, putting infer the first assertion of the lemma. Further by (10) and (11) and we get the last assertion of the lemma.

Proposition 3
For every x ∈ X F the weak solution x of the problem y, x (t, y)) + G x (t, y, x(t, y)), belongs to X F .
Further by (10) we get and hence .
Remark 5 Proposition 3 apparently may suggest that it is more convenient to apply a suitable fixed point theorem in order to get the existence of weak solutions to problem (1). Indeed, the above mentioned proposition states that the map H assigning to x ∈ X F a weak solution x ∈ X F of (27) has the property H(X F ) ⊂ X F . However this is only starting point in fixed point theory. In order to proceed with the so called topological method we must prove that the map H and set H(X F ) possess suitable properties. However the assumptions (As) do not imply directly (if at all) neither that H is contraction nor that H(X F ) is convex or relatively compact. We have chosen variational approach which ensures not only the existence of solutions but also certain variational properties of solutions which are absolutely new in that case.
Proof It is clear, by Lemma 1, that {x j tt } from X tt has a subsequence converging weakly in L2 (0, T ; H −1 (Ω)). By the same lemma corresponding sequence {x j } is also weakly (up to some subsequence) converging in L2 (0, T ; H 1 (Ω)) to some element x. By the definition of X F and uniqueness of weak solutions of (27) we infer that {x j tt } has a subsequence weakly convergent to x tt . The last, the fact that {x j t } is bounded in L2 (0, T ; L 2 (Ω)) and a known theorem imply that for this subsequent {x j t } is convergent strongly in L2 (0, T ; L 2 (Ω)) to x t .
We observe that functional J F is well defined on X F . Moreover we have Proof Note that by the construction of X F it is a ball in U, we have also that U ⊂ H 1 ((0, T )× Ω). It is easy to see that X F is weakly compact in H 1 ((0, T ) × Ω) and that the functionals y, x(t, y))dydt, x →  H 1 ((0, T ) × Ω) is strongly convergent in L 2 ((0, T ) × Ω)). Hence we infer the assertions of the Lemma.

Lemma 4 The functional J F attains its maximum in X
Proof By definition of the set X F and Lemma 3 we see that the functional J F is bounded in X F . We denote by {x j } a maximizing sequence for J F in X F . This sequence has a subsequence which we denote again by {x j } converging weakly in L2 (0, T ; H 1 (Ω)) and strongly in L2 (0, T ; L2 (Ω)), hence also strongly in L2 ((0, T ) × Ω; R) to a certain element x ∈ U . Moreover {x j } is also convergent almost everywhere. Thus by the construction of the set X F and uniqueness of weak solutions of (27) we observe that x ∈ X F . Hence To consider properly the dual action functional let us put and W 1F y = q ∈ L2 0, T ; L2 (Ω) n : there exists p ∈ W 1F t such that p t (·, ·) − div q(·, ·) ∈ L2 (0, T ; L2 (Ω)) .
and define a set on which a dual functional will be considered. Definition of X Fd We say that an element ( p, q) ∈ W 1F t × W 1F y belongs to X Fd provided that there exists x ∈ X F such that for a.e. t ∈ (0, T ) with q (t, ·) = ∇x (t, ·) .
We observe that neither X F nor X Fd is a linear space. Thus even standard calculations using convexity arguments are rather difficult. What helps us is a special structure of the sets X F and X Fd which despite their nonlinearity makes these calculations possible. Now we may state the main result of the paper which is the following existence theorem.

Variational Principle: The Case l = F x − G x
We state the necessary conditions. We observe that due to the construction of the set X F and Lemma 2 it follows that a maximizing sequence in X F for the functional J F may be assumed to be weakly convergent in L2 (0, T ; H 1 (Ω)) and strongly in L2 (0, T ; L2 (Ω)).
By the definitions of J F , J F D , relations (37), (36) and the Fenchel-Young inequality it follows that , p (T, ·)) L 2 (Ω) Therefore we get that Letx ∈ X F denote the element corresponding to p accordingly to the definition of the set X Fd i.e.x t = p . It is clear that then J F (x) ≤ J F (x). We observe that by the Fenchel-Young equality we have . y, x (t, y)) dydt y, x (t, y))) dydt Therefore by (37), (36) and standard convexity arguments Hence by the Fenchel-Young transform we obtain that p (t, y) = x t (t, y) i.e. (33).

Variational Principle: The Case l = −F x − G x
A similar theorem is true for the problem and corresponding to it functional defined on U with same hypotheses (As) and the set X F . Really, Lemmas 1-4 are still valid as sign of F does not change their proofs, also the proof of the above theorem does not change. Hence we get for (39) the following theorem.

Proof of Stability: Theorem 3
In this section we prove some stability results in our nonlinear case i.e. Theorem 3. To this effect let us assume we are given a sequence {x 0 n } in H 1 (Γ ) converging to x 0 in H 1 (Γ ).
Let {x Γ n } be a sequence of solutions to (5) corresponding to {x 0 n }. Then by Proposition 2 and its proof we know that J t Γ (x Γ n ) = 0, t ∈ [0, T ], n = 1, 2, . . . and that J t Γ (x) ≥ 0 for all x ∈ U t , t ∈ [0, T ] (moreover J t Γ is weakly lower semicontinuous in U t , t ∈ [0, T ]). Therefore the sequence {x Γ n } is bounded in H 1 (Γ ) and hence it has a subsequence (which we shall denote again by {x Γ n }) converging weakly in H 1 (Γ ) to some x Γ ∈ H 1 (Γ ) and J t Γ (x Γ ) = 0 on U t , t ∈ [0, T ]. The last implies that x Γ is a solution to (5) corresponding to x 0 .
Next let us assume that {x 0 n }, {x 1 n } are given sequences in H 1 (Γ ), L2 (Ω), respectively and such that (Ω) and x 0 n (y) = x Γ n (0, y) on Γ , n = 1, 2, . . . . Thus we assume hypotheses (HΓ ), (As) to be considered and satisfied. In consequence, all the assertions of Theorem 4 are true for all n with estimations independent on n. Because, we have nonlinear problem we cannot expect the same type of continuous dependence as in [13] for linear case. First note that for each x 0 n , x 1 n there exists a solution x n ∈ X F ⊂ U to (1) determined by x 0 n and x 1 n , n = 1, 2, . . . . Therefore, choosing suitable subsequence, let x ∈ X F be a weak limit in H 1 ((0, T ) × Ω) and strong in L2 ((0, T ) × Ω) of {x n } . We know also that for all n, J F (x n ) = J F D p n , q n where p n , q n denote the pair of sequences corresponding to {x n } and satisfying x nt (t, ·) = p n (t, ·) , ∇x n (t, ·) = q n (t, ·) , − p nt (t, ·) + div q n (t, ·) − l (t, y, x n (t, ·)) = 0.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.