NEW ALTERNATING DIRECTION FINITE ELEMENT SCHEME FOR NONLINEAR PARABOLIC EQUATION

A new alternating direction (AD) finite element (FE) scheme for 3-dimensional nonlinear parabolic equation and parabolic integro-differential equation is studied. By using AD, the 3-dimensional problem is reduced to a family of single space variable problems, calculation work is simplified; by using FE, high accuracy is kept; by using various techniques for priori estimate for differential equations such as inductive hypothesis reasoning, the difficulty arising from the nonlinearity is treated. For both FE and ADFE schemes, the convergence properties are rigorously demonstrated, the optimal H1-and L2-norm space estimates and the 0((△t)2) estimate for time variable are obtained.     

AD GALERKIN ANALYSIS FOR NONLINEAR PSEUDO—HYPERBOLIC EQUATIONS

AD(Alternating direction)Galerkin schemes for d-dimensional nonlinear pseudo-hyperbolic equations are studied.By using patch approximation technique,AD procedure is realized,and calculation,work is simplified.By using Galerkin approach,highly computational accuracy is kept.By using various priori estimate techniques for differential equations,difficulty coming form non-linearity is treated,and optimal H^1 and L^2 convergence prop-erties are demonstrated.Moreover,although all the existed AD Galerkin schemes using patch approximation are limited to have only one order accuracy in time increment,yet the schemes formulated in this paper have second order accuracy in it.This implies an essential advancement in AD Galerkin aualysis.     

AN A.D.I.GALERKIN METHOD FOR NONLINEAR PARABOLIC INTEGRO-DIFFERENTIAL EQUATION USING PATCH APPROXIMATION

An A. D. I. Galerkin scheme for three-dimensional nonlinear parabolic integro-differen-tial equation is studied. By using alternating-direction, the three-dimensional problem is reduced to a family of single space variable problems, the calculation is simplified; by using a local approxima-tion of the coefficients based on patches of finite elements, the coefficient matrix is updated at each time step; by using Ritz-Volterra projection, integration by part and other techniques, the influence coming from the integral term is treated; by using inductive hypothesis reasoning, the difficulty coming from the nonlinearity is treated. For both Galerkin and A. D. I. Galerkin schemes the con-vergence properties are rigorously demonstrated, the optimal H~1-norm and L~2-norm estimates are obtained.     

一类非线性对流扩散问题的FDSD预测校正格式

1．引言由Hughes和Brooks门提出,并经Johnson等人［‘－‘1发展的流线扩散法（Streamline－DiffusionMetho人以下简称SD方法）是求解对流占优扩散问题（包括纯双曲问题）的一种有效的数值方法．由于良好的数值稳定性及其高阶收敛率,SD方法已广泛地应用于计算流体等诸多科学工程计算．然而,传统的sD方法利用时一空有限元求解发展型问题,导致对高维问题工作量过于庞大;其编程实现较复杂,对非线性问题也不便进行线性化处理．为使SD方法能够较简便地应用于高维和非线性问题,孙撒问提出了仅对空间域作有限元离散,而对时间域作差分…     

AN ADI GALERKIN METHOD WITH MOVING FINITE ELEMENT SPACES FOR A CLASS OF SECOND-ORDER HYPERBOLIC EQUATIONS

1　 IntroductionADI Galerkin methods were first formulated for the solution of nonlinear parabolic andlinear second-order hyperbolic problems on rectangular regions by Douglas and Dupont[1 ] .These methods combine alternating-direction method and Galerkin finite element methodtogether.So,they have the advantage of reducing the solution of a multidimensional problemto the solution of sets of independent one-dimensional problems,decreasing the amount ofcalculation,natural parallelism and highe…     

带有记忆项和阻尼项的非线性双曲型方程的有限元方法

This paper studies the finite element method for some nonlinear hyperbolic partial differential equations with memory and dampling terms.A Crank-Nicolson approximation for this kind of equations is presented.By using the elliptic Ritz-Volterra projection,the analysis of the error estimates for the finite element numerical solutions and the optimal H1-norm error estimate are demonstrated.     

一类带弱奇异核非线性偏积分微分方程的全离散有限元

1引言我们将研究下面一类带弱奇异核非线性偏积分微分方程的数值解:u_t-▽·(a(u)▽u)-integral from n=0 to tβ(t-s)△u(s)ds=f(u),x∈Ω,t∈(?),(1.1) u(·,t)=0,x∈(?)Ω,t∈J,(1.2) u(·,0)=v(x),x∈Ω,(1.3)其中Ω为平面上的凸角域,J=(0,T],α和f为R上的光滑函数,满足0相似文献   

三维非线性抛物型积分-微分方程的A.D.I.Galerkin方法及其分析

