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

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

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

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

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

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

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

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

基于变分原理的二维热传导方程差分格式

研究二维热传导方程的差分数值模拟.用变分原理在不规则结构网格上建立热流通量形式的差分格式.将热流通量作为未知函数求泛函极值,并与温度函数联立求解.克服通常九点格式用插值方法计算网格边界上的热传导系数和网格结点上的温度所引入的误差.     

多块结构网格上的Kershaw扩散格式

Kershaw格式是在四边形结构网格上求解扩散方程的-种经典格式.基于对Kershaw格式中"流"的深入理解,将其拓展到包含非结构点的多块结构网格,分别推导退化非结构点和强化非结构点情况下的Kershaw格式,拓展的Kershaw格式满足流连续条件.三个数值算例的计算结果与精确解吻合得很好,表明将Kershaw格式拓展到多块结构网格的正确性和有效性.     

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

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

多流管网格上九点格式的节点插值方法（英文）

多流管方法是二维多介质辐射流体力学数值模拟中一类常用的求解方法,它采用Lagrange-Euler混合型四边形网格,称为多流管网格。通常其网格品质高于一般的四边形网格。在这类网格上,可以利用网格特性对九点扩散格式中的节点插值方法进行改进。本文利用调和平均点和梯度离散构造的方法提出几种节点插值方法。并给出数值实验,说明现有应用程序中的节点插值方法损失精度,而新的节点插值方法能够使得九点格式在多流管网格上具有二阶精度。     

多流管网格上九点格式的节点插值方法

多流管方法是二维多介质辐射流体力学数值模拟中一类常用的求解方法,它采用Lagrange-Euler混合型四边形网格,称为多流管网格.通常其网格品质高于一般的四边形网格.在这类网格上,可以利用网格特性对九点扩散格式中的节点插值方法进行改进.本文利用调和平均点和梯度离散构造的方法提出几种节点插值方法.并给出数值实验,说明现...     

四阶紧致差分格式对静力模式中垂直网格计算频散性的影响

为了说明四阶紧致差分格式在大气和海洋数值模式中的潜在价值,提出一种通用方法,推导静力线性斜压适应方程组在微分和差分情况下的频散关系,水平尺度分100 km,10 km和1 km三种情况,从频率、水平群速和垂直群速方面,对采用二阶中央差和四阶紧致差分格式情况下,非跳点网格(N网格)、Lorenz网格(L网格)、Charney-Phillips网格(CP网格)、Lorenz时间跳点网格(LTS网格)和Charney-Phillips时间跳点网格(CPTS网格)的计算特性进行比较,发现采用高精度的四阶紧致差分格式总体上可以明显减少上述三种水平尺度波动在N网格、CP网格、L网格和CPTS网格上的频率、水平群速和垂直群速误差,但对LTS网格,采用四阶紧致差分格式,会使得计算水平群速和垂直群速误差变大.     

The finite volume scheme preserving extremum principle for diffusion equations on polygonal meshes

We construct a new nonlinear finite volume scheme for diffusion equation on polygonal meshes and prove that the scheme satisfies the discrete extremum principle. Our scheme is locally conservative and has only cell-centered unknowns. Numerical results are presented to show how our scheme works for preserving discrete extremum principle and positivity on various distorted meshes.     

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.     

非结构四边形网格下的一类保对称有限体元格式

针对定常扩散问题,在非结构四边形网格下,通过选取特殊的控制体和有限体元空间,建立两种保对称有限体元格式,在拟一致网格剖分下,当扩散系数光滑时,证明有限体元解函数在L²和H¹范数下均具有饱和误差阶.数值实验验证理论结果的正确性,同时表明新格式对扭曲大变形四边形网格、间断系数问题具有较强的适应性.在正交网格下,第二种格式对流(flux)函数在单元中心点的值还具有超逼近性.     

An efficient and accurate reconstruction algorithm for the formulation of cell-centered diffusion schemes

