Search for Articles:

Contents

ISSN: xxxx-xxxx (Online) xxxx-xxxx (Print)

\(L^{\infty}\)- error estimate for impulse control QVIs with vanishing control cost

Messaoud Boulbrachene1
1Department of Mathematics, Sultan Qaboos University. P.O. Box 36, Muscat 123, Oman
Copyright © Messaoud Boulbrachene. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

This paper deals with the finite element approximation of the elliptic impulse control quasi-variational inequality (QVI), when the impulse control cost goes to zero. By means of the concepts of subsolutions for QVIs and a Lipschitz dependence property with respect to the impulse cost, an \(L^{\infty}\) error estimate is derived for both the impulse control QVI and the correponding asymptotic problem.

Keywords: Quasi variational inequalities, impulse control cost, finite elements, subsolutions, \(L^{\infty}\)- error estimate

1. Introduction

This paper deals with the standard finite element approximation in the \(L^{\infty}\)– norm of the impulse quasi-variational inequality (QVI): find \(u_{k}\in H^{1}(\Omega)\) such that \[\begin{cases} a(u_{k},v-u_{k})\geqq(f,\text{ }v-u_{k})\;&\forall v\in H^{1}(\Omega ),\\ u_{k}\leq k+Mu_{k}\,,\,v\leq k+Mu_{k}\,,&u_{k}\geq0, \end{cases} \label{CQVI} \tag{1}\] as the impulse cost \(k\) goes to zero.

Here \(\Omega\) is a bounded convex domain of \(\mathbb{R}^{N},\) \(N\geq1,\) with smooth boundary, \((.,.)\) corresponds to the inner product in \(L^{2}(\Omega)\), \(\ f\geq0\) is a function in \(L^{\infty}(\Omega),\) \(a(.\,,.)\) is the bilinear form \[a(u\,,v)=\int_{\Omega}\nabla u\nabla v(x)+cuvdx\text{, }\label{BF} \tag{2}\] where \(c\) is a positive constant, and \(M\) is the nonlinear operator, defined from \(L^{\infty}\left( \Omega\right)\) into itself, by \[M\varphi(x)=\inf\varphi(x+\xi)\text{; }x\in\Omega\text{, }\xi\geq0\text{; }x+\xi\in\bar{\Omega},\label{IMP} \tag{3}\] and \(k\) is a positive constant.

Let \(\Omega\) be decomposed into triangles and let \(\tau_{h}\) denote the set of all those elements; \(h\) \(>0\) is the mesh size.We assume that the family \(\tau_{h}\) is regular and quasi-uniform; let also \(\mathbb{V}_{h}\) denote the finite element space consisting of continuous piecewise linear functions, \(\{\varphi_{s}\},\,s=1,2,…m(h)\) the basis of \(\mathbb{V}_{h}\), and \(r_{h}\) the usual corresponding Lagrange interpolation operator.

