High order generalized upwind schemes and numerical solution of singular perturbation problems

High even order generalizations of the traditional upwind method are introduced to solve second order ODE-BVPs without recasting the problem as a first order system. Both theoretical analysis and numerical comparison with central difference schemes of the same order show that these new methods may avoid typical oscillations and achieve high accuracy. Singular perturbation problems are taken into account to emphasize the main features of the proposed methods. AMS subject classification (2000)  65L10, 65L12, 65L50     

On the solution of Maxwell's equations in polygonal domains

This paper is concerned with the structure of the singular and regular parts of the solution of time‐harmonic Maxwell's equations in polygonal plane domains and their effective numerical treatment. The asymptotic behaviour of the solution near corner points of the domain is studied by means of discrete Fourier transformation and it is proved that the solution of the boundary value problem does not belong locally to H² when the boundary of the domain has non‐acute angles. A splitting of the solution into a regular part belonging to the space H², and an explicitly described singular part is presented. For the numerical treatment of the boundary value problem, we propose a finite element discretization which combines local mesh grading and the singular field methods and derive a priori error estimates that show optimal convergence as known for the classical finite element method for problems with regular solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

Recognition of the Finite Simple Groups F 4(2m) by Spectrum

The spectrum of a finite group is the set of its element orders. A finite group G is said to be recognizable by spectrum, if every finite group with the same spectrum as G is isomorphic to G. The purpose of the paper is to prove that for every natural m the finite simple Chevalley group F ₄(2 ^m ) is recognizable by spectrum.     

Vassiliev invariants for braids on surfaces

We show that Vassiliev invariants separate braids on a closed oriented surface, and we exhibit a universal Vassiliev invariant for these braids in terms of chord diagrams labeled by elements of the fundamental group of the surface.

     

Mass‐lumping or not mass‐lumping for eigenvalue problems

In this article we analyze the effect of mass‐lumping in the linear triangular finite element approximation of second‐order elliptic eigenvalue problems. We prove that the eigenvalue obtained by using mass‐lumping is always below the one obtained with exact integration. For singular eigenfunctions, as those arising in non convex polygons, we prove that the eigenvalue obtained with mass‐lumping is above the exact eigenvalue when the mesh size is small enough. So, we conclude that the use of mass‐lumping is convenient in the singular case. When the eigenfunction is smooth several numerical experiments suggest that the eigenvalue computed with mass‐lumping is below the exact one if the mesh is not too coarse. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 653–664, 2003     

A boundary element method for the calculation of noise barrier insertion loss in the presence of atmospheric turbulence

Atmospheric turbulence is an important factor that limits the amount of attenuation a barrier can provide in the outdoor environment. It is therefore important to develop a reliable method to predict its effect on barrier performance. The boundary element method (BEM) has been shown to be a very effective technique for predicting barrier insertion loss in the absence of turbulence. This paper develops a simple and efficient modification of the BEM formulation to predict the insertion loss of a barrier in the presence of atmospheric turbulence. The modification is based on two alternative methods: (1) random realisations of log-amplitude and phase fluctuations of boundary sources and (2) de-correlation of source coherence using the mutual coherence function (MCF). An investigation into the behaviours of these two methods is carried out and simplified forms of the methods developed. Some systematic differences between the predictions from the methods are found. When incorporated into the BEM formulation, the method of random realisations and the method of MCF de-correlation provide predictions that agree well with predictions by the parabolic equation method and by the scattering cross-section method on a variety of thin barrier configurations.     

局部有限超拓扑的局部紧性

本文讨论了赋予局部有限拓扑的非空闲子集超空间的局部紧性.主要结果是:X正则,则其闭子集超空间局部紧当且仅当X可表示成一个紧空间与一个离散空间的拓扑和.     

Inner bounds for the spectrum of quasinormal operators

A linear operator in a separable Hilbert space is called a quasinormal one if it is a sum of a normal operator and a compact one. In the paper, bounds for the spectrum of quasinormal operators are established. In addition, the lower estimate for the spectral radius is derived. Under some restrictions, that estimate improves the well-known results. Applications to integral operators and matrices are discussed. Our results are new even in the finite-dimensional case.

     

用数值方法分析微孔激光器出射光的近场特性

微孔激光器作为应用于近场光信息存储系统中的一种新型光源，它的出射光斑的近场特性对于近场光存储是十分重要的。针对纳米孔径运用角谱进行Fox-Li数值迭代，得到不同孔径微孔激光器的基模光强分布，然后运用二维非线性时域有限差分法分析微孔激光器出射端即微孔金属膜的近场光学性质，模拟计算了不同孔径和厚度的微孔金属膜的光强近场分布，从应用于近场光存储的角度，给出反映其近场光学特性的相关数据。发现由于TM模式下金属存在局域表面等离子增强效应，使得其出射强度比TE模式高一个数量级，从而更适于作为实际中近场光存储系统和原理试验的光学头。