带移动奇异源热方程的移动网格方法

研究了一个带若干奇异源热方程的数值求解,其源的移动由一个常微分方程描述.基于移动观察区域和区域分解思想提出了一个移动网格预估校正算法.网格方程可自然的通过并行高效求解,算法避免了跳跃信息[u]的计算而使物理方程的离散格式变得非常简单,且仍保持了空间上的二阶收敛性.数值例子验证了算法的收敛性和高效性,并模拟了非线性源函数带来的爆破现象.     

椭圆外区域上Helmholtz问题的耦合法

以椭圆外区域上Helmholtz方程为例,研究一种带有椭圆人工边界的自然边界元与有限元耦合法,给出了耦合变分问题的适定性及误差分析并给出数值例子.理论分析及数值结果表明,用方法求解椭圆外问题是十分有效的.为求解具有长条型内边界外Helmholtz问题提供了一种很好的数值方法.     

求解修正的Helmholtz方程关于非光滑区域问题

研究了修正的(纯虚波数)Helmholtz方程在阻尼边界条件下,求解含单个角点的闭区域问题.通过采用单双层混合位势来表示其解,进而对其角型区域进行求解.最后,通过数值例子来说明此方法的可行性与可靠性.     

单调迭代结合虚拟区域法求解非线性障碍问题

讨论了二阶半线性椭圆方程障碍问题的数值求解问题.用单调迭代算法求解障碍问题,并用改进的虚拟区域法求解相关的不规则区域上具有Dirichlet边界条件的椭圆方程.在计算过程中,传统的有限元离散会导致用扩展区域规则网格计算不规则物体边界上积分的困难.为了克服此困难,给出了一种新的基于有限差分的算法,从而使得偏微分快速算法可用.算法结构简单,易于编程实现.对有扩散和增长障碍的logistic人口模型数值模拟说明算法可行且高效.     

Schwarz 算法的 Lions 框架与异步并行算法的收敛性证明

经典的Schwarz算法,早在1870年就提出了,是求解不规则椭圆型方程的交替法.本世纪苏联学者等又在变分框架下论证了收敛性.近年来以 Schwarrz算法为基础的区域分解算法,发展十分活跃,当前该方法正与并行算法,预处理,快速直接解,多水平及多层网格技术结合,成为计算数学领域内最有前途一个分支.虽然表面看 Schwarz 算法不是并行的,但康立山等打破了分解区域为两子域的     

求解一类Helmholtz方程Cauchy问题的中心差分正则化方法

Helmholtz方程Cauchy问题是严重不适定问题,本文我们在一个带形区域上考虑了一类Helmholtz方程Cauchy问题:已知Cauchy数据u(0,y)=g(y),在区间0＜x＜1上求解.我们用半离散的中心差分方法得到了这一问题的正则化解,给出了正则化参数的选取规则,得到了误差估计.     

区域分解算法——偏微数值解新技术(Ⅲ)

<正> 3.1经典Sohwar二交替方法 重迭型区域分解算法是以Schwarz交替法为理论依据.1869年德国数学家H.A.Schwarz首次用交替方法论证两个相互重迭区域和集上的Laplaoe方程Diriohlet问题解的存在性,稍后Neumann注意到这一思想可以用于求解两个相互覆盖区域的Diriohlet     

一类各向异性外问题的重叠型区域分解算法

本文以椭圆外调和问题的自然边界归化为基础，提出了求解各向异性常系数椭圆方程的一种重叠型区域分解算法，并分析了算法的收敛性及收敛速度．理论分析及数值实验表明，该方法对于求解各向异性外问题非常有效．     

一种求解修正的Helmholtz方程Cauchy问题的数值方法

