A nonoverlapping domain decomposition method for nonlinear physical processes

A coupling technique for a nonoverlapping domain decomposition method is investigated in this paper. We test this method on a nonlinear heat conduction process taking place in a multi-chip module which has many geometrically structured subdomains. The problem of interface condition between two adjacent subdomains is tackled by incorporating a defect equation which is solved iteratively. We evaluate two defect equations on this scheme, i.e. the difference of normal derivative and the residual of the heat conduction equation itself. Both equations lead to systems of nonlinear equations which are solved by means of quasi-Newton methods with an adaptive alpha rate instead of a Jacobian matrix. The simulation suggests that the second defect equation is much more accurate. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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

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

Projection solution methods for one nonlinear singular integral equation

In this paper we consider a general projection method for the solution of a nonlinear singular integral equation and its applications in the method of orthogonal polynomials, the subdomains method, and the collocation method.     

Multidomain direct method and its applicability to the local time-step concept

A new multidomain direct (noniterative) method for solving boundary-value problems is presented. Using this method, the solution is expressed by a linear combination of auxiliary functions and unknowns which pertain to the boundaries of a subdomain. This approach enables us to solve problems independently and exactly, without any iterations between subdomains. As a consequence, different types of equations and coordinate systems may be considered in different subdomains. Moreover, different boundary conditions and variable (in space) time steps may be imposed on the subdomains as well. Applications are given for initial boundary-value problems with known analytical solutions, including a highly nonlinear porosity equation. © 1994 John Wiley & Sons, Inc.     

多孔介质驱动问题的混合元最小二乘特征有限元方法及其收敛分析

１．引言多孔介质二相驱动问题的数学模型是偶合的非线性偏微分方程组的初边值问题．该问题可转化为压力方程和浓度方程［１－４］．浓度方程一般是对流占优的对流扩散方程,它的对流速度依赖于比浓度方程的扩散系数大得多的Ｆａｒｃｙ速度．因此Ｄａｒｃｙ速度的求解精度直接影响着浓度的求解精度．为了提高速度的求解精度,７０年代Ｐ．Ａ．Ｒａｖｉａｔ和Ｊ．Ｍ．Ｔｈｏｍａｓ提出混合有限元方法［５］．Ｊ．ＤｏｕｇｌａｓＪｒ,Ｔ．Ｆ．Ｒｕｓｓｅｌｌ,Ｒ．Ｅ．Ｅｗｉｎｇ,Ｍ．Ｆ．Ｗｈｅｅｌｅｒ［１］－［４］,［９］,［１２］袁…     

Fractional Crank–Nicolson–Galerkin finite element scheme for the time‐fractional nonlinear diffusion equation

This article presents a finite element scheme with Newton's method for solving the time‐fractional nonlinear diffusion equation. For time discretization, we use the fractional Crank–Nicolson scheme based on backward Euler convolution quadrature. We discuss the existence‐uniqueness results for the fully discrete problem. A new discrete fractional Gronwall type inequality for the backward Euler convolution quadrature is established. A priori error estimate for the fully discrete problem in L²(Ω) norm is derived. Numerical results based on finite element scheme are provided to validate theoretical estimates on time‐fractional nonlinear Fisher equation and Huxley equation.     

Solving the 2-D elliptic Monge-Ampère equation by a Kansa’s method

In this paper, a Kansa’s method is designed to solve numerically the Monge-Ampère equation. The primitive Kansa’s method is a meshfree method which applying the combination of some radial basis functions (such as Hardy’s MQ) to approximate the solution of the linear parabolic, hyperbolic and elliptic problems. But this method is deteriorated when is used to solve nonlinear partial differential equations. We approximate the solution in some local triangular subdomains by using the combination of some cubic polynomials. Then the given problems can be computed in each subdomains independently. We prove the stability and convergence of the new method for the elliptic Monge-Ampère equation. Finally, some numerical experiments are presented to demonstrate the theoretical results.     

多孔介质中可压缩可混溶驱动问题的有限体积元法

