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.     

Solution of the hyperbolic mild‐slope equation using the finite volume method

A finite volume solver for the 2D depth‐integrated harmonic hyperbolic formulation of the mild‐slope equation for wave propagation is presented and discussed. The solver is implemented on unstructured triangular meshes and the solution methodology is based upon a Godunov‐type second‐order finite volume scheme, whereby the numerical fluxes are computed using Roe's flux function. The eigensystem of the mild‐slope equations is derived and used for the construction of Roe's matrix. A formulation that updates the unknown variables in time implicitly is presented, which produces a more accurate and reliable scheme than hitherto available. Boundary conditions for different types of boundaries are also derived. The agreement of the computed results with analytical results for a range of wave propagation/transformation problems is very good, and the model is found to be virtually paraxiality‐free. Copyright © 2003 John Wiley & Sons, Ltd.     

一类雨量预报准确性评价方法

给出了一种利用网格点的预报雨量来模拟区域内任一点预报雨量的方法,建立了预报雨量函数,进而得出预报雨量与实测雨量的整体偏方差的数学模型;利用QPF方法建立出准确率、预报效率、空报率、漏报率等评价预报方法优劣指标的数学模型,并对两种预报方法得出的数据进行了计算,得出上述模型中的数值如下:方案方法总体偏方差预报准确率漏报率空报率预报效率方法1 0.99987×10679.431%11.463%25.621%62.256%方法2 1.00518×10657.509%26.976%31.705%58.936%比较这些指标,可得第一种预报结果较好.     

ROBUSTNESS OF AN UPWIND FINITE DIFFERENCE SCHEME FOR SEMILINEAR CONVECTION-DIFFUSION PROBLEMS WITH BOUNDARY TURNING POINTS

We consider a singularly perturbed semilinear convection-diffusion problem with a boundary layer of attractive turning-point type. It is shown that its solution can be decomposed into a regular solution component and a layer component. This decomposi-tion is used to analyse the convergence of an upwinded finite difference scheme on Shishkin meshes.     

Analysis of a finite-difference scheme for a singularly perturbed problem with two small parameters

We study a model linear convection-diffusion-reaction problem where both the diffusion term and the convection term are multiplied by small parameters ε_d and ε_c, respectively. Depending on the size of the parameters the solution of the problem may exhibit exponential layers at both end points of the domain. Sharp bounds for the derivatives of the solution are derived using a barrier-function technique. These bounds are applied in the analysis of a simple upwind-difference scheme on Shishkin meshes. This method is established to be almost first-order convergent, independently of the parameters ε_d and ε_c.     

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.     

THE GENERALIZED MAXIMUM ANGLE CONDITION FOR THE Q1 ISOPARAMETRIC ELEMENT

We consider the quadrilateral Q1 isoparametric element and establish an optimal error estimate in H^1 norm for the interpolation operator under a weaker mesh condition which admits anisotropic quadrilaterals and allows the quadrilateral to become a regular triangle in the sense of maximum angle condition [5, 11].     

曲率流的参数化有限元逼近

许多物理现象可以在数学上描述为受曲率驱动的自由界面运动,例如薄膜和泡沫的演变、晶体生长,等等.这些薄膜和界面的运动常依赖于其表面曲率,从而可以用相应的曲率流来描述,其相关自由界面问题的数值计算和误差分析一直是计算数学领域中的难点.参数化有限元法是曲率流的一类有效计算方法,已经能够成功模拟一些曲面在几类基本的曲率流下的演化过程.本文重点讨论曲率流的参数化有限元逼近,它的产生、发展和当前的一些挑战.     

曲边区域上的多边形网格间断有限元离散及其多重网格算法

本文利用多边形网格上的间断有限元方法离散二阶椭圆方程,在曲边区域上,采用多条直短边逼近曲边的以直代曲的策略,实现了高阶元在能量范数下的最优收敛.本文还将这一方法用于带曲边界面问题的求解,同样得到高阶元的最优收敛.此外我们还设计并分析了这一方法的\linebreakW-cycle和Variable V-cycle多重网格预条件方法,证明当光滑次数足够多时,多重网格预条件算法一致收敛.最后给出了数值算例,证实该算法的可行性并验证了理论分析的结果.