ABSTRACT
A nonlinear boundary value problem related to an equation of Kirchhoff type is considered. The existence of multiple positive solutions is proved through Avery-Peterson Fixed Point Theorem. A numerical method based on Levenberg-Marquadt algorithm combined with a heuristic process is present in order to align numerical and theoretical aspects.
Keywords:
Multiple solution; Kirchhoff Equation; numerical solutions
RESUMO
Um problema de valor de fronteira associado é uma equação do tipo Kirchhoff é considerado. Resultados de existência de múltiplas soluções são provados utilizando o teorema de Avery-Peterson. Um método numérico baseado no algoritmo de Levenberg-Maquardt combinado com um processo heurístico é apresentado com finalidade de elucidar aspectos numéricos e teóricos.
Palavras-chave:
Múltiplas soluções; equação de Kirchhoff; soluções numéricas
1 INTRODUCTION
In this paper we present a study on second order equation Kirchhoff problem, given by
where and are continuous maps.
Variations of (1.1) can be related to stationary state of Kirchhoff equation 88 G.R. Kirchhoff. “Vorlesungen u¨ber mathematische physik: mechanik”, volume 1. Teubner (1876). . For example, as elucidated in 1111 T.F. Ma, E.S. Miranda & M.B. de Souza Cortes. A nonlinear differential equation involving reflection of the argument. Arch. Math.(Brno), 40(1) (2004), 63-68. considering argument reflection and later in 11 C. Alves, F. Corrêa & T.F. Ma. Positive solutions for a quasilinear elliptic equation of Kirchhoff type. Computers & Mathematics with Applications, 49(1) (2005), 85-93. with more general boundary conditions a stationary state of the Kirchhoff equation of kind of:
can be associated with their respectives equations. This kind of equation commonly appear in the context of free vibrations in elastic strings, consequently, they form a relevant object of study.
Many authors have studied problems related to (1.2). For a wide spectrum of study, we recommend 22 P. Amster & M.C. Mariani. A fixed point operator for a nonlinear boundary value problem. Journal of mathematical analysis and applications, 266(1) (2002), 160-168., 33 A. Arosio & S. Panizzi. On the well-posedness of the Kirchhoff string. Transactions of the American Mathematical Society, 348(1) (1996), 305-330.), (66 M. Chipot & J.F. Rodrigues. On a class of nonlocal nonlinear elliptic problems. ESAIM: Mathematical Modelling and Numerical Analysis, 26(3) (1992), 447-467.), (1010 T.F. Ma. Existence results for a model of nonlinear beam on elastic bearings. Applied Mathematics Letters, 13(5) (2000), 11-15. and more recently 55 G. Caristi, S. Heidarkhani & A. Salari. Variational Approaches to Kirchhoff-Type Second-Order Impulsive Differential Equations on the Half-Line. Results in Mathematics, 73(1) (2018), 44.. Concerning the problem (1.1) we can observe by 1212 A.L.M. Martinez, E.V. Castelani, J. Da Silva & W.V.I. Shirabayashi. A note about positive solutions for an equation of Kirchhoff type. Applied Mathematics and Computation, 218(5) (2011), 2082-2090. that studies with theoretical aspects were developed using Banach’s Fixed Point Theorem or Leray-Schauder Alternative combined with Krasnoselskii’s Theorem. Consequently, it is natural to ask about existence results by using Avery-Peterson Theorem 44 R. Avery & A.C. Peterson. Three positive fixed points of nonlinear operators on ordered Banach spaces. Computers & Mathematics with Applications, 42(3-5) (2001), 313-322.. In this sense, this work complement previous results of existence with the cited theorem, and this study is present in Section 2.
Numerical studies related to second order equations, normally, are presented as illustration of the existence Banach’s Fixed Point Theorem applied to equation but no strategy is given to ilustraste more general results (like the existence results provided by Avery-Peterson Theorem). The common sense suggests that optimization methods combined with heuristic processes can provide good results to find multiple numerical solutions. As consequence, a method based on Levenberg-Maquardt 1313 J.J. Moré. The Levenberg-Marquardt algorithm: implementation and theory. In “Numerical analysis”. Springer (1978), pp. 105-116.), (99 M.I.A. Lourakis. A brief description of the Levenberg-Marquardt algorithm implemented by levmar. Foundation of Research and Technology, 4(1) (2005), 1-6. is shown in Section 3 and comparisons with other methods using classical strategies are established. Still, in Section 3, a simple but effective heuristic process is introduced in order to find multiple numerical solutions. Final remarks are given in Section 4.
2 MULTIPLE SOLUTIONS
Let be the Banach space of the continuously differentiable functions in [0, 1]. Let us consider with the norm
We begin this section by observing that the solutions of (1.1) can be written as:
where G is the Green’s function
Considering the result that we will show, we must note some properties of the function G. In fact, we have that
then G satisfies:
Let m be a constant in [0,1/2]. Thus, we can obtain the inequalities
and
Thus, concerning the representation given in (2.1) we define the operator as:
where M and f are assumed to be continuous functions and there is a constant such that for all y in the domain of M.
We claim that T is continuous and completely continuous. Continuity follows immediately from the Lebesgue dominated convergence theorem and the fact that
with . To show complete continuity we will use the Arzela-Ascoli’s theorem. Let be bounded, in other words, there exists with for each . Now if we have
where is determined by the bounded set and functions q, f and M. It is easy to check that . Then imply that T Ω is a bounded equicontinuous family on [0, 1].
Consequently the Arzela-Ascoli theorem implies is completely continuous.
Consider the following hypotheses:
(A1) There are positive constants d, A, B such that:
Next, we define a cone P by
Remark 1. Denoting
we can extract some properties related to F. Since
from (A1), we find that F is concave, and since, we can conclude that
For the purpose of this work we need to introduce the main tools.
Avery-Peterson theorem.
Now, we need to consider the convex sets
and the closed set
Theorem 2. Let P be a cone in a real Banach space X. Let γ and θ nonnegative continuous convex functionals on P, α be a nonnegative continuous concave functional on P, and ψ be a nonnegative continuous functional on P satisfying for , such that for some positive numbers µ and d,
for all . Suppose
is completely continuous and there exist positive numbers a, b, c with a < b, such that
Then T has at least three distinct fixed points in.
In our main result (given by Theorem 2) we will show that the Problem 1.1 has at least three positive solutions.
Theorem 3.Suppose that the hypothesis (A1) is satisfied. Suppose, in addition, that there exist a, such that f satisfies the following conditions:
Then, Problem (1.1) has at least three positive solutions.
Proof. We will apply Avery-Peterson theorem (as stated in 44 R. Avery & A.C. Peterson. Three positive fixed points of nonlinear operators on ordered Banach spaces. Computers & Mathematics with Applications, 42(3-5) (2001), 313-322.). Let us consider T and P as defined before. Furthermore, we need define the following functionals
Therefore, we need to verify that
is completely continuous and there exist positive numbers b and c with , such that
Using (A1) we have if , and obtain:
Therefore T applies in .
Now, we consider
and
Clearly, we have . Let us demonstrate (2.8). Using (A3) we obtain
Let us demonstrate (2.9). Let with . Then
Now, let us show (2.10). Thus, let u ∈ R(γ, ψ;a;d) with ψ(u) = a. From (A1)-(A2) we have,
Applying Avery-Peterson theorem we obtain the result.
Example 2.1. This example shows a problem that has at least three positive solutions, according to Theorem 3. Suppose that
and
Taking and we find
(A1) for
(A2) for
where
(A3) for
where
Thus, from Theorem 3, the problem has at least three positive solutions.
3 NUMERICAL SOLUTIONS
In most studies, numerical solutions are obtained by fixed point methods, according to 1212 A.L.M. Martinez, E.V. Castelani, J. Da Silva & W.V.I. Shirabayashi. A note about positive solutions for an equation of Kirchhoff type. Applied Mathematics and Computation, 218(5) (2011), 2082-2090.. More specifically, an iterative sequence based on operator given by equation (2.1) define the method. The basic idea of the proposed method is to use the Levenberg-Maquardt algorithm 1313 J.J. Moré. The Levenberg-Marquardt algorithm: implementation and theory. In “Numerical analysis”. Springer (1978), pp. 105-116.. This type of algorithm arises frequently in the context of optimization related to the data adjustment problem 77 C.T. Kelley. “Iterative methods for optimization”, volume 18. Siam (1999).. The Levenberg-Marquardt algorithm can be easily used for the solution of non-linear systems. Algorithm 1 briefly describes the ideas to solve the Problem (1.1).
The motivation for the Algorithm 1 is the fact that fixed point methods are tendentious to find solutions in which the operator T is a contraction and consequently, chosen an initial approximation u 0 we have generally two possibilities: the method converges to solution given by Banach’s Theorem or the method diverges. Anyway, if we have multiple solutions (as in Theorem 3), we can try to find these solutions. For this reason, our proposed algorithm, using an appropriated initial approximation, allows to find others solutions. So the development of an heuristic to find a better initial approximation is relevant.
3.1 An heuristic procedure for initial approximations
We know the solutions that we are looking for must be concave or convex and will satisfy the condition . Thus, approaches by parabolas are reasonable ways of approaching the solution. In this sense, the heuristic procedure proposed in this paper consists in generating parables about initial points as follows:
where nr is a random number in . For practical purposes, the proposed procedure is defined by Algorithm 2.
Naturally, this procedure returns multiple answers. So we need to establish a way to compare solutions in order to distinguish them. Note that the magnitude of the solutions may be different. In this sense, we say that the numerical solutions u* and u** are equivalent if
is satisfied.
Additional tests are run by using an algorithm based on the Banach Fixed Point Theorem. Considering some hypothesis about the functions f and M, we can demonstrate, using the Banach Fixed Point Theorem (details in 1212 A.L.M. Martinez, E.V. Castelani, J. Da Silva & W.V.I. Shirabayashi. A note about positive solutions for an equation of Kirchhoff type. Applied Mathematics and Computation, 218(5) (2011), 2082-2090.), the local convergence of algorithms that use the iterative sequence
This procedure is described by Algorithm 3, which is a version of the algorithm proposed in 1212 A.L.M. Martinez, E.V. Castelani, J. Da Silva & W.V.I. Shirabayashi. A note about positive solutions for an equation of Kirchhoff type. Applied Mathematics and Computation, 218(5) (2011), 2082-2090..
For practical purposes, the procedure is formalized by Algorithm 4, which is equivalent to Algorithm 2 and it calculates possible multiple solutions.
3.2 Examples
The examples below show how the Algorithm 2 can be promissor in order to find multiple solutions. We run Algorithm 2 and 4 with N = 10. For Algorithm 3, we consider as stopping criteria:
and for Algorithm 1 we consider .
Example 3.1. Consider Problem (1.1) defined by
and
The solutions are and . In this example we consider d = 10 and n = 20 points for the spaced mesh. The numerical results are compared with the exact solution.
Running Algorithm 2, the procedure converged to the solution in 9 out of 10 times in which the Algorithm 1 was called. After 7 iterations of Algorithm 1, six initializations converged to the solution u and the obtained accuracy is . After 4 iterations of Algorithm 1, three initializations converged to the solution and the obtained accuracy is . The divergence means that the maximum number of iterations, i = 50, was extrapolated. At the point where the execution stopped, the solution was near of u.
On the other hand, running Algorithm 4, the procedure converged to the solution ū (the less norm) in 8 out of 10 times in which the Algorithm 3 was called and reached the required accuracy after k iterations, with . The two divergences mean that Algorithm 3 stopped running since the function M became null and then there was a divide by zero. For accuracy or order < 10−4, the algorithm diverges in most of times, due to divide by zero.
Example 3.2. Consider the Problem (1.1) defined by
and
The solutions are and . In this example, we consider d = 10 and n = 20 points for the spaced mesh. The numerical results are compared with the exacts solutions.
Running Algorithm 2, the procedure converged to the solution in 10 out of 10 times in which the Algorithm 1 was called. After 5 iterations of Algorithm 1, five initializations converged to the solution u and the obtained accuracy is . After 5 iterations of Algorithm 1, three initializations converged to the solution ū and the obtained accuracy is . After 12 iterations of Algorithm 1, two initializations converged to the solution u and the obtained accuracy is .
On the other hand, running Algorithm 4, the procedure converged to the solution null (the less norm) in 10 out of 10 times in which the Algorithm 3 was called and reached the required accuracy after k iterations, with .
We have conducted an additional test, which run Algorithm 2 using the functions defined in Example 1. In this example, we defined n = 20 and d = 55. Using the criterion established in (3.1), we obtain that the 3 solutions are different. These results illustrate the result of existence given by Theorem 3. Figure 1 shows a graphical representation of these solutions.
Running Algorithm 4, the procedure converged to the solution of less norm in 9 out of 10 times in which the Algorithm 3 was called (according to the criteria described in 3.1) and reached the required accuracy after 5 iterations.
As a conclusion, Algorithm 4 have found only the solution of less norm in all the run examples. On the other hand, Algorithm 2 have found, to the same run examples, all the possible solutions. Therefore, Algorithm 2 is more robust and suitable to determine multiple solutions.
4 FINAL REMARKS
Considering the numerical aspects of this work, we present a new algorithm (Algorithm 1) and a new heuristic that allows obtain multiple solutions. The method has robustness to solve second order equation Kirchhoff problem. Of course, the cost of this robustness is a slightly higher cost of computer processing, especially when we compare the method presented in last section with fixed-point methods (using the operator given by (7)).
Fortunately, the level of processing is not absurd even considering finer meshes. Furthermore, the method can be adapted for parallel programming and consequently, new features can be exploited in the computational field.
REFERENCES
-
1C. Alves, F. Corrêa & T.F. Ma. Positive solutions for a quasilinear elliptic equation of Kirchhoff type. Computers & Mathematics with Applications, 49(1) (2005), 85-93.
-
2P. Amster & M.C. Mariani. A fixed point operator for a nonlinear boundary value problem. Journal of mathematical analysis and applications, 266(1) (2002), 160-168.
-
3A. Arosio & S. Panizzi. On the well-posedness of the Kirchhoff string. Transactions of the American Mathematical Society, 348(1) (1996), 305-330.
-
4R. Avery & A.C. Peterson. Three positive fixed points of nonlinear operators on ordered Banach spaces. Computers & Mathematics with Applications, 42(3-5) (2001), 313-322.
-
5G. Caristi, S. Heidarkhani & A. Salari. Variational Approaches to Kirchhoff-Type Second-Order Impulsive Differential Equations on the Half-Line. Results in Mathematics, 73(1) (2018), 44.
-
6M. Chipot & J.F. Rodrigues. On a class of nonlocal nonlinear elliptic problems. ESAIM: Mathematical Modelling and Numerical Analysis, 26(3) (1992), 447-467.
-
7C.T. Kelley. “Iterative methods for optimization”, volume 18. Siam (1999).
-
8G.R. Kirchhoff. “Vorlesungen u¨ber mathematische physik: mechanik”, volume 1. Teubner (1876).
-
9M.I.A. Lourakis. A brief description of the Levenberg-Marquardt algorithm implemented by levmar. Foundation of Research and Technology, 4(1) (2005), 1-6.
-
10T.F. Ma. Existence results for a model of nonlinear beam on elastic bearings. Applied Mathematics Letters, 13(5) (2000), 11-15.
-
11T.F. Ma, E.S. Miranda & M.B. de Souza Cortes. A nonlinear differential equation involving reflection of the argument. Arch. Math.(Brno), 40(1) (2004), 63-68.
-
12A.L.M. Martinez, E.V. Castelani, J. Da Silva & W.V.I. Shirabayashi. A note about positive solutions for an equation of Kirchhoff type. Applied Mathematics and Computation, 218(5) (2011), 2082-2090.
-
13J.J. Moré. The Levenberg-Marquardt algorithm: implementation and theory. In “Numerical analysis”. Springer (1978), pp. 105-116.
Publication Dates
-
Publication in this collection
Sep-Dec 2018
History
-
Received
15 Mar 2018 -
Accepted
25 June 2018