The discrete counterpart of (1) consists of seeking \(u_{k,h} \in\mathbb{V}_{h}\) such that \[\left\{ \begin{array} [c]{l} a(u_{k,h},v-u_{k,h})\geqq(f,v-u_{k,h})\;\forall v\in\mathbb{V}_{h},\\ v\leq k+r_{h}Mu_{k,h},\text{ }u_{k,h}\leq k+r_{h}Mu_{k,h}. \end{array} \right. \label{QVIH} \tag{4}\]

In this paper, under a \(W^{1,\infty}\)– uniform regularity assumption, we derive error estimate for the QVI (1) \[\left\Vert u_{k}-u_{k,h}\right\Vert _{\infty}\leq Ch,\] and for the asymptotic problem \[\left\Vert u_{0}-u_{0,h}\right\Vert _{\infty}\leq Ch,\label{MR0} \tag{5}\] where \(u_{0}\) is the solution of the continuous asymptotic QVI (as \(k\rightarrow0\)), that is \[\left\{ \begin{array} [c]{l} a(u_{0},v-u_{0})\geqq(f,\text{ }v-u_{0})\;\forall v\in H^{1}(\Omega ),\\ u_{0}\leq Mu_{0}\,,\,v\leq Mu_{0}\,\qquad ,\text{ }u_{0}\geq0, \end{array} \right. \label{ASYMP} \tag{6}\] (see [1]), \(u_{0,h}\) is the solution of the discrete asymptotic QVI (as \(k\rightarrow0\)), that is, \[\left\{ \begin{array} [c]{l} a(u_{0,h},v-u_{0,h})\geqq(f,\text{ }v-u_{0,h})\;\forall v\in H^{1} (\Omega),\\ u_{0,h}\leq r_{h}Mu_{0,h}\,,\,v\leq r_{h}Mu_{0,h}\text{, }u_{0,h} \geq0,\, \end{array} \right. \label{ASYMPH} \tag{7}\] and \(C\) is a constant independent of both \(h\) and \(k\).

To achieve this, we develop an approach that, in both the continuous and discrete settings, combines two key ingredients: a characterization of the QVI solution as the upper envelope of the set of subsolutions, and a Lipschitz continuity property with respect to the impulse cost.

It’s worth mentioning that the approximation approach developped in this paper is completely different from the one in [2].

The remainder of the paper is organized as follows. §2 recalls several qualitative properties of the solution to the continuous QVI (1), establishes its Lipschitz dependence with respect to the impulse cost, and introduces the continuous asymptotic QVI. §3 presents the corresponding discrete theory: assuming a discrete maximum principle, we derive analogous qualitative properties for the solution of the discrete QVI (4), and state the discrete asymptotic QVI. Finally, in §4, we give a detailed analysis of the approximation, and convergence order.

2. Preliminaries

2.1. Elliptic variational inequality (VI)

Let \(g\in L^{\infty}\left( \Omega\right)\) and \(\psi\in W^{1,\infty}\left( \Omega\right)\) such that \(\partial\psi/\partial n\geq0\) on \(\partial\Omega\). The following problem is called an elliptic variational inequality (VI): find \(\omega=\sigma(\psi)\in H^{1}(\Omega)\) such that \[\left\{ \begin{array} [c]{l} a(\omega,v-\omega)\geqq(g,v-\omega)\,\forall v\in H^{1}(\Omega),\\ v\leq\psi\,,\,\,\omega\leq\psi\,. \end{array} \right. \label{VI} \tag{8}\]

Theorem 1. [3] Let \(\psi\) and \(\tilde{\psi}\) in \(W^{1,\infty}(\Omega)\).Then, we have \[\left\Vert \omega-\tilde{\omega}\right\Vert _{\infty}\leq\left\Vert \psi-\tilde{\psi}\right\Vert _{\infty},\]

The discrete analog of VI (8) is the VI: find \(\omega_{h}=\sigma _{h}(\psi)\in\mathbb{V}_{h}\) such that \[\left\{ \begin{array} [c]{l} a(\omega_{h},v-\omega_{h})\geqq(g,v-\omega_{h})\,\,\,\forall v\in \mathbb{V}_{h},\\ v\leq r_{h}\psi\,,\,\,\omega_{h}\leq r_{h}\psi.\, \end{array} \right. \label{VIG} \tag{9}\]

We assume that the matrix with generic coefficients \(a\,(\varphi_{l} ,\varphi_{s})\,\,\)is an M-Matrix.This is the so-called discrete maximum principle (d.m.p) [4].

Theorem 2. Under the d.m.p, we have \[\left\Vert \omega_{h}-\tilde{\omega}_{h}\right\Vert _{\infty}\leq\left\Vert \psi-\tilde{\psi}\right\Vert _{\infty}.\]

Theorem 3. [5] Let \(\psi\in W^{1,\infty}\left( \Omega\right)\) and \(\omega\in W^{1,\infty}\left( \Omega\right)\), then \[\left\Vert \omega-\omega_{h}\right\Vert _{\infty}\leq Ch.\]

3. The continuous QVI

Next, we give some qualitative results for the solution of the QVI (1). Indeed, starting from \(\hat{u}^{0}\in\) \(H^{1}(\Omega),\)the unique solution of the equation  \[a(\hat{u}^{0},v)=(f,v)\text{ }\forall v\in H^{1}(\Omega). \label{EQ} \tag{10}\]

We define the continuous sequence \(\ \left( \hat{u}_{k}^{n}\right)\) such that \(\hat{u}_{k}^{n}\) \(\in H^{1}(\Omega)\) solves the variational inequality (VI) \[\left\{ \begin{array} [c]{l} a(\hat{u}_{k}^{n},v-\hat{u}_{k}^{n})\geq(f,v-\hat{u}_{k}^{n})\text{ }\forall v\in H^{1}(\Omega),\\ \hat{u}_{k}^{n}\leq k+M\hat{u}_{k}^{n-1},v\leq k+M\hat{u}_{k}^{n-1}. \end{array} \right. \label{ITER1} \tag{11}\]

Similarly, starting from \(\check{u}^{0}=0\), we define the continuous sequence \(\left( \check{u}_{k}^{n}\right)\) such that \(\check{u}_{k}^{n}\in H^{1}(\Omega)\) solves the variational inequality (VI) \[\left\{ \begin{array} [c]{l} a(\check{u}_{k}^{n},v-\check{u}_{k}^{n})\geq(f,v-\check{u}_{k}^{n})\text{ }\forall v\in H^{1}(\Omega),\\ \check{u}_{k}^{n}\leq k+M\check{u}_{k}^{n-1},v\leq k+M\check{u}_{k}^{n-1}. \end{array} \right. \label{ITER2} \tag{12}\]

Note that due to standard comparison results in elliptic variational inequalities, the sequences \(\left( \hat{u}_{k}^{n}\right)\) and \(\left( \check{u}_{k}^{n}\right)\) are decreasing and increasing, respectively.

Theorem 4. [6] Assume that \(f\geq f_{0}>0\). Then, the sequences\(\left( \hat {u}_{k}^{n}\right)\) and \(\left( \check{u}_{k}^{n}\right)\) defined in (11) and (12), respectively, converge decreasingly and increasingly to the unique solution \(u_{k}\) of QVI. (1).

Definition 1. [6] A function \(w_{k}\) \(\in H^{1}(\Omega)\) is said to be a subsolution for QVI (1) if \[\left\{ \begin{array} [c]{l} a(w_{k},v)\,\leq\,\,(f,v)\,\,\,\forall v\in H^{1}(\Omega),\text{ } v\geq0,\\ w_{k}\leq k+Mw_{k},\text{ \ }v\leq k+Mw_{k}. \end{array} \right. \label{CSUB2} \tag{13}\]

Let \(\mathbb{X}\) denote the set of such subsolutions.

Theorem 5. [6] The solution of QVI (1) is the maximum element of \(\ \mathbb{X}.\)

Notation 1. Let \(k,\) \(\tilde{k}\) be two positive parameters and \(u_{k}=\partial(k)\), \(u_{\tilde{k}}=\) \(\partial(\tilde{k})\) be the corresponding solutions to QVI (1), respectively.

Lemma 1 (Monotonicity property). \(k\geq\) \(\tilde{k}\) implies \(u_{k}\geq u_{\tilde{k}}\).

Proof. We proceed by induction. For \(n=0\), since \(M\) is nondecreasing and \(k\geq\) \(\tilde{k}\), we have \(k+M\hat{u}_{k}^{0}\geq\tilde{k}+M\hat{u}_{k}^{0},\) which implies that \(\hat{u}_{k}^{1}\geq\) \(\hat{u}_{\tilde{k}}^{1}.\)Let us now assume that \(\hat{u}_{k}^{n-1}\geq\) \(\hat{u}_{\tilde{k}}^{n-1}\). Then, as \(k\geq\) \(\tilde{k}\), we have \(k+M\hat{u}_{k}^{n-1}\geq\tilde{k}+M\) \(\hat{u}_{\tilde {k}}^{n-1}\), and making use of standard comparison results in elliptic VIs, we get \(\hat{u}_{k}^{n}\geq\hat{u}_{\tilde{k}}^{n}.\) Hence, due to Theorem 1, passing to the limit, as \(n\rightarrow\infty,\) it follows that \(u_{k}\geq u_{\tilde{k}.}\) ◻

Theorem 6 (Lipschitz property). Under conditions of Lemma 1,we have \[\left\Vert u_{k}-u_{\tilde{k}}\right\Vert _{\infty}\leq\left\vert k-\tilde {k}\right\vert .\label{CLIPK} \tag{14}\]

Proof. Set \[\Phi=\left\vert k-\tilde{k}\right\vert .\]

Then, since  \[\begin{aligned} k & \leq\tilde{k}+\left\vert k-\tilde{k}\right\vert \\ & \leq\tilde{k}+\Phi, \end{aligned}\] applying Lemma 1, we get \[\partial(\,k)\leq\partial(\,\tilde{k}+\Phi). \]

But, \(\partial(\,\tilde{k}+\Phi)\) being the solution of QVI with source term \(f\) and impulse cost \(\tilde{k}+\Phi\) , and \(\partial(\tilde{k})+\Phi\) being the solution of QVI with source term \(f+c.\Phi\) and impulse cost \(\tilde {k}+\Phi\), we have, \(f+c\Phi\geq f\), and hence, by standard comparison result for elliptic QVIs, \(\partial(\,\tilde{k}+\Phi)\leq\partial(\tilde{k})+\Phi\). Thus, \[\partial(\,k)\leq\partial(\,\tilde{k})+\Phi.\]

In a similar way we can also get \[\partial(\,\tilde{k})\leq\partial(k)+\Phi,\] which completes the proof. ◻

Corollary 1. we have \[\left\Vert u_{k}-u_{0}\right\Vert _{\infty}\leq k.\label{Rate} \tag{15}\]

Proof. It suffices to take the limit on (14), as \(\tilde{k}\rightarrow0.\) ◻

4. The discrete QVI

4.1. Qualitative properties

Next, we shall give analog qualitative properties for the discrete QVI (4).Their respective proofs will be omitted as they are similar to those of their continuous analogues. Indeed, starting from \(u_{h}^{0}\), the unique solution to the equation

\[a(\hat{u}_{h}^{0},v)=(f,v)\text{ }\forall v\in\mathbb{V}_{h},\] we define the continuous sequence \(\left( \hat{u}_{k,h} ^{n}\right)\) such that \(\hat{u}_{k,h}^{n}\) \(\in\mathbb{V}_{h}\) solves the discrete variational inequality (VI) \[\left\{ \begin{array} [c]{l} a(\hat{u}_{k,h}^{n},v-u_{k,h}^{n})\geq(f,v-\hat{u}_{k,h}^{n})\text{ }\forall v\in\mathbb{V}_{h},\\ \hat{u}_{k,h}^{n}\leq k+r_{h}M\hat{u}_{k,h}^{n-1},v\leq k+r_{h}M\hat{u} _{k,h}^{n-1}. \end{array} \right. \label{DUNUP} \tag{16}\]

Likewise, starting from \(\check{u}_{h}^{0}=0\), we define the continuous sequence \((\check{u}_{k,h}^{n})\) such that \(\check {u}_{k,h}^{n}\in\mathbb{V}_{h}\) solves the discrete variational inequality (VI) \[\left\{ \begin{array} [c]{l} a(\check{u}_{k,h}^{n},v-\check{u}_{h}^{n})\geq(f,v-\check{u}_{k,h}^{n})\text{ }\forall v\in\mathbb{V}_{h},\\ \check{u}_{k,h}^{n}\leq k+r_{h}M\check{u}_{k,n-1h}^{n},v\leq k+r_{h}M\check {u}_{k,n-1h}^{n}. \end{array} \right. \label{DUND} \tag{17}\]

Note that, due to standard comparison results in discrete elliptic variational inequalities, the sequences \(\left( \hat{u}_{h}^{n}\right)\) and \((\check {u}_{h}^{n})\) are decreasing and increasing, respectively.

Theorem 7.[2] Let the d.m.p hold. Then, the sequences \(\left( \hat{u}_{k,h} ^{n}\right)\) and \((\check{u}_{k,h}^{n})\) defined in (16)and (17), respectively, converge decreasingly and increasingly to the unique solution \(u_{h}\) of QVI (4).

Definition 2. A function \(w_{h}\in\) \(\mathbb{V}_{h}\) is said to be a subsolution for the QVI (4) if \[\left\{ \begin{array} [c]{l} a(w_{k,h},v)\,\leq\,\,(f,v)\,\,\,\forall\varphi_{s},\\ w_{k,h}\leq k+r_{h}Mw_{h},\text{ }v\leq k+r_{h}Mw_{h}. \end{array} \right. \tag{18}\]

Let \(\mathbb{X}_{h}\) denote the set of such subsolutions.

Theorem 8. [2] The solution of QVI (4) is the maximum element of \(\ \mathbb{X}_{h}.\)

Notation 2. Let \(k,\tilde{k}\) be two positive parameters and \(u_{k,h}=\partial_{h}(k),\) \(u_{\tilde{k},h}=\partial_{h}(\,\tilde{k})\) be the corresponding solutions to (4), respectively.

Lemma 2(Monotonicity property).Let the d.m.p hold.Then\(\,\,k\geq\) \(\tilde {k}\) implies that \(u_{k,h}\geq u_{\tilde{k},h}.\)

Theorem 9(Lipschitz property).Uder conditions of Lemma 2, we have \[\left\Vert u_{k,h}-u_{\tilde{k},h}\right\Vert _{\infty}\leq\left\vert k-\tilde{k}\right\vert .\label{DLIPK} \tag{19}\]

Corollary 2. We have \[\left\Vert u_{k},_{h}-u_{0,h}\right\Vert_{\infty}\leq k. \tag{20}\]

5. Error analysis

This section is devoted to demonstrate the proposed approximation method converges in the \(L^{\infty}\left( \Omega\right)\) norm. The proof stands on the construction of an appropriate subsolution for each of the continuous and discrete QVIs, respectively. From now on, \(C\) will denote a constant independent of both \(h\) and \(k\)).

5.1. Error estimate for the QVI (1)

Theorem 10. We have \[\left\Vert u_{k}-u_{k,h}{}\right\Vert _{\infty}\leq Ch.\]

The proof stands on the following two fundamental lemmas.

Lemma 3(Construction of a discrete subsolution). Assume that \[\max\left\{ \left\Vert Mu_{k}\right\Vert _{W^{1,\infty}\left( \Omega\right) }\text{, }\left\Vert u_{k}\right\Vert _{W^{1,\infty}\left( \Omega\right) }\right\} \leq C.\label{REGUL 1} \tag{21}\]

Then, there exists a function \(\alpha_{k,h}\) such that: \[\left\{ \begin{array} [c]{l} \alpha_{k,h}\leq u_{k,h},\qquad \text{and}\ \\ \left\Vert \alpha_{k,h}-u_{k}\right\Vert _{\infty}\leq Ch.\ \end{array} \right.\]

Proof. Let us first, just for the sake of simplicity, adopt the notations: \[u=u_{k},\text{ }u_{h}=\ u_{k,h}.\]

Now, let \(\bar{u}_{h}\) be the solution of the discrete variational inequality (VI) \[\left\{ \begin{array} [c]{l} a(\bar{u}_{h},v-\bar{u}_{h})\geqq(f,v-\bar{u}_{h})\,\,\,\forall v\in \mathbb{V}_{h},\\ \bar{u}_{h}\leq k+r_{h}Mu,\,v\leq k+r_{h}Mu. \end{array} \right. \label{AUXVI1} \tag{22}\] where \(u\) is the solution of (1).Then, under \(W^{1,\infty}\)-regularity (21)of \(u\) and \(Mu\), using standard \(L^{\infty}\) error estimate for VIs [5], we get \[\left\Vert u-\bar{u}_{h}\right\Vert _{\infty}\leq Ch.\label{EREST1} \tag{23}\]

Now, since \(\bar{u}_{h}\) is solution to the discrete VI (22), it is also a discrete subsolution, that is, \[\left\{ \begin{array} [c]{l} a(\bar{u}_{h},\varphi_{s})\leq(f,\varphi_{s})\,\,\,\forall s=1,2,…,m(h),\\ \bar{u}_{h}\leq k+r_{h}Mu. \end{array} \right.\]

On the other hand, since \(M\) and \(r_{h}\) are both Lipschitz with constant equal to \(1\), using (23), yields \[\begin{aligned} \bar{u}_{h} & \leq k+r_{h}Mu-r_{h}M\bar{u}_{h}+r_{h}M\bar{u}_{h} \,\\ & \leq\left\Vert k+r_{h}Mu-r_{h}M\bar{u}_{h}\right\Vert _{\infty}+r_{h} M\bar{u}_{h}\,\\ & \leq k+r_{h}M\bar{u}_{h}\,+\left\Vert u-\bar{u}_{h}\right\Vert _{\infty}\\ & \leq k+Ch\,+r_{h}M\bar{u}_{h}\,. \end{aligned}\]

So, \(\bar{u}_{h}\) is a subsolution for the QVI with impulse cost \(\tilde{k}=k+Ch\)

Let \(\bar{U}_{h}\) be the solution of such a QVI, that is, \(\bar{U}_{h}\) \(=\partial_{h}(k+Ch).\)Then, as \(u_{h}=\partial_{h}(k),\) making use of Theorem 9, we get \[\begin{aligned} \left\Vert u_{h}-\bar{U}_{h}\right\Vert _{\infty} & \leq C\left\vert k-(k+Ch)\right\vert \\ & \leq Ch. \end{aligned}\]

Hence, due to Theorem 8, we have \[\bar{u}_{h}\leq\bar{U}_{h}\leq u_{h}+Ch,\] and putting \[\alpha_{h}=\bar{u}_{h}-Ch,\] we get \[\alpha_{h}\leq u_{h}.\]

Finally, using (23), we obtain \[\begin{aligned} \left\Vert \alpha_{h}-u\right\Vert _{\infty} & \leq\left\Vert \bar{u} _{h}-u\right\Vert _{\infty}+Ch\\ & \leq Ch. \end{aligned}\] ◻

Lemma 4(Construction of a continuous subsolution).Assume that \[\max\left\{ \left\Vert Mu_{k,h}\right\Vert _{W^{1,\infty}\left( \Omega\right) }\text{, }\left\Vert u_{k,h}\right\Vert _{W^{1,\infty}\left( \Omega\right) }\right\} \leq C. \tag{24}\]

Then, there exists a function \(\beta_{k}^{(h)}\) such that: \[\left\{ \begin{array} [c]{l} \beta_{k}^{(h)}\leq u_{k},\qquad \text{and}\\ \left\Vert \beta_{k}^{(h)}-u_{k,h}\right\Vert _{\infty}\leq Ch. \end{array} \right.\]

Proof. Let \(\bar{u}\), be the solution of the continuous variational inequality \[\left\{ \begin{array} [c]{l} a(\bar{u},v-\bar{u})\geqq(f,v-\bar{u})\,\forall v\in H^{1}\left( \Omega\right),\\ \bar{u}\leq k+Mu_{h},\,v\leq k+Mu_{h}, \end{array} \right. \label{AUXVI2} \tag{25}\] with obstacle \(\psi=k+Mu_{h}\). It is clear that \(u_{h}=\sigma_{h}(k+Mu_{h})\) is nothing but the approximation of \(\bar{u}=\sigma(k+Mu_{h})\). Hence, as \(\left\Vert \psi\right\Vert _{W^{1,\infty}\left( \Omega\right) }\leq C\), making use of standard \(L^{\infty}\) – error estimate for VIs [5], we get, \[\left\Vert \bar{u}-u_{h}\right\Vert _{\infty}\leq Ch. \label{ERBETA} \tag{26}\]

Now, as the solution of VI (25) is a subsolution, we have \[\left\{ \begin{array} [c]{l} a(\bar{u},v)\leq(f,v)\,\,\,\forall v\in H^{1}(\Omega)\text{, }v>0,\\ \bar{u}\leq k+Mu_{h}\,,\,\,\,v\leq k+Mu_{h}, \end{array} \right.\] but, as \[\begin{aligned} \bar{u} & \leq k+Mu_{h}-M\bar{u}+M\bar{u}\\ & \leq k+\left\Vert Mu_{h}-M\bar{u}\right\Vert _{\infty}+M\bar{u} \\ & \leq k+\left\Vert u_{h}-\bar{u}\right\Vert _{\infty}+M\bar{u}\\ & \leq k+Ch+M\bar{u},\, \end{aligned}\] it follows that \(\bar{u}\) is a subsolution for the QVI with impulse cost \(\tilde{k}=k+Ch\) .

Denote by \(\bar{U}\) the solution of such a QVI, that is, \(\bar{U} =\partial\left( k+Ch\right)\). Then, as \(u=\partial(k)\), making use of Theorem 6, we have \[\begin{aligned} \left\Vert u-\omega\right\Vert _{\infty} & \leq C\left\vert k-(k+Ch)\right\vert \\ & \leq Ch \end{aligned}\]

Hence, making use of Theorem 5, we obtain \[\bar{u}\leq\bar{U}\leq u+Ch.\] Now, putting \[\beta^{(h)}=\bar{u}-Ch\text{,}\] we clearly have \[\beta^{(h)}\leq u.\]

Finally, using (26), we obtain \[\begin{aligned} \left\Vert \beta^{(h)}-u_{h}\right\Vert _{\infty} & \leq\left\Vert \bar {u}-Ch-u_{h}\right\Vert _{\infty}\\ & \leq Ch, \end{aligned}\] which completes the proof. ◻

Now, combining Lemmas 3 and 4, we are in position to derive the main result, we are in a position to prove Theorem 10.

Proof. Indeed, making use of both Lemmas 3 and 4, we have \[\begin{aligned} u_{k,h} & \leq\beta_{k}^{(h)}+Ch\\ & \leq u_{k}+Ch\\ & \leq\alpha_{k,{\large h}}+Ch \end{aligned}\]

Thus \[\left\Vert u_{k}-u_{k,h}{}\right\Vert _{\infty}\leq Ch,\] which completes the proof. ◻

Corollary 3. We have \[\text{ }\left\Vert u_{0}-u_{0,h}{}\right\Vert _{\infty}\leq Ch.\]

Proof. Using Theorem 10, and Corollaries 1 and 2, we get \[\begin{aligned} \left\Vert u_{0}-u_{0,h}{}\right\Vert _{\infty}&\leq \left\Vert u_{0} -u_{k}\right\Vert _{\infty}+\left\Vert u_{k}-u_{k,h}\right\Vert {}_{\infty }+\left\Vert u_{k,h}-u_{0,h}\right\Vert {}_{\infty}\\ & \leq Ck+Ch+Ck. \end{aligned}\]

Thus, passing to the limit, as \(k\rightarrow0\), the desired result follows. ◻

6. Conclusion

We have derived maximum error estimates of the standard finite element approximation of elliptic impulse control quasi variational inequalities (QVIs) with vanishing control cost.For this purpose, we developped a method that stands, in both the continuous settings, on intrinsic qualitative properties of the QVI like the continuous dependence with respect to the control cost and the characterization of the solution of the QVI as the least upper bound of the set of subsolutions.

References

  1. Perthame, B. (1988).Vanishing impulse cost in the quasi-variational inequality for ergodic impulse control. Asymptotic Analysis, 1(1), 13-21.

  2. Cortey-Dumont, P. (1980). Approximation numérique d’une inéquation quasi variationnelle liée à des problèmes de gestion de stock. RAIRO. Analyse Numérique, 14(4), 335-346.

  3. Bensoussan, A., & Lions, J. L. (2011). Applications of Variational Inequalities in Stochastic Control (Vol. 12). Elsevier.

  4. Lu, C., Huang, W., & Qiu, J. (2014). Maximum principle in linear finite element approximations of anisotropic diffusion–convection–reaction problems. Numerische Mathematik, 127(3), 515-537.

  5. Cortey-Dumont, P. (1985). On the finite element approximation in the \(L^{\infty}\) norm of variational inequalities with nonlinear operators. Numerish Mathematik, 47(1), 45-57.

  6. Bensoussan, A., & Lions, J. L. (1984). Impulse Control and Quasi-Variational Inequalities. Gauthiers-Villars, Paris.