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.     

Numerical method of solution to loaded nonlocal boundary value problems for ordinary differential equations

A numerical method is suggested for solving systems of nonautonomous loaded linear ordinary differential equations with nonseparated multipoint and integral conditions. The method is based on the convolution of integral conditions into local ones. As a result, the original problem is reduced to an initial value (Cauchy) problem for systems of ordinary differential equations and linear algebraic equations. The approach proposed is used in combination with the linearization method to solve systems of loaded nonlinear ordinary differential equations with nonlocal conditions. An example of a loaded parabolic equation with nonlocal initial and boundary conditions is used to show that the approach can be applied to partial differential equations. Numerous numerical experiments on test problems were performed with the use of the numerical formulas and schemes proposed.     

The method of fundamental solutions with eigenfunction expansion method for nonhomogeneous diffusion equation

In this article we describe a numerical method to solve a nonhomogeneous diffusion equation with arbitrary geometry by combining the method of fundamental solutions (MFS), the method of particular solutions (MPS), and the eigenfunction expansion method (EEM). This forms a meshless numerical scheme of the MFS‐MPS‐EEM model to solve nonhomogeneous diffusion equations with time‐independent source terms and boundary conditions for any time and any shape. Nonhomogeneous diffusion equation with complex domain can be separated into a Poisson equation and a homogeneous diffusion equation using this model. The Poisson equation is solved by the MFS‐MPS model, in which the compactly supported radial basis functions are adopted for the MPS. On the other hand, utilizing the EEM the diffusion equation is first translated to a Helmholtz equation, which is then solved by the MFS together with the technique of the singular value decomposition (SVD). Since the present meshless method does not need mesh generation, nodal connectivity, or numerical integration, the computational effort and memory storage required are minimal as compared with other numerical schemes. Test results for two 2D diffusion problems show good comparability with the analytical solutions. The proposed algorithm is then extended to solve a problem with irregular domain and the results compare very well with solutions of a finite element scheme. Therefore, the present scheme has been proved to be very promising as a meshfree numerical method to solve nonhomogeneous diffusion equations with time‐independent source terms of any time frame, and for any arbitrary geometry. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

NUMERICAL SOLUTIONS OF PARABOLIC PROBLEMS ON UNBOUNDED 3-D SPATIAL DOMAIN

In this paper,the numerical solutions of heat equation on 3-D unbounded spatial do-main are considered. n artificial boundary Γ is introduced to finite the computationaldomain.On the artificial boundary Γ,the exact boundary condition and a series of approx-imating boundary conditions are derived,which are called artificial boundary conditions.By the exact or approximating boundary condition on the artificial boundary,the originalproblem is reduced to an initial-boundary value problem on the bounded computationaldomain,which is equivalent or approximating to the original problem.The finite differencemethod and finite element method are used to solve the reduced problems on the finitecomputational domain.The numerical results demonstrate that the method given in thispaper is effective and feasible.     

The spectral methods for parabolic Volterra integro-differential equations

In this paper we study the numerical solutions to parabolic Volterra integro-differential equations in one-dimensional bounded and unbounded spatial domains. In a bounded domain, the given parabolic Volterra integro-differential equation is converted to two equivalent equations. Then, a Legendre-collocation method is used to solve them and finally a linear algebraic system is obtained. For an unbounded case, we use the algebraic mapping to transfer the problem on a bounded domain and then apply the same presented approach for the bounded domain. In both cases, some numerical examples are presented to illustrate the efficiency and accuracy of the proposed method.     

The numerical solution of a periodic parabolic problem subject to a nonlinear boundary condition

The numerical solution by finite differences of a periodic parabolic problem subject to a nonlinear boundary condition is considered. It is shown that Newton's method can be used to solve the nonlinear equations provided a suitable initial approximation is known, and a method for constructing this first approximation is given.     

非均匀Reissner板弯曲的精确元法

本文在阶梯折算法和精确解析法的基础上,提出构造有限元的新方法——精确元法.该方法不用变分原理,可适用于任意变系数正定和非正定偏微分方程.利用该方法,得到Reissuer板弯曲的一个非协调单元,它具有十五个自由度.由于节点位移参数仅含有挠度和转角,因此处理任意边界条件非常容易.文中给出证明,位移和内力均收敛于精确解.由精确元法所得到的单元不仅能用于厚板,也可用于薄板.文末给出四个算例.算例表明,利用本文的方法,可获得满意的结果,并有较高的数值精度.     