研究三维非线性抛物型积分-微分方程的A.D.I.Galerkin方法.通过交替方向,化三维为一维,简化计算;通过Galerkin法,保持高精度.成功处理了Volterra项的影响;对所提Galerkin及A.D.I.Galerkin格式给出稳定性和收敛性分析,得到最佳H1和L2模估计.     

一类非线性反应扩散方程组的有限元分析

菸屏抗娴耐ǘ嘶故侵苟耍蟛饬康亩急匦胧撬堑闹芯叮皇堑ヒ恢芯丁Ｂ菸屏抗娴闹芯叮ǔ２捎萌敕ú饬俊？墒牵敕ú獬隽恐挡环螱B／T14791－93螺纹中径的定义。三针法直接测出的是螺纹量规的单一中径。按GB／T14791－93要求,螺纹的中径是一假想圆柱的直径,该圆柱的母线通过牙型上沟槽和凸起宽度相等的地方。一般应该用三针法测出单一中径,再在工具显微镜上测出螺距的实际偏差值,然后按下式计算出螺纹量规的中径氖当平獾拇嬖谖ㄒ恍约捌湮蟛钭钣牛葉（１）模和最优Ｌ~（２）模估计结果．在　２中,建立（１）的交替…     

非线性分数阶反应扩散方程组的间断时空有限元方法

利用时间间断空间连续的时空有限元方法构造了空间分数阶反应扩散方程组的可以逐时间层求解的全离散格式.在时间离散区间上,采用Radau积分公式,将插值理论与有限元理论相结合,给出了全离散格式解的存在唯一性结果,并证明了所给格式是无条件稳定的,进而详细给出最优阶L~∞(L~2)模误差估计过程.最后用数值算例验证了理论分析的正确性.     

GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS

