Skip to main content

Positive blow-up solutions of nonlinear models from real world dynamics

Abstract

In this paper, we investigate the structure and properties of the set of positive blow-up solutions of the differential equation ( t k v ( t ) ) = t k h(t,v(t)), t(0,T]R, where k(1,). The differential equation is studied together with the boundary conditions lim t 0 + v(t)=, v(T)=0. We specify conditions for the data function h which guarantee that the set of all positive solutions to the above boundary value problem is nonempty. Further properties of the solutions are discussed and results of numerical simulations are presented.

MSC:34B18, 34B16, 34A12.

1 Introduction

In this paper, we investigate the structure and properties of the set of positive blow-up solutions of the differential equation

( t k v ( t ) ) = t k h ( t , v ( t ) ) ,t(0,T]R,
(1)

where k(1,).

Models in the form of (1) arise in many applications. Among others, they occur in the study of phase transitions of Van der Waals fluids [13], in population genetics, where they characterize the spatial distribution of the genetic composition of a population [4, 5], in the homogeneous nucleation theory [6], in relativistic cosmology for particles which can be treated as domains in the universe [7], and in the nonlinear field theory, in particular, in context of bubbles generated by scalar fields of the Higgs type in Minkowski spaces [8].

Here, we assume that h is positive and satisfies the Carathéodory conditions on [0,T]×[0,). We define a positive solution of (1) as a function v which satisfies (1) for a.e. t[0,T], is positive on (0,T), and has absolutely continuous first derivative on each compact subinterval in (0,T].

According to Lemma 4, if v is a positive solution of (1), then either

lim t 0 + v(t)[0,), lim t 0 + v (t)=0,

or

lim t 0 + v(t)=, lim t 0 + v (t)=.

In the literature, bounded solutions of (1) have been widely investigated; see e.g. [913]. Such solutions are characterized by the initial condition lim t 0 + v (t)=0. In contrast to this, some real problems lead to the investigation of unbounded solutions which are characterized by the condition lim t t v(t)= for some t [0,T] and which are called blow-up solutions. We refer to [1416]. Here, we are interested in blow-up solutions of (1), where t =0. In particular, (1) will be considered together with the boundary conditions

lim t 0 + v(t)=,v(T)=0.
(2)

In this case, we speak about a positive solution of problem ( 1 ), ( 2 ). Let us denote by Z the set of all positive solutions to (1), (2). Moreover, let R={vZ: t k v is bounded on (0,T]} and R c ={vR: v (T)=c} for c0. Since v (T)0 for each vZ, it is obvious that R= c 0 R c .

Our main goal is to find conditions for the data function h in (1), which guarantee that the set R c is nonempty for each c0 and then to investigate the properties of this set. For example, we prove that the difference of any two functions in R c , c0, retains its sign on [0,T), and that there exist minimal and maximal solutions v c , min , v c , max R c for each c0, cf. Theorem 5. If the interior of the set {(t,x) R 2 :0tT, v c , min (t)x v c , max (t)} is nonempty, we show that this interior is fully covered by ordered graphs of other functions belonging to R c for each c0; see Theorem 6. Finally, in Theorem 7, the existence of a positive constant α 0 such that α 0 = lim t 0 + t k 1 v 0 , min (t) is shown and all functions v from are uniquely characterized by the condition

lim t 0 + t k 1 v(t)=α,α[ α 0 ,).

If we denote such v by v α and define w α (t)= t k 1 v(t) for t(0,T], and w α (0)=α, then we find that the graphs of these functions do not intersect, cf. Theorem 8, and that for each α > α 0 , the set { w α : α 0 α α } is compact in C[0,T]; see Theorem 9.

The study of a structure of positive solutions to other types of ordinary differential equations can be found for example in [1719].

Notation

Let us denote by C[0,T] the Banach space of functions continuous on [0,T] equipped with the maximum norm

x =max { | x ( t ) | : t [ 0 , T ] } .

Similarly, C 1 [0,T] means the Banach space of functions having a continuous first derivative on [0,T] with the corresponding maximum norm x + x . By L 1 [0,T] we denote the set of functions which are Lebesgue integrable on [0,T]. Moreover, A C 1 [0,T] is the set of functions whose first derivative is absolutely continuous on [0,T], while A C loc 1 (0,T] is the set of functions having absolutely continuous first derivative on each compact subinterval of (0,T]. We say that h(t,x) satisfies the Carathéodory conditions on [0,T]×[0,), if the following three conditions hold:

  1. (i)

    The function h(,x):[0,T]R is measurable for all x[0,).

  2. (ii)

    The function h(t,):[0,)R is continuous for a.e. t[0,T].

  3. (iii)

    For each compact set U[0,) there exists a function m U L 1 [0,T] such that

    |h(t,x)| m U (t)for a.e.  t [ 0 , T ]  and all  x U .

For functions satisfying above conditions, we use the notation hCar([0,T]×[0,)).

Structure of the paper

The paper is organized as follows. In Section 2 we discuss properties of the solutions of the auxiliary Dirichlet problem (3), (4). We recapitulate previous results from [20] and also present new results in Theorems 1, 2, and 3. The main results of the paper can be found in Section 3, where we describe a relation between solutions of problem (3), (4) and blow-up solutions of problem (1), (2); see Theorem 4. Using this relation and the results of Section 2, we obtain various interesting properties of blow-up solutions; see Theorems 5 to 9. Section 4 contains three examples illustrating the theoretical findings. Final remarks and open problems are formulated in Section 5.

2 Auxiliary results

In this section we consider the auxiliary singular differential equation,

u (t)+ a t u (t) a t 2 u(t)=f ( t , u ( t ) ) ,
(3)

where a<1. For k:=a and v(t):= t k u(t), (3) becomes a special case of (1) and therefore, results obtained for (3) apply for (1). For the further analysis, we assume that f satisfies the following conditions:

(H1) fCar([0,T]×[0,)),

(H2) 0<f(t,x) for a.e. t[0,T] and all x[0,),

(H3) f(t,x) is increasing in x for a.e. t[0,T] and

lim x 1 x 0 T f(t,x)dt=0.

We now study (3) subject to the boundary conditions

u(0)=0,u(T)=0,
(4)

and require that

u (T)=c,c0,
(5)

holds. We call a function u:[0,T]R a positive solution of the Dirichlet problem (3), (4) if uA C 1 [0,T], u>0 on (0,T), u satisfies the boundary conditions (4), and (3) holds for a.e. t[0,T]. Clearly, for each positive solution u of (3), (4), there exists c0 such that (5) is satisfied.

We now denote by S the set of all positive solutions of problem (3), (4), and S c ={uS: u (T)=c}, where c0.

In the following lemma we cite those results from [20] which will be used in the analysis of problem (1), (2).

