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

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

一阶双曲问题间断有限元的后验误差分析

一阶双曲问题的有限元后验误差估计至今没有得到很好的解决.本文对d维区域上一阶双曲问题的k次间断有限元逼近提出了一种新的后验误差分析方法, 进而建立了间断有限元解在DG范数下(强于L₂范数)基于误差余量型的后验误差估计. 数值计算验证了本文理论分析的有效性. 本文方法也适用于其他变分问题有限元逼近的后验误差分析.     

DG油田特高含水期层系组合优化方法研究

DG油田X区块经过多年开发,可采储量变少,储层非均质性严重,后续开采难以维持,研究DG油田层系优化组合,改变现状达到可持续开采的目的是重中之重.根据灰色关联及层系聚类法对DG油田进行合理,科学的开发层系划分.利用数理统计和油藏工程方法,结合数值模拟软件进行数值模拟,填补了层系组合的指标体系的相关空白,建立一整套层系组合优化的经济检验指标,形成一整套成熟的层系优化组合方法.研究表明,运用此方法在目标油田可以科学、快速的给出层系组合优化方案,同时组合的结果符合油田实际情况,可显著增加经济效益改善油田开发效果.     

紧致DG模和Gorenstein DG代数

证明同调有界的连通微分分次代数（简称为DG代数）上的紧致DG模的ampli-tude与基代数的amplitude的差恰为该DG模的投射维数.由此可得非平凡的正则DG代数是同调无界的.对正则DG代数A,若它的同调代数H（A）是分次Koszul代数,则证明H（A）有有限的整体维数;如果把条件减弱为A是Koszul DG代数,则给出了一个H（A）的整体维数为无限的例子.对一般的正则DG代数A,给出了其为Gorenstein DG代数的一些等价刻画.对同调有限维的连通DG代数A,证明由紧致对象全体构成的三角范畴Dc（A）和Dc（Aop）存在Auslander-Reiten三角当且仅当A和Aop都是Gorenstein DG代数.当A是非平凡的正则DG代数,且H（A）是局部有限维时,Dc（A）不存在Auslander-Reiten三角.对正则DG代数A,转而讨论了Auslander-Reiten三角在Dlbf（A）以及Dlbf（Aop）上的存在性.     

一个非平凡的Calabi-Yau DG代数

证明例1中的DG代数不仅是Koszul,同调光滑DG代数,而且还是一个Calabi-Yau DG代数.该例子说明一个Calabi-Yau DG代数的同调分次代数不一定具有Calabi-Yau性质,甚至可能不是同调光滑的;另外,该例子还说明一个Calabi-Yau DG代数忘掉微分后得到的分次代数不一定是分次Calabi-Yau代数.     

连通微分分次代数的整体维数

首次把有理同伦论中的同伦不变量-锥长度(cone length)引入到微分分次(简记为DG)同调代数中,定义了连通DG代数上DG模的锥长度.连通DG代数A的左(右)整体维数定义为所有DGA-模(Aop-模)的锥长度的上确界.在一些特殊情形下,发现连通.DG代数A的左(右)整体维数与H(A)的整体维数有着密切的关系.任意一个连通分次代数,如果将它视为微分为O的连通DG代数,其左(右)整体维数与其作为连通分次代数的整体维数是一致的.因此该定义是连通分次代数整体维数的一种推广形式.证明A的整体维数足三角范畴D(A)以及Dc(A)的维数的一个上界.当A是正则DG代数时,给出了A的左(右)整体维数的一个有限上界.     

Duality for cochain DG algebras

This paper develops a duality theory for connected cochain DG algebras,with particular emphasis on the non-commutative aspects.One of the main items is a dualizing DG module which induces a duality between the derived categories of DG left-modules and DG right-modules with finitely generated cohomology.As an application,it is proved that if the canonical module k=A/A≥1 has a semi-free resolution where the cohomological degree of the generators is bounded above,then the same is true for each DG module with finitely generated cohomology.     

DG polynomial algebras and their homological properties