1. IntroductionIn the numerical simulation of the Navier-Stokes equations one encounters three seriousdifficulties in the case of large Reynolds numbers f the treatment of the incomPressibility con-dition divu = 0, the treatment of the noIilinear terms and the large time integration. For thetreatment of the incoInPressibility condition, one use the penalty method in the case of finiteelemellts [1--2l and for the treatmen of the noulinar terms and the large tfor integration, oneuse the nonlin…     

非线性抛物方程混合有限元方法的高精度分析

采用双线性元及零阶Raviart-Thomas元（Q₁₁+Q₁₀×Q₀₁）对非线性抛物方程讨论了一种H¹-Galerkin混合有限元方法.提出一个线性化的二阶格式,利用数学归纳法有技巧的导出了原始变量u在H¹（Ω）模意义下及流量p=▽u在L²（Ω）模意义下的O（h²+τ²）阶超逼近性质.引入一个有关初始点的时间离散方程,并利用其得到了▽ ·在L²（Ω）模意义下的O（h²+τ²）阶的超逼近结果.同时利用插值后处理技巧得到整体超收敛.最后,数值算例结果验证了理论分析（其中,h是剖分参数,τ是时间步长）.     

THE UPWIND FINITE ELEMENT SCHEME AND MAXIMUM PRINCIPLE FOR NONLINEAR CONVECTION-DIFFUSION PROBLEM

In this paper, a kind of partial upwind finite element scheme is studied for twodimensional nonlinear convection-diffusion problem. Nonlinear convection term approximated by partial upwind finite element method considered over a mesh dual to the triangular grid, whereas the nonlinear diffusion term approximated by Galerkin method. A linearized partial upwind finite element scheme and a higher order accuracy scheme are constructed respectively. It is shown that the numerical solutions of these schemes preserve discrete maximum principle. The convergence and error estimate are also given for both schemes under some assumptions. The numerical results show that these partial upwind finite element scheme are feasible and accurate.     

ERROR ESTIMATES FOR TWO-SCALE COMPOSITE FINITE ELEMENT APPROXIMATIONS OF NONLINEAR PARABOLIC EQUATIONS

We study spatially semidiscrete and fully discrete two-scale composite finite element method for approximations of the nonlinear parabolic equations with homogeneous Dirich-let boundary conditions in a convex polygonal domain in the plane.This new class of finite elements,which is called composite finite elements,was first introduced by Hackbusch and Sauter[Numer.Math.,75(1997),pp.447-472]for the approximation of partial differential equations on domains with complicated geometry.The aim of this paper is to introduce an efficient numerical method which gives a lower dimensional approach for solving par-tial differential equations by domain discretization method.The composite finite element method introduces two-scale grid for discretization of the domain,the coarse-scale and the fine-scale grid with the degrees of freedom lies on the coarse-scale grid only.While the fine-scale grid is used to resolve the Dirichlet boundary condition,the dimension of the finite element space depends only on the coarse-scale grid.As a consequence,the resulting linear system will have a fewer number of unknowns.A continuous,piecewise linear composite finite element space is employed for the space discretization whereas the time discretization is based on both the backward Euler and the Crank-Nicolson methods.We have derived the error estimates in the L∞(L2)-norm for both semidiscrete and fully discrete schemes.Moreover,numerical simulations show that the proposed method is an efficient method to provide a good approximate solution.     

Navier-Stokes方程的最佳有限元非线性Galerkin算法

1．引言对于Navier－Stokes方程有限元数值求解方面的研究已有很多的文章和专著,多数是采用有限元Galerkin算法,例见文献[1-4]．然而,由于Navier-Stokes方程在大雷诺数时有其强的非线性性和对时间土的长期依赖性,用计算机求解Navier－Stokes方程在速度和容量方面是难以承受的．为了克服这些困难,最近人们提出了有限元非线性Galerkin算法,见文献卜8］,然而这种算法只是在某一有限时刻之后具有好的收敛速度,在初始时刻的某一区间不能达到好的收敛速度．本文应用Taylor展开技术导出了数值求解二维非定常Navier－Stokes方程的最佳…     

轴对称非线性磁场分布的有限元分析

一、引言 轴对称电磁结构在轴对称电流作用下的磁场分布问题本是三维问题,利用轴对称关系可将三维麦克斯韦方程化为形如下面(1)式的二维非线性椭圆型方程.它与标准的平面磁场方程有所不同.而且v(B)不是B的单增函数.实际工程中的大型圆柱变压器、平     

STRONG CONVERGENCE OF A FULLY DISCRETE FINITE ELEMENT METHOD FOR A CLASS OF SEMILINEAR STOCHASTIC PARTIAL DIFFERENTIAL EQUATIONS WITH MULTIPLICATIVE NOISE

This paper develops and analyzes a fully discrete finite element method for a class of semilinear stochastic partial differential equations(SPDEs)with multiplicative noise.The nonlinearity in the diffusion term of the SPDEs is assumed to be globally Lipschitz and the nonlinearity in the drift term is only assumed to satisfy a one-sided Lipschitz condition.These assumptions are the same ones as the cases where numerical methods for general nonlinear stochastic ordinary differential equations(SODEs)under"minimum assumptions"were studied.As a result,the semilinear SPDEs considered in this paper are a direct generalization of these nonlinear SODEs.There are several difficulties which need to be overcome for this generalization.First,obviously the spatial discretization,which does not appear in the SODE case,adds an extra layer of difficulty.It turns out a spatial discretization must be designed to guarantee certain properties for the numerical scheme and its stiffness matrix.In this paper we use a finite element interpolation technique to discretize the nonlinear drift term.Second,in order to prove the strong convergence of the proposed fully discrete finite element method,stability estimates for higher order moments of the H1-seminorm of the numerical solution must be established,which are difficult and delicate.A judicious combination of the properties of the drift and diffusion terms and some nontrivial techniques is used in this paper to achieve the goal.Finally,stability estimates for the second and higher order moments of the L2-norm of the numerical solution are also difficult to obtain due to the fact that the mass matrix may not be diagonally dominant.This is done by utilizing the interpolation theory and the higher moment estimates for the H1-seminorm of the numerical solution.After overcoming these difficulties,it is proved that the proposed fully discrete finite element method is convergent in strong norms with nearly optimal rates of convergence.Numerical experiment results are also presented to validate the theoretical results and to demonstrate the efficiency of the proposed numerical method.     

A coupling method based on new MFE and FE for fourth-order parabolic equation

In this article, a coupling method of new mixed finite element (MFE) and finite element (FE) is proposed and analyzed for fourth-order parabolic partial differential equation. First, the fourth-order parabolic equation is split into the coupled system of second-order equations. Then, an equation is solved by finite element method, the other equation is approximated by the new mixed finite element method, whose flux belongs to the square integrable space replacing the classical H(div;Ω) space. The stability for fully discrete scheme is derived, and both semi-discrete and fully discrete error estimates are obtained. Moreover, the optimal a priori error estimates in L ² and H ¹-norm for both the scalar unknown u and the diffusion term γ and a priori error estimate in (L ²)²-norm for its flux σ are derived. Finally, some numerical results are provided to validate our theoretical analysis.     

非线性抛物型方程的变网格有限元法

Three linear finite element schemes with moving mesh are established for a class of nonlinear parabolic equations.Optimal L^2-norm and energy norm error estimates are derived under certain conditions.     

NONLINEAR GALERKIN METHOD FOR NAVIER-STOKES EQUATIONS WITH STREAM-VORTICITY FORM

I Let D be a launder domain in Rz with Lipehitz continuous bodare I'. ac consider the stream-vortidty form of hmedependent Navier-stokes Muahons describing the 'flow Of a ~ incompreSSible nuid confined in Dwhere to and & are vorticity and stream function. BecaUSe the action (1) doeS not include the differentialten z, ~ tie ~ con&tion of' dab at In tthe Paper, we give a ho element n~ Galerkin ~, acs ~ is ~ an tab finite element spaceS X. and X* for the aPPmxhaation of the ~ ac v~ fUncti…