本文研究了无界带形区域Ω={(x,y)|0相似文献   

半圆域内的二维线性椭圆偏微分方程（英文）

本文研究了半圆域内的二维线性椭圆偏微分方程.利用Fokas提出的求解凸多边形区域内的线性椭圆偏微分方程的变换方法,我们改进了这个方法来研究半圆域内Laplace方程,修改Helmholtz方程和Helmholtz方程的解,并且导出了这些方程解的积分表达式,讨论了Helmholtz方程的广义Dirichlet到Neumann映射.     

A Source Transfer Domain Decomposition Method for Helmholtz Equations in Unbounded Domain Part II: Extensions

In this paper we extend the source transfer domain decomposition method (STDDM) introduced by the authors to solve the Helmholtz problems in two-layered media, the Helmholtz scattering problems with bounded scatterer, and Helmholtz problems in 3D unbounded domains. The STDDM is based on the decomposition of the domain into non-overlapping layers and the idea of source transfer which transfers the sources equivalently layer by layer so that the solution in the final layer can be solved using a PML method defined locally outside the last two layers. The details of STDDM is given for each extension. Numerical results are presented to demonstrate the efficiency of STDDM as a preconditioner for solving the discretization problem of the Helmholtz problems considered in the paper.     

A residual-based a posteriori error estimator for a two-dimensional fluid–solid interaction problem

In this paper, we develop an a posteriori error analysis of a mixed finite element method for a fluid–solid interaction problem posed in the plane. The media are governed by the acoustic and elastodynamic equations in time-harmonic regime, respectively, and the transmission conditions are given by the equilibrium of forces and the equality of the normal displacements of the solid and the fluid. The coupling of primal and dual-mixed finite element methods is applied to compute both the pressure of the scattered wave in the linearized fluid and the elastic vibrations that take place in the elastic body. The finite element subspaces consider continuous piecewise linear elements for the pressure and a Lagrange multiplier defined on the interface, and PEERS for the stress and rotation in the solid domain. We derive a reliable and efficient residual-based a posteriori error estimator for this coupled problem. Suitable auxiliary problems, the continuous inf-sup conditions satisfied by the bilinear forms involved, a discrete Helmholtz decomposition, and the local approximation properties of the Clément interpolant and Raviart–Thomas operator are the main tools for proving the reliability of the estimator. Then, Helmholtz decomposition, inverse inequalities, and the localization technique based on triangle-bubble and edge-bubble functions are employed to show the efficiency. Finally, some numerical results confirming the reliability and efficiency of the estimator are reported.     

A nonconforming domain decomposition approximation for the Helmholtz screen problem with hypersingular operator

We present and analyze a nonconforming domain decomposition approximation for a hypersingular operator governed by the Helmholtz equation in three dimensions. This operator appears when considering the corresponding Neumann problem in unbounded domains exterior to open surfaces. We consider small wave numbers and low‐order approximations with Nitsche coupling across interfaces. Under appropriate assumptions on mapping properties of the weakly singular and hypersingular operators with Helmholtz kernel, we prove that this method converges almost quasioptimally, that is, with optimal orders reduced by an arbitrarily small positive number. Numerical experiments confirm our error estimate. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 125–141, 2017     

Orthogonal Helmholtz decomposition in arbitrary dimension using divergence-free and curl-free wavelets

We present tensor-product divergence-free and curl-free wavelets, and define associated projectors. These projectors enable the construction of an iterative algorithm to compute the Helmholtz decomposition of any vector field, in wavelet domain. This decomposition is localized in space, in contrast to the Helmholtz decomposition calculated by Fourier transform. Then we prove the convergence of the algorithm in dimension two for any kind of wavelets, and in larger dimension for the particular case of Shannon wavelets. We also present a modification of the algorithm by using quasi-isotropic divergence-free and curl-free wavelets. Finally, numerical tests show the validity of this approach for a large class of wavelets.     

二维Helmholtz方程外问题基于自然边界归化的非重叠型区域分解算法

１．引言 设是平面光滑闭曲线,是以为边界的外部区域,考虑二维Ｈｅｌｍｈｏｌｔｚ方程外Ｎｅｕｍａｎｎ问题并在无穷远处满足Ｓｏｍｍｅｒｆｅｌｄ辐射条件其中 ｉ＝是区域的边界 的外法线方向,即指向由 包围的内部区域． ｋ在许多情况下（例如约化波动方程）是实数,在另一些情况下则是纯虚数,本文仅讨论ｋ为纯虚数的情况,且不失一般性,可设Ｉｍ（ｋ）＞０． 用某些数值方法求解线性抛物型方程或线性双曲型方程的初边值问题时,可能间接地导致求解Ｈｅｌｍｈｏｌｔ。方程的外问题［１０,１１;１２;１３］．例如,用自然边界无法求…     

Preconditioned iterative methods on sparse subspaces

When some rows of the system matrix and a preconditioner coincide, preconditioned iterations can be reduced to a sparse subspace. Taking advantage of this property can lead to considerable memory and computational savings. This is particularly useful with the GMRES method. We consider the iterative solution of a discretized partial differential equation on this sparse subspace. With a domain decomposition method and a fictitious domain method the subspace corresponds a small neighborhood of an interface. As numerical examples we solve the Helmholtz equation using a fictitious domain method and an elliptic equation with a jump in the diffusion coefficient using a separable preconditioner.     

A posteriori error estimation for the dual mixed finite element method of the elasticity problem in a polygonal domain

In this article, we propose a residual based reliable and efficient error estimator for the new dual mixed finite element method of the elasticity problem in a polygonal domain, introduced by M. Farhloul and M. Fortin. With the help of a specific generalized Helmholtz decomposition of the error on the strain tensor and the classical decomposition of the error on the gradient of the displacements, we show that our global error estimator is reliable. Efficiency of our estimator follows by using classical inverse estimates. The lower and upper error bounds obtained are uniform with respect to the Lamé coefficient λ, in particular avoiding locking phenomena. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

椭圆外区域上Helmholtz问题的自然边界元法

本文研究椭圆外区域上Helmholtz方程边值问题的自然边界元法.利用自然边界归化原理,获得该问题的Poisson积分公式及自然积分方程,给出了自然积分方程的数值方法.由于计算的需要,我们详细地讨论了Mathieu函数的计算方法(当0相似文献   

Boundary value problems for asymptotically homogeneous elliptic second order operators

In this paper we show the Dirichlet and Neumann problems over exterior regions have unique solutions in certain weighted Sobolev spaces. Two applications are given: (1) The Dirichlet problem for semi-linear operators, and (2) a Helmholtz decomposition for vector fields on exterior regions.