In this paper, we introduce and study differential graded(DG for short) polynomial algebras. In brief, a DG polynomial algebra A is a connected cochain DG algebra such that its underlying graded algebra A~# is a polynomial algebra K[x_1, x_2,..., x_n] with |xi| = 1 for any i ∈ {1, 2,..., n}. We describe all possible differential structures on DG polynomial algebras, compute their DG automorphism groups, study their isomorphism problems, and show that they are all homologically smooth and Gorenstein DG algebras. Furthermore, it is proved that the DG polynomial algebra A is a Calabi-Yau DG algebra when its differential ?_A≠ 0 and the trivial DG polynomial algebra(A, 0) is Calabi-Yau if and only if n is an odd integer.     

基于改进的模糊综合评判方法的聚表剂性能评价

聚表剂凭借其良好的增黏性,抗盐性以及乳化特性成为陆相油田潜在的优势驱替剂.以华鼎Ⅰ、海博BI和DG聚表剂为研究对象,分析了影响聚表剂溶液黏度性质的单因素指标,基于改进的模糊综合评判方法对三种聚表剂的粘度性能进行优选;通过并联岩心模型驱油实验对优选结果进行验证.结果表明华鼎Ⅰ的粘度性能最优,DG聚表剂的粘度性能最差;并联岩心驱油实验中以改善剖面能力为验证指标与模糊评价结果相吻合;以采出程度为验证指标与模糊评价结果不吻合,原因在于不同聚表剂之间的洗油能力存在差异.综合来看,改进的模糊数学评判方法对于聚表剂的优选具有一定的参考意义.     

三角形的一条性质及其应用

在三角形中,有如下一条常用的性质:图1如图1,P为△ABC内任一点,射线AP、BP、CP分别交BC、CA、AB于点D、E、F,EF交AP于点G.则AP·DG AG·DP=AP·DG AD·PG=2.证明如图1所示,由面积关系可得AG PG=S△AEF S△PEF=S△AEF S△PAF·S△PAF S△PEF=EB PB·AC EC=S△EBC S△PBC·S△ABC S△EBC=S△ABC S△PBC=AD PD.于是AG·PD=AD·PG=(AP+PD)(AP-AG)=AP2+AP·PD-AP·AG-AG·PD=AP(AP+PD-AG)-AG·PD=AP·DG-AG·PD,即AP·DG=2·AG·PD.所以AP·DG AG·PD=2.同理AP·DG AD·PG=2.故AP·DG AG·PD=AP·DG AD·PG=2.注(1)此处的证明是联想到“A、G、P、D P交点为     

A posteriori error estimate for discontinuous Galerkin finite element method on polytopal mesh

In this work, we derive a posteriori error estimates for discontinuous Galerkin finite element method on polytopal mesh. We construct a reliable and efficient a posteriori error estimator on general polygonal or polyhedral meshes. An adaptive algorithm based on the error estimator and DG method is proposed to solve a variety of test problems. Numerical experiments are performed to illustrate the effectiveness of the algorithm.     

Mixed Discontinuous Galerkin Time-Stepping Method for Linear Parabolic Optimal Control Problems

In this paper, we discuss the mixed discontinuous Galerkin (DG) finite element approximation to linear parabolic optimal control problems. For the state variables and the co-state variables, the discontinuous finite element method is used for the time discretization and the Raviart-Thomas mixed finite element method is used for the space discretization. We do not discretize the space of admissible control but implicitly utilize the relation between co-state and control for the discretization of the control. We derive a priori error estimates for the lowest order mixed DG finite element approximation. Moveover, for the element of arbitrary order in space and time, we derive a posteriori $L^2(0, T ;L^2(Ω))$ error estimates for the scalar functions, assuming that only the underlying mesh is static. Finally, we present an example to confirm the theoretical result on a priori error estimates.     

Analysis of the parareal approach based on discontinuous Galerkin method for time-dependent Stokes equations

This paper analyzes a parareal approach based on discontinuous Galerkin (DG) method for the time-dependent Stokes equations. A class of primal discontinuous Galerkin methods, namely variations of interior penalty methods, are adopted for the spatial discretization in the parareal algorithm (we call it parareal DG algorithm). We study three discontinuous Galerkin methods for the time-dependent Stokes equations, and the optimal continuous in time error estimates for the velocities and pressure are derived. Based on these error estimates, the proposed parareal DG algorithm is proved to be unconditionally stable and bounded by the error of discontinuous Galerkin discretization after a finite number of iterations. Finally, some numerical experiments are conducted which confirm our theoretical results, meanwhile, the efficiency of the parareal DG algorithm can be seen through a parallel experiment.     

