Abstract
In this work we present a unified error analysis for abstract space discretizations of wave-type equations with nonlinear quasi-monotone operators. This yields an error bound in terms of discretization and interpolation errors that can be applied to various equations and space discretizations fitting in the abstract setting. We use the unified error analysis to prove novel convergence rates for a non-conforming finite element space discretization of wave equations with nonlinear acoustic boundary conditions and illustrate the error bound by a numerical experiment.
Similar content being viewed by others
1 Introduction
In this paper we present a unified error analysis for abstract non-conforming space discretizations of nonlinear wave-type equations with quasi-monotone operators. The unified error analysis was introduced in [15, 16] for linear wave equations and extended in [18] to semilinear problems. It is an abstract framework in which wave equations as well as a variety of spatial discretizations are considered as evolution equations in Hilbert spaces. Using such an abstract framework allows to derive an abstract error bound in terms of approximation properties of the space discretization method. This error bound can then be used to prove convergence rates for all space discretizations of wave-type equations which fit into the abstract setting. This was demonstrated in [16] for isoparametric finite element discretizations of wave equations with linear acoustic boundary conditions and a dG discretization of Maxwell equations.
The aim of this paper is to extend the unified error analysis to nonlinear evolution equations with quasi-monotone operators. As a specific application, we use this theory to prove error bounds for a non-conforming finite element discretization of the wave equation with nonlinear acoustic boundary conditions. This is a generalization of the results in the thesis [20].
Acoustic boundary conditions were first mentioned in [5]. Since then, many papers studied their properties, wellposedness, and stability, and they are still in the focus of current research, cf.[4, 10, 11, 21, 23, 26] and references therein.
However, there are only very few numerical papers considering these boundary conditions. We are aware of [16] and [17]. In these papers space discretizations for wave equations with linear acoustic boundary conditions were derived and analyzed in the energy and the \(L^2\)-norm, respectively. In the present paper, we now consider the space discretization of nonlinear acoustic boundary conditions as proposed in [13, 14, 28], and extend the results from [16] to this case.
Since acoustic boundary conditions include derivatives on the boundary, they are usually posed on domains with smooth boundaries. A common choice to discretize such problems are isoparametric finite elements. Since this involves to approximate the boundary of the domain, the discretization becomes non-conforming. Unfortunately, this makes the error analysis much more involved since the exact and the numerical solution are not defined on the same domain which causes errors in the bilinear form ot the weak form.
We derive an order p isoparametric finite element discretization of the wave equation with nonlinear acoustic boundary conditions and show that it fits into the setting of the unified error analysis. Using the abstract error bound we prove order p convergence of the method in the energy norm, where we tackle the appearing approximation errors stemming from the domain approximation or interpolation by known error bounds from [7, 8]. A major difficulty lies in the discretization of the nonlinearities, since this must be done in such a way that it preserves the quasi-monotonicity to ensure the stability of the numerical scheme. Furthermore, the discretization errors of the nonlinearities have to be bounded. While both is straightforward for conforming discretizations, it turns out to be much more involved in the non-conforming case.
We are not aware of any other results in this direction, neither of such a general error analysis for non-conforming space discretizations of nonlinear wave-type equations, nor of results concerning the discretization of wave equations with nonlinear acoustic boundary conditions. Nevertheless, we mention the following works going in the same direction. In [9], a full discretization in an abstract framework similar to the one used in this paper was considered. But only a conforming space discretization was analyzed and no error bounds but only weak convergence of the discretization was shown. For quasilinear equations, a related framework was introduced in [19, 22], covering quasilinear wave and Maxwell equations. However, the error analysis in this work relies on properties of quasilinear operators that cannot be used for nonlinear acoustic boundary conditions and in general for equations with maximal quasi-monotone operators.
This paper is structured as follows. In Sect. 2 we introduce the wave equation with nonlinear acoustic boundary conditions with a corresponding finite element space discretization and state an error bound of the spatial discretization. We then present in Sect. 3 the unified error analysis for nonlinear first-order evolution equations and use the results in Sect. 4 to analyze nonlinear second-order wave-type equations. As main results we derive abstract error bounds for the space discretizations. Finally, in Sect. 5 we use these abstract bounds of the unified error analysis to prove the space discretization error bound for the wave equations with nonlinear acoustic boundary conditions and illustrate it with some numerical experiments.
2 The wave equation with nonlinear acoustic boundary conditions
In this section we present the analytical framework for the wave equation with acoustic boundary conditions and a suitable finite element space discretization. Additionally, we present a space discretization error bound which we will prove by application of the unified error analysis in Sect. 5.
2.1 Problem statement and analytical framework
Let \(\varOmega \subset \mathbb {R}^n, n=2,3\), be a bounded domain with \(C^2\)-boundary \(\varGamma \) and outer normal vector \(\textbf{n}\). We consider the acousitc wave equation with non-local reacting acoustic boundary conditions in the following form: seek \(u:[0,T]\times \varOmega \rightarrow \mathbb {R}, \delta :[0,T]\times \varGamma \rightarrow \mathbb {R}\) satisfying
Here \(\varDelta _\varGamma \) denotes the Laplace–Beltrami operator an \(\varGamma \).
Remark 2.1
It is possible to include nonlinear forcing terms \(F_{\varOmega }(\textbf{x},u)\) and \(F_{\varGamma }(\textbf{x},\delta )\) at the right-hand side of (1a) and (1b), respectively. This was considered in [20] for the wave equation with kinetic boundary conditions and such terms can be treated similarly for the acoustic boundary conditions. We omit this here for the sake of a clearer presentation.
We make the following assumptions on the coefficients and nonlinearities in (1).
Assumption 2.2
-
a)
The constants satisfy \(c_\varOmega ,c_\varGamma , \mu > 0, \quad k_\varOmega ,k_\varGamma \ge 0, \quad d, \rho \in \mathbb {R}\).
-
b)
The function \(\theta \in C(\mathbb {R};\mathbb {R})\) satisfies \(\theta (0)=0\) and is strictly monotonically increasing with
$$\begin{aligned} \left( \theta (\xi _1) - \theta (\xi _2)\right) (\xi _1 - \xi _2)\ge \theta _0 \left| \xi _1 - \xi _2\right| ^2, \qquad \xi _1,\xi _2 \in \mathbb {R}, \end{aligned}$$for some \(\theta _0>0\). Further, there exist
$$\begin{aligned} 1\le \zeta {\left\{ \begin{array}{ll} < \infty , &{} n=2, \\ \le 3, &{} n= 3, \end{array}\right. } \end{aligned}$$(2)and a constant \(C>0\) such that for all \(\xi \in \mathbb {R}\)
$$\begin{aligned} \left| \theta (\xi )\right| \le C (1 + \left| \xi \right| ^{\zeta }). \end{aligned}$$(3) -
c)
The function \(\eta :\mathbb {R}\rightarrow \mathbb {R}\) is globally Lipschitz continuous and satisfies \(\eta (0) = 0\). We then have that \(\tilde{\eta }\) defined via \(\tilde{\eta }(\xi ) = \eta (\xi ) - \frac{\rho }{c_\varOmega } \xi \) is also Lipschitz continuous and denote the Lipschitz constant of \(\tilde{\eta }\) by \(L_{\eta }\).
-
d)
The inhomogeneities satisfy \(f_{\varOmega }\in W_{loc}^{1,1}([0,\infty );C({\overline{\varOmega }}))\) and \(f_{\varGamma }\in W_{loc}^{1,1}([0,\infty );C(\varGamma ))\).
Weak formulation To prove wellposedness and derive a finite element discretization, we now present a weak formulation of the wave equation with acoustic boundary conditions (1). We make use of the densely embedded Hilbert spaces
where
Note that for \(k\ge 3\), the spaces \(H^k(\varGamma )\) require more boundary regularity to be well-defined, e.g., \(\varGamma \in C^k\), which denotes that \(\varGamma \) is a \(C^k\) boundary.
By multiplying (1a) and (1b) with test functions defined on \(\varOmega \) and \(\varGamma \), respectively, applying integration by parts and inserting the nonlinear coupling (1c), we obtain the the weak formulation of (1): seek \(\textbf{u}= \left[ u,\delta \right] ^\intercal \in C^2([0,T];H)\cap C^1([0,T];V)\) satisfying
where for \(\textbf{v}= \left[ v,z \right] ^\intercal , {\varvec{\varphi }}= \left[ \varphi , \psi \right] ^\intercal \in V\) we have
Note that \(m\) is an inner product on H and \({{\tilde{a}}}{:}{=}a+ m\) is an inner product on V.
Remark 2.3
Assumption 2.2 ensures that (4) is globally wellposed, we comment on this in Sect. 4.1, cf.Corollary 4.4, and Sect. 5.
2.2 Finite element space discretization
For the space discretization of (1) we consider the bulk-surface finite element method from [7] which was also used in [16] to discretize the wave equation with linear acoustic boundary conditions. We give a brief introduction of the finite element spaces and refer to [7] for further details on the bulk-surface finite element method.
The bulk-surface finite element method Let \(\varGamma \in C^{p+1}\) for some \(p\ge 1\) and let \({\mathcal {T}}_h^\varOmega \) be a consistent and quasi-uniform mesh consisting of isoparametric elements \(K\) of degree p which discretizes \(\varOmega \). By h we denote the maximal mesh width of \({\mathcal {T}}_h^\varOmega \). The discretized domain is then given by
and its boundary by \(\varGamma _h= \partial \varOmega _h\). The bulk and the surface finite element space of order p are then defined by
respectively. Here, \(\widehat{K}\) denotes the reference triangle with corresponding polynomial space \(\mathbb {P}_p(\widehat{K})\) of order p, and \(F_{K}\) is the transformation from \(\widehat{K}\) to \(K\). Note that by construction we have \(v_h\big |_{{\varGamma _h}}\in V_{h,p}^\varGamma \) for all \(v_h\in V_{h,p}^\varOmega \).
As approximation space for V we set \(V_h= V_{h,p}^\varOmega \times V_{h,p}^\varGamma \). Note that, since \(\varOmega _h\) is only an approximation of \(\varOmega \), we have \(V_h \nsubseteq V\), i.e., the discretization is non-conforming. Hence, to relate functions in \(V_h\) with functions in V, in [7], for \(\textbf{v}_h= \left[ v_h,\vartheta _h \right] ^\intercal \in V_h\) and \(h<h_0\) sufficiently small, a lifted version
was constructed. By \(I_{h,\varOmega }:C({\overline{\varOmega }})\rightarrow V_{h,p}^\varOmega \) and \(I_{h,\varGamma }:C(\varGamma ) \rightarrow V_{h,p}^\varGamma \) we denote the order p nodal interpolation operators in \(\varOmega \) and on \(\varGamma \), respectively, and set for \(\textbf{v}= \left[ v,\vartheta \right] ^\intercal \in V\)
The spatially discretized equation We now state the finite element discretization of (1). For this, let
be an elementwise defined quadrature formula that approximates the integral \(\int _{\varGamma _h}\cdot \, \textrm{d}\textbf{s}\). We require that the quadrature formula has positive weights and is of order greater than 2p, s.t.polynomials up to degree 2p are integrated exactly and we have for all \(z_h,\psi _h \in V_{h,p}^\varGamma \)
For \(\textbf{v}_h= \left[ v_h,z_h \right] ^\intercal , \varvec{\mathbf {\varphi }_h}= \left[ \varphi _h, \psi _h \right] ^\intercal \in V_h\) we define
Then, the spatial discretization of (1) is given by: seek \(\textbf{u}_h:[0,T]\rightarrow V_h\) s.t.
Remark 2.4
The use of the quadrature formulas instead of the interpolation in the definition of the discretized nonlinearity \({\mathcal {D}}_h\) is required to prove that \({\mathcal {D}}_h\) is quasi-monotone, cf.Lemma 5.3.
To prove an error bound of the discretization we pose the following assumptions on the exact solution and the data:
Assumption 2.5
-
a)
Let \(T>0\). For the inhomogeneities and the nonlinearities in (1) we assume the additional regularity
$$\begin{aligned} f_{\varOmega }\in L^\infty \big ([0,T];H^{\max \lbrace 2 , p \rbrace }(\varOmega )\big ),{} & {} \quad f_{\varGamma }\in L^\infty \big ([0,T];H^{\max \lbrace 2 , p \rbrace }(\varGamma )\big ), \end{aligned}$$(9a)$$\begin{aligned} \theta ,\eta \in C^{\max \lbrace 2 , p \rbrace }(\mathbb {R};\mathbb {R}).{} & {} \end{aligned}$$(9b)Furthermore, we assume that the strong solution \(u,\delta \) of (1) satisfies on [0, T]
$$\begin{aligned} u&\in L^\infty \big ([0,T];H^{p+1}(\varOmega )\big ),&u'&\in L^\infty \big ([0,T]; H^{p+2}(\varOmega ) \cap W^{p+1,\infty }(\varOmega )\big ),\\ u''&\in L^\infty \big ([0,T];H^{\max \lbrace 2 , p \rbrace }(\varOmega )\big ),\\ \delta&\in L^\infty \big ([0,T];H^{p+1}(\varGamma )\big ),&\delta '&\in L^\infty \big ([0,T]; H^{p+1}(\varGamma )\cap W^{p,\infty }(\varGamma )\big ),\\ \delta ''&\in L^\infty \big ([0,T];H^{\max \lbrace 2 , p \rbrace }(\varGamma )\big ). \end{aligned}$$ -
b)
Let the discrete initial values satisfy
$$\begin{aligned} \left\| \textbf{u}^{\,0}_h- I_h\textbf{u}^{\,0}\right\| _{\mathbb {H}^1} + \left\| \textbf{v}^{\,0}_h- I_h\textbf{v}^{\,0}\right\| _{\mathbb {H}^0} \le C_{\text {iv}} h^p \end{aligned}$$with a constant \(C_{\text {iv}}\) independent of h.
As main theorem, we state the following error bound for the finite element discretization of the wave equation with nonlinear acoustic boundary conditions.
Theorem 2.6
Let Assumption 2.2 be satisfied and \(\textbf{u}= \left[ u,\delta \right] ^\intercal \) be the solution of (1) on [0, T]. Further, let Assumption 2.5 be satisfied and let \(\textbf{u}_h= \left[ u_h,\delta _h \right] ^\intercal \) be the spatial approximation of \(\textbf{u}\), obtained with the bulk-surface finite element method of order p. Then, there exists some \(h_0>0\) s.t. the error bound
holds true for all \(h<h_0\), where \(c' = \frac{1}{\mu } \Big ( \frac{L_{\eta }^2c_\varOmega }{4\theta _0} - d \Big )\) and C is a constant independent of h.
In the next two sections we will now present a general theory for the error analysis of non-conforming space discretizations which we then use to proof Theorem 2.6 in Sect. 5.
3 Abstract space discretizations of first-order evolution equations with monotone operators
In this section we present the unified error analysis for abstract space discretizations of first-order evolution equations with maximal monotone operators. This generalizes the results from [16] and [18] for linear and semilinear equations, respectively. The results of this section are part of the dissertation [20].
We first present the continuous equation and the corresponding abstract space discretization, before we prove an error bound.
3.1 Analytical setting
Let X be a Hilbert space with scalar product \({\big (}{\cdot , \cdot }{\big )_{X}}\) in which we consider the evolution equation
In the following, we omit the t arguments in evolution equations. We pose the following classical assumptions to ensure that (10) is wellposed.
Assumption 3.1
-
a)
The nonlinear operator \(\mathcal {S}:D(\mathcal {S})\rightarrow X\) is quasi-monotone and maximal, i.e., there is a \(c_{\text {qm}}> 0\) s.t.
$$\begin{aligned} {\big (}{\mathcal {S}(y)-\mathcal {S}(z), y-z}{\big )_{X}} \ge -c_{\text {qm}}\left\| y-z\right\| _{X}^{2}\qquad \text {for all }y,z \in D(\mathcal {S}), \end{aligned}$$and there exists some \(\lambda >c_{\text {qm}}\) s.t.. \({\text {range}}(\lambda + \mathcal {S}) = X\).
-
b)
The inhomogeneity satisfies \(g\in W_{loc}^{1,1}([0,\infty );X)\).
The following wellposedness result can, e.g., be found in [25, Corollary IV.4.1].
Theorem 3.2
Let Assumption 3.1 hold true. Then, the evolution equation (10) is globally wellposed, i.e., (10) has a unique strong solution \(x\in C([0,\infty );X)\) which satisfies \(x(t)\in D(\mathcal {S})\) for all \(t\in [0,\infty )\), \(x(0)=x^0\), and (10a) is satisfied for almost all \(t\in [0,\infty )\).
We further state the following stability result which is essential for the latter error analysis.
Theorem 3.3
Let Assumption 3.1 be satisfied and for \(T>0\) and \(i=1,2\) let \(x_i\) be the strong solutions of
with \(g_i \in W^{1,1}([0,T];X)\). Then for all \(t \in [0,T]\)
Proof
The result can be derived with energy estimates similar to [25, Theorem IV.4.1A]. \(\square \)
3.2 Abstract space discretization
We now present an abstract space discretization of the evolution equation (10). Let \((X_h)_h\) be a family of finite dimensional vector spaces with scalar products \({\big (}{\cdot , \cdot }{\big )_{X_h}}\), where h is a discretization parameter, e.g., the maximal mesh width of a finite element discretization. For all \(X_h\in (X_h)_h\) we seek an approximations \(x_h\in X_h\) to the solution x of (10). Therefore, let \(\mathcal {S}_h\) and \(g_h\) be approximations of \(\mathcal {S}\) and \(g\), respectively, which satisfy the following assumptions similar to Assumption 3.1.
Assumption 3.4
-
a)
The nonlinear operator \(\mathcal {S}_h:X_h\rightarrow X_h\) is quasi-monotone, i.e., there is a \(\widehat{c}_{\text {qm}}> 0\) independent of h s.t.
$$\begin{aligned} {\big (}{\mathcal {S}_h(y_h)-\mathcal {S}_h(z_h), y_h-z_h}{\big )_{X_h}} \ge -\widehat{c}_{\text {qm}}\Vert {y_h-z_h}\Vert ^2_X\qquad \text {for all }y_h,z_h\in X_h. \end{aligned}$$(11) -
b)
The inhomogeneity satisfies \(g_h\in W_{loc}^{1,1}([0;\infty );X_h)\).
The discretized evolution equation is then given by
Since these assumptions are similar to the continuous case, we obtain by Theorem 3.2 that (12) is globally wellposed.
In the following we introduce a framework for the error analysis of the abstract space discretization that is similar to the linear case presented in [16]. To cover non-conforming space discretizations where \(X_h\nsubseteq X\), as they appear in Sect. 2, we make the following assumptions to relate the discrete and the continuous problem.
Assumption 3.5
-
a)
There exists a lift operator \({\mathcal {L}}_h\in \mathcal {L}(X_h,X)\) which satisfies
$$\begin{aligned} {}{\Vert }{{\mathcal {L}}_hy_h}{\Vert }{_X} \le \widehat{C}_X\Vert {y_h}\Vert X_h\qquad \text {for all }y_h\in X_h\end{aligned}$$(13)for some constant \(\widehat{C}_X>0\) independent of h. The adjoint of the lift operator \({{\mathcal {L}}}_h^*\in \mathcal {L}(X,X_h)\) is defined via
$$\begin{aligned} {\big (}{{{\mathcal {L}}}_h^*y, y_h}{\big )_{X_h}} = {\big (}{y, {\mathcal {L}}_hy_h}{\big )_{X}}, \qquad \text {for all }y\in X, y_h\in X_h. \end{aligned}$$ -
b)
Let \(Z \hookrightarrow X\) be a densely embedded subspace of X on which a reference operator \(J_h\in \mathcal {L}(Z;X_h)\) is defined which satisfies
$$\begin{aligned} {}{\Vert }{J_h}{\Vert }{_{X_h \leftarrow Z}} \le \widehat{C}_{J_h}\end{aligned}$$for some constant \(\widehat{C}_{J_h}>0\) independent of h.
As the term \({\mathcal {L}}_hJ_h-{\text {I}}\) appears in the error bound in Theorem 3.7, the reference operator should be chosen such that the approximation \({\mathcal {L}}_hJ_hz \approx z\) (\(z\in Z\)) is of the order of convergence which one wants to proof for the spatial discretization. It could, e.g., be an interpolation or a projection operator. One possible example for a lift operator is the lift defined in (6). We will consider this in Sect. 5.
The space discretization error bound is given in terms of the following terms:
Definition 3.6
(Remainder and error terms)
-
a)
The remainder of the nonlinear monotone operator is given by
$$\begin{aligned} R_h:D(\mathcal {S})\cap Z \rightarrow X_h,\qquad R_h(z){:}{=}{{\mathcal {L}}}_h^*\mathcal {S}(z) - \mathcal {S}_h\left( J_hz \right) . \end{aligned}$$(14) -
b)
We define the error term
$$\begin{aligned} \begin{aligned} E_h(t) =\;&\Vert x_{h}^{0}-J_{h} X^{0}\Vert x_{h}(+) t{}{\big \Vert }{\mathrm e^{-\widehat{c}_{\text {qm}}\cdot }({{\mathcal {L}}}_h^*- J_h)x'}{\big \Vert }{_{L^\infty \left( [0,t];X_h \right) }} \\&+ t{}{\Vert }{\mathrm e^{-\widehat{c}_{\text {qm}}\cdot }R_h(x)}{\Vert }{_{L^\infty \left( [0,t];X_h \right) }} + t{}{\Vert }{\mathrm e^{-\widehat{c}_{\text {qm}}\cdot }{{\mathcal {L}}}_h^*g- g_h}{\Vert }{_{L^\infty \left( [0,t];X_h \right) }}. \end{aligned} \end{aligned}$$(15)
We now can state and prove an error bound of the abstract space discretization, cf.[20, Thm. 2.10].
Theorem 3.7
Let Assumptions 3.1, 3.4, and 3.5 be satisfied and x be the strong solution of (10) on [0, T] with \(x,x'\in L^\infty ([0,T];Z)\). Furthermore, let \(x_h\) be the solution of (12) on [0, T]. Then, for all \(t\in [0,T]\) the lifted discrete solution satisfies the error bound
Proof
We split the error via \({\mathcal {L}}_hx_h(t)- x(t) = {\mathcal {L}}_he_h+ ({\mathcal {L}}_hJ_h- {\text {I}})x(t)\), where
is the discrete error. The full error can thus be bounded by
and we further investigate the discrete error. By applying the adjoint lift to (10a) we obtain
Adding \(J_hx', \mathcal {S}_h(J_hx),\) and \(g_h\) on both sides yields
where
Under Assumption 3.4, the stability estimate from Theorem 3.3 holds also true in the discrete case with \(\widehat{c}_{\text {qm}}\) instead of \(c_{\text {qm}}\). Hence, we obtain by Theorem 3.3 applied to (12) and (18) the following bound for the discrete error
where we used (19) and (14). Together with (17), we finally obtain (16).\(\square \)
In the following section we will use this result to derive error bounds for second-order nonlinear wave-type equations.
4 Abstract space discretizations of second-order evolution equations with nonlinear damping
In this section we apply the theory of Sect. 3 to second-order evolution equations. As in the previous section, we first introduce the continuous problem and then present and analyze the abstract space discretization. This is a generalization of the linear unified error analysis introduced in [16] and also an extension of the framework considered in the dissertation [20] which does not cover the acoustic boundary conditions with nonlinear coupling from Sect. 2, cf.Remark 4.2 and Sect. 5.
4.1 Analytical setting
Let V, H be Hilbert spaces es and let V be densely embedded in H. We consider the following variational equation, which is typical for a weak formulation of a second-order partial differential equation. Seek \(u \in C^2([0,T];H)\cap C^1([0,T];V)\) with
To ensure the wellposedness of (21) we pose the following assumptions.
Assumption 4.1
-
a)
The bilinear form \(m:H\times H \rightarrow \mathbb {R}\) is a scalar product on H with induced complete norm \({}{\Vert }{\cdot }{\Vert }{_{m}}{\cdot }\) In the following, we equip H with \(m\).
-
b)
The bilinear form \(a:V \times V \rightarrow \mathbb {R}\) is symmetric and there exists a constant \(c_G\ge 0\) s.t.
$$\begin{aligned} {{\tilde{a}}}{:}{=}a+ c_Gm\end{aligned}$$is a scalar product on V with induced complete norm \({}{\Vert }{\cdot }{\Vert }{_{{{\tilde{a}}}}}\). From now on, we equip V with \({{\tilde{a}}}\).
-
c)
The nonlinearity \({\mathcal {D}}\in C(V;V^*)\) satisfies \({\mathcal {D}}(0) = 0\) and is quasi-monotone, i.e., there is a constant \(\beta _{\text {qm}}\ge 0\) s.t.
$$\begin{aligned} {\langle }{{\mathcal {D}}(v)-{\mathcal {D}}(w),v-w}{\rangle _{V^{*}\times V}} \ge -\beta _{\text {qm}}\Vert v-w\Vert _m\qquad \text {for all }v,w \in V. \end{aligned}$$ -
d)
The inhomogeneity satisfies \(f\in W_{loc}^{1,1}([0,\infty );H)\).
We denote by \(C_{H,V}\) the embedding constant of V into H, i.e.,
Formulation as evolution equation We identify H with its dual space \( H^*\) to obtain the Gelfand triple
with dense embeddings. We thus have for all \(v\in V, w\in H\)
To reformulate (21) as an evolution equation, we define the operator \({\mathcal {A}}\in \mathcal {L}(V,V^*)\) associated to \(a\) via
Then, we can rewrite (21) equivalently as an evolution equation in \(V^*\): Seek \(u \in C^2([0,T];H)\cap C^1([0,T];V)\) satisfying
Note that (25) implicitly contains the condition
due to \(u'',f\in H\).
Remark 4.2
In [20], the stricter assumption \({\mathcal {D}}\in C(V;H)\) was posed. However, this does not cover the acoustic boundary conditions with nonlinear coupling (1c) as we will see in Sect. 5, cf.Remark 5.2.
First-order formulation We rewrite (25) into a first-order formulation in the framework of Sect. 3.1. For this let \(u' = v\) and we define
with
Then, (25) is equivalent to the first-order evolution equation (10).
In the following we show that the assumptions of Sect. 3.1 are satisfied. The subsequent lemma is a slight extension of [20, Lemma 2.14].
Lemma 4.3
The nonlinear operator \(\mathcal {S}\) is maximal and quasi-monotone with constant
and \(D(\mathcal {S})\) is dense in X.
Proof
We start by proving the quasi-monotonicity. For \(x_1 = \left[ u_1, v_1 \right] ^\intercal , x_2 = \left[ u_2, v_2 \right] ^\intercal \in D(\mathcal {S})\) we calculate by using Assumption 4.1, (23), and the definitions of \(\mathcal {S}\) and \({\mathcal {A}}\)
In the next step we prove the maximality and proceed similar as in the proof of [27, Theorem 4.1]. We have to show that there exists a \(\lambda > 0\) such that for every \(h = \left[ h_1,h_2 \right] ^\intercal \in X = V \times H\) there exists a solution \(x = \left[ v, w \right] ^\intercal \in D(\mathcal {S})\) of the stationary problem \((\lambda + \mathcal {S})x = h\) or equivalently
By solving (28a) for v and plugging it into (28b) we obtain
We thus investigate the operator \(T = \lambda + \frac{1}{\lambda }{\mathcal {A}}+ {\mathcal {D}}\in C(V;V^*)\) which can be decomposed via \(T = T_1 + T_2\) with
For
we then have that T is monotone as the sum of monotone operators. Further, we have for all \(v\in V\)
where we used that \(T_1\) is coercive due to the choice of \(\lambda \), and \(T_2\) is monotone with \(T_2(0) = 0\). Thus, T is coercive, i.e.
We apply [3, Corollary 2.3] stating that continuous, monotone, and coercive operators from a reflexive Banach space to its dual space are surjective. This yields the existence of a solution \(v\in V\) of (29) and thus also of a solution \(x = \left[ v, w \right] ^\intercal \in V\times V\) of (28). We further obtain by (28b) \(x\in D(\mathcal {S})\) since
The density of \(D(\mathcal {S})\) in X follows from the maximality and the quasi-monotonicity of \(\mathcal {S}\) and \(\mathcal {S}(0)=0\), cf [25, Prop. I.4.2].\(\square \)
Corollary 4.4
Assumption 4.1 implies that the first-order formulation of (25) satisfies Assumption 3.1.
Proof
By Lemma 4.3 we have that Assumption 3.1 a) is satisfied. Assumption 3.1 b) is directly implied by Assumption 4.1 d).\(\square \)
By Theorem 3.2 we then directly obtain the wellposedness of (21).
Corollary 4.5
Let Assumption 4.1 hold true and let \(\left[ u^0,v^0 \right] ^\intercal \in D(\mathcal {S})\), i.e., \(u^0,v^0\in V\) with \({\mathcal {A}}u^0+ {\mathcal {D}}(v^0) \in H\). Then, (21) is globally wellposed, i.e., there exists a unique strong solution \(\left[ u,v \right] ^\intercal \in C([0,\infty );V\times H)\).
4.2 Space discretization
We consider a family \((V_h)_h\) of finite dimensional vector spaces related to a discretization parameter h and the following discretized version of (21) in \(V_h\in (V_h)_h\): seek \(u_h\in C^2([0,T];V_h)\) with
Here, \(m_{h}, a_h, {\mathcal {D}}_h,\) and \(f_h\) are approximations of the corresponding continuous counterparts.
We pose the following assumptions similar to Assumption 4.1.
Assumption 4.6
All constants in the following statements are independent of h.
-
a)
The bilinear form \(m_{h}\) is a scalar product on \(V_h\). We denote \(V_h\) equipped with this scalar product \(m_{h}\) by \(H_h\) and the induced norm by \({}{\Vert }{\cdot }{\Vert }{_{m_{h}}}\).
-
b)
The bilinear form \(a_h:V_h\times V_h\rightarrow \mathbb {R}\) is symmetric and there exists a constant \(\widehat{c}_{G}\ge 0\) s.t.
$$\begin{aligned} {{\tilde{a}}}_h{:}{=}a_h+ \widehat{c}_{G}m_{h}\end{aligned}$$is a scalar product on \(V_h\) with induced norm \({}{\Vert }{\cdot }{\Vert }{_{{{\tilde{a}}}_h}}\). In the following, we equip \(V_h\) with \({{\tilde{a}}}_h\).
-
c)
The nonlinearity \({\mathcal {D}}_h\in C(V_h;H_h)\) satisfies \({\mathcal {D}}_h(0)=0\) and is continuous and quasi-monotone with constant \(\widehat{\beta }_{\text {qm}}\).
-
d)
The inhomogeneity satisfies \(f_h\in W_{loc}^{1,1}([0,\infty );H_h)\).
-
e)
There exists a constant \(\widehat{C}_{H,V}>0\) s.t.
$$\begin{aligned} {}{\Vert }{v_h}{\Vert }{_{m_{h}}}\le \widehat{C}_{H,V}{}{\Vert }{v_h}{\Vert }{_{{{\tilde{a}}}_h}} \qquad \text {for all }v_h\in V_h. \end{aligned}$$(31)
The operator \({\mathcal {A}}_h\in \mathcal {L}(V_h;V_h)\) related to \(a_h\) is defined via
We then can reformulate (30) as an evolution equation in \(V_h\):
Analogously to the continuous equation we rewrite (32) in a first-order formulation and therefore define \(X_h= V_h\times H_h\). With
(32) is then of the form (12).
Corollary 4.7
Assumption 4.6 implies that the first-order formulation of (32) satisfies Assumption 3.4. Furthermore, (11) holds true with \(\widehat{c}_{\text {qm}}= \frac{1}{2} \widehat{c}_{G}\widehat{C}_{H,V}+ \widehat{\beta }_{\text {qm}}\).
Proof
Since the setting in the discrete case from Assumption 4.6 is similar to the continuous one from Assumption 4.1 with constants independent of h, the proof of Lemma 4.3 transfers directly to the discrete case. \(\square \)
Similar to the first-order case, we require the existence of suitable operators to relate continuous and discrete functions of the abstract non-conforming space discretization.
Assumption 4.8
-
a)
There exists a lift operator \({{\mathcal {L}}_h^V\in \mathcal {L}(V_h;V)}\) satisfying
$$\begin{aligned} \begin{aligned} \Vert {\mathcal {L}}_h^Vv_h\Vert _m \le \widehat{C}_H{}{\Vert }{v_h}{\Vert }{_{m_{h}}}, \qquad {}{\Vert }{{\mathcal {L}}_h^Vv_h}{\Vert }{_{{{\tilde{a}}}}} \le \widehat{C}_V{}{\Vert }{v_h}{\Vert }{_{{{\tilde{a}}}_h}}, \end{aligned} \end{aligned}$$(34)for all \(v_h\in V_h\) with constants \(\widehat{C}_H, \widehat{C}_V> 0\) independent of h.
-
b)
There exists an interpolation operator \(I_h\in \mathcal {L}(Z^V;V_h)\), defined on a dense subspace \(Z^V\) of V, which satisfies
$$\begin{aligned} {}{\Vert }{I_h}{\Vert }{_{H_h\leftarrow Z^V}} \le \widehat{C}_{I_h}\end{aligned}$$(35)with a constant \(\widehat{C}_{I_h}> 0\) independent of h.
To apply the results of Sect. 3.2, we now define the first-order reference and lift operator.
Definition 4.9
-
a)
The adjoint lift operators \({\mathcal {L}}_h^{V* }:V \rightarrow V_h\) and \({\mathcal {L}}_h^{H*}:H \rightarrow H_h\) w.r.t. the scalar products of V and H are defined via
$$\begin{aligned} \begin{aligned} {m_{h}\big (}{{\mathcal {L}}_h^{H*}v, w_h}{\big )}&{:}{=}{m\big (}{v,{\mathcal {L}}_h^Vw_h}{\big )}{} & {} \text {for all } v \in H, w_h\in H_h,\\ {{{\tilde{a}}}_h\big (}{{\mathcal {L}}_h^{V* }v, w_h}{\big )}&{:}{=}{{{\tilde{a}}}\big (}{v, {\mathcal {L}}_h^Vw_h}{\big )}{} & {} \text {for all } v \in V, w_h\in V_h. \end{aligned} \end{aligned}$$(36) -
b)
We define the first-order lift operator \({\mathcal {L}}_h:X_h\rightarrow X\) by
$$\begin{aligned} {\mathcal {L}}_h\begin{bmatrix} v_h\\ w_h \end{bmatrix} {:}{=}\begin{bmatrix} {\mathcal {L}}_h^Vv_h\\ {\mathcal {L}}_h^Vw_h \end{bmatrix}. \end{aligned}$$ -
c)
We define the first-order reference operator \({J_h:Z \rightarrow X_h}\) by
$$\begin{aligned} J_h\begin{bmatrix} v \\ w \end{bmatrix} {:}{=}\begin{bmatrix} {\mathcal {L}}_h^{V* }v \\ I_hw \end{bmatrix} \end{aligned}$$(37)on \(Z = V \times Z^V\overset{d}{\hookrightarrow }X\).
Lemma 4.10
The first-order lift and reference operators from Definition 4.9 satisfy Assumption 3.5 with \(\widehat{C}_X= \max \lbrace \widehat{C}_V, \widehat{C}_H\rbrace \) and \(\widehat{C}_{J_h}= \max \lbrace \widehat{C}_V,\widehat{C}_{I_h}\rbrace \).
Proof
This is a direct consequence of Assumption 4.8. \(\square \)
In the following we now bound the first-order remainder term which is for \(z=\left[ v,w \right] ^\intercal \in D(\mathcal {S}) \cap Z\) given by
To do so, we use the following error terms in the scalar products, which are for \(v_h, w_h\in V_h\) defined via
We obtain the following bound for the remainder term, cf.[20, Lem. 2.23]
Lemma 4.11
Let Assumption 4.1 and 4.6 be satisfied. Then, for \(z=\left[ v,w \right] ^\intercal \in D(\mathcal {S}) \cap Z\), the remainder of the monotone operator can be bounded by
i.e., against errors in the scalar products, interpolation errors, and the discretization error of the nonlinear operator.
Proof
The proof works similar to the proof of [16, Lemma 4.7] and relies on the identity
where \({\big (}{\cdot , \cdot }{\big )_{X_h}}\) is the scalar product on \(X_h\). Thus, let \(y_h= \left[ \varphi _h, \psi _h \right] ^\intercal \in X_h\) with \(\Vert y_h\Vert X_h=1\). By (38) we obtain
and we bound the first two summands separately. To bound the first one, we use (39), (34), and \({}{\Vert }{\varphi _h}{\Vert }{_{{{\tilde{a}}}_h}}\le 1\) to obtain
By using the definitions of \({{\tilde{a}}},{{\tilde{a}}}_h\), \({}{\Vert }{\psi _h}{\Vert }{_{m_{h}}}\le 1\) and (39), (22), (34), (31), we bound the second summand in (41) via
Similar to (42), we further estimate
We finally obtain the assertion by collecting all terms. \(\square \)
We are now in the position to prove the following error bound which is a generalization of [20, Thm. 2.24], cf. Remark 4.2. It is applicable to all equations and space discretizations fitting in the abstract framework of this section. In the dissertation [20], it was used to prove novel convergence rates for the wave equation with nonlinear kinetic boundary conditions and in this paper, we use it to prove Theorem 2.6.
Theorem 4.12
Let Assumptions 4.1, 4.6, and 4.8 be satisfied and u be the strong solution of (25) on [0, T] with \(u, u', u'' \in L^\infty ([0,T];Z^V)\). Further, let \(u_h\) be the semidiscrete solution of (32) on [0, T]. Then, for all \(t \in [0,T]\), the lifted semidiscrete solution satisfies the error bound
with a constant C that is independent of h and t. The other constants are given by
and the abstract space discretization errors
Proof
By Corollaries 4.4, 4.7, and Lemma 4.10, we have that the first-order formulations of (25) and (32) satisfy all assumptions of Theorem 3.7.
By applying Theorem 3.7 and employing the error bound (16), we obtain
with
In the remaining proof, we bound the different terms against \(E_{h,i}, i = 1, \ldots , 4\). For the remainder term we apply the bound (40) and obtain for all \(t\in [0,T]\)
By the definitions of \(J_h\) and \({{\mathcal {L}}}_h^*\) we further have for the discretization errors of the initial values and the inhomogeneity
The reference error can be decomposed for all \(t\in [0,T]\) via
where we have similar to (42)
In the same way, we finally bound
\(\square \)
Having this abstract theory at hand, we can now return to the wave equation with nonlinear acoustic boundary conditions from Sect. 2 and give the proof of Theorem 2.6 in the next section.
5 Numerical analysis of wave equations with nonlinear acoustic boundary conditions
In this section we will use the unified error analysis for second-order equations from Sect. 4 to prove the error bound from Theorem 2.6. We start by verifying that all assumptions are satisfied.
Lemma 5.1
Let Assumption 2.2 be satisfied. Then, with the definitions in (5), Assumption 4.1 is satisfied with \(\beta _{\text {qm}}= \frac{1}{\mu } \Big ( \frac{L_{\eta }^2c_\varOmega }{4\theta _0} - d \Big ), c_G= 1,\) and \(C_{H,V}= 1\).
Proof
We clearly have that \(m\) is a scalar product on H and that \({{\tilde{a}}}{:}{=}a+ m\) is a scalar product on V. Further, Assumption 2.2 d) directly implies Assumption 4.1 d).
Thus it remains to prove Assumption 4.1 c). By Assumption 2.2 a), b), c) and (5a), (5c) we obtain
This proves the the quasi-monotonicity of \({\mathcal {D}}\).
In the next step we show \({\mathcal {D}}\in C(V;V^*)\). We emphasize that the trace inequality
holds true for \(q = \zeta +1\) with \(\zeta \) from the growth condition (3), cf.[24, Thms. 2.4.2 and 2.4.6]. For \(\textbf{v}_1 = \left[ v_1,z_1 \right] ^\intercal , \textbf{v}_2 = \left[ v_2,z_2 \right] ^\intercal , {\varvec{\varphi }}= \left[ \varphi , \psi \right] ^\intercal \in V\) with \({}{\Vert }{{\varvec{\varphi }}}{\Vert }{_{{{\tilde{a}}}}} = 1\) this yields together with the Hölder and the Minkowski inequalities and the global Lipschitz continuity of \(\eta \)
We hence obtain
By the trace inequality (45), the growth condition (3), the relation \(\zeta = q-1\), and [12, Theorem 2] we further have \(v \mapsto \theta (v) \in C(H^1(\varOmega );L^\frac{q}{q-1}(\varGamma ))\). This yields
which proves \({\mathcal {D}}\in C(V, V^*)\). \(\square \)
Remark 5.2
It is not possible to prove the stronger condition \({\mathcal {D}}\in C(V,H)\). This is due to the fact that the calculation (46) strongly relies on \(\varphi \in H^1(\varOmega )\) and is not possible for a test function \(\varphi \in L^2(\varOmega )\).
Lemma 5.1 ensures that the weak formulation (4) of (1) fits in the setting of Sect. 4.1 and, hence, is locally wellposed by Corollary 4.4.
We now prove, that the bulk-surface finite element space discretization from Sect. 2.2 fits into the abstract setting of Sect. 4.2.
Lemma 5.3
Let Assumption 2.2 hold true. Then, the bulk-surface finite element space discretization of (1) satisfies Assumption 4.6 with \(\widehat{\beta }_{\text {qm}}= \frac{1}{\mu } \Big ( \frac{L_{\eta }^2 c_\varOmega }{4\theta _0} - d \Big ), \widehat{c}_{G}= 1,\) and \(\widehat{C}_{H,V}= 1\).
Proof
Since \(a_h\) and \(m_{h}\) are defined as in continuous case, Assumption 4.6 a) and b) are satisfied. Assumption 4.6 d) follows from Assumption 2.2 d) and the continuity of the interpolation operator.
It remains to prove Assumption 4.6 c). To show the quasi-monotonicity, we proceed analogously to the proof in the continuous case from Lemma 5.1 and obtain by (8c)
where we used that the quadrature formula has positive weights and satisfies (7).
Finally, \({\mathcal {D}}_h\) is continuous, since \(V_h\) is a finite dimensional space and, thus, convergence in \(V_h\) implies uniform pointwise convergence and especially convergence in all quadrature nodes. \(\square \)
To prove an error bound for the semidiscretization, we apply the theory of Sect. 4.2 and therefore have to specify the operators from Assumption 4.8.
Definition 5.4
-
a)
The lift operator \({\mathcal {L}}_h^V\in \mathcal {L}(V_h;V)\) is defined via
$$\begin{aligned} {\mathcal {L}}_h^V\left[ v_h,z_h \right] ^\intercal {:}{=}\left[ v_h^\ell ,z_h^\ell \right] ^\intercal \qquad \text {for all }\left[ v_h,z_h \right] ^\intercal \in V_h\end{aligned}$$with \(v_h^\ell \) from (6).
-
b)
We set \(Z^V{:}{=}H^2(\varOmega )\times H^2(\varGamma )\subset C({\overline{\varOmega }})\times C(\varGamma )\).
-
c)
We define the interpolation operator via \(I_h\left[ v_h,z_h \right] ^\intercal {:}{=}\left[ I_{h,\varOmega }v_h,I_{h,\varGamma }z_h \right] ^\intercal \).
Our error analysis relies on the following properties of the lift and the interpolation operators.
First of all, there exist element-wise norm equivalences related to the lift, which were shown in [8, Lemmas 5.3 and 7.3].
Lemma 5.5
There exist \(C_{\varOmega ,\varOmega _h}> c_{\varOmega ,\varOmega _h}>0\), \(C_{\varGamma ,\varGamma _h}>c_{\varGamma ,\varGamma _h}>0\) independent of \(h<h_0\) sufficiently small s.t.for all \(v_h\in V_{h,p}^\varOmega \), \(\vartheta _h\in V_{h,p}^\varGamma \), \(k=0,1,\ldots ,p+1\), and \(K_\varOmega \in {\mathcal {T}}_h\), \(K_\varGamma \in {\mathcal {T}}_h^\varGamma \) and for all \(h<h_0\) sufficiently small we have
where \(K_\varOmega ^\ell = G_h(K_\varOmega ), K_\varGamma ^\ell = G_h(K_\varGamma )\). By construction, the lift additionally preserves the \(L^\infty \) norm, i.e.,
Further, we have the following bounds of the geometric errors stemming from the domain approximation (cf.[7, proof of Lemma 6.2]).
Lemma 5.6
For \(u_h, \varphi _h \in V_{h,p}^\varOmega \) and \(\vartheta _h, \psi _h\in V_{h,p}^\varGamma \), the following bounds hold true for all \(h<h_0\) sufficiently small:
The nodal interpolation satisfy the following error bounds, which follow from [8, Theorem 4.28, Theorem 5.9] for the bulk and [8, Theorem 6.24, Theorem 7.10] for the surface interpolation, respectively.
Lemma 5.7
Let \(1\le k \le p\) and \(h<h_0\) sufficiently small.
-
a)
Globally, the interpolation operators satisfy for all \(v \in H^{k+1}(\varOmega )\), and \(\vartheta \in H^{k+1}(\varGamma )\) the error bounds
$$\begin{aligned} \left\| v-(I_{h,\varOmega }v)^\ell \right\| _{L^2(\varOmega )} + h\left\| v-(I_{h,\varOmega }v)^\ell \right\| _{H^1(\varOmega )}&\le C h^{k+1}\left\| v\right\| _{H^{k+1}(\varOmega )} , \end{aligned}$$(49a)$$\begin{aligned} \left\| \vartheta -(I_{h,\varGamma }\vartheta )^\ell \right\| _{L^2(\varGamma )} + h\left\| \vartheta -(I_{h,\varGamma }\vartheta )^\ell \right\| _{H^1(\varGamma )}&\le C h^{k+1}\left\| \vartheta \right\| _{H^{k+1}(\varGamma )} , \end{aligned}$$(49b)with a constant C independent of h.
-
b)
Locally, on each element \(K_\varOmega \in {\mathcal {T}}_h^\varOmega \), \(K_\varGamma \in {\mathcal {T}}_h^\varGamma \), the interpolation operators satisfy for all \(0\le r\le k\) and all \(v \in H^{k+1}(K_\varOmega ^\ell ), \vartheta \in H^{k+1}(K_\varGamma ^\ell )\), the error bounds
$$\begin{aligned} \left\| v-(I_{h,\varOmega }v)^\ell \right\| _{H^r(K_\varOmega ^\ell )}&\le C h^{k+1-r}\left\| v\right\| _{H^{k+1}(K_\varOmega ^\ell )} , \end{aligned}$$(50a)$$\begin{aligned} \left\| \vartheta -(I_{h,\varGamma }\vartheta )^\ell \right\| _{H^r(K_\varGamma ^\ell )}&\le C h^{k+1-r}\left\| \vartheta \right\| _{H^{k+1}(K_\varGamma ^\ell )} , \end{aligned}$$(50b)with a constant C independent of h.
-
c)
Locally, on each element \(K_\varOmega \in {\mathcal {T}}_h^\varOmega \), \(K_\varGamma \in {\mathcal {T}}_h^\varGamma \), and for every \(v_h\in H^{k+1}(K_\varOmega ), \vartheta _h\in H^{k+1}(K_\varGamma )\), the \(L^\infty \) error bounds
$$\begin{aligned} \left\| v_h-I_{h,\varOmega }v_h^\ell \right\| _{L^\infty (K_\varOmega )}&\le C h^{k+1}\left\| v_h\right\| _{W^{k+1,\infty }(K_\varOmega )} , \end{aligned}$$(51a)$$\begin{aligned} \left\| \vartheta _h-I_{h,\varGamma }\vartheta _h^\ell \right\| _{L^\infty (K_\varGamma )}&\le C h^{k+1}\left\| \vartheta _h\right\| _{W^{k+1,\infty }(K_\varGamma )} \end{aligned}$$(51b)hold true with a constant C independent of h.
The following lemma is a direct consequence of Lemmas 5.6 and 5.7.
Lemma 5.8
The operators defined in Definition 5.4 satisfy Assumption 4.8 with
where \(C_{\varOmega ,\varOmega _h}\) and \(C_{\varGamma ,\varGamma _h}\) are given in (47).
We are now in the position to prove the error bound of the space discretization. Proof of Theorem 2.6. We apply Theorem 4.12. By Lemmas 5.1, 5.3, and 5.8 we have that all assumptions are satisfied and we have to bound the space discretization error terms \(E_{h,i}\) in (44). In the following we always assume that \(h<h_0\) is sufficiently small that all the bounds we apply are valid.
The terms \(E_{h,1},E_{h,3}\) and \(E_{h,4}\) also appeared in the linear case and were bounded under Assumption 2.5 in [16, Proof of Thm 5.3] by order \(h^p\).
It thus remains to bound \(E_{h,2}\). For \(t \in [0,T]\) and \(\textbf{v}= \textbf{u}'(t) = \left[ v,z \right] ^\intercal \in V\) we calculate
We bound the different terms separately. Let \({}{\Vert }{{\varvec{\varphi }}_h}{\Vert }{_{m_{h}}} = {}{\Vert }{\left[ \varphi _h, \psi _h \right] ^\intercal }{\Vert }{_{m_{h}}} = 1\). We then have
For the first summand on the right-hand side of (53) we have by the continuity of the lift operator and the interpolation error (49b)
The second summand can be bounded using the geometric error estimate (48c) by
To bound the third summand on the right-hand side of (53) we use that for the nodal interpolation we have \(I_{h,\varGamma }\theta (v) = I_{h,\varGamma }\theta ((I_{h,\varGamma }v)^\ell )\in V_{h,p}^\varGamma \) and that the order of quadrature formula is greater than 2p to obtain
where we denote by \(\sigma (\varGamma _h)\) the measure of \(\varGamma _h\) and we used the \(L^\infty \) interpolation error bound (51b). The terms \(\left\| \theta (I_{h,\varGamma }v)\right\| _{W^{ p,\infty }(F)}\) are bounded due since we have \(v \in W^{p+1,\infty }(\varOmega )\), \(v \mapsto v\big |_{{\varGamma }}\in C(W^{p+1,\infty }(\varOmega );W^{p,\infty }(\varGamma ))\), \(I_{h,\varGamma }\in C(W^{p,\infty }(\varGamma );W^{p,\infty }(\varGamma _h))\), and \(\theta \in C^{p}\), cf.(9b). In total we obtain in (53)
Similarly we obtain
but since z is defined on \(\varGamma \) we only require the regularity \(z\in W^{p,\infty }(\varGamma )\) in contrast to \(v \in W^{p+1,\infty }(\varOmega )\).
To bound the third term in (52), we make use of the classical inverse estimate \(\left\| v_h\right\| _{H^1(\varOmega _h)} \le Ch^{-1}\left\| v_h\right\| _{L^2(\varOmega _h)}\) (cf.[6, Lem. 4.5.3]), (49b), and the trace inequality \(\left\| v\right\| _{H^{p+1}(\varGamma )} \le C\left\| v\right\| _{H^{p+2}(\varOmega )}\) to obtain
Similarly, we have
but we do not need the inverse estimate in this case and hence the smoothness \(z\in H^{p+1}(\varGamma )\) is sufficient. In total we obtain \(E_{h,2}\le Ch^p\) and
Theorem 4.12 gives then the desired result. \(\square \)
6 Numerical experiment
In this section we illustrate Theorem 2.6 with a numerical experiment.
Let \(\varOmega = B_1(0) \subset \mathbb {R}^2\) be the unit disc and in (1) we set
where \(r = r(\textbf{x}) = \textbf{x}^2_1 + \textbf{x}_2^2\). Then, Assumption 2.2 is satisfied and the exact solution of (1) is given by
Error from (54) at \(t=0.7\) for the test example. The dashed plots are straight lines of slope 1, 2 and 3, respectively
We implemented the experiments in the C++ finite element library deal.ii, cf.[1, 2]. The code which was used for the numerical is available at https://doi.org/10.5445/IR/1000150159. For the time integration we use the implicit midpoint rule with time step size \(\approx 2\cdot 10^{-4}\) and solve the arising nonlinear systems with the simplified Newton method. For the spatial discretization we use the bulk-surface finite element methods of orders \(p = 1,2,3\).
We consider the error
instead of the error from Theorem 2.6 since the computation of the lift is quite laborious. We evaluated the integrals with a quadrature rule of degree 2p, so that the quadrature error is negligible. The restriction of u to \(\varOmega _h\) is possible for this example since we have \(\varOmega _h\subset \varOmega \).
In Fig. 1 the error \(\textbf{E}(t)\) is plotted against the maximal mesh width h. We chose \(t=0.7\) since it keeps distance to the roots of \(\sin (2\pi t)\). We observe that the error converges with order p as predicted by Theorem 2.6 which indicates that our proven convergence rates are optimal. Note that for order \(p=3\) and small h, the plot approaches the plateau of the time discretization error.
References
Arndt, D., Bangerth, W., Blais, B., Fehling, M., Gassmöller, R., Heister, T., Heltai, L., Köcher, U., Kronbichler, M., Maier, M., Munch, P., Pelteret, J.P., Proell, S., Simon, K., Turcksin, B., Wells, D., Zhang, J.: The deal.II library, version 9.3. J. Numer. Math. 28(3), 171–186 (2021)
Arndt, D., Bangerth, W., Davydov, D., Heister, T., Heltai, L., Kronbichler, M., Maier, M., Pelteret, J.P., Turcksin, B., Wells, D.: The deal.II finite element library: design, features, and insights. Comput. Math. Appl. 81, 407–422 (2021). https://doi.org/10.1016/j.camwa.2020.02.022
Barbu, V.: Nonlinear Differential Equations of Monotone Types in Banach Spaces. Springer Monographs in Mathematics, Springer, New York (2010). https://doi.org/10.1007/978-1-4419-5542-5
Beale, J.T.: Spectral properties of an acoustic boundary condition. Indiana Univ. Math. J. 25(9), 895–917 (1976). https://doi.org/10.1512/iumj.1976.25.25071
Beale, J.T., Rosencrans, S.I.: Acoustic boundary conditions. Bull. Am. Math. Soc. 80, 1276–1278 (1974). https://doi.org/10.1090/S0002-9904-1974-13714-6
Brenner, S.C., Scott, L.R.: The Mathematical Theory of Finite Element Methods. Texts in Applied Mathematics, vol. 15, 3rd edn. Springer, New York (2008). https://doi.org/10.1007/978-0-387-75934-0
Elliott, C.M., Ranner, T.: Finite element analysis for a coupled bulk-surface partial differential equation. IMA J. Numer. Anal. 33(2), 377–402 (2013). https://doi.org/10.1093/imanum/drs022
Elliott, C.M., Ranner, T.: A unified theory for continuous-in-time evolving finite element space approximations to partial differential equations in evolving domains. IMA J. Numer. Anal. (2020). https://doi.org/10.1093/imanum/draa062.Draa062
Emmrich, E., Šiška, D., Thalhammer, M.: On a full discretisation for nonlinear second-order evolution equations with monotone damping: construction, convergence, and error estimates. Found. Comput. Math. 15(6), 1653–1701 (2015). https://doi.org/10.1007/s10208-014-9238-4
Frota, C.L., Vicente, A.: Uniform stabilization of wave equation with localized internal damping and acoustic boundary condition with viscoelastic damping. Z. Angew. Math. Phys. 69(3), 24 (2018). https://doi.org/10.1007/s00033-018-0977-y
Gal, C.G., Goldstein, G.R., Goldstein, J.A.: Oscillatory boundary conditions for acoustic wave equations, pp. 623–635 (2003). https://doi.org/10.1007/s00028-003-0113-z. (Dedicated to Philippe Bénilan)
Goldberg, H., Kampowsky, W., Tröltzsch, F.: On Nemytskij operators in \(L_p\)-spaces of abstract functions. Math. Nachr. 155, 127–140 (1992). https://doi.org/10.1002/mana.19921550110
Graber, P.J.: Uniform boundary stabilization of a wave equation with nonlinear acoustic boundary conditions and nonlinear boundary damping. J. Evol. Equ. 12(1), 141–164 (2012)
Graber, P.J., Said-Houari, B.: On the wave equation with semilinear porous acoustic boundary conditions. J. Differ. Equ. 252(9), 4898–4941 (2012). https://doi.org/10.1016/j.jde.2012.01.042
Hipp, D.: A unified error analysis for spatial discretizations of wave-type equations with applications to dynamic boundary conditions. Ph.D. thesis, Karlsruher Institut für Technologie (KIT) (2017). https://doi.org/10.5445/IR/1000070952
Hipp, D., Hochbruck, M., Stohrer, C.: Unified error analysis for nonconforming space discretizations of wave-type equations. IMA J. Numer. Anal. 39(3), 1206–1245 (2019). https://doi.org/10.1093/imanum/dry036
Hipp, D., Kovács, B.: Finite element error analysis of wave equations with dynamic boundary conditions: \(L^2\) estimates. IMA J. Numer. Anal. 41(1), 683–728 (2021). https://doi.org/10.1093/imanum/drz073
Hochbruck, M., Leibold, J.: Finite element discretization of semilinear acoustic wave equations with kinetic boundary conditions. Electron. Trans. Numer. Anal. 53, 522–540 (2020). https://doi.org/10.1553/etna_vol53s522
Hochbruck, M., Maier, B.: Error analysis for space discretizations of quasilinear wave-type equations. IMA J. Numer. Anal. (2021). https://doi.org/10.1093/imanum/drab073.Drab073
Leibold, J.: A unified error analysis for the numerical solution of nonlinear wave-type equations with application to kinetic boundary conditions. Ph.D. thesis, Karlsruher Institut für Technologie (KIT) (2021). https://doi.org/10.5445/IR/1000130222
Ma, T.F., Souza, T.M.: Pullback dynamics of non-autonomous wave equations with acoustic boundary condition. Differ. Integr. Equ. 30(5–6), 443–462 (2017)
Maier, B.: Error analysis for space and time discretizations of quasilinear wave-type equations. Ph.D. thesis, Karlsruher Institut für Technologie (KIT) (2020). https://doi.org/10.5445/IR/1000120935. https://publikationen.bibliothek.kit.edu/1000120935
Mugnolo, D., Vitillaro, E.: The wave equation with acoustic boundary conditions on non-locally reacting surfaces. arXiv preprint arXiv:2105.09219 (2021)
Nečas, J.: Direct Methods in the Theory of Elliptic Equations. Springer Monographs in Mathematics. Springer, Heidelberg (2012). https://doi.org/10.1007/978-3-642-10455-8. (Translated from the 1967 French original by Gerard Tronel and Alois Kufner, Editorial coordination and preface by Šárka Nečasová and a contribution by Christian G. Simader)
Showalter, R.E.: Monotone Operators in Banach Space and Nonlinear Partial Differential Equations. Mathematical Surveys and Monographs, vol. 49. American Mathematical Society, Providence (1997)
Vicente, A., Frota, C.L.: On a mixed problem with a nonlinear acoustic boundary condition for a non-locally reacting boundaries. J. Math. Anal. Appl. 407(2), 328–338 (2013)
Vitillaro, E.: On the wave equation with hyperbolic dynamical boundary conditions, interior and boundary damping and source. Arch. Ration. Mech. Anal. 223(3), 1183–1237 (2017). https://doi.org/10.1007/s00205-016-1055-2
Wu, J.: Well-posedness for a variable-coefficient wave equation with nonlinear damped acoustic boundary conditions. Nonlinear Anal. 75(18), 6562–6569 (2012). https://doi.org/10.1016/j.na.2012.07.032
Acknowledgements
I thank Benjamin Dörich and Marlis Hochbruck for the careful reading of this manuscript and helpful discussions.
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Project-ID 258734477 - SFB 1173.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Leibold, J. A unified error analysis for nonlinear wave-type equations with application to acoustic boundary conditions. Numer. Math. 152, 907–936 (2022). https://doi.org/10.1007/s00211-022-01326-8
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00211-022-01326-8