- Research
- Open access
- Published:
A free boundary problem arising in the ecological models with N-species
Boundary Value Problems volume 2015, Article number: 4 (2015)
Abstract
This paper is concerned with the one-dimensional free boundary problem for quasilinear reaction-diffusion systems arising in the ecological models with N-species, where some of the species are made up of two separated groups and the mankind’s influence is taken into account. In the problem under consideration, there are n free boundaries, the coefficients of the equations are allowed to be discontinuous on the free boundaries and the reaction functions are mixed quasimonotone. The aim is to show the local existence of the solutions for the free boundary problem by the fixed point method, and the global existence and uniqueness of the solutions for the corresponding diffraction problem by the approximation and estimate methods.
1 Introduction
For some animal species, different groups may live in different habitats separated by free boundaries (see [1]). For the case of one-species and one-dimensional, Kim and Lin [1] proposed a model to describe the evolution of the free boundary. In this paper, we consider the one-dimensional ecological models with N-species where the mankind’s influence is taken into account. Assume that there are N species sharing the same habitat, the one-dimensional domain \((0, d)\). For some \(1\leq n \leq N\), assume further that each of the n species is divided into two groups, whereas each of the rest of \(N- n\) species has a single group.
To describe the free boundary problem, we introduce some notations:
-
\(Q_{T} :=(0, d)\times(0, T]\) is the open rectangle (and upper horizontal line).
-
\(\bar{Q}_{T}\) is the entire closed rectangle.
-
\(S_{T}\) is made up of the closed left and right boundaries of \(Q_{T}\).
-
\(\Gamma ^{l}_{T}:= \{(x,t): x=\varphi ^{l}(t),0\leq t \leq T \}\) is the free boundary for species l, \(l\in\{ 1,\ldots,n\}\).
-
\(Q^{l-}_{T}:= \{(x,t):0< x<\varphi ^{l}(t),0<t\leq T \}\), \(Q^{l+}_{T}:= \{(x,t):\varphi ^{l}(t)< x<d,0<t\leq T \}\) are the subdomains of \(Q_{T}\) separated by \(\Gamma _{T}^{l}\), \(l\in\{1,\ldots,n\}\).
-
\(Q_{k,T}:= \{(x,t):\varphi ^{k-1}(t)< x<\varphi ^{k}(t),0<t\leq T \}\) is the region between consecutive boundaries, \(k\in\{1,\ldots,n+1\}\), where \(\varphi ^{0}(t)\equiv0\), \(\varphi ^{n+1}(t)\equiv d\) for \(0\leq t\leq T\).
Let \(u^{l}=u^{l}(x,t)\) be the population density of the lth species. For any \(l\in\{1,\ldots,n\}\), the diffusion coefficients of the lth species in \(Q^{l-}_{T}\), \(Q^{l+}_{T}\) are positive constants \(\hat {a}^{l}\) and \(\tilde {a}^{l}\), respectively. For any \(l\in\{n+1,\ldots,N\}\), the diffusion coefficient of the lth species in the whole domain \(Q_{T}\) is positive constant \(a^{l}\). Then Fick’s law shows that the flux of the population density can be represented as
Assume that for each \(l\in\{1,\ldots,N\}\), harvesting rate or stocking rate for the lth species depends on its density and flux. Using the principle of conservation, we conclude that the vector functions \(\mathbf {u}=(u^{1},\ldots,u^{N})\) and \(\boldsymbol {\varphi }=(\varphi ^{1},\ldots, \varphi ^{n})\) are governed by quasilinear reaction-diffusion equations
where \(\hat {b}^{i}(u^{i})u^{i}_{x}\), \(\tilde {b}^{i}(u^{i})u^{i}_{x}\) (\(i=1,\ldots,n\)) and \(b^{j}(u^{j})u^{j}_{x}\) (\(j=n+1,\ldots,N\)) represent the mankind’s influence, and \(G^{l}(x,t,\mathbf {u})\) (\(l=1,\ldots,N\)) are the reaction functions defined by (see the example in Section 5)
Here and below,
On the parabolic boundary, u is required to satisfy the condition
For each \(l\in\{1,\ldots,n\}\), on the free boundary \(\Gamma ^{l}_{T}\), assume that the density and the flux are continuous (see [1]). By the conservation law of population, the rate of increase of population \(\varphi ^{l}_{t}(t)u^{l-}(\varphi ^{l}(t),t)\) is equal to the flux \(-\hat {a}^{l} u^{l-}_{x}(\varphi ^{l}(t),t) \), where \(u^{l-}(\varphi ^{l}(t),t)\) and \(u^{l+}(\varphi ^{l}(t),t)\) represent the limits of \(u^{l}\) from left and right for the space variable, respectively. Then u and φ are required to satisfy the free boundary conditions
Condition (1.4) is used to describe the rates of change of the free boundaries.
The motivation to deal with problem (1.1a)-(1.4) are some real processes. For example, consider 4-species system with two predators and two preys, where two predator species comprise two groups of wolves and two groups of lions, two prey species comprise a herd of sheep and a herd of horses, and four species share the same habitat. In each of predator species, the two groups live in different habitats separated by a free boundary and have different diffusion coefficients and biological habits. Then this ecological system can be described by problem (1.1a)-(1.4).
The aim of this paper is to investigate: (i) the global existence and uniqueness of the solutions for diffraction problem (1.1a)-(1.3); (ii) the local existence of the solutions for problem (1.1a)-(1.4).
The free boundary problems often appear in different fields, such as physics, ecology and chemistry. They have been investigated extensively in the literature (see [1–13] and the references therein). For the free boundary problems of single equations, the linear problems describing the flow of fluids in a one-dimensional porous medium were studied by Kamynin [2] and Evans and Friedman [3–5], and the quasilinear problem was studied by Liang [6]. The work in [1] is concerned with the other quasilinear problem describing ecological model with one-species. For the free boundary problems of parabolic systems, Mimura et al. [7], Ling et al. [8] and Lin [9] investigated the problems arising from two-species models in ecology, where the number of the free boundaries is one and the conditions on the free boundary are Stefan conditions. In this paper, there are n free boundaries in problem (1.1a)-(1.4) and the vector functions \(\mathbf{g}_{i}(\mathbf {u}):=(g^{1}_{i}(\mathbf {u}),\ldots,g^{N}_{i}(\mathbf {u}))\), \(i=1,\ldots,n\), are mixed quasimonotone.
Since for any given \(\boldsymbol {\varphi }(t)\), problem (1.1a)-(1.3) is a diffraction problem, we must first investigate the corresponding diffraction problem. In the study of diffraction problem, \(\Gamma ^{i}_{T}\), \(i=1,\ldots,n\), are called inner boundaries. The papers in [14–16] are concerned with the diffraction problems of parabolic systems with the inner boundaries in the special form \(\{x:\varphi _{k}(x)=0\}\times[0,T]\), \(k=1,\ldots,K\), where \(\varphi _{k}(x)\in C^{2+\alpha _{0}}\). In this paper, by the approximation and estimate methods, we show the existence and uniqueness of the solutions for diffraction problem (1.1a)-(1.3) with the inner boundaries in the form \(\Gamma ^{i}_{T}: x=\varphi ^{i}(t)\), \(i=1,\ldots,n\), where \(\varphi ^{i}(t)\in W^{1}_{2}(0,T)\).
The rest of this paper is organized as follows. In Section 2, we state the notations, definitions, hypotheses and main results. Section 3 is devoted to investigating the global existence and uniqueness of the solutions for the corresponding diffraction problem (1.1a)-(1.3). Section 4 is concerned with the local existence of solutions for the free boundary problem (1.1a)-(1.4) by Schauder’s fixed point theorem. Finally, in Section 5 we give an example.
2 The notations, hypotheses and main results
Fist, let us define some basic function spaces (see [17])
equipped with the norms
respectively, where \(D\subset\mathbb{R}\) is an arbitrary interval.
\(V_{2}(Q_{T})\) is the Banach space consisting of all elements of \(W_{2}^{1,0}(Q_{T})\) having a finite norm
\(V_{2}^{1,0}(Q_{T})\) is the Banach space consisting of all elements of \(V_{2}(Q_{T})\) that are continuous in t in the norm of \(L^{2}(0,d)\) with a finite norm
is the set of all functions in \(W_{2}^{1,1}(Q_{T})\) that vanish on \(S_{T}\) in the sense of trace. We write \(H^{1}(0,d):=W^{1}_{2}(0,d)\).
\(C^{\alpha }(\bar{Q}_{T})\) is the spaces of Hölder continuous functions in \(\bar{Q}_{T}\) with exponent \(\alpha \in(0, 1)\).
The product function spaces with N-components of \(W_{2}^{1,1}(Q_{T})\), \(V_{2}(Q_{T})\), \(V_{2}^{1,0}(Q_{T})\) and \(C^{\alpha }(\bar{Q}_{T})\) are denoted by \(\mathcal{W}_{2}^{1,1}(\mathcal{Q}_{T})\), \(\mathcal{V}_{2}(\mathcal{Q}_{T})\), \(\mathcal{V}_{2}^{1,0}(\mathcal{Q}_{T})\) and \(\mathcal{C}^{\alpha }(\bar{\mathcal{Q}}_{T})\), respectively, and the product function spaces with n-components of \(W^{1}_{p}(0,T)\), \(L^{p}(0,T)\) and \(C^{\alpha }([0,T])\) are denoted by \(\mathbf {W}^{1}_{p}(0,T)\), \(\mathbf {L}^{p}(0,T)\) and \(\mathbf {C}^{\alpha }([0,T])\), respectively. We note that the norm of \(\boldsymbol {\varphi }(t)\) in \(\mathbf {L}^{p}([0,T])\) is defined by
Let us recall the following definition. Write the vector u in the split form
and rewrite the functions \(f^{l}\) as
where \(\rho^{l}\), \(\omega^{l}\) are two nonnegative integers such that \(\rho^{l}+\omega^{l}=N-1\) and \([\mathbf {u}]_{\sigma}\) denotes a vector with σ number of component of u. The split form of u varies with respect to l and is determined by the quasimonotone property of \(f^{l}\) given as follows.
Definition 2.1
(see [16])
The vector function \(\mathbf{f}(\cdot, \mathbf {u}):=(f^{1}(\cdot, \mathbf {u}),\ldots,f^{N}(\cdot, \mathbf {u}))\) is said to be mixed quasimonotone in \(\mathcal{J}\subset\mathbb{R}^{N}\) with index vector \(([\mathbf {u}]_{\rho^{1}};\ldots;[\mathbf {u}]_{\rho^{N}})\) if for each \(l=1,\ldots,N\), there exist nonnegative integers \(\varrho^{l}\), \(\omega^{l}\), satisfying the relation
such that \(f^{l}(\cdot, u^{l},[\mathbf {u}]_{\rho^{l}},[\mathbf {u}]_{\omega^{l}})\) is nondecreasing in \([\mathbf {u}]_{\varrho^{l}}\) and is nonincreasing in \([\mathbf {u}]_{\omega^{l}}\) for all \(\mathbf {u}\in\mathcal{J}\).
Throughout the paper we make the following hypotheses:
(H)
-
(i)
\(d^{i}\), \(i=1,\ldots,n\), are constants satisfying the relation \(0< d^{1}<\cdots<d^{n}<d\). \(\psi^{l}(x,t)\), \(l=1,\ldots,N\), are all positive functions, and they possess the following properties: for some constants \(\alpha _{0}\in(0,1)\) and \(\chi_{0}\in(0,\min\{d^{1}, d-d^{n}\})\),
$$\begin{aligned}& \psi^{l}(x,t)\in C^{\alpha _{0}} \bigl([-\chi_{0}, d+\chi _{0}]\times[0,T] \bigr)\cap W^{1,1}_{2} \bigl((- \chi_{0}, d+\chi_{0})\times(0,T) \bigr), \\& \psi^{l}(x,t)\in C^{1+\alpha _{0}} \bigl([-\chi_{0}, \chi_{0}]\times[0,T] \bigr)\cap C^{1+\alpha _{0}} \bigl([d- \chi_{0}, d+\chi_{0}]\times[0,T] \bigr), \\& \quad l=1,\ldots,N, \\& \psi^{i}(x,0)\in H^{1}(0,d)\cap C^{2+\alpha _{0}} \bigl( \bigl[0,d^{i} \bigr] \bigr)\cap C^{2+\alpha _{0}} \bigl( \bigl[d^{i},d \bigr] \bigr),\quad i=1,\ldots,n, \\& \psi^{j}(x,0)\in C^{2+\alpha _{0}} \bigl([0,d] \bigr),\quad j=n+1, \ldots,N, \end{aligned}$$and
$$\begin{aligned} \hat {a}^{i}\psi^{i-}_{x} \bigl(d^{i},0 \bigr)=\tilde {a}^{i}\psi^{i+}_{x} \bigl(d^{i},0 \bigr),\quad i=1,\ldots,n. \end{aligned}$$ -
(ii)
There exist constant vectors \(\mathbf{M}:=(M^{1},\ldots ,M^{N})\) and \(\mathbf{m}:=(m^{1},\ldots,m^{N})\), \(\textbf{0}<\mathbf{m}\leq\mathbf{M}\), such that for each \(l=1,\ldots,N\),
$$ \left \{ \begin{array}{@{}l} g^{l}_{k}( M^{l},[\mathbf{M}]_{\varrho^{l}},[\mathbf{m}]_{\omega^{l}} )\leq 0,\\ g^{l}_{k}( m^{l},[\mathbf{m}]_{\varrho^{l}},[\mathbf{M}]_{\omega^{l}} )\geq0,\quad k=1,\ldots,n+1,\\ m^{l}< \psi^{l}(x,t)<M^{l} \quad((x,t)\in S_{T}\cup\{\Omega \times \{0\} \} ), \end{array} \right . $$(2.1)where the component \([\mathbf{\mathbf {u}}]_{\varrho^{l}}\) is independent of k. Let
The vector functions \(\mathbf{g}_{k}(\mathbf {u})=(g^{1}_{k}(\mathbf {u}),\ldots,g^{N}_{k}(\mathbf {u}))\), \(k=1,\ldots,n+1\), are mixed quasimonotone in with the same index vector \(([\mathbf {u}]_{\rho^{1}};\ldots;[\mathbf {u}]_{\rho^{N}})\).
-
(iii)
\(\hat {b}^{i}(u^{i}),\tilde {b}^{i}(u^{i}), b^{j}(u^{j})\in C^{1+\alpha _{0}}(\mathbb{R})\), \(i=1,\ldots,n\), \(j=n+1,\ldots,N\), and , \(k=1,\ldots,n+1\), \(l=1,\ldots,N\). There exists a positive nondecreasing function \(\mu(\theta)\) for \(\theta\in [0,+\infty)\) such that
$$ \left \{ \begin{array}{@{}l} |\hat {b}^{i}(u^{i}),\tilde {b}^{i}(u^{i})|\leq\mu(|u^{i}|),\quad i=1,\ldots,n,\\ |b^{j}(u^{j})|\leq\mu(|u^{j}|),\quad j=n+1,\ldots,N. \end{array} \right . $$(2.2)
Hypothesis (H)(i) implies that there exists a positive constant \(\mu_{1}\) such that
where \(\delta_{0}:=\min_{i=1,\ldots,n+1}(d^{i}-d^{i-1})\), \(d^{0}=0\) and \(d^{n+1}=d\). Hypothesis (H)(ii) is used to guarantee the existence of the coupled weak upper and lower solutions for the approximation problem of (1.1a)-(1.3). (H)(iii) shows that on \(Q_{T}\) the coefficients of the equations are allowed to be discontinuous, but on \(Q_{k,T}\) the reaction functions are required to be independent of \((x,t)\) and have Hölder continuous derivatives with respect to u, and on \(Q^{l-}_{T}\) and \(Q^{l+}_{T}\) the remaining coefficients of the lth equation are also required to be independent of \((x,t)\) and have Hölder continuous derivatives with respect to u.
For convenience, let
and
Definition 2.2
For a given \(\boldsymbol {\varphi }(t)\in \mathbf {W}^{1}_{2}(0,T)\), a vector function \(\mathbf {u}(x,t)\) is said to be a solution of diffraction problem (1.1a)-(1.3) corresponding to \(\boldsymbol {\varphi }(t)\) if (i) For each \(l=1,\ldots,N\), \(u^{l}\in C^{\alpha ,\alpha /2}(\bar{Q}_{T})\) for some \(0<\alpha <1\), \(u^{l}\in V_{2}^{1,0}(Q_{T})\), \(A^{l}(x,t) u^{l}_{x}\in V_{2}(Q_{T})\), \(u^{l}_{t}\in L^{2}(Q_{T})\), and \(u^{l}\in C^{2+\alpha ,1+\alpha /2}(\bar{Q}')\) for any given \(k\in\{1,\ldots,n+1\}\) and any given subdomain \(Q', \bar{Q}'\subset Q_{k,T}\); (ii) u satisfies pointwise the equations in (1.1a) and (1.1b) for \((x,t)\in Q_{k,T}\) (\(k=1,\ldots,n+1\)) and the parabolic boundary condition (1.2), and satisfies the inner boundary condition (1.3) for almost all \(t\in[0,T]\); (iii) For any vector function and any \(\tau\in[0,T]\), the following equalities hold:
Definition 2.3
\(\{\mathbf {u}(x,t), \boldsymbol {\varphi }(t)\}\) is said to be a solution of problem (1.1a)-(1.4) if (i) u is a solution of diffraction problem (1.1a)-(1.3) corresponding to \(\boldsymbol {\varphi }(t)\); (ii) \(\boldsymbol {\varphi }(t)\) is in \(\mathbf {C}^{2+\alpha }((0,T])\), and (1.4) holds for \(t\in(0,T]\).
Throughout the paper, \(\alpha _{1}\), \(\alpha _{2}\), \(C(\cdots)\) and \(C_{j}(\cdots)\) \((j=1,2, \ldots)\) denote constants depending only on T, d, \(d^{i}\), \(1/m^{i}\), \(M^{i}\) (\(i=1,\ldots,n\)), \(\mu(|\mathbf{M}|)\), \(\mu_{1}\), and the quantities appearing in parentheses. Constant C in different expressions may be different.
The main results in this paper are the following theorems.
Theorem 2.1
Let hypothesis (H) be satisfied.
(i) If φ is in \(\mathbf {W}^{1}_{2}(0, T)\) and satisfies
for some \(\delta>0\), then the corresponding diffraction problem (1.1a)-(1.3) has a unique solution u in . Moreover, the following estimates hold:
(ii) If \(\boldsymbol {\varphi }_{1}\), \(\boldsymbol {\varphi }_{2}\) are in \(\mathbf {W}^{1}_{2}(0, T)\) and satisfy (2.5), and if \(\mathbf {u}_{1}\), \(\mathbf {u}_{2}\) are the solutions in corresponding to \(\boldsymbol {\varphi }_{1}\) and \(\boldsymbol {\varphi }_{2}\), respectively, then
(iii) If φ is in \(\mathbf {C}^{2}((0, T])\) and u is the solution in , then for any given \(t_{1}\in(0,T)\) and any given \(l\in\{1,\ldots,n\}\), \(u^{l}_{x}(x,t)\) and \(u^{l}_{t}(x,t)\) are continuous with respect to x in \([\varphi ^{l}(t)-\delta/4,\varphi ^{l}(t)]\) and in \([\varphi ^{l}(t),\varphi ^{l}(t)+\delta/4]\) for almost all \(t\in[t_{1},T]\), and \(u^{l}_{xt}\) is square integrable in \(Q^{l}(\delta,t_{1})\), where \(Q^{l}(\delta,t_{1}):=\{(x,t):x\in(\varphi ^{l}(t)-\delta/4, \varphi ^{l}(t)+\delta/4),t\in(t_{1},T)\}\).
Theorem 2.2
Let hypothesis (H) hold. Then there exist \(T>0\) and a pair of vector functions \(\mathbf {u}(x,t)\), \(\boldsymbol {\varphi }(t)\), such that \(\{\mathbf {u}(x,t), \boldsymbol {\varphi }(t)\}\) is a solution to free boundary problem (1.1a)-(1.4).
3 The corresponding diffraction problem
Let \(\boldsymbol {\varphi }(t)\) be given and let (2.5) hold for some \(\delta>0\). The aim of this section is to prove Theorem 2.1.
3.1 An approximation problem
We will construct an approximation problem of (1.1a)-(1.3). We first construct some approximation functions.
For an arbitrary \(\varepsilon >0\), choose \(\boldsymbol {\varphi }_{\varepsilon }=(\varphi ^{1}_{\varepsilon }(t),\ldots, \varphi ^{n}_{\varepsilon }(t))\in \mathbf {C}^{2}([0, T])\) such that \(\boldsymbol {\varphi }_{\varepsilon }(t)\to \boldsymbol {\varphi }(t)\) in \(\mathbf {W}^{1}_{2}(0, T)\) as \(\varepsilon \to0\). Then, for small enough ε,
Let \(s_{\varepsilon }=s_{\varepsilon }(\theta)\) be a smooth function with values between 0 and 1 such that \(|\mathrm{d}s_{\varepsilon }(\theta)/\mathrm{d}\theta|\leq C/\varepsilon \) for all \(\theta\in\mathbb{R}\), \(s_{\varepsilon }(\theta)=1\) for \(\theta\leq0\) and \(s_{\varepsilon }(\theta)=0\) for \(\theta\geq \varepsilon \), and let
Define
and
Then according to hypothesis (H)(ii) and (iii), it follows that the vector function \(\mathbf{G}_{\varepsilon }(\cdot, \mathbf {u})=(G^{1}_{\varepsilon }(\cdot, \mathbf {u}),\ldots,G^{N}_{\varepsilon }(\cdot, \mathbf {u}))\) is mixed quasimonotone in with index vector \(([\mathbf {u}]_{\rho^{1}};\ldots;[\mathbf {u}]_{\rho^{N}})\), and
where \(\nu:=\min_{i=1,\ldots,n, j=n+1,\ldots,N}\{\hat {a}^{i},\tilde {a}^{i},a^{j}\}\). The definition of function \(s_{\varepsilon }(\theta)\) implies that
In addition, it is obvious from (2.5) that
where \(\bar{Q}_{T}\), \(\bar{Q}_{\tau,T}\) are the closure of \({Q}_{T}\) and \({Q}_{\tau,T}\), respectively. Thus by (3.9) and the definition of functions \(z_{\varepsilon ,k}(x,t)\), an argument similar to the one used in [16, Lemma 3.2] shows that
and
We next construct the approximation functions of \(\psi^{l}(x,t)\). Let \(\omega(|x|)\) be a sufficiently smooth nonnegative function such that \(\omega(|x|)=0\) for \(|x|\geq1\) and \(\int_{|x|\leq 1}\omega(x)\,\mathrm{d}x=1\), and let \(\lambda =\lambda (x)\) be a sufficiently smooth nonnegative function taking values in \([0,1]\) such that \(\lambda (x)=0\) for \(\delta\leq x\leq d-\delta\), \(\lambda (x)=1\) for \(x\leq\delta/2\) or \(x\geq d-\delta/2\), and \(\lambda _{x}(x)\leq C/\delta\) for all \(x\in \mathbb{R}\). Define
Then hypothesis (H)(i) and [17, Chapter II] imply that
and
and (2.1) and (2.3) imply that
In addition, for \(\varepsilon <\delta/8\),
Employing the above approximation functions, we consider the following approximation problem:
Lemma 3.1
Problem (3.16) has a unique classical solution \(\mathbf {u}_{\varepsilon }\) in , and the following estimates hold:
where constants \(\alpha _{2}\) and C are independent of ε.
Proof
In [15], by using the method of upper and lower solutions, together with the associated monotone iterations and various estimates, we investigated the existence and uniqueness of the global piecewise classical solutions of the quasilinear parabolic system with discontinuous coefficients and continuous delays under various conditions including mixed quasimonotone property of reaction functions. The same problem was also discussed for the system with continuous coefficients without time-delay.
It is obvious that problem (3.16) is the special case of [15, problem (1.1)] without discontinuous coefficients and time delays. Hypothesis (H)(ii) shows that \(\tilde{\mathbf {u}}=\mathbf{M}\), \(\hat{\mathbf {u}}=\mathbf{m}\) are the coupled weak upper and lower solutions of (3.16) in the sense of [15, Definition 2.2]. By (3.4)-(3.6) and (3.11)-(3.14), we conclude from [15, Theorem 4.1] that problem (3.16) has a unique classical solution \(\mathbf {u}_{\varepsilon }\) in . Furthermore, using (3.1), (3.5), (3.6) and (3.12)-(3.14), the proof similar to that of [16, Lemma 3.3] shows that estimates (3.17) and (3.18) hold.
To prove (3.19), we first fix \(l\in\{1,\ldots,n\}\). In view of (3.7), (3.8), (3.10) and (3.15), we find that \(u^{l}_{\varepsilon }\) is the solution of the following problems for single equation:
and
By (3.22), (3.23), (3.12) and (3.14), the proof similar to that of [17, Chapter VI, Lemma 3.1] gives (3.19) for \(l\in\{1,\ldots,n\}\). The similar argument shows that (3.19) holds for \(l\in\{n+1,\ldots,N\}\).
We next prove (3.20). For any fixed \(l\in\{1,\ldots,n\}\), let
and let
By a direct computation we have
Then (3.16), (3.2), (3.3) and (3.10) imply that the function \(q=q(y,t)\) satisfies
where \(E^{l}_{T}:=\{(y, t):-\varphi ^{l}_{\varepsilon }(t)< y< d-\varphi ^{l}_{\varepsilon }(t), 0<t\leq T\}\).
A double integration by parts gives
Thus multiplying the equation in (3.25) by \(q_{t}\), integrating it on \(E^{l}_{\tau}\) and using (3.5), (3.6), (3.24) and (3.26), we find that
Furthermore, by (3.12)-(3.14), (3.19) and Cauchy’s inequality, we deduce that for any \(\sigma>0\),
Choosing \(\sigma=\nu/4\), we get
Consequently,
This, together with Gronwall’s inequality (see [17, Chapter II, Lemma 5.5]), implies that
Hence we deduce from (3.25), (3.27), (3.5) and (3.6) that
which, together with (3.24), yields (3.20) for \(l\in \{1,\ldots,n\}\).
For any fixed \(l\in\{n+1,\ldots,N\}\), we consider the equality
A similar argument gives (3.20) for \(l\in\{n+1,\ldots,N\}\). Therefore, (3.20) holds for all \(l\in\{1,\ldots,N\}\).
It remains to prove (3.21). For each \(l\in\{1,\ldots,N\}\), since \(A^{l}_{\varepsilon }(x,t)u^{l}_{\varepsilon x}\) is in \(V_{2}(Q_{T})\), then (3.20) and [17, Chapter 2, formula (3.8)] show that
Thus (3.21) holds. □
Lemma 3.2
Let \(\mathbf {u}_{\varepsilon _{1}}(x,t)\), \(\mathbf {u}_{\varepsilon _{2}}(x,t)\) be the solutions in for problem (3.16) corresponding to \(\boldsymbol {\varphi }_{\varepsilon _{1}}(t)\) and \(\boldsymbol {\varphi }_{\varepsilon _{2}}(t)\), respectively. Then
where \(I_{3,0}:=\sum_{l=1}^{N}\|\psi^{l}_{\varepsilon _{1}}(x,0)-\psi^{l}_{\varepsilon _{2}}(x,0)\| _{H^{1}(0,d)}^{2} +\|\boldsymbol {\varphi }_{\varepsilon _{1}}-\boldsymbol {\varphi }_{\varepsilon _{2}}\|_{\mathbf {L}^{2}(0,T)}+\varepsilon _{1}+\varepsilon _{2}\).
Proof
Let \(\mathbf {w}=(w^{1},\ldots,w^{N})=\mathbf {u}_{\varepsilon _{1}}-\mathbf {u}_{\varepsilon _{2}}\), and let \(l\in\{1,\ldots,N\}\) be fixed. We see from (3.16) that
In view of (3.15), we find \(w^{l}(0,t)=w^{l}(d,t)=0\). Multiplying the equation in (3.29) by \(w^{l}\) and integrating by parts on \(Q_{\tau}\), we deduce that, for any \(\tau\in[0,T]\),
Let us estimate \(I^{l}_{3,1}\), \(I^{l}_{3,2}\) and \(I^{l}_{3,3}\). Since (3.2), (3.3) and (3.10) imply that
and
then it follows from Cauchy’s inequality that, for any \(\sigma>0\),
and
where
According to the definition of function \(s_{\varepsilon }(\theta)\), we see that \(s_{\varepsilon _{1}} (x-\varphi ^{l}_{\varepsilon _{1}}(t))-s_{\varepsilon _{2}}(x-\varphi ^{l}_{\varepsilon _{2}}(t))=0\) if \(x<\min\{\varphi ^{l}_{\varepsilon _{1}}(t), \varphi ^{l}_{\varepsilon _{2}}(t)\}\) or \(x>\max\{\varphi ^{l}_{\varepsilon _{1}}(t)+\varepsilon _{1},\varphi ^{l}_{\varepsilon _{2}}(t)+\varepsilon _{2}\}\). Thus by (3.21) we have
and
Summing equality (3.30) with respect to l from \(l=1\) to \(l=N\), using (3.31)-(3.35) and Minkowski’s inequality, and choosing \(\sigma=\nu/8\), we then conclude that
This, together with Gronwall’s inequality, yields
Combining the two inequalities above and (3.21) leads us to estimate (3.28). □
3.2 The solutions of the diffraction problem
Proof of Theorem 2.1
We divide the proof into three steps.
Step 1. We prove the global existence of the solutions. Let us discuss the behavior of the solution \(\mathbf {u}_{\varepsilon }\) associated with \(\boldsymbol {\varphi }_{\varepsilon }(t)\) by Theorem 2.1 as \(\varepsilon \to0\).
We first see from (3.16) that for any \(\tau\in[0,T]\) and any vector function ,
Furthermore, according to estimates (3.17), (3.20), (3.28) and the Arzela-Ascoli theorem, we conclude that there exists a subsequence (we retain the same notation for it) \(\{\mathbf {u}_{\varepsilon }\}\) such that
Thus u is in \(\mathcal{C}^{\alpha _{1},\alpha _{1}/2}(\bar{\mathcal {Q}}_{T})\), u satisfies the parabolic condition (1.2), and estimates (2.6) and (2.7) hold.
We next show that for each \(l=1,\ldots,N\), the sequences \(\{A^{l}_{\varepsilon }(x,t)u^{l}_{\varepsilon x}\}\), \(\{B^{l}_{\varepsilon }(x,t,u^{l}_{\varepsilon })u^{l}_{\varepsilon x}\}\), \(\{G_{\varepsilon }^{l}(x,t,\mathbf {u}_{\varepsilon })\}\) converge in \(L^{2}(Q_{T})\) to \(A^{l}(x,t)u^{l}_{x}\), \(B^{l}(x,t,u^{l}) u^{l}_{x}\) and \(G^{l}(x,t,\mathbf {u})\), respectively. Since
and
then it follows from (2.2), (3.3) and Lebesgue dominated convergence theorem that
and
and from (3.10) that
The similar argument shows that \(\|A^{l}_{\varepsilon }(x,t)u^{l}_{\varepsilon x}-A^{l}(x,t) u^{l}_{x}\|_{L^{2}(Q_{T})}\to0\) for each \(l=1,\ldots,N\).
Based on the above arguments for sequences \(\{A^{l}_{\varepsilon }(x,t)u^{l}_{\varepsilon x}\}\), \(\{B^{l}_{\varepsilon }(x,t,u^{l}_{\varepsilon })u^{l}_{\varepsilon x}\}\), \(\{G_{\varepsilon }^{l}(x,t,\mathbf {u}_{\varepsilon })\}\) and \(\{u^{l}_{\varepsilon t}\}\), by letting \(\varepsilon \to0\), we conclude from (3.36) that (2.4) holds.
We also see from (3.20) that there exists a subsequence \(\{\mathbf {u}_{\varepsilon '}\}\) (denoted by \(\{\mathbf {u}_{\varepsilon }\}\) still) such that for each \(l=1,\ldots,N\), \(\{(A_{\varepsilon }(x,t)u^{l}_{\varepsilon x})_{x}\}\) converge weakly in \(L^{2}(Q_{T})\) to \(\varpi^{l}(x,t)\). Recalling that \(A_{\varepsilon }(x,t)u^{l}_{\varepsilon x}\to A(x,t) u^{l}_{x}\) in \(L^{2}(Q_{T})\), we deduce \(\varpi_{l}=(A(x,t) u^{l}_{x})_{x}\). These, together with (3.20), imply that \(u^{l}_{xx}\in L^{2}(Q^{l-}_{T})\cap L^{2}(Q^{l+}_{T})\) for each \(l=1,\ldots,n\), \(u^{l}_{xx}\in L^{2}(Q_{T})\) for each \(l=n+1,\ldots,N\), \(A^{l}(x,t) u^{l}_{x}\in V_{2}(Q_{T})\) for each \(l=1,\ldots,N\), and (2.8) holds. Thus (2.4) implies that u satisfies the equations in (1.1a) and (1.1b) for almost all \((x,t)\in Q_{T}\) and the inner boundary condition (1.3) for almost all \(t\in[0,T]\) (see [17, Chapter 3, Section 13]). As we have done in the derivation of (3.21), estimate (2.8) yields (2.9).
For fixed \(l\in\{1,\ldots,N\}\) and \(k\in\{1,\ldots,n+1\}\), \(u^{l}\) satisfies the linear equation
where
and
Then for any subdomains \(Q'\) and \(Q''\) satisfying \(\bar{Q}\subset Q''\) and \(\bar{Q}''\subset Q_{k,T}\), we have \(\underline{b}^{l}(x,t), \underline{g}^{l}(x,t)\in C^{\alpha _{1},\alpha _{1}/2}(\bar{Q}'')\). The parabolic regularity theory shows that \(u^{l}\in C^{2+\alpha _{1},1+\alpha _{1}/2}(\bar{Q}')\). Hence u satisfies pointwise the equations in (1.1a) and (1.1b) for \((x,t)\in Q_{k,T}\). Consequently, u is a solution in of problem (1.1a)-(1.3) and estimates (2.6)-(2.9) hold.
Step 2. In what follows, we will show that the solution in for problem (1.1a)-(1.3) is unique and estimates (2.10) and (2.11) hold.
Let \(\mathbf {u}_{1}\), \(\mathbf {u}_{2}\) be the solutions in corresponding to \(\boldsymbol {\varphi }_{1}\) and \(\boldsymbol {\varphi }_{2}\), respectively. Set \(\mathbf {w}=\mathbf {u}_{1}-\mathbf {u}_{2}\). Then . We choose \(\boldsymbol {\eta }=\mathbf {w}\) in (2.4) to find
and
For each \(l=1,\ldots,N\), by a subtraction of the above equations for \(i=1,2\), we conclude that
and
Then
Setting \(\epsilon=1/2\) and summing the above inequalities with respect to l from \(l=1\) to \(l=N\), we have
Again by Gronwall’s inequality we deduce (2.10), which, together with [17, Chapter 2, formula (3.8)], gives (2.11). Therefore the solution in for problem (1.1a)-(1.3) associated with \(\boldsymbol {\varphi }(t)\) is unique.
Step 3. For \(\boldsymbol {\varphi }\in \mathbf {C}^{2}((0,T])\), we discuss the regularity of u.
For any fixed \(l\in\{1,\ldots,n\}\), let
Then v satisfies
We will use the result of [18] to obtain the regularity of v. To do this, we need the estimate of \(\operatorname{ess\,sup}_{t_{1}\leq t\leq T}\int_{-\delta/2}^{\delta/2}(v_{t})^{2}\, \mathrm{d}y\) for any fixed \(t_{1}\in(0,T)\). Let \(\xi=\xi(y,t)\) be a smooth function with values between 0 and 1 such that \(\xi=0\) for \(|y|\geq \delta\) or \(t\leq t_{1}/2\), \(\xi=1\) for \((y,t)\in[-\delta/2,\delta/2 ]\times[t_{1},T]\) and \(|\xi_{y}|+|\xi_{t}|\leq C(1/\delta,1/t_{1})\), and let
For any small enough Δt, consider the equality \(\int_{t_{1}/2}^{\tau} \int_{-\delta}^{\delta}[v_{t}-\mathcal {H}(y,t,v)]_{(t)}v_{(t)}(\xi)^{2}\,\mathrm{d}x\,\mathrm{d}t\). By employing the formula of integration by parts and [17, Chapter II, formula (4.7)], we get
where \(v_{y(t)}=(v_{y})_{(t)}\). Some tedious computation and Cauchy’s inequality yield
We choose \(\epsilon=\nu_{0}/4\) and employ (2.8) to find
As we have done in the derivation of (3.20), by Gronwall’s inequality we get
Consequently,
By [17, Chapter II, Lemma 4.11], this inequality implies that
Hence \(u^{l}_{xt}\in\L^{2}(Q^{l}(\delta,t_{1}))\). Using (3.37), hypothesis (H)(iii) and [18, Theorem 1.1], we deduce that \(v_{y}(y,t)\) is continuous with respect to y in \([-\delta/4,0]\) and in \([0, \delta/4]\) for almost all \(t\in[t_{1},T]\), and \(v_{t}\) is continuous in \([-\delta/4,\delta/4]\times[t_{1},T]\). Since \(u^{l}_{x}(x,t)=v_{y}(x-\varphi ^{l}(t),t)\) and \(u^{l}_{t}(x,t)=v_{t}(x-\varphi ^{l}(t),t)-v_{y}(x-\varphi ^{l}(t),t)\varphi ^{l}_{t}(t)\), then for almost all \(t\in[t_{1},T]\), \(u^{l}_{x}(x,t)\) and \(u^{l}_{t}(x,t)\) are continuous with respect to x in \([\varphi ^{l}(t)-\delta/4,\varphi ^{l}(t)]\) and in \([\varphi ^{l}(t),\varphi ^{l}(t)+\delta/4]\). □
The following corollary follows directly from Theorem 2.1.
Corollary 3.3
Assume that \(\boldsymbol {\varphi }_{m}(t)\in \mathbf {C}^{2}([0,T])\), \(m=1,2,\ldots\) , satisfy (2.5) and the sequence \(\{\boldsymbol {\varphi }_{m}(t)\}\) converges in \(\mathbf {W}^{1}_{2}(0,T)\) to \(\boldsymbol {\varphi }(t)\). If \(\mathbf {u}_{m}\), u are the solutions in of (1.1a)-(1.3) corresponding to \(\boldsymbol {\varphi }_{m}\) and φ, respectively, then there exists a subsequence (we retain the same notation for it) \(\{\mathbf {u}_{m}\}\) such that
4 The free boundary problem
The goal of this section is to prove the local existence of solutions for problem (1.1a)-(1.4).
In this section, let
Then, for \(\mathbf {r}(t)\in\mathbb{B}\), we have
Thus \(\varphi ^{i}(t)=d^{i}+\int_{0}^{t}r^{i}(t)\,\mathrm{d}t\), \(i=1,\ldots,n\), satisfy
and
For \(\mathbf {r}(t)\in\mathbb{B}\), define
where \(\varphi ^{i}(t)=d^{i}+\int_{0}^{t}r^{i}(\tau)\,\mathrm {d}\tau\) (\(i=1,\ldots,n\)) and u is the solution in of problem (1.1a)-(1.3) associated with \(\boldsymbol {\varphi }(t)=(\varphi ^{1}(t),\ldots, \varphi ^{n}(t))\).
Proof of Theorem 2.2
By a slight modification, we use the methods and the framework of [3, 5] to prove Theorem 2.2. We only sketch the main steps and omit most of the detailed proofs.
Step 1. We will show that \(\mathcal{G}(\mathbb{B})\subseteq \mathbb{B}\).
Employing (4.1), (4.2), (2.9) and Hölder’s inequality, we get
and
Consequently, \(\mathcal{G}(\mathbb{B})\subseteq\mathbb{B}\).
Step 2. By (2.11) and a slight modification of the proof of [3, Proposition 3.4], we can prove that \(\mathcal{G}:\mathbb{B} \to\mathbb{B}\) is continuous and \(\mathcal{G}(\mathbb{B})\subset\mathbb{B}\) is precompact in \(\mathbf {L}^{3}(0,T)\).
Step 3. It follows from Schauder’s fixed point theorem that \(\mathcal{G}:\mathbb{B}\to\mathbb{B}\) has at least one fixed point \(\mathbf {r}(t)\). Then the curves \(\varphi ^{i}(t)=d^{i}+\int_{0}^{t}r^{i}(\tau)\,\mathrm{d}\tau\), \(i=1,\ldots,n\), and the vector function \(\mathbf {u}(x,t)\) associated with \(\boldsymbol {\varphi }(t)\) by Theorem 2.1 satisfy (1.4) for almost all \(t\in[0,T]\). Moreover, (1.4), (2.6) and (2.9) imply that \(\boldsymbol {\varphi }(t)\) is in \(\mathbf {W}^{1}_{4}(0,T)\).
Step 4. Using the regularity result of Theorem 2.1(iii) and the approximation result of Corollary 3.3, and calculating \(\int_{0}^{T}\frac{\mathrm{d}}{\mathrm{d}t}\int_{0}^{d} (A^{l}(x,t))^{3}(u^{l}_{x})^{4}\,\mathrm{d}x\,\mathrm{d}t\), we can conclude from the proof similar to that of [5, Lemma 2.5] that
In addition, for each \(l=n+1,\ldots,N\), \(u^{l}\) satisfies the linear equation
Estimate (4.3) yields \(\underline{f}^{l}(x,t)\in L^{4}(Q_{T})\). Thus, according to the local parabolic \(L^{p}\) estimates, we get \(u^{l}\in W^{2,1}_{4}(Q')\) for any interior domain \(Q'\) of \(Q_{T}\). This implies that
for some \(\alpha \in(0,1)\).
Step 5. By (4.3), (4.4) and hypothesis (H)(iii), a similar argument as that for [5, Theorem 2.6] shows that for each \(l=1,\ldots,n\), \(\varphi ^{l}(t)\in C^{2+\alpha }((0,T])\). Then (1.4) holds for \(t\in(0, T]\) and \(\{\mathbf {u}(x,t), \boldsymbol {\varphi }(t)\}\) is a local solution of problem (1.1a)-(1.4). □
Remark 4.1
Theorem 2.2 concerns the local existence of solutions for (1.1a)-(1.4). For the case \(n=1\), the global existence of solutions of (1.1a)-(1.4) can be proved by the same argument as that in [4]. For the case \(n>1\), the difficulty to prove the global existence of solutions may be that the free boundaries can intersect with each other.
5 An example
In this section, we give an example satisfying the conditions in hypothesis (H).
Consider the Lotka-Volterra reaction-diffusion model with two predators \(u^{1}\), \(u^{2}\), and one prey \(u^{3}\), where each of predators \(u^{1}\), \(u^{2}\) has two groups and the prey \(u^{3}\) is native and has only one group. Then the reaction-diffusion functions are defined by
where \(\gamma^{l}_{k}\), \(\beta ^{l}_{k}\), \(\varsigma^{l}_{k}\), \(\sigma^{l}_{k}\) (\(l,k=1,2\)), \(\gamma^{3}\), \(\beta ^{3}\), \(\varsigma^{3}_{3}\) and \(\sigma^{3}_{3}\) are all positive constants. For the mankind’s influence, we assume that
where \(\hat{\chi}^{l}\), \(\tilde{\chi}^{l}\) (\(l=1,2\)) and \(\chi^{3}\) are constants. Thus \(\mathbf {u}=(u^{1},u^{2},u^{3})\) and \(\boldsymbol {\varphi }=(\varphi ^{1},\varphi ^{2})\) are governed by system (1.1a)-(1.4), where \(n=2\), \(N=3\) and the coefficients are replaced by (5.1), (5.2).
Assume that
and
We will look for \(\mathbf{M}=(M^{1},M^{2},M^{3})\) and \(\mathbf{m}=(m^{1},m^{2},m^{3})\) satisfying (2.1). This requires that
and
Choose \(m^{1}=m^{2}=m^{3}\) and \(M^{1}=M^{2}\). It is obvious that (5.6) and (5.7) hold if
and
By (5.3)-(5.5), we can first choose \(m^{1}=m^{2}=m^{3}\) small enough such that
and
Furthermore, choose \(M^{1}\), \(M^{2}\) and \(M^{3}\) such that
and
The relations (5.8)-(5.10) are guaranteed by the above M and m. Hence M and m satisfy (5.6)-(5.8); and consequently, M and m satisfy (2.1). Moreover, \(\mathbf{g}_{k}(\mathbf {u})=(g^{1}_{k}(\mathbf {u}),g^{2}_{k}(\mathbf {u}),g^{3}_{k}(\mathbf {u}))\), \(k=1,2,3\), are mixed quasimonotone in with the same index vector \(((u^{3});(u^{3});[\mathbf {u}]_{0})\). The above argument shows that the conditions in hypothesis (H) can be satisfied. Then we get the following corollary from Theorem 2.2.
Corollary 5.1
In problem (1.1a)-(1.4), let \(n=2\), \(N=3\), and let the coefficients be given by (5.1), (5.2). Assume that the conditions in (5.3)-(5.5) and hypothesis (H)(i) are satisfied. Then there exist \(T>0\) and a pair of functions \(\mathbf {u}(x,t)\), \(\boldsymbol {\varphi }(t)\) such that \(\{\mathbf {u}(x,t), \boldsymbol {\varphi }(t)\}\) is a solution of (1.1a)-(1.4).
References
Kim, K, Lin, Z: A free boundary problem for a parabolic system describing an ecological model. Nonlinear Anal., Real World Appl. 10, 428-436 (2009)
Kamynin, LI: On the linear Verigin problem. Dokl. Akad. Nauk SSSR 150, 1210-1213 (1963) (in Russian)
Evans, LC: A free boundary problem: the flow of two immiscible fluids in a one-dimensional porous medium, I. Indiana Univ. Math. J. 26, 914-932 (1977)
Evans, LC: A free boundary problem: the flow of two immiscible fluids in a one-dimensional porous medium, II. Indiana Univ. Math. J. 27, 93-111 (1978)
Evans, LC, Friedman, A: Regularity and asymptotic behavior of two immiscible fluids in a one-dimensional porous medium. J. Differ. Equ. 31, 366-391 (1979)
Liang, J: The one-dimensional quasilinear Verigin problem. J. Partial Differ. Equ. 4, 74-96 (1991)
Mimura, M, Yamada, Y, Yotsutani, S: A free boundary problem in ecology. Jpn. J. Appl. Math. 2, 151-186 (1985)
Ling, Z, Tang, Q, Lin, Z: A free boundary problem for two-species competitive model in ecology. Nonlinear Anal., Real World Appl. 11, 1775-1781 (2010)
Lin, Z: A free boundary problem for a predator prey model. Nonlinearity 20, 1883-1892 (2007)
Guan, Z: A class of free boundary problems to quasi-linear parabolic equations. Acta Math. Appl. Sin. 15, 197-205 (1999)
Yi, FH: Global classical solution of Muskat free boundary problem. J. Math. Anal. Appl. 288, 442-461 (2003)
Merz, W, Rybka, P: A free boundary problem describing reaction-diffusion problems in chemical vapor infiltration of pyrolytic carbon. J. Math. Anal. Appl. 292, 571-588 (2004)
Yi, F: One dimensional combustion free boundary problem. Glasg. Math. J. 46, 63-75 (2004)
Boyadjiev, G, Kutev, N: Diffraction problems for quasilinear reaction-diffusion systems. Nonlinear Anal. 55, 905-926 (2003)
Tan, QJ: Systems of quasilinear parabolic equations with discontinuous coefficients and continuous delays. Adv. Differ. Equ. 2011, Article ID 925173 (2011). doi:10.1155/2011/925173
Tan, QJ, Pan, CY: Diffraction problems for quasilinear parabolic systems with boundary intersecting interfaces. Bound. Value Probl. 2013, Article ID 99 (2013)
Ladyzenskaya, OA, Solonnikov, VA, Ural’ceva, NN: Linear and Quasilinear Equations of Parabolic Type. Am. Math. Soc., Providence (1968)
Tan, QJ, Leng, ZJ: The regularity of generalized solutions for the n-dimensional quasi-linear parabolic diffraction problems. Nonlinear Anal. 69, 4624-4642 (2008)
Acknowledgements
The author would like to thank the reviewers and the editors for their valuable suggestions and comments. The work was supported by the research fund of Department of Education of Sichuan Province and the research fund of Chengdu Normal University (CS14ZD01).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The author declares that they have no competing interests.
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
Tan, QJ. A free boundary problem arising in the ecological models with N-species. Bound Value Probl 2015, 4 (2015). https://doi.org/10.1186/s13661-014-0267-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-014-0267-3