In this paper, we discuss the existence of positive solutions to the singular Dirichlet boundary value problems (BVPs) for ordinary differential equations (ODEs) of the form
where . The nonlinearity may be singular for the space variables and/or . Moreover, since , the differential operator on the left-hand side of the differential equation is singular at . Sufficient conditions for the existence of positive solutions of the above BVPs are formulated and asymptotic properties of solutions are specified. The theory is illustrated by numerical experiments computed using the open domain MATLAB code bvpsuite, based on polynomial collocation.
MSC: 34B18, 34B16, 34A12.
Keywords:singular ordinary differential equation of the second order; time singularities; space singularities; positive solutions; existence of solutions; polynomial collocation
In the present work, we analyze the existence of positive solutions to the singular Dirichlet BVP,
For such functions, we use the notation . Moreover, may become singular when the space variables x and/or y vanish, which means that may become unbounded for and a.e. and all , and/or it may be unbounded for and a.e. and all . Finally, since , Eq. (1a) has a singularity of the first kind at the time variable because
The differential operator on the left-hand side of Eq. (1a) can be equivalently written as and, after the substitution , it takes the form , which arises in numerous important applications. Operators of such type were studied in phase transitions of Van der Waals fluids [1-4], in population genetics, especially in models for the spatial distribution of the genetic composition of a population [5,6], in the homogeneous nucleation theory , in relativistic cosmology for description of particles which can be treated as domains in the universe , and in the nonlinear field theory , in particular, when describing bubbles generated by scalar fields of Higgs type in the Minkowski spaces .
The aim of this paper is to study the case which is fundamentally different from the case . The latter setting was studied in [11,12], where the structure and properties of the set of all positive solutions to (1a) and (1b) were investigated (the cardinality of this set is a continuum).
In the sequel, we introduce the basic notation and state the preliminary results required in the analysis of problem (1a) and (1b). Here, we focus our attention on the case and prove the existence of at least one positive solution of (1a) and (1b). In contrast to [11,12], we consider the more general situation in which f may be also singular at . This means that we have to deal with the following additional difficulties.
Let u be a positive solution of problem (1a) and (1b), where has a singularity at . Then there exists such that , and hence f is unbounded in a neighborhood of the point . Unfortunately, we do not know the exact position of and therefore, it is not possible to construct a universal Lebesgue integrable majorant for all functions , where are positive solutions of a sequence of auxiliary regular problems. Consequently, the Lebesgue dominated convergence theorem is not applicable and we have to use arguments based on the Vitali convergence theorem instead; see Lemma 2. Another tool used in the proofs is a combination of regularization and sequential techniques with the Leray-Schauder nonlinear alternative.
The investigation of singular Dirichlet BVPs has a long history and a lot of methods for their analysis are available. One of the most important ones is the topological degree method providing various fixed point theorems and existence alternative theorems; see, e.g., Lemma 1. For more information on the topological degree method and its application to numerous BVPs, including Dirichlet problems, we refer the reader to the monographs by Mawhin [13-15].
Throughout this paper, we work with the following conditions on the function f in (1a):
we denote the norms in and , respectively. denotes the set of functions whose first derivative is absolutely continuous on , while is the set of functions having absolutely continuous first derivative on each compact subinterval of . We use the symbol to denote the Lebesgue measure of ℳ.
As a first step in the analysis, we investigate auxiliary regular BVPs of the form
To show the solvability of problem (7a) and (7b), we use the following alternative of Leray-Schauder type which follows from [, Theorem 5.1].
The paper is organized as follows. In Section 2, we collect auxiliary results used in the subsequent analysis. Section 3 is devoted to the study of limit properties of solutions to Eq. (7a). In Section 4, we investigate auxiliary regular problems associated with the singular problem (1a) and (1b). We show their solvability and describe properties of their solutions. An existence result for the singular problem (1a) and (1b) is given in Section 5. Finally, in Section 6, we illustrate the theoretical findings by means of numerical experiments.
In this section, auxiliary statements necessary for the subsequent analysis are formulated.
Hence, cf. (10),
we conclude from the relation
The last equality contradicts (13). Consequently, (12) holds and the result follows. □
Proof Since (cf. (10))
3 Limit properties of solutions to Eq. (7a)
Here, we investigate asymptotic properties of solutions of (7a). We also provide a related integral equation this solution satisfies.
We now discuss the linear Euler differential equation
Let H be the function given in Lemma 3. By Lemma 3(ii), H can be extended on with and −H satisfies (16) for a.e. . Therefore, each function which satisfies Eq. (16) a.e. on has the form for , with some . By assumption we know that satisfies (16) a.e. on , and therefore there exist such that , . Since by assumption on , we have . Consequently, the function u can be extended on the interval in the class and (15) holds on . □
4 Auxiliary regular problems
In order to prove the solvability of problem (7a) and (7b), we first have to investigate the problem
Lemma 7Let (H1) hold. Thenuis a solution of problem (17a) and (17b) if and only ifuis a solution of the integral equation
Let u be a solution of Eq. (18) in . Then . Hence, Lemma 3(ii) (with ρ replaced by ) guarantees that and u is a solution of Eq. (17a). Moreover, . Consequently, u is a solution of problem (17a) and (17b) which completes the proof. □
The following results provide bounds for solutions of problem (17a) and (17b).
and therefore, we have
We are now in the position to prove the existence result for problem (7a) and (7b).
Proof Let S be a positive constant in Lemma 8 and let us define
It follows from Lemma 5 and (9) that
we have by the Lebesgue dominated convergence theorem. Hence, by (28), is a continuous operator. We now show that the set is relatively compact in . It follows from and bounded in that there exists such that
Then by Lemma 5 and (9), the inequalities
holds for a.e. and all (cf. (9)). Consequently, the set is equicontinuous on . Hence, the set is relatively compact in by the Arzelà-Ascoli theorem. As a result, is a compact operator and the condition (i) follows.
Due to the fact that by Lemma 7 any fixed point u of the operator is a solution of problem (17a) and (17b), Lemma 8 guarantees that u satisfies inequality (20). Therefore, has property (ii). Consequently, by Lemmas 1 and 8, for each , problem (7a) and (7b) has a solution u satisfying estimates (19)-(21). □
By Lemma 9, the inequalities
5 The existence result for BVP (1a) and (1b)
This section is devoted to the main result on the existence of positive solutions to the original BVP (1a) and (1b).
Theorem 1Let (H1)-(H3) hold. Then problem (1a) and (1b) has at least one positive solution.
Proof By Lemma 9, for each , problem (7a) and (7b) has a solution satisfying inequalities (29)-(31), where S is a positive constant and . Moreover, by Lemma 10, the sequence is uniformly integrable on . We now prove that is equicontinuous on . Since is a fixed point of the operator given in (27), the equality
follows. By Lemma 4 (for ), the sequence is equicontinuous on . Since the sequence is uniformly integrable on , the sequence is equicontinuous on . Hence, it follows from (33), that is equicontinuous on . We summarize: is bounded in and is equicontinuous on . Also, . Using appropriate subsequences, if necessary, we can assume, by the Arzelà-Ascoli theorem and the Bolzano-Weierstrass theorem, that is convergent in and is convergent in ℝ. Let and . With in (29)-(31), we conclude
it follows from Lemma 2 that
6 Numerical simulations
where is a parameter. We can use the above formulation because problem (1a) and (1b) is solvable for f satisfying the assumptions of Theorem 1 and, therefore, solutions of problem (1a) and (1b) can be computed as solutions of problem (36a) and (36b) using the proper value depending on f. The values are provided for given f in Examples 1 and 2, below.
The reason for changing the boundary conditions from (1b) to (36b) is that the differential equation (36a) subject to (1b) is not well posed; see . However, to enable successful numerical treatment, well-posedness of the model is crucial. This property means that Eq. (36a) subject to proper boundary conditions has at least a locally unique solution,a and this solution depends continuously on the problem data. The well-posedness of the problem is important for two reasons. First of all, it allows to express errors in the solution of the analytical problem in terms of modeling errors and data errors (all measured via appropriate norms). Therefore, when the errors in the data become smaller due to more precise modeling or smaller measurement inaccuracies, the errors in the solution will decrease. The second reason is that the well-posedness decides if the numerical simulation will be at all successful. If the analytical problem is ill-posed, then the inevitable round-off errors can become extremely magnified and fully spoil the accuracy of the approximation.
By , problem (37), (36b) is well posed and therefore it is suitable for the numerical treatment. To see this, we need to look at a general solution of the homogeneous equation
6.1 MATLAB Code bvpsuite
To illustrate the analytical results discussed in the previous section, we solved numerically examples of the form (36a) and (36b) using a MATLAB™ software package bvpsuite 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 two to eight. Collocation has been investigated in the 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 the 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. 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 webcite. For further information, see . This software proved useful for the approximation of numerous singular BVPs important for applications; see, e.g., [3,9,27,28].
Before dealing with two nonlinear models specified in Sections 6.3 and 6.4, we have to compute numerical solutions for a simpler linearc model of the form
where a was chosen as . Since in this case the exact solution is given, , the value is available, , respectively. In Figure 1, the numerical solutions of BVPs (39a) and (39b) are shown. They will be used as starting values for the numerical solution of Examples 1 and 2; see Sections 6.3 and 6.4, respectively. All numerical results have been obtained using collocation with five Gaussian collocation points on an equidistant grid (justified by a very simple solution structure) with the step size 0.01.
Figure 1. Problem (39a) and (39b): Numerical solutions for different values ofa.
6.3 Example 1
We first investigate the following problem:
The nonlinearity f in (40a) has the form
It follows from Theorem 1 that there exists at least one value of such that the related solution u of problem (40a) and (40b) with is positive on with . Using formula (35), we now determine an interval containing all admissible values of c.
Let u be a solution of problem (1a) and (1b) with f from (41). Then by (35), we obtain
In order to solve the nonlinear problem (40a) and (40b), we first have to solve a series of auxiliary problems for parameter-dependent differential equations
We begin the calculations with and increase its value gradually until we arrive at ; cf. (40a). In each step we use the solution of the previous problem to solve the next one. The aim is to find a good starting value for both the solution u and the value before solving the BVP (40a) and (40b), i.e., find the final value of such that .
After the last step, we have solved problem (40a) and (40b) subject to boundary conditions , . In this case, the value of was not small enough to consider it a reasonable approximation for . Therefore, we use a shooting idea combined with a bisection strategy to find a better value for . The complete numerical results for Example 1 can be found in Table 1 and Figure 2.
Figure 2. Problem (40a) and (40b): Numerical solutions for different values ofa. Values of.
Table 1. Problem (40a) and (40b): Complete data of the numerical simulation for different values ofa
6.4 Example 2
The above approach has been also accordingly applied for Example 2. Here, we consider the problem
The right-hand side f in Eq. (46a) now reads
Theorem 1 guarantees the existence of at least one such that a solution u of problem (46a) and (46b) is positive on and holds. We now again determine an interval containing all such values of c. Let u be a solution of problem (1a) and (1b) with f given in (47). Inequality (19) yields
and hence by (35),
For Example 2, the auxiliary ODE is constructed using ODE (39a),
After the last step, we have solved BVP (46a) and (46b) subject to boundary conditions , , but also, in this case, the value of is too large and we have to find a better value for . The complete numerical results for Example 2 can be found in Table 2 and Figure 3.
In the present article, we deal with the existence of positive solutions to the singular Dirichlet problem of the form
where , and the nonlinearity may be singular at the space variables and/or . The main result for the existence of positive solutions of the above BVP is Theorem 1. It is illustrated by numerical simulations using the MATLAB code bvpsuite, based on polynomial collocation. For the successful numerical treatment, the above problem has to be reformulated to obtain its well-posed form
Here, it is only known that , where can be specified depending on functions f arising in Examples 1 and 2. Now, a simple shooting method combined with the bisection idea is used to find c in such a way that .
The authors declare that they have no competing interests.
IR and SS contributed to the analytical part of the work and AS and EBW contributed to its numerical part. All authors read and approved the final version of the manuscript.
Dedicated to Jean Mawhin on the occasion of his 70th birthday.
This research was supported by the grant Matematické modely a struktury, PrF-2012-017. The authors thank the referees for suggestions which improved the paper.
This BVP can have more than one solution, but they may not lay close together.
The required smoothness of higher derivatives is related to the order of the used collocation method.
The nonlinear term in f has been omitted; see (40a) and (46a).
Gouin, H, Rotoli, G: An analytical approximation of density profile and surface tension of microscopic bubbles for Van der Waals fluids. Mech. Res. Commun.. 24, 255–260 (1997). Publisher Full Text
Kitzhofer, G, Koch, O, Lima, P, Weinmüller, E: Efficient numerical solution of the density profile equation in hydrodynamics. J. Sci. Comput.. 32, 411–424 (2007). Publisher Full Text
Fischer, RA: The wave of advance of advantageous genes. Ann. Eugen.. 7, 355–369 (1937). Publisher Full Text
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.. 11(1), 112–132 (2013). Publisher Full Text
Agarwal, RP, O’Regan, D, Rachůnková, I, Staněk, S: Two-point higher-order BVPs with singularities in phase variables. Comput. Math. Appl.. 46, 1799–1826 (2003). Publisher Full Text
Rachůnková, I, Staněk, S: A singular boundary value problem for odd-order differential equations. J. Math. Anal. Appl.. 291, 741–756 (2004). Publisher Full Text
Weinmüller, EB: On the boundary value problems of ordinary second order differential equations with a singularity of the first kind. SIAM J. Math. Anal.. 15, 287–307 (1984). Publisher Full Text
de Hoog, F, Weiss, R: Collocation methods for singular boundary value problems. SIAM J. Numer. Anal.. 15, 198–217 (1978). Publisher Full Text
Weinmüller, EB: Collocation for singular boundary value problems of second order. SIAM J. Numer. Anal.. 23, 1062–1095 (1986). Publisher Full Text
Budd, C, Koch, O, Weinmüller, E: From nonlinear PDEs to singular ODEs. Appl. Numer. Math.. 56, 413–422 (2006). Publisher Full Text
Rachůnková, I, Pulverer, G, Weinmüller, E: A unified approach to singular problems arising in the membrane theory. Appl. Math.. 55, 47–75 (2010). Publisher Full Text