On recent developments in the theory of abstract differential equations with fractional derivatives

This note is motivated from some recent papers treating the problem of the existence of a solution for abstract differential equations with fractional derivatives. We show that the existence results in [Agarwal et al. (2009) [1], Belmekki and Benchohra (2010) [2], Darwish et al. (2009) [3], Hu et al. (2009) [4], Mophou and N’Guérékata (2009) [6] and [7], Mophou (2010) [8] and [9], Muslim (2009) [10], Pandey et al. (2009) [11], Rashid and El-Qaderi (2009) [12] and Tai and Wang (2009) [13]] are incorrect since the considered variation of constant formulas is not appropriate. In this note, we also consider a different approach to treat a general class of abstract fractional differential equations.     

Extending the Newton-Kantorovich hypothesis for solving equations

The famous Newton-Kantorovich hypothesis (Kantorovich and Akilov, 1982 [3], Argyros, 2007 [2], Argyros and Hilout, 2009 [7]) has been used for a long time as a sufficient condition for the convergence of Newton’s method to a solution of an equation in connection with the Lipschitz continuity of the Fréchet-derivative of the operator involved. Here, using Lipschitz and center-Lipschitz conditions, and our new idea of recurrent functions, we show that the Newton-Kantorovich hypothesis can be weakened, under the same information. Moreover, the error bounds are tighter than the corresponding ones given by the dominating Newton-Kantorovich theorem (Argyros, 1998 [1]; [2] and [7]; Ezquerro and Hernández, 2002 [11]; [3]; Proinov 2009, 2010 [16] and [17]).Numerical examples including a nonlinear integral equation of Chandrasekhar-type (Chandrasekhar, 1960 [9]), as well as a two boundary value problem with a Green’s kernel (Argyros, 2007 [2]) are also provided in this study.     

On the Secant method

We present a new semilocal convergence analysis for the Secant method in order to approximate a locally unique solution of a nonlinear equation in a Banach space setting. Our analysis is based on the weaker center-Lipschitz concept instead of the stronger Lipschitz condition which has been ubiquitously employed in other studies such as Amat et al. (2004)  [2], Bosarge and Falb (1969)  [9], Dennis (1971)  [10], Ezquerro et al. (2010)  [11], Hernández et al. (2005, 2000)   and , Kantorovich and Akilov (1982)  [14], Laasonen (1969)  [15], Ortega and Rheinboldt (1970)  [16], Parida and Gupta (2007)  [17], Potra (1982, 1984–1985, 1985)  ,  and , Proinov (2009, 2010)   and , Schmidt (1978) [23], Wolfe (1978)  [24] and Yamamoto (1987)  [25] for computing the inverses of the linear operators. We also provide lower and upper bounds on the limit point of the majorizing sequences for the Secant method. Under the same computational cost, our error analysis is tighter than that proposed in earlier studies. Numerical examples illustrating the theoretical results are also given in this study.     

Linearly perturbed polyhedral normal cone mappings and applications

Under a mild regularity assumption, we derive an exact formula for the Fréchet coderivative and some estimates for the Mordukhovich coderivative of the normal cone mappings of perturbed polyhedra in reflexive Banach spaces. Our focus point is a positive linear independence condition, which is a relaxed form of the linear independence condition employed recently by Henrion et al. (2010) [1], and Nam (2010) [3]. The formulae obtained allow us to get new results on solution stability of affine variational inequalities under linear perturbations. Thus, our paper develops some aspects of the work of Henrion et al. (2010) [1] Nam (2010) [3] Qui (in press) [12] and Yao and Yen (2009) [6] and [7].     

Langevin type limiting processes for adaptive MCMC

Adaptive Markov Chain Monte Carlo (AMCMC) is a class of MCMC algorithms where the proposal distribution changes at every iteration of the chain. In this case it is important to verify that such a Markov Chain indeed has a stationary distribution. In this paper we discuss a diffusion approximation to a discrete time AMCMC. This diffusion approximation is different when compared to the diffusion approximation as in Gelman et al. [5] where the state space increases in dimension to ∞. In our approach the time parameter is sped up in such a way that the limiting process (as the mesh size goes to 0) approaches to a non-trivial diffusion process.     

Radial basis function collocation method for the numerical solution of the two-dimensional transient nonlinear coupled Burgers’ equations

