A new finite element realization of the perfectly matched layer method for Helmholtz scattering problems on polygonal domains in two dimensions

In this paper we propose a new finite element realization of the Perfectly Matched Layer method (PML-method). Our approach allows to deal with a wide class of polygonal domains and with certain types of inhomogeneous exterior domains. Among the covered inhomogeneities are open waveguide structures playing an essential role in integrated optics. We give a detailed insight into implementation aspects. Numerical examples show exponential convergence behavior to the exact solution with the thickness of the PML sponge layer.     

Minimal area problems for functions with integral representation

We study the minimization problem for the Dirichlet integral in some standard classes of analytic functions. In particular, we solve the minimal areaa ₂-problem for convex functions and for typically real functions. The latter gives a new solution to the minimal areaa ₂-problem for the classS of normalized univalent functions in the unit disc. Supported by NSF grant DMS-0412908.     

Well‐conditioned boundary integral formulations for high‐frequency elastic scattering problems in three dimensions

We construct and analyze a family of well‐conditioned boundary integral equations for the Krylov iterative solution of three‐dimensional elastic scattering problems by a bounded rigid obstacle. We develop a new potential theory using a rewriting of the Somigliana integral representation formula. From these results, we generalize to linear elasticity the well‐known Brakhage–Werner and combined field integral equation formulations. We use a suitable approximation of the Dirichlet‐to‐Neumann map as a regularizing operator in the proposed boundary integral equations. The construction of the approximate Dirichlet‐to‐Neumann map is inspired by the on‐surface radiation conditions method. We prove that the associated integral equations are uniquely solvable and possess very interesting spectral properties. Promising analytical and numerical investigations, in terms of spherical harmonics, with the elastic sphere are provided. Copyright © 2014 John Wiley & Sons, Ltd.     

Bang-bang property for Bolza problems in two dimensions

Consider the following Bolza problem: $$\begin{gathered} \min \int {h(x,u) dt,} \hfill \\ \dot x = F(x) + uG(x), \hfill \\ \left| u \right| \leqslant 1, x \in \Omega \subset \mathbb{R}^2 , \hfill \\ x(0) = x_0 , x(1) = x_1 . \hfill \\ \end{gathered} $$ We show that, under suitable assumptions onF, G, h, all optimal trajectories are bang-bang. The proof relies on a geometrical approach that works for every smooth two-dimensional manifold. As a corollary, we obtain existence results for nonconvex optimization problems.     

A volume integral equation method for periodic scattering problems for anisotropic Maxwell's equations

This paper presents a volume integral equation method for an electromagnetic scattering problem for three-dimensional Maxwell's equations in the presence of a biperiodic, anisotropic, and possibly discontinuous dielectric scatterer. Such scattering problem can be reformulated as a strongly singular volume integral equation (i.e., integral operators that fail to be weakly singular). In this paper, we firstly prove that the strongly singular volume integral equation satisfies a Gårding-type estimate in standard Sobolev spaces. Secondly, we rigorously analyze a spectral Galerkin method for solving the scattering problem. This method relies on the periodization technique of Gennadi Vainikko that allows us to efficiently evaluate the periodized integral operators on trigonometric polynomials using the fast Fourier transform (FFT). The main advantage of the method is its simple implementation that avoids for instance the need to compute quasiperiodic Green's functions. We prove that the numerical solution of the spectral Galerkin method applied to the periodized integral equation converges quasioptimally to the solution of the scattering problem. Some numerical examples are provided for examining the performance of the method.     

On free boundary problems in two dimensions

Let L be a linear elliptic operator in two dimensions with analytic coefficients and of second order, andu(x, y) a solution of Lu=0 in a simply connected domain ω with rectifiable boundary Γ. Suppose ψ(x, y) analytic on ω∪Γ and L ψ≠0 there.H is shown that ifu and ψ coincide with first derivatives on an open portion Γ₀ of Γ, then Γ₀ permits the representation λ=x (θ),y=y (θ) withx(θ),y(θ)analytic functions of a real parameter θ.     

A heuristic for solving large bin packing problems in two and three dimensions

The more-dimensional bin packing problem (BPP) considered here requires packing a set of rectangular-shaped items into a minimum number of identical rectangular-shaped bins. All items may be rotated and the guillotine cut constraint has to be respected. A straightforward heuristic is presented that is based on a method for the container loading problem following a wall-building approach and on a method for the one-dimensional BPP. 1,800 new benchmark instances are introduced for the two-dimensional and three-dimensional BPP. The instances include more than 1,500 items on average. Applied to these very large instances, the heuristic generates solutions of acceptable quality in short computation times. Moreover, the influence of different instance parameters on the solution quality is investigated by an extended computational study.     