A posteriori error estimates of discontinuous Galerkin methods for non-standard Volterra integro-differential equations

^** Email: jingtang{at}lsec.cc.ac.cn^*** Email: hermann{at}math.mun.ca In this paper we establish a posteriori error estimates forthe discontinuous Galerkin (DG) method applied to linear, semilinearand non-standard (non-linear) Volterra integro-differentialequations. We also present an analysis of the DG method withquadrature for the memory term. Numerical experiments basedon three integro-differential equations are used to illustratevarious aspects of the error analysis.     

Accuracy Enhancement of Discontinuous Galerkin Method for Hyperbolic Systems

We study the enhancement of accuracy, by means of the convolution post-processing technique, for discontinuous Galerkin(DG) approximations to hyperbolic problems. Previous investigations have focused on the superconvergence obtained by this technique for elliptic, time-dependent hyperbolic and convection-diffusion problems. In this paper, we demonstrate that it is possible to extend this post-processing technique to the hyperbolic problems written as the Friedrichs' systems by using an upwind-like DG method. We prove that the $L_2$-error of the DG solution is of order $k+1/2$, and further the post-processed DG solution is of order $2k+1$ if $Q_k$-polynomials are used. The key element of our analysis is to derive the $(2k+1)$-order negative norm error estimate. Numerical experiments are provided to illustrate the theoretical analysis.     

Numerical approximation of the Oldroyd‐B model by the weighted least‐squares/discontinuous Galerkin method

We consider a numerical method for the Oldroyd‐B model of viscoelastic fluid flows by a combination of the weighted least‐squares (WLS) method and the discontinuous Galerkin (DG) finite element method. The constitutive equation is decoupled from the momentum and continuity equations, and the approximate solution is computed iteratively by solving the Stokes problem and a linearized constitutive equation using WLS and DG, respectively. An a priori error estimate for the WLS/DG method is derived and numerical results supporting the estimate are presented. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

A POSTERIORI ERROR ESTIMATES FOR A MODIFIED WEAK GALERKIN FINITE ELEMENT APPROXIMATION OF SECOND ORDER ELLIPTIC PROBLEMS WITH DG NORM

In this paper,we derive a residual based a posteriori error estimator for a modified weak Galerkin formulation of second order elliptic problems.We prove that the error estimator used for interior penalty discontinuous Galerkin methods still gives both upper and lower bounds for the modified weak Galerkin method,though they have essentially different bilinear forms.More precisely,we prove its reliability and efficiency for the actual error measured in the standard DG norm.We further provide an improved a priori error estimate under minimal regularity assumptions on the exact solution.Numerical results are presented to verify the theoretical analysis.     

A posteriori error estimates for discontinuous Galerkin method to the elasticity problem

This work concerns with the discontinuous Galerkin (DG) method for the time‐dependent linear elasticity problem. We derive the a posteriori error bounds for semidiscrete and fully discrete problems, by making use of the stationary elasticity reconstruction technique which allows to estimate the error for time‐dependent problem through the error estimation of the associated stationary elasticity problem. For fully discrete scheme, we make use of the backward‐Euler scheme and an appropriate space‐time reconstruction. The technique here can be applicable for a variety of DG methods as well.     

Adaptive discontinuous Galerkin methods in multiwavelets bases

We use a multiwavelet basis with the Discontinuous Galerkin (DG) method to produce a multi-scale DG method. We apply this Multiwavelet DG method to convection and convection-diffusion problems in multiple dimensions. Merging the DG method with multiwavelets allows the adaptivity in the DG method to be resolved through manipulation of multiwavelet coefficients rather than grid manipulation. Additionally, the Multiwavelet DG method is tested on non-linear equations in one dimension and on the cubed sphere.     

Interior penalty discontinuous Galerkin method for Maxwell's equations: Energy norm error estimates

We develop the symmetric interior penalty discontinuous Galerkin (DG) method for the time-dependent Maxwell equations in second-order form. We derive optimal a priori error estimates in the energy norm for smooth solutions. We also consider the case of low-regularity solutions that have singularities in space.