人工边界方法

人工边界方法是数值求解无界区域上偏微分方程的一类有效计算方法。本文以二维Poisson方程的外问题和声波方程为例,介绍人工边界方法的基本思想和核心技术。     

圆外区域求解Helmholtz方程

1引言许多科学和工程计算问题都可以归结为无界区域上的偏微分方程边值问题．而求解椭圆方程边值问题的常用技术是有限元方法,可是对于无界区域,在用有限元方法求解时,往往遇到困难．最简单的办法显然是直接略去区域的无界部分求解,但这样做或者导致过低的计算精度,或者要付出很高的计算代价．边界归化,即将求解偏微分方程边值问题转化为边界积分方程,是求解某些无界区域问题的强有力的手段．自70年代以来,有限元和     

数值解多维问题的外推与组合技术的若干新进展

本文综述近年来数值解多维问题的外推与组合技术的新进展，内容包括分裂外推及其在偏微分方程、多堆积分方程、多维数值积分中的应用；Ｃ．Ｚｅｎｇｅｒ的稀疏网格法与组合求解技术；以及解边界积分方程的组合方法，本文通过算例表明这些方法是非常有效的，是解多维问题的钥匙。     

Winkler地基上固支薄板自由振动问题的准Green函数方法

将准Green函数方法应用于求解Winkler地基上固支薄板的自由振动问题.即利用问题的基本解和边界方程构造一个准Green函数,这个函数满足了问题的齐次边界条件.采用Green公式,将Winkler地基上固支薄板自由振动问题的振型控制微分方程化为第二类Fredholm积分方程.通过边界方程的适当选择,积分方程核的奇异性被克服了.数值算例表明,该方法具有较高的精度,是一种有效的数学方法.     

新的三维力学GELD正演和反演算法

在本文中 ,我们提出了新的整体积分和局部微分GILD的力学正演和反演方法 .我们建立了弹性和塑性力学的体积分微分方程 .我们证明了这个体积分方程和伽辽金虚功原理等价 .新的GILD方法是基于这个体积分微分方程 .GL方法是进一步的发展 ,GL是一种整体场和局部场相互作用的全新方法 .在这个方法中 ,仅仅需要解 3× 3或者 6 × 6的局部小矩阵 .特别是 ,用GL方法求解无限域的偏微分方程时 ,不需要任何人工边界 ,不需要任何吸收边界条件和不需要任何边界积分方程 .新的三维力学GILD正演和反演算法已被应用研究奈米材料的力学性质的模拟计算 .我们获得非常好的奈米材料的力学变形的超拉力的力学性质 .我们提出了新的奈米地球物理新概念和发现了GILD数值量子     

边界元方法——一个老而新的发展中的数值方法

<正> 一、前言物理问题的数学模型的形成,一般有几种不同的形式,可以直接表示成偏微分方程的形式,也可表现为区域上的变分形式,或者归结为边界上的某个积分方程的形式.这些不同的数学形式在理论上是等价的,但在实践中不等效,它们分别导致有限差分法、有限单元法和边界元方法等不同的数值方法.     

一维热传导方程移动边界识别的计算方法

构造了一种正则化的积分方程方法来由Cauchy数据确定一维热传导方程的移动边界．在将区域延拓至规则区域后,通过Fourier方法将问题转化为一个第一类Volterra积分方程．然后分别用Lavrentiev正则化方法以及Tikhonov正则化方法将不稳定的第一类Volterra积分方程转化为适定的第二类积分方程,并分别将积分方程转化为常微分方程组,并用Runge—Kutta方法数值求解,以及直接离散来求解．最后通过自由边界上的条件得到数值的移动边界．通过一些数值试验表明此方法是有效可行的,并且给出的方法无需迭代,数值计算较简单．     

一类半线性对流扩散问题特征-有限体积法H1模误差估计

1引言有限体积法是由Baliga和Patankar提出的一种数值求解偏微分方程,特别是物理学中保持守恒律方程的有效方法．由于其运用原方程的体积积分公式和有限控制体积来离散方程．使方程在控制体积上保持守恒律这一重要的物理特性,自出现以来,有了很大的发展([2-4],[10])．特征线方法([1],[8],[9])则是一种非常适合求解对流占优扩散方程的数值     

一类带弱奇异核偏积分微分方程空间半离散的Legendre-Galerkin谱方法

本文给出了数值求解一类偏积分微分方程的空间半离散Legendre-Galerkin谱方法;证明了解的稳定性及收敛性。     

二阶和四阶椭圆型偏微分方程的镜像基本解方法及应用

将基本解方法推广到二阶和四阶椭圆型偏微分方程的对称问题,在边界上不需要处理奇异积分.通过坐标变换,将一般二阶和四阶椭圆型偏微分方程化为目前研究较为成熟的调和或双调和方程.再根据镜像法构造出适合对称条件的基本解函数,简化了计算,且不影响计算的精度.通过数值计算结果可以看出,利用镜像技术构造出的基本解,前期准备数据少,可保持精度,是一种有效的数值方法.     