A parallel method for time-discretization of parabolic problems based on contour integral representation and quadrature

We treat the time discretization of an initial-value problem for a homogeneous abstract parabolic equation by first using a representation of the solution as an integral along the boundary of a sector in the right half of the complex plane, then transforming this into a real integral on the finite interval , and finally applying a standard quadrature formula to this integral. The method requires the solution of a finite set of elliptic problems with complex coefficients, which are independent and may therefore be done in parallel. The method is combined with spatial discretization by finite elements.

     

An adaptive fast direct solver for boundary integral equations in two dimensions

We describe an algorithm for the rapid direct solution of linear algebraic systems arising from the discretization of boundary integral equations of potential theory in two dimensions. The algorithm is combined with a scheme that adaptively rearranges the parameterization of the boundary in order to minimize the ranks of the off-diagonal blocks in the discretized operator, thus obviating the need for the user to supply a parameterization r of the boundary for which the distance ‖r(s)−r(t)‖ between two points on the boundary is related to their corresponding distance |s−t| in the parameter space. The algorithm has an asymptotic complexity of , where N is the number of nodes in the discretization. The performance of the algorithm is illustrated with several numerical examples.     

Recurrent dimensions of quasi-periodic solutions for nonlinear evolution equations

In this paper we introduce recurrent dimensions of discrete dynamical systems and we give upper and lower bounds of the recurrent dimensions of the quasi-periodic orbits. We show that these bounds have different values according to the algebraic properties of the frequency and we investigate these dimensions of quasi-periodic trajectories given by solutions of a nonlinear PDE.

     

Exponentially fitted shape functions for advection-dominated flow problems in two dimensions

We sketch out the basic properties of three novel exponentially fitted shape functions that generalize to two-dimensional triangular elements the one-dimensional functions used in the well-known Scharfetter-Gummel method. This scheme is widely employed in the finite element approximation of the drift-diffusion equations arising in semiconductor device modelling.     

关于代数多项式的一个积分表示式

研究了代数多项式导数的Bernstein不等式和Markov不等式．通过代数多项式导数的一个积分表示式,给出这两个著名不等式以及它们的离散形式的证明．     

Volume singular integral equations for problems of scattering on three-dimensional dielectric structures

We consider volume singular integral equations describing problems of scattering of electromagnetic waves in bounded three-dimensional dielectric structures. We analyze the equations mathematically. Existence and uniqueness theorems are presented, and the spectrum of the integral operators is studied.     

A new efficient DQ algorithm for the solution of elliptic problems in higher dimensions

Two new efficient algorithms are developed to approximate the derivatives of sufficiently smooth functions. The new techniques are based on differential quadrature method with quartic B-spline bases as test functions. To obtain the weighting coefficients of differential quadrature method (DQM), we use the midpoints of a uniform partition mixed with near-boundary grid points. This enables us to obtain the weighting coefficients without adding the new extra relations. By obtaining the error bounds, it is proved that the method in its classic form is non-optimal. Then, some new weighting coefficients are constructed to obtain higher accuracy. By obtaining the error bounds, it is proved that the new algorithm is superconvergent. Afterwards, by defining some new symbols, we find a way to approximate the partial derivatives of multivariate functions. Also, some approximations are constructed to the mixed derivatives of multivariate functions. Finally, the applicability of the methods is examined by solving some well-known problems of partial differential equations. Some examples of 2D and 3D biharmonic, Poisson, and convection-diffusion equations are solved and compared to the existing methods to show the efficiency of the proposed algorithms.     

Error estimates for a discretized Galerkin method for a boundary integral equation in two dimensions

We present a priori and a posteriori estimates for the error between the Galerkin and a discretized Galerkin method for the boundary integral equation for the single layer potential on the square plate. Using piecewise constant finite elements on a rectangular mesh we study the error coming from numerical integration. The crucial point of our analysis is the estimation of some error constants, and we demonstrate that this is necessary if our methods are to be used. After the determination of these constants we are in the position to prove invertibility and quasioptimal convergence results for our numerical scheme, if the chosen numerical integration formulas are sufficiently precise. © 1992 John Wiley & Sons, Inc.