A fast finite difference method for biharmonic equations on irregular domains and its application to an incompressible Stokes flow

Biharmonic equations have many applications, especially in fluid and solid mechanics, but is difficult to solve due to the fourth order derivatives in the differential equation. In this paper a fast second order accurate algorithm based on a finite difference discretization and a Cartesian grid is developed for two dimensional biharmonic equations on irregular domains with essential boundary conditions. The irregular domain is embedded into a rectangular region and the biharmonic equation is decoupled to two Poisson equations. An auxiliary unknown quantity Δu along the boundary is introduced so that fast Poisson solvers on irregular domains can be used. Non-trivial numerical examples show the efficiency of the proposed method. The number of iterations of the method is independent of the mesh size. Another key to the method is a new interpolation scheme to evaluate the residual of the Schur complement system. The new biharmonic solver has been applied to solve the incompressible Stokes flow on an irregular domain.      

Numerical treatment of degenerate diffusion equations via Feller's boundary classification,and applications

A numerical method is devised to solve a class of linear boundary‐value problems for one‐dimensional parabolic equations degenerate at the boundaries. Feller theory, which classifies the nature of the boundary points, is used to decide whether boundary conditions are needed to ensure uniqueness, and, if so, which ones they are. The algorithm is based on a suitable preconditioned implicit finite‐difference scheme, grid, and treatment of the boundary data. Second‐order accuracy, unconditional stability, and unconditional convergence of solutions of the finite‐difference scheme to a constant as the time‐step index tends to infinity are further properties of the method. Several examples, pertaining to financial mathematics, physics, and genetics, are presented for the purpose of illustration. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Weak Galerkin finite element method for valuation of American options

We introduce a weak Galerkin finite element method for the valuation of American options governed by the Black-Scholes equation. In order to implement, we need to solve the optimal exercise boundary and then introduce an artificial boundary to make the computational domain bounded. For the optimal exercise boundary, which satisfies a nonlinear Volterra integral equation, it is resolved by a higher-order collocation method based on graded meshes. With the computed optimal exercise boundary, the front-fixing technique is employed to transform the free boundary problem to a one- dimensional parabolic problem in a half infinite area. For the other spatial domain boundary, a perfectly matched layer is used to truncate the unbounded domain and carry out the computation. Finally, the resulting initial-boundary value problems are solved by weak Galerkin finite element method, and numerical examples are provided to illustrate the efficiency of the method.     

Numerical Solution of Partial Differential Equations in Arbitrary Shaped Domains Using Cartesian Cut-Stencil Finite Difference Method. Part I: Concepts and Fundamentals

A new finite difference (FD) method, referred to as "Cartesian cut-stencil FD", is introduced to obtain the numerical solution of partial differential equations on any arbitrary irregular shaped domain. The 2^nd-order accurate two-dimensional Cartesian cut-stencil FD method utilizes a 5-point stencil and relies on the construction of a unique mapping of each physical stencil, rather than a cell, in any arbitrary domain to a generic uniform computational stencil. The treatment of boundary conditions and quantification of the solution accuracy using the local truncation error are discussed. Numerical solutions of the steady convection-diffusion equation on sample complex domains have been obtained and the results have been compared to exact solutions for manufactured partial differential equations (PDEs) and other numerical solutions.     

数学物理中的总体和局部场GL方法

