基于Voronoi网格的扩散方程差分格式

在Voronoi网格上利用一种基于回路积分法的有限体积法构造扩散方程的的差分格式.在这种特殊的网格上离散扩散方程比通常在四边形网格上离散的格式要简单,不会引进角点未知量,提高了对网格边上的流的离散精度,及差分格式整体精度.这种Voronoi网格上的扩散计算也可以与单元中心流体力学计算耦合.数值算例表明这种格式比四边形网格上的格式精度高,且能更好的应对网格扭曲情形.     

基于二阶格式的有限体积保正格式

从二阶线性格式出发，通过对法向通量进行重构，得到非线性两点通量，获得四面体网格上的单元中心型有限体积保正格式.该格式适用于求解间断和各向异性扩散系数问题.无需假设辅助未知量非负，避免了辅助未知量计算出负时"遇负置零"的人为处理方式；并且证明该格式在每个非线性Picard迭代步具有强保正性，即当源项和边界条件非负时，线性化格式的非平凡解是严格大于零的.数值算例验证该格式具有二阶收敛性且是保正的.     

基于WENO重构的高阶移动网格动理学格式

提出一种基于WENO重构的高阶(至少三阶)移动网格动理学格式.利用流体力学方程的积分形式得到移动网格上离散格式,再利用自适应移动网格方法移动网格,进而得到网格速度,利用WENO重构得到高阶插值多项式,最后使用时间方向上精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高精度、高分辨率的特点.     

Non-negative mixed finite element formulations for a tensorial diffusion equation

We consider the tensorial diffusion equation, and address the discrete maximum–minimum principle of mixed finite element formulations. In particular, we address non-negative solutions (which is a special case of the maximum–minimum principle) of mixed finite element formulations. It is well-known that the classical finite element formulations (like the single-field Galerkin formulation, and Raviart–Thomas, variational multiscale, and Galerkin/least-squares mixed formulations) do not produce non-negative solutions (that is, they do not satisfy the discrete maximum–minimum principle) on arbitrary meshes and for strongly anisotropic diffusivity coefficients.     

基于移动网格的熵稳定格式

提出一种基于移动网格的熵稳定格式求解双曲型守恒律方程.该方法利用等分布原理得到新的网格分布,基于守恒型插值公式计算新的网格上的物理量,使用熵稳定数值通量和三阶强稳定Runge-Kutta时间推进方法得到下一时刻的数值解.数值算例表明该格式不仅能有效提高解在间断处的分辨率,而且能消除可能产生的伪振荡.     

Discrete Maximum Principle for Poisson Equation with Mixed Boundary Conditions Solved by $hp$-FEM

We present a proof of the discrete maximum principle (DMP) for the 1D Poisson equation $−u''=f$ equipped with mixed Dirichlet-Neumann boundary conditions. The problem is discretized using finite elements of arbitrary lengths and polynomial degrees ($hp$-FEM). We show that the DMP holds on all meshes with no limitations to the sizes and polynomial degrees of the elements.     

An improved monotone finite volume scheme for diffusion equation on polygonal meshes

We construct a new nonlinear monotone finite volume scheme for diffusion equation on polygonal meshes. The new scheme uses the cell-edge unknowns instead of cell-vertex unknowns as the auxiliary unknowns in order to improve the accuracy of monotone scheme. Our scheme is locally conservative and has only cell-centered unknowns. Numerical results are presented to show how our scheme works for preserving positivity on various distorted meshes. Specially, numerical results show that the new scheme is robust, and more accurate than the existing monotone scheme on some kinds of meshes.     

九点格式中网格边上切向流的计算——节点自适应加权格式

针对网格扭曲的不同情形,直接考虑网格边上切向流的离散.基于扩散方程法向流连续的条件,给出离散法向流的构造,导出扭曲网格上九点计算格式中网格边上离散切向流的表达式,从而推导出加权系数的计算公式,适应于各种扭曲的网格.数值结果表明,与九点格式中节点量简单加权的方法相比,基于网格边离散切向流的节点自适应加权九点格式的精度有明显改进,迭代求解次数减少,计算效率明显提高.     

大长宽比扭曲网格上辅助网格差分方法

讨论抛物型方程的离散差分格式的构造,对九点差分格式进行了适用范围的讨论,并在此基础上提出辅助网格差分方法,用于处理因网格长宽比大且扭曲较大的网格引起的计算精度与计算效率降低的问题,该方法从守恒方程出发,将九点差分格式应用于按某种合适的方式进行重分之后的网格上,减少由于网格正则性差以及网格节点上的物理量采用周围网格量的加权平均等原因所引起的计算误差,得到一个新的但其解仍然逼近原来网格上的物理量的方程组.所构造的方法便于实施,且更适合于对实际物理模型的模拟,能比较好地适应流体大变形导致的网格扭曲,数值试验表明它有较好的数值精度和稳定性.     

