%% R.E. Showalter: Chapter 6 \chapter{Second Order Evolution Equations} \section{Introduction} % 1 \setcounter{equation}{0} \setcounter{theorem}{0} We shall find well-posed problems for evolution equations which contain the second order time derivative of the solution. These arise, for example, when we attempt to use the techniques of the preceding chapters to solve a Cauchy problem for the {\it wave equation\/} \begin{equation}\label{eq611} \partial_t^2 u(x,t) - \Delta_n u(x,t) = F(x,t)\ . \end{equation} The corresponding abstract problem will contain the second order evolution equation \begin{equation}\label{eq612} u''(t) + \A u(t) = f(t)\ , \end{equation} where $\A$ is an operator which contains $-\Delta_n$ in some sense. Wave equations with damping or friction occur in practice, e.g., the {\it telegraphists equation\/} $$\partial_t^2 u(x,t) + R\cdot\partial_t u(x,t) -\Delta_n u(x,t) = F(x,t)\ ,$$ so we shall add terms to \eqn{612} of the form $\B u'(t)$. Finally, certain models in fluid mechanics lead to equations, for example, \begin{equation}\label{eq613} \partial_t^2 (\Delta_n u(x,t)) +\partial_n^2 u(x,t)=0\ ,\qquad x=(x_1,\ldots,x_n)\ , \end{equation} which contain spatial derivatives in the terms with highest (= second) order time derivatives. These motivate us to consider abstract evolution equations of the form \begin{equation}\label{eq614} \C u''(t) +\B u'(t) +\A u(t) = f(t)\ ,\qquad t>0\ . \end{equation} We consider in Section 2 equations of the form \eqn{614} in which $\C$ is invertible; this situation is similar to that of Section~V.2, so we call \eqn{614} a {\it regular equation\/} then. The equation \eqn{613} is known as Sobolev's equation, so we call \eqn{614} a {\it Sobolev equation\/} when $\C$ is invertible and both $\C^{-1}\B$ and $\C^{-1}\A$ are continuous. This situation is studied in Section~3 and is the analogue of (first-order) pseudoparabolic problems. Section~4 will be concerned with \eqn{614} when $\C$ is {\it degenerate\/} in the sense of Section~V.4. Such equations arise, for example from a system described by appropriately coupled wave and heat equations $$\begin{array}{ll} \partial_t^2 u(x,t)-\Delta_n u(x,t) =0\ ,&\qquad x\in G_1\ ,\\ \noalign{\vskip6pt} \partial_tu (x,t)-\Delta_n u(x,t)=0\ ,&\qquad x\in G_2\ . \end{array}$$ Here the operator $\C$ is multiplication by the characteristic function of $G_1$ and $\B$ is multiplication by the characteristic function of $G_2$. $G_1$ and $G_2$ are disjoint open sets whose closures intersect in an $(n-1)$-dimensional manifold or {\it interface\/}. Additional examples will be given in Section~5. \section{Regular Equations} % 2 \setcounter{equation}{0} Let $V$ and $W$ be Hilbert spaces with $V$ a dense subspace of $W$ for which the injection is continuous. Thus, we identify $W'\le V'$ by duality. Let $\A \in \L(V,V')$ and $\C\in \L(W,W')$ be given. Suppose $D(B) \le V$ and $B:D(B)\to V'$ is linear. If $u_0 \in V$, $u_1\in W$ and $f\in C((0,\infty),W')$ are given, we consider the problem of finding $u\in C([0,\infty),V)\cap C^1((0,\infty), V)\cap C^1 ([0,\infty),W) \cap C^2 ((0,\infty),W)$ such that $u(0)=u_0$, $u'(0)= u_1$, and \begin{equation}\label{eq621} \C u''(t) + B u'(t) + \A u(t) = f(t) \end{equation} for all $t>0$. Note that for any such solution of \eqn{621} we have $u'(t)\in D(B)$ and $Bu'(t)+\A u(t)\in W'$ for all $t>0$. We shall solve \eqn{621} by reducing it to a first order equation on a product space and then applying the results of Section~V.2. The idea is to write \eqn{621} in the form $$ \left(\begin{array}{cc} \A &0\\ 0&\C \end{array}\right) \left(\begin{array}{c} u\\ u'\end{array}\right)' + \left(\begin{array}{cc} 0&-\A\\ \A&B\end{array}\right) \left(\begin{array}{c} u\\ u'\end{array}\right) = \left(\begin{array}{c} 0\\ f(t)\end{array}\right)\ .$$ Define $V_m=V\times W$, the product Hilbert space with scalar-product given by $$([x_1,x_2],[y_1,y_2])_{V_m} = (x_1,y_1)_V + (x_2,y_2)_W\ ,\qquad [x_1,x_2],[y_1,y_2]\in V\times W\ .$$ We have then $V'_m = V'\times W'$, and we define $\M\in \L(V_m,V'_m)$ by $$\M([x_1,x_2]) = [\A x_1,\C x_2]\ ,\qquad [x_1,x_2]\in V_m\ .$$ Define $D= \{[x_1,x_2]\in V\times D(B) :\A x_1+Bx_2\in W'\}$ and $L\in L(D,V'_m)$ by $$L([x_1,x_2]) = [-\A x_2,\A x_1 +Bx_2]\ ,\qquad [x_1,x_2]\in D\ .$$ If $u(\cdot)$ is a solution of \eqn{621}, then the function defined by $w(t) = [u(t),u'(t)]$, $t\ge 0$, satisfies the following: $w\in C([0,\infty),V_m)\cap C^1 ((0,\infty),V_m)$, $w(0) = [u_0,u_1]\in V_m$, and \begin{equation}\label{eq622} \M w'(t) + Lw(t) = [0,f(t)]\ ,\qquad t>0\ . \end{equation} This is precisely the situation of Section V.2, so we need only to find conditions on the data in \eqn{621} so that Theorems~\ref{thm5-2A} or \ref{thm5-2B} of Chapter V are applicable. This leads to the following. \begin{theorem}\label{thm6-2A} Let $V$ and $W$ be Hilbert spaces with $V$ dense and continuously imbedded in $W$. Assume $\A \in \L(V,V')$ and $\C \in \L(W,W')$ are the Riesz maps of $V$ and $W$, respectively, and let $B$ be linear from the subspace $D(B)$ of $V$ into $V'$. Assume that $B$ is monotone and that $\A + B + \C:D(B)\to V'$ is surjective. Then for every $f\in C^1([0,\infty),W')$ and $u_0 \in V$, $u_1\in D(B)$ with $\A u_0 + Bu_1 \in W'$, there exists a unique solution $u(t)$ of \eqn{621} (on $t\ge 0$) with $u(0)=u_0$ and $u'(0)=u_1$. \end{theorem} \proof Since $\A$ and $\C$ are Riesz maps of their corresponding spaces, we have \begin{eqnarray*} \M([x_1,x_2])([y_1,y_2]) &=& \A x_1(y_1) +\C x_2(y_2)\\ &=& (x_1,y_1)_V + (x_2,y_2)_W\\ &=& ([x_1,x_2],[y_1,y_2])_{V_m}\ ,\qquad [x_1,x_2],[y_1,y_2] \in V_m\ , \end{eqnarray*} so $\M$ is the Riesz map of $V_m$. Also we have for $[x_1,x_2]\in D$ $$L([x_1,x_2])([y_1,y_2]) = -\A x_2 (y_1)+ (\A x_1 + Bx_2)(y_2) \ ,\qquad [y_1,y_2] \in V_m\ ,$$ hence, $L([x_1,x_2])([x_1,x_2]) = -\overline{\A x_1(x_2)} + \A x_1(x_2) + Bx_2(x_2)$ since $\A$ is symmetric. From this we obtain $$\Re L ([x_1,x_2]) ([x_1,x_2]) = \Re B x_2 (x_2)\ ,\qquad [x_1,x_2]\in D\ ,$$ so $B$ being monotone implies $L$ is monotone. Finally, if $f_1\in V'$ and $f_2\in W'$, then we can find $x_2\in D(B)$ such that $(\A + B+\C)x_2 = f_2-f_1$. Setting $x_1=x_2 + \A^{-1} f_1\in V$, we have a pair $[x_1,x_2]\in D$ for which $(\M +L)[x_1,x_2] = [f_1,f_2]$. (Note that $\A x_1 + Bx_2 = f_2-\C x_2 \in W'$ as required.) Thus $\M +L$ is a surjection of $D$ onto $V'_m$. Theorem \ref{thm5-2A} of Chapter~V asserts the existence of a solution $w(t) = [u(t),v(t)]$ of \eqn{622}. Since $\A$ is a norm-preserving isomorphism, $v(t) = u'(t)$ and the result follows. A special case of Theorem \ref{thm6-2A} that occurs frequently in applications is that $D(B) =V$ and $B=\B \in \L(V,V')$. Then one needs only to verify that $\B$ is monotone, for then $\A+\B+\C$ is $V$-coercive, hence surjective. Furthermore, in this case we may define $\L\in \L(V_\ell,V'_\ell)$ and $V_\ell = V\times V$ by $$\L([x_1,x_2])([y_1,y_2]) = -\A x_2 (y_1) +(\A x_1 +\B x_2) (y_2)\ ,\qquad [x_1,x_2],[y_1,y_2]\in V_\ell\ .$$ Thus, Theorem \ref{thm5-2B} of Chapter V applies if we can show that $\L(\cdot)(\cdot)$ is $V_\ell$-elliptic. Of course we need only to verify that $(\lambda \M+\L) (\cdot)(\cdot)$ is $V_\ell$-elliptic for some $\lambda>0$ (Exercise V.2.3), and this leads us to the following. \begin{theorem}\label{thm6-2B} Let $\A$ and $\C$ be the Riesz maps of the Hilbert spaces $V$ and $W$, respectively, where $V$ is dense and continuously imbedded in $W$. Let $\B \in \L(V,V')$ and assume $\B +\lambda\C$ is $V$-elliptic for some $\lambda>0$. Then for every H\"older continuous $f:[0,\infty)\to W'$, $u_0\in V$ and $u_1\in W$, there is a unique solution $u(t)$ of \eqn{621} on $t>0$ with $u(0)=u_0$ and $u'(0)=u_1$. \end{theorem} Theorem \ref{thm6-2B} applies to evolution equations of second order which are parabolic, i.e., those which can be solved for more general data $u_0$, $u_1$ and $f(\cdot)$, and whose solutions are smooth for all $t>0$. Such problems occur when energy is strongly dissipated; we give examples below. The situation in which energy is conserved is described in the following result. We leave its proof as an exercise, as it is a direct consequence of either Theorem \ref{thm6-2B} above or Section IV.5. \begin{theorem}\label{thm6-2C} In addition to the hypotheses of Theorem \ref{thm6-2A}, assume that $\Re Bx(x)=0$ for all $x\in D(B)$ and that both $\A+B+\C$ and $\A-B+\C$ are surjections of $D(B)$ onto $V'$. Then for every $f\in C^1(\RR,W')$ and $u_0\in V$, $u_1\in D(B)$ with $\A u_0+ Bu_1\in W'$, there exists a unique solution of \eqn{621} on $\RR$ with $u(0)=u_0$ and $u'(0)=u_1$. \end{theorem} We shall describe how Theorems \ref{thm6-2A} and \ref{thm6-2C} apply to an {\it abstract wave equation\/}. Examples will be given afterward. Assume we are given Hilbert spaces $V\le H$, and $B$, and a linear surjection $\gamma :V\to B$ with kernel $V_0$ such that $\gamma$ factors into an isomorphism of $V/V_0$ onto $B$, the injection $V\hookrightarrow H$ is continuous, $V_0$ is dense in $H$, and $H$ is identified with its dual $H'$ by the Riesz map. We thereby obtain continuous injections $V_0\hookrightarrow H\hookrightarrow V'_0$ and $V\hookrightarrow H\hookrightarrow V'$. Let $a_1:V\times V\to \KK$ and $a_2:B\times B\to \KK$ be continuous symmetric sesquilinear forms and define $a:V\times V\to\KK$ by \begin{equation}\label{eq623} a(u,v) = a_1(u,v) + a_2(\gamma(u),\gamma(v))\ ,\qquad u,v\in V\ . \end{equation} Assume $a(\cdot,\cdot)$ is $V$-elliptic. Then $a(\cdot,\cdot)$ is a scalar-product on $V$ which gives an equivalent norm on $V$, so we hereafter consider $V$ with this scalar-product, i.e., $(u,v)_V\equiv a(u,v)$ for $u,v\in V$. The form \eqn{623} will be used to prescribe an abstract boundary value problem as in Section~III.3. Thus, we define $A:V\to V'_0$ by $$Au (v) = a_1(u,v)\ ,\qquad u\in V\ ,\ v\in V_0$$ and $D_0 = \{u\in V:Au\in H\}$. Then Theorem III.\ref{thm3-2C} gives the abstract boundary operator $\partial_1\in L(D_0,B')$ for which $$a_1(u,v)- (Au,v)_H = \partial_1u(\gamma v)\ ,\qquad u\in D_0\ ,\ v\in V\ .$$ We define $D=\{ u\in V:\A u\in \H\}$, where $\A$ is the Riesz map of $V$ given by $$\A u (v)= a(u,v)\ ,\qquad u,v\in V\ ,$$ and $\A_2 :B\to B'$ is given by $$\A_2 \varphi (\psi) = a_2 (\varphi,\psi)\ ,\qquad\varphi,\psi\in B\ .$$ Then, we recall from Corollary III.\ref{cor3-3B} that $u\in D$ if and only if $u\in D_0$ and $\partial_1 u+\A_2 (\gamma u) =0$. Let $(\cdot,\cdot)_W$ be a scalar-product on $H$ whose corresponding norm is equivalent to that of $(\cdot,\cdot)_H$, and let $W$ denote the Hilbert space consisting of $H$ with the scalar-product $(\cdot,\cdot)_W$. Then the Riesz map $\C$ of $W$ satisfies $\C\in\L(H)$ (and $\C^{-1}\in \L(H)$). Suppose we are also given an operator $\B\in\L(V,H)$ which is monotone (since $H\le V'$). \begin{theorem}\label{thm6-2D} Assume we are given the Hilbert spaces $V$, $H$, $B$, $V_0$, $W$ and linear operators $\gamma$, $\partial_1$, $\A_2$, $A$, $\A$, $\B$ and $\C$ as above. Then for every $f\in C^1([0,\infty),H)$, $u_0\in D$ and $u_1\in V$, there is a unique solution $u(\cdot)$ of \eqn{621}, and it satisfies \begin{equation}\label{eq624} \left. \begin{array}{ll} \C u''(t) + \B u'(t) + Au(t) =f(t)\ ,&\qquad t\ge0\ ,\\ \noalign{\vskip6pt} u(t)\in V\ ,\quad \partial_1 u(t)+\A_2 \gamma(u(t))=0\ ,&\qquad t\ge 0\ ,\\ \noalign{\vskip6pt} u(0)=u_0\ ,\quad u'(0)=u_1\ ,&\end{array}\right\} \end{equation} \end{theorem} \proof Since $\A +\B +\C \in\L(V,V')$ is $V$-elliptic, it is surjective so Theorem \ref{thm6-2A} (with $B=\B$) asserts the existence of a unique solution. Also, since each of the terms $\C u''(t)$, $\B u'(t)$ and $f(t)$ of the equation \eqn{621} are in $H$, it follows that $\A u(t)\in H$ and, hence, $u(t)\in D$. This gives the middle line in \eqn{624}. \qed In each of our examples below, the first line in \eqn{624} will imply an abstract wave equation, possibly with damping, and the second line will imply boundary conditions. \subsection{} % 2.1 Let $G$ be open in $\RR^n$ and take $H=L^2(G)$. Let $\rho \in L^\infty (G)$ satisfy $\rho (x)\ge c>0$ for $x\in G$, and define $$(u,v)_W\equiv \int_G \rho (x) u(x)\overline{v(x)}\,dx\ ,\qquad u,v\in H\ .$$ Then $\C$ is just multiplication by $\rho (\cdot)$. Suppose further that $\partial G$ is a $C^1$ manifold and $\Gamma$ is a closed subset of $\partial G$. We define $V= \{v\in H^1(G) : \gamma_0 (v) (s)=0$, a.e., $s\in\Gamma\}$, $\gamma=\gamma_0|_V$ and, hence, $V_0 = H_0^1(G)$ and $B$ is the range of $\gamma$. Note that $B\hookrightarrow L^2(\partial G\sim \Gamma)\hookrightarrow B'$. We define $$a_1 (u,v) = \int_G \nabla u\cdot\nabla \bar v\,dx\ ,\qquad u,v\in V\ ,$$ and it follows that $A= -\Delta_n$ and $\partial_1$ is the normal derivative $${\partial u\over\partial \nu} = \nabla u\cdot\nu$$ on $\partial G$. Let $\alpha \in L^\infty (\partial G)$ satisfy $\alpha (s) \ge0$, a.e. $s\in \partial G$, and define $$a_2 (\varphi,\psi) = \int_{\partial G\sim\Gamma} \alpha (s)\varphi (s) \overline{\psi (s)}\,ds\ ,\qquad \varphi,\psi \in B\ .$$ Then $\A_2$ is multiplication by $\alpha(\cdot)$. Assume that for each $t\in [0,T]$ we are given $F(\cdot,t)\in L^2(G)$, that $\partial_t F(x,t)$ is continuous in $t$ for almost every $x\in G$, and $|\partial_tF(x,t)|\le g(x)$ for some $g\in L^2(G)$. It follows that the map $t\mapsto F(\cdot,t)\equiv f(t)$ belongs to $C^1([0,T],L^2(G))$. Finally, let $U_0 (\cdot)\in D$ (see below) and $U_1(\cdot)\in V$ be given. Then, if $u(\cdot)$ denotes the solution of \eqn{624} it follows from Theorem IV.\ref{thm4-7A} that we can construct a function $U\in L^2(G\times [0,T])$ such that $U(\cdot,t)=u(t)$ in $L^2(G)$ for each $t\in [0,T]$ and this function satisfies the partial differential equation \begin{equation}\label{eq625} \rho (x)\partial_t^2 U(x,t) -\Delta_n U(x,t) = F(x,t)\ ,\qquad x\in G\ ,\ 0\le t\le T \end{equation} and the initial conditions $$U(x,0) = U_0 (x)\ ,\quad \partial_t U(x,0) = U_1(x)\ ,\quad \hbox{a.e. }\ x\in G\ .$$ Finally, from the inclusion $u(t)\in D$ we obtain the boundary conditions for $t\ge0$ \begin{equation}\label{eq626} \left. \begin{array}{ll} U(s,t) = 0\ ,&\qquad \hbox{a.e. }\ s\in \Gamma\ \hbox{ and}\\ \noalign{\vskip6pt} \ds {\partial U(s,t)\over \partial\nu} +\alpha (s)U(s,t)=0\ , &\qquad\hbox{a.e. }\ s\in \partial G\sim\Gamma\ .\end{array}\right\} \end{equation} The first equation in \eqn{626} is the boundary condition of {\it first type\/}. The second is the boundary condition of {\it second type\/} where $\alpha (s)=0$ and of {\it third type\/} where $\alpha (s)>0$. (Note that $U_0$ necessarily satisfies the conditions of \eqn{626} with $t=0$ and that $U_1$ satisfies the first condition in \eqn{626}. If $F(\cdot,t)$ is given as above but for each $t\in [-T,T]$, then Theorem \ref{thm6-2C} (and Theorem III.\ref{thm3-7A}) give a solution of \eqn{625} on $G\times [-T,T]$. \subsection{} % 2.2. In addition to all the data above, suppose we are given $R(\cdot) \in L^\infty (G)$ and a vector field $\mu (x)=(\mu_1(x),\ldots, \mu_n(x))$, $x\in G$, with each $\mu_j \in C^1(\bar G)$. We define $\B \in \L(V,H)$ (where $V\le H^1(G)$ and $H=L^2(G))$ by \begin{equation}\label{eq627} \B u (v) = \int_G \left( R(x)u(x)+ {\partial u(x)\over\partial \mu}\right) \overline{v(x)}\,dx \end{equation} the indicated directional derivative being given by $${\partial u(x)\over\partial\mu} \equiv \sum_{j=1}^n \partial_j u(x) \mu_j (x)\ .$$ From the Divergence Theorem we obtain $$2\Re \int_G {\partial u(x)\over\partial \mu}\,\overline{u(x)}\,dx + \int_G \biggl( \sum_{j=1}^n \partial_j \mu_j(x)\biggr) |u(x)|^2\,dx = \int_{\partial G} (\mu\cdot\nu)|u(x)|^2\,ds\ ,$$ where $\mu\cdot\nu = \sum_{j=1}^n \mu_j (s)\nu_j (s)$ is the indicated euclidean scalar-product. Thus, $\B$ is monotone if $$\begin{array}{rcll} \ds -\left(\tfrac12\right)\sum_{j=1}^n\partial_j\mu_j (x)+\Re \{R(x)\}&\ge&0\ , &\qquad x\in G\\ \noalign{\vskip6pt} \mu (s)\cdot\nu (s)&\ge&0\ ,&\qquad s\in \partial G\sim\Gamma\ .\end{array}$$ The first equation represents friction or energy dissipation distributed throughout $G$ and the second is friction distributed over $\partial G$. Note that these are determined by the divergence of $\mu$ and the normal component of $\mu$, respectively. If $u(\cdot)$ is a solution of \eqn{624} and the corresponding $U(\cdot,\cdot)$ is obtained as before from Theorem IV.\ref{thm4-7A}, then $U(\cdot,\cdot)$ is a generalized solution of the initial-boundary value problem $$\cases{ \ds \rho (x)\partial_t^2 U(x,t) + R(x)\partial_t U(x,t) + \partial_t {\partial U(x,t)\over\partial \mu} - \Delta_n U(x,t) = F(x,t)\ ,&\cr \noalign{\vskip6pt} \hfill x\in G\ ,\ t\ge 0&\cr \noalign{\vskip6pt} U(s,t) =0\ ,\qquad \hbox{a.e. }\ s\in \Gamma\ ,&\cr \noalign{\vskip6pt} \ds {\partial U (s,t)\over\partial \nu} +\alpha (s) U(s,t)=0\ ,\qquad\hbox{a.e. } \ s\in \partial G\sim\Gamma&\cr \noalign{\vskip6pt} U(x,0) = U_0(x)\ ,\quad \partial_t U(x,0) = U_1(x)&\cr}$$ One could similarly solve problems with the fourth boundary condition, oblique derivatives, transition conditions on an interface, etc., as in Section III.4. We leave the details as exercises. We now describe how Theorem \ref{thm6-2B} applies to an {\it abstract viscoelasticity equation\/}. \begin{theorem}\label{thm6-2E} Assume we are given the Hilbert spaces $V$, $H$, $B$, $V_0$, $W$ and linear operators $\gamma$, $\partial_1$, $\A_2$, $A$, $\A$, $\B$ and $\C$ as in Theorem \ref{thm6-2D}. Then for every $f:[0,\infty)\to H$ which is H\"older continuous, $u_0\in V$ and $u_1\in H$, there is a unique solution $u(t)$ of \eqn{621} with $B= \B+ \varep \A$ and $\varep>0$. This solution satisfies \begin{equation}\label{eq628} \left.\begin{array}{ll} \C u''(t) + (\B +\varep A) u'(t) +Au(t)=f(t)\ ,&\qquad t>0\ ,\\ \noalign{\vskip6pt} u(t) \in V\ ,&\qquad t\ge0\ ,\\ \noalign{\vskip6pt} \partial_1 (\varep u'(t)+u(t)) +\A_2\gamma (\varep u'(t) +u(t))=0\ , &\qquad t>0\ ,\\ \noalign{\vskip6pt} u(0)=u_0\ ,\quad u'(0)=u_1\ .&\end{array}\right\} \end{equation} \end{theorem} \proof This follows immediately from $$\Re Bx(x) \ge \varep \A x(x) = \varep \|x\|_V^2\ ,\qquad x\in V\ ,$$ (since $\B$ is monotone) and the observation that $\varep u'(t)+u(t)\in D$ for $t>0$. \subsection{} %2.3 Let all spaces and operators be chosen just as in Section 2.1 above. Suppose $U_0\in V$, $U_1\in H$ and $f(t) = F(t,\cdot)$, $t\ge0$, where $F(\cdot,\cdot)$ is given as in Theorem IV.\ref{thm4-7C}. Then we obtain a generalized solution of the initial-boundary value problem \begin{equation}\label{eq629} \left.\begin{array}{l} \rho (x)\partial_t^2U(x,t)-\varep \partial_t\Delta_n U(x,t) - \Delta_n U(x,t) = F(x,t)\ ,\\ \hfill\hbox{a.e. } x\in G\ ,\ t>0\ ,\\ \noalign{\vskip6pt} U(s,t)=0\ ,\qquad \hbox{a.e. }s\in\Gamma\ ,\ t\ge0\ ,\\ \noalign{\vskip6pt} \ds {\partial\over\partial\nu} (\varep \partial_t U(s,t) +U(s,t)) +\alpha (s) (\varep \partial_t U(s,t) +U(s,t)) =0\ ,\\ \hfill\hbox{a.e. } s\in\partial G\sim\Gamma\ ,\ t>0\ ,\\ \noalign{\vskip6pt} U(x,0)= U_0(x) \ ,\quad \partial_t U(x,0)= U_1(x)\ ,\qquad x\in G\ . \end{array}\right\} \end{equation} In certain applications the coefficient $\varep>0$ corresponds to {\it viscosity\/} in the model and it distinguishes the preceding parabolic problem from the corresponding hyperbolic problem in Section~2.1. Problems with viscosity result in very strong damping effects on solutions. Dissipation terms of lower order like \eqn{627} could easily be added to the system \eqn{629}, and other types of boundary conditions could be obtained. \section{Sobolev Equations} % 3 \setcounter{equation}{0} We shall give sufficient conditions for a certain type of evolution equation to have either a weak solution or a strong solution, a situation similar to that for pseudoparabolic equations. The problems we consider here have the strongest operator as the coefficient of the term in the equation with the second order derivative. \begin{theorem}\label{thm6-3A} Let $V$ be a Hilbert space and $\A, \B, \C \in \L(V,V')$. Assume that the sesquilinear form corresponding to $\C$ is $V$-elliptic. Then for every $u_0,u_1\in V$ and $f\in C(\RR,V)$ there is a unique $u\in C^2(\RR,V)$ such that \begin{equation}\label{eq631} \C u'' (t) +\B u'(t) +\A u(t) = f(t)\ ,\qquad t\in \RR\ , \end{equation} and $u(0)= u_0$, $u'(0)=u_1$. \end{theorem} \proof The change of variable $v(t)\equiv e^{-\lambda t} u(t)$ gives an equivalent problem with $\A$ replaced by $\A+\lambda\B +\lambda^2\C$, and this last operator is $V$-coercive if $\lambda$ is chosen sufficiently large. Hence, we may assume $\A$ is $V$-elliptic. If we define $\M$ and $\L$ as in Section~2.2, then $\M$ is $V\times V \equiv V_m$-elliptic, and Theorem V.\ref{thm5-3A} then applies to give a solution of \eqn{622}. The desired result then follows. A solution $u\in C^2(\RR,V)$ of \eqn{631} is called a {\it weak solution\/}. If we are given a Hilbert space $H$ in which $V$ is continuously imbedded and dense, we define $D(C) = \{v\in V:\C v\in H\}$ and $C=\C|_{D(C)}$. The corresponding restrictions of $\B$ and $\A$ to $H$ are denoted similarly. A (weak) solution $u$ of \eqn{631} for which each term belongs to $H$ at each $t\in \RR$ is called a {\it strong solution\/}, and it satisfies \begin{equation}\label{eq632} C u'' (t) + Bu' (t) +Au(t) = f(t)\ ,\qquad t\in \RR\ . \end{equation} \begin{theorem}\label{thm6-3B} Let the Hilbert space $V$ and operators $\A$, $\B$, $\C$ be given as in Theorem \ref{thm6-3A}. Let the Hilbert space $H$ and corresponding operators $A$, $B$, $C$ be defined as above, and assume $D(C)\subset D(A)\cap D(B)$. Then for every pair $u_0\in D(A)$, $u_1\in D(C)$, and $f\in C(\RR,H)$, there is a unique strong solution $u(\cdot)$ of \eqn{632} with $u(0)=u_0$, $u'(0)=u_1$. \end{theorem} \proof We define $M[x_1,x_2] = [Ax_1,Cx_2]$ on $D(A)\times D(C) = D(M)$ and $L[x_1,x_2] = [-Ax_2,Ax_1+Bx_2]$ on $D(A)\times D(A) \cap D(B)$ and apply Theorem V.\ref{thm5-3B}. \begin{corollary}\label{cor6-3C} In the situation of Theorem \ref{thm6-3B}, the weak solution $u(\cdot)$ is strong if and only if, for some $t_0\in \RR$, $u(t_0)\in D(A)$ and $u'(t_0)\in D(C)$. \end{corollary} We give an example which includes the classical Sobolev equation from fluid mechanics and an evolution equation of the type used to describe certain vibration problems. Let $G$ be open in $\RR^n$ and suppose that $\partial G$ is a $C^1$ manifold and that $\Gamma$ is a closed subset of $\partial G$. Let $V= \{v\in H^1(G): \gamma v(s)=0$, a.e.\ $s\in \Gamma\}$ and $$\C u(v) = (u,v)_{H^1(G)}\ ,\qquad u,v\in V\ .$$ Suppose $a_j(\cdot) \in L^\infty (G)$ for $j=1,2,\ldots,n$, and define $$\A u (v) = \sum_{j=1}^n \int_G a_j (x)\partial_j u(x)\partial_j \overline{v(x)}\,dx\ ,\qquad u,v\in V\ .$$ Let the functions $t\mapsto F(\cdot,t):\RR\to L^2(G)$ and $t\mapsto g(\cdot,t):\RR\to L^2(\partial G)$ be continuous and define $f\in C(\RR,V')$ by $$f(t) (v) =\int_G F(x,t)\overline{v(x)}\,dx + \int_{\partial G} g(s,t)\gamma\overline{v(s)}\,ds\ ,\qquad v\in V\ .$$ Then for each pair $U_0,U_1\in V$, we obtain from Theorems \ref{thm6-3A} and IV.\ref{thm4-7A} a unique generalized solution of the problem \begin{equation}\label{eq633} \left.\begin{array}{l} \ds \partial_t^2 U(x,t)-\Delta_n \partial_t^2U(x,t) -\sum_{j=1}^n \partial_j (a_j(x)\partial_j U(x,t)) = F(x,t)\ ,\\ \hfill x\in G\ ,\ t>0\ ,\\ \noalign{\vskip6pt} U(s,t)=0\ ,\qquad s\in\Gamma\ ,\\ \noalign{\vskip6pt} \ds\partial_\nu \partial_t^2 U(s,t) +\sum_{j=1}^n a_j(s)\partial_jU(s,t) = g(s,t)\ ,\qquad s\in \partial G\sim\Gamma\ ,\\ \noalign{\vskip6pt} U(x,0)=U_0(x)\ ,\quad \partial_t U(x,0)= U_1(x)\ .\end{array}\right\} \end{equation} In the special case of $a_j\equiv 0$, $1\le j\le n-1$, and $a_n(x)\equiv 1$, the partial differential equation in \eqn{633} is {\it Sobolev's equation\/} which describes inertial waves in rotating fluids. Terms due to temperature gradients will give \eqn{633} with $a_j(x)\equiv a>0$, $1\le j\le n-1$, and $a_n(x)\equiv 1$. Finally, if $a_j(x)\equiv a>0$, $1\le j\le n$, then the partial differential equation in \eqn{633} is {\it Love's equation\/} for longitudinal vibrations with lateral inertia. Suppose now that $g\equiv 0$ in the above, hence, $f\in C(\RR,H)$, where $H=L^2(G)$. If we assume $\Gamma=\partial G$, hence, $V=H_0^1(G)$, then $D(C) = H_0^1 (G)\cap H^2(G) \subset D(A)$, so Theorem \ref{thm6-3B} gives a smoother solution of \eqn{633} whenever $U_0,U_1\in D(C)$. If instead we assume $a_j(x)\equiv a$, $1\le j\le n$, then $D(C)=D(A)$, and Theorem \ref{thm6-3B} gives a smoother solution of \eqn{633} whenever $U_0,U_1\in D(C)$. Similar problems containing dissipation effects can easily be added, and we leave these to the exercises. In particular, there is motivation to consider problems like \eqn{633} with viscosity. \section{Degenerate Equations} % 4 \setcounter{equation}{0} We shall consider evolution equations of the form \eqn{621} wherein the leading operator $\C$ may not necessarily be the Riesz map of a Hilbert space. In particular, certain applications lead to \eqn{621} with $\C$ being symmetric and monotone. Our plan is to first solve a first order system like \eqn{622} by using one of Theorems V.\ref{thm5-4A} or V.\ref{thm5-4B}. Then the first and second components will be solutions (of appropriate modifications) of \eqn{621}. Also we shall obtain well-posed problems for a first order evolution equation in which the leading operator is not necessarily symmetric. (The results of Section~V.4 do not apply to such a situation.) \subsection{} % 4.1 Let $\A$ be the Riesz map of a Hilbert space $V$ to its dual $V'$. Let $\C\in\L (V,V')$ and suppose its sesquilinear form is symmetric and non-negative on $V$. Then it follows (cf., Section~V.4) that $x\mapsto \C x(x)^{1/2}$ is a seminorm on $V$; let $W$ denote the corresponding seminorm space. Finally, suppose $D(B)\le V$ and $B\in L(D(B),V')$ are given. Now we define $V_m$ to be the product $V\times W$ with the seminorm induced by the symmetric and non-negative sesquilinear form $$m(x,y) = \A x_1 (y_1) + \C x_2 (y_2)\ ,\qquad x,y\in V_m\equiv V\times W\ .$$ The identity $\M x(y) = m(x,y)$, $x,y\in V_m$, defines $\M\in\L(V_m,V'_m)$. Finally we define $D\equiv\{[x_1,x_2]\in V\times D(B): \A x_1 + Bx_2\in W'\}$ and the linear map $L:D\to V'_m$ by $$L[x_1,x_2] = [-\A x_2, \A x_1 + Bx_2]\ .$$ We shall apply Theorem V.\ref{thm5-4A} to obtain the following result. \begin{theorem}\label{thm6-4A} Let $\A$ be the Riesz map of the Hilbert space $V$ and let $W$ be the seminorm space obtained from a symmetric and monotone $\C \in \L(V,V')$. Let $D(B)\le V$ and $B\in L(D(B),V')$ be monotone. Assume $B+\C$ is strictly monotone and $\A+B+\C: D(B)\to V'$ is a surjection. Let $f\in C^1([0,\infty),W')$ and $g\in C^1([0,\infty),V')$. If $V_m$ and $D$ are the spaces denoted above, then for every pair $[u_0,u_1]\in D$ there exists a unique function $w(\cdot):[0,\infty)\to D$ such that $\M w(\cdot) \in C^1([0,\infty),V'_m)$, $\M w(0)=\M[u_0,u_1]$, and \begin{equation}\label{eq641} (\M w)'(t) + Lw (t) = [-g(t),f(t)]\ ,\qquad t\ge 0\ . \end{equation} \end{theorem} \proof We need to verify that the hypotheses of Theorem V.\ref{thm5-4A} are valid in this situation. First note that $K(\M)\cap D= \{[0,x_2]:x_2\in D(B)$, $Bx_2\in W'$, $\C x_2=0\}$. But if $y\in D(B)$ with $By\in W'$, then there is a $K\ge0$ such that $$|By (x)| \le K|\C x(x)|^{1/2}\ ,\qquad x\in V\ ,$$ hence, $|By(y)| \le K|\C y (y)|^{1/2} = 0$ if $\C y=0$. Thus, we have shown that $$\Re (B+\C) x_2 (x_2) = 0\ ,\qquad x= [0,x_2]\in K(\M) \cap D\ ,$$ so $B+\C$ being strictly monotone implies that $K(\M) \cap D= \{[0,0]\}$. Finally, just as in the proof of Theorem \ref{thm6-2A}, it follows from $\A+B+\C$ being surjective that $\M+L$ is surjective, so all the hypotheses of Theorem V.\ref{thm5-4A} are true. Let $w(\cdot)$ be the solution of \eqn{641} from Theorem \ref{thm6-4A} and set $w(t) = [u(t),v(t)]$ for each $t\ge0$. If we set $g\equiv0$ and eliminate $v(\cdot)$ from the system \eqn{641}, then we obtain an equivalent second order evolution equation which $u(\cdot)$ satisfies and, thereby, the following result. \begin{corollary}\label{cor6-4B} Let the spaces and operators be given as in Theorem \ref{thm6-4A}. For every $f\in C^1([0,\infty),W')$ and every pair $u_0\in V$, $u_1\in D(B)$ with $\A u_0+Bu_1 \in W'$ there exists a unique $u(\cdot)\in C^1([0,\infty), V)$ such that $\C u'(\cdot)\in C^1 ([0,\infty),W')$, $u(0)=u_0$, $\C u'(0)= \C u_1$, and for each $t\ge0$, $u'(t)\in D(B)$, $\A u(t)+Bu'(t)\in W'$, and \begin{equation}\label{eq642} (\C u'(t))' + Bu' (t) + \A u(t) = f(t)\ . \end{equation} \end{corollary} Similarly, the function $v(\cdot)$ obtained from a solution of \eqn{641} satisfies a second order equation. \begin{corollary}\label{cor6-4C} Let the spaces and operators be given as in Theorem \ref{thm6-4A}. If $F\in C([0,\infty),W')$, $g\in C^1([0,\infty),V')$, $u_1\in D(B)$ and $U_2\in W'$, then there exists a unique $v(\cdot):[0,\infty)\to D(B)$ such that $\C v (\cdot) \in C^1 ([0,\infty),W')$, $(\C v)'+Bv(\cdot) \in C^1([0,\infty),V')$, $\C v(0) = \C u_1$, $(\C v'+Bv)(0)=U_2+Bu_1$, and for each $t\ge0$, \begin{equation}\label{eq643} ((\C v)' (t) + Bv(t))' + \A v(t) = F(t) + g(t)\ . \end{equation} \end{corollary} \proof Given $F(\cdot)$ as above, define $f(\cdot) \in C^1([0,\infty),W')$ by $f(t) = \int_0^t F$. With $u_1$ and $U_2$ as above, there is a unique $u_0\in V$ for which $\A u_0 = -Bu_1-U_2$. Thus, $\A u_0 + Bu_1 \in W'$ so Theorem \ref{thm6-4A} gives a unique $w(\cdot)$ as indicated. Letting $w(t) \equiv [u(t),v(t)]$ for $t\ge0$, we have immediately $v(t)\in D(B)$ for $t\ge0$, $\C v\in C^1([0,\infty),W')$ and $\C v(0) = \C u_1$. The second line of \eqn{641} shows $$(\C v)' + Bv = f-\A u\in C^1([0,\infty),V')$$ and the choice of $u_0$ above gives $(\C v)' (0)+Bv(0)=U+Bu_1$. Eliminating $u(\cdot)$ from \eqn{641} gives \eqn{643}. This establishes the existence of $v(\cdot)$. The uniqueness result follows by defining $u(\cdot)$ by the second line of \eqn{641} and then noting that the function defined by $w(t)\equiv [u(t),v(t)]$ is a solution of \eqn{641}. Finally, we record the important special case of Corollary \ref{cor6-4C} that occurs when $\C=0$. This leads to a well-posed problem for a first order equation whose leading operator is not necessarily symmetric. \begin{corollary}\label{cor6-4D} Let the spaces $V$, $D(B)$ and operators $B$, $\A$ be given as in Theorem \ref{thm6-4A} but with $\C=0$, hence, $W'=\{0\}$. Then for every $g\in C^1([0,\infty),V')$ and $u_1\in D(B)$, there exists a unique $v:[0,\infty)\to D(B)$ such that $Bv(\cdot)\in C^1([0,\infty),V')$, $Bv(0)=Bu_1$, and for each $t\ge0$, \begin{equation}\label{eq644} (Bv)' (t) + \A v(t) = g(t)\ . \end{equation} \end{corollary} \subsection{} % 4.2 Each of the preceding results has a parabolic analogue. We begin with the following. \begin{theorem}\label{thm6-4E} Let $\A$ be the Riesz map of the Hilbert space $V$ and let $W$ be the seminorm space obtained from a symmetric and monotone $\C \in\L(V,V')$. Let $\B \in \L(V,V')$ be monotone and assume that $\B+\lambda\C$ is $V$-elliptic for some $\lambda>0$. Then for every pair of H\"older continuous functions $f:[0,\infty)\to W'$, $g:[0,\infty)\to V'$ and each pair $u_0\in V$, $U_1\in W'$, there exists a unique function $w:[0,\infty)\to V_m$ such that $\M w(\cdot)\in C([0,\infty),V'_m) \cap C^1 ((0,\infty),V'_m)$, $\M w(0) = [\A u_0,U_1]$, and for all $t>0$, $$(\M w)'(t) + \L w(t) = [-g(t),f(t)]\ ,$$ where $\L\in \L (V\times V,V'\times V')$ is defined by $\L[x_1,x_2] = [-\A x_2,\A x_1 + \B x_2]$, and $\M$ is given as in Theorem \ref{thm6-4A}. \end{theorem} \proof By introducing a change-of-variable, if necessary, we may replace $\L$ by $\lambda \M +\L$. Since for $x\equiv [x_1,x_2]\in V\times V$ we have $$\Re(\lambda\M + \L)x(x)=\lambda \A x_1(x_1) + (\B+\lambda \C) x_2(x_2)\ ,$$ we may assume $\L$ is $V\times V$-elliptic. The desired result follows from Theorem V.\ref{thm5-4B}. \begin{corollary}\label{cor6-4F} Let the spaces and operators be given as in Theorem \ref{thm6-4E}. For every H\"older continuous $f:[0,\infty)\to W'$, $u_0\in V$ and $U_1\in W'$, there exists a unique $u(\cdot) \in C([0,\infty),V)\cap C^1((0,\infty),V)$ such that $\C u'(\cdot) \in C((0,\infty),W')\cap C^1((0,\infty),W')$, $u(0)=u_0$, $\C u'(0)=U_1$, and \begin{equation}\label{eq645} (\C u'(t))' + \B u'(t) + \A u(t) = f(t)\ ,\qquad t>0\ . \end{equation} \end{corollary} \begin{corollary}\label{cor6-4G} Let the spaces and operators be given as in Theorem \ref{thm6-4E}. Suppose $F: (0,\infty)\to W'$ is continuous at all but a finite number of points and for some $p>1$ we have $\int_0^T \|F(t)\|_{W'}^p\,dt <\infty$ for all $T>0$. If $g:[0,\infty)\to V'$ is H\"older continuous, $u_1\in V$ and $U_2\in V'$, then there is a unique function $v(\cdot):[0,\infty)\to V$ such that $\C v\in C([0,\infty),W')\cap C^1((0,\infty),W')$, $(\C v)'+\B v\in C([0,\infty),V')$ and is continuously differentiable at all but a finite number of points, $\C v(0)=\C u_1$, $(\C v' +\B v)(0)=U_2+ \B u_1$, and \begin{equation}\label{eq646} ((\C v)'(t) +\B v(t))' + \A v(t) = F(t) +g(t) \end{equation} at those points at which the derivative exists. \end{corollary} \proof Almost everything follows as in Corollary \ref{cor6-4C}. The only difference is that we need to note that with $F(\cdot)$ as given above, the function $f(t) =\int_0^t F$ satisfies \begin{eqnarray*} \| f(t)-f(\tau)\|_{W'} &\le & \int_\tau^t \|F\|_{W'} \le |t-\tau|^{1/q} \biggl( \int_\tau^t \|F\|_{W'}^p\biggr)^{1/p}\\ \noalign{\vskip6pt} &\le &|t-\tau|^{1/q} \biggl( \int_0^T \|F\|_{W'}^p\biggr)^{1/p}\ , \qquad 0\le \tau\le t\le T\ ,\end{eqnarray*} where $1/q = 1-1/p \ge0$. Hence, $f$ is H\"older continuous. \section{Examples} % 5 \setcounter{equation}{0} We shall illustrate some applications of our preceding results by various examples of initial-boundary value problems. In each such example below, the operator $\A$ will correspond to one of the elliptic boundary value problems described in Section~III.4, and we refer to that section for the computations as well as occasional notations. Our emphasis here will be on the {\it types\/} of operators that can be chosen for the remaining coefficients in either of \eqn{642} or \eqn{643}. We begin by constructing the operator $\A$ from the abstract boundary value problem of Section~III.3. Let $V$, $H$ and $B$ be Hilbert spaces and $\gamma:V\to B$ a linear surjection with kernel $V_0$, and assume $\gamma$ factors into a norm-preserving isomorphism of $V/V_0$ onto $B$. Assume the injection $V\hookrightarrow H$ is continuous, $V_0$ is dense in $H$, and $H$ is identified with $H'$. Then we obtain the continuous injections $V_0\hookrightarrow H\hookrightarrow V'_0$ and $V\hookrightarrow H\hookrightarrow V'$ and $$(f,v)_H = f(v)\ ,\qquad f\in H\ ,\ v\in V\ .$$ Let $a_1 :V\times V\to \KK$ and $a_2 :B\times B\to \KK$ be continuous, sesquilinear and symmetric forms and define \begin{equation}\label{eq651} a(u,v)\equiv a_1 (u,v) + a_2 (\gamma u,\gamma v)\ ,\qquad u,v\in V\ . \end{equation} We shall assume $a(\cdot,\cdot)$ is $V$-elliptic; thus, $a(\cdot,\cdot)$ is a scalar-product on $V$ whose norm is equivalent to the original one on $V$. Hereafter, we shall take $a(\cdot,\cdot)$ as the scalar-product on $V$; the corresponding Riesz map $\A \in \L(V,V')$ is given by $$\A u(v) = a(u,v)\ ,\qquad u,v\in V\ .$$ Similarly, we define $A\in \L(V,V'_0)$ by \begin{equation}\label{eq652} Au(v) = a_1(u,v) \ ,\qquad u\in V\ ,\ v\in V_0\ , \end{equation} Let $D_0\equiv \{u\in V:Au\in H\}$, and denote by $\partial \in L(D_0,B')$ the abstract Green's operator constructed in Theorem III.\ref{thm3-2C} and characterized by the identity \begin{equation}\label{eq653} a_1(u,v) - (Au,v)_H = \partial u(\gamma(v))\ ,\qquad u\in D_0\ ,\ v\in V\ . \end{equation} Finally, we denote by $\A_2 \in \L (B,B')$ the operator given by $$\A_2 \varphi (\psi) = a_2(\varphi,\psi)\ ,\qquad \varphi,\psi\in B\ .$$ It follows from \eqn{651}, \eqn{652} and \eqn{653} that \begin{equation}\label{eq654} \A u(v) - (Au,v)_H = (\partial u+\A_2 (\gamma u))(\gamma v)\ ,\qquad u\in D_0\ ,\ v\in V\ , \end{equation} and this identity will be used to characterize the weak or variational boundary conditions below. Let $c:H\times H\to \KK$ be continuous, non-negative, sesquilinear and symmetric; define the monotone $\C \in \L(H)$ by $$\C u(v) = c(u,v)\ ,\qquad u,v\in H\ ,$$ where $\C u\in H$ follows from $H'=H$. Note that the inclusion $W'\subset H$ follows from the continuity of the injection $H\hookrightarrow W$, where $W$ is the space $H$ with seminorm induced by $c(\cdot,\cdot)$. Finally let $\B\in\L(V,H)$ be a given monotone operator $$\Re \B u(v) \ge 0\ ,\qquad u\in V\ ,\ v\in H\ ,$$ and assume $\C +\B$ is strictly monotone: $$(\C +\B)u(u) = 0\ \hbox{ only if }\ u=0\ .$$ \begin{theorem}\label{thm6-5A} Let the Hilbert spaces and operators be given as above. For every $f\in C^1([0,\infty),W')$ and every pair $u_0,u_1\in V$ with $\A u_0 +\B u_1 \in W'$, there exists a unique $u\in C^1([0,\infty),V)$ such that $\C u'\in C^1([0,\infty),W')$, $u(0)=u_0$, $\C u'(0)=\C u_1$, and for each $t\ge0$, \begin{eqnarray} &&(\C u'(t))' + \B u'(t) +Au (t) = f(t)\ ,\label{eq655}\\ \noalign{\vskip6pt} &&u(t)\in D_0\subset V\ ,\quad \partial u(t) +\A_2\gamma(u(t))=0\ . \label{eq656}\end{eqnarray} \end{theorem} \proof The existence and uniqueness of $u(\cdot)$ follows from Corollary \ref{cor6-4B}. With $\C$ and $\B$ as above \eqn{642} shows that $\A u(t)\in H$ for all $t\ge0$, so \eqn{656} follows from Corollary III.\ref{cor3-3B}. (Cf. \eqn{654}.) To be sure, the pair of equations \eqn{655}, \eqn{656}, is equivalent to \eqn{642}. \qed We illustrate Theorem \ref{thm6-5A} in the examples following in Sections~5.1 and 5.2. \subsection{} %5.1 Let $G$ be open in $\RR^n$, $H=L^2(G)$, $\Gamma\subset \partial G$ and $V=\{v\in H^1(G): \gamma_0 (v)(s)=0$, a.e.\ $s\in \Gamma\}$. Let $p\in L^\infty (G)$ with $p(x) \ge0$, $x\in G$, and define \begin{equation}\label{eq657} c(u,v) = \int_G p(x)u(x)\overline{v(x)}\, dx\ ,\qquad u,v\in H\ . \end{equation} Then $\C$ is multiplication by $p$ and $W'=\{p^{1/2} v:v\in L^2(G)\}$. Let $R\in L^\infty (G)$ and the real vector field $\mu(x)= (\mu_1(x),\ldots, \mu_n(x))$ be given with each $\mu_j \in C^1(\bar G)$; assume $$\begin{array}{rcll} \ds -\left(\tfrac12\right) \sum_{j=1}^n \partial_j \mu_j (x) + \Re \{R(x)\} &\ge&0\ ,&\qquad x\in G\ ,\\ \noalign{\vskip6pt} \ds \left(\tfrac12\right) \mu(s)\cdot\nu (s)&\ge&0\ ,&\qquad s\in\partial G \sim \Gamma\ .\end{array}$$ Then $\B \in \L(V,H)$ given by \eqn{627} is monotone. Furthermore, we shall assume $$p(x) -\left(\tfrac12\right) \sum_{j=1}^n \partial_j \mu_j (x) + \Re \{R(x)\} >0\ ,\qquad x\in G\ ,$$ and this implies $\C+\B$ is strictly-monotone. Let $a_0,a_{ij}\in L^\infty (G)$, $1\le i,j\le n$, and assume $a_0(x) \ge0$, $a_{ij}(x)= \overline{a_{ji}(x)}$, $x\in G$, and that \begin{equation}\label{eq658} a(u,v) \equiv \int_G \biggl\{ \sum_{i,j=1}^n a_{ij} (x)\partial_i u(x) \partial_j \overline{v(x)} + a_0 (x) u(x) \overline{v(x)}\biggr\}\,dx \end{equation} is $V$-coercive (cf.\ Section III.5). Then \eqn{658} is a scalar product on $V$ whose norm is equivalent to that of $H^1(G)$ on $V$. Let $F(\cdot,t)\in L^2(G)$ be given for each $t\ge0$ such that $t\mapsto F(\cdot,t)$ belongs to $C^1([0,\infty),L^2(G))$ (cf.\ Section 2.1). Then $f(t)\equiv p^{1/2} F(\cdot,t)$ defines $f\in C^1([0,\infty),W')$. Finally, let $U_0,U_1\in V$ satisfy $\A U_0 +\B U_1\in W'$. (This can be translated into an elliptic boundary value problem.) Using Theorem IV.\ref{thm4-7A}, we can obtain a (measurable) function $U(\cdot,\cdot)$ on $G\times [0,\infty)$ which is a solution of the initial-boundary value problem \begin{eqnarray} &&\partial_t (p(x)\partial_t U(x,t)) + R(x)\partial_t U(x,t) + {\partial\over\partial\mu} (\partial_t U(x,t)) \label{eq659}\\ \noalign{\vskip6pt} &&\qquad\qquad -\sum_{j=1}^n \partial_j a_{ij}(x) \partial_i U(x,t) + a_0(x) U(x,t) \nonumber\\ \noalign{\vskip6pt} &&\qquad\qquad = p^{1/2} (x)F(x,t)\ , \qquad x\in G\ ,\ t\ge0\ ;\nonumber\\ \noalign{\vskip6pt} &&\left.\begin{array}{ll} U(s,t) =0\ ,&\qquad s\in\Gamma\ ,\\ \noalign{\vskip6pt} \ds {\partial U(s,t)\over\partial \nu_A}=0\ , &\qquad s\in\partial G\sim\Gamma\ ;\end{array} \right\} \label{eq6510}\\ \noalign{\vskip6pt} &&\left.\begin{array}{ll} U(x,0)= U_0(x)\ ,&\\ \noalign{\vskip6pt} p(x)\partial_tU(x,0) = p(x)U_1(x)\ ,&\qquad x\in G\ .\end{array}\right\} \label{eq6511} \end{eqnarray} We refer to Section III.4.1 for notation and computations involving the operators associated with the form \eqn{658}. The partial differential equation \eqn{659} is of mixed hyperbolic-parabolic type. Note that the initial conditions \eqn{6511} imposed on the solution at $x\in G$ depend on whether $p(x)>0$ or $p(x)=0$. Also, the equation \eqn{659} is satisfied at $t=0$, thereby imposing a compatibility condition on the initial data $U_0,U_1$. Finally, we observe that \eqn{6510} contains the boundary condition of {\it first type\/} along $\Gamma$ and the boundary condition of {\it second type\/} on $\partial G\sim\Gamma$. \subsection{} % 5.2 Let $H$ and $\C$ be given as in Section 5.1; let $V=H^1(G)$ and define $\B$ by \eqn{627} with $\mu\equiv 0$ and assume $$\begin{array}{rcll} \Re\{ R(x)\}&\ge&0\ ,&\qquad p(x)\ge0\ ,\\ \noalign{\vskip6pt} p(x)+\Re \{R(x)\}&>&0\ ,&\qquad x\in G\ ,\end{array}$$ as before. Define $$\begin{array}{rcll} a_1(u,v)&=&\int_G \nabla u\cdot\nabla \bar v&\qquad u,v\in V\ ,\\ \noalign{\vskip6pt} a_2(\varphi,\psi)&=&\int_{\partial G} \alpha (s) \varphi (s) \overline{\psi (x)}\, ds\ ,&\qquad \varphi,\psi\in L^2(\partial G)\end{array}$$ where $\alpha\in L^\infty (\partial G)$, $\alpha (s) \ge0$, a.e.\ $s\in\partial G$. Then $\A_2$ is multiplication by $\alpha$. We assume that $a(\cdot,\cdot)$ given by \eqn{651} is $V$-coercive (cf. Corollary III.\ref{cor3-5E}). With $F(\cdot,\cdot)$, $U_0$, and $U_1$ as above, we obtain a unique generalized solution of the problem \begin{eqnarray} &&\partial_t (p(x)\partial_t U(x,t)) + R(x)\partial_t U(x,t) - \Delta_n U(x,t) \label{eq6512}\\ &&\hskip1in = p^{1/2} (x)F(x,t)\ ,\qquad x\in G\ ,\ t\ge0\ ,\nonumber\\ \noalign{\vskip6pt} &&{\partial U(s,t)\over\partial \nu} + \alpha (s)U(s,t)=0\ , \qquad\qquad s\in\partial G\ ,\ t\ge0\ ,\label{eq6513}\end{eqnarray} and \eqn{6511}. We note that at those $x\in G$ where $p(x)>0$, \eqn{6512} is a (hyperbolic) wave equation and \eqn{6511} specifies initially $U$ and $\partial_t U$, whereas at those $x\in G$ where $p(x)=0$, \eqn{6512} is a homogeneous (parabolic) diffusion equation and only $U$ is specified initially. The condition \eqn{6513} is the boundary condition of {\it third type\/}. If we choose $V= \{v\in H^1(G):\gamma_0(v)$ is constant$\}$ as in Section III.4.2 and prescribe everything else as above, then we obtain a solution of \eqn{6512}, \eqn{6511} and the boundary condition of {\it fourth type\/} \begin{equation}\label{eq6514} \left.\begin{array}{l} U(s,t)=h(t)\ ,\qquad s\in \partial G\ ,\\ \noalign{\vskip6pt} \ds \int_{\partial G} {\partial U(s,t)\over\partial\nu}\, ds + \int_{\partial G} \alpha (s)\,ds \cdot h(t)=0\ .\end{array}\right\} \end{equation} Note that $h(\cdot)$ is an unknown in the problem. Boundary value problems with {\it periodic boundary conditions\/} can be put in the form of \eqn{6514}. \subsection{} % 5.3 Let $H= L^2(G)$, $V= H^1(G)$, and define $\A$ as in Section 5.2. Set $\B \equiv 0$ and define $$c(u,v) = \int_G p(x)u(x)\overline{v(x)}\,dx + \int_{\partial G} \sigma (s) u(s)\overline{v(s)}\,ds\ ,\qquad u,v\in V$$ when $p\in L^\infty (G)$ and $\sigma \in L^\infty (\partial G)$ satisfy $p(x)>0$, $x\in G$, and $\sigma (s)\ge0$, $s\in\partial G$. Let $t\mapsto F(\cdot,t)$ be given in $C^1 ([0,\infty),L^2(G))$ and $t\mapsto g(\cdot,t)$ be given in $C^1([0,\infty),L^2(\partial G))$; then define $f\in C^1([0,\infty),W')$ by $$f(t)(v) = \int_G p^{1/2} (x)F(x,t)\overline{v(x)}\,dx + \int_{\partial G}\!\!\! \sigma^{1/2} (s)g(s,t)\overline{v(s)}\,ds\ , \quad v\in V\ ,\ t\ge0\ .$$ Let $U_0,U_1\in V$ with $\A U_0\in W'$. (This last inclusion is equivalent to requiring $\Delta_nU_0=p^{1/2}H$ for some $H\in L^2(G)$ and $\partial_\nu U_0 +\alpha U_0 = \sigma^{1/2}h$ for some $h\in L^2(\partial G)$.) Then Corollary \ref{cor6-4B} applies to give a unique solution $u$ of \eqn{642} with initial conditions. From this we obtain a solution $U$ of the problem $$\cases{\partial_t (p(x)\partial_t U(x,t)) - \Delta_n U(x,t) = p^{1/2}(x)F(x,t)\ ,&\cr \noalign{\vskip6pt} \hfill x\in G\ ,\ t\ge 0\ ,&\cr \noalign{\vskip6pt} \partial_t (\sigma (s)\partial_t U(s,t)) + \partial_\nu U(s,t) +\alpha (s)U(s,t)= \sigma^{1/2} (s) g(s,t)\ ,&\cr \noalign{\vskip6pt} \hfill s\in\partial G\ ,\ t\ge0\ ,&\cr \noalign{\vskip6pt} U(x,0)= U_0(x)\ ,\quad \partial_tU(x,0)= U_1(x)\ .&\cr}$$ The boundary condition is obtained {\it formally\/} since we do not know $\Delta_nU(\cdot,t)\in L^2(G)$ for all $t>0$; hence, \eqn{653} is not directly applicable. Such boundary conditions arise in models of vibrating membranes (or strings) with boundaries (or ends) loaded with a mass distribution, thereby introducing an inertia term. Such problems could also contain mass distributions (or point loads) on internal regions. Similarly, internal or boundary damping can be included by appropriate choices of $\B$, and we illustrate this in the following example. \subsection{} % 5.4 Let $H$, $V$, $\A$ and $\C$ be given as in Section 5.3. Assume $R\in L^\infty(G)$, $r\in L^\infty (\partial G)$ and that $\Re \{R(x)\} \ge0$, $x\in G$, $\Re \{r(s)\}\ge0$, $s\in\partial G$. We define $\B\in \L(V,V')$ by $$\B u(v)= \int_G R(x)u(x)\overline{v(x)}\,dx + \int_{\partial G} r(s)u(s)\overline{v(s)}\,ds\ ,\qquad u,v\in V\ .$$ We need only to assume $p(x) + \Re \{R(x)\}>0$ for $x\in G$; then Corollary \ref{cor6-4C} is applicable. Let $t\mapsto F_1(\cdot,t)$ in $C([0,\infty),L^2(G))$, $t\mapsto G_1(\cdot,t)$ in $C^1([0,\infty),L^2(\partial G))$, and $t\mapsto G_2(t)$ in $C^1([0,\infty),L^2(G))$ be given. We then define $F\in C([0,\infty),W')$ and $g\in C^1([0,\infty),V')$ by \begin{eqnarray*} F(t) &=&p^{1/2} F_1(\cdot,t)\ ,\cr \noalign{\vskip6pt} g(t)(v)&=& \int_{\partial G} \sigma^{1/2} (s) G_1(s,t) \overline{v(s)}\,ds + \int_G G_2 (x,t)\overline{v(x)}\,dx \ ,\qquad v\in V\ .\end{eqnarray*} If $U_1\in V$ and $V_1\in L^2(G)$, and $V_2\in L^2(\partial G)$, then $U_2\in W'$ is defined by $$U_2 (v) = \int_G p^{1/2} (x) V_1 (x)\overline{v(x)}\,dx + \int_{\partial G}\sigma^{1/2} (s)V_2 (s)\overline{v(s)}\,ds\ ,\qquad v\in V\ ,$$ and Corollary \ref{cor6-4C} gives a generalized solution of the following problem: $$\cases{\partial_t^2 (p(x)U(x,t)) + \partial_t (R(x)U(x,t))-\Delta_n U(x,t)&\cr \noalign{\vskip6pt} \hfill = p^{1/2} (x)F_1(x,t) +G_2(x,t)\ ,\qquad x\in G\ ,&\cr \noalign{\vskip6pt} \partial_t^2 (\sigma (s)U(s,t)) +\partial_t (r(s)U(s,t)) + \partial_\nu U(s,t) + (s)U(s,t) &\cr \noalign{\vskip6pt} \hfill = \sigma^{1/2}(s)G_1(s,t)\ ,\ \ \qquad s\in \partial G\ ,\ t>0\ ,&\cr \noalign{\vskip6pt} p(x)U(x,0) = p(x)U_1(x)\ ,\cr \noalign{\vskip6pt} \sigma(s) U(s,0) = \sigma (s)U_1(s)\ ,\qquad s\in \partial G\cr \noalign{\vskip6pt} \partial_t(p(x)U(x,0)) + R(x) U(x,0) = p^{1/2} (x)V_1 (x)\ ,&\cr \noalign{\vskip6pt} \partial_t (\sigma (s)U(s,0)) + r(s) U(s,0) = \sigma^{1/2}(s)V_2(s)\ .&\cr}$$ The right side of the partial differential equation could contain singularities in $x$ as well. When $\Re \{R(x)\}>0$ in $G$, the preceding problem with $p\equiv 0$ and $\sigma\equiv 0$ is solved by Corollary \ref{cor6-4D}. Similarly one can obtain generalized solutions to boundary value problems containing partial differential equations of the type \eqn{633}; that is, equations of the form \eqn{659} plus the fourth-order term $-\partial_t (\Delta_n \partial_t U(x,t))$. Finally, we record an abstract parabolic boundary value problem which is solved by using Corollary \ref{cor6-4F}. Such problems arise in classical models of linear viscoelasticity (cf. \eqn{629}). \begin{theorem}\label{thm6-5B} Let the Hilbert spaces and operators be given as in Theorem \ref{thm6-5A}, except we do not assume $\B +C$ is strictly monotone. If $\varep>0$, $f:[0,\infty)\to W'$ is H\"older continuous, $u_0\in V$ and $U_1\in W'$, there exists a unique $u\in C([0,\infty),V)\cap C^1([0,\infty), V)$ such that $\C u'\in C([0,\infty),W')\cap C^1((0,\infty),W')$, $u(0)=u_0$, $Cu'(0)=U_1$, and \eqn{655}, \eqn{656} hold for each $t>0$. \end{theorem} \exercises \begin{description} \item[1.1.] Use the separation-of-variables technique to obtain a series representation for the solution $u$ of \eqn{611} with $u(0,t)=u(\pi,t)=0$, $u(x,0)=u_0(x)$ and $\partial_t u(x,0)=u_1(x)$. \item[1.2.] Repeat the above for the viscoelasticity equation $$\partial_t^2 u-\varep \partial_t \Delta_n u-\Delta_n u= F(x,t)\ , \qquad \varep >0\ .$$ \item[1.3.] Compare the convergence rates of the two series solutions obtained above. \medskip \item[2.1.] Explain the identification $V'_m=V'\times W'$ in Section 2.1. \item[2.2.] Use Theorem \ref{thm6-2A} to prove Theorem \ref{thm6-2C}. \item[2.3.] Use the techniques of Section 2 to deduce Theorem \ref{thm6-2C} from IV.5. \item[2.4.] Verify that the function $f$ in Section 2.1 belongs to $C^1([0,T],L^2(G))$. \item[2.5.] Use Theorem \ref{thm6-2A} to construct a solution of \eqn{625} satisfying the fourth boundary condition. Repeat for each of the examples in Section III.4. \item[2.6.] Add the term $\int_{\partial G} r(s) u(s)\overline{v(s)}\,ds$ to \eqn{627} and find the initial-boundary value problem that results. \item[2.7.] Show that Theorem \ref{thm6-2A} applies to appropriate problems for the equation $$\partial_t^2 u(x,t) + \partial_x^3 \partial_t u(x,t)- \partial_x^2 u(x,t) = F(x,t)\ .$$ \item[2.8.] Find some well-posed problems for the equation $$\partial_t^2 u(x,t) + \partial_x^4 u(x,t) = F(x,t)\ .$$ \item[3.1.] Complete the proofs of Theorem \ref{thm6-3B} and Corollary \ref{cor6-3C}. \item[3.2.] Verify that \eqn{633} is the characterization of \eqn{631} with the given data. \medskip \item[4.1.] Use Corollary \ref{cor6-4D} to solve the problem $$\begin{array}{l} \partial_t\partial_x u(x,t)-\partial_x^2 u(x,t) = F(x,t)\\ \noalign{\vskip6pt} u(0,t) = cu(1,t)\\ \noalign{\vskip6pt} u(x,0)= u_0(x)\end{array}$$ for $|c|\le 1$, $c\ne1$. \item[4.2.] For each of the Corollaries of Section 4, give an example which illustrates a problem solved by that Corollary only. \medskip \item[5.1.] In the proof of Theorem \ref{thm6-5A}, verify that \eqn{642} is equivalent to the pair \eqn{655}, \eqn{656}. \item[5.2.] In Section 5.1, show $\C+\B$ is strictly monotone, give sufficient conditions for \eqn{658} to be $V$-elliptic, and characterize the condition $\A U_0+ \B U_1 \in W'$ as requiring that $U_0$ satisfy an elliptic boundary value problem (cf. Section 5.3). \item[5.3.] In Section 5.2, give sufficient conditions for $a(\cdot,\cdot)$ to be $V$-elliptic. \item[5.4.] Show the following problem with periodic boundary conditions is well-posed: $\partial_t^2 u-\partial_x^2 u = F(x,t)$, $u(x,0) = u_0(x)$, $\partial_t u(x,0) = u_1(x,0)$, $u(0,t) = u(1,t)$, $\partial_x u(0,t) = \partial_x u(1,t)$. Generalize this to higher dimensions. \item[5.5.] A vibrating string loaded with a point mass $m$ at $x=\tfrac12$ leads to the following problem: $\partial_t^2 u= \partial_x^2 u$, $u(0,t)= u(1,t)=0$, $u(x,0)= u_0(x)$, $\partial_t u(x,0) = u_1(x)$, $u((\tfrac12)^-,t) = u((\tfrac12)^+,t)$, $m\partial_t^2 u(\tfrac12,t) = \partial_xu((\tfrac12)^+,t) - \partial_x u((\tfrac12)^-,t)$. Use the methods of Section 5.3 to show this problem is well-posed. \end{description}