电报方程的时间间断时空有限元方法

为同时高精度逼近速度和位移,利用时间间断的时空有限元与降阶的思想,对一类电报方程的初边值问题建立一种时间间断时空有限元格式.利用有限差分方法与有限元方法相结合的技巧,证明了格式的稳定性和收敛性,得到了速度的L∞(L2)模和位移的L∞(H1)模最优误差估计.最后用数值算例验证了理论分析结果和所提算法的有效性.     

发展型方程的时间间断时空有限元方法

时空有限元方法通过统一时间和空间变量,克服了传统有限元方法对时间作差分离散引起的时间上的低精度,不但具有时、空高精度,而且在无结构网格上耗散特性好、无条件稳定,成为解决时间依赖问题的有效方法.本文利用抛物问题给出时间允许间断而空间连续的时空有限元方法的基本概念和过程,给出抛物型方程、积分-微分方程、双曲方程、Sobolev方程和其他高阶方程的算例,验证方法的精度和稳定性,并综合评价时间间断时空有限元方法目前的发展现状和应用前景.     

Sobolev 方程的$H^1$-Galerkin混合有限元方法

对Sobolev方程采用H1-Galerkin混合有限元方法进行数值模拟．给出了一维空间中该方法的半离散和全离散格式及其最优误差估计;并将该方法推广到二维和三维空间．与H1-Galerkin有限元方法相比,该方法不仅降低了对有限元空间的连续性要求;而且与传统的混合有限元方法具有相同的收敛阶,但其有限元空间的选取却不需要满足LBB相容条件．数值例子将进一步说明该方法的可行性与有效性．     

Sobolev型方程有限元的后处理方法

1 引言 考虑下述Sobolev型方程的混合问题 (a) (b) (c) 其中Ω为R~2中具有边界的矩形域,a,b,f,u_o。为适当光滑且有界的已知函数,a(x,t)有正下界a_*. Sobolev型方程是重要的数学物理方程之一,文[1]导出了问题(1.1)标准有限元方法的最优L_2(2≤P<∞)估计.本文研究矩形剖分上的双k次有限元方法,用插值算子对近似解进行后处理,仅增加极少工作量,使整体收敛性提高一阶.本文证明了误差及任意阶时间导 进行后处理,仅增加极少工作量,使整体收敛性提高一阶.本文证明了误差及任意阶时间导数的H~1,W~(1,∞),L_p和L_∞的超收敛估计.若采用文[2]的预处理方法构造最优剖分,可将本文结果推广到一般区域(仍超收敛1/2阶).这样,采用低次有限元可获得高阶精度,从而大大节省了计算量.     

热传导方程的间断Galerkin数值解法

本文对具有周期边界的热传导方程采用间断Galerkin(DG)方法给出数值求解方法,并利用傅里叶分析,对数值解进行L∞-误差估计,以一次分段多项式为例,得到半离散格式的误差估计.     

Sobolev方程的连续时空有限元方法

本文研究Sobolev方程的连续时空有限元方法.首先建立Sobolev方程的连续时空有限元格式,然后证明了解的存在唯一性和稳定性并给出连续时空有限元解各种范数下的误差估计.最后给出数值算例来验证理论分析的正确性,并进一步说明本文所建立的格式关于时间可以得到比传统有限元方法更高的精度.     

半线性Sobolev方程的H~1-Galerkin混合有限元方法

利用H~1-Galerkin混合有限元方法研究了一维半线性Sobolev方程,得到了半离散解的最优阶误差估计,优点是不需验证LBB相容性条件.     

流固耦合方程组间断Galerkin方法的探索

主要通过对复杂接触表面问题以及流固耦合方程组中边界间断问题的分析,探讨其间断Galerkin方法的有限元计算.保留有限元线性离散的计算优势,有效地弱化了边界间断对流场中速度的影响,得到流固耦合方程组的空间半离散有限元格式,为数值计算提供了有力的理论支撑.     

一种推广的对流扩散方程的局部化间断Galerkin方法

研究求解一种产生于径向渗流问题的推广的对流扩散方程的局部化间断Galerkin方法,对一般非线性情形证明了方法的L^2稳定性;对线性情形证明了,当方法取有限元空间为κ次多项式空间时,数值解逼近的L^∞（0,T;L^2)模的误差阶为κ。     

对流反应扩散方程的稳定化时间间断时空有限元解的误差估计

在流线迎风Petrov-Galerkin(SUPG)稳定化有限元数值格式的基础上,结合时间方向的变分离散,构造对流反应扩散方程的稳定化时间间断时空有限元格式.该类格式在工程上有一些数值模拟应用,但相关文献没有看到类似数值格式的理论证明.本文以Radau点为节点,构造时间方向的Lagrange插值多项式,证明了稳定化有限...     

On discontinuous Galerkin finite element method for singularly perturbed delay differential equations

