# New-Type Solutions of the Modified Fischer-Kolmogorov Equation

Authors: P. Drábek
Journal: Abstract and Applied Analysis
 (2011)

---

## Abstract

We prove the existence of new-type solutions of the modified
Fischer-Kolmogorov equation with slow/fast diffusion and with possibly nonsmooth double-well potential. We show that a certain relation
between the rate of the diffusion and the smoothness of the potential may originate new type solutions which do not occur in the classical Fischer-Kolmogorov equation. The main focus of this paper is to show the sensitivity of the mathematical modelling with respect to the chosen form of the diffusion term and the shape of the double-well potential.


---

## Body

## 1. Introduction

The Fischer-Kolmogorov (FK) equation(1.1)∂u∂t=∂2u∂x2+u-u3
was originally proposed in 1937 to model the interaction of dispersal and fitness in biological populations; see [1]. In connection with the modelling of an active pulse transmission line simulating a nerve axon, it is called the Newell-Whitehead equation or the Nagumo equation; see [2]. The stationary solutions of FK equation are solutions of the second-order equation(1.2)u′′+u-u3=0.This equation has a first integral, and we can define the energy functional(1.3)E(u)=def-12(u′)2+14(1-u2)2,
which is constant along solutions of (1.2).Since (1.2) is equivalent to the first-order system(1.4)u′=v,v′=-u+u3,
solutions of (1.2) can be depicted as orbits of the system (1.4) in the (u,u′)-plane; see [1, page 11]. Complete classification of solutions of (1.2) is given in [1, pages 10–12].Let us focus on the solutionsu(x) of (1.2) with zero energy(1.5)E(u(x))=0,x∈R.
It is easy to see that the only solutions of (1.2) satisfying (1.5) are(1)
constant solutions:u(x)≡-1 and u(x)≡1,(2)
kinks(called also domain walls): heteroclinic orbits that connect the saddle points (u,u′)=(±1,0). These solutions describe transition layers between the two stable uniform states −1 and +1 and have explicit form
(1.6)u(x)=±tanh(x-x02),
where x0∈ℝ is arbitrary but fixed number. We may assume, without loss of generality, that x0=0 due to the fact that (1.2) is autonomous.In contrast with some fourth-order problems, the dynamical system (1.4) associated with the second-order FK equation does not admit any homoclinic orbits, and so (1.2) does not possess any solutions which are called pulses (see [1, page 12]). In this paper, we want to show that this fact is a consequence of the combined effect of the linear diffusion term and the C2-regularity of the double-well potential (1.7)W(u)=def14(1-u2)2.It may well happen in practice that the diffusion is slower than the linear one. Then besides kinks, also manynew types of nonconstant stationary solutions on zero energy level may occur. These solutions might be of the shape similar to the kinks mentioned above, but the shape of solutions might be of pulse type as well. Moreover, one can also get periodic solutions or any random C1-transition between the stable states −1 and +1; see the next section. This may happen even in the case when the double-well potential is a smooth, C∞-function. Similar phenomenon occurs also when the diffusion is expressed by the linear Laplace operator and the double-well potential lacks its C2-regularity at the global minimizers.On the other hand, if the slow diffusion is compensated by higher regularity of the double-well potential, then the kinks are the only nonconstant  stationary solutions on zero energy level. In general, we show that a certain relation between the rate of the diffusion and the smoothness of the potential may originate new type solutions which do not occur in the classical Fischer-Kolmogorov equation or if such a relation does not hold true.Themain focus of this paper is to show the sensitivity of the mathematical modelling with respect to the chosen form of the diffusion term and the shape of the double-well potential as well.


## 2. Modified FK Equation

