Comparison between the convergence rates of the Chebyshev method and the related (2,2)‐step methods

An optimal Chebyshev method for solving A x = b , where all the eigenvalues of the real and non‐symmetric matrix A are located in the open right half plane, is dependent on an optimal ellips∂Ω^* such that the spectrum of A is contrained in Ω^*, the closed interior of the ellipse. The relationship between the convergence rates of the Chebyshev method and the closely related (2,2)‐step iterative methods are studied. (2,2)‐step iterative methods are faster than an optimal Chebyshev method under certain conditions. A numerical example illustrates such an improvement of a (2,2)‐step iterative method. Copyright © 2000 John Wiley & Sons, Ltd.     

Strong convergence of the split‐step θ‐method for stochastic age‐dependent population equations

In this paper, we constructed the split‐step θ (SSθ)‐method for stochastic age‐dependent population equations. The main aim of this paper is to investigate the convergence of the SS θ‐method for stochastic age‐dependent population equations. It is proved that the proposed method is convergent with strong order 1/2 under given conditions. Finally, an example is simulated to verify the results obtained from the theory, and comparative analysis with Euler method is given, the results show the higher accuracy of the SS θ‐method. Copyright © 2014 John Wiley & Sons, Ltd.     

Determination of space‐dependent coefficients from temperature measurements using the conjugate gradient method

In this article, we consider coefficient identification problems in heat transfer concerned with the determination of the space‐dependent perfusion coefficient and/or thermal conductivity from interior temperature measurements using the conjugate gradient method (CGM). We establish the direct, sensitivity and adjoint problems and the iterative CGM algorithm which has to be stopped according to the discrepancy principle in order to reconstruct a stable solution for the inverse problem. The Sobolev gradient concept is introduced in the CGM iterative algorithm in order to improve the reconstructions. The numerical results illustrated for both exact and noisy data, in one‐ and two‐dimensions for single or double coefficient identifications show that the CGM is an efficient and stable method of inversion.     

On convergence and complexity of the modified forward‐backward method involving new linesearches for convex minimization

In optimization theory, convex minimization problems have been intensively investigated in the current literature due to its wide range in applications. A major and effective tool for solving such problem is the forward‐backward splitting algorithm. However, to guarantee the convergence, it is usually assumed that the gradient of functions is Lipschitz continuous and the stepsize depends on the Lipschitz constant, which is not an easy task in practice. In this work, we propose the modified forward‐backward splitting method using new linesearches for choosing suitable stepsizes and discuss the convergence analysis including its complexity without any Lipschitz continuity assumption on the gradient. Finally, we provide numerical experiments in signal recovery to demonstrate the computational performance of our algorithm in comparison to some well‐known methods. Our reports show that the proposed algorithm has a good convergence behavior and can outperform the compared methods.     

Finite element analysis and approximation of Burgers’‐Fisher equation

In this article, the authors present finite element analysis and approximation of Burgers’‐Fisher equation. Existence and uniqueness of weak solution is proved by Galerkin's finite element method for non‐smooth initial data. Next, a priori error estimates of semi‐discrete solution in norm, are derived and the convergence of semi‐discrete solution is established. Then, fully discretization of the problem is done with the help of Euler's backward method. The nonlinearity is removed by lagging it to previous known level. The scheme is found to be convergent. Positivity of fully discrete solution is discussed, and bounds on time step are discovered for which the solution preserves its positivity. Finally, numerical experiments are performed on some examples to demonstrate the effectiveness of the scheme. The proposed scheme found to be fast, easy and accurate.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1652–1677, 2017     

A nonparametric approach to 3D shape analysis from digital camera images — I

This article for the first time develops a nonparametric methodology for the analysis of projective shapes of configurations of landmarks on real 3D objects from their regular camera pictures. A fundamental result in computer vision, emulating the principle of human vision in space, claims that, generically, a finite 3D configuration of points can be retrieved from corresponding configurations in a pair of camera images, up to a projective transformation. Consequently, the projective shape of a 3D configuration can be retrieved from two of its planar views, and a projective shape analysis can be pursued from a sample of images. Projective shapes are here regarded as points on projective shape manifolds. Using large sample and nonparametric bootstrap methodology for extrinsic means on manifolds, one gives confidence regions and tests for the mean projective shape of a 3D configuration from its 2D camera images. Two examples are given: an example of testing for accuracy of a simple manufactured object using mean projective shape analysis, and a face identification example. Both examples are data driven based on landmark registration in digital images.     

Some new nonlinear wave solutions and their convergence for the (2+1)‐dimensional Broer–Kau–Kupershmidt equation