The discrete variational derivative method based on discrete differential forms

As is well known, for PDEs that enjoy a conservation or dissipation property, numerical schemes that inherit this property are often advantageous in that the schemes are fairly stable and give qualitatively better numerical solutions in practice. Lately, Furihata and Matsuo have developed the so-called “discrete variational derivative method” that automatically constructs energy preserving or dissipative finite difference schemes. Although this method was originally developed on uniform meshes, the use of non-uniform meshes is of importance for multi-dimensional problems. On the other hand, the theories of discrete differential forms have received much attention recently. These theories provide a discrete analogue of the vector calculus on general meshes. In this paper, we show that the discrete variational derivative method and the discrete differential forms by Bochev and Hyman can be combined. Applications to the Cahn–Hilliard equation and the Klein–Gordon equation on triangular meshes are provided as demonstrations. We also show that the schemes for these equations are H¹-stable under some assumptions. In particular, one for the nonlinear Klein–Gordon equation is obtained by combination of the energy conservation property and the discrete Poincaré inequality, which are the temporal and spacial structures that are preserved by the above methods.     

Analysis of the monotonicity conditions in the mimetic finite difference method for elliptic problems

The maximum principle is one of the most important properties of solutions of partial differential equations. Its numerical analog, the discrete maximum principle (DMP), is one of the most difficult properties to achieve in numerical methods, especially when the computational mesh is distorted to adapt and conform to the physical domain or the problem coefficients are highly heterogeneous and anisotropic. Violation of the DMP may lead to numerical instabilities such as oscillations and to unphysical solutions such as heat flow from a cold material to a hot one. In this work, we investigate sufficient conditions to ensure the monotonicity of the mimetic finite difference (MFD) method on two- and three-dimensional meshes. These conditions result in a set of general inequalities for the elements of the mass matrix of every mesh element. Efficient solutions are devised for meshes consisting of simplexes, parallelograms and parallelepipeds, and orthogonal locally refined elements as those used in the AMR methodology. On simplicial meshes, it turns out that the MFD method coincides with the mixed-hybrid finite element methods based on the low-order Raviart–Thomas vector space. Thus, in this case we recover the well-established conventional angle conditions of such approximations. Instead, in the other cases a suitable design of the MFD method allows us to formulate a monotone discretization for which the existence of a DMP can be theoretically proved. Moreover, on meshes of parallelograms we establish a connection with a similar monotonicity condition proposed for the Multi-Point Flux Approximation (MPFA) methods. Numerical experiments confirm the effectiveness of the considered monotonicity conditions.     

辐射流体力学的分子动理学自适应加密方法(英文)

通过构造新的平衡分布函数和结合分区自适应网格加密方法,对不带扩散项的平衡辐射流体力学方程,构造二阶的分子动理学BGK-AMR格式.一方面在关心的计算区域中局部加密计算网格,提高计算精度的同时大大节省了计算网格数量和计算时间;另一方面,不同于已有的参数强耦合平衡分布函数,新构造的平衡分布函数中各参数不相互依赖,简化了辐射流体力学分子动理学格式的计算.一维和二维的数值算例显示了格式的性能.     

二维扩散方程的局部二阶线性保极值节点计算方法及应用

对一般四边形网格设计一种优化的节点控制体, 并构造了一种扩散方程的保极值二阶收敛的局部线性节点计算格式(优化控制体节点格式, VOC格式)。在网格不出现异常节点的情况下, 证明VOC格式是保极值、线性精确和二阶收敛的。而且在均匀的矩形网格上, 修正的逆距离加权格式与VOC格式等价, 从而对间断系数问题也是局部二阶收敛的。VOC格式可以用于单元中心型线性扩散格式和保正格式的节点值计算。数值算例表明对扭曲网格上的间断系数问题, VOC格式是二阶收敛的。采用VOC格式计算节点值的线性九点格式具有线性精确性和二阶收敛性, 采用VOC格式的保正格式也具有二阶收敛性。     

适应强扭曲网格的扩散方程时空二阶精度计算

以全局支撑算子方法为基础,通过引入面通量,构造了具有局部模板点的时空二阶精度格式。对于大变形扭曲网格,格式采用法向修正技术和合理的单元角体积计算方法,可以保持通量的精确性。算例表明该格式在非凸网格上能够精确获得线性解; 在非光滑网格上可以达到时空二阶精度; 能够较好地保持对称性; 并适合于三维非结构网格上的求解。     