Error estimates for shooting methods in two-point boundary value problems for second-order equations

This paper is concerned with a procedure for estimating the global discretization error arising when a boundary value problem for a system of second order differential equations is solved by the simple shooting method, without transforming the original problem in an equivalent first order problem. Expressions of the global discretization error are derived for both linear and nonlinear boundary value problems, which reduce the error estimation for a boundary value problem to that for an initial value problem of same dimension. The procedure extends to second order equations a technique for global error estimation given elsewhere for first order equations. As a practical result the accuracy of the estimates for a second order problem is increased compared with the estimates for the equivalent first order problem.     

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.

     

Interpolation error-based a posteriori error estimation for two-point boundary value problems and parabolic equations in one space dimension

Summary. I derive a posteriori error estimates for two-point boundary value problems and parabolic equations in one dimension based on interpolation error estimates. The interpolation error estimates are obtained from an extension of the error formula for the Lagrange interpolating polynomial in the case of symmetrically-spaced interpolation points. From this formula pointwise and seminorm a priori estimates of the interpolation error are derived. The interpolant in conjunction with the a priori estimates is used to obtain asymptotically exact a posteriori error estimates of the interpolation error. These a posteriori error estimates are extended to linear two-point boundary problems and parabolic equations. Computational results demonstrate the convergence of a posteriori error estimates and their effectiveness when combined with an hp-adaptive code for solving parabolic systems. Received April 17, 2000 / Revised version received September 25, 2000 / Published online May 30, 2001     

An augmented Galerkin procedure for the boundary integral method applied to mixed boundary value problems

Here we present the asymptotic error analysis for the boundary element approximation of the direct boundary integral equations for the plane mixed boundary value problem of the Laplacian. The boundary elements are defined by B-splines for the smooth parts of the boundary charges and additional singular functions at the collision points. The asymptotic error estimates include estimates for the stress intensity factors which occur as additional unknowns to be computed within the Galerkin scheme. The numerical analysis is based on the uniqueness of the problem, a coerciveness inequality, the triangular principal part and an extended shift theorem of the boundary integral operators.     

FINITE ELEMENT NONLINEAR GALERKIN COUPLING METHOD FOR THE EXTERIOR STEADY NAVIER-STOKES PROBLEM

1.IntroductionNonlinearGalerkinmethodsaremultilevelschemesforthedissipativeevolutionpartialdifferentialequations.Theycorrespondtothesplittingsoftheunknownu:u=y z)wherethecomponentsareofdifferentorderofmagnitudewithrespecttoaparameterrelatedtothespati...     

细菌模型的非协调有限元分析

将非协调元应用于描述细菌传播的反应扩散方程组的初边值问题.借助单元的一些特性和非协调误差估计技巧,分别在半离散和全离散有限元格式下,研究了其数值解与精确解的误差估计,得到了最优的误差估计以及超逼近结果.     

Meshless Galerkin algorithms for boundary integral equations with moving least square approximations

In this paper, we first give error estimates for the moving least square (MLS) approximation in the Hk norm in two dimensions when nodes and weight functions satisfy certain conditions. This two-dimensional error results can be applied to the surface of a three-dimensional domain. Then combining boundary integral equations (BIEs) and the MLS approximation, a meshless Galerkin algorithm, the Galerkin boundary node method (GBNM), is presented. The optimal asymptotic error estimates of the GBNM for three-dimensional BIEs are derived. Finally, taking the Dirichlet problem of Laplace equation as an example, we set up a framework for error estimates of the GBNM for boundary value problems in three dimensions.     

A posterior error estimates for the nonlinear grating problem with transparent boundary condition

The nonlinear grating problem is modeled by Maxwell's equations with transparent boundary conditions. The nonlocal boundary operators are truncated by taking sufficiently many terms in the corresponding expansions. A finite element method with the truncation operators is developed for solving the nonlinear grating problem. The two posterior error estimates are established. The a posterior error estimate consists of two parts: finite element discretization error and the truncation error of the nonlocal boundary operators. In particular, the truncation error caused by truncation operations is exponentially decayed when the parameter N is increased. Numerical experiment is included to illustrate the efficiency of the method. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1101–1118, 2015     

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 Solutions of Harmonic and Biharmonic Canonical Integral Equations in Interior or Exterior Circular Domains

Elliptic boundary-value problems can be reduced to integral equations on the boundary by many different ways. The canonical reduction, suggested by Prof. Feng Kang, is a natural and direct approach of boundary reduction. This paper gives the numerical method for solving harmonic and biharmonic canonical integral equations in interior or exterior circular domains, together with their convergence and error estimates. Using the theory of distributions, the difficulty caused by the singularities of integral kernel is overcome. Results of several numerical calculations verify the theoretical estimates.