Original ArticlesA boundary shape function iterative method for solving nonlinear singular boundary value problems
Introduction
The singular boundary value problem (BVP) occurs in many areas for the applications in science and engineering, including thermal explosions [3], tumor growth models [1], electroosmotic flows [4], [5], [6], modeling of heat sources in human head [13], oxygen diffusion [26], physiology [15], electrohydrodynamic flow of a fluid [36], viscous flows [40], science and engineering [32], as well as various physical models [37].
The current paper designs a novel numerical method to find the singular solution of nonlinear BVP, exactly satisfying the prescribed boundary conditions on the boundaries of a unit interval. When a singular point happens at the boundary, it might be a difficult task to solve this sort problem. Many computational methods have been developed for solving the nonlinear BVPs of the singular type [2], [12], [15], [19], [29], [30], [33], [34].
The boundary shape function method (BSFM) was first introduced by Liu and Chang [24] to find the periodic solution of the nonlinear jerk equations, and by Liu et al. [25] to solve the optimal control problems of nonlinear Duffing oscillators. Later, the methodology of the BSFM was extended to solve the BVP with multipoint boundary conditions [22] and the singularly perturbed BVP with Robin boundary conditions [23]. The present problem with nonlinear BVP and with a singular point at the boundary is more difficult to be solved.
We arrange the paper as follows. First, we give a brief sketch of the shooting method to solve the nonlinear singular BVP in Section 2, which is recast to be an initial value problem (IVP) but with unknown initial conditions. In Section 3, we introduce the boundary shape function (BSF) for satisfying the specified Robin boundary conditions automatically. In Section 4, we take the BSF as the solution, such that we can exactly transform the nonlinear singular BVP to an IVP for the new variable with known initial conditions, and from it the fulfillment of the Robin boundary conditions is automatic for the original variable. The iterative algorithm is developed to find the unknown terminal values of the new variable, which can be iteratively determined with a fast convergence speed. In Section 5, we solve some numerical examples by the proposed iterative algorithm based on the BSFM and compare to that obtained from the shooting method. In Section 6, using the BSFM we solve the Blasius equation and a nonlinear singular differential equation of a pressurized spherical membrane with a strong singularity. Finally, the conclusions are drawn in Section 7.
Section snippets
Shooting method
In the fields of science and engineering, some axisymmetrical diffusion problems in the steady state can be modeled by a class of two-point BVPs [8], [9]: where and . If , we encounter a nonlinear singular boundary value problem. With certain conditions on , the existence and uniqueness of the solution of Eqs. (1)–(3) have been discussed in [27], [28], [31].
We
Boundary shape function
The following results are required to construct the boundary shape function [23].
Theorem 1 are shape functions and satisfy
Proof Refer [23] for the proof of the existence of and . □
Theorem 2 Suppose that is a free function. With Eqs. (13) and (14) for and and there exists a boundary shape function which automatically
The boundary shape function method
In Theorem 2, we can replace by as a boundary shape function to automatically satisfy the Robin boundary conditions (2), (3). is defined by Eq. (15) and thus the variable transformation from to a new variable is given by which is inserted into Eq. (1) to obtain a new ordinary differential equation (ODE) for : We place in to stress that they are unknown parameters in the
Numerical tests
We compare the shooting method in Section 2 and the BSFM in Section 4 with the following examples.
Two practical nonlinear singular problems
In the application of the BSFM to solve practical nonlinear singular problems, we may encounter the situation which is more complex than that in Eq. (1). In such problems the singularities with and and at the boundary point are more difficult to be treated than the singularity appeared in Eq. (1) at the boundary point . We solve two such problems in the below by using the BSFM.
Conclusions
For the second-order nonlinear singular boundary value problems, the boundary shape functions, exactly satisfying the given Robin boundary conditions, are derived. Because the singularity is appeared on the boundary, the numerical method is designed to exactly and automatically preserve the Robin boundary conditions. In the shooting method, we have applied the fictitious time integration method to solve the residual equation for satisfying the right-boundary condition. The present paper derived
References (40)
- et al.
Mathematical model of tumor growth IV: Effect of necrotic core
Math. Biosci.
(1989) Electroosmotic flow in slit microchannel containing salt-free solution
Eur. J. Mech. B Fluids
(2010)Electroosmotic flow in a dissimilarly charged slit microchannel containing salt-free solution
Eur. J. Mech. B Fluids
(2012)Taylor series method for solving a class of nonlinear singular boundary value problems arising in applied science
Appl. Math. Comput.
(2014)Convergence of variational iteration method applied to two-point diffusion problems
Appl. Math. Model.
(2016)Numerical solutions of the classical Blasius flat-plate problem
Appl. Math. Comput.
(2005)Integral methods for computing solutions of a class of singular two-point boundary value problems
Appl. Math. Comput.
(2002)The distribution of heat sources in the human head: A theoretical consideration
J. Theoret. Biol.
(1975)An iterative solution for rotationally symmetric nonlinear membrane problems
Int. J. Non-Linear Mech.
(1966)- et al.
A novel approach for the solution of a class of singular boundary value problems arising in physiology
Math. Comput. Modelling
(2010)
Modified Adomian decomposition method and computer implementation for solving singular boundary value problems arising in various physical problems
Comput. Chem. Eng.
The Lie-group shooting method for boundary-layer problems with suction/injection/ reverse flow conditions for power-law fluids
Int. J. Non-Linear Mech.
Boundary shape function method for nonlinear BVP, automatically satisfying prescribed multipoint boundary conditions
Bound. Value Probl.
The periods and periodic solutions of nonlinear jerk equations solved by an iterative algorithm based on a shape function method
Appl. Math. Lett.
A re-examination of oxygen diffusion in a spherical cell with Michaelis–Menten nonlinear oxygen uptake kinetics
J. Theoret. Biol.
On a class of weakly regular singular two-point boundary value problems, II
J. Differential Equations
On a class of regular singular two-point boundary value problems
J. Math. Anal. Appl.
On the convergence of a finite difference method for a class of singular boundary value problems arising in physiology
J. Comput. Appl. Math.
On the convergence of a fourth-order method for a class of singular boundary value problems
J. Comput. Appl. Math.
Existence-uniqueness results for a class of singular boundary value problems arising in physiology
Nonlinear Anal. RWA
Cited by (11)
An explicit representation of the three-point Hermite interpolant for the numerical solution of singular boundary value problems
2022, Applied Numerical MathematicsCitation Excerpt :Wavelet Galerkin method (WGM) was proposed by Nosrati Sahlan and Hashemizadeh [51] to approximate the solution of a class of practical SBVPs. Furthermore, the simplified reproducing kernel method (SRKM) [32], Haar wavelet collocation method (HWCM) [53], B-spline functions method (BSM) [4], boundary shape function method (BSFM) [28], He's variational iteration method (HVIM) [20], neural networks [49,50], and many other numerical techniques have been recently utilized to solve different types of nonlinear SBVPs [18,25,26,40,48,56,58,59,61,62]. In addition to the analytical importance of the Hermite interpolation, the widespread utilization of this interpolation method in solving a variety of differential equations demonstrates its practical importance [16,17,27,29,46,55].
Boundary shape function iterative method for nonlinear second-order boundary value problems with nonlinear boundary conditions
2022, Mathematics and Computers in SimulationAn adaptive pair of one-step hybrid block Nyström methods for singular initial-value problems of Lane–Emden–Fowler type
2022, Mathematics and Computers in SimulationCitation Excerpt :Many researchers have widely studied the model in (1) and provided different strategies to solve it. To indicate some of the existing methods to give numerical solutions to the SSIVPs given in (1) and related problems, we can mention the implicit Euler method proposed by Koch et al. [7], the homotopy-perturbation technique by [5], the pseudospectral method presented by Mehrpouya [10], the explicit one-step strategies reported by Kutniv et al. [8], the collocation method by Bhrawy and Alofi [4], the Haar Wavelet resolution technique in [24], the heuristic approach in [22], an iterative method in [9], the quasilinearization method by Singh et al. [23], the hybrid block Nyström methods in [20], the analytical methods presented in [26], the Adomian decomposition method in [6], the spline techniques in [16,17] or a compact finite difference method in [18] and the references enclosed in those manuscripts. Another possibility to deal with (1) is to transform the differential equation in an equivalent system of first order differential equations and use any of the methods for this kind of problems, as the one in [11], although this is usually more expensive computationally.
EXISTENCE AND STABILITY OF SOLUTIONS FOR A COUPLED HADAMARD TYPE SEQUENCE FRACTIONAL DIFFERENTIAL SYSTEM ON GLUCOSE GRAPHS
2024, Journal of Applied Analysis and ComputationWavelets based physics informed neural networks to solve non-linear differential equations
2023, Scientific ReportsSolving nonlinear boundary value problems by a boundary shape function method and a splitting and linearizing method
2023, International Journal of Nonlinear Sciences and Numerical Simulation