各向异性扩散方程的高精度算法

针对多介质各向异性扩散方程,本文设计了一种非结构多边形网格高精度有限体积计算格式.为了能适应网格大变形,在构造格式框架时除了用到单元中心量外还引入了节点量作为中间变量,并通过推广孪生逼近算法于各向异性扩散系数情形消除节点量,使算法回归于单元中心量计算流程.数值算例表明,该方法能较好适应大变形网格及间断系数各向异性扩散方程计算.     

《计算数学》2021年第43卷第1期摘要

袁光伟.非正交网格上满足极值原理的扩散格式[J].计算数学,2021,43(1)1-16.摘要:构造了非正交网格上扩散方程新的非线性单元中心型有限体积格式,证明了该格式满足离散极值原理.且在适当条件下具有强制性、以及在离散H1范数下解的有界性和一阶收敛性.     

非结构任意多边形网格辐射扩散方程有限体积格式

本文基于非结构任意多边形网格体系,给出了求解辐射扩散方程的中心型有限体积格式,格式中出现的网格节点未知量由相邻的网格中心未知量加权给出,综合考虑网格几何及扩散系数的影响,给出了节点未知量的一种加权方式,数值实验表明格式在各种非结构网格上具有较强的适应性.     

基于非结构网格的Gas-Kinetic方法

为解决带有复杂几何边界条件的高速流体计算问题,提出基于非结构网格的Gas-Kinetic方法.对于二维非结构网格,以三角形网格作为计算单元,形成在该网格控制单元中物理量导数求解的新方法.通过物理量导数得到在控制体积元边界上的通量,然后用每个计算时间步中求出的边界通量和控制体积元中的物理量,求出下一计算时间步所需的新物理量,依次进行计算直到计算结果收敛为止.采用NACA0012翼型进行数值计算验证,结果表明该方法简单高效,适用于低速和高速流体的计算.     

可互溶液体的真实感模拟

液体互溶混合现象在日常生活中非常常见,但是由于不同液体之间的交互扩散过程非常复杂,对于这类现象的真实感模拟非常困难.为此,提出一种基于体积函数的方法来模拟液体之间的互溶扩散现象.该方法中,流体方程采用二阶精度的有限体积法求解,通过在每个计算网格单元中跟踪记录每种液体组分在网格单元内所占的体积比例,来模拟不同液体之间的交互扩散过程;并利用基于八叉树的自适应网格细分算法对计算过程进行加速.实验结果表明,与已有方法相比,文中方法更加稳定,更容易获得细节丰富的液体交互模拟结果.     

基于物理的波浪实时模拟方法

为解决广域范围内波浪状态的实时计算和可视化问题,结合无粘滞流体力学的物理模型,提出了一种以三角形为控制单元的有限体积简化算法。该算法的优点在于：以不规则边界区域的非规则三角网格为基础,通过简化通量向量分裂方法获得三角形控制单元的边界数值通量,能快速逼近二维浅水方程的解进而模拟非规则边界浅水的实时流动。实验结果显示,所提方法能在符合现实世界物理规律的前提下较好地实现大规模波浪的实时可视化模拟。     

三维非结构网格上求解双曲型守恒律方程的一类三阶精度有限体积格式

考虑标量双曲型守恒律方程,对三维非结构四面体网格给出了一类满足局部极值原理的三阶精度有限体积格式.方法的主要思想是时间和空间分开处理,时间离散采用三阶TVD Runge-Kutta方法;对空间,在每一个四面体单元上基于最小二乘原理构造一个二次多项式,结合数值解光滑探测器和梯度限制器,使其在光滑区域具有高阶精度,在间断区域满足局部极值原理.该格式具有间断分辨能力高,编程实现简便,计算速度快等优点.典型算例的数值试验表明,该格式是有效的.     

一类Lagrange坐标系下的ENO有限体积格式

本文首先从积分形式的二维Lagrange流体力学方程组出发,使用ENO高阶插值多项式,推广了四边形结构网格下的一阶有限体积格式,构造得到了一类结构网格下的高精度有限体积格式.该格式针对单介质问题具有良好的计算效果,同时在处理多介质问题时,不会产生物质界面附近强烈的震荡.结合有效的守恒重映方法,用ALE方法进行数值模拟,得到了预期的效果.     

二维非结构网格的非振荡有限体积方法

