The Two-dimensional Electromagnetic Scattering from Periodic Chiral Structures and Its Finite Element Approximation

In this paper, we consider the electromagnetic scattering from periodic chiral structures. The structure is periodic in one direction and invariant in another direction. The electromagnetic fields in the chiral medium are governed by the Maxwell equations together with the Drude-Born-Fedorov equations. We simplify the problem to a two-dimensional scattering problem and we show that for all but possibly a discrete set of wave numbers, there is a unique quasi-periodic weak solution to the diffraction problem. The diffraction problem can be solved by finite element method. We also establish uniform error estimates for the finite element method and the error estimates when the truncation of the nonlocal transparent boundary operators takes place.     

Numerical analysis of diffraction by periodic structures: TM polarization

Summary. The diffraction of a time harmonic wave in TM polarization by some periodic inhomogeneous material is studied in this paper. The diffraction problem is modeled by a generalized Helmholtz equation with transparent boundary conditions. Existence and uniqueness are proved for the model problem. Uniform error estimates for the finite element approximation of the solution are established. Error estimates are also obtained when the truncation of the nonlocal boundary operators takes place. Received July 24, 1994 / Revised version received August 24, 1995     

Error analysis of the DtN-FEM for the scattering problem in acoustics via Fourier analysis

In this paper, we are concerned with the error analysis for the finite element solution of the two-dimensional exterior Neumann boundary value problem in acoustics. In particular, we establish explicit priori error estimates in H¹ and L²- norms including both the effect of the truncation of the DtN mapping and that of the numerical discretization. To apply the finite element method (FEM) to the exterior problem, the original boundary value problem is reduced to an equivalent nonlocal boundary value problem via a Dirichlet-to-Neumann (DtN) mapping represented in terms of the Fourier expansion series. We discuss essential features of the corresponding variational equation and its modification due to the truncation of the DtN mapping in appropriate function spaces. Numerical tests are presented to validate our theoretical results.     

Discontinuous Galerkin method for nonlinear convection-diffusion problems with mixed Dirichlet-Neumann boundary conditions

The paper is devoted to the analysis of the discontinuous Galerkin finite element method (DGFEM) applied to the space semidiscretization of a nonlinear nonstationary convection-diffusion problem with mixed Dirichlet-Neumann boundary conditions. General nonconforming meshes are used and the NIPG, IIPG and SIPG versions of the discretization of diffusion terms are considered. The main attention is paid to the impact of the Neumann boundary condition prescribed on a part of the boundary on the truncation error in the approximation of the nonlinear convective terms. The estimate of this error allows to analyse the error estimate of the method. The results obtained represent the completion and extension of the analysis from V. Dolej?í, M. Feistauer, Numer. Funct. Anal. Optim. 26 (2005), 349–383, where the truncation error in the approximation of the nonlinear convection terms was proved only in the case when the Dirichlet boundary condition on the whole boundary of the computational domain was considered.     

Error estimates of the finite element method for the exterior Helmholtz problem with a modified DtN boundary condition

A priori error estimates in the H¹- and L²-norms are established for the finite element method applied to the exterior Helmholtz problem, with modified Dirichlet-to-Neumann (MDtN) boundary condition. The error estimates include the effect of truncation of the MDtN boundary condition as well as that of discretization of the finite element method. The error estimate in the L²-norm is sharper than that obtained by the author [D. Koyama, Error estimates of the DtN finite element method for the exterior Helmholtz problem, J. Comput. Appl. Math. 200 (1) (2007) 21-31] for the truncated DtN boundary condition.     

Nonlinear fourth-order elliptic equations with nonlocal boundary conditions

This paper is concerned with a class of fourth-order nonlinear elliptic equations with nonlocal boundary conditions, including a multi-point boundary condition in a bounded domain of Rn. Also considered is a second-order elliptic equation with nonlocal boundary condition, and the usual multi-point boundary problem in ordinary differential equations. The aim of the paper is to show the existence of maximal and minimal solutions, the uniqueness of a positive solution, and the method of construction for these solutions. Our approach to the above problems is by the method of upper and lower solutions and its associated monotone iterations. The monotone iterative schemes can be developed into computational algorithms for numerical solutions of the problem by either the finite difference method or the finite element method.     

