A Vertex-Centered and Positivity-Preserving Finite Volume Scheme for Two-Dimensional Three-Temperature Radiation Diffusion Equations on General Polygonal Meshes

Two-dimensional three-temperature (2-D 3-T) radiation diffusion equations are widely used to approximately describe the evolution of radiation energy within a multimaterial system and explain the exchange of energy among electrons, ions and photons. In this paper, we suggest a new positivity-preserving finite volume scheme for 2-D 3-T radiation diffusion equations on general polygonal meshes. The vertex unknowns are treated as primary ones for which the finite volume equations are constructed. The edge-midpoint and cell-centered unknowns are used as auxiliary ones and interpolated by the primary unknowns, which makes the final scheme a pure vertex-centered one. By comparison, most existing positivity-preserving finite volume schemes are cell-centered and based on the convex decomposition of the co-normal. Here, the co-normal decomposition is not convex in general, leading to a fixed stencil of the flux approximation and avoiding a certain search algorithm on complex grids. Moreover, the new scheme effectively alleviates the numerical heat-barrier issue suffered by most existing cell-centered or hybrid schemes in solving strongly nonlinear radiation diffusion equations. Numerical experiments demonstrate the second-order accuracy and the positivity of the solution on various distorted grids. For the problem without analytic solution, the contours of the numerical solutions obtained by our scheme on distorted meshes accord with those on smooth quadrilateral meshes.     

Alternating Direction Finite Volume Element Methods for Three-Dimensional Parabolic Equations

<正>This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes,one is analogous to Douglas finite difference scheme with second-order splitting error,the other two schemes have third-order splitting error,and the last one is an extended LOD scheme.The L~2 norm and H~1 semi-norm error estimates are obtained for the first scheme and second one,respectively.Finally,two numerical examples are provided to illustrate the efficiency and accuracy of the methods.     

双曲型积分-微分方程的有限体积元方法

本文研究了双曲型积分—微分方程的有限体积元方法，利用基于有限体积元的Ritz—Volterra投影的逼近性质，得到了半离散有限体积元解的最优阶L2，H^1，L∞和W^1,∞模误差估计．     

Two-grid Methods for Finite Volume Element Approximations of Nonlinear Sobolev Equations

In this article, two-grid methods are studied for solving nonlinear Sobolev equation using the finite volume element method. The methods are based on one coarse grid space and one fine grid space. The nonsymmetric and nonlinear iterations are only executed on the coarse grid (with grid size H), and the fine grid solution (with grid size h) can be obtained in a single symmetric and linear step. The optimal H¹ error estimates are presented for the proposed methods, which show that the two-grid methods achieve optimal approximation as long as the mesh sizes satisfy h = 𝒪(H³|ln H|). As a result, solving such a large class of nonlinear Sobolev equations will not be much more difficult than solving one linearized equation.     

Meshfree Finite Volume Element Method for Constrained Optimal Control Problem Governed by Random Convection Diffusion Equations

In this paper, we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients. There are two contributions of this paper. Firstly, we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space, which is competitive for high-dimensional random inputs. Secondly, the a priori error estimates are derived for the state,the co-state and the control variables. Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.     

拟线性对流占优扩散方程的扩展特征混合有限元方法数值模拟

讨论了拟线性对流占优扩散问题的数值模拟.对对流部分采用特征线格式进行离散,以消除流动锋线前沿的数值弥散现象,保证格式的稳定性;而对扩散部分采用扩展混合有限元方法,同时逼近未知函数,未知函数的梯度及伴随向量函数.理论分析和数值算例表明,此方法是稳定的,具有最优L2逼近精度.     

对流扩散方程的混合时间间断时空有限元方法

构造并分析二阶对流扩散方程的混合时间间断时空有限元格式．利用混合有限元方法将二阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散低阶方程．证明数值解的稳定性、存在唯一性和收敛性．最后通过数值结果验证该算法的有效性和可行性．     

定常Navier—Stokes方程的Petrov—Galerkin方法

研究了定常Navier－Stokes方程的四种Petrov－Galerkin有限元方法：PG1，PG2，SD和GLS．它们都是稳定的，避免了经典混合方法中必要的Babuska－Brezzi条件．给出了各种方法有限元解的存在性、唯一性和唯一解的误差估计．     