为了求解物理化学生物材料和金融中的微分方程,提出了一种总体(Global)和局部(Local)场方法.微分方程的求解区域可以是有限域,无限域,或具曲面边界的部分无限域.其无限域包括有限有界不均匀介质区域.其不均匀介质区域被分划为若干子区域之和.在这含非均匀介质的无限区域,将微分方程的解显式地表示为在若干非均匀介质子区域上和局部子曲面的积分的递归和.把正反算的非线性关系递归地显式化.在无限均匀区域,微分方程的解析解被称为初始总体场.微分方程解的总体场相继地被各个非均匀介质子区域的局部散射场所修正.这种修正过程是一个子域接着另个子域逐步相继地进行的.一旦所有非均匀介质子区域被散射扫描和有限步更新过程全部完成后,微分方程的解就获得了.称其为总体和局部场的方法,简称为GL方法.GL方法完全地不同于有限元及有限差方法,GL方法直接地逐子域地组装逆矩阵而获得解.GL方法无需求解大型矩阵方程,它克服了有限元大型矩阵解的困难.用有限元及有限差方法求解无限域上的微分方程时,人为边界及其上的吸收边界条件是必需的和困难的,人为边界上的吸收边界条件的不精确的反射会降低解的精确度和毁坏反算过程.GL方法又克服了有限元和有限差方法的人为边界的困难.GL方法既不需要任何人为边界又不需要任何吸收边界条件就可以子域接子域逐步精确地求解无限域上的微分方程.有限元和有限差方法都仅仅是数值的方法,GL方法将解析解和数值方法相容地结合起来.提出和证明了三角的格林函数积分方程公式.证明了当子域的直经趋于零时,波动方程的GL方法的数值解收敛于精确解.GL方法解波动方程的误差估计也获得了.求解椭圆型,抛物线型,双曲线型方程的GL模拟计算结果显示出我们的GL方法具有准确,快速,稳定的许多优点.GL方法可以是有网,无网和半网算法.GL方法可广泛应用在三维电磁场,三维弹塑性力学场,地震波场,声波场,流场,量子场等方面.上述三维电磁场等应用领域的GL方法的软件已经由作者研制和发展了。     

A layer adaptive B‐spline collocation method for singularly perturbed one‐dimensional parabolic problem with a boundary turning point

In this article, we develop a parameter uniform numerical method for a class of singularly perturbed parabolic equations with a multiple boundary turning point on a rectangular domain. The coefficient of the first derivative with respect to x is given by the formula a₀(x, t)x^p, where a₀(x, t) ≥ α > 0 and the parameter p ∈ [1,∞) takes the arbitrary value. For small values of the parameter ε, the solution of this particular class of problem exhibits the parabolic boundary layer in a neighborhood of the boundary x = 0 of the domain. We use the implicit Euler method to discretize the temporal variable on uniform mesh and a B‐spline collocation method defined on piecewise uniform Shishkin mesh to discretize the spatial variable. Asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular component. These bounds are applied in the convergence analysis of the proposed scheme on Shishkin mesh. The resulting method is boundary layer resolving and has been shown almost second‐order accurate in space and first‐order accurate in time. It is also shown that the proposed method is uniformly convergent with respect to the singular perturbation parameter ε. Some numerical results are given to confirm the predicted theory and comparison of numerical results made with a scheme consisting of a standard upwind finite difference operator on a piecewise uniform Shishkin mesh. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1143–1164, 2011     

抛物型初边值问题的有限元与边界积分耦合的离散化及其误差分析

1．引言边界元方法是近二十几年来迅速发展起来的一类新的偏微分方程的数值方法．它的独特之处是将空间的维数降低一维,从而倍受工程技术人员的青睐,并在工程技术与计算数学领域得到越来越广泛的重视和研究．对椭圆型问题,边界元方法的理论与应用研究已取得丰硕成果;对发展型问题,近年来在理论方面的研究也已取得重要进展［6－11］．但边界元方法难以处理非均质问题,而有限元对各类问题及各种区域具有较好的适应性,将两者结合起来可充分发挥各自的优点．文山提出了一种抛物方程初边值问题的有限元与边界积分的耦合方法,其主要思想是…     

Domain Decomposition,Optimal Control of Systems Governed by Partial Differential Equations,and Synthesis of Feedback Laws

We present an iterative domain decomposition method for the optimal control of systems governed by linear partial differential equations. The equations can be of elliptic, parabolic, or hyperbolic type. The space region supporting the partial differential equations is decomposed and the original global optimal control problem is reduced to a sequence of similar local optimal control problems set on the subdomains. The local problems communicate through transmission conditions, which take the form of carefully chosen boundary conditions on the interfaces between the subdomains. This domain decomposition method can be combined with any suitable numerical procedure to solve the local optimal control problems. We remark that it offers a good potential for using feedback laws (synthesis) in the case of time-dependent partial differential equations. A test problem for the wave equation is solved using this combination of synthesis and domain decomposition methods. Numerical results are presented and discussed. Details on discretization and implementation can be found in Ref. 1.     