1.引言 自从1983年Harten提出了TVD格式后,高分辨率有限差分方法(TVD,ENO等)在计算流体力学领域已经得到了广泛的应用,并取得了很好的计算效果,但对几何形状非常复杂的计算问题,有限差分方法有一定的局限性.非结构网格有限体积方法可以计算任何几何形状的二,三维问题,所以对非结构网格以及有限体积方法的研究越来越受到人们的重     

二维非结构网格浅水波方程的二阶精度非振荡有限体积方法

考虑浅水波方程,对二维非结构网格给出了一种非振荡有限体积方法.该方法的主要思想是在每一个三角形单元上采用最小二乘的思想构造一个重构函数,而时间离散采用二步TVD Runge- Kutta方法.最后用该格式对二维溃坝问题进行了数值试验,得到了满意的结果.     

Uniform Second Order Convergence of a Complete Flux Scheme on Unstructured 1D Grids for a Singularly Perturbed Advection–Diffusion Equation and Some Multidimensional Extensions

The accurate and efficient discretization of singularly perturbed advection–diffusion equations on arbitrary 2D and 3D domains remains an open problem. An interesting approach to tackle this problem is the complete flux scheme (CFS) proposed by G. D. Thiart and further investigated by J. ten Thije Boonkkamp. For the CFS, uniform second order convergence has been proven on structured grids. We extend a version of the CFS to unstructured grids for a steady singularly perturbed advection–diffusion equation. By construction, the novel finite volume scheme is nodally exact in 1D for piecewise constant source terms. This property allows to use elegant continuous arguments in order to prove uniform second order convergence on unstructured one-dimensional grids. Numerical results verify the predicted bounds and suggest that by aligning the finite volume grid along the velocity field uniform second order convergence can be obtained in higher space dimensions as well.     

非正规四边形网格二维抛物方程的有限体积差分方法

扩散方程的数值模拟是计算流体力学和数值热传导问题中的一个重要的基础性课题。 热传导数值计算中,需要计算各种非线性的扩散方程,扩散方程的数值模拟是各种线性、非线性的流体力学方程数值计算的基础,研究扩散方程的高精度,高效率和守恒的数值方法,     

Non-unified Compact Residual-Distribution Methods for Scalar Advection–Diffusion Problems

This paper solves the advection–diffusion equation by treating both advection and diffusion residuals in a separate (non-unified) manner. An alternative residual distribution (RD) method combined with the Galerkin method is proposed to solve the advection–diffusion problem. This Flux-Difference RD method maintains a compact-stencil and the whole process of solving advection–diffusion does not require additional equations to be solved. A general mathematical analysis reveals that the new RD method is linearity preserving on arbitrary grids for the steady-state advection–diffusion equation. The numerical results show that the flux difference RD method preserves second-order accuracy on various unstructured grids including highly randomized anisotropic grids on both the linear and nonlinear scalar advection–diffusion cases.     

A unifying grid approach for solving potential flows applicable to structured and unstructured grid configurations

In this study, an efficient numerical method is proposed for unifying the structured and unstructured grid approaches for solving the potential flows. The new method, named as the “alternating cell directions implicit - ACDI”, solves for the structured and unstructured grid configurations equally well. The new method in effect applies a line implicit method similar to the Line Gauss Seidel scheme for complex unstructured grids including mixed type quadrilateral and triangle cells. To this end, designated alternating directions are taken along chains of contiguous cells, i.e. ‘cell directions’, and an ADI-like sweeping is made to update these cells using a Line Gauss Seidel like scheme. The algorithm makes sure that the entire flow field is updated by traversing each cell twice at each time step for unstructured quadrilateral grids that may contain triangular cells. In this study, a cell-centered finite volume formulation of the ACDI method is demonstrated. The solutions are obtained for incompressible potential flows around a circular cylinder and a forward step. The results are compared with the analytical solutions and numerical solutions using the implicit ADI and the explicit Runge-Kutta methods on single-and multi-block structured and unstructured grids. The results demonstrate that the present ACDI method is unconditionally stable, easy to use and has the same computational performance in terms of convergence, accuracy and run times for both the structured and unstructured grids.     

Application of Second Order TVD and WENO Schemes in Internal Aerodynamics

We deal with the comparison of several finite volume TVD schemes and finite difference ENO schemes and we describe a second order finite volume WENO scheme which was developed for the case of general unstructured meshes. The proposed second order WENO reconstruction is much simpler than the original ENO scheme introduced in [Harten and Chakravarthy 1991]. Moreover, the proposed WENO method is very easily extendible for unstructured meshes in 3D. All above mentioned schemes are applied for the solution of 2D and 3D transonic flows in the turbines and channels and the numerical solution is compared to experimental results or to the results obtained by other authors.     