Lemma 1 Let (H1)-(H3) hold. Then the following statements hold:

  1. (a)

    For each c0 the set S c is nonempty and there exist functions u c , min , u c , max S c such that u c , min (t)u(t) u c , max (t) for t[0,T] and u S c .

  2. (b)

    If c 1 > c 2 0, u i S c i , i=1,2, then u 1 (t)> u 2 (t) for t(0,T).

  3. (c)

    If c0, u i S c , i=1,2, and u 1 ( t 0 )> u 2 ( t 0 ) for some t 0 (0,T), then either u 1 (t)> u 2 (t) for t(0,T) or there exists t ( t 0 ,T) such that u 1 (t)> u 2 (t) for t(0, t ) and u 1 (t)= u 2 (t) for t[ t ,T].

  4. (d)

    For each t 0 (0,T) and each A> u 0 , max ( t 0 ) there exists uS satisfying u( t 0 )=A.

  5. (e)

    S c is a one-point set for each c[0,)Γ, where Γ[0,) is at most countable.

  6. (f)

    For each 0KQ<, the set K c Q S c is compact in C 1 [0,T].

  7. (g)

    If c0, then u S c if and only if it is a solution of the equation

    u(t)=t c T a + 1 | a + 1 | ( T a 1 t a 1 ) +t t T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) ds
    (6)

in the set C 1 [0,T].

We now formulate new results which complete those from [20]. We first establish a relation between S 0 and the set {(t,x) R 2 :0tT, u 0 , min (t)x u 0 , max (t)} if its interior is nonempty. This question was a short time ago an open problem [[20], Remark 4.4]. We note that the relation between S c and the set {(t,x) R 2 :0tT, u c , min (t)x u c , max (t)} with c>0 having nonempty interior is described in Lemma 1(d).

Theorem 1 Let (H1)-(H3) hold. Let us assume that there exists t 0 (0,T) such that u 0 , min ( t 0 )< u 0 , max ( t 0 ). Then for any Q( u 0 , min ( t 0 ), u 0 , max ( t 0 )) there exists u S 0 satisfying u( t 0 )=Q.

Proof Step 1. Auxiliary Dirichlet problem.

Choose Q( u 0 , min ( t 0 ), u 0 , max ( t 0 )). Consider (3) subject to the Dirichlet conditions

u( t 0 )=Q,u(T)=0.
(7)

We claim that there exists a solution v to problem (3), (7) such that

u 0 , min (t)v(t) u 0 , max (t)for  t [ t 0 , T ] .
(8)

We show this result utilizing the method of lower and upper functions. It follows from (H1) that there exists φ L 1 [ t 0 ,T] such that

| a t 2 x+f(t,x)|φ(t)for a.e. t[ t 0 ,T] and  u 0 , min (t)x u 0 , max (t).
(9)

Let

W= t 0 T φ(t)dt,S=max { u 0 , min , u 0 , max , ( T t 0 ) | a | W } +1,

and let χ:R[0,1] be given as

χ(x)={ 1 if  | x | S , 2 | x | S if  S < | x | < 2 S , 0 if  | x | 2 S .

We now consider the auxiliary differential equation

v (t)=χ ( v ( t ) ) ( a t v ( t ) + a t 2 v ( t ) + f ( t , v ( t ) ) ) .
(10)

It is not difficult to verify that χ(y)|y|S for yR. Hence, cf. (9),

χ(y)| a t y+ a t 2 x+f(t,x)|χ(y) | a y | t +χ(y)φ(t) | a | S t +φ(t)
(11)

for a.e. t[ t 0 ,T] and all u 0 , min (t)x u 0 , max (t), yR. Since χ( u 0 , min (t))=1 and χ( u 0 , max (t))=1 for t[ t 0 ,T], u 0 , min (T)= u 0 , max (T)=0, u 0 , min ( t 0 )<Q< u 0 , max ( t 0 ), and u 0 , min , u 0 , max solve (3) on [0,T], we conclude that u 0 , min and u 0 , max are lower and upper functions of problem (3), (7) (see e.g. [21] or [22]). This fact together with (11) implies the existence of a solution v to problem (3), (7) satisfying (8), cf. [[21], Lemma 3.7]. Moreover,

u 0 , min ( t ) T t v ( t ) T t u 0 , max ( t ) T t ,t[ t 0 ,T),

and taking the limit tT we obtaina u 0 , min (T) v (T) u 0 , max (T), which together with u 0 , min (T)= u 0 , max (T)=0 gives v (T)=0.

We now prove that | v |<S on [ t 0 ,T]. The proof is indirect. Let us assume that there exists ξ[ t 0 ,T) such that | v (ξ)|=S and | v (t)|<S for t(ξ,T]. Then χ( v (t))=1 on this interval, and therefore

v (t)= a t v (t)+ a t 2 v(t)+f ( t , v ( t ) ) for a.e. t[ξ,T].

From (8) and (9) we conclude that

| a t 2 v(t)+f ( t , v ( t ) ) |φ(t)for a.e. t[ξ,T].

Since

v ( t ) = v ( T ) v ( t ) = t T v ( s ) d s = a t T v ( s ) s d s + t T ( a s 2 v ( s ) + f ( s , v ( s ) ) ) d s ,

we have

| v (t)||a| t T | v ( s ) | s ds+ t T φ(s)ds.

In particular,

| v (t)||a| t T | v ( s ) | s ds+Wfor t[ξ,T].

By the Gronwall lemma,

| v (t)|Wexp ( | a | t T d s s ) =W ( T t ) | a | ,t[ξ,T].

Therefore

S=| v (ξ)|W ( T ξ ) | a | ,

which contradicts

W ( T ξ ) | a | W ( T t 0 ) | a | S1.

Consequently, | v (t)|<S for t[ t 0 ,T], and so χ( v ) = 1 on this interval. Thus, v is a solution of problem (3), (7).

Step 2. Continuation of the solution v.

It follows from the arguments given in Step 1 that v is a solution of problem (3), (7) on [ t 0 ,T], v (T)=0 and u 0 , min ( t 0 )<v( t 0 )< u 0 , max ( t 0 ) because v( t 0 )=Q( u 0 , min ( t 0 ), u 0 , max ( t 0 )). It is easy to verify that the equality

( t a + 2 ( v ( t ) t ) ) = t a + 1 f ( t , v ( t ) )
(12)

is satisfied for a.e. t[ t 0 ,T]. We now integrate the last equality two times over [t,T][ t 0 ,T] and have (note that v(T)= v (T)=0)

v(t)=t t T s a 2 ( s T ξ a + 1 f ( ξ , v ( ξ ) ) d ξ ) ds,t[ t 0 ,T].

Let u be a solution of problem (3), (7) on an interval J which is a left-continuation of v. Let us assume that u is not continuable. Let γ=inf{t:tJ}. Then 0γ< t 0 and (12) with v replaced by u holds for a.e. tJ. The integration now yields

u(t)=t t T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) ds,tJ,
(13)

and we claim that

u 0 , min (t)<u(t)< u 0 , max (t)for t(γ, t 0 ].
(14)