A nonsymmetric discontinuous Galerkin FEM with interior penalties has been applied to one-dimensional singularly perturbed problem with a constant negative shift. Using higher order polynomials on Shishkin-type layer-adapted meshes, a robust convergence has been proved in the corresponding energy norm. Numerical experiments support theoretical findings.     

A space–time discontinuous Galerkin method for linear convection-dominated Sobolev equations

This article presents a space–time discontinuous Galerkin (DG) finite element method for linear convection-dominated Sobolev equations. The finite element method has basis functions that are continuous in space and discontinuous in time, and variable spatial meshes and time steps are allowed. In the discrete intervals of time, using properties of the Radau quadrature rule, eliminates the restriction to space–time meshes of convectional space–time Galerkin methods. The existence and uniqueness of the approximate solution are proved. An optimal priori error estimate in L^∞(H¹) is derived. Numerical experiments are presented to confirm theoretical results.     

A discontinuous Galerkin finite element method for time dependent partial differential equations with higher order derivatives

In this paper, we develop a new discontinuous Galerkin (DG) finite element method for solving time dependent partial differential equations (PDEs) with higher order spatial derivatives. Unlike the traditional local discontinuous Galerkin (LDG) method, the method in this paper can be applied without introducing any auxiliary variables or rewriting the original equation into a larger system. Stability is ensured by a careful choice of interface numerical fluxes. The method can be designed for quite general nonlinear PDEs and we prove stability and give error estimates for a few representative classes of PDEs up to fifth order. Numerical examples show that our scheme attains the optimal -th order of accuracy when using piecewise -th degree polynomials, under the condition that is greater than or equal to the order of the equation.

     

Connections between discontinuous Galerkin and nonconforming finite element methods for the Stokes equations

We study a discontinuous Galerkin finite element method (DGFEM) for the Stokes equations with a weak stabilization of the viscous term. We prove that, as the stabilization parameter γ tends to infinity, the solution converges at speed γ^?1 to the solution of some stable and well‐known nonconforming finite element methods (NCFEM) for the Stokes equations. In addition, we show that an a posteriori error estimator for the DGFEM‐solution based on the reconstruction of a locally conservative H(div, Ω)‐tensor tends at the same speed to a classical a posteriori error estimator for the NCFEM‐solution. These results can be used to affirm the robustness of the DGFEM‐method and also underline the close relationship between the two approaches. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A weak Galerkin finite element method for the stokes equations

This paper introduces a weak Galerkin (WG) finite element method for the Stokes equations in the primal velocity-pressure formulation. This WG method is equipped with stable finite elements consisting of usual polynomials of degree k≥1 for the velocity and polynomials of degree k?1 for the pressure, both are discontinuous. The velocity element is enhanced by polynomials of degree k?1 on the interface of the finite element partition. All the finite element functions are discontinuous for which the usual gradient and divergence operators are implemented as distributions in properly-defined spaces. Optimal-order error estimates are established for the corresponding numerical approximation in various norms. It must be emphasized that the WG finite element method is designed on finite element partitions consisting of arbitrary shape of polygons or polyhedra which are shape regular.     

A weak Galerkin finite element method for the Oseen equations

In this paper, a weak Galerkin finite element method for the Oseen equations of incompressible fluid flow is proposed and investigated. This method is based on weak gradient and divergence operators which are designed for the finite element discontinuous functions. Moreover, by choosing the usual polynomials of degree i ≥ 1 for the velocity and polynomials of degree i ? 1 for the pressure and enhancing the polynomials of degree i ? 1 on the interface of a finite element partition for the velocity, this new method has a lot of attractive computational features: more general finite element partitions of arbitrary polygons or polyhedra with certain shape regularity, fewer degrees of freedom and parameter free. Stability and error estimates of optimal order are obtained by defining a weak convection term. Finally, a series of numerical experiments are given to show that this method has good stability and accuracy for the Oseen problem.     

Space–time discontinuous Galerkin finite element method for shallow water flows

A space–time discontinuous Galerkin (DG) finite element method is presented for the shallow water equations over varying bottom topography. The method results in nonlinear equations per element, which are solved locally by establishing the element communication with a numerical HLLC flux. To deal with spurious oscillations around discontinuities, we employ a dissipation operator only around discontinuities using Krivodonova's discontinuity detector. The numerical scheme is verified by comparing numerical and exact solutions, and validated against a laboratory experiment involving flow through a contraction. We conclude that the method is second order accurate in both space and time for linear polynomials.     

Extrapolation discontinuous Galerkin method for ultraparabolic equations

Ultraparabolic equations arise from the characterization of the performance index of stochastic optimal control relative to ultradiffusion processes; they evidence multiple temporal variables and may be regarded as parabolic along characteristic directions. We consider theoretical and approximation aspects of a temporally order and step size adaptive extrapolation discontinuous Galerkin method coupled with a spatial Lagrange second-order finite element approximation for a prototype ultraparabolic problem. As an application, we value a so-called Asian option from mathematical finance.