Extending the applicability of Newton’s method for a class of boundary value problems using the shooting method
Introduction
In many scientific and engineering problems the main objective is to solve an equation. Equivalently, the principal aim may be to find a good local approximation x* for the equation:using the determined Fréchet–differentiable operator f defined in a nonempty, open and convex domain E of a Banach space T1 with correspondence in a Banach space T2.
We have focused, in this work, on problems such as the following:with
Although there are equations appearing in these type of problems that can be solved by applying an analytical way, the solution is often approximated numerically through an iterative method due to its high level of complexity. Normally, Newton’s method is used due to its easily, efficient and simple implementation.
A common technique used to solve problems involving differential equations with boundary conditions is the so called shooting method. This method acquires this name because of the desire, to arrive to the boundary final point, starting from the point on the frontier boundary, or in case of it not being exact, to approximate it as much as possible. Other methods can also be used to obtain solutions. The main complication with this method is that in order to ensure convergence, the steps involved must be deeply studied. There is also the possibility that it may not even be convergent.
The type of convergence usually analyzed is the semilocal type, which is focused on the information around a starting point in order to turn on conditions of convergence to the iterative process towards the sought solution. There is a plethora of studies that analyze Newton’s method as in Amat et al. [1], Argyros and Hilout [3], Argyros and Magreñán [4], Burden et al. [6], Kantorovich and Akilov [7], Kelley [9], 18 [10], Parhi and Gupta [11], Rheinboldt [12], Sánchez [13], Traub [14].
The problem appearing in this method and in others is that the starting point is needed, so that it approximates the solution, being such starting point unknown at first. Once conditions are given on this point, there are specific conditions that can ensure convergence.
The layout of this work is: in Section 2 the shooting method of resolution for problems with border values is detailed. A study of the general convergence of Newton’s method that appears in these problems is performed as well depending on whether the function involved is increasing or not. A semilocal and local convergence analysis is presented for Newton’s method involving Banach space valued operators once the method is defined, in Section 3. The convergence criteria are improved. Finally, in Section 4, examples of the new developed theory are included.
Section snippets
Shooting method
Consider the boundary value problemwhere is a real valued function and is an arbitrary function. The shooting method is a usual way to approximate a solution of a boundary value problem like Eq. (2.1). Here, we fix an initial aproximation of the initial slope, α1, and we compute the solution to obtain (see Fig. 1).
We have several techniques if the value of is larger than desired. One method, for example requires
Convergence analysis
Consider H: E ⊂ T1 → T2 to be continuously differentiable in the Fréchet sense.
We present the semi–local followed by the local convergence analysis of Newton’s method Definition 3.1 We say that operator H′ satisfies the center–Lipschitz condition on E, if there exists L0 > 0 such that for each v ∈ E, and some fixed v0 ∈ ESuppose that μ > 0 is such that . Define the set . Definition 3.2 We say that operator H′ satisfies the restricted–Lipschitz condition
Numerical examples
In this section we show several examples where the previous conditions are analyzed. Example 4.1 As an academic and motivational example, let us set . Consider function ϕ asUsing (4.1), we get and . Notice that for each The Kantorovich convergence criterion (3.12) is not satisfied, sinceHowever, our condition (3.9) gives Example 4.2 Let and .
References (14)
- et al.
A modified Chebyshev’s iterative method with at least sixth order of convergence
Appl. Math. Comput.
(2008) A new tool to study real dynamics: the convergence plane
Appl. Math. Comput.
(2014)Convergence and Applications of Newton–Type Iterations
(2008)- et al.
Computational Methods in Nonlinear analysis
(2013) - et al.
A Contemporary Study of Iterative Method
(2018) - et al.
Computer Methods for Ordinary Differential Equations and Differential Algebraic Eccuations
(1998) - et al.
Numerical Analysis
(2016)
Cited by (5)
Modelling study on freezing process of water droplet on inclined cold plate surface with droplet dynamic behavior considered
2022, International Journal of Heat and Mass TransferCitation Excerpt :Inputting appropriate value of P0, the solution of Eq. (15) is obtained numerically by using the Runge-Kutta method [35] under the initial conditions of Eq. (17a) and (17b). P0 is corrected repeatedly with Newton's method [36] to satisfy Eq. (17c) and (17d). Finally, the droplet profile under spreading type can be obtained.
Adaptive simulations to pressure distribution for creeping motion of Carreau nanofluid with variable fluid density effects: Physiological applications
2022, Thermal Science and Engineering ProgressAiming at methods' wider adoption: Applicability determinants and metrics
2021, Computer Science ReviewCitation Excerpt :This is especially because the search engines transformed the term into its declinational forms including such popular words as ‘application’ or ‘applied’. Moreover, the rough analysis of results revealed that the majority of them concerned discussions on the application of a proposed solution to a particular domain [13–26]. The results of the core search are summarised in Tables 1 and 2.
Concentration-dependent electrical and thermal conductivity effects on magnetoHydrodynamic Prandtl nanofluid in a divergent–convergent channel: Drug system applications
2022, Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering