An a priori error estimate for the finite element modelling of electromagnetic waves interacting with a periodic diffraction grating

An a priori error estimate using a so called α,β‐ periodic transformation to study electromagnetic waves in a periodic diffraction grating is derived. It has been reported for single scattering that there is an instability in numerical methods for high wavenumbers. To address this problem, the analytical solution of the scattering problem when the domain is scatterer free and an unknown function called the α,β‐quasi periodic solution are used to transform the associated Helmholtz problem. The well‐posedness of the resulting continuous problem is analysed before approximating its solution using a finite element discretization. To guarantee the uniqueness of this approximate solution, an a priori error estimate is derived. Finally, numerical results are presented that suggest that the α,β‐quasi periodic method converges at a far lower number of degrees of freedom than the α,0‐quasi periodic method reported previously; especially for high wavenumbers. This is particularly true when the incident wave only undergoes a small perturbation because of the presence of the scatterer. Copyright © 2012 John Wiley & Sons, Ltd.     

Existence Uniqueness and Decay of Solution for Fractional Boussinesq Approximation

The Boussinesq approximation finds more and more frequent use in geological practice. In this paper, the asymptotic behavior of solution for fractional Boussinesq approximation is studied. After obtaining some a priori estimates with the aid of commutator estimate, we apply the Galerkin method to prove the existence of weak solution in the case of periodic domain. Meanwhile, the uniqueness is also obtained. Because the results obtained are independent of domain, the existence and uniqueness of the weak solution for Cauchy problem is also true. Finally, we use the Fourier splitting method to prove the decay of weak solution in three cases respectively.     

A two-sided Arnoldi algorithm with stopping criterion and MIMO selection procedure

In this paper we introduce a two-sided Arnoldi method for the reduction of high order linear systems and we propose useful extensions, first of all a stopping criterion to find a suitable order for the reduced model and secondly, a selection procedure to significantly improve the performance in the multi-input multi-output (MIMO) case. One application is in micro-electro-mechanical systems (MEMS). We consider a thermo-electric micro thruster model, and a comparison between the commonly used Arnoldi algorithm and the two-sided Arnoldi is performed.     

关于连续素数上取值的Legendre符号和的估计

利用Karatsuba的方法,研究了关于连续素数上取值的Legendre符号和∑p≤N(p q)的估计问题,得到了当q0.75+ω相似文献   

Estimating complex covariance by observing two variables at a time

The estimation of covariance matrices is central in array signal processing systems. This note addresses complex covariance estimation for the situation, where the complex data are available only as independent pairwise sets (observations) corresponding to individual elements of the matrix. The formulation for the empirical estimate and the normal maximum likelihood estimate is developed for the general case of different sample sizes for each observation. The approach allows, for example, the estimate of the p by p covariance matrix of a p-port sensor array from a two-port measurement instrument.     

关于“序列号”方法的讨论

对于无放回抽样和有放回抽样两种情形,研究了"序列号"方法,给出并讨论了某些相关的估计量的优良性.     

Hardy spaces H_L~p(R~n) associated with operators satisfying k-Davies-Gaffney estimates

Let L be a one-to-one operator of type ω having a bounded H∞ functional calculus and satisfying the k-Davies-Gaffney estimates with k ∈ N. In this paper, the authors introduce the Hardy space HLp(Rn) with p ∈ (0, 1] associated with L in terms of square functions defined via {e-t2kL}t>0 and establish their molecular and generalized square function characterizations. Typical examples of such operators include the 2k-order divergence form homogeneous elliptic operator L1 with complex bounded measurable coefficients and the 2k-order Schrdinger type operator L2 := (-Δ)k + Vk, where Δ is the Laplacian and 0≤V∈Llock(Rn). Moreover, as an application, for i ∈ {1, 2}, the authors prove that the associated Riesz transform ▽k(Li-1/2) is bounded from HLip (Rn) to Hp(Rn) for p ∈ (n/(n + k), 1] and establish the Riesz transform characterizations of HL1p (Rn) for p ∈ (rn/(n + kr), 1] if {e-tL1 }t>0 satisfies the Lr-L2 k-off-diagonal estimates with r ∈ (1, 2]. These results when k := 1 and L := L1 are known.     

A finite difference scheme for two-dimensional semiconductor device of heat conduction on composite triangular grids

The momentary state of a semiconductor device of heat conduction is described by a system of four nonlinear partial differential equations. One elliptic equation is for the electrostatic, two parabolic equations are for the electron concentration and the hole concentration, and one heat exchange equation is for the temperature. According to the necessary of practical numerical simulations and based on the balance equation, finite difference schemes for two-dimensional transient behavior of a semiconductor device of heat conduction on composite triangular grids are constructed. Studying their stability and convergence properties, the error estimate in the energy norm is obtained. Finally, a numerical example is given.     

Natural hp‐BEM for the electric field integral equation with singular solutions

We apply the hp ‐version of the boundary element method (BEM) for the numerical solution of the electric field integral equation (EFIE) on a Lipschitz polyhedral surface Γ. The underlying meshes are supposed to be quasi‐uniform triangulations of Γ, and the approximations are based on either Raviart‐Thomas or Brezzi‐Douglas‐Marini families of surface elements. Nonsmoothness of Γ leads to singularities in the solution of the EFIE, severely affecting convergence rates of the BEM. However, the singular behavior of the solution can be explicitly specified using a finite set of functions (vertex‐, edge‐, and vertex‐edge singularities), which are the products of power functions and poly‐logarithmic terms. In this article, we use this fact to perform an a priori error analysis of the hp ‐BEM on quasi‐uniform meshes. We prove precise error estimates in terms of the polynomial degree p, the mesh size h, and the singularity exponents. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

Adaptive meshless Galerkin boundary node methods for hypersingular integral equations

Adaptive refinement techniques are developed in this paper for the meshless Galerkin boundary node method for hypersingular boundary integral equations. Two types of error estimators are derived. One is a perturbation error estimator that is formulated based on the difference between numerical solutions obtained using two consecutive nodal arrangements. The other is a projection error estimator that is formulated based on the difference between the numerical solution itself and its projection. These error estimators are proven to have an upper and a lower bound by the constant multiples of the exact error in the energy norm. A localization scheme is presented to accomodate the non-local property of hypersingular integral operators for the needed computable local error indicators. The convergence of the adaptive meshless techniques is verified theoretically. To confirm the theoretical results and to show the efficiency of the adaptive techniques, numerical examples in 2D and 3D with high singularities are provided.