This paper examines the numerical solution of the transient nonlinear coupled Burgers’ equations by a Local Radial Basis Functions Collocation Method (LRBFCM) for large values of Reynolds number (Re) up to 10³. The LRBFCM belongs to a class of truly meshless methods which do not need any underlying mesh but works on a set of uniform or random nodes without any a priori node to node connectivity. The time discretization is performed in an explicit way and collocation with the multiquadric radial basis functions (RBFs) are used to interpolate diffusion-convection variable and its spatial derivatives on decomposed domains. Five nodded domains of influence are used in the local support. Adaptive upwind technique [1] and [54] is used for stabilization of the method for large Re in the case of mixed boundary conditions. Accuracy of the method is assessed as a function of time and space discretizations. The method is tested on two benchmark problems having Dirichlet and mixed boundary conditions. The numerical solution obtained from the LRBFCM for different value of Re is compared with analytical solution as well as other numerical methods [15], [47] and [49]. It is shown that the proposed method is efficient, accurate and stable for flow with reasonably high Reynolds numbers.     

Stability in cellular neural networks with a piecewise constant argument

In this paper, by using the concept of differential equations with piecewise constant arguments of generalized type [1], [2], [3] and [4], a model of cellular neural networks (CNNs) [5] and [6] is developed. The Lyapunov-Razumikhin technique is applied to find sufficient conditions for the uniform asymptotic stability of equilibria. Global exponential stability is investigated by means of Lyapunov functions. An example with numerical simulations is worked out to illustrate the results.     

Kantorovich-type semilocal convergence analysis for inexact Newton methods

We provide a new semilocal convergence analysis for generating an inexact Newton method converging to a solution of a nonlinear equation in a Banach space setting. Our analysis is based on our idea of recurrent functions. Our results are compared favorably to earlier ones by others and us (Argyros (2007, 2009) [5] and [6], Argyros and Hilout (2009) [7], Guo (2007) [15], Shen and Li (2008) [18], Li and Shen (2008) [19], Shen and Li (2009) [20]). Numerical examples are provided to show that our results apply, but not earlier ones [15], [18], [19] and [20].     

A parameter-uniform numerical method for time-dependent singularly perturbed differential-difference equations

A numerical study is made for solving a class of time-dependent singularly perturbed convection–diffusion problems with retarded terms which often arise in computational neuroscience. To approximate the retarded terms, a Taylor’s series expansion has been used and the resulting time-dependent singularly perturbed differential equation is approximated using parameter-uniform numerical methods comprised of a standard implicit finite difference scheme to discretize in the temporal direction on a uniform mesh by means of Rothe’s method and a B-spline collocation method in the spatial direction on a piecewise-uniform mesh of Shishkin type. The method is shown to be accurate of order O(M⁻¹ + N⁻² ln³N), where M and N are the number of mesh points used in the temporal direction and in the spatial direction respectively. An extensive amount of analysis has been carried out to prove the uniform convergence with respect to the singular perturbation parameter. Numerical results are given to illustrate the parameter-uniform convergence of the numerical approximations. Comparisons of the numerical solutions are performed with an upwind and midpoint upwind finite difference scheme on a piecewise-uniform mesh to demonstrate the efficiency of the method.     

Fixed point theorems for some generalized multivalued nonexpansive mappings

In this paper, we introduce a condition on multivalued mappings which is a multivalued version of condition (C_λ) defined by Garcia-Falset et al. (2011) [3]. It is shown here that some of the classical fixed point theorems for multivalued nonexpansive mappings can be extended to mappings satisfying this condition. Our results generalize the results in Lim (1974), Lami Dozo (1973), Kirk and Massa (1990), Garcia-Falset et al. (2011), Dhompongsa et al. (2009) and Abkar and Eslamian (2010) [4], [5], [6], [3], [7] and [8] and many others.     

Laplace–Beltrami enhancement for unstructured two-dimensional meshes having dendritic elements and boundary node movement

The Laplace–Beltrami mesh enhancement algorithm of Hansen et al. ,  and  has been implemented and broadened to include meshes containing dendritic elements and allowing for boundary node movement. This implementation operates on an unstructured two-dimensional mesh by forming an equivalent weak statement using finite element interpolation, assembly, and solution ideas to iteratively place those nodes allowed to move. Moving boundary nodes are constrained to follow the boundary geometry described as a Wilson–Fowler spline (e.g., [3, Section 2.1.3.1]). Implementation details concerning the element basis set modifications, the metric tensor for dendritic element treatment and boundary node movement are presented. Laplacian (e.g., [6]) enhancement is included as a special case. Results are presented which illustrate the algorithm for three test problems.     

Convergence and certain control conditions for hybrid viscosity approximation methods

Very recently, Yao, Chen and Yao [20] proposed a hybrid viscosity approximation method, which combines the viscosity approximation method and the Mann iteration method. Under the convergence of one parameter sequence to zero, they derived a strong convergence theorem in a uniformly smooth Banach space. In this paper, under the convergence of no parameter sequence to zero, we prove the strong convergence of the sequence generated by their method to a fixed point of a nonexpansive mapping, which solves a variational inequality. An appropriate example such that all conditions of this result are satisfied and their condition β_n→0 is not satisfied is provided. Furthermore, we also give a weak convergence theorem for their method involving a nonexpansive mapping in a Hilbert space.     