A new reconstruction algorithm is proposed for constructing cell-centered diffusion schemes on distorted meshes. Its main feature is that edge unknowns are defined at certain balance points, the locations of which depend on the diffusion coefficient and the skewness of grid cells, so as to obtain a two-point reconstruction stencil. Implementing the new algorithm for the approximation of gradients, we extend the IDC (improved deferred correction) scheme, which was proposed by Traoré et al. [P. Traoré, Y. Ahipo, C. Louste, A robust and efficient finite volume scheme for the discretization of diffusive flux on extremely skewed meshes in complex geometries, J. Comput. Phys. 228 (2009) 5148–5159], to handle diffusion problems with discontinuous coefficients. Numerical results demonstrate the accuracy and efficiency of the extended scheme.     

各向异性扩散问题Kershaw格式的守恒保正修复算法

针对各向异性扩散方程Kershaw格式的数值解在正交网格及扭曲网格上计算出负的现象,给出一种守恒的保正修复算法（CENZ）,该算法对简单遇负置零（ENZ）方法进行改进,使修复后的数值解不仅具有非负性,而且保持法向通量的局部守恒性.数值算例表明,该方法不受计算网格类型和扩散系数各向异性比的限制,可用于对任意违背单调性（或保正性）的有限体积格式数值解的修复.     

A Positivity-Preserving Finite Volume Scheme for Nonequilibrium Radiation Diffusion Equations on Distorted Meshes

In this paper, we propose a new positivity-preserving finite volume scheme with fixed stencils for the nonequilibrium radiation diffusion equations on distorted meshes. This scheme is used to simulate the equations on meshes with both the cell-centered and cell-vertex unknowns. The cell-centered unknowns are the primary unknowns, and the element vertex unknowns are taken as the auxiliary unknowns, which can be calculated by interpolation algorithm. With the nonlinear two-point flux approximation, the interpolation algorithm is not required to be positivity-preserving. Besides, the scheme has a fixed stencil and is locally conservative. The Anderson acceleration is used for the Picard method to solve the nonlinear systems efficiently. Several numerical results are also given to illustrate the efficiency and strong positivity-preserving quality of the scheme.     

Linearity preserving nine-point schemes for diffusion equation on distorted quadrilateral meshes

In this paper, we employ the so-called linearity preserving method, which requires that a difference scheme should be exact on linear solutions, to derive a nine-point difference scheme for the numerical solution of diffusion equation on the structured quadrilateral meshes. This scheme uses firstly both cell-centered unknowns and vertex unknowns, and then the vertex unknowns are treated as a linear combination of the surrounding cell-centered unknowns, which reduces the scheme to a cell-centered one. The weights in the linear combination are derived through the linearity preserving approach and can be obtained by solving a local linear system whose solvability is rigorously discussed. Moreover, the relations between our linearity preserving scheme and some existing schemes are also discussed, by which a generalized multipoint flux approximation scheme based on the linearity preserving criterion is suggested. Numerical experiments show that the linearity preserving schemes in this paper have nearly second order accuracy on many highly skewed and highly distorted structured quadrilateral meshes.     

Compact finite volume schemes on boundary-fitted grids

The paper focuses on the development of a framework for high-order compact finite volume discretization of the three-dimensional scalar advection–diffusion equation. In order to deal with irregular domains, a coordinate transformation is applied between a curvilinear, non-orthogonal grid in the physical space and the computational space. Advective fluxes are computed by the fifth-order upwind scheme introduced by Pirozzoli [S. Pirozzoli, Conservative hybrid compact-WENO schemes for shock–turbulence interaction, J. Comp. Phys. 178 (2002) 81] while the Coupled Derivative scheme [M.H. Kobayashi, On a class of Padé finite volume methods, J. Comp. Phys. 156 (1999) 137] is used for the discretization of the diffusive fluxes.Numerical tests include unsteady diffusion over a distorted grid, linear free-surface gravity waves in a irregular domain and the advection of a scalar field. The proposed methodology attains high-order formal accuracy and shows very favorable resolution characteristics for the simulation of problems with a wide range of length scales.