A posteriori H1 error estimation for a singularly perturbed reaction diffusion problem on anisotropic meshes

The paper deals with a singularly perturbed reaction diffusionmodel problem. The focus is on reliable a posteriori error estimatorsfor the H¹ seminorm that can be applied to anisotropic finiteelement meshes. A residual error estimator and a local problemerror estimator are proposed and rigorously analysed. They arelocally equivalent, and both bound the error reliably. Threemodifications of these estimators are introduced and discussed. Much attention is given to the performance of the error estimatorin numerical experiments. This helps to identify those estimatorsthat are suitable for practical applications.     

A mixed hp FEM for the approximation of fourth‐order singularly perturbed problems on smooth domains

We consider fourth‐order singularly perturbed problems posed on smooth domains and the approximation of their solution by a mixed Finite Element Method on the so‐called Spectral Boundary Layer Mesh. We show that the method converges uniformly, with respect to the singular perturbation parameter, at an exponential rate when the error is measured in the energy norm. Numerical examples illustrate our theoretical findings.     

Wavelet multiresolution interpolation Galerkin method for nonlinear boundary value problems with localized steep gradients

The wavelet multiresolution interpolation for continuous functions defined on a finite interval is developed in this study by using a simple alternative of transformation matrix. The wavelet multiresolution interpolation Galerkin method that applies this interpolation to represent the unknown function and nonlinear terms independently is proposed to solve the boundary value problems with the mixed Dirichlet-Robin boundary conditions and various nonlinearities, including transcendental ones, in which the discretization process is as simple as that in solving linear problems, and only common two-term connection coefficients are needed. All matrices are independent of unknown node values and lead to high efficiency in the calculation of the residual and Jacobian matrices needed in Newton’s method, which does not require numerical integration in the resulting nonlinear discrete system. The validity of the proposed method is examined through several nonlinear problems with interior or boundary layers. The results demonstrate that the proposed wavelet method shows excellent accuracy and stability against nonuniform grids, and high resolution of localized steep gradients can be achieved by using local refined multiresolution grids. In addition, Newton’s method converges rapidly in solving the nonlinear discrete system created by the proposed wavelet method, including the initial guess far from real solutions.

     

On the Order Reduction Approach for Singularly Perturbed Optimal Control Systems

The order reduction method for singularly perturbed optimal control systems consists of employing the system obtained while setting the small parameter to be zero. In many situations the differential-algebraic system thus obtained indeed provides an appropriate approximation to the singularly perturbed problem with a small parameter. In this paper we establish that if relaxed controls are allowed then the answer to the question whether or not this method is valid depends essentially on one simple parameter: the dimension of the fast variable, denoted n. More specifically, if n=1 then the order reduction method is indeed applicable, while if n>1 then the set of singularly perturbed optimal control systems for which it is not applicable is dense (in the L norm).     

Schrödinger Operators with Empty Singularly Continuous Spectra

Let H be a semibounded perturbation of the Laplacian H ₀ in L ²( ^d ). For an admissible function sufficient conditions are given for the completeness of the scattering system (H), (H ₀). If is the exponential function and if e^{– H} is an integral operator we denote the kernel of the difference D = e^{– H} – e^{– H} ₀ by D (x, y), > 0. The singularly continuous spectrum of H is empty if^d dx ^d dy |D(x,y)| (1 + |y|²)< for some > 1. This result is applied to potential perturbations and to perturbations by imposing Dirichlet boundary conditions.     

Concentrating solutions of some singularly perturbed elliptic equations

We study singularly perturbed elliptic equations arising from models in physics or biology, and investigate the asymptotic behavior of some special solutions. We also discuss some connections with problems arising in differential geometry.      

Simple eigenvalues for the Steklov problem in a domain with a small hole. A functional analytic approach

Let be a bounded open domain of . Let denote the outward unit normal of . We assume that the Steklov problem Δu = 0 in and on has a simple eigenvalue of . Then we consider an annular domain obtained by removing from a small‐cavity size of ε > 0, and we show that under proper assumptions there exists a real valued and real analytic function defined in an open neighborhood of (0,0) in and such that is a simple eigenvalue for the Steklov problem Δu = 0 in and on for all ε > 0 small enough, and such that . Here denotes the outward unit normal of , and δ_2,2 ≡ 1 and δ_2,n ≡ 0 if n ≥ 3. Then related statements have been proved for corresponding eigenfunctions. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical solutions based on Chebyshev collocation method for singularly perturbed delay parabolic PDEs

The main purpose of this work is to provide a numerical approach for the delay partial differential equations based on a spectral collocation approach. In this research, a rigorous error analysis for the proposed method is provided. The effectiveness of this approach is illustrated by numerical experiments on two delay partial differential equations. Copyright © 2013 John Wiley & Sons, Ltd.     

Uniform convergence analysis of hybrid numerical scheme for singularly perturbed problems of mixed type

In this article, we consider a class of singularly perturbed mixed parabolic‐elliptic problems whose solutions possess both boundary and interior layers. To solve these problems, a hybrid numerical scheme is proposed and it is constituted on a special rectangular mesh which consists of a layer resolving piecewise‐uniform Shishkin mesh in the spatial direction and a uniform mesh in the temporal direction. The domain under consideration is partitioned into two subdomains. For the spatial discretization, the proposed scheme is comprised of the classical central difference scheme in the first subdomain and a hybrid finite difference scheme in the second subdomain, whereas the time derivative in the given problem is discretized by the backward‐Euler method. We prove that the method converges uniformly with respect to the perturbation parameter with almost second‐order spatial accuracy in the discrete supremum norm. Numerical results are finally presented to validate the theoretical results.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1931–1960, 2014