We consider modified FK equation(2.1)∂u∂t=∂∂x(|∂u∂x|p-2∂u∂x)+u(1-u2)α-1,
where α>1 and p>1 are real parameters. Obviously, for α=p=2, (2.1) reduces to FK equation. The one-dimensional p-Laplacian (∂/∂x)(|∂u/∂x|p-2(∂u/∂x)) corresponds to the fast or slow diffusion if p∈(1,2) or p∈(2,∞), respectively. The double-well potential (2.2)Wα(u)=def12α(1-u2)α
keeps the original shape of W=W2, but it exhibits the lack of C2-smoothness at the global minimizers ±1 if α∈(1,2).A counterpart of (1.2) is the quasilinear elliptic second-order equation (2.3)(|u′|p-2u′)′+u(1-u2)α-1=0,
and the energy functional is given by (2.4)Eα,p(u)=def-1p|u′|p+12α(1-u2)α.
It is clear that u(x)≡±1 are constant solutions of (2.3) which correspond to the zero energy ℰα,p(u(x))=0, x∈ℝ. We observe that for any solution corresponding to the zero energy we have u′(x)=0 if and only if either u(x)=-1 or else u(x)=+1. It then follows that a solution corresponding to the zero energy is either monotone or constant (equal to either −1 or +1) on suitable subintervals of ℝ. In particular, we also observe that a bounded monotone solution of (2.3) corresponding to the zero energy must be increasing or decreasing transition between two stable uniform states −1 and +1.Letu(x) be such a transition. Since (2.3) is autonomous, we may assume that u(0)=0. Without loss of generality, we may also assume that u′(0)=(p/2α)1/p due to the uniqueness theorem for the initial value problem, see, for example, [3]. Then there is x0, 0<x0≤∞ such that 1>u(x)>0, u′(x)>0, x∈(0,x0). Hence, there exists a positive inverse function x(u) of u(x) on (0,x0). We then introduce new variables t=u and z(t)=(u′(x(t)))p, compare with [1, page 53]. Upon the substitution into the first integral (2.5)Eα,p(u)=-1p|u′|p+12α(1-u2)α=0,
we obtain (2.6)-1pz(t)+12α(1-t2)α=0,
that is, (2.7)z(t)=p2α(1-t2)α.
Hence, (2.8)x(=x(u))=∫0u(x)dt(z(t))1/p=(2αp)1/p∫0u(x)dt(1-t2)α/p,x∈[0,x0).
In particular, for x↗x0, u(x)↗1, and so we get (2.9)x0=(2αp)1/p∫01dt(1-t2)α/p.
It then follows from (2.9) that  (1)
x0=∞ if α≥p,(2)
0<x0<∞ if α<p.In the first case, we thus have strictly increasing solutionu(x) on [0,∞) such that u(x)↗1 as x↗∞. We extend u to (-∞,0) as u(x)=-u(-x), x∈(-∞,0), and by the uniqueness theorem, we get a global solution of (2.3) which is defined on the entire ℝ. It is actually an increasing odd kink of the shape similar to that given by u(x)=tanh(x/2) (which occurs in the case α=p=2). Similarly, we obtain a decreasing odd kink, and all other bounded nonconstant  solutions differ from these two kinks just modulo translations; see Figure 1.Figure 1We want to emphasize that the situation is completely different in the second case. A solutionu(x) of (2.3) on [0,x0) is given implicitly by (2.8). We have u(0)=0, u′(0)=(p/2α)1/p, and we can define u(x0)=1, u′(x0)=0. We then extend u(x) as follows: u(x)=-u(-x), x∈[-x0,0), u(x)=-1, x∈(-∞,-x0), u(x)=1, x∈(x0,∞). We obtain a global solution of (2.3) which is of kink type; see Figure 2. It differs from the previous increasing odd kink by the property that the transition from −1 to +1 occurs on the finite interval [-x0,x0].Figure 2In the language of the dynamical system(2.10)u′=|v|p′-2v,v′=u(1-u2)α-1,p′=p/(p-1), associated with (2.3), this solution is a heteroclinic orbit that connects the saddle points (u,u′)=(±1,0). In contrast with the case α≥p, this orbit connects the equilibria in a “finite time.”Notice that we have a plenty of other possibilities how to extendu(x) given by (2.8) from [0,x0] to get a global solution of (2.3). In particular, we can extend u(x) as follows: u(x)=-u(-x), x∈[-x0,0), u(x)=u(2x0-x), x∈(x0,3x0], u(x)=-1, x∈(-∞,-x0)∪(3x0,∞); see Figure 3. We obtain a global solution of (2.3) which is of pulse type. It differs from the pulses defined in [1] by the property that the solution leaves one stable uniform state (-1), it touches the other one (+1), and it returns back to the original stable state (-1). Both transitions occur on finite intervals of the same length. Such a solution is a homoclinic orbit of (2.10) which emanates from (-1,0), passes through (+1,0), and returns back to (-1,0) in a “finite time.”Figure 3Similarly, we can construct infinitely many (even continua) of stationary solutions of (2.3) on zero energy level. Some of them might be periodic; the other ones are just random transitions between uniform steady states (-1) and (+1); see Figures 4 and 5.Figure 4Figure 5


## 3. Discussions

We would like to point out that the same results hold true if the double-well potentialWα is replaced by more general nonnegative even function F(u) with a positive local maximum at 0 and two global minima at ±1, F(±1)=0, F strictly concave in (-ξ,ξ), strictly convex in (-∞,-ξ)∪(ξ,∞) for some ξ∈(0,1), F∈C2 with possible exception at ±1, where we might have F∈C1 only. Then the analogues of (2.8) and (2.9) read as (3.1)x=∫0u(x)dt(F(t))1/p,x∈[0,x0),x0=∫01dt(F(t))1/p,
respectively, and we have either  (1)
x0=∫01(dt/(F(t))1/p)=∞,   or else,(2)
0<x0=∫01(dt/(F(t))1/p)<∞.The only bounded nonconstant  solutions in the first case are kinks. The existence of new-type solutions in the second case is then related to the nonuniqueness of the solution of the initial value problem(3.2)(|u′|p-2u′)′-F′(u)=0,u(x0)=±1,u′(x0)=0.It is clear from Section2 that the existence of new-type solutions depends on the relation between the rate of the diffusion and the regularity of the double-well potential at its global minimizers. In particular, equations (3.3)∂u∂t=∂2u∂x2+u(1-u2)α-1
with 1<α<2 and (3.4)∂u∂t=∂∂x(|∂u∂x|p-2∂u∂x)+u-u3
with p>2 exhibit the existence of new-type stationary solutions.The generalization of these results to reaction-diffusion type equations in higher space dimensions is far from being obvious. For example, in the radial symmetric case, the counterpart of (2.3) is not an autonomous equation. Hence, we cannot use similar calculations with a first integral, and continua of new solutions cannot be constructed simply by graphing transitions between steady stable states as in one-dimensional case. Nevertheless, the author conjectures that even in higher space dimensions new-type solutions occur under suitable conditions on p and/or α. To prove this fact and to classify these solutions might be an interesting open problem.



---

*Source: 247619-2011-06-14.xml*