On the h-adaptive coupling of FE and BE for viscoplastic and elasto-plastic interface problems

This paper presents a posteriori error estimates for the symmetric finite element and boundary element coupling for a nonlinear interface problem: A bounded body with a viscoplastic or plastic material behaviour is surrounded by an elastic body. The nonlinearity is treated by the finite element method while large parts of the linear elastic body are approximated using the boundary element method. Based on the a posteriori error estimates we derive an algorithm for the adaptive mesh refinement of the boundary elements and the finite elements. Its implementation is documented and numerical examples are included.     

Error estimates of the DtN finite element method for the exterior Helmholtz problem

A priori error estimates are established for the DtN (Dirichlet-to-Neumann) finite element method applied to the exterior Helmholtz problem. The error estimates include the effect of truncation of the DtN boundary condition as well as that of the finite element discretization. A property of the Hankel functions which plays an important role in the proof of the error estimates is introduced.     

A problem with wells for the steady diffusion equation

This paper deals with a problem with wells for which nonlocal boundary conditions are given. It is shown that the problem is equivalent to a mixed problem without wells. For this formulation, an error estimate of a mixed finite element method in the 2D case is studied.     

A posteriori error analysis for linearization of nonlinear elliptic problems and their discretizations

The paper is devoted to a posteriori quantitative analysis for errors caused by linearization of non-linear elliptic boundary value problems and their finite element realizations. We employ duality theory in convex analysis to derive computable bounds on the difference between the solution of a non-linear problem and the solution of the linearized problem, by using the solution of the linearized problem only. We also derive computable bounds on differences between finite element solutions of the nonlinear problem and finite element solutions of the linearized problem, by using finite element solutions of the linearized problem only. Numerical experiments show that our a posteriori error bounds are efficient.     

Finite Difference Methods for the Heat Equation with a Nonlocal Boundary Condition

We consider the numerical solution by finite difference methods of the heat equation in one space dimension, with a nonlocal integral boundary condition, resulting from the truncation to a finite interval of the problem on a semi-infinite interval. We first analyze the forward Euler method, and then the $θ$-method for $0 < θ ≤ 1$, in both cases in maximum-norm, showing $O(h^2 + k)$ error bounds, where $h$ is the mesh-width and $k$ the time step. We then give an alternative analysis for the case $θ = 1/2$, the Crank-Nicolson method, using energy arguments, yielding a $O(h^2$ + $k^{3/2}$) error bound. Special attention is given the approximation of the boundary integral operator. Our results are illustrated by numerical examples.     

Finite element approximation of nonlocal parabolic problem

In this article, we present a finite element scheme combined with backward Euler method to solve a nonlocal parabolic problem. An important issue in the numerical solution of nonlocal problems while using Newton's method is related to its structure. In fact differently from the local case where the Jacobian matrix is sparse and banded, in the nonlocal case the Jacobian matrix is dense and computations are much more onerous compared to that for differential equations. In order to avoid this difficulty, we use the technique given by Gudi (SIAM J Numer Anal 50 (2012), 657–668) for elliptic nonlocal problem of Kirchhoff type. We discuss the well‐posedness of the weak formulation at continuous as well as at discrete levels. We also derive a priori error estimates for semidiscrete and fully discrete formulations in L² and H¹ norms. Results based on the usual finite element method are provided to confirm the theoretical estimates. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 786–813, 2017     

一个扩散问题的自然边界元法与有限元法组合

本文讨论由Helmholtz方程描述的扩散问题的自然边界元法与有限元法的组合．取一个圆作为公共边界，用Fourier展开建立边界积分方程，将无界区域上的问题化为有界区域上的非局部边值问题．在变分方程中公共边界上的未知量只包含函数本身而不包含其法向导数，从而减少了未知数的数目，并且边界元剐度矩阵只有极少量不同的元素，有利于数值计算．这种组台方法优越于建立在直接边界元法基础上的组合方法．文中证明了变分解的唯一性，数值解的收敛性和误差估计．最后讨论了数值技术并给出一个算倒．     

A Posteriori Error Estimates for Axisymmetric and Nonlinear Problems

