- Research
- Open access
- Published:
Singularly perturbed second order semilinear boundary value problems with interface conditions
Boundary Value Problems volume 2015, Article number: 47 (2015)
Abstract
In this paper we study a class of singularly perturbed interface boundary value problems with discontinuous source terms. We first establish a lemma of lower-upper solutions by using the Schauder fixed point theorem. By the method of boundary functions and the lemma of lower-upper solutions we obtain the existence, asymptotic estimates, and uniqueness of the solution with boundary and interior layers for the proposed problem.
1 Introduction
We consider the interface boundary value problem with discontinuous source terms
where ε is a small and positive parameter, A and B being given constants, and
The functions \(f_{1}\) and \(f_{2}\) are smooth enough on \([a,d]\times\mathbb {R}\) and on \([d,b]\times\mathbb{R}\), respectively, and \(f_{1}(d,u) \neq f_{2}(d,u)\) for \(u\in\mathbb{R}\). We denote by \([w](d)\) the jump of a function w at the point d, i.e. \([w](d)=w(d^{+})-w(d^{-})\). The solution of the problem and its first order derivative have a jump at \(x=d\in(a,b)\). The original problem can be regarded as the coupling of the left problem and the right problem.
In recent years, the boundary value problems with interface conditions have appeared in applications such as heat transfer in layers composite material [1], the one dimensional metal-oxide-semiconductor structure [2], and population genetics [3]. The study of boundary value problems with interface conditions has attracted much attention, especially in numerical aspects; see [2, 4–8] and references therein. For instance, de Falco and O’Riordan [2] considered singularly perturbed reaction-diffusion equations with discontinuous data and interface conditions, where one of the problems was given as follows: Find \(u_{\varepsilon}\in C^{0}[0,1]\cap C^{2}((0,d)\cup(d,1))\) such that
with the diffusion coefficient
They suggested a parameter-uniform method based on piecewise-uniform Shishkin meshes to solve the problem above.
In this paper, we are devoted to the study of the existence, uniqueness, and asymptotics of the singularly perturbed interface boundary value problem (1.1)-(1.3), whose solution exhibits an interior layer due to the discontinuity of the source term. By using the Schauder fixed point theorem, we first establish a lower and upper solutions lemma which is an extension of classical theory of lower and upper solutions (see [9], for instance). By the method of boundary functions (see [10, 11], for example) and the lemma of lower-upper solutions we obtain the existence, asymptotic estimates, and uniqueness of the solution with boundary and interior layers for the proposed problem.
The remainder of this paper is organized as follows. In Section 2 we establish the lemma of lower-upper solutions for a class of two-point boundary value problems with interface conditions by the Schauder fixed point theorem, which will be used to prove our main result. With the asymptotic expansions and the lemma of lower and upper solutions established in Section 2, the asymptotic estimates, existence, and uniqueness of the solution for the problem (1.1)-(1.3) are obtained in Section 3.
2 Basic lemmas
For \(j=1, 2\), we define the vector spaces \(Q^{j}[a,b]\) of functions by
We define \(\widetilde{Q}^{2}[a,b]\) by
where \(\Lambda=\{d,x_{1},x_{2},\ldots,x_{k}\}\) and k is a given positive integer.
It is easily verified that \(Q^{1}[a,b]\) is a Banach space endowed with the norm \(\|u\|=\|u\|_{\infty}+\|u'\|_{\infty}\). The following lemma is a generalized Arzelà -Ascoli theorem on families of functions in \(Q^{1}[a,b]\).
Lemma 2.1
Assume that a bounded set E in \(Q^{1}[a,b]\) is piecewise equicontinuous, that is, \(E|_{[a,d)}=\{u|_{[a,d)}: u\in E\}\) is equicontinuous on \([a,d)\) and \(E|_{(d,b]}=\{u|_{(d,b]}: u\in E\}\) is equicontinuous on \((d,b]\). Then E is a relatively compact subset of \(Q^{1}[a,b]\).
Proof
The proof follows almost the same lines as that of the classical Arzelà -Ascoli theorem (see, for instance, [12]), noting that for \(u\in Q^{1}[a,b]\), u, and \(u'\) both have left and right limits at \(x=d\). Thus, details are omitted here. □
Let us consider the two-point boundary value problem with interface conditions
where the constants A, B, \(c_{1}\), \(c_{2}\), and the function \(f(x,u)\) are given in (1.1)-(1.3). Obviously, all solutions to (2.1)-(2.3) belong to \(Q^{2}[a,b]\).
Definition 2.1
We say that a function \(\alpha\in\widetilde{Q}^{2}[a,b]\) is a lower solution of the problem (2.1)-(2.3) if
We say that a function \(\beta\in\widetilde{Q}^{2}[a,b]\) is an upper solution of the problem (2.1)-(2.3) if
Lemma 2.2
Assume that α and β are lower and upper solutions of the problem (2.1)-(2.3) such that \(\alpha\leq\beta\). Then the problem (2.1)-(2.3) has at least one solution \(u\in Q^{1}[a,b]\) such that for all \(x\in[a,d)\cup(d,b]\),
Proof
Let us consider the modified problem
where \(\gamma(x,u)\) is defined by
The homogeneous problem of (2.6) and (2.8) can be divided into the following two initial value problems:
and
Obviously, the problems (2.10)-(2.11) and (2.12)-(2.13) have general solutions
and
respectively.
Consider the modified problem (2.6)-(2.8). Let us write the boundary value problem as an integral equation
where
and
If \(x\in[a,d)\), the integral equation becomes
If \(x\in(d,b]\), the integral equation becomes
In order that the solutions of the two problems satisfy the interface conditions (2.7) at \(x=d\), we must have
Substituting the integral equations into (2.17) and (2.18), and considering (2.11) and (2.13), we must have
Consequently, we obtain
Define an operator \(T: Q^{1}[a,b]\rightarrow Q^{1}[a,b]\) as follows:
Observe that \(f(x,\gamma(x,u))-\gamma(x,u): ([a,d)\cup (d,b] )\times\mathbb{R}\rightarrow\mathbb{R}\) is uniformly bounded in \(u\in Q^{1}[a,b]\). It follows that from Lemma 2.1 that the set \(T (Q^{1}[a,b] )\) is a relatively compact subset of \(Q^{1}[a,b]\). Moreover, T is continuous. Hence, it follows from the Schauder fixed point theorem that T has at least one fixed point \(u(x)\in Q^{1}[a,b]\).
We are now ready to prove that each solution \(u(x)\) of the problem (2.6)-(2.8) satisfies \(\alpha(x)\leq u(x)\leq\beta (x)\) for \(x\in[a,d)\cup(d,b]\).
Let us first prove that \(u(x)\leq\beta(x)\). Suppose, on the contrary, that the function \(h(x)=u(x)-\beta(x)\) has a positive maximum at some \(x_{0}\in[a,d)\cup(d,b]\). First, if \(x_{0}\neq d^{\pm},x_{i}^{\pm}\), we have \(h(x_{0})>0\), \(h'(x_{0})=0\), \(h''(x_{0})\leq0\), then
which yields a contradiction. Second, if \(x_{0}=d^{-}\) or \(x_{0}=d^{+}\), we have
so we can see that the function \(h(x)=u(x)-\beta(x)\) is continuous. Then we have
But
which yields also a contradiction. Third, if \(x_{0}=x_{i}^{-}\) or \(x_{0}=x_{i}^{+}\), we can similarly reach a contradiction. Therefore, we prove that \(u(x)\leq\beta(x)\) for all \(x\in[a,d)\cup(d,b]\).
In a similar way, we can prove that \(\alpha(x)\leq u(x)\) for all \(x\in [a,d)\cup(d,b]\).
Therefore, the solution of (2.6)-(2.8) is also that of (2.1)-(2.3) and satisfies \(\alpha(x)\leq u(x)\leq\beta (x)\) for \(x\in[a,d)\cup(d,b]\). The proof of the lemma is completed. □
Lemma 2.3
Assume that function \(f(x,u)\) is given in (1.1)-(1.3), and \(f(x,u)\) is strictly increasing with respect to u. Then the problem (2.1)-(2.3) has at most one solution.
Proof
Assume u and w are distinct solutions of the given problem (2.1)-(2.3). Without loss of generality, we can assume that there is some \(x\in[a,d)\cup(d,b]\) such that \(u(x)>w(x)\). Let
Since
and
we have \(h(x)\in C^{1}[a,b]\). Owing to \(h(a)=h(b)=0\), \(h(x)\) achieves the positive maximum at some \(x_{0}\in(a,b)\).
First, if \(x_{0}\neq d\), we have
However, on the other hand,
This is a contradiction.
Second, if \(x_{0}=d\), noting that \(h(a)=0\) and \(h(d)>0\), there exists some \(x^{*}\in[a,d)\) such that \(h(x)>0\) for \(x\in[x^{*},d]\) and \(h'(x^{*})\geq0\). However, on the other hand,
which yields a contradiction. Then the proof is complete. □
3 Main results
In this section, we are interested in the asymptotic behavior of solution with respect to the small parameter ε, as well as the existence and uniqueness for the problem (1.1)-(1.3). For the sake of simplicity, we only consider the approximation of zero order.
Because of discontinuity at \(x=d\), the original problem (1.1)-(1.3) can be viewed as the coupling of the left problem
and the right problem
where \(\gamma_{L}(\varepsilon)\), \(\gamma_{R}(\varepsilon)\) are constants dependent on ε which will be determined later.
We first make the following basic assumptions.
-
(H1)
The functions \(f_{1}\in C^{2} ([a,d]\times \mathbb{R} )\), \(f_{2}\in C^{2} ([d,b]\times\mathbb{R} )\), and
$$\begin{aligned}& f_{1}(d,u)\neq f_{2}(d,u), \quad \mbox{for } u\in\mathbb{R}, \\& \frac{\partial f_{1}}{\partial u}(x,u)\geq\sigma_{0}>0,\quad \mbox {for } (x,u) \in[a,d]\times\mathbb{R}, \\& \frac{\partial f_{2}}{\partial u}(x,u)\geq\sigma_{0}>0,\quad \mbox {for } (x,u) \in[d,b]\times\mathbb{R}. \end{aligned}$$ -
(H2)
The reduced problem \(f(x,u)=0\) has a solution
$$ u(x)=\left \{ \begin{array}{l@{\quad}l} \varphi(x), &x\in[a,d), \\ \psi(x), &x\in(d,b], \end{array} \right . $$such that \(\varphi(x)\in C^{2}[a,d]\) and \(\psi(x)\in C^{2}[d,b]\).
Let us construct the formal asymptotic solution of the left problem (3.1)-(3.2). We seek the formal solution of the form
where
Putting (3.5) into (3.1) and equating the coefficients of like powers of ε we can obtain a series of recursive equations determining these coefficients. From (H2) it follows that \(u_{0}^{(-)}(x)=\varphi(x)\), \(x\in[a,d)\). The left boundary layer term \(Q_{0}^{(-)}(\tau)\) solves
subject to the boundary conditions
The right boundary layer term \(V_{0}^{(-)}(\eta)\) solves
subject to the boundary conditions
Therefore, the formal approximation of zero order for the left problem (3.1)-(3.2) reads
In a completely similar way, we can get the formal approximation of zero order for the right problem (3.3)-(3.4), where \(Q_{0}^{(+)}(\eta)\) solves
and \(V_{0}^{(+)}(\xi)\) solves
The following two lemmas are concerned with the asymptotic behavior of the boundary layer terms for the left problem, whose proofs are essentially similar to that of Lemma 3.2 in [13] and thus are omitted here.
Lemma 3.1
Under the assumptions (H1) and (H2), for sufficiently small \(\varepsilon >0\), the problem (3.6) and (3.7) has a solution \(Q_{0}^{(-)}(\tau)\) satisfying the following estimate:
Lemma 3.2
Under the assumptions (H1) and (H2), for sufficiently small \(\varepsilon >0\) the problem (3.8) and (3.9) has a solution \(V_{0}^{(-)}(\eta)\) satisfying the following estimate:
In a similar way, we can prove that \(Q_{0}^{(+)}(\eta)\) and \(V_{0}^{(+)}(\xi)\) satisfy the following estimates:
and
In order that the solutions of the two problems satisfy the interface conditions (1.2) at \(x=d\), we must have
By substituting the formal solutions into (3.14) we have
In view of (3.9), multiplying (3.8) by \(2\frac {dV_{0}^{(-)}}{d\eta}\) and integrating from −∞ to 0 we obtain
where
Analogously, it follows from (3.10)-(3.11) that
where
Noting that \(V_{0}^{(-)}(\eta) \frac{dV_{0}^{(-)}}{d\eta}>0\) for \(\eta \leq0\) and \(Q_{0}^{(+)}(\eta)\frac{dQ_{0}^{(+)}}{d\eta}<0\) for \(\eta \geq0\) we have from (3.17) and (3.18)
Putting (3.19)-(3.20) into (3.14) yields
Substituting the formal solutions into (3.15) we have
Noting that \(\gamma_{0}^{(+)}\), \(\gamma_{0}^{(-)}\) can be determined by (3.21) and (3.22). We suppose that:
-
(H3)
The system of equations (3.21) and (3.22) has a solution \((\gamma_{0}^{(+)}, \gamma_{0}^{(-)})\).
Likewise, \(\gamma_{i}^{(\pm)}\) (\(i\geq1\)) can be determined recursively.
Thus, we obtain the formal asymptotic solution of the problem (1.1)-(1.3). Now we are in a position to state our main result.
Theorem 3.1
Let the conditions (H1)-(H3) hold. Then for sufficiently small \(\varepsilon >0\) the boundary value problem (1.1)-(1.3) has a unique solution \(u(x,\varepsilon )\in Q^{1}[a,b]\) satisfying:
where
Proof
In order to prove our main result we need to construct suitable upper and lower solutions. To this end, we need to make some necessary preparations.
Firstly, from the construction of asymptotic solution and the assumption (H1)-(H3) we know that there is a positive constant M such that
and, moreover,
Now we select the barrier functions as follows:
where
and
with \(x_{1}= \frac{a+d}{2}\) and \(x_{2}= \frac {b+d}{2}\), and positive constants \(\lambda_{i}\) (\(i=1,2,3,4\)) to be determined later.
It is easily verified that \(Q_{1}^{(-)}(\tau)\), \(V_{1}^{(-)}(\eta )\), \(Q_{1}^{(+)}(\eta)\), and \(V_{1}^{(+)}(\xi)\) have the following properties:
-
(i)
\(Q_{1}^{(-)}(\tau)\), \(V_{1}^{(-)}(\eta)\), \(Q_{1}^{(+)}(\eta)\) and \(V_{1}^{(+)}(\xi )\) are nonnegative functions;
-
(ii)
$$ {\biggl.\frac{dV_{1}^{(-)}}{d\eta}\biggr|_{\eta =0}}=-\frac{M}{4\sigma_{0}}< 0,\qquad {\biggl.\frac{dQ_{1}^{(+)}}{d\eta}\biggr|_{\eta=0}}=\frac {M}{4\sigma_{0}}>0; $$
-
(iii)
\(V_{1}^{(-)}(\eta)\) and \(V_{1}^{(+)}(\xi)\) are solutions of the equation
$$- \frac{d^{2}y(t)}{dt^{2}}+\sigma_{0}y(t)+Mt{\mathrm{e}}^{\sqrt {\sigma_{0}}t}=0, $$and \(Q_{1}^{(-)}(\tau)\), \(Q_{1}^{(+)}(\eta)\) are solutions of the equation
$$- \frac{d^{2}y(t)}{dt^{2}}+\sigma_{0}y(t)-Mt{\mathrm{e}}^{-\sqrt {\sigma_{0}}t}=0. $$
In order to make sure the functions α and β are continuous at \(x_{1}\), \(x_{2}\) and meet
the constants \(\lambda_{i}\) (\(i=1,2,3,4\)) that we select must satisfy
With (3.23) and the definitions of α and β, we easily verify that
Using the above properties (i)-(iii) and the assumptions (H1)-(H3), if \(x\in(a,x_{1})\) we have
provided ε is small enough and \(\lambda_{1}\) is large enough, where \(0\geq\theta_{1}, \theta_{2}\geq1\), and
If \(x\in(x_{1},d)\) we have
provided ε is small enough and \(\lambda_{2}\) is large enough, where \(0\geq\theta_{3},\theta_{4}\geq1\), and
Similarly, we have for ε small enough
We can prove in a similar way that for sufficiently small \(\varepsilon>0\),
It follows from Lemma 2.2 and Lemma 2.3 that the boundary value problem (1.1)-(1.3) has a unique solution \(u(x,\varepsilon)\in Q^{1}[a,b]\) such that
The proof is completed. □
Finally, we present an example to illustrate our results. Consider an interface two-point boundary value problem:
where
The left reduced problem \(f_{1}(x,u)=u-x=0\) has a unique smooth solution
and the right reduced problem \(f_{2}(x,u)=u+x=0\) has a unique smooth solution
The problems (3.6)-(3.7), (3.8)-(3.9), (3.10)-(3.11), and (3.12)-(3.13) for the boundary layer terms have the solutions
Putting (3.24)-(3.27) into (3.21) and (3.16), we can see that
By Theorem 3.1 the above problem has the unique solution
where
References
Mikhailov, M, Özisik, M: Unified Analysis and Solutions of Heat and Mass Diffusion. Dover, New York (1994)
de Falco, C, O’Riordan, E: Interior layers in a reaction-diffusion equation with a discontinuous diffusion coefficient. Int. J. Numer. Anal. Model. 7, 444-461 (2010)
Chen, CK: A fixed interface boundary value problem for differential equations: a problem arising from population genetics. Dyn. Partial Differ. Equ. 3, 199-208 (2006)
Aitbayev, R: Existence and uniqueness for a two-point interface boundary value problem. Electron. J. Differ. Equ. 2013, 242 (2013)
Farrell, PA, Hegarty, AF, Miller, JJH, O’Riordan, E, Shishkin, GI: Global maximum norm parameter-uniform numerical method for a singularly perturbed convection-diffusion problem with discontinuous convection coefficient. Math. Comput. Model. 40, 1375-1392 (2004)
Farrell, PA, O’Riordan, E, Shishkin, GI: A class of singularly perturbed quasilinear differential equations with interior layers. Math. Comput. 78, 103-127 (2009)
Huang, Z: Tailored finite point method for the interface problem. Netw. Heterog. Media 4, 91-106 (2009)
Loubenets, A, Ali, T, Hanke, M: Highly accurate finite element method for one-dimensional elliptic interface problems. Appl. Numer. Math. 59, 119-134 (2009)
Coster, CD, Habets, P: Two-Point Boundary Value Problems: Lower and Upper Solutions. Elsevier, New York (2006)
Vasiléva, AB, Butuzov, VF, Kalachev, LV: The Boundary Function Method for Singular Perturbation Problems. SIAM Studies in Applied Mathematics. SIAM, Philadelphia (1995)
de Jager, EM, Jiang, F: The Theory of Singular Perturbations. North-Holland, Amsterdam (1996)
Zeidler, E: Applied Functional Analysis: Applications to Mathematical Physics. Springer, New York (1995)
Xie, F: On a class of singular boundary value problems with singular perturbation. J. Differ. Equ. 252, 2370-2387 (2012)
Acknowledgements
The authors are grateful to the referees for their valuable comments. The second author was supported by the Natural Science Foundation of Shanghai (No. 15ZR1400800), in part by the National Natural Science Foundation of China (No. 11371087) and by the Fundamental Research Funds for the Central Universities.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.
Rights and permissions
Open Access This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.
About this article
Cite this article
Lin, H., Xie, F. Singularly perturbed second order semilinear boundary value problems with interface conditions. Bound Value Probl 2015, 47 (2015). https://doi.org/10.1186/s13661-015-0309-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-015-0309-5