On the numerical solution of singular two‐point boundary value problems: A domain decomposition homotopy perturbation approach

This paper reports a modified homotopy perturbation algorithm, called the domain decomposition homotopy perturbation method (DDHPM), for solving two‐point singular boundary value problems arising in science and engineering. The essence of the approach is to split the domain of the problem into a number of nonoverlapping subdomains. In each subdomain, a method based on a combination of HPM and integral equation formalism is implemented. The boundary condition at the right endpoint of each inner subdomain is established before deriving an iterative scheme for the components of the solution series. The accuracy and efficiency of the DDHPM are demonstrated by 4 examples (2 nonlinear and 2 linear). In comparison with the traditional HPM, the proposed domain decomposition HPM is highly accurate.     

Existence to singular boundary value problems with sign changing nonlinearities using an approximation method approach

This paper studies the existence of solutions to the singular boundary value problem

A block monotone iterative method for numerical solutions of nonlinear elliptic boundary value problems

The aim of this article is to develop a new block monotone iterative method for the numerical solutions of a nonlinear elliptic boundary value problem. The boundary value problem is discretized into a system of nonlinear algebraic equations, and a block monotone iterative method is established for the system using an upper solution or a lower solution as the initial iteration. The sequence of iterations can be computed in a parallel fashion and converge monotonically to a maximal solution or a minimal solution of the system. Three theoretical comparison results are given for the sequences from the proposed method and the block Jacobi monotone iterative method. The comparison results show that the sequence from the proposed method converges faster than the corresponding sequence given by the block Jacobi monotone iterative method. A simple and easily verified condition is obtained to guarantee a geometric convergence of the block monotone iterations. The numerical results demonstrate advantages of this new approach. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A modified truncated singular value decomposition method for discrete ill‐posed problems

Truncated singular value decomposition is a popular method for solving linear discrete ill‐posed problems with a small to moderately sized matrix A. Regularization is achieved by replacing the matrix A by its best rank‐k approximant, which we denote by A_k. The rank may be determined in a variety of ways, for example, by the discrepancy principle or the L‐curve criterion. This paper describes a novel regularization approach, in which A is replaced by the closest matrix in a unitarily invariant matrix norm with the same spectral condition number as A_k. Computed examples illustrate that this regularization approach often yields approximate solutions of higher quality than the replacement of A by A_k.Copyright © 2014 John Wiley & Sons, Ltd.     

A sub-supersolution approach for Neumann boundary value problems with gradient dependence

Existence and location of solutions to a Neumann problem driven by an nonhomogeneous differential operator and with gradient dependence are established developing a non-variational approach based on an adequate method of sub-supersolution. The abstract theorem is applied to prove the existence of finitely many positive solutions or even infinitely many positive solutions for a class of Neumann problems.     

An upper and lower solution approach for singular boundary value problems with sign changing non‐linearities

We obtain via Schauder's fixed point theorem new results for singular second‐order boundary value problems where our non‐linear term f(t,y,z) is allowed to change sign. In particular, our problem may be singular at y=0, t=0 and/or t=1. Copyright © 2002 John Wiley & Sons, Ltd.     

A Riemannian inexact Newton‐CG method for stochastic inverse singular value problems

In this article, we consider the stochastic inverse singular value problem (ISVP) of constructing a stochastic matrix from the prescribed realizable singular values. We propose a Riemannian inexact Newton‐CG method with various choices of forcing terms for solving the stochastic ISVP. We show the proposed method converges linearly or superlinearly for different forcing terms under some assumptions. We also extend the proposed method to the case of prescribed entries. Finally, we report some numerical results to demonstrate the effectiveness of the proposed method. MOS SUBJECT CLASSIFICATION 65F18; 65F15; 15A18; 65K05; 90C26; 90C48     

Existence of positive solutions for nonlinear singular boundary value problem with p‐Laplacian

In this paper, we study the existence of positive solutions for the following Sturm–Liouville‐like four‐point singular boundary value problem (BVP) with p‐Laplacian where ?_p(s)=|s|^p?2 s, p>1, f is a lower semi‐continuous function. Using the fixed‐point theorem of cone expansion and compression of norm type, the existence of positive solution and infinitely many positive solutions for Sturm–Liouville‐like singular BVP with p‐Laplacian are obtained. Copyright © 2009 John Wiley & Sons, Ltd.     

Self‐adjoint singularly perturbed boundary value problems: an adaptive variational approach

This study is intended to provide a modified variational algorithm for the numerical solution of a class of self‐adjoint singularly perturbed boundary value problem, which is equally applicable to other classes of problems. The principle of the method lies in the introduction of a mixed piecewise domain decomposition and manipulating the variational iterative approach for tackling this class of problems. The uniform convergence of the technique to the exact solution is demonstrated. Numerical results, computational comparisons, suitable error measures and illustrations are provided to testify efficiently and demonstrate the convergence, efficiency and applicability of the method. Copyright © 2012 John Wiley & Sons, Ltd.     

Existence of solutions to nonlinear Neumann boundary value problems with <Emphasis Type="Italic">p</Emphasis>-Laplacian operator and iterative construction