We show inequality (14) indirectly. Let us assume that (14) does not hold. Then there exists ν(γ, t 0 ) such that

u 0 , min (t)<u(t)< u 0 , max (t)for t(ν, t 0 ]

and either u(ν)= u 0 , max (ν) or u(ν)= u 0 , min (ν). Assume that u(ν)= u 0 , max (ν).

Since, by (H3), f(t,u(t))<f(t, u 0 , max (t)) for a.e. t[ν, t 0 ] and f(t,u(t))f(t, u 0 , max (t)) for a.e. t[ t 0 ,T], we have

u ( ν ) = ν ν T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) d s < ν ν T s a 2 ( s T ξ a + 1 f ( ξ , u 0 , max ( ξ ) ) d ξ ) d s = u 0 , max ( ν ) ,

which is not possible. The case u(ν)= u 0 , min (ν) can be discussed analogously. Hence, (14) holds.

Suppose that γ>0. Then J=(γ,T] and since u is bounded on (γ,T], it follows that lim sup t γ + | u (t)|=. The integration of the equality

u (t)= a t u (t)+ a t 2 u(t)+f ( t , u ( t ) ) for a.e. t(γ,T]

over [t,T](γ,T] gives

u (t)=a t T u ( s ) s ds+ t T ( a s 2 u ( s ) + f ( s , u ( s ) ) ) ds,t(γ,T].

Since | a t 2 u(t)+f(t,u(t))|p(t) for a.e. t[γ,T], where p L 1 [γ,T], we have

| u (t)||a| t T | u ( s ) | s ds+ γ T p(s)ds,t(γ,T].

Using the Gronwall lemma, we deduce that for t(γ,T],

| u (t)| ( T t ) | a | γ T p(s)ds ( T γ ) | a | γ T p(s)ds

holds. This is a contradiction.

Therefore γ=0, and then (14) yields lim t 0 + u(t)=0. Consequently J=[0,T], and the assertion of the theorem follows from (13). □

In the next corollary we extend the statement (d) from Lemma 1 to a large set of A values.

Corollary 1 For each t 0 (0,T) and each A u 0 , min ( t 0 ) there exists uS satisfying u( t 0 )=A.

Proof The result follows immediately from Lemma 1(d) and Theorem 1. □

Remark 1 Corollary 1 says that the set U={(t,x) R 2 :t[0,T],x u 0 , min (t)} is covered by the graphs of functions from S, that is,

U= { ( t , u ( t ) ) : t [ 0 , T ] , u S } .

By Lemma 1(b), (c) we know that functions from S are uniquely determined by the values −c of their derivatives at the right end point t=T only in the case that S c is a singleton set for each c0. Since we cannot uniquely determine all functions from S via their derivatives at t=T, see Lemma 1(e), we discuss their derivatives at the singular point t=0.

Lemma 2 Let (H1)-(H3) hold. Let u S c , c0. Then

u (0)= c | a + 1 | + 0 T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) ds.
(15)

Proof It follows from Lemma 1(g) that (6) holds for t[0,T]. Since u(0)=0, we have (note that a+1<0)

u ( 0 ) = lim t 0 + u ( t ) u ( 0 ) t = lim t 0 + ( c T a + 1 | a + 1 | ( T a 1 t a 1 ) + t T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) d s ) = c | a + 1 | + 0 T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) d s ,

and (15) follows. □

Corollary 2 Let uS. Then

u (0)= u ( T ) | a + 1 | + 0 T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) ds.

Proof The result follows from Lemma 2, since u S u ( T ) . □

Corollary 3 Let u,vS, uv. Then u (0) v (0).

Proof Since uv, there exists t 0 (0,T) such that u( t 0 )v( t 0 ). We can assume that for instance u( t 0 )<v( t 0 ). Then u<v on (0, t 0 ] and uv on ( t 0 ,T] by Lemma 1(b)(c). Hence, u (T) v (T) and

0 T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) d s < 0 T s a 2 ( s T ξ a + 1 f ( ξ , v ( ξ ) ) d ξ ) d s ,

which together with Corollary 2 gives u (0)< v (0). □

Corollary 4 Let u S c , c0. Then

u (0)> c | a + 1 | .
(16)

In particular, u 0 , min (0)>0.

Proof Inequality (16) follows from Lemma 2 and the fact that f(t,u(t))>0 for a.e. t[0,T] by (H2). □

Corollary 5 Let uS. Then

0 T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) ds= 1 | a + 1 | 0 T f ( s , u ( s ) ) ds.

Proof The integration of (3) over [ε,T](0,T] gives

u (T) u (ε)+a ε T u ( t ) t dta ε T u ( t ) t 2 dt= ε T f ( t , u ( t ) ) dt.

Using integration by parts we obtain

ε T u ( t ) t 2 dt= u ( ε ) ε + ε T u ( t ) t dt,

and, therefore,

u (T) u (ε)a u ( ε ) ε = ε T f ( t , u ( t ) ) dt.

Taking the limit ε0, we have

u (T)+|a+1| u (0)= 0 T f ( t , u ( t ) ) dt.

On the other hand, it follows from Corollary 2 that

u (T)+|a+1| u (0)=|a+1| 0 T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) ds.

Combining the above two equalities yields the result. □

Lemma 3 Let (H1)-(H3) hold. Let c0 and S c be not a singleton set. Then for each α[ u c , min (0), u c , max (0)] there exists a unique u S c such that u (0)=α. Consequently, [ u c , min (0), u c , max (0)]={ u (0):u S c }.

Proof Let u S c , u c , min u u c , max . Then it follows from Lemma 1(a) and (H3) that

0 T s a 2 ( s T ξ a + 1 f ( ξ , u c , min ( ξ ) ) d ξ ) d s < 0 T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) d s < 0 T s a 2 ( s T ξ a + 1 f ( ξ , u c , max ( ξ ) ) d ξ ) d s .

Hence, by Lemma 2,

u c , min (0)< u (0)< u c , max (0)for u S c , u c , min u u c , max .

Since S c is a compact set in C 1 [0,T] by Lemma 1(f), the set I={ u (0):u S c }[ u c , min (0), u c , max (0)] is closed. In fact, let { β n }I, β n = u n (0), where u n S c , and let lim n β n =β. Then there exist a subsequence { β k n } of { β n } and u S c such that lim n u k n =u in C 1 [0,T]. In particular, β k n = u k n (0) u (0) as n. Therefore, β= u (0)I.

It remains to prove that I=[ u c , min (0), u c , max (0)]. Assume that the equality does not hold. Then, from the structure of bounded and closed subsets of the existence of an open interval (α,β)[ u c , min (0), u c , max (0)]I, α,βI, follows. Let α= u α (0) and β= u β (0), where u α , u β S c . Due to Lemma 1(c), there exists t 0 (0,T) such that u α < u β on (0, t 0 ], u α u β on ( t 0 ,T]. Choose τ( u α ( t 0 ), u β ( t 0 )). By Corollary 1, there exists v S c such that v( t 0 )=τ. Then u α <v< u β on (0, t 0 ] and u α v u β on ( t 0 ,T]. Consequently, u α (0)< v (0)< u β (0) on ( t 0 ,T], that is, v (0)(α,β), which is not possible. □