Majorizing functions and two-point Newton-type methods

The semi-local convergence of a Newton-type method used to solve nonlinear equations in a Banach space is studied. We also give, as two important applications, convergence analyses of two classes of two-point Newton-type methods including a method mentioned in [5] and the midpoint method studied in [1], [2] and [12]. Recently, interest has been shown in such methods [3] and [4].     

Convergence analysis of the LDG method applied to singularly perturbed problems

Considering a two‐dimensional singularly perturbed convection–diffusion problem with exponential boundary layers, we analyze the local discontinuous Galerkin (DG) method that uses piecewise bilinear polynomials on Shishkin mesh. A convergence rate O(N^‐1 lnN) in a DG‐norm is established under the regularity assumptions, while the total number of mesh points is O(N²). The rate of convergence is uniformly valid with respect to the singular perturbation parameter ε. Numerical experiments indicate that the theoretical error estimate is sharp. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

The impact of smooth W-grids in the numerical solution of singular perturbation two-point boundary value problems

This paper develops a semi-analytic technique for generating smooth nonuniform grids for the numerical solution of singularly perturbed two-point boundary value problems. It is based on the usual idea of mapping a uniform grid to the desired nonuniform grid. We introduce the W-grid, which depends on the perturbation parameter ? ? 1. For problems on [0, 1] with a boundary layer at one end point, the local mesh width h_i = x_i+1 − x_i, with 0 = x₀ < x₁ < ? < x_N = 1, is condensed at either 0 or 1. Two simple 2nd order finite element and finite difference methods are combined with the new mesh, and computational experiments demonstrate the advantages of the smooth W-grid compared to the well-known piecewise uniform Shishkin mesh. For small ?, neither the finite difference method nor the finite element method produces satisfactory results on the Shishkin mesh. By contrast, accuracy is vastly improved on the W-grid, which typically produces the nominal 2nd order behavior in L², for large as well as small values of N, and over a wide range of values of ?. We conclude that the smoothness of the mesh is of crucial importance to accuracy, efficiency and robustness.     

Analytic solution of a class of fractional differential equations with variable coefficients by operatorial methods

In this note we show the analytic solution of a class of fractional differential equations with variable coefficients by using operatorial methods. Taking inspiration from previous papers by Dattoli et al. [4], [5] and [6] about spectral properties of Laguerre derivative, we here generalize some of their results to fractional evolution equations. Besides that, we have two interesting generalized examples. One is about telegraph equation with time dependent coefficient. The other, that could be of some interest for realistic applications, is the fractional diffusion with a space-dependent diffusion coefficient.     

Performance enhancement and validation of the open-source software for modeling of ship infrared signatures (OSMOSIS)

Designing stealth technologies for modern warships requires methods for signatures reduction (radar, infrared, etc). First, we have to model these signatures. In Lapierre et al. (2006), Lapierre et al. (2007) [3] and [4], we proposed an infrared (IR) ship signature simulation software (OSMOSIS) that can manage parametric emissivity. OSMOSIS comprises a smart computation time and memory management. Evaluating the effectiveness of an IR signature reduction method implies simulating the evolution of the IR signature with time. This requires a huge computation time. In this paper, we propose a novel hierarchical meshing strategy leading to an important reduction of the computation time. The time complexity of the algorithm is shown to be , where N is the requested number of facets of the mesh of the ship’s surface. Then, we validate OSMOSIS on temperature measurements performed on an “L-shape” object.     

A computational modeling of the behavior of the two-dimensional reaction–diffusion Brusselator system

This paper studies a meshfree technique for the numerical solution of the two-dimensional reaction–diffusion Brusselator system along with Dirichlet and Neumann boundary conditions. Combination of collocation method using the radial basis functions (RBFs) with first order accurate forward difference approximation is employed for obtaining meshfree solution of the problem. Different types of RBFs are used for this purpose. The method is shown to converge to the only equilibrium point of the system. Performance of the proposed method is successfully tested in terms of various error norms. In the case of non-availability of exact solution, performance of the new method is compared with the results obtained from the existing methods [7] and [8]. The elementary stability analysis is established theoretically and is also supported by numerical results.     

A unifying theorem for Newton’s method on spaces with a convergence structure

We present a semilocal convergence theorem for Newton’s method (NM) on spaces with a convergence structure. Using our new idea of recurrent functions, we provide a tighter analysis, with weaker hypotheses than before and with the same computational cost as for Argyros (1996, 1997, 1997, 2007) [1], [2], [3] and [5], Meyer (1984, 1987, 1992) [13], [14] and [15]. Numerical examples are provided for solving equations in cases not covered before.