一种简单的精确捕捉接触间断的黎曼求解器

基于Godunov型数值格式的有限体积法是求解双曲型守恒律系统的主流方法,其中用来计算界面数值通量的黎曼求解器在很大程度上决定了数值格式在计算中的表现。单波的Rusanov求解器和双波的HLL求解器具有简单、高效和鲁棒性好等优点,但是在捕捉接触间断时耗散太大。全波的HLLC格式能够精确捕捉接触间断,但是在计算中出现的激波不稳定现象限制了其在高马赫数流动问题中的应用。本文利用双曲正切函数和五阶WENO格式来重构界面两侧的密度值,并且结合边界变差下降算法来减小Rusanov格式耗散项中的密度差,从而提高格式对于接触间断的分辨率。研究表明,相比于全波的HLLC求解器,本文构造的黎曼求解器不仅具有更高的接触分辨率,而且还具有更好的激波稳定性。     

基于非结构网格求解二维浅水方程的高精度有限体积方法

采用HLL格式,在三角形非结构网格下采用有限体积离散,建立了求解二维浅水方程的高精度的数值模型.本文采用多维重构和多维限制器的方法来获得高精度的空间格式以及防止非物理振荡的产生,时间离散采用三阶Runge-Kutta法以获得高阶的时间精度.基于三角形网格,底坡源项采用简单的斜底模型离散,为保证计算格式的和谐性,对经典的HLL格式计算的数值通量中的静水压力项进行了修正.算例证明本文提出的方法的和谐性并具有高精度的间断捕捉能力和稳定性.     

一种求解Euler方程的新型高阶精度数值方法

提出了一种求解Euler方程的新型高阶精度数值方法.该数值方法基于一种新的矢通量分裂格式,将矢通量项分裂成压力通量项和对流通量项.与传统矢通量分裂格式相比,新的矢通量分裂格式能够更好地捕捉特征场内的中间特征波,从而增强格式的分辨率.同时,为了提高这种矢通量分裂格式的空间精度,我们在近似求解压力通量项黎曼问题时对界面处的独立物理变量进行高阶插值.在时间步上,采用显式最优的三阶龙格-库塔方法进行推进.数值试验表明,与传统数值方法相比,本文提出的新方法同时具有高精度和高分辨率的优点.     

基于混合型WENO格式数值模拟四类Riemann问题

本文基于三阶WENO格式和三阶WENO-Z格式，利用有限体积法研究了同一格式在不同方向使用不同的求解器，以及在不同方向采用相同精度的不同格式数值模拟4类Riemann问题，分析各向异性对数值计算结果的影响．数值结果表明，无论是使用不同的格式，还是使用不同的求解器进行数值模拟，都会导致数值结果不同程度地失去图像对称性．     

一种健壮的混合Roe黎曼求解器

使用Roe格式计算多维流动问题时,在强激波附近会出现数值激波不稳定现象。带有剪切粘性的HLLEC格式不仅可以捕捉接触间断,而且表现出很好的稳定性。混合Roe格式和HLLEC格式来消除数值激波不稳定性。在强激波附近,通过激波面法向和网格界面法向的夹角来定义开关函数,使得数值通量在激波面横向切换成HLLEC格式。在其余地方,数值通量依然使用Roe格式来计算。数值试验表明,混合格式不仅消除了Roe格式的数值激波不稳定性,还最大程度地减少了HLLEC格式所带来的剪切耗散,保留了Roe格式高分辨率的优点。     

基于制造解的非结构二阶有限体积离散格式的精度测试与验证

随着计算机技术的飞速进步,计算流体力学得到迅猛发展,数值计算虽能够快速得到离散结果,但是数值结果的正确性与精度则需要通过严谨的方法来进行验证和确认.制造解方法和网格收敛性研究作为验证与确认的重要手段已经广泛应用于计算流体力学代码验证、精度分析、边界条件验证等方面.本文在实现标量制造解和分量制造解方法的基础上,通过将制造解方法精度测试结果与经典精确解(二维无黏等熵涡)精度测试结果进行对比,进一步证实了制造解精度测试方法的有效性,并将两种制造解方法应用于非结构网格二阶精度有限体积离散格式的精度测试与验证,对各种常用的梯度重构方法、对流通量格式、扩散通量格式进行了网格收敛性精度测试.结果显示,基于Green-Gauss公式的梯度重构方法在不规则网格上会出现精度降阶的情况,导致流动模拟精度严重下降,而基于最小二乘(least squares)的梯度重构方法对网格是否规则并不敏感.对流通量格式的精度测试显示,所测试的各种对流通量格式均能达到二阶精度,且各方法精度几乎相同;而扩散通量离散中界面梯度求解方法的选择对流动模拟精度有显著影响.     

Level set方法在自适应Cartesian网格上的应用