Since functions from S belong to A C 1 [0,T], we have u (0)< for each uS. Corollary 4 yields u 0 , min (0)>0. Let us denote

J:=[ u 0 , min (0),)(0,).
(17)

Lemma 3 implies that functions from S can be uniquely determined by the values of their derivatives at the singular point t=0; see Theorem 2.

Theorem 2 Let (H1)-(H3) hold. Then there exists a unique uS satisfying u (0)=α if and only if αJ.

Proof We first show

{ u ( 0 ) : u S } =J.
(18)

It follows from Lemma 3 that

{ u ( 0 ) : u 0 c K S c } = [ u 0 , min ( 0 ) , u K , max ( 0 ) ]

for each K0, where we set u c , min =u= u c , max if S c is a singleton set and S c ={u}. In view of Corollary 4, u (0)>c/|a+1| for u S c . Consequently, (18) follows.

Let us now choose αJ. Then there exists uS satisfying u (0)=α. The uniqueness of u follows from Corollary 3. □

For αJ, we denote by u α the unique element of S satisfying u α (0)=α.

Theorem 3 Let (H1)-(H3) hold. Assume that { α n }J is a convergent sequence and lim n α n =α. Then lim n u α n = u α in C 1 [0,T].

Proof Let K=sup{ u α n (T):nN}. Then K< because α n = u α n (0)> u α n (T)/|a+1| for nN by Corollary 4. Since { u α n } 0 c K S c and 0 c K S c is compact in C 1 [0,T] by Lemma 1(f), the sequence { u α n } is relatively compact in C 1 [0,T]. Let { u α n } be a subsequence of { u α n } which is convergent in C 1 [0,T], and let lim n u α n =u. Then u S c 0 for a c 0 [0,K] and u (0)= lim n u α n (0)= lim n α n =α. Therefore, u= u α and hence any subsequence { u α n } of { u α n } converging in C 1 [0,T] has the same limit u α . Consequently, { u α n } is convergent in C 1 [0,T] and u α is its limit. □

3 Blow-up solutions and their properties

This section contains the main results of the paper. First, we present a lemma which describes how positive solutions of (1) may behave at the singular point t=0.

Lemma 4 Let k(1,) and let the function hCar([0,T]×(0,)) be positive. Let v be a positive solution of (1). Then either

lim t 0 + v(t)[0,), lim t 0 + v (t)=0,
(19)

or

lim t 0 + v(t)=, lim t 0 + v (t)=.
(20)

Proof By Corollary 3.5 in [12], if sup{|v(t)|+| v (t)|:t(0,T]}<, then v satisfies (19). Now, assume that

sup { | v ( t ) | + | v ( t ) | : t ( 0 , T ] } =.
(21)

Then v has to satisfy

sup { | v ( t ) | : t ( 0 , T ] } =,
(22)

because otherwise, if (22) was not true, then sup{|v(t)|:t(0,T]}<, contradicting (21). Since h is positive, (1) indicates that the function t k v (t) is increasing on (0,T] and hence there exists lim t 0 + t k v (t). The next part of the proof is divided into three cases.

  1. (i)

    Let us assume that

    lim t 0 + t k v (t)=0.
    (23)

Since t k v (t) is increasing on (0,T], we deduce t k v (t)>0 and v (t)>0 for t(0,T]. Consequently, v is increasing on (0,T]. This together with v>0 and vA C loc 1 (0,T] yields sup{v(t):t(0,T]}<. Therefore, h(t,v(t)) L 1 [0,T]. Integrating (1) and using (23), we obtain

t k v ( t ) = 0 t s k h ( s , v ( s ) ) d s , t ( 0 , T ] , 0 < v ( t ) = 1 t k 0 t s k h ( s , v ( s ) ) d s < 0 T h ( s , v ( s ) ) d s < , t ( 0 , T ] ,

which is a contradiction to (22).

  1. (ii)

    Assume that

    lim t 0 + t k v (t)=>0.

Then there exists t 0 (0,T) such that

v (t)>/2 t k for t(0, t 0 ].
(24)

Hence, v is increasing on (0, t 0 ] and there exists lim t 0 + v(t)0.

By integration, we obtain from (24)

v( t 0 )v(t)> ( t 0 1 k t 1 k ) 2 ( 1 k ) ,t(0, t 0 ]

and, by virtue of k>1 and >0, we arrive at

lim t 0 + v(t)v( t 0 ) lim t 0 + ( t 0 1 k t 1 k ) 2 ( 1 k ) =,

which again is a contradiction.

  1. (iii)

    Assume that

    lim t 0 + t k v (t)=<0.
    (25)

Then there exists t 0 (0,T) such that

v (t)</2 t k for t(0, t 0 ].
(26)

Hence, v is decreasing on (0, t 0 ] and there exists lim t 0 + v(t)(0,]. By integration, we obtain from (26)

v( t 0 )v(t)< ( t 0 1 k t 1 k ) 2 ( 1 k ) ,t(0, t 0 ]

and, since k>1 and <0, we have

lim t 0 + v(t)v( t 0 ) lim t 0 + ( t 0 1 k t 1 k ) 2 ( 1 k ) =.

In addition (25) yields lim t 0 + v (t)= and (20) follows. □

Now, we investigate the existence and properties of blow-up solutions of (1), for the case that the function h has the form

h(t,x)=ψ(t)+g(t,x).

In particular, we study the equation

v (t)+ k t v (t)=ψ(t)+g ( t , v ( t ) ) ,
(27)

where k(1,) and ψ, g satisfy the conditions

( H 1 0 ) ψ L 1 [0,T] and ψ>0 for a.e. t[0,T],

( H 2 0 ) gCar([0,T]×[0,)) and g(t,x) is increasing in x for a.e. t[0,T],

( H 3 0 ) 0g(t,x)ϕ(x) for a.e. t[0,T] and all x[0,), where ϕC[0,), ϕ is increasing and

lim x ϕ ( x ) x =0.

Together with (27) we consider the conditions

lim t 0 v(t)=,v(T)=0,
(28)

and

v (T)=c.
(29)

We define a positive solution to problem (27), (28) as a function vA C loc 1 (0,T] such that v>0 on (0,T), v satisfies the boundary conditions (28), and (27) holds for a.e. t[0,T].

Define a set Z by

Z= { v A C loc 1 ( 0 , T ] : v  is a positive solution of (27), (28) } .
(30)

Clearly, for each vZ there exists c0 such that (29) holds.

Lemma 5 Let ( H 1 0 )-( H 3 0 ) hold. Let vZ. Then v <0 on (0,T).

Proof Since ( t k v ) = t k ( v +(k/t) v ), the integration of (27) over the subinterval [t,T](0,T] gives

