基于半平面上自然边界归化的无界区域上的Schwarz交替法及其离散化

1．引言由于并行技术的不断发展，人们越来越重视区域分解法的研究．对于闭曲线o外部的无界区域Ω上的椭圆边值问题，近年来基于自然边界归化理论[2,5,6,10]，提出了无界区域上的一类重叠型和不重叠型区域分解算法［11，12。13］，即将无界区域Ω分解为一个很小的有界区域Ω1和一个圆外无界区域Ω2，在Ω1和Ω2上分别有限元法和自然边界元法交替求解．其中，对于连续情形的重叠型区域分解法可利用投影理论得到意义下的几何收敛性［11］．对于连续情形和离散情形的重叠型区域分解法还可利用极值原理证明在最大模意义下的几何收敛性113］本文…     

无界区域上Stokes问题的自然边界元与有限元耦合法

§1.引言 对于用有限元方法求解平面有界区域上的Stokes问题,国内外已有大量工作,例如可见[2]、[9]及其所引文献.但对无界区域上的这一问题,由于区域的无界性给有限元方法带来了困难,边界元方法及边界元与有限元的耦合法便显示其优越性.本文提出用自然边界元与有限元的耦合法求解无界区域上的Stokes问题.这一耦合法早在作者以前的工作中被应用于求解调和问题、重调和问题和平面弹性问题,但将它用于求解     

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

１．引言 二十年来,自然边界元法已在椭圆问题求解方面取得了许多研究成果。它可以直接用来解决圆内（外）区域、扇形区域、球内（外）区域及半平面区域等特殊区域上的椭圆边值问题［１,２,５］,也可以结合有限元法求解一般区域上的椭圆边值问题,例如基于自然边界归化的耦合算法及区域分解算法就是处理断裂区域问题及外问题的一种有效手段［２－４,６］。 人们在设计求解外问题的耦合算法或者区域分解算法时,通常选取圆周或球面作人工边界。但对具有长条型内边界的外问题,以圆周或球面作人工边界显然并非最佳选择,它将会导致大量的…     

无界区域上基于自然边界归化的一种重叠型区域分解法及其离散化

1．引言由于科学技术的迅猛发展，人们遇到许多大规模科学和工程计算问题．随着并行计算机的出现和应用，并行技术越来越得到人们的重视和研究．区域分解法成为并行计算和处理这类问题的主要方法之一．但是，对于无界区域上的椭圆边值问题，因进行区域分解后至少有一个区域仍为无界区域，故仅应用通常的区域分解算法求解是不够的．由于边界归化是处理无界区域问题的有效手段，通常采用边界元和有限元耦合的方法求解此类问题IZ，6。8。121．或片什适当的人工边界并在此边界上加近似边界条件，再在有限区域应用有限元方法求解【人习．近年来…     

无界区域Stokes 问题非重叠型区域分解算法及其收敛性

本文研究无界区域Stokes方程外问题的利用有限元法和自然边界归化的非蕈叠型区域分解算法,此方法对无界区域Stokes问题非常有效.给出连续和离散情形的D-N算法及其收敛性分析,得到算法收敛的充要条件及充分条件,并得到最优的松弛因子和压缩因子,最后给出数值算例予以验证.     

曲边有限元和自然边界元求解外边值问题的耦合法

1　引　　言边界元与有限元耦合法在科学和工程计算中有着独特的作用 .由于区域的无限性给人们常用的有限元方法带来困难 ,边界元方法又难以独立处理非线性和非均质的问题以及具有不规则边界的区域上的问题 ,而两者相结合却可以克服各自的缺点 ,故边界元与有限元耦合法在处理一般区域问题特别是无界区域问题时便得到科学与工程界的青睐 ,获得了比较广泛的应用 .自然边界元方法并不引入新的变量 ,属于直接边界元方法[2 ] [8] .它保持能量不变和原边值问题的许多有用性质 ,例如双线性型的对称性和强制性 ,从而自然积分方程的解的存在唯一性及…     

抛物型初边值问题的自然积分方程及其数值解法