The auxiliary space method and optimal multigrid preconditioning techniques for unstructured grids

An abstract framework ofauxiliary space method is proposed and, as an application, an optimal multigrid technique is developed for general unstructured grids. The auxiliary space method is a (nonnested) two level preconditioning technique based on a simple relaxation scheme (smoother) and an auxiliary space (that may be roughly understood as a nonnested coarser space). An optimal multigrid preconditioner is then obtained for a discretized partial differential operator defined on an unstructured grid by using an auxiliary space defined on a more structured grid in which a furthernested multigrid method can be naturally applied. This new technique makes it possible to apply multigrid methods to general unstructured grids without too much more programming effort than traditional solution methods. Some simple examples are also given to illustrate the abstract theory and for instance the Morley finite element space is used as an auxiliary space to construct a preconditioner for Argyris element for biharmonic equations. Some numerical results are also given to demonstrate the efficiency of using structured grid for auxiliary space to precondition unstructured grids.     

Spectral Difference Method for Unstructured Grids II: Extension to the Euler Equations

An efficient, high-order, conservative method named the spectral difference method has been developed recently for conservation laws on unstructured grids. It combines the best features of structured and unstructured grid methods to achieve high-computational efficiency and geometric flexibility; it utilizes the concept of discontinuous and high-order local representations to achieve conservation and high accuracy; and it is based on the finite-difference formulation for simplicity. The method is easy to implement since it does not involve surface or volume integrals. Universal reconstructions are obtained by distributing solution and flux points in a geometrically similar manner for simplex cells. In this paper, the method is further extended to nonlinear systems of conservation laws, the Euler equations. Accuracy studies are performed to numerically verify the order of accuracy. In order to capture both smooth feature and discontinuities, monotonicity limiters are implemented, and tested for several problems in one and two dimensions. The method is more efficient than the discontinuous Galerkin and spectral volume methods for unstructured grids.     

Parallel adaptive solution for two dimensional 3-T energy equation on UG

Two-dimensional (2D) energy equation coupled with three temperatures such as electron, ion and photon is widely used to approximately describe the evolution of radiation energy across multiple materials and to study the exchange of energy among electrons, ions and photons for numerical research on laser-driven implosion of a fuel capsule in inertial confinement fusion experiments. The numerical solution of such equations is always fascinating because of its strongly nonlinear phenomena and strongly discontinuous interfaces. Using the UG framework, this paper successfully solves such equations on 2D unstructured grids with a fully implicit finite volume discretization scheme and parallel adaptive multigrid. Significant numerical results using 32 processors are given and analyzed.     

DIMEX Runge–Kutta finite volume methods for multidimensional hyperbolic systems

We propose a class of finite volume methods for the discretization of time-dependent multidimensional hyperbolic systems in divergence form on unstructured grids. We discretize the divergence of the flux function by a cell-centered finite volume method whose spatial accuracy is provided by including into the scheme non-oscillatory piecewise polynomial reconstructions. We assume that the numerical flux function can be decomposed in a convective term and a non-convective term. The convective term, which may be source of numerical stiffness in high-speed flow regions, is treated implicitly, while the non-convective term is always discretized explicitly. To this purpose, we use the diagonally implicit–explicit Runge–Kutta (DIMEX-RK) time-marching formulation. We analyze the structural properties of the matrix operators that result from coupling finite volumes and DIMEX-RK time-stepping schemes by using M-matrix theory. Finally, we show the behavior of these methods by some numerical examples.     

Toward a higher order unsteady finite volume solver based on reproducing kernel methods

During the last decades, research efforts are headed to develop high order methods on CFD and CAA to reach most industrial applications (complex geometries) which need, in most cases, unstructured grids. Today, higher-order methods dealing with unstructured grids remain in infancy state and they are still far from the maturity of structured grids-based methods when solving unsteady cases. From this point of view, the development of higher order methods for unstructured grids become indispensable. The finite volume method seems to be a good candidate, but unfortunately it is difficult to achieve space flux derivation schemes with very high order of accuracy for unsteady cases. In this paper we propose, a high order finite volume method based on Moving Least Squares approximations for unstructured grids that is able to reach an arbitrary order of accuracy on unsteady cases. In order to ensure high orders of accuracy, two strategies were explored independently: (1) a zero-mean variables reconstruction to enforce the mean order at the time derivative and (2) a pseudo mass matrix formulation to preserve the residuals order.