T k v (T) t k v (t)= t T s k ( ψ ( s ) + g ( s , v ( s ) ) ) ds.

Assume that v ( t 0 )=0 for some t 0 (0,T). Then

T k v (T)= t 0 T s k ( ψ ( s ) + g ( s , v ( s ) ) ) ds,

which is not possible, since v (T)0 and t 0 T s k (ψ(s)+g(s,v(s)))ds>0. □

Moreover, let us define sets

R= { v Z : t k v  is bounded on  ( 0 , T ] }
(31)

and

R c = { v R : v ( T ) = c } for c0.
(32)

It is obvious that R= c 0 R c . Using the substitution

a=k,u(t)= t a v(t)= t k v(t)for t(0,T],
(33)

we can rewrite (27) and obtain the form

u (t)+ a t u (t) a t 2 u(t)= t a ( ψ ( t ) + g ( t , t a x ) ) .
(34)

We now introduce a function f,

f(t,x)= t a ( ψ ( t ) + g ( t , t a x ) ) for a.e. t[0,T] and all x[0,).

Under conditions ( H 1 0 )-( H 3 0 ), the function f satisfies (H1)-(H3); see the proof of [[20], Theorem 5.1]. Therefore, the results of Section 2 hold for problem (34), (4). As in Section 2, we define the sets S and S c for solutions of (34).

The following result is the key-stone to the analysis of the structure of the set and describes the relation between the sets S and .

Theorem 4 Let ( H 1 0 )-( H 3 0 ) hold and c0. Let us assume that (33) holds. Then v R c if and only if u S c 0 and c 0 = T k c.

Proof () Let v R c . Then (28), (32), and (33) provide the following properties: uA C loc 1 (0,T], u (T)= c 0 = T k c, and u is bounded on (0,T]. It follows from the equality

( t k v ( t ) ) = t k ( ψ ( t ) + g ( t , v ( t ) ) ) for a.e. t[0,T]

that

( t a ( t a u ( t ) ) ) = t a ( ψ ( t ) + g ( t , t a u ( t ) ) ) =f ( t , u ( t ) ) for a.e. t[0,T].

Hence,

u (t)+ a t u (t) a t 2 u(t)=f ( t , u ( t ) ) for a.e. t[0,T].

We now argue as in the proof of [[20], Lemma 3.3] and arrive at

u(t)=t c 0 T a + 1 | a + 1 | ( T a 1 t a 1 ) +t t T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) ds
(35)

for t(0,T]. Since u is bounded on (0,T], we have f(t,u(t)) L 1 [0,T]. Then, by [[20], Lemma 2.1],

0 t a 1 t T s a + 1 f ( s , u ( s ) ) d s 0 T f ( s , u ( s ) ) d s , 0 t T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) d s 1 | a + 1 | 0 T f ( s , u ( s ) ) d s ,

for t(0,T], and hence, sup{|u(t)|+| u (t)|:t(0,T]}<. Now, [[23], Corollary 1] guarantees that u can be extended on [0,T] with uA C 1 [0,T] such that the equality (35) holds for t[0,T]. Consequently, u S c 0 .

() Let u S c 0 and c 0 = T k c. Then it follows from (4), (5), and (33) that vA C loc 1 (0,T] and v (T)=c, lim t 0 + t k v(t)=u(0)=0. Hence, t k v is bounded on (0,T] and

( t a ( t a u ( t ) ) ) =f ( t , u ( t ) ) for a.e. t[0,T]

yields ( t k v ( t ) ) =f(t, t k v(t)), that is,

v (t)+ k t v(t)= t k f ( t , t k v ( t ) ) =ψ(t)+g ( t , v ( t ) ) for a.e. t[0,T].

Since u S c 0 , the equality (35) holds on [0,T], and consequently,

t a u(t) t a + 1 t T s a 2 ( s T ξ a + 1 f ( ξ , u ( ξ ) ) d ξ ) ds,t(0,T].

Hence, note that a+1<0, lim t 0 t a u(t)=, and so lim t 0 v(t)=. As a result, we have v R c . □

Now, we are in the position to provide results on the solvability of problem (27), (28) and formulate the properties of its solutions.

Theorem 5 Let ( H 1 0 )-( H 3 0 ) hold. Then the following statements hold:

  1. (a)

    For each c0 the set R c is nonempty and there exist v c , min , v c , max R c such that v c , min (t)v(t) v c , max (t) for t(0,T] and v R c .

  2. (b)

    If c 1 > c 2 0, v i R c i , i=1,2, then v 1 (t)> v 2 (t) for t(0,T).

  3. (c)

    If c0, v i R c , i=1,2, and v 1 ( t 0 )> v 2 ( t 0 ) for some t 0 (0,T), then either v 1 (t)> v 2 (t) for t(0,T) or there exists t ( t 0 ,T) such that v 1 (t)> v 2 (t) for t(0, t ) and v 1 (t)= v 2 (t) for t[ t ,T].

  4. (d)

    R c is a singleton set for each c[0,)Γ, where Γ[0,) is at most countable.

Proof The result follows by combining results from Lemma 1(a), (b), (c), and (e) with those from Theorem 4. □

If c0, c 0 = T k c, and if (33) holds, then Theorems 4 and 5 yield

v c , min (t)= t k u c 0 , min (t), v c , max (t)= t k u c 0 , max (t),t(0,T].
(36)

The next theorem shows that the set

V= { ( t , x ) R 2 : t ( 0 , T ] , x v 0 , min ( t ) }

is covered by graphs of the functions from .

Theorem 6 Let ( H 1 0 )-( H 3 0 ) hold. Then, for each t 0 (0,T) and each Q v 0 , min ( t 0 ), there exists vR satisfying v( t 0 )=Q. In particular, if for some c0 and t ˜ (0,T) the inequality v c , min ( t ˜ )< v c , max ( t ˜ ) holds, then for each Q ˜ ( v c , min ( t ˜ ), v c , max ( t ˜ )) there exist v R c satisfying v( t ˜ )= Q ˜ .

Proof Choose t 0 (0,T) and Q v 0 , min ( t 0 ). Define A:= t 0 k Q. Then, using (33) and (36), we deduce u 0 , min ( t 0 )= t 0 k v 0 , min ( t 0 )A. Therefore, by Corollary 1, there exists uS such that u( t 0 )=A. Consequently, (33) and Theorem 4 give vR and v( t 0 )=Q. The last statement follows from Theorem 5(c). □

By Corollary 4, there exists α 0 >0 such that

α 0 = u 0 , min (0)= lim t 0 + u 0 , min ( t ) t ,
(37)

and hence, by (33) and (36) with c= c 0 =0,

α 0 = lim t 0 + t k 1 v 0 , min (t).
(38)

Using constants from the interval [ α 0 ,), we can uniquely determine all functions in .

