Abstract
In this paper a standard numerical method with piecewise linear interpolation on Shishkin mesh is suggested to solve singularly perturbed boundary value problem for second order ordinary delay differential equations with discontinuous convection coefficient and source term. An error estimate is derived by using the supremum norm and it is of almost first order convergence. Numerical results are provided to illustrate the theoretical results.
Keywords
Singularly perturbed problem; Convection–diffusion problem; Discontinuous convection coefficient; Shishkin mesh; Delay
2010 Mathematics Subject Classification
65L10; 65L11; 65L12
1. Introduction
Singularly perturbed ordinary differential equations with a delay are ordinary differential equations in which the highest derivative is multiplied by a small parameter and involving at least one delay term. Such type of equations arises frequently from the mathematical modelling of various practical phenomena, for example, in the modelling of the human pupil-light reflex [14], the study of bistable devices [4] and variational problems in control theory [10], etc. It is important to develop suitable numerical methods to solve singularly perturbed differential equations with a delay, whose accuracy does not depend on the parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon } , that is the methods are uniformly convergent with respect to the parameter.
In the past, only very few people had worked in the area Numerical Methods to Singularly Perturbed Delay Differential Equation(SPDDE). But in the recent years, there has been growing interest in this area. The authors of [12], [6], [16], [1] and [2] suggested some numerical methods for singularly perturbed delay differential equations with continuous data. Recently few authors in [20], [21] and [17] suggested some numerical method for singularly perturbed delay differential equations with discontinuous data.
In the present paper, as mentioned in the above abstract, motivated by the works of [7], [3] and [13], we consider the following singularly perturbed boundary value problem (2.1) for second order ordinary delay differential equations with discontinuous convection coefficient and suggest a parameter uniform numerical method. It is proved that this method is uniformly convergent of order Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle O(N^{-1}ln^2N)} .
The present paper is organized as follows. In Section 2, the problem of study with discontinuous data is stated. Existence of the solution to the problem is established in Section 3. A maximum principle of the DDE is established in Section 4. Further a stability result is derived. Analytical results of the problem are derived in Section 5. The present numerical method is described in Section 6 and an error estimate is derived in Section 7. Section 8 presents numerical results.
2. Statement of the problem
Through out the paper, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle C,C_1}
denote generic positive constants independent of the singular perturbation parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon }
and the discretization parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle N}
of the discrete problem. Further, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle I_N}
denotes Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \lbrace 0,1,\ldots ,N\rbrace }
. The supremum norm is used for studying the convergence of the numerical solution to the exact solution to a singular perturbation problem: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Vert u\Vert _{\Omega }=sup_{x\in \Omega }\vert u(x)\vert } .
Motivated by the works of [8], [3] and [13], we consider the following BVP for SPDDE.
Find Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u\in Y=C^0(\overline{\Omega })\cap C^1(\Omega )\cap C^2(\Omega ^{{_\ast}})}
such that
|
(2.1) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle 0<\epsilon \ll 1} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a,f}
are sufficiently smooth and bounded in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega ^{{_\ast}}}
. The function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle b}
is a sufficiently smooth function on Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \overline{\Omega }}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega =(0,2)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \overline{\Omega }=[0,2]} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega ^{{_\ast}}=\Omega ^{-}\cup \Omega ^+} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega ^{-}=(0,1)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega ^+=(1,2)}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi }
is smooth on Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle [-1,0]}
.
The above problem (2.1) is equivalent to
|
(2.2) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u(1-)}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u(1+)}
denote the left and right limits of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u}
at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=1}
, respectively.
3. Existence result
For the reader’s convenience some known results are briefly reported on this section and in Section 4. They can be used here with some modifications.
Theorem 3.1.
The problem (2.1) has a solution Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u\in C^0(\bar{\Omega })\cap C^1(\Omega )\cap C^2(\Omega ^{{_\ast}})} .
Proof.
The proof is by construction. Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle y_1}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle y_2}
be particular solutions of the DDEs,
|
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle y_1=\phi (x),x\in [-1,0]} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a_1,a_2\in C^2(\bar{\Omega })}
with the above properties.
Consider the function
|
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi _1} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi _2}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi _3}
are the solutions of the following problems, respectively:
|
and
|
It is easy to see that the above function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle y}
satisfies the differential equation (2.1) and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u(0)=y(0)} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u(2)=y(2)}
. Using the similar arguments given in [7, Theorem 1], and [19, Theorems 2,3] one can prove the existence of the solution. □
Note: For the existence of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi _i,i=1,2,3}
one may refer to [18] and [5].
4. Stability result
Theorem 4.1 Maximum Principle.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle w\in C^0(\overline{\Omega })\cap C^2(\Omega ^{{_\ast}})} be any function satisfying Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle w(0)\geq 0} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle w(2)\geq 0} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Pw(x)\geq 0,\forall x\in \Omega ^{{_\ast}}} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle w^{{'}}(1+)-w^{{'}}(1-)=[w^{{'}}](1)\leq 0} . Then Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle w(x)\geq 0,\forall x\in \overline{\Omega }} .
In the following we use the function
|
(4.1) |
Proof.
Using the above function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle s}
and the procedure adopted in [20, Theorem 3.1], one can prove this theorem. □
Corollary 4.2 Stability Result.
For any Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u\in Y} we have
|
(4.2) |
Proof.
Using the barrier function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \psi ^{\pm }(x)=CC_1s(x)\pm u(x),x\in \overline{\Omega }} , where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle C_1=max\lbrace \vert u(0)\vert ,\vert u(2)\vert ,sup_{\xi \in \Omega ^{{_\ast}}}\vert Pu(\xi )\vert \rbrace }
and the procedure adopted in [20, Theorem 3.2], we can prove this corollary. □
Note: An immediate consequence of the Corollary 4.2 is that, the solution of the BVP (2.1) is unique.
5. Analytical results
Theorem 5.1.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u} be the solution of the problem (2.1), then we have the following bounds
|
Proof.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x\in \Omega ^{-}} . Then we have,
|
Integrating (2.2) from 0 to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x}
we get,
|
Therefore,
|
By the mean value theorem there exists a Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle z\in (0,\epsilon )}
such that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert \epsilon u^{{'}}(z)\vert \leq 2\Vert u\Vert _{\overline{\Omega }}}
. Therefore Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon \vert u^{{'}}(0)\vert \leq C(\Vert u\Vert _{\overline{\Omega }}+\Vert f\Vert _{\Omega }+\Vert \phi \Vert _{[-1,0]})} . Hence,
|
Similarly one can show that, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon \vert u^{{'}}(x)\vert \leq C,x\in \Omega ^+} .
From (2.2) it is easy to show that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Vert u^{(k)}\Vert _{\Omega ^{{_\ast}}}\leq C\epsilon ^{-k},k=2,3} . Hence the proof. □
To derive uniform error estimates, we need sharper bounds on the derivatives of the solution Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u} . We derive these using the following decomposition of the solution into smooth and singular components Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u(x)=v(x)+w(x)}
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v}
can be written in the form Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v=v_0+\epsilon v_1+\epsilon ^2v_2}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v_0,v_1}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v_2}
are defined respectively to be the solutions of the following problems:
Find Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v_0\in C^0(\Omega ^{{_\ast}})\cap C^1(\Omega ^{{_\ast}})}
such that
|
(5.1) |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v_1\in C^0(\Omega ^{{_\ast}})\cap C^1(\Omega ^{{_\ast}}\cup \lbrace 2\rbrace )}
such that
|
(5.3) |
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v_2\in Y^{{_\ast}}}
such that
|
(5.5) |
We assume that, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Vert v_0^{{''}}\Vert _{\Omega ^{{_\ast}}}\leq C}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Vert v_1^{{'''}}\Vert _{\Omega ^{{_\ast}}}\leq C}
.
Thus the smooth component Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v}
satisfies the following:
find Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v\in C^0(\Omega ^{{_\ast}}\cup \lbrace 0,2\rbrace )\cap C^2(\Omega ^{{_\ast}})}
such that
|
(5.7) |
Further Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle w}
satisfies the problem, that is, find Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle w\in C^0(\Omega ^{{_\ast}}\cup \lbrace 0,2\rbrace )\cap C^2(\Omega ^{{_\ast}})}
such that
|
(5.8) |
Note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v+w=u\in Y^{{_\ast}}} .
Theorem 5.2.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle w} be the solutions of the regular and singular components of the solution Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u} . Then
|
Proof.
Integrating the differential equation (5.1)–(5.4) separately on Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega ^{-}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega ^+}
, we get Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Vert v_i\Vert \leq C,i=0,1}
and by the stability result we have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Vert v_2\Vert \leq C}
. Therefore Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Vert v\Vert _{\Omega ^{{_\ast}}}\leq C} . Similarly one can prove that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Vert v^{(k)}\Vert _{\Omega ^{{_\ast}}}\leq C(1+\epsilon ^{2-k}),k=0,1,2,3} .
Note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert w(x)\vert \leq \vert u(x)\vert +\vert v(x)\vert } . From the stability result we have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert u(1)\vert \leq C} . Further, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert v(1)\vert \leq C} . Therefore Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert w(1)\vert \leq \eta }
(say). Now consider the barrier function
|
It is easy to check that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi _1^{\pm }(0)=\eta exp(\frac{-\alpha }{\epsilon })\pm w(0)\geq 0,\varphi _1^{\pm }(1)=\eta \pm w(1)\geq 0} .
Applying the result given in [9, Theorem 2.1] on Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle [0,1]} , we get Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi _1^{\pm }(x)\geq 0} .
Consider the barrier function
|
It is easy to see that, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi _2^{\pm }(1)=C_1(\epsilon +1-\epsilon exp(\frac{-\alpha }{\epsilon }))\pm w(1)\geq 0} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi _2^{\pm }(2)=C_1(\epsilon +exp(\frac{-\alpha }{\epsilon })-\epsilon )\pm w(2)\geq 0} . Again applying the result given in [9, Theorem 2.1] on Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle [1,2]} , then we get Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi _2^{\pm }(x)\geq 0} . Using the procedure adopted in [7, Lemma 4], one can prove the rest of this theorem. □
Note: From the above theorem it is easy to show that
|
(5.9) |
6. Discrete problem
In this section, mesh selection strategy, namely piecewise uniform mesh (Shishkin mesh), is explained. Also upwind finite difference scheme with piecewise linear interpolation on Shishkin mesh for the problem (2.1) is described.
6.1. Mesh selection strategy
Since the BVP (2.1) exhibits strong interior layers at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=1}
and a weak boundary layer at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=2}
, we choose a piecewise uniform Shishkin mesh on Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle [0,2]} . For this we divide the interval Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle [0,2]}
into five subintervals, namely Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega _1=[0,1-\tau _1]}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega _2=[1-\tau _1,1]} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega _3=[1,1+\tau _2]} ,Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega _4=[1+\tau _2,2-\tau _2]} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Omega _5=[2-\tau _2,2]} , where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \tau _1=min\lbrace 0.5,\frac{2\epsilon lnN}{\alpha }\rbrace ,\tau _2=min\lbrace 0.25,\frac{2\epsilon lnN}{\alpha }\rbrace } . Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h_1=4N^{-1}(1-\tau _1)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h_2=4N^{-1}\tau _1} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h_3=8N^{-1}\tau _2} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h_4=4N^{-1}(1-2\tau _2)} . The mesh Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \overline{\Omega }^N=\lbrace x_0,x_1,\ldots ,x_N\rbrace }
is defined by
|
6.2. A finite difference scheme for (2.2)
On Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \overline{\Omega }^N} , we define the following scheme for the BVP (2.2):
|
(6.1) |
where
|
6.3. Discrete stability result
Lemma 6.1 Discrete Maximum Principle.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Z(x_i)} be a mesh function satisfying Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Z(x_0)\geq 0,Z(x_N)\geq 0} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle P^NZ(x_i)\geq 0,i\in I_N\setminus \lbrace 0,N/2,N\rbrace } and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle (D^+-D^{-})Z(x_{N/2})=[DZ](x_{N/2})\leq 0} . Then Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Z(x_i)\geq 0,\forall x_i\in \overline{\Omega }^N} .
Proof.
Define Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle s(x_i)=\lbrace \begin{array}{c} \frac{3}{2}+\frac{x_i}{2}\mbox{,}\quad x_i\in [0,1]\cap \overline{\Omega }^N\mbox{,}\\ 3-x_i\mbox{,}\quad x_i\in [1,2]\cap \overline{\Omega }^N\mbox{.} \end{array}}
Note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle P^Ns(x_i)>0,\forall x_i\in \Omega ^{{_\ast}}\cap \overline{\Omega }^N}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle [Ds](x_{N/2})<0}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle s(x_i)>0,\forall x_i\in \overline{\Omega }^N}
.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \mu ^{{_\ast}}=max\lbrace \frac{-Z(x_i)}{s(x_i)}:x_i\in \overline{\Omega }^N\rbrace } . Then there exists Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x_i^{{_\ast}}\in \overline{\Omega }^N}
such that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Z(x_i^{{_\ast}})+\mu ^{{_\ast}}s(x_i^{{_\ast}})=0}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Z(x_i)+\mu ^{{_\ast}}s(x_i)\geq 0,\forall x_i\in \overline{\Omega }^N}
. Therefore the mesh function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle (Z+\mu ^{{_\ast}}s)}
attains its minimum at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x_i=x_i^{{_\ast}}}
. Suppose the theorem does not hold true, then Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \mu ^{{_\ast}}>0} .
Case (i)Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle :(x_i^{{_\ast}}\in \Omega ^{-}\cap \overline{\Omega }^N)}
|
It is a contradiction.
Case (ii)Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle :(x_i^{{_\ast}}\in \Omega ^+\cap \overline{\Omega }^N)}
|
It is a contradiction.
Case (iii)Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle :(x_i^{{_\ast}}=x_{N/2})}
|
It is a contradiction. Hence the proof of the theorem. □
Lemma 6.2.
For any mesh function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U(x_i)} we have
|
Proof.
One can easily prove this lemma by using Lemma 6.1 and the discrete barrier function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi ^{\pm }(x_i)=CC_1s(x_i)\pm U(x_i),x_i\in \overline{\Omega }^N} , where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \begin{array}{l}\displaystyle C_1=max\lbrace \vert U(x_0)\vert ,\vert U(x_N)\vert \\\displaystyle max_{j\in I_N\setminus \lbrace 0,N/2,N\rbrace }P^NU(x_j)\rbrace \end{array}} . □
Analogous to the continuous function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u} , we decompose the numerical solution Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U(x_i)}
defined by (6.1)–(6.3) as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U(x_i)=V(x_i)+W(x_i)}
, where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle V(x_i)}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle W(x_i)}
satisfy the following:
|
(6.4) |
and
|
(6.5) |
Theorem 6.3.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U(x_i)} be the numerical solution of (2.2) defined by (6.1)–(6.3) and further let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle V(x_i)} be the numerical solution of (5.7) given by (6.4). Then,
|
Proof.
Consider a mesh function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \begin{array}{l}\displaystyle \varphi ^{\pm }(x_i)=C_1N^{-1}[s(x_i)+\eta (x_i)]+C_1\psi (x_i)\vert U(x_{N/2})-\\\displaystyle -V(x_{N/2})\vert \pm (U(x_i)-V(x_i)),i\in I_N\end{array}}
where
|
It is easy to see that, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi ^{\pm }(x_0)\geq 0}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi ^{\pm }(x_N)\geq 0}
for a suitable Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle C_1>0}
. Further,
|
Note that, for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle i\in I_N\setminus \lbrace 0,N/2,N\rbrace } , we have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle P^N(U(x_i)-V(x_i))=0} .
Hence Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle P^N\varphi ^{\pm }(x_i)\geq 0,i\in I_N\setminus \lbrace 0,N/2,N\rbrace }
by a proper choice of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle C_1}
.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x_i=x_{N/2}} , then Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \begin{array}{l}\displaystyle [D]\varphi ^{\pm }(x_i)=-C_1\frac{7N^{-1}}{2}-C_12\vert U(x_{N/2})-\\\displaystyle -V(x_{N/2})\vert \pm [[D]U(x_i)-[D]V(x_i)]\leq 0\end{array}} , by (6.2), (6.4) and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \begin{array}{l}\displaystyle (x_{i+1}-\\\displaystyle -x_{i-1})\vert \delta ^2v(x_i)\vert \leq max_{[x_{i-1},x_{i+1}]}\vert v^{''}(x)\vert N^{-1}\end{array}} [15, page 52]. Then by Lemma 6.1, we have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi ^{\pm }(x_i)\geq 0,\forall i\in I_N} . Hence the proof. □
7. Error analysis
In this section we derive an error estimate for the numerical solution obtained by the scheme (6.1)–(6.3) for the problem (2.1).
Lemma 7.1.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v} be the solution of the problem (5.7) and let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle V(x_i)} be its numerical solution defined by (6.4). Then, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert v(x_i)-V(x_i)\vert \leq CN^{-1},i\in I_N} .
Proof.
Now,
|
Since Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert v^I(x_i)-v(x_i-1)\vert \leq CN^{-2}}
[11], then Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert P^N(v(x_i)-V(x_i))\vert \leq CN^{-1},i\in I_N\setminus \lbrace 0,N/2,N\rbrace }
. Then by Lemma 6.2, we have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert v(x_i)-V(x_i)\vert \leq CN^{-1},i\in I_N} . Hence the proof. □
Lemma 7.2.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle w} be the solution to the problem (5.8) and let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle W(x_i)} be its numerical solution defined by (6.5). If Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon \leq CN^{-1}} , then we have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert w(x_i)-W(x_i)\vert \leq CN^{-1}ln^2N} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle i\in I_N} .
Proof.
Note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert w(x_i)-W(x_i)\vert \leq \vert u(x_i)-U(x_i)\vert +\vert v(x_i)-V(x_i)\vert } . Then by Eq. (5.9), Theorem 6.3 and Lemma 7.1, we have
|
Therefore
|
(7.1) |
Now consider a mesh function
|
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \tau =min\lbrace \tau _1,\tau _2\rbrace } . From (7.1), it is clear that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi ^{\pm }(x_{N/4})\geq 0}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi ^{\pm }(x_{5N/8})\geq 0}
for a suitable choice of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle C_1>0}
.
|
Note that,
|
Also note that, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert w^I(x_i)-w(x_i-1)\vert \leq CN^{-1}}
[11]. Further, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert P^N(w(x_i)-W(x_i))\vert \leq C_2\epsilon ^{-2}N^{-1}}
, where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle C_2>0}
a constant independent of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon }
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle N}
.
Therefore, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle P^N\varphi ^{\pm }(x_i)\geq 0,i\in \lbrace N/4+1,\ldots ,N/2-1,N/2+1,\ldots ,5N/8-1\rbrace } . Then by the Lemma 6.1, we have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert w(x_i)-W(x_i)\vert \leq CN^{-1}ln^2N,i=N/4+1,\ldots ,5N/8-1} . Hence the proof. □
Theorem 7.3.
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u} be the solution of the problem (2.2), Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U(x_i)} be its numerical solution defined by (6.1)–(6.3). Then Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \vert u(x_i)-U(x_i)\vert \leq CN^{-1}ln^2N,i\in I_N} .
Proof.
The desired estimate follows from the fact that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u=v+w,U=V+W}
and from the Lemma 7.1 and Lemma 7.2. □
8. Numerical examples
In this section, three examples are given to illustrate the numerical method discussed in this paper. We use the double mesh principle to estimate the error and compute the experiment rate of convergence in our computed solutions for all problems. For this we put Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle D_{\epsilon }^M=max_{0\leq i\leq M}\vert U_i^M-U_{2i}^{2M}\vert } , where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U_i^M}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle U_{2i}^{2M}}
are the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle i}
th components of the numerical solutions on meshes of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle M}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle 2M}
points respectively. We compute the uniform error and rate of convergence as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle D^M=max_{\epsilon }D_{\epsilon }^M\mbox{and}\quad p^M=log_2(\frac{D^M}{D^{2M}})}
. For the following examples the numerical results are presented for the values of perturbation parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon \in \lbrace 2^{-27},2^{-12},\cdots ,2^{-6}\rbrace } .
Example 8.1.
|
(8.1) |
Table 1 presents the values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle D^N}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle p^N}
for this problem. Fig. 1 and Fig. 2 represent the numerical solution and the maximum point wise error for this problem, respectively.
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon } | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle N}
(Number of mesh points) | ||||||
|---|---|---|---|---|---|---|---|
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \downarrow } | 16 | 32 | 64 | 128 | 256 | 512 | 1024 |
| 2−6 | 6.9037e−2 | 4.4673e−2 | 3.9152e−2 | 2.4284e−2 | 1.6478e−2 | 1.0510e−2 | 6.3683e−3 |
| 2−7 | 7.1024e−2 | 4.5980e−2 | 3.3034e−2 | 2.7566e−2 | 1.6621e−2 | 1.0508e−2 | 6.3483e−3 |
| 2−8 | 7.2001e−2 | 4.6628e−2 | 3.3467e−2 | 2.3775e−2 | 1.8534e−2 | 1.0523e−2 | 6.3774e−3 |
| 2−9 | 7.2486e−2 | 4.6951e−2 | 3.3683e−2 | 2.3943e−2 | 1.6407e−2 | 1.1663e−2 | 6.3716e−3 |
| 2−10 | 7.2727e−2 | 4.7112e−2 | 3.3790e−2 | 2.4026e−2 | 1.6488e−2 | 1.0447e−2 | 6.9885e−3 |
| 2−11 | 7.2847e−2 | 4.7192e−2 | 3.3844e−2 | 2.4068e−2 | 1.6528e−2 | 1.0484e−2 | 6.3334e−3 |
| 2−12 | 7.2907e−2 | 4.7232e−2 | 3.3871e−2 | 2.4089e−2 | 1.6548e−2 | 1.0502e−2 | 6.3508e−3 |
| 2−13 | 7.2937e−2 | 4.7253e−2 | 3.3884e−2 | 2.4099e−2 | 1.6558e−2 | 1.0511e−2 | 6.3595e−3 |
| 2−14 | 7.2952e−2 | 4.7263e−2 | 3.3891e−2 | 2.4104e−2 | 1.6563e−2 | 1.0516e−2 | 6.3638e−3 |
| 2−15 | 7.2960e−2 | 4.7268e−2 | 3.3894e−2 | 2.4107e−2 | 1.6566e−2 | 1.0518e−2 | 6.3660e−3 |
| 2−16 | 7.2964e−2 | 4.7270e−2 | 3.3896e−2 | 2.4108e−2 | 1.6567e−2 | 1.0519e−2 | 6.3671e−3 |
| 2−17 | 7.2966e−2 | 4.7271e−2 | 3.3897e−2 | 2.4109e−2 | 1.6567e−2 | 1.0520e−2 | 6.3676e−3 |
| 2−18 | 7.2967e−2 | 4.7272e−2 | 3.3897e−2 | 2.4109e−2 | 1.6568e−2 | 1.0520e−2 | 6.3679e−3 |
| 2−19 | 7.2967e−2 | 4.7272e−2 | 3.3897e−2 | 2.4109e−2 | 1.6568e−2 | 1.0520e−2 | 6.3680e−3 |
| 2−20 | 7.2967e−2 | 4.7272e−2 | 3.3898e−2 | 2.4109e−2 | 1.6568e−2 | 1.0520e−2 | 6.3681e−3 |
| 2−21 | 7.2967e−2 | 4.7273e−2 | 3.3898e−2 | 2.4109e−2 | 1.6568e−2 | 1.0520e−2 | 6.3681e−3 |
| 2−22 | 7.2967e−2 | 4.7273e−2 | 3.3898e−2 | 2.4109e−2 | 1.6568e−2 | 1.0520e−2 | 6.3682e−3 |
| 2−23 | 7.2967e−2 | 4.7273e−2 | 3.3898e−2 | 2.4109e−2 | 1.6568e−2 | 1.0520e−2 | 6.3682e−3 |
| 2−24 | 7.2967e−2 | 4.7273e−2 | 3.3898e−2 | 2.4109e−2 | 1.6568e−2 | 1.0520e−2 | 6.3682e−3 |
| 2−25 | 7.2967e−2 | 4.7273e−2 | 3.3898e−2 | 2.4109e−2 | 1.6568e−2 | 1.0520e−2 | 6.3682e−3 |
| 2−26 | 7.2967e−2 | 4.7273e−2 | 3.3898e−2 | 2.4109e−2 | 1.6568e−2 | 1.0520e−2 | 6.3681e−3 |
| 2−27 | 7.2967e−2 | 4.7273e−2 | 3.3898e−2 | 2.4110e−2 | 1.6568e−2 | 1.0521e−2 | 6.3682e−3 |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle D^N} | 7.2967e−2 | 4.7273e−2 | 3.9152e−2 | 2.7566e−2 | 1.8534e−2 | 1.1663e−2 | 6.9885e−3 |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle p^N} | 6.2625e−1 | 2.7192e−1 | 5.0620e−1 | 5.7270e−1 | 6.6825e−1 | 7.3887e−1 | – |
|
|
|
Fig. 1. Numerical solution of the problem stated in Example 8.1. |
|
|
|
Fig. 2. Maximum point wise error for the problem stated in Example 8.1. |
Example 8.2.
|
(8.2) |
Table 2 presents the values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle D^N}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle p^N}
for this problem. Fig. 3 and Fig. 4 represent the numerical solution and the maximum point wise error for this problem, respectively.
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon } | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle N}
(Number of mesh points) | ||||||
|---|---|---|---|---|---|---|---|
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \downarrow } | 16 | 32 | 64 | 128 | 256 | 512 | 1024 |
| 2−6 | 4.4334e−2 | 2.6901e−2 | 2.4542e−2 | 1.4287e−2 | 9.6205e−3 | 6.1830e−3 | 3.7665e−3 |
| 2−7 | 4.5728e−2 | 2.7859e−2 | 1.9444e−2 | 1.6827e−2 | 9.7456e−3 | 6.1893e−3 | 3.7555e−3 |
| 2−8 | 4.6412e−2 | 2.8332e−2 | 1.9796e−2 | 1.3927e−2 | 1.1337e−2 | 6.2057e−3 | 3.7820e−3 |
| 2−9 | 4.6751e−2 | 2.8566e−2 | 1.9970e−2 | 1.4069e−2 | 9.5858e−3 | 7.1177e−3 | 3.7784e−3 |
| 2−10 | 4.6920e−2 | 2.8683e−2 | 2.0057e−2 | 1.4140e−2 | 9.6544e−3 | 6.1463e−3 | 4.2869e−3 |
| 2−11 | 4.7004e−2 | 2.8742e−2 | 2.0101e−2 | 1.4176e−2 | 9.6887e−3 | 6.1776e−3 | 3.7477e−3 |
| 2−12 | 4.7046e−2 | 2.8771e−2 | 2.0122e−2 | 1.4193e−2 | 9.7058e−3 | 6.1932e−3 | 3.7624e−3 |
| 2−13 | 4.7067e−2 | 2.8785e−2 | 2.0133e−2 | 1.4202e−2 | 9.7144e−3 | 6.2010e−3 | 3.7697e−3 |
| 2−14 | 4.7077e−2 | 2.8792e−2 | 2.0139e−2 | 1.4207e−2 | 9.7186e−3 | 6.2049e−3 | 3.7734e−3 |
| 2−15 | 4.7083e−2 | 2.8796e−2 | 2.0141e−2 | 1.4209e−2 | 9.7208e−3 | 6.2068e−3 | 3.7752e−3 |
| 2−16 | 4.7085e−2 | 2.8798e−2 | 2.0143e−2 | 1.4210e−2 | 9.7219e−3 | 6.2078e−3 | 3.7761e−3 |
| 2−17 | 4.7087e−2 | 2.8799e−2 | 2.0143e−2 | 1.4211e−2 | 9.7224e−3 | 6.2083e−3 | 3.7766e−3 |
| 2−18 | 4.7087e−2 | 2.8799e−2 | 2.0144e−2 | 1.4211e−2 | 9.7227e−3 | 6.2086e−3 | 3.7768e−3 |
| 2−19 | 4.7088e−2 | 2.8800e−2 | 2.0144e−2 | 1.4211e−2 | 9.7228e−3 | 6.2087e−3 | 3.7769e−3 |
| 2−20 | 4.7088e−2 | 2.8800e−2 | 2.0144e−2 | 1.4211e−2 | 9.7229e−3 | 6.2087e−3 | 3.7770e−3 |
| 2−21 | 4.7088e−2 | 2.8800e−2 | 2.0144e−2 | 1.4211e−2 | 9.7229e−3 | 6.2088e−3 | 3.7770e−3 |
| 2−22 | 4.7088e−2 | 2.8800e−2 | 2.0144e−2 | 1.4211e−2 | 9.7229e−3 | 6.2088e−3 | 3.7770e−3 |
| 2−23 | 4.7088e−2 | 2.8800e−2 | 2.0144e−2 | 1.4211e−2 | 9.7229e−3 | 6.2088e−3 | 3.7770e−3 |
| 2−24 | 4.7088e−2 | 2.8800e−2 | 2.0144e−2 | 1.4211e−2 | 9.7229e−3 | 6.2088e−3 | 3.7770e−3 |
| 2−25 | 4.7088e−2 | 2.8800e−2 | 2.0144e−2 | 1.4211e−2 | 9.7229e−3 | 6.2088e−3 | 3.7770e−3 |
| 2−26 | 4.7088e−2 | 2.8800e−2 | 2.0144e−2 | 1.4211e−2 | 9.7229e−3 | 6.2088e−3 | 3.7770e−3 |
| 2−27 | 4.7088e−2 | 2.8800e−2 | 2.0144e−2 | 1.4211e−2 | 9.7227e−3 | 6.2090e−3 | 3.7772e−3 |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle D^N} | 4.7088e−2 | 2.8800e−2 | 2.4542e−2 | 1.6827e−2 | 1.1337e−2 | 7.1177e−3 | 4.2869e−3 |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle p^N} | 7.0929e−1 | 2.3083e−1 | 5.4448e−1 | 5.6974e−1 | 6.7152e−1 | 7.3149e−1 | – |
|
|
|
Fig. 3. Numerical solution of the problem stated in Example 8.2. |
|
|
|
Fig. 4. Maximum point wise error for the problem stated in Example 8.2. |
Example 8.3.
|
(8.3) |
Table 3 presents the values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle D^N}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle p^N}
for this problem.
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon } | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle N}
(Number of mesh points) | ||||||
|---|---|---|---|---|---|---|---|
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \downarrow } | 16 | 32 | 64 | 128 | 256 | 512 | 1024 |
| 2−6 | 6.8150e−2 | 3.7767e−2 | 2.8472e−2 | 1.5380e−2 | 9.2391e−3 | 5.2438e−3 | 2.9335e−3 |
| 2−7 | 6.9961e−2 | 3.8976e−2 | 2.3702e−2 | 1.7956e−2 | 9.3354e−3 | 5.3363e−3 | 3.0515e−3 |
| 2−8 | 7.0863e−2 | 3.9595e−2 | 2.4127e−2 | 1.4829e−2 | 1.0949e−2 | 5.3842e−3 | 3.0680e−3 |
| 2−9 | 7.1312e−2 | 3.9913e−2 | 2.4339e−2 | 1.4992e−2 | 9.2106e−3 | 6.2932e−3 | 3.0871e−3 |
| 2−10 | 7.1537e−2 | 4.0072e−2 | 2.4446e−2 | 1.5074e−2 | 9.2781e−3 | 5.3106e−3 | 3.5991e−3 |
| 2−11 | 7.1649e−2 | 4.0151e−2 | 2.4499e−2 | 1.5115e−2 | 9.3119e−3 | 5.3405e−3 | 3.0583e−3 |
| 2−12 | 7.1705e−2 | 4.0191e−2 | 2.4525e−2 | 1.5135e−2 | 9.3288e−3 | 5.3554e−3 | 3.0725e−3 |
| 2−13 | 7.1733e−2 | 4.0211e−2 | 2.4539e−2 | 1.5145e−2 | 9.3372e−3 | 5.3629e−3 | 3.0795e−3 |
| 2−14 | 7.1747e−2 | 4.0221e−2 | 2.4545e−2 | 1.5150e−2 | 9.3415e−3 | 5.3666e−3 | 3.0831e−3 |
| 2−15 | 7.1754e−2 | 4.0226e−2 | 2.4549e−2 | 1.5153e−2 | 9.3436e−3 | 5.3685e−3 | 3.0848e−3 |
| 2−16 | 7.1758e−2 | 4.0228e−2 | 2.4550e−2 | 1.5154e−2 | 9.3446e−3 | 5.3694e−3 | 3.0857e−3 |
| 2−17 | 7.1760e−2 | 4.0229e−2 | 2.4551e−2 | 1.5155e−2 | 9.3452e−3 | 5.3699e−3 | 3.0862e−3 |
| 2−18 | 7.1760e−2 | 4.0230e−2 | 2.4552e−2 | 1.5155e−2 | 9.3454e−3 | 5.3701e−3 | 3.0864e−3 |
| 2−19 | 7.1761e−2 | 4.0230e−2 | 2.4552e−2 | 1.5155e−2 | 9.3456e−3 | 5.3702e−3 | 3.0865e−3 |
| 2−20 | 7.1761e−2 | 4.0230e−2 | 2.4552e−2 | 1.5155e−2 | 9.3456e−3 | 5.3703e−3 | 3.0865e−3 |
| 2−21 | 7.1761e−2 | 4.0231e−2 | 2.4552e−2 | 1.5155e−2 | 9.3457e−3 | 5.3703e−3 | 3.0866e−3 |
| 2−22 | 7.1761e−2 | 4.0231e−2 | 2.4552e−2 | 1.5155e−2 | 9.3457e−3 | 5.3703e−3 | 3.0866e−3 |
| 2−23 | 7.1761e−2 | 4.0231e−2 | 2.4552e−2 | 1.5155e−2 | 9.3457e−3 | 5.3703e−3 | 3.0866e−3 |
| 2−24 | 7.1761e−2 | 4.0231e−2 | 2.4552e−2 | 1.5155e−2 | 9.3457e−3 | 5.3704e−3 | 3.0866e−3 |
| 2−25 | 7.1761e−2 | 4.0231e−2 | 2.4552e−2 | 1.5155e−2 | 9.3457e−3 | 5.3704e−3 | 3.0866e−3 |
| 2−26 | 7.1761e−2 | 4.0231e−2 | 2.4552e−2 | 1.5155e−2 | 9.3457e−3 | 5.3703e−3 | 3.0866e−3 |
| 2−27 | 7.1761e−2 | 4.0231e−2 | 2.4552e−2 | 1.5155e−2 | 9.3455e−3 | 5.3703e−3 | 3.0868e−3 |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle D^N} | 7.1761e−2 | 4.0231e−2 | 2.8472e−2 | 1.7956e−2 | 1.0949e−2 | 6.2932e−3 | 3.5991e−3 |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle p^N} | 8.3491e−1 | 4.9876e−1 | 6.6508e−1 | 7.1360e−1 | 7.9898e−1 | 8.0616e−1 | – |
9. Discussion
A BVP for one type of SPDDEs is considered. To obtain an approximate solution to this type of problem, an upwind finite difference scheme with piecewise linear interpolation on Shishkin mesh is presented. The method is shown to be of almost first order convergence. This is very much reflected on the numerical results (Table 1, Table 2 and Table 3). Also Fig. 1 and Fig. 3 represent that the model problems stated in Example 8.1 and Example 8.2 exhibit strong interior layers at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=1}
and a weak boundary layer at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=2}
. Fig. 2 and Fig. 4 represent the maximum point wise error for the numerical solutions. Further these Fig. 2 and Fig. 4 represent the uniform convergence of the numerical method presented in this paper. The authors of [7] have considered second order ordinary differential equations with discontinuous convection coefficient with different signs on different subdomains. The solution to the problem considered in [7] exhibits strong interior layers at an interior point. Whereas the problem considered in this paper exhibits strong interior layers at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=1}
and weak boundary layer at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=2}
(see Theorem 5.2). This is due to the presence of the delay term with the differential equation. Therefore, to accommodate these interior layers and boundary layer in numerical solution, the Shishkin mesh Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \bar{\Omega }^N}
has been constructed in Section 6.1. In [7], the authors have suggested a uniform numerical method without interpolation, whereas the finite difference method with interpolation is needed in this paper, since the point Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x_i-1,i>N/2}
need not be a mesh point.
References
- [1] G.M. Amiraliyev, E. Cimen; Numerical method for a singularly perturbed convection–diffusion problem with delay; Appl. Math. Comput., 216 (2010), pp. 2351–2359
- [2] A.R. Ansari, S.A. Baker, G.I. Shishkin; A parameter-robust finite difference method for singularly perturbed delay parabolic partial differential equations; J. Comput. Appl. Math., 205 (2007), pp. 552–566
- [3] Z. Cen; A hybrid difference scheme for a singularly perturbed convection–diffusion problem with discontinuous convection coefficient; Appl. Math. Comput., 169 (2005), pp. 689–699
- [4] M.W. Derstine, H.M. Gibbs, F.A. Hopf, D.L. Kaplan; Bifurcation gap in a hybrid optical system; Phys. Rev. A, 26 (1982), pp. 3720–3722
- [5] P.W. Eloe, Youssef, N. Raffoul, C. Christopher, Tisdell; Existence, uniqueness and constructive results for delay differential equations; Electron. J. Differential Equations, 2005 (121) (2005), pp. 1–11
- [6] F. Erdogan; An exponentially fitted method for singularly perturbed delay differential equation; Adv. Differ. Equ., 2009 (2009), pp. 1–9 http://dx.doi.org/10.1155/2009/781579
- [7] P.A. Farrell, A.F. Hegarty, J.J.H. Miller, E. O’Riordan, G.I. Shishkin; Global maximum norm parameter-uniform numerical method for a singularly perturbed convection–diffusion problems with discontinuous convection coefficient; Math. Comput. Modelling, 40 (2004), pp. 1375–1392
- [8] P.A. Farrell, A.F. Hegarty, J.J.H. Miller, E. O’Riordan, G.I. Shishkin; Singularly perturbed convection diffusion problems with boundary and weak interior layers; J. Comput. Appl. Math., 166 (2004), pp. 133–151
- [9] P.A. Farrell, A.F. Hegarty, J.J.H. Miller, E. O’Riordan, G.I. Shishkin; Robust Computational Techniques for Boundary Layers; Chapman Hall/CRC, Boca Raton (2000)
- [10] V. Glizer; Asymptotic analysis and solution of a finite-horizon Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h_{\infty }} control problem for singularly-perturbed linear systems with small state delay; J. Optim. Theory Appl., 117 (2003), pp. 295–325
- [11] M.K. Jain, S.R.K. Iyengar, R.K. Jain; Numerical Methods for Scientific and Engineering Computation; New Age International Publisher (1993)
- [12] M. Kadalbajoo, K. Sharma; Numerical treatment of boundary value problems for second order singularly perturbed delay differential equations; Comput. Appl. Math., 24 (2) (2005), pp. 151–172
- [13] C.G. Lange, R.M. Miura; Singularly perturbation analysis of boundary-value problems for differential-difference equations; SIAM. J. Appl. Math., 42 (3) (1982), pp. 502–530
- [14] A. Longtin, J. Milton; Complex oscillations in the human pupil light reflex with mixed and delayed feedback; Math. Biosci., 90 (1988), pp. 183–199
- [15] J.J.H. Miller, E. O’Riordan, G.I. Shishkin; Fitted Numerical Methods for Singular Perturbation Problems: Error Estimates in the Maximum Norm for Linear Problems in One and Two Dimensions; World Scientific Publishing Co. Pte. Ltd, Singapore (2012)
- [16] J. Mohapatra, S. Natesan; Uniform convergence analysis of finite difference scheme for singularly perturbed delay differential equation on an adaptively generated grid; Numer. Math.:Theory Methods Appl., 3 (1) (2010), pp. 1–22
- [17] S. Nicaise, C. Xenophontos; Robust approximation of singularly perturbed delay differential equations by the hp finite element method; Comput. Methods Appl. Math., 13 (1) (2013), pp. 21–37
- [18] S.K. Ntouyas, G.Y. Sficas, P.C. Tsamato; An existence principle for boundary value problems for second order functional differential equations; Nonlinear Anal,Theory Methods Appl., 20 (3) (1993), pp. 215–222
- [19] M.H. Protter, H.F. Weinberger; Maximum Principle in Differential Equations; Prentice Hall, New York (1967)
- [20] V. Subburayan, N. Ramanujam; Asymptotic initial value technique for singularly perturbed convection diffusion delay problems with boundary and weak interior layers; Appl. Math. Lett., 25 (12) (2012), pp. 2272–2278
- [21] V. Subburayan, N. Ramanujam; An initial value technique for singularly perturbed convection–diffusion problems with a delay; J. Optim. Theory Appl., 158 (1) (2013), pp. 234–250 http://dx.doi.org/10.1007/s10957-012-0200-9
Document information
Published on 07/10/16
Licence: Other
Share this document
claim authorship
Are you one of the authors of this document?