有界区域上多孔介质中可压缩可混溶驱动问题由两个非线性抛物型方程耦合而成：压力方程和饱和度方程均是抛物型方程．运用有限体积元法对两个方程进行数值分析,给出了全离散有限体积元格式,并通过详细的理论分析,得到了近似解与原问题真解的最优H^1模误差估计。     

三维两相渗流驱动问题迎风区域分裂显隐差分法

对三维两相渗流驱动问题提出了两种迎风区域分裂显隐差分格式．压力方程采用了七点差分格式,为了能达到实际并行计算的要求,对饱和度方程采用了迎风区域分裂差分法,内边界处和各子区域分别对应显隐格式．得到了离散l2模收敛性分析,最后给出数值试验,支撑了理论分析结果．     

A mixed and discontinuous Galerkin finite volume element method for incompressible miscible displacement problems in porous media

The incompressible miscible displacement problem in porous media is modeled by a coupled system of two nonlinear partial differential equations, the pressure‐velocity equation and the concentration equation. In this article, we present a mixed finite volume element method for the approximation of pressure‐velocity equation and a discontinuous Galerkin finite volume element method for the concentration equation. A priori error estimates in L^∞(L²) are derived for velocity, pressure, and concentration. Numerical results are presented to substantiate the validity of the theoretical results. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR A FORWARD-BACKWARD HEAT EQUATION

A space-time finite element method,discontinuous in time but continuous in space, is studied to solve the nonlinear forward-backward heat equation. A linearized technique is introduced in order to obtain the error estimates of the approximate solutions. And the numerical simulations are given.     

A parallel two-level finite element variational multiscale method for the Navier–Stokes equations

A combination method of two-grid discretization approach with a recent finite element variational multiscale algorithm for simulation of the incompressible Navier–Stokes equations is proposed and analyzed. The method consists of a global small-scale nonlinear Navier–Stokes problem on a coarse grid and local linearized residual problems in overlapped fine grid subdomains, where the numerical form of the Navier–Stokes equations on the coarse grid is stabilized by a stabilization term based on two local Gauss integrations at element level and defined by the difference between a consistent and an under-integrated matrix involving the gradient of velocity. By the technical tool of local a priori estimate for the finite element solution, error bounds of the discrete solution are estimated. Algorithmic parameter scalings are derived. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the method.     

EXTRAPOLATION AND A-POSTERIORI ERROR ESTIMATORS OF PETROV-GALERKIN METHODS FOR NON-LINEAR VOLTERRA INTEGRO-DIFFERENTIAL EQUATIONS

1. IntroductionThe pmpme of tabs Paper is to show that the ~ardson edrapolation can be used toenhance the nUmerical solutions generated by a cab of Petrov-Gaierkin lhate element methodsfor the nonlinear VOlterra integrO-chrential equation (VIDE):where j = j(t,y): I x R --+ R and k = k(t,8,g): D x R - R (with D:= {(t,8): 0 S & S t ST}) denote given hmctions.Throughout tab paperl it will always be assumed that the problem (1.1) possesses a piquesolution y E C'(I), namely, the given hmc…     

A new mixed scheme based on variation of constants for Sobolev equation with nonlinear convection term

A new mixed scheme which combines the variation of constants and the H ¹-Galerkin mixed finite element method is constructed for nonlinear Sobolev equation with nonlinear convection term. Optimal error estimates are derived for both semidiscrete and fully discrete schemes. Finally, some numerical results are given to confirm the theoretical analysis of the proposed method.     

Local and parallel finite element algorithm based on the partition of unity method for the incompressible MHD flow

Based on the partition of unity method (PUM), a local and parallel finite element method is designed and analyzed for solving the stationary incompressible magnetohydrodynamics (MHD). The key idea of the proposed algorithm is to first solve the nonlinear system on a coarse mesh, divide the globally fine grid correction into a series of locally linearized residual problems on some subdomains derived by a class of partition of unity, then compute the local subproblems in parallel, and obtain the globally continuous finite element solution by assembling all local solutions together by the partition of unity functions. The main feature of the new method is that the partition of unity provide a flexible and controllable framework for the domain decomposition. Finally, the efficiency of our theoretical analysis is tested by numerical experiments.     

