A numerical solution to nonlinear second order three-point boundary value problems in the reproducing kernel space
Introduction
Boundary value problems for nonlinear differential equations arise in the mathematical modeling of viscoelastic and inelastic flows, deformation of beams and plate deflection theory [1], [2], [3]. Multi-point nonlinear boundary value problems, which take into account the boundary data at intermediate points of the interval under consideration, have been addressed by many authors [4], [5], [6], [7], [8]. Despite of the large amount of works which are done on the theoretical aspects of these kind of equations (please see [4], [5], [6], [7], [9], [10], [11] and references therein), few works are available on the numerical analysis. Authors of [12] solved numerically multi-point boundary value problems by using the Sinc-collocation method. In [13], [14], [15], [16], [17], [18], [19], [20], finite difference methods have been proposed for the numerical solution of various nonlocal boundary value problem. Adomian decomposition method [21] and He’s variational iteration method [22] are employed for solving multi-point boundary value problems.
In [7], the existence and uniqueness of solutions of the following second-order three-point boundary value problems have been studied by the monotone iterative methodwhere .
In [8], the authors investigated the numerical solutions of singular second order three-point boundary value problems as follows:where or and . Besides, the authors also imposed many restrictive conditions on and .
In this work, a numerical method will be given for Eqs. (1.1), (1.2). It is worth mentioning that we only need to assume that for Eq. (1.2). Without loss of generality, we may put . The proposed numerical method depends on a new reproducing kernel space. The construction of reproducing kernel space is innovative. Its reproducing kernel function has a lot of good properties that are beneficial to calculation, especially, the new reproducing kernel space satisfies all boundary conditions.
The rest of the paper is organized as follows. A new reproducing kernel space for solving problem Eq. (1.1) is constructed in Section 2; In Section 3, based on the reproducing kernel function and the minimum-point of function, the numerical algorithm is discussed. In Section 4, some numerical examples are shown to verify the effect of our method. Also a conclusion is given in Section 5.
Section snippets
A constructive method for the reproducing kernel space
In order to solve Eq. (1.1), a reproducing kernel space is defined byThe inner product and norm of are defined by Lemma 2.1 The function space is a reproducing kernel space.
The proof can be found in [23]. It is very important to obtain the representation of reproducing kernel function, since it is the base of our algorithm. By Lemma 2.1, there exists a reproducing
Numerical algorithm
In this section, we shall explain how to obtain approximate solution of Eq. (1.1).
Let and . Since operator is bounded, we can rewrite Eq.(1.1) asChoosing a countable dense subset and defining bythen is a complete system of the space (see [23]). Furthermore, we obtain an orthogonal basis of by Gram–Schmidt process, such that Theorem 3.1 If is
Numerical examples
Here some numerical examples are studied to demonstrate the accuracy of the present method. Results obtained by the method are compared with analytical solution of each example and are found to be in good agreement with each other. Example 1 Consider the following nonlinear second-order differential equation:with an exact solution . Applying our method and taking the number of nodes , the relative errors are shown
Conclusion
In this work, we construct a novel reproducing kernel space and give the way to express the reproducing kernel function. A numerical algorithm is presented based on the properties of reproducing kernel function and the minimum-point of function. Through numerical experiments, we demonstrated the efficiency and superiority of our proposed algorithm.
Acknowledgments
The authors appreciate the constructive comments and suggestions provided from the kind referees and editor. This work was supported by Youth Foundation of Heilongjiang Province under Grant QC2010036 and Fundamental Research Funds for the Central Universities under Grant No. HIT.NSRIF.2009050.
References (23)
- et al.
Dynamic boundary value problems of the second-order: Bernstein–Nagumo conditions and solvability
Nonlinear Anal.
(2007) - et al.
Positive solutions to a second order three-point boundary value problem
J. Math. Anal. Appl.
(2003) - et al.
The method of lines for solution of the one-dimensional wave equation subject to an integral conservation condition
Comput. Math. Appl.
(2008) Efficient techniques for the second-order parabolic equation subject to nonlocal specifications
Appl. Numer. Math.
(2005)The one-dimensional heat equation subject to a boundary integral specification
Chaos Soliton Fract.
(2007)Numerical techniques for a parabolic equation subject to an overspecified boundary condition
Appl. Math. Comput.
(2002)Soil settlement effects on buildings
Vibrations of Shells and Plates
(1993)- et al.
Nonlinear nonlocal boundary value problems and perturbed Hammerstein integral equations
Proc. Edinburgh. Math. Soc.
(2006) - et al.
Nonlocal boundary value problems of arbitrary order
J. London Math. Soc.
(2009)
Existence and uniqueness of solutions of second-order three-point boundary value problems with upper and lower solutions in the reversed order
Nonlinear Anal.
Cited by (32)
Ten non-polynomial cubic splines for some classes of Fredholm integral equations
2022, Ain Shams Engineering JournalAn iterative kernel based method for fourth order nonlinear equation with nonlinear boundary condition
2018, Communications in Nonlinear Science and Numerical SimulationAn effective computational method for solving linear multi-point boundary value problems
2018, Applied Mathematics and ComputationCitation Excerpt :In the past twenty years, the problem of existence and uniqueness of the solutions for the MBVPs has been extensively investigated; see [5–13] and the references therein. Meanwhile, quite a few methods have been proposed to handle the numerical solutions of the MBVPs, such as the shooting method [4,14], the finite difference method [15], the Adomian decomposition method [16], the variational iteration method [17], the homotopy analysis method [18], the Sinc-collocation method [19], the optimal homotopy asymptotic method [20], the shifted Jacobi spectral method [21], the reproducing kernel method and its modifications [22–31], and the Pade approximant [32]. More recently, an effective algorithm based mainly on differential transform method was proposed to solve the general n order MBVPs [33].
An iterative reproducing kernel method in Hilbert space for the multi-point boundary value problems
2018, Journal of Computational and Applied MathematicsReproducing kernel method for the numerical solution of the Brinkman–Forchheimer momentum equation
2017, Journal of Computational and Applied MathematicsCitation Excerpt :In 1950, Aronszajn produced a systematic reproducing kernel space method based on the Bergman’s works [9,10]. Since then, many researchers in various fields of science and engineering have used the reproducing kernel space method [11–22]. Based on our best knowledge, Cui and Lin are pioneers in linear and nonlinear numerical analysis using the reproducing kernel method (RKM) [23–25].
A novel method for solving a class of second order nonlinear differential equations with finitely many singularities
2015, Applied Mathematics LettersCitation Excerpt :In [4], Umer Saeed and Mujeeb ur Rehman discussed a class of nonlinear boundary value problems by wavelet-Galerkin and quasilinearization. For these references, one can see [5–9]. The numerical treatment of the second order nonlinear differential equation presents some major computational difficulties.