Monotone finite volume schemes for diffusion equations on unstructured triangular and shape-regular polygonal meshes

We consider a non-linear finite volume (FV) scheme for stationary diffusion equation. We prove that the scheme is monotone, i.e. it preserves positivity of analytical solutions on arbitrary triangular meshes for strongly anisotropic and heterogeneous full tensor coefficients. The scheme is extended to regular star-shaped polygonal meshes and isotropic heterogeneous coefficients.     

Arbitrary Lagrangian–Eulerian finite-element method for computation of two-phase flows with soluble surfactants

A finite-element scheme based on a coupled arbitrary Lagrangian–Eulerian and Lagrangian approach is developed for the computation of interface flows with soluble surfactants. The numerical scheme is designed to solve the time-dependent Navier–Stokes equations and an evolution equation for the surfactant concentration in the bulk phase, and simultaneously, an evolution equation for the surfactant concentration on the interface. Second-order isoparametric finite elements on moving meshes and second-order isoparametric surface finite elements are used to solve these equations. The interface-resolved moving meshes allow the accurate incorporation of surface forces, Marangoni forces and jumps in the material parameters. The lower-dimensional finite-element meshes for solving the surface evolution equation are part of the interface-resolved moving meshes. The numerical scheme is validated for problems with known analytical solutions. A number of computations to study the influence of the surfactants in 3D-axisymmetric rising bubbles have been performed. The proposed scheme shows excellent conservation of fluid mass and of the total mass of the surfactant.     

反应流体的移动网格动理学格式

在移动网格上构造一种反应流的动理学格式.首先利用BGK模型推导含化学反应的流体力学方程组,并利用其积分形式构造移动网格上离散格式,再利用自适应移动网格方法得到网格速度,最后利用时间精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.一维与二维的数值实验表明这种格式同时具有高精度、高分辨率的特点.     

Geometric formulations and variational integrators of discrete autonomous Birkhoff systems

The variational integrators of autonomous Birkhoff systems are obtained by the discrete variational principle. The geometric structure of the discrete autonomous Birkhoff system is formulated. The discretization of mathematical pendulum shows that the discrete variational method is as effective as symplectic scheme for the autonomous Birkhoff systems.     

基于非结构四边形网格的WENO格式

基于非结构四边形网格发展求解双曲守恒律的三阶加权基本无振荡（WENO）格式.针对任意非结构四边形网格选取重构模板，并给出基于线性多项式的三阶线性重构.但对于一般的非结构四边形网格，会出现非常大的线性权和负权，使得非线性重构的WENO格式对光滑问题也不稳定.本文给出一个处理非常大的线性权的优化重构方法，对优化后得到的负线性权采用分裂方法进行处理.对于非线性权，提出一种考虑局部网格和物理量间断的新光滑度量因子.采用优化重构方法和新的非线性权，当前的三阶WENO格式在质量很差的网格上也具有很好的稳定性.理论的三阶精度在数值精度测试算例中得到验证，同时一范数和无穷范数的误差绝对值不依赖于网格质量；具有强间断的数值结果证明了当前格式的有效性.     

A central WENO scheme for solving hyperbolic conservation laws on non-uniform meshes

This paper proposes a new WENO procedure to compute multi-scale problems with embedded discontinuities, on non-uniform meshes.In a one-dimensional context, the WENO procedure is first defined on a five-points stencil and designed to be fifth-order accurate in regions of smoothness. To this end, we define a finite-volume discretization in which we consider the cell averages of the variable as the discrete unknowns. The reconstruction of their point-values is then ensured by a unique fifth-order polynomial. This optimum polynomial is considered as a symmetric and convex combination, by ideal weights, of four quadratic polynomials.The symmetric nature of the resulting interpolation has an important consequence: the choice of ideal weights has no influence on the accuracy of the discretization. This advantage enables to formulate the interpolation for non-uniform meshes. Following the methodology of the classic WENO procedure, non-oscillatory weights are then calculated from the ideal weights.We adapt this procedure for the non-linear weights to maintain the theoretical convergence properties of the optimum reconstruction, whatever the problem considered.The resulting scheme is a fifth-order WENO method based on central interpolation and TVD Runge–Kutta time-integration. We call this scheme the CWENO5 scheme.Numerical experiments in the scalar and the 1D Euler cases make it possible to check and to validate the options selected. In those experiments, we emphasize the resolution power of the method by computing test cases that model realistic aero-acoustic problems. Finally, the new algorithm is directly extended to bi-dimensional problems.