Piecewise homotopy analysis method and convergence analysis for formally well-posed initial value problems

In this paper, we apply the Schwarz waveform relaxation (SWR) method to the one-dimensional Schrödinger equation with a general linear or a nonlinear potential. We propose a new algorithm for the Schrödinger equation with time-independent linear potential, which is robust and scalable up to 500 subdomains. It reduces significantly computation time compared with the classical algorithms. Concerning the case of time-dependent linear potential or the nonlinear potential, we use a preprocessed linear operator for the zero potential case as a preconditioner which leads to a preconditioned algorithm. This ensures high scalability. In addition, some newly constructed absorbing boundary conditions are used as the transmission conditions and compared numerically.     

A two-level additive Schwarz method for the Morley nonconforming element approximation of a nonlinear biharmonic equation

In this paper, we consider the well known Morley nonconformingelement approximation of a nonlinear biharmonic equation whichis related to the well-known two-dimensional Navier–Stokesequations. Firstly, optimal energy and H¹-norm estimates areobtained. Secondly, a two-level additive Schwarz method is presentedfor the discrete nonlinear algebraic system. It is shown thatif the Reynolds number is sufficiently small, the two-levelSchwarz method is optimal, i.e. the convergence rate of theSchwarz method is independent of the mesh size and the numberof subdomains.     

Schwarz preconditioners for the spectral element discretization of the steady Stokes and Navier-Stokes equations

Summary. The - spectral element discretization of the Stokes equation gives rise to an ill-conditioned, indefinite, symmetric linear system for the velocity and pressure degrees of freedom. We propose a domain decomposition method which involves the solution of a low-order global, and several local problems, related to the vertices, edges, and interiors of the subdomains. The original system is reduced to a symmetric equation for the velocity, which can be solved with the conjugate gradient method. We prove that the condition number of the iteration operator is bounded from above by , where C is a positive constant independent of the degree N and the number of subdomains, and is the inf-sup condition of the pair -. We also consider the stationary Navier-Stokes equations; in each Newton step, a non-symmetric indefinite problem is solved using a Schwarz preconditioner. By using an especially designed low-order global space, and the solution of local problems analogous to those decribed above for the Stokes equation, we are able to present a complete theory for the method. We prove that the number of iterations of the GMRES method, at each Newton step, is bounded from above by . The constant C does not depend on the number of subdomains or N, and it does not deteriorate as the Newton iteration proceeds. Received March 2, 1998 / Revised version received October 12, 1999 / Published online March 20, 2001     

An element agglomeration nonlinear additive Schwarz preconditioned Newton method for unstructured finite element problems

This paper extends previous results on nonlinear Schwarz preconditioning (Cai and Keyes 2002) to unstructured finite element elliptic problems exploiting now nonlocal (but small) subspaces. The nonlocal finite element subspaces are associated with subdomains obtained from a non-overlapping element partitioning of the original set of elements and are coarse outside the prescribed element subdomain. The coarsening is based on a modification of the agglomeration based AMGe method proposed in Jones and Vassilevski 2001. Then, the algebraic construction from Jones, Vassilevski and Woodward 2003 of the corresponding non-linear finite element subproblems is applied to generate the subspace based nonlinear preconditioner. The overall nonlinearly preconditioned problem is solved by an inexact Newton method. A numerical illustration is also provided.This work was performed under the auspices of the U.S. Department of Energy by the University of California Lawrence Livermore National Laboratory: contract/grant number: W-7405-Eng-48. The contribution of the second author was also partially supported by Polish Scientific Grant 2/P03A/005/24.     

杂交有限元的区域分解法

近几年来,由于并行计算机的迅速发展,求解椭圆型微分方程的区城分解法又引起了人们的重视.这一方法的基本思想是把求解区域分解成许多子区域,每个子区域上用一台计算机求解.这种方法适应并行机的需要,是用并行机解大型椭圆型偏微分方程的一