Transmission problems for Maxwell's equations with weakly Lipschitz interfaces

We prove sufficient conditions on material constants, frequency and Lipschitz regularity of interface for well posedness of a generalized Maxwell transmission problem in finite energy norms. This is done by embedding Maxwell's equations in an elliptic Dirac equation, by constructing the natural trace space for the transmission problem and using Hodge decompositions for operators d and δ on weakly Lipschitz domains to prove stability. We also obtain results for boundary value problems and transmission problems for the Hodge–Dirac equation and prove spectral estimates for boundary singular integral operators related to double layer potentials. Copyright © 2005 John Wiley & Sons, Ltd.     

Convergence analysis of a hp-finite element approximation of the time-harmonic Maxwell equations with impedance boundary conditions in domains with an analytic boundary

We consider a nonconforming hp -finite element approximation of a variational formulation of the time-harmonic Maxwell equations with impedance boundary conditions proposed by Costabel et al. The advantages of this formulation is that the variational space is embedded in H¹ as soon as the boundary is smooth enough (in particular it holds for domains with an analytic boundary) and standard shift theorem can be applied since the associated boundary value problem is elliptic. Finally in order to perform a wavenumber explicit error analysis of our problem, a splitting lemma and an estimation of the adjoint approximation quantity are proved by adapting to our system the results from Melenk and Sauter obtained for the Helmholtz equation. Some numerical tests that illustrate our theoretical results are also presented. Analytic regularity results with bounds explicit in the wavenumber of the solution of a general elliptic system with lower order terms depending on the wavenumber are needed and hence proved.     

Two-sided a posteriori error bounds for electro-magnetostatic problems

This paper is concerned with the derivation of computable and guaranteed upper and lower bounds of the difference between exact and approximate solutions of a boundary value problem for static Maxwell equations. Our analysis is based upon purely functional argumentation and does not invoke specific properties of the approximation method. For this reason, the estimates derived in the paper at hand are applicable to any approximate solution that belongs to the corresponding energy space. Such estimates (also called error majorants of the functional type) have been derived earlier for elliptic problems. Bibliography: 24 titles.     

线性椭圆问题最低次混合有限元方法的最优最大模误差估计

本文对线性椭圆问题的最低次混合元方法提出了构造混合元空间的充分条件,并建立了新的插值算子.据此得到了混合元解,伴随向量函数及其散度的最优最大模误差估计.     

Graded mesh refinement and error estimates of higher order for DGFE solutions of elliptic boundary value problems in polygons

Error estimates for DGFE solutions are well investigated if one assumes that the exact solution is sufficiently regular. In this article, we consider a Dirichlet and a mixed boundary value problem for a linear elliptic equation in a polygon. It is well known that the first derivatives of the solutions develop singularities near reentrant corner points or points where the boundary conditions change. On the basis of the regularity results formulated in Sobolev–Slobodetskii spaces and weighted spaces of Kondratiev type, we prove error estimates of higher order for DGFE solutions using a suitable graded mesh refinement near boundary singular points. The main tools are as follows: regularity investigation for the exact solution relying on general results for elliptic boundary value problems, error analysis for the interpolation in Sobolev–Slobodetskii spaces, and error estimates for DGFE solutions on special graded refined meshes combined with estimates in weighted Sobolev spaces. Our main result is that there exist a local grading of the mesh and a piecewise interpolation by polynoms of higher degree such that we will get the same order O (h^α) of approximation as in the smooth case. © 2011 Wiley Periodicals, Inc. Numer Mehods Partial Differential Eq, 2012     

Finite Element Approximation of Elliptic Dirichlet Optimal Control Problems

We develop a priori error analysis for the finite element Galerkin discretization of elliptic Dirichlet optimal control problems. The state equation is given by an elliptic partial differential equation and the finite dimensional control variable enters the Dirichlet boundary conditions. We prove the optimal order of convergence and present a numerical example confirming our results.     

阻尼Sine-Gordon方程的H1-Galerkin混合元方法数值解

利用H1-Galerkin混合有限元方法讨论阻尼Sine-Gordon方程,得到一维情况下半离散和全离散格式的最优阶误差估计,并且推广应用到二维和三维情况,而且不用验证LBB相容性条件.     

伪双曲型积分-微分方程的H~1-Galerkin混合元法误差估计