Mortar Upwind Finite Volume Element Method with Crouzeix-Raviart Element for Parabolic Convection Diffusion Problems

In this paper,we study the semi-discrete mortar upwind finite volume element method with the Crouzeix-Raviart element for the parabolic convection diffusion problems. It is proved that the semi-discrete mortar upwind finite volume element approximations derived are convergent in the H~1-and L~2-norms.     

非定常线性化Navier-Stokes方程的非协调流线扩散有限元法分析

对非定常线性化Navier-Stokes方程提出了非协调流线扩散有限元方法．用向后Euler格式离散时间,用流线扩散法处理扩散项带来的非稳定性．速度采用不连续的分片线性逼近,压力采用分片常数逼近．得到了离散解的存在唯一性以及在一定范数意义下离散解的稳定性和误差估计．     

Constraint-Preserving Hybrid Finite Element Methods for Maxwell’s Equations

Foundations of Computational Mathematics - Maxwell’s equations describe the evolution of electromagnetic fields, together with constraints on the divergence of the magnetic and electric flux...     

Cubic Finite Volume Methods for Second Order Elliptic Equations with Variable

In this paper, we present a finite volume framework for second order elliptic equations with variable coefficients based on cubic Hermite element. We prove the optimal H^1 norm error estimates. A numerical example is given at the end to show the feasibility of the method.     

一维定常型对流占优扩散方程的一类迎风有限体积格式

本文针对一维定常型对流占优扩散方程提出了一类迎风有限体积格式 .该格式对对流项具有二阶精度 ,对扩散项保持一阶精度 ,符合对流占优扩散问题强对流、弱扩散的特点 .     

Sobolev-Volterra投影与积分微分方程有限元数值分析

本文提出一类称之为Sobolev-Volterra投影的有限元投影,研究了有关性质并将之应用于伪抛物型积分微分方程有限元方法、伪双曲型积分微分方程有限元方法以及三维伪双曲型积分微分方程交替方向有限元方法的数值分析.     

A Nonconforming Finite Element Method of Streamline Diffusion Type for the Incompressible Navier-Stokes Equations

A nonconforming finite element method of streamline diffusion type for solving the stationary and incompressible Navier-Stokes equation is considered. Velocity field and pressure field are approximated by piecewise linear and piecewise constant functions, respectively. The existence of solutions of the discrete problem and the strong convergence of a subsequence of discrete solutions are established. Error estimates are presented for the uniqueness case.     

二维空间中分数阶扩散方程的变网格有限元方法

针对二维空间分数阶偏微分方程,给出了一个变网格全离散有限元格式,并得到了相应最优误差估计.其主要思想是对空间变量采用有限元离散,对时间交量采用差分,但不同时刻的有限元网格可以不同.这对于没计相应的自适应算法是十分有益的.     

二次变系数椭圆方程的三次有限体积元方法

In this paper, we present a finite volume framework for second order elliptic equations with variable coefficients based on cubic Hermite element. We prove the optimal H1 norm error estimates. A numerical example is given at the end to show the feasibility of the method.     

Streamline Diffusion Methods for Operator Equations

To solve a class of operator equations numerically, some general streamline diffusion methods with satisfactory convergence properties are presented in this paper. It is proved that the approximation accuracy is only half a power of $h$, the mesh size, from being optimal when these methods are applied to mixed problems and convection-diffusion problems.     

非线性抛物型方程组的二次有限体积元方法及其误差估计

讨论基于三角形网格的二维非线性抛物型方程组的有限体积元方法,其中试探函数空间为二次Lagrange元,检验函数空间为分片常数函数空间,对问题的全离散格式证明了最优的能量模误差估计。最后给出一个相关数值算例以验证格式的有效性。     

Superconvergence of H 1-Galerkin Mixed Finite Element Methods for Second-Order Elliptic Equations

We derive superconvergence result for H ¹-Galerkin mixed finite element method for second-order elliptic equations over rectangular partitions. Compared to standard mixed finite element procedure, the method is not subject to the Ladyzhenskaya–Bab?ska–Brezzi (LBB) condition and the approximating finite element spaces are allowed to be of different polynomial degrees. Superconvergence estimate of order 𝒪(h ^k+3), where k ≥ 1 is the order of the approximating polynomials employed in the Raviart–Thomas elements, is established for the flux via a postprocessing technique.