１．引言数值求解无界区域的偏微分方程,自然的处理方式是削去区域的无界部分,即引入一条适当的人工边界ｒ。,将原问题的求解限制在一个适当的有界区域Ｄ内,这样必须在人工边界上引入适当边界的条件．于是很自然地导致这样一个问题：'是否存在一个人工边界条件,使得在这边界条件下,原问题在区域Ｄ内所求得的数值解与原无界区域的解在Ｄ上的限制是完全一致的？"这里我们的着眼点是寻求与原无界区域问题等价的数学形式,以便于数值求解．因为边界元方法可以将区域内的问题转化到区域的边界上去处理,经典的边界元方法常被应用．七十年代…     

无界区域抛物方程自然边界元方法

本文应用自然边界元方法求解无界区域抛物型初边值问题。首先将控制方程对时间进行离散化,得到关于时间步长离散化的椭圆型问题。通过Fourier展开,导出相应问题的自然积分方程和Poisson积分公式。研究了自然积分算子的性质,并讨论了自然积分方程的数值解法,最后给出数值例子。从而解决了抛物型问题的自然边界归化和自然边界元方法。     

圆外区域求解Helmholtz方程

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

无穷凹角区域各向异性问题的自然边界元与有限元耦合法

本文研究无穷凹角区域上一类各向异性问题的自然边界元与有限元耦合法.利用自然边界归化原理,获得圆弧或椭圆弧人工边界上的自然积分方程,给出了耦合的变分形式及其数值方法,以及逼近解的收敛性和误差估计,最后给出了数值例子,以示方法的可行性和有效性.     

An efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations

In this paper, we propose an efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations. Firstly, we carry out a truncation from a three‐dimensional unbounded domain to a bounded spherical domain. By using spherical coordinate transformation and spherical harmonic expansion, we transform the original problem into a series of one‐dimensional eigenvalue problem that can be solved effectively. Secondly, we introduce a weighted Sobolev space to treat the singularity in the effective potential. Using the property of orthogonal polynomials in weighted Sobolev space, the error estimate for the approximate eigenvalues and corresponding eigenfunctions are proved. Error estimates show that our numerical method can achieve spectral accuracy for approximate eigenvalues and eigenfunctions. Finally, we give some numerical examples to demonstrate the efficiency of our algorithms and the correctness of the theoretical results.     

求解Fisher型方程的混合Legendre-球面调和谱方法

主要研究两同心球所界球形区域上偏微分方程的谱方法,建立了与区域形状相适应的混合Legendre-球面调和正交逼近的部分结果,在此基础上提出了数值求解两同心球所界球形区域上Fisher型方程的混合Legendre-球面调和谱格式,并分别给出了格式的收敛性及相关的数值结果.     

两同心球间旋转流动类Lorenz方程组的分歧分析

0引言两个同心旋转球之间的流动又称为球Couette流动.作为一个简单的模型,研究它能够为揭示流动失稳转捩至湍流这一重大理论课题的规律提供线索;同时,由于球Couette流动更象全球大气流动,研究它也能成为研究大气物理提供一个粗略的模型,为这一方面     

Furtivity and Masking Problems in Time-Dependent Electromagnetic Obstacle Scattering

In this paper, we consider furtivity and masking problems in time-dependent three-dimensional electromagnetic obstacle scattering. That is, we propose a criterion based on a merit function to minimize or to mask the electromagnetic field scattered by a bounded obstacle when hit by an incoming electromagnetic field and, with respect to this criterion, we drive the optimal strategy. These problems are natural generalizations to the context of electromagnetic scattering of the furtivity problem in time-dependent acoustic obstacle scattering presented in Ref. 1. We propose mathematical models of the furtivity and masking time-dependent three-dimensional electromagnetic scattering problems that consist in optimal control problems for systems of partial differential equations derived from the Maxwell equations. These control problems are approached using the Pontryagin maximum principle. We formulate the first-order optimality conditions for the control problems considered as exterior problems defined outside the obstacle for systems of partial differential equations. Moreover, the first-order optimality conditions derived are solved numerically with a highly parallelizable numerical method based on a perturbative series of the type considered in Refs. 2–3. Finally, we assess and validate the mathematical models and the numerical method proposed analyzing the numerical results obtained with a parallel implementation of the numerical method in several experiments on test problems. Impressive speedup factors are obtained executing the algorithms on a parallel machine when the number of processors used in the computation ranges between 1 and 100. Some virtual reality applications and some animations relative to the numerical experiments can be found in the website http://www.econ.unian.it/recchioni/w10/.     