采用基于自适应Cartesian网格的level set方法对多介质流动问题进行数值模拟。采用基于四叉树的方法来生成自适应Cartesian网格。采用有限体积法求解Euler方程,控制面通量的计算采用HLLC（Hartern, Lax, van Leer, Contact）近似黎曼解方法。level set方程也采用有限体积法求解,采用Lax-Friedchs方法计算通量,通过窄带方法来减少计算量,界面的处理采用ghost fluid方法。Runge-Kutta显式时间推进,时间、空间都是二阶精度。对两种不同比热比介质激波管问题进行数值模拟,其结果和精确解吻合;对空气/氦气泡相互作用等问题进行模拟,取得令人满意的结果。     

求解可压缩流的二维通量分裂格式

传统的一维通量分裂格式在计算界面数值通量时,只考虑网格界面法向的波系。采用传统的TV格式分别求解对流通量和压力通量。通过求解考虑了横向波系影响的角点数值通量来构造一种真正二维的TV通量分裂格式。在计算一维数值算例时,该格式与传统的TV格式具有相同的数值通量计算公式,因此其保留了传统的TV格式精确捕捉接触间断和膨胀激波的优点。在计算二维算例时,该格式比传统的TV格式具有更高的分辨率;在计算二维强激波问题时,消除了传统TV格式的非物理现象,表现出更好的鲁棒性;此外,该格式大大提高了稳定性CFL数,从而具有更高的计算效率。因此,本文方法是一种精确、高效并且具有强鲁棒性的数值方法,在可压缩流的数值模拟中具有广阔的应用前景。     

一类格心型ALE有限体积格式方法

现在国内外流行的ALE有限体积格式基本上都基于交错网榕进行格式的离散.该类格武在进行重映时,速度、密度和能量需要分别进行重映计算,效率较低,而且由于速度在网格角点.而密度、能量在网格中心,重映时会出现动能和内能不协调现泉.本文在巳有格心型Lagrange有限体积格式研究的基础上,结合Abgrall R.等关于榕心型格式下的网格角点速度的计算方法,利用最小二乘法进行高阶插值多项式重构,构造了一类新的格心型的高精度Lagrangian有限体积格式,并结合有效的高精度ENO守恒重映方法,获得了一类格心型的高精度ALE有限体积格式.数值试验的结果说明本文的格式是有效的,高精度的,收敛的,并且避免了物理量的不协调现象.     

管系内一维可压缩流动数值建模与特性分析

针对复杂管系内可压缩流体,基于有限体积法,采用HLLC(Harten-Lax-vanLeerContact)格式和黎曼求解器构建了有限控制体数值离散方法,引入虚拟节点用于连接有限控制体,借助虚拟节点给出控制体之间数值通量的计算格式,发展了一种管道内一维流动数值建模方法。针对含有分支管路的管系,在管道连接部位构建了分支管路拟一维流动数值计算模型。基于所发展的一维流动数值方法,建立了变径管道和含60°分支管道内流动计算模型,验证了该方法的收敛性和有效性;基于虚拟节点的数值格式处理变径管激波问题具有一定精度优势。研究了变径管和分支管模型中可压缩流体激波、稀疏波等的传播机理,分析了管径对相邻支管压力的影响,为工程管路设计提供了参考。     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

A Godunov‐type fractional semi‐implicit method based on staggered grid for dam‐break modeling

A semi‐implicit finite volume model based upon staggered grid is presented for solving shallow water equation. The model employs a time‐splitting scheme that uses a predictor–corrector method for the advection term. The fluxes are calculated based on a Riemann solver in the prediction step and a downwind scheme in the correction step. A simple TVD scheme is employed for shock capturing purposes in which the Minmond limiter is used for flux functions. As a consequence of using staggered grid, an ADI method is adopted for solving the discretized equations for 2‐D problems. Several 1‐D and 2‐D flows have been modeled with satisfactory results when compared with analytical and experimental test cases. The model is also capable of simulating supercritical as well as subcritical flow. Copyright © 2010 John Wiley & Sons, Ltd.     

A simple high‐resolution advection scheme

A simple, robust, mass‐conserving numerical scheme for solving the linear advection equation is described. The scheme can estimate peak solution values accurately even in regions where spatial gradients are high. Such situations present a severe challenge to classical numerical algorithms. Attention is restricted to the case of pure advection in one and two dimensions since this is where past numerical problems have arisen. The authors' scheme is of the Godunov type and is second‐order in space and time. The required cell interface fluxes are obtained by MUSCL interpolation and the exact solution of a degenerate Riemann problem. Second‐order accuracy in time is achieved via a Runge–Kutta predictor–corrector sequence. The scheme is explicit and expressed in finite volume form for ease of implementation on a boundary‐conforming grid. Benchmark test problems in one and two dimensions are used to illustrate the high‐spatial accuracy of the method and its applicability to non‐uniform grids. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of the Jacobian evaluation on Newton's solution of the Euler equations