<正>1引言考虑如下一类具有Lipschitz连续边界(?)Ω的凸有界区域Ω上的伪双曲型积分微分方程其中Ω(?)R~d,(d=1,2,3)J=(0,T],对于固定的T,0T∞,函数0a_0≤a(x,t)≤     

A plasmastatic model of the galathea-belt magnetic trap

The mathematical apparatus of plasmastatics, which includes the MHD equilibrium equation and steady-state Maxwell equations, is reduced, in two-dimensional problems arising due to symmetry, to a single scalar second-order elliptic equation with a nonlinear right-hand side known as the Grad-Shafranov equation. In this paper, we numerically solve a series of boundary value problems for this equation that model equilibrium plasma configurations in the magnetic field of the belt-like galathea trap in a cylinder with two plasma embedded conductors. The mathematical model is outlined, the results of calculations of the magnetic field and plasma pressure in the cylinder depending on the parameters of the problem are presented, and the main integral characteristics of the trap are calculated. The existence and uniqueness of the solution is discussed; the limiting values of the maximal pressure at which there exists a solution of the equilibrium problem are found.     

HILBERT SPACES ASSOCIATED WITH SECOND ORDER ELLIPTIC OPERATORS AND APPLICATIONS TO BOUNDARY VALUE PROBLEMS AND VARIATIONAL INEQUALITIES

Abstract

Boundary value problems and variational inequalities, associated with second order elliptic operators, will be studied in a Hilbert space framework. In this space, functions will have (at least) locally square integrable derivatives of order up to two. Also the conormal derivative, extended by continuity, will be square integrable on the boundary of the region considered. Criteria for approximating elements of the Hilbert space by smooth functions will be given and thus closed convex sets, associated with inequalities on the boundary, exist.

The idea of the present approach originated from the method suggested by Lions and Magenes, for putting some regular elliptic problems in the variational setting. The differential equation is multiplied by Qv, with Q some operator and v a function and the result is integrated as required.     

Error bounds for finite element solutions of mildly nonlinear elliptic boundary value problems

Summary The equivalence in a Hilbert space of variational and weak formulations of linear elliptic boundary value problems is well known. This same equivalence is proved here for mildly nonlinear problems where the right hand side of the differential equation involves the solution function. A finite element approximation to the solution of the weak problem ina finite dimensional subspace of the original Hilbert space is defined. An inequality bounding the error in this approximation over all functions of the space is derived, and in particular this holds for an interpolant to the weak solution. Thus this inequality, together with previously known, interpolation error bounds, produces a bound on the finite element solution to this nonlinear problem. An example of a mildly nonlinear Poisson problem is given.     

Multiscale asymptotic expansion and finite element methods for the mixed boundary value problems of second order elliptic equation in perforated domains