We propose and examine the primal and dual finite element method for solving an axially symmetric elliptic problem with mixed boundary conditions. We derive an a posteriori error estimate and generalize the method used for a nonlinear elliptic problem. Finally, an a posteriori error estimate for a nonlinear parabolic problem based on the concept of hierarchical finite element basis functions is introduced.     

三维热传导型半导体问题的特征混合元方法和分析

本文研究三维热传导型半导体态问题的特征混合元方法及其理论分析,其数学模型是一类非线性偏微分方程的初边值问题,对电子位势方程提出混合元逼近,对电子,空穴浓度方程笔挺表限元逼近;对热传导方程采用对时间向后差分的Ｇａｌｅｒｋｉｎ逼近,应用微分方程先验估计理论和技巧得到了最优阶Ｌ＾２误差估计。     

Spectral analysis of a high-order Hermitian algorithm for structural dynamics

The spectral analysis of an efficient step-by-step direct integration algorithm for the structural dynamic equation is presented. The proposed algorithm is formulated in terms of two Hermitian finite difference operators of fifth-order local truncation error and it is unconditionally stable with no numerical damping presenting a fourth-order truncation error for period dispersion (global error). In addition, although it is in competition with higher-order algorithms presented in the literature, the computational effort is similar to that of the classical second-order Newmark’s method. The numerical application for nonlinear structural dynamic problems is also considered.     

The coupling of boundary integral and finite element methods for nonmonotone nonlinear problems ∗

In this paper, we apply the coupling of the boundary integral and finite element methods to study the weak solvability of certain nonmonotone nonlinear exterior boundary value problems. In order to convert the original exterior problem into an equivalent nonlocal boundary value problem on a finite region, we employ two different approaches based on the use of one and two integral equations on the coupling boundary. Existence of a solution for the associated weak formulation, and convergence properties of the corresponding Galerkin approximations are deduced from fundamental results in nonlinear functional analysis. Indeed, the main arguments of our proofs are based on a surjectivity theorem for mappings of type (S) and on the Fredholm alternative for nonlinear A-proper mappings.     

三维热传导型半导体问题的交替方向特征有限元方法及理论分析

１．引言 三维热传导型半导体器件瞬态问题的数学模型由四个非线性偏微分方程描述 ［１,２］．工程研究中一般考虑绝流边条件,由于绝流条件可以看作一反射条件来处理、为了数值分析方便,我们在此考虑三维周期问题： 其中,　＝［０,１］３,未知函数是电子位势　;电子,空穴浓度ｅ,ｐ;温度函数Ｔ．方程（１,１）－（１.４）中出现的系数均有正的上下界,且是　周期的． ａ＝Ｑ／ε,Ｑ,ε分别表示电子负荷和介电系数,均为正常数．Ｎ（ｘ）是给定的函数．Ｄｓ（ｘ）为扩散系数,μｓ（ｘ）为迁移率,ｓ＝ｅ,Ｐ．Ｒ（ｅ,ｐ,Ｔ）…     

Error estimates for the finite element approximation of linear elastic equations in an unbounded domain

In this paper we present error estimates for the finite element approximation of linear elastic equations in an unbounded domain. The finite element approximation is formulated on a bounded computational domain using a nonlocal approximate artificial boundary condition or a local one. In fact there are a family of nonlocal approximate boundary conditions with increasing accuracy (and computational cost) and a family of local ones for a given artificial boundary. Our error estimates show how the errors of the finite element approximations depend on the mesh size, the terms used in the approximate artificial boundary condition, and the location of the artificial boundary. A numerical example for Navier equations outside a circle in the plane is presented. Numerical results demonstrate the performance of our error estimates.

     

Interface problem in holonomic elastoplasticity

The three-dimensional interface problem with the homogeneous Lamé system in an unbounded exterior domain and holonomic material behaviour in a bounded interior Lipschitz domain is considered. Existence and uniqueness of solutions of the interface problem are obtained rewriting the exterior problem in terms of boundary integral operators following the symmetric coupling procedure. The numerical approximation of the solutions consists in coupling of the boundary element method (BEM) and the finite element method (FEM). A Céa-like error estimate is presented for the discrete solutions of the numerical procedure proving its convergence.