Newton's method is developed for solving the 2‐D Euler equations. The Euler equations are discretized using a finite‐volume method with upwind flux splitting schemes. Both analytical and numerical methods are used for Jacobian calculations. Although the numerical method has the advantage of keeping the Jacobian consistent with the numerical residual vector and avoiding extremely complex analytical differentiations, it may have accuracy problems and need longer execution time. In order to improve the accuracy of numerical Jacobians, detailed error analyses are performed. Results show that the finite‐difference perturbation magnitude and computer precision are the most important parameters that affect the accuracy of numerical Jacobians. A method is developed for calculating an optimal perturbation magnitude that can minimize the error in numerical Jacobians. The accuracy of the numerical Jacobians is improved significantly by using the optimal perturbation magnitude. The effects of the accuracy of numerical Jacobians on the convergence of the flow solver are also investigated. In order to reduce the execution time for numerical Jacobian evaluation, flux vectors with perturbed flow variables are calculated only for neighbouring cells. A sparse matrix solver that is based on LU factorization is used. Effects of different flux splitting methods and higher‐order discretizations on the performance of the solver are analysed. Copyright © 2005 John Wiley & Sons, Ltd.     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

Numerical Solution of the Polymer System by Front Tracking

The paper describes the application of front tracking to the polymer system, an example of a nonstrictly hyperbolic system. Front tracking computes piecewise constant approximations based on approximate Riemann solutions and exact tracking of waves. It is well known that the front tracking method may introduce a blowup of the initial total variation for initial data along the curve where the two eigenvalues of the hyperbolic system are identical. It is demonstrated by numerical examples that the method converges to the correct solution after a finite time, and that this time decreases with the discretization parameter.For multidimensional problems, front tracking is combined with dimensional splitting, and numerical experiments indicate that large splitting steps can be used without loss of accuracy. Typical CFL numbers are in the range 10–20, and comparisons with Riemann free, high-resolution methods confirm the high efficiency of front tracking.The polymer system, coupled with an elliptic pressure equation, models two-phase, three-component polymer flooding in an oil reservoir. Two examples are presented, where this model is solved by a sequential time stepping procedure. Because of the approximate Riemann solver, the method is non-conservative and CFL numbers must be chosen only moderately larger than unity to avoid substantial material balance errors generated in near-well regions after water breakthrough. Moreover, it is demonstrated that dimensional splitting may introduce severe grid orientation effects for unstable displacements that are accentuated for decreasing discretization parameters.     

Numerical scheme for accurately capturing gas migration described by 1D multiphase drift flux model

A numerical scheme is designed and implemented to solve a simplified set of equations modeling 1-D multi-phase flow based on drift flux model in an isothermal setup with phase dissolution. The difficulty in obtaining the analytical Jacobian of the fluxes leads to the difficulty in obtaining an efficient linearized Riemann solver which in turn affects the accuracy in capturing the contact wavefront/gas migration. To address this issue a fully explicit second order finite volume solver based on flux corrected transport (FCT) is implemented. The choice of variables used for limiting the fluxes affects the amount of numerical diffusion and an appropriate choice of the gradient in volume fraction is used. Practical test cases while drilling in the oil and gas industry, of gas injection inside a well annulus and shut-in of the vertical well are presented. The results conclude that the FCT solver is better and efficient for accurately capturing gas migration for multi-phase models with phase behaviour involving slip velocities given by algebraic relations.     

A novel multidimensional solution reconstruction and edge‐based limiting procedure for unstructured cell‐centered finite volumes with application to shallow water dynamics

On unstructured meshes, the cell‐centered finite volume (CCFV) formulation, where the finite control volumes are the mesh elements themselves, is probably the most used formulation for numerically solving the two‐dimensional nonlinear shallow water equations and hyperbolic conservation laws in general. Within this CCFV framework, second‐order spatial accuracy is achieved with a Monotone Upstream‐centered Schemes for Conservation Laws‐type (MUSCL) linear reconstruction technique, where a novel edge‐based multidimensional limiting procedure is derived for the control of the total variation of the reconstructed field. To this end, a relatively simple, but very effective modification to a reconstruction procedure for CCFV schemes, is introduced, which takes into account geometrical characteristics of computational triangular meshes. The proposed strategy is shown not to suffer from loss of accuracy on grids with poor connectivity. We apply this reconstruction in the development of a second‐order well‐balanced Godunov‐type scheme for the simulation of unsteady two‐dimensional flows over arbitrary topography with wetting and drying on triangular meshes. Although the proposed limited reconstruction is independent from the Riemann solver used, the well‐known approximate Riemann solver of Roe is utilized to compute the numerical fluxes, whereas the Green–Gauss divergence formulation for gradient computations is implemented. Two different stencils for the Green–Gauss gradient computations are implemented and critically tested, in conjunction with the proposed limiting strategy, on various grid types, for smooth and nonsmooth flow conditions. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of dam‐ and dyke‐break hydrodynamics on dynamically adaptive quadtree grids

Flooding due to the failure of a dam or dyke has potentially disastrous consequences. This paper presents a Godunov‐type finite volume solver of the shallow water equations based on dynamically adaptive quadtree grids. The Harten, Lax and van Leer approximate Riemann solver with the Contact wave restored (HLLC) scheme is used to evaluate interface fluxes in both wet‐ and dry‐bed applications. The numerical model is validated against results from alternative numerical models for idealized circular and rectangular dam breaks. Close agreement is achieved with experimental measurements from the CADAM dam break test and data from a laboratory dyke break undertaken at Delft University of Technology. Copyright © 2004 John Wiley Sons, Ltd.