By using some results of pseudo-monotone operator, we discuss the existence and uniqueness of the solution of one kind nonlinear Neumann boundary value problems involving the p-Laplacian operator. We also construct an iterative scheme converging strongly to this solution.     

Extremal solutions for p‐Laplacian boundary value problems with the right‐handed Riemann‐Liouville fractional derivative

The paper is concerned with the solvability for several nonlinear boundary value problems of fractional p‐Laplacian differential equation involving the right‐handed Riemann‐Liouville derivative. By applying monotone iterative technique, lower and upper solutions method and the Banach fixed point theorem, sufficient conditions for existence and uniqueness of extremal solutions are obtained and they extend existing results. At last, two examples are provided to illustrate the results.     

An efficient analytic approach for solving two‐point nonlinear boundary value problems by homotopy analysis method

In this paper, we use homotopy analysis method (HAM) to solve two‐point nonlinear boundary value problems that have at least one solution. The new approach provides the solution in the form of a rapidly convergent series with easily computable components using symbolic computation software. The scheme shows importance of choice of convergence‐control parameter ? to guarantee the convergence of the solutions of nonlinear differential equations. This scheme is tested on three nonlinear exactly solvable differential equations. Two of the examples are practical in science and engineering. The results demonstrate reliability, simplicity and efficiency of the algorithm developed. Copyright © 2011 John Wiley & Sons, Ltd.     

Positive solutions of non‐positone Dirichlet boundary value problems with singularities in the phase variables

The paper presents existence results for positive solutions of the differential equations x ″ + μh (x) = 0 and x ″ + μf (t, x) = 0 satisfying the Dirichlet boundary conditions. Here μ is a positive parameter and h and f are singular functions of non‐positone type. Examples are given to illustrate the main results. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A numerical study of the SVD–MFS solution of inverse boundary value problems in two‐dimensional steady‐state linear thermoelasticity

We study the reconstruction of the missing thermal and mechanical data on an inaccessible part of the boundary in the case of two‐dimensional linear isotropic thermoelastic materials from overprescribed noisy measurements taken on the remaining accessible boundary part. This inverse problem is solved by using the method of fundamental solutions together with the method of particular solutions. The stabilization of this inverse problem is achieved using several singular value decomposition (SVD)‐based regularization methods, such as the Tikhonov regularization method (Tikhonov and Arsenin, Methods for solving ill‐posed problems, Nauka, Moscow, 1986), the damped SVD and the truncated SVD (Hansen, Rank‐deficient and discrete ill‐posed problems: numerical aspects of linear inversion, SIAM, Philadelphia, 1998), whilst the optimal regularization parameter is selected according to the discrepancy principle (Morozov, Sov Math Doklady 7 (1966), 414–417), generalized cross‐validation criterion (Golub et al. Technometrics 22 (1979), 1–35) and Hansen's L‐curve method (Hansen and O'Leary, SIAM J Sci Comput 14 (1993), 1487–503). © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 168–201, 2015     

A layer adaptive B‐spline collocation method for singularly perturbed one‐dimensional parabolic problem with a boundary turning point

In this article, we develop a parameter uniform numerical method for a class of singularly perturbed parabolic equations with a multiple boundary turning point on a rectangular domain. The coefficient of the first derivative with respect to x is given by the formula a₀(x, t)x^p, where a₀(x, t) ≥ α > 0 and the parameter p ∈ [1,∞) takes the arbitrary value. For small values of the parameter ε, the solution of this particular class of problem exhibits the parabolic boundary layer in a neighborhood of the boundary x = 0 of the domain. We use the implicit Euler method to discretize the temporal variable on uniform mesh and a B‐spline collocation method defined on piecewise uniform Shishkin mesh to discretize the spatial variable. Asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular component. These bounds are applied in the convergence analysis of the proposed scheme on Shishkin mesh. The resulting method is boundary layer resolving and has been shown almost second‐order accurate in space and first‐order accurate in time. It is also shown that the proposed method is uniformly convergent with respect to the singular perturbation parameter ε. Some numerical results are given to confirm the predicted theory and comparison of numerical results made with a scheme consisting of a standard upwind finite difference operator on a piecewise uniform Shishkin mesh. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1143–1164, 2011     

A numerical approach for solving weakly singular partial integro‐differential equations via two‐dimensional‐orthonormal Bernstein polynomials with the convergence analysis

In this paper, we develop an efficient matrix method based on two‐dimensional orthonormal Bernstein polynomials (2D‐OBPs) to provide approximate solution of linear and nonlinear weakly singular partial integro‐differential equations (PIDEs). First, we approximate all functions involved in the considerable problem via 2D‐OBPs. Then, by using the operational matrices of integration, differentiation, and product, the solution of Volterra singular PIDEs is transformed to the solution of a linear or nonlinear system of algebraic equations which can be solved via some suitable numerical methods. With a small number of bases, we can find a reasonable approximate solution. Moreover, we establish some useful theorems for discussing convergence analysis and obtaining an error estimate associated with the proposed method. Finally, we solve some illustrative examples by employing the presented method to show the validity, efficiency, high accuracy, and applicability of the proposed technique.