Theorem 7 Let ( H 1 0 )-( H 3 0 ) hold and let α 0 be as in (38). Then the following statements hold:

  1. (a)

    For each vR there exists α[ α 0 ,) such that

    lim t 0 + t k 1 v(t)=α, lim t 0 + t k v (t)=(1k)α.
    (39)
  2. (b)

    For each α[ α 0 ,) there exists a unique vR satisfying (39).

Proof (a) Choose vR. Using (33), (37), and Theorem 4, we have t k v=uS and

lim t 0 + t k 1 v(t)= lim t 0 u ( t ) t = u (0).
(40)

By Lemma 1(c), u (0) α 0 . Hence, if we denote u (0)=α, the first condition in (39) follows. Moreover,

lim t 0 + t k v (t)= lim t 0 + t k ( t k u ( t ) ) = lim t 0 + ( k u ( t ) t + u ( t ) ) =(1k) u (0),
(41)

and therefore the second condition in (39) holds.

  1. (b)

    Choose α[ α 0 ,). Theorem 2 guarantees that there exists a unique uS satisfying u (0)=α. Using (33) and Theorem 4, we conclude t k u=vR. Then v satisfies (40) and (41) which results in (39). It remains to prove that v is unique. Assume that there exists a function wR, wv, such that w satisfies (39). Let z= t k w. Then zS and z (0)= lim t 0 z(t)/t= lim t 0 t k 1 w(t)=α. Consequently, by Theorem 2, we arrive at z=u and w=v, which is a contradiction. □

Let α 0 be given in (38) and choose α[ α 0 ,). Keeping Theorem 7(b) in mind, there exists a unique vR satisfying (39). We denote such v by v α and define a function w α :[0,T]R as