THE OVERLAPPING DOMAIN DECOMPOSITION METHODFOR HARMONIC EQUATION OVER EXTERIOR THREE-DIMENSIONAL DOMAIN

1.IntroductionManyscientificandengineeringproblemscallbereducedtoexteriorboundaryvalueproblemsofpartialdifferentialequations.Althoughthenumericalmethodstosolveboundaryvalueproblems,suchasthefiniteelementmethodandthefinitedifferencemethod,areveryeffectiveonboundeddomains,yetweoftenfinditdifficulttousethemtodealwithunboundedproblems.Theboundaryreductionisaforcefulmeanstosolvesomeproblemsoverunboundeddomains.Amongmanyboundaryreductions,thenaturalboundaryreductionfoundedbyK.FengandD.H.Yuhassomed…     

Generalized combined field integral equations for the iterative solution of the three-dimensional Maxwell equations

The construction of new second-kind Fredholm integral equations for the numerical solution of problems of high-frequency electromagnetic scattering by a perfect conductor is proposed. These formulations are characterized by some eigenvalue clusterings. They are especially well adapted to Krylov subspace iterative solvers. Their derivation is based on the incorporation of a sufficiently accurate approximation of the exact operator linking the Cauchy data of the scattering boundary-value problem to the classical integral relations. This operator is related to the concept of the On-Surface Radiation Condition (OSRC). These formulations can be considered as a natural generalization of the well-known Brakhage–Werner and combined field integral equations. The efficiency of the approach is established through an analytical and numerical study in the spherical case.     

Stable Boundary Conditions and Discretization for $P_N$ Equations

A solution to the linear Boltzmann equation satisfies an energy bound, which reflects a natural fact: The energy of particles in a finite volume is bounded in time by the energy of particles initially occupying the volume augmented by the energy transported into the volume by particles entering the volume over time. In this paper, we present boundary conditions (BCs) for the spherical harmonic $(P_N)$ approximation, which ensure that this fundamental energy bound is satisfied by the $P_N$ approximation. Our BCs are compatible with the characteristic waves of $P_N$ equations and determine the incoming waves uniquely. Both, energy bound and compatibility, are shown on abstract formulations of $P_N$ equations and BCs to isolate the necessary structures and properties. The BCs are derived from a Marshak type formulation of BC and base on a non-classical even/odd-classification of spherical harmonic functions and a stabilization step, which is similar to the truncation of the series expansion in the $P_N$ method. We show that summation by parts (SBP) finite difference on staggered grids in space and the method of simultaneous approximation terms (SAT) allows to maintain the energy bound also on the semi-discrete level.     

Forced Vibration of Three-Layered Spherical and Ellipsoidal Shells under Axisymmetric Loads

A variant of vibration theory for three-layered shells of revolution under axisymmetric loads is elaborated by applying independent kinematic and static hypotheses to each layer, with account of transverse normal and shear strains in the core. Based on the Reissner variational principle for dynamic processes, equations of nonlinear vibrations and natural boundary conditions are obtained. The numerical method proposed for solving initial boundary-value problems is based on the use of integrodifferential approach for constructing finite-difference schemes with respect to spatial and time coordinates. Numerical solutions are obtained for dynamic deformations of open three-layered spherical and ellipsoidal shells, over a wide range of geometric and physical parameters of the core, for different types of boundary conditions. A comparative analysis is given for the results of investigating the dynamic behavior of three-layered shells of revolution by the equations proposed and the shell equations of Timoshenko and Kirhhoff-Love type, with the use of unified hypotheses across the heterogeneous structure of shells.     

On numerical density approximations of solutions of SDEs with unbounded coefficients

We study a numerical method to compute probability density functions of solutions of stochastic differential equations. The method is sometimes called the numerical path integration method and has been shown to be fast and accurate in application oriented fields. In this paper we provide a rigorous analysis of the method that covers systems of equations with unbounded coefficients. Working in a natural space for densities, L ¹, we obtain stability, consistency, and new convergence results for the method, new well-posedness and semigroup generation results for the related Fokker-Planck-Kolmogorov equation, and a new and rigorous connection to the corresponding probability density functions for both the approximate and the exact problems. To prove the results we combine semigroup and PDE arguments in a new way that should be of independent interest.