In this paper, we will discuss the mixed boundary value problems for the second order elliptic equation with rapidly oscillating coefficients in perforated domains, and will present the higher-order multiscale asymptotic expansion of the solution for the problem, which will play an important role in the numerical computation . The convergence theorems and their rigorous proofs will be given. Finally a multiscale finite element method and some numerical results will be presented. This work is Supported by National Natural Science Foundation of China (grant # 10372108, # 90405016), and Special Funds for Major State Basic Research Projects( grant # TG2000067102)     

Error estimates of H^1-Galerkin mixed finite element method for Schrodinger equation

An H1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.     

A general mixed covolume framework for constructing conservative schemes for elliptic problems

We present a general framework for the finite volume or covolume schemes developed for second order elliptic problems in mixed form, i.e., written as first order systems. We connect these schemes to standard mixed finite element methods via a one-to-one transfer operator between trial and test spaces. In the nonsymmetric case (convection-diffusion equation) we show one-half order convergence rate for the flux variable which is approximated either by the lowest order Raviart-Thomas space or by its image in the space of discontinuous piecewise constants. In the symmetric case (diffusion equation) a first order convergence rate is obtained for both the state variable (e.g., concentration) and its flux. Numerical experiments are included.

     

On the absence of eigenvalues of Maxwell and Lamé systems in exterior domains

The arguments showing non‐existence of eigensolutions to exterior‐boundary value problems associated with systems—such as the Maxwell and Lamé system—rely on showing that such solutions would have to have compact support and therefore—by a unique continuation property—cannot be non‐trivial. Here we will focus on the first part of the argument. For a class of second order elliptic systems it will be shown that L²‐solutions in exterior domains must have compact support. Both the asymptotically isotropic Maxwell system and the Lamé system with asymptotically decaying perturbations can be reduced to this class of elliptic systems. Copyright © 2004 John Wiley & Sons, Ltd.     

A power penalty approach to a mixed quasilinear elliptic complementarity problem

In this paper, a power penalty approximation method is proposed for solving a mixed quasilinear elliptic complementarity problem. The mixed complementarity problem is first reformulated as a double obstacle quasilinear elliptic variational inequality problem. A nonlinear elliptic partial differential equation is then defined to approximate the resulting variational inequality by using a power penalty approach. The existence and uniqueness of the solution to the partial differential penalty equation are proved. It is shown that, under some mild assumptions, the sequence of solutions to the penalty equations converges to the unique solution of the variational inequality problem as the penalty parameter tends to infinity. The error estimates of the convergence of this penalty approach are also derived. At last, numerical experimental results are presented to show that the power penalty approximation method is efficient and robust.

     

A POSTERIORI ERROR ESTIMATES OF FINITE ELEMENT METHOD FOR PARABOLIC PROBLEMS

1IntroductionThebase0fadaPtivecomputing0ffiniteelementmethodisap0steri0rierr0restimates.I.Babuskaisthepioneerinthisfields.Manytechniquesaredevel0pedtoobtainaposteri0rierrorestimators.See[1-3,7-8,19-201.Theyaremainlybased0nthejumps0fthederiva-tivesontheboundariesoftl1eelel11elltandtheresidualintheelemellts.Recelltresultssh0wthatthereareveryclosedrelatiollsbetweellasymptoticexactap0steri0rierrorestimatesandsuperc0nvergence-SeealsoQ.Linetal.[11-13],andChen-Huang['].Therehasbeenmuchprogressill…     

A Fast Algorithm for Two-Dimensional Elliptic Problems

In this paper, we extend the work of Daripa et al. [14–16,7] to a larger class of elliptic problems in a variety of domains. In particular, analysis-based fast algorithms to solve inhomogeneous elliptic equations of three different types in three different two-dimensional domains are derived. Dirichlet, Neumann and mixed boundary value problems are treated in all these cases. Three different domains considered are: (i) interior of a circle, (ii) exterior of a circle, and (iii) circular annulus. Three different types of elliptic problems considered are: (i) Poisson equation, (ii) Helmholtz equation (oscillatory case), and (iii) Helmholtz equation (monotone case). These algorithms are derived from an exact formula for the solution of a large class of elliptic equations (where the coefficients of the equation do not depend on the polar angle when written in polar coordinates) based on Fourier series expansion and a one-dimensional ordinary differential equation. The performance of these algorithms is illustrated for several of these problems. Numerical results are presented.     

On the use of characteristic variables in viscoelastic flow problems

The system of partial differential equations governing the flow of an upper converted Maxwell fluid is known to be of mixed elliptic–hyperbolic type. The hyperbolic nature of the constitutive equation requires that, where appropriate, inflow conditions are prescribed in order to obtain a well-posed problem. Although there are three convective derivatives in the constitutive equation there are only two characteristic quantities whichare transported along the streamlines. These characteristicquantities are identified. A spectral element method is describedin which continuity of the characteristic variables is usedto couple the extra stress components between contiguous elements.The continuity of the characteristic variables is treated asa constraint on the constitutive equation. These conditionsdo not necessarily impose continuity on the extra-stress components.The velocity and pressure follow from the doubly constrainedweak formulation which enforces a divergence-free velocityfield and irrotational polymeric stress forces. This meansthat both the pressure and the extra-stress tensor are discontinuous.Numerical results are presented to demonstrate this procedure.The theory is applied to the upper convected Maxwell modelwith vanishing Reynolds number. No regularization techniques such as streamline upwind Petrov Galerkin (SUPG), elastic viscous split stress (EVSS) or explicitly elliptic momentum equation(EEME) are used.     

Stability of the space identification problem for the elliptic‐telegraph differential equation

The present paper is devoted to study the space identification problem for the elliptic‐telegraph differential equation in Hilbert spaces with the self‐adjoint positive definite operator. The main theorem on the stability of the space identification problem for the elliptic‐telegraph differential equation is proved. In applications, theorems on the stability of three source identification problems for one dimensional with nonlocal conditions and multidimensional elliptic‐telegraph differential equations are established.