w α (t):={ t k 1 v α ( t ) if  t ( 0 , T ] , α if  t = 0 .
(42)

Then w α C[0,T]. We now specify further properties of functions w α .

Theorem 8 Let ( H 1 0 )-( H 3 0 ) hold and let α 0 be from (38). Let α,β[ α 0 ,) and α<β. Then either w α (t)< w β (t) for t[0,T) or there exists t (0,T) such that w α (t)< w β (t) for t[0, t ) and w α (t)= w β (t) for t[ t ,T].

Proof According to (42), there exist c 1 , c 2 [0,), v α R c 1 , and v β R c 2 such that w α (t)= t k 1 v α (t), w β (t)= t k 1 v β (t) for t(0,T]. Also, w α , w β C[0,T] and, since w α (0)=α<β= w β (0), it follows that w α < w β in a right neighborhood of t=0. Therefore, there exists t 0 (0,T) such that v α ( t 0 )< v β ( t 0 ) and hence, c 1 c 2 . (Note that if c 1 > c 2 , then Theorem 5(b) yields v α (t)> v β (t) for t(0,T), which is not possible.) The result now follows from Theorem 5(b) for c 1 < c 2 and from Theorem 5(c) for c 1 = c 2 . □

Corollary 6 Let { α n }[ α 0 ,) be a convergent sequence and let us denote lim n α n =α. Then lim n w α n = w α in C[0,T].

Proof The proof is indirect. Assume that the statement of the corollary does not hold. Then there exist ε, { t n }[0,T] and a subsequence of { α n }, we denote it again by { α n }, such that

| w α n ( t n ) w α ( t n )|εfor nN.
(43)

Let N 1 ={nN: α n α} and N 2 ={nN: α n α}. Let N 1 be countable. Then there exists a decreasing subsequence { α n } n N of { α n } n N 1 . By Theorem 8, the sequence { w α n } is not increasing on [0,T] and w α n w α on this interval. Using (33) and (42), we obtain w α n (t)= u α n (t)/t, w α (t)= u α (t)/t. Since lim n u α n = u α in C 1 [0,T] it follows from Theorem 3, note that w α n (0)= α n , that w α (0)=α and that lim n w α n (t)= w α (t) for t[0,T]. This fact together with the monotonicity of { w α n } and w α C[0,T] gives lim n w α n = w α in C[0,T], which contradicts (43). Hence lim n w α n = w α in C[0,T]. □

For α > α 0 , we now define a set W ,

W = { w α  from  ( 42 )  with  v α R  satisfying  ( 39 ) : α 0 α α } .
(44)

Theorem 9 Let ( H 1 0 )-( H 3 0 ) hold and let α 0 be from (38). For each α > α 0 the set W defined by (44) is compact in C[0,T].

Proof Consider α[ α 0 , α ]. By Theorem 8 we have

w α 0 (t) w α (t) w α (t),t[0,T].

According to (33), (42), and Theorem 4 we have

u α (t)=t w α (t), u α (0)= w α (0)=α, u α S,
(45)
u α 0 (t) u α (t) u α (t),t[0,T].
(46)

Let us denote

c α := u α (T)and S := α 0 α α S c α .

Conditions (4) and (46) imply 0 c α 0 c α c α . Consequently, using Theorem 4, we obtain

w α W if and only if u α S for α [ α 0 , α ] .
(47)

Consider a sequence { w α n } W . Then there exists a sequence { u α n } S with u α n =t w α n . By Lemma 1(f), the set S is compact in C 1 [0,T]. Therefore, there exists a subsequence { u α k } which converges in C 1 [0,T] to u ¯ S . In particular, lim k u α k (0)= u ¯ (0). Let us denote α ¯ := u ¯ (0), then we have lim k α k = α ¯ [ α 0 , α ]. This yields, by Corollary 6, lim k w α k = w α ¯ in C[0,T] and w α ¯ W . Thus we showed that for any sequence in W , there exists a uniformly converging subsequence on [0,T] with a limit in W . □

4 Examples

Using the open domain Matlab Code bvpsuite, we numerically simulate three model problems in order to illustrate theoretical statements made above.

4.1 Matlab Code bvpsuite

The MATLABTM software package bvpsuite is designed to solve BVPs in ODEs and differential algebraic equations. The solver routine is based on a class of collocation methods whose orders may vary from 2 to 8. Collocation has been investigated in context of singular differential equations of first and second order in [24, 25], respectively. This method could be shown to be robust with respect to singularities in time and retains its high convergence order in case that the analytical solution is appropriately smooth. The code also provides an asymptotically correct estimate for the global error of the numerical approximation. To enhance the efficiency of the method, a mesh adaptation strategy is implemented, which attempts to choose grids related to the solution behavior, in such a way that the tolerance is satisfied with the least possible effort. The error estimate procedure and the mesh adaptation work dependably provided that the solution of the problem and its global error are appropriately smooth.b The code and the manual can be downloaded from http://www.math.tuwien.ac.at/~ewa. For further information see [26]. This software proved useful for the approximation of numerous singular BVPs important for applications; see e.g. [9, 2729].

Example 1 The first example is used to comment on Theorem 5 and to show that if a function g(t,x) in (27) is increasing in x, the graphs of positive solutions of problem (27), (28) cannot intersect. Here, we choose

k=2,T=1,ψ(t)=0.01,g(t,x)= x 2 / 3 ,

and thus problem (27), (28) becomes

v (t)+ 2 t v =0.01+ v 2 3 (t),v(0)=,v(1)=0.
(48)

Using substitution (33), which has the form

a=2,u(t)= t 2 v(t),
(49)

we transform (48) onto the form, cf. (34),

u (t) 2 t u (t)+ 2 t 2 u(t)= t 2 ( 0.01 + ( t 2 u ( t ) ) 2 3 ) ,

and we solve the terminal value problem given by

u (t) 2 t u (t)+ 2 t 2 u(t)=0.01 t 2 + t 2 3 u 2 3 (t),u(1)=0, u (1)=c,
(50)

where c=0,1,2,3, and 4. Then, using (49), we recalculate these solution to their blow-up form and obtain different positive solutions v of problem (48). Note that their graphs are ordered after c, cf. Figure 1.

Figure 1
figure 1

Example 1. Solutions u(t) of (50) (first) and v(t) of (48) (second) for different values of c, c=0 orange, c=1 green, c=2 red, c=3 dark blue, c=4 magenta. We used a collocation method of order 8 based on Gaussian points and grid adaptation strategy to satisfy the tolerances TO L a b s =TO L r e l = 10 7 .

Example 2 We designed this example to see if the condition ( H 2 0 ) is essential for statement (b) of Theorem 5. To this aim, we choose a function g(t,x) which is decreasing in x in (27). Let

k=2,T=1,ψ(t)=0.01,g(t,x)= 1 1 + x .

Then problem (27), (28) has the form

v (t)+ 2 t v =0.01+ 1 1 + v ( t ) ,v(0)=,v(1)=0.
(51)

Using substitution (49), we rewrite (51) and have

u (t) 2 t u (t)+ 2 t 2 u(t)= t 2 ( 0.01 + 1 1 + t 2 u ( t ) ) .

The resulting terminal value problem reads

u (t) 2 t u (t)+ 2 t 2 u(t)=0.01 t 2 + t 4 t 2 + u ( t ) ,u(1)=0, u (1)=c,
(52)

where c=0,1,2,3, and 4. Then, using (49), we obtain different positive blow-up solutions of problem (51). However, their graphs are again ordered; see Figure 2. Hence, condition ( H 2 0 ) seems not to be necessary for statement (b) of Theorem 5 to hold. On the other hand, the question, if another model with a decreasing nonlinearity g exists, where solutions cross each other and are not ordered after c, remains open.

Figure 2
figure 2

Example 2. Solutions u(t) of (52) (first) and v(t) of (51) (second) for different values of c, c=0 orange, c=1 green, c=2 red, c=3 dark blue, c=4 magenta. We used a collocation method of order 8 based on Gaussian points and grid adaptation strategy to satisfy the tolerances TO L a b s =TO L r e l = 10 8 .

Example 3 Since, Example 2 did not provide the solution crossings, we constructed another model problem, where in (27) the function g(t,x) oscillates in x. Let

k=2,T=1,ψ(t)=0.01,g(t,x)=10 ( 1 + sin ( 10 t 2 x ) ) ;

then problem (27), (28) has the form

v (t)+ 2 t v =10.01+10sin ( 10 t 2 v ( t ) ) ,v(0)=,v(1)=0.
(53)

Using (49), we again transform (53) and solve the terminal value problem

u (t) 2 t u (t)+ 2 t 2 u(t)= t 2 ( 10.01 + 10 sin ( 10 u ( t ) ) ) ,u(1)=0, u (1)=c,
(54)

where c=0,1,2,3, and 4. Then, using (49), we calculate the related positive blow-up solutions of problem (53). As can be seen in Figure 3, their graphs are again ordered after c and the question of solution crossings remains open even if we consider oscillating nonlinearities.

Figure 3
figure 3

Example 3. Solutions u(t) of (54) (first) and v(t) of (53) (second) for different values of c, c=0 orange, c=1 green, c=2 red, c=3 dark blue, c=4 magenta. We used a collocation method of order 8 based on Gaussian points and grid adaptation strategy to satisfy the tolerances TO L a b s =TO L r e l = 10 9 .

To show the advantage of the grid adaptation, we carried out the following additional test. For c=4, we evaluated the solution u at the point t=0.04 and solved the boundary value problem

v (t)+ 2 t v =10.01+10sin ( 10 t 2 v ( t ) ) ,v(0.04):= u ( 0.04 ) 0.04 2 ,v(1)=0.
(55)

In Figure 4, solutions u(t) of (54) and v(t) of (55), together with the final grids are shown on the interval [0,1] (left) and [0.04,1] (right). Due to the different smoothness of u and v, their accuracy strongly varies. While the absolute and relative errors of u (or rather their estimates) on the grid with 52 grid points is 10−13, the absolute and relative errors of v on the grid with exactly the same number of points is only 10−9. Note that the code has automatically adapted the location of the grid points to correctly reflect the solution behavior.

Figure 4
figure 4

Example 3, c=4. Solutions u(t) of (54) and v(t) of (55), shown on the interval [0,1] (first) and [0.04,1] (second). We used a collocation method of order 8 based on Gaussian points and grid adaptation strategy.

5 Conclusions

We have described the set of all positive solutions of problem (1), (2), where h has the form h(t,x)=ψ(t)+g(t,x) and assumptions ( H 1 0 )-( H 3 0 ) hold. By Theorem 5 we know that is nonempty and that for each c0 there exists at least one function vR fulfilling v (T)=c. In addition, graphs of functions from do not intersect and has a minimal element v 0 , min .

If we choose an arbitrary c >0 and denote R ={vR: c v (T)0}, then there exists a maximal element v c , max in R . Clearly v 0 , min (T)=0 and v c , max (T)= c . According to Theorem 6, the interior of the set {(t,x) R 2 :0tT, v 0 , min (t)x v c , max (t)} is fully covered by graphs of functions from . Finally, we deduce from Theorems 7-9 that the set { t k 1 v:v R } is compact in C[0,T].

Example 1 illustrates the results of Section 3 and hence, according to ( H 2 0 ), we have chosen the increasing function g(x)= x 2 / 3 in (48). Figure 1 shows ordered graphs of solutions.

In contrast to this, in Example 2 and Example 3, we have chosen the decreasing function g(x)=1/(1+x) in (51) and the non-monotonous function g(t,x)=10(1+sin(10 t 2 x)) in (53), respectively. We can see on Figure 2 and Figure 3 that the graphs of solutions are ordered in both cases. But to prove such order of solutions for non-increasing g is an open problem. On the other hand, the question of the construction of problem (27), (28) whose positive solutions cross each other remains open, as well.

Endnotes

aNote that u 0 , min (T)= u 0 , max (T)=v(T)=0.

bThe required smoothness of higher derivatives is related to the order of the used collocation method.

References

  1. Bongiorno V, Scriven LE, Davis HT: Molecular theory of fluid interfaces. J. Colloid Interface Sci. 1967, 57: 462-475.

    Article  Google Scholar 

  2. Gouin H, Rotoli H: An analytical approximation of density profile and surface tension of microscopic bubbles for van der Waals fluids. Mech. Res. Commun. 1997, 24: 255-260. 10.1016/S0093-6413(97)00022-0

    Article  Google Scholar 

  3. van der Waals JD, Kohnstamm R: Lehrbuch der Thermodynamik, Vol. 1. Maas & von Suchtelen, Leipzig; 1908.

    Google Scholar 

  4. Fife PC Lecture Notes in Biomathematics 28. In Mathematical Aspects of Reacting and Diffusing Systems. Springer, New York; 1979.

    Chapter  Google Scholar 

  5. Fischer RA: The wave of advance of advantageous genes. Ann. Eugenics 1937, 7: 355-369. 10.1111/j.1469-1809.1937.tb02153.x

    Article  Google Scholar 

  6. Abraham FF: Homogeneous Nucleation Theory. Academic Press, New York; 1974.

    Google Scholar 

  7. Linde AP: Particle Physics and Inflationary Cosmology. Harwood Academic, Chur; 1990.

    Book  Google Scholar 

  8. Derrick GH: Comments on nonlinear wave equations as models for elementary particles. J. Math. Phys. 1965, 5: 1252-1254.

    Article  MathSciNet  Google Scholar 

  9. Kitzhofer G, Koch O, Lima P, Weinmüller EB: Efficient numerical solution of the density profile equation in hydrodynamics. J. Sci. Comput. 2007, 32: 411-424. 10.1007/s10915-007-9141-0

    Article  MathSciNet  Google Scholar 

  10. Lima PM, Chemetov NV, Konyukhova NB, Sukov AI: Analytical-numerical investigation of bubble-type solutions of nonlinear singular problems. J. Comput. Appl. Math. 2006, 189: 260-273. 10.1016/j.cam.2005.05.004

    Article  MathSciNet  Google Scholar 

  11. Rachůnková I, Koch O, Pulverer G, Weinmüller EB: On a singular boundary value problem arising in the theory of shallow membrane caps. J. Math. Anal. Appl. 2007, 332: 523-541. 10.1016/j.jmaa.2006.10.006

    Article  MathSciNet  Google Scholar 

  12. Rachůnková I, Staněk S, Weinmüller EB, Zenz M: Neumann problems with time singularities. Comput. Math. Appl. 2010, 60: 722-733. 10.1016/j.camwa.2010.05.019

    Article  MathSciNet  Google Scholar 

  13. Rachůnková I, Tomeček J: Bubble-type solutions of nonlinear singular problems. Math. Comput. Model. 2010, 51: 658-669. 10.1016/j.mcm.2009.10.042

    Article  Google Scholar 

  14. Astashova I: On power and non-power asymptotic behavior of positive solutions to Emden-Fowler type higher-order equations. Adv. Differ. Equ. 2013, 2013(220):1-15.

    Google Scholar 

  15. Karls M, Mohammed A: Integrability of blow-up solutions to some non-linear differential equations. Electron. J. Differ. Equ. 2004, 33: 1-8.

    MathSciNet  Google Scholar 

  16. Kiguradze IT, Chanturia TA: Asymptotic Properties of Solutions of Nonautonomous Ordinary Differential Equations. Kluwer Academic, London; 1993.

    Book  Google Scholar 

  17. Jaroš J, Kusano T: On black hole solutions of second order differential equations with a singular nonlinearity in the differential operator. Funkc. Ekvacioj 2000, 43: 491-509.

    Google Scholar 

  18. Jaroš J, Kusano T: On white hole solutions of a class of nonlinear ordinary differential equations of the second order. Funkc. Ekvacioj 2002, 45: 319-339.

    Google Scholar 

  19. Tanigawa T: On the structure of positive solutions of a class of fourth order nonlinear differential equations. Ann. Mat. Pura Appl. 2006, 185: 521-536. 10.1007/s10231-005-0165-5

    Article  MathSciNet  Google Scholar 

  20. Rachůnková I, Staněk S: Properties of the set of positive solutions to Dirichlet boundary value problems with time singularities. Cent. Eur. J. Math. 2013, 11: 112-132. 10.2478/s11533-012-0047-1

    Article  MathSciNet  Google Scholar 

  21. Kiguradze IT, Shekhter BL: Singular boundary value problems for second-order ordinary differential equations. J. Sov. Math. 1988, 43: 2340-2417.

    Article  Google Scholar 

  22. Lomtatidze A, Malaguti L: On a two-point boundary value problem for the second order differential equations with singularities. Nonlinear Anal. 2003, 52: 1553-1567. 10.1016/S0362-546X(01)00148-1

    Article  MathSciNet  Google Scholar 

  23. Rachůnková I, Spielauer A, Staněk S, Weinmüller EB: Positive solutions of nonlinear Dirichlet BVPs in ODEs with time and space singularities. Bound. Value Probl. 2013, 2013(6):1-22.

    Google Scholar 

  24. de Hoog F, Weiss R: Collocation methods for singular boundary value problems. SIAM J. Numer. Anal. 1978, 15: 198-217. 10.1137/0715013

    Article  MathSciNet  Google Scholar 

  25. Weinmüller EB: Collocation for singular boundary value problems of second order. SIAM J. Numer. Anal. 1986, 23: 1062-1095. 10.1137/0723074

    Article  MathSciNet  Google Scholar 

  26. Kitzhofer G, Pulverer G, Simon C, Koch O, Weinmüller EB: The new MATLAB solver bvpsuite for the solution of singular implicit BVPs. J. Numer. Anal. Ind. Appl. Math. 2010, 5: 113-134.

    Google Scholar 

  27. Budd C, Koch O, Weinmüller EB: From nonlinear PDEs to singular ODEs. Appl. Numer. Math. 2006, 56: 413-422. 10.1016/j.apnum.2005.04.012

    Article  MathSciNet  Google Scholar 

  28. Hammerling R, Koch O, Simon C, Weinmüller EB: Numerical solution of eigenvalue problems in electronic structure computations. J. Comp. Physiol. 2010, 181: 1557-1561.

    Google Scholar 

  29. Rachůnková I, Pulverer G, Weinmüller EB: A unified approach to singular problems arising in the membrane theory. Appl. Math. 2010, 55: 47-75. 10.1007/s10492-010-0002-z

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

This research was supported by the grants IGA PrF_2013_013 and IGA-PrF_2014028. The authors are grateful to the referees for useful comments and suggestions which improved the paper.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Irena Rachůnková.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The authors read and approved the final draft. IR and SS contributed to the analytical part of the paper and JG and EW contributed to the numerical part of the paper.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gschwindl, J., Rachůnková, I., Staněk, S. et al. Positive blow-up solutions of nonlinear models from real world dynamics. Bound Value Probl 2014, 121 (2014). https://doi.org/10.1186/1687-2770-2014-121

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-121

Keywords