SUPERCONVERGENCE OF GRADIENT RECOVERY SCHEMES ON GRADED MESHES FOR CORNER SINGULARITIES

For the linear finite element solution to the Poisson equation, we show that supercon- vergence exists for a type of graded meshes for corner singularities in polygonal domains. In particular, we prove that the L^2-projection from the piecewise constant field △↓UN to the continuous and piecewise linear finite element space gives a better approximation of △↓U in the Hi-norm. In contrast to the existing superconvergence results, we do not assume high regularity of the exact solution.     

High order finite volume methods for singular perturbation problems

In this paper we establish a high order finite volume method for the fourth order singular perturbation problems.In conjunction with the optimal meshes,the numerical solutions resulting from the method have optimal convergence order.Numerical experiments are presented to verify our theoretical estimates.     

A higher-order macroscopic model for pedestrian flows

This paper develops a higher-order macroscopic model of pedestrian crowd dynamics derived from fluid dynamics that consists of two-dimensional Euler equations with relaxation. The desired directional motion of pedestrians is determined by an Eikonal-type equation, which describes a problem that minimizes the instantaneous total walking cost from origin to destination. A linear stability analysis of the model demonstrates its ability to describe traffic instability in crowd flows. The algorithm to solve the macroscopic model is composed of a splitting technique introduced to treat the relaxation terms, a second-order positivity-preserving central-upwind scheme for hyperbolic conservation laws, and a fast-sweeping method for the Eikonal-type equation on unstructured meshes. To test the applicability of the model, we study a challenging pedestrian crowd flow problem of the presence of an obstruction in a two-dimensional continuous walking facility. The numerical results indicate the rationality of the model and the effectiveness of the computational algorithm in predicting the flux or density distribution and the macroscopic behavior of the pedestrian crowd flow. The simulation results are compared with those obtained by the two-dimensional Lighthill-Whitham-Richards pedestrian flow model with various model parameters, which further shows that the macroscopic model is able to correctly describe complex phenomena such as “stop-and-go waves” observed in empirical pedestrian flows.     

A brief survey on numerical methods for solving singularly perturbed problems

In the present paper, a brief survey on computational techniques for the different classes of singularly perturbed problems is given. This survey is a continuation of work performed earlier by the first author and contains the literature of the work done by the researchers during the years 2000-2009. However some older important relevant papers are also included in this survey. We also mentioned those papers which are not surveyed in the previous survey papers by the first author of this paper, see [Appl. Math. Comput. 30 (1989) 223-259, 130 (2002) 457-510, 134 (2003) 371-429] for details. Thus this survey paper contains a surprisingly large amount of literature on singularly perturbed problems and indeed can serve as an introduction to some of the ideas and methods for the singular perturbation problems.     

Finite element simulation of compressible particle-laden gas flows

A macroscopic two-fluid model of compressible particle-laden gas flows is considered. The governing equations are discretized by a high-resolution finite element method based on algebraic flux correction. A multidimensional limiter of TVD type is employed to constrain the local characteristic variables for the continuous gas phase and conservative fluxes for a suspension of solid particles. Special emphasis is laid on the efficient computation of steady state solutions at arbitrary Mach numbers. To avoid stability restrictions and convergence problems, the characteristic boundary conditions are imposed weakly and treated in a fully implicit manner. A two-way coupling via the interphase drag force is implemented using operator splitting. The Douglas-Rachford scheme is found to provide a robust treatment of the interphase exchange terms within the framework of a fractional-step solution strategy. Two-dimensional simulation results are presented for a moving shock wave and for a steady nozzle flow.     

Morley type non‐C0 nonconforming rectangular plate finite elements on anisotropic meshes

In this article, two Morley type non‐C⁰ nonconforming rectangular finite elements are discussed to numerically solve the fourth order plate bending problem under anisotropic meshes. The optimal anisotropic interpolation error and consistency error estimates are obtained by using some novel approaches. Some numerical tests are given to confirm the theoretical analysis. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

The adaptive immersed interface finite element method for elliptic and Maxwell interface problems

We propose self-adaptive finite element methods with error control for solving elliptic and electromagnetic problems with discontinuous coefficients. The meshes in the methods do not need to fit the interfaces. New error indicators are introduced to control the error due to non-body-fitted meshes. Flexible h-adaptive strategies are developed, which can be systematically extended to a large class of interface problems. Extensive numerical experiments are performed to support the theoretical results and to show the competitive behavior of the adaptive algorithm even for interfaces involving corner or tip singularities.     

High contrast multi-pass Fabry-Perot interferometer

A technique for improving the contrast ratio of a Fabry-Perot spectrometer is presented. Using triplepasses of radiation through a mesh-reflector Fabry-Perot it is shown how the contrast ratio may be improved by orders of magnitude. Experimental verification has been demonstrated and the application of the technique to specific cases of interest discussed briefly. Overall improvement in finesse as well as contrast ratio has been achieved with only small loss in transmission.     

Locally Supported,Piecewise Polynomial Biorthogonal Wavelets on Nonuniform Meshes

In this paper, biorthogonal wavelets are constructed on nonuniform meshes. Both primal and dual wavelets are locally supported, continuous piecewise polynomials. The wavelets generate Riesz bases for the Sobolev spaces (H^s) for (|s| < 3/2). The wavelets at the primal side span standard Lagrange finite element spaces.