A Smooth Regularization of the Projection Formula for Constrained Parabolic Optimal Control Problems

We present a smooth, that is, differentiable regularization of the projection formula that occurs in constrained parabolic optimal control problems. We summarize the optimality conditions in function spaces for unconstrained and control-constrained problems subject to a class of parabolic partial differential equations. The optimality conditions are then given by coupled systems of parabolic PDEs. For constrained problems, a non-smooth projection operator occurs in the optimality conditions. For this projection operator, we present in detail a regularization method based on smoothed sign, minimum and maximum functions. For all three cases, that is, (1) the unconstrained problem, (2) the constrained problem including the projection, and (3) the regularized projection, we verify that the optimality conditions can be equivalently expressed by an elliptic boundary value problem in the space-time domain. For this problem and all three cases we discuss existence and uniqueness issues. Motivated by this elliptic problem, we use a simultaneous space-time discretization for numerical tests. Here, we show how a standard finite element software environment allows to solve the problem and, thus, to verify the applicability of this approach without much implementation effort. We present numerical results for an example problem.     

求解带非线性边界条件的反应扩散方程数值解的组合迭代方法

1　引　　言我们首先考虑如下抛物型方程ｕｔ-ＤΔｕ =ｆ(ｘ ,ｔ ,ｕ) (ｔ∈ ( 0 ,Ｔ],ｘ∈Ω ) ｕ/ ν+ βｕ =ｇ(ｘ ,ｔ ,ｕ) (ｔ∈ ( 0 ,Ｔ],ｘ∈ Ω )ｕ(ｘ ,0 ) =ψ(ｘ) (ｘ∈Ω )( 1 .1 )其中Ｔ为正常数 ,Ω 是ＲＰ 空间的有界区域 记ＱＴ=Ω × ( 0 ,Ｔ],ＳＴ= Ω × ( 0 ,Ｔ],假设在ＱＴ上Ｄ≡ｄ(ｘ ,ｔ) >0 ,在ＳＴ 上β≡β(ｘ ,ｔ)≥ 0 又设 ｆ(ｘ ,ｔ,ｕ) ,ｇ(ｘ ,ｔ,ｕ)为关于ｕ的非线性函数 ,且对ｘ ,ｔ各参数满足Ｈ¨ｏｌｄｅｒ连续条件 将 ( 1 .1 )离散化之后我们得到相应的有限差分系统 ,当 ｇ(ｘ ,ｔ,ｕ)为ｕ的线性…     

Analysis of a domain decomposition method for the nearly elastic wave equations based on mixed finite element methods

This paper concerns the numerical solution of the nearly elasticwave equations with the first-order absorbing boundary condition;these equations describe the motion of a nearly elastic solidin the frequency domain. Two mixed finite elements, the Johnson-Mercierelement and the Arnold-Douglas-Gupta element, are adapted andanalyzed for the problem. The resulting mixed finite elementequations are complex-valued and are neither Hermitian nor definite.As a result, most standard iterative methods fail to convergefor the systems. To solve the mixed finite element equations,a parallelizable domain decomposition iterative method is proposed.The convergence of the method is demonstrated and a rate ofconvergence of the form 1 - Ch is derived. These results arevalid for the case when the original domain is decomposed intosubdomains which consist of an individual element associatedwith the above two mixed finite elements.     

DISCRETE APPROXIMATIONS FOR SINGULARLY PERTURBED BOUNDARY VALUE PROBLEMS WITH PARABOLIC LAYERS,Ⅱ

l)ThisworkwassupportedbyNWOthroughgrantIBo7-3Go12.BOUNDAarv^LUEPRoBLEMFORELLIPTICEQUMIONwiTHMIXEDBOUNDAavCONDITION1.IntroductionInthispedwesketchavarietyofspecialmethodswhichareusedforconstructinge-unifornilyconvergelltschemes-WeshaJldemonstrateamethodwhichachieveshaprovedaccuracyforsolvingsingularlyperturbedb0undaryvalueproblemforeiliPicequatiouswithparabolicboundarylayers-InSecti0n4weshallintroduceanaturalclass,B,oftritefferenceschemes,inwhich(bytheabovementi0nedaP…