We use the bifurcation method of dynamical systems to study the (2+1)‐dimensional Broer–Kau–Kupershmidt equation. We obtain some new nonlinear wave solutions, which contain solitary wave solutions, blow‐up wave solutions, periodic smooth wave solutions, periodic blow‐up wave solutions, and kink wave solutions. When the initial value vary, we also show the convergence of certain solutions, such as the solitary wave solutions converge to the kink wave solutions and the periodic blow‐up wave solutions converge to the solitary wave solutions. Copyright © 2014 John Wiley & Sons, Ltd.     

Block‐centered upwind multistep difference method and convergence analysis for numerical simulation of oil reservoir

The three‐dimensional displacement of two‐phase flow in porous media is a preliminary problem of numerical simulation of energy science and mathematics. The mathematical model is formulated by a nonlinear system of partial differential equations to describe incompressible miscible case. The pressure is defined by an elliptic equation, and the concentration is defined by a convection‐dominated diffusion equation. The pressure generates Darcy velocity and controls the dynamic change of concentration. We adopt a conservative block‐centered scheme to approximate the pressure and Darcy velocity, and the accuracy of Darcy velocity is improved one order. We use a block‐centered upwind multistep method to solve the concentration, where the time derivative is approximated by multistep method, and the diffusion term and convection term are treated by a block‐centered scheme and an upwind scheme, respectively. The composite algorithm is effective to solve such a convection‐dominated problem, since numerical oscillation and dispersion are avoided and computational accuracy is improved. Block‐centered method is conservative, and the concentration and the adjoint function are computed simultaneously. This physical nature is important in numerical simulation of seepage fluid. Using the convergence theory and techniques of priori estimates, we derive optimal estimate error. Numerical experiments and data show the support and consistency of theoretical result. The argument in the present paper shows a powerful tool to solve the well‐known model problem.     

A study of integrability and symmetry for the (p + 1)th Boltzmann equation via Painlevé analysis and Lie‐group method

In this paper,we applied the Painlevé property test on Krook‐Wu model of the nonlinear Boltzmann equation (p = 1). As a result, by using Bäcklund transformation, we obtained three solutions two of them were known earlier, while the third one is new and more general, we have also two reductions one of them is Abel's equation. Also, Lie‐group method is applied to the (p + 1)th Boltzmann equation. The complete Lie algebra of infinitesimal symmetries is established. Three nonequivalent sub‐algebraic of the complete Lie algebra are used to investigate similarity solutions and similarity reductions in the form of nonlinear ordinary equations for (p + 1)th Boltzmann equation; we obtained two general solutions for (p + 1)th Boltzmann equation and new solutions for Krook‐Wu model of Boltzmann equation (p = 1). Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling and analysis of system reliability using phase‐type distribution closure properties

Phase‐type distribution closure properties are utilized to devise algorithms for generating reliability functions of systems with basic structures. These structures include series, parallel, K‐out‐of‐N, and standby structures with perfect/imperfect switch. The algorithms form a method for system reliability modeling and analysis based on the relationship between the system lifetime and component lifetimes for general structures. The proposed method is suitable for functional system reliability analysis, which can produce reliability functions of systems with independent components instead of only system reliability values. Once the system reliability function is obtained, other reliability measures such as the system's hazard function and mean time to failure can be obtained efficiently using only matrix algebra. Dimensional and numerical comparisons with computerized symbolic processing are also presented to show the superiority of the proposed method.     

Numerical analysis of a finite volume scheme for a seawater intrusion model with cross‐diffusion in an unconfined aquifer

We consider a degenerate parabolic system modeling the flow of fresh and saltwater in a porous medium in the context of seawater intrusion. We propose and analyze a finite volume scheme based on two‐point flux approximation with upwind mobilities. The scheme preserves at the discrete level the main features of the continuous problem, namely the nonnegativity of the solutions, the decay of the energy and the control of the entropy and its dissipation. Based on these nonlinear stability results, we show that the scheme converges toward a weak solution to the problem. Numerical results are provided to illustrate the behavior of the model and of the scheme.     

A new technique of using homotopy analysis method for solving high‐order nonlinear differential equations

In this paper, a new technique of homotopy analysis method (HAM) is proposed for solving high‐order nonlinear initial value problems. This method improves the convergence of the series solution, eliminates the unneeded terms and reduces time consuming in the standard homotopy analysis method (HAM) by transform the nth‐order nonlinear differential equation to a system of n first‐order equations. Second‐ and third‐ order problems are solved as illustration examples of the proposed method. Copyright © 2010 John Wiley & Sons, Ltd.