基于有限体积水锤方程的Godunov格式离散

基于有限体积法建立了管道瞬变流的离散格式,采用特征分解技术计算界面通量,并通过重构和通量限制建立二阶精度的TVD格式。此格式保证了质量和动量的守恒性,物理意义明确,计算速度快,适用范围广。通过Riemann问题算例、水锤实例和各种阀门组合情况下的管网水锤实例验证了格式具有高分辨率、无虚假振荡和对克朗数灵敏度低等优点。     

An exact multiphase Riemann solver for compressible cavitating flows

This paper presents a new one-fluid method for simulating formation and the collapse of cavitation regions in water during an isothermal process. In this method, the fluid phase changes are included into the wave pattern of an exact Riemann solver. The model behavior is assessed by comparing the numerical results with the other numerical models for several 1D Riemann problems. One-dimensional water hammer problems with vast creation and collapsing of cavitation zones are simulated as well—and the numerical results are compared to experimental results. The new model results are in very good agreement with accepted results reported in the literature. The presented results clearly show that the new model is able to capture the various behaviors of water during the phase change in the saturation dome and the vapor state, which was neglected in previous studies. Finally, the new model is adopted to an ALE method on an adaptive triangular grid to simulate an underwater explosion phenomenon inside a rigid cylinder—and the results are compared with other simulations.     

Explicit solution to the exact Riemann problem and application in nonlinear shallow‐water equations

The Riemann solver is the fundamental building block in the Godunov‐type formulation of many nonlinear fluid‐flow problems involving discontinuities. While existing solvers are obtained either iteratively or through approximations of the Riemann problem, this paper reports an explicit analytical solution to the exact Riemann problem. The present approach uses the homotopy analysis method to solve the nonlinear algebraic equations resulting from the Riemann problem. A deformation equation defines a continuous variation from an initial approximation to the exact solution through an embedding parameter. A Taylor series expansion of the exact solution about the embedding parameter provides a series solution in recursive form with the initial approximation as the zeroth‐order term. For the nonlinear shallow‐water equations, a sensitivity analysis shows fast convergence of the series solution and the first three terms provide highly accurate results. The proposed Riemann solver is implemented in an existing finite‐volume model with a Godunov‐type scheme. The model correctly describes the formation of shocks and rarefaction fans for both one and two‐dimensional dam‐break problems, thereby verifying the proposed Riemann solver for general implementation. Copyright © 2007 John Wiley & Sons, Ltd.     

固定网格下的特征线法求解溃坝问题

溃坝问题是典型的非线性双曲方程的Riemann问题,其数值求解的难点在于对间断面的捕捉以及避免间断面处在数值计算过程中产生数值色散,因而为求解此问题所产生的各种数值计算方法的优劣也体现在这两个方面。本文针对溃坝问题提出一种新的计算方法。该方法基于对偶变量推导的浅水波方程,根据方程的特点,从方程的特征值和黎曼不变量出发,采用高精度的激波捕捉方法计算黎曼不变量的位置随时间的变化,然后映射至不随时间变化的固定网格。根据黎曼不变量的位置,采用保形分段三次Hermite插值将物理量映射至网格节点。计算结果显示,该方法不仅操作简单,计算量小,而且结果准确。     

A general approximate‐state Riemann solver for hyperbolic systems of conservation laws with source terms

An approximate‐state Riemann solver for the solution of hyperbolic systems of conservation laws with source terms is proposed. The formulation is developed under the assumption that the solution is made of rarefaction waves. The solution is determined using the Riemann invariants expressed as functions of the components of the flux vector. This allows the flux vector to be computed directly at the interfaces between the computational cells. The contribution of the source term is taken into account in the governing equations for the Riemann invariants. An application to the water hammer equations and the shallow water equations shows that an appropriate expression of the pressure force at the interface allows the balance with the source terms to be preserved, thus ensuring consistency with the equations to be solved as well as a correct computation of steady‐state flow configurations. Owing to the particular structure of the variable and flux vectors, the expressions of the fluxes are shown to coincide partly with those given by the HLL/HLLC solver. Computational examples show that the approximate‐state solver yields more accurate solutions than the HLL solver in the presence of discontinuous solutions and arbitrary geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

Multigrid third‐order least‐squares solution of Cauchy–Riemann equations on unstructured triangular grids

In this paper, a multigrid algorithm is developed for the third‐order accurate solution of Cauchy–Riemann equations discretized in the cell‐vertex finite‐volume fashion: the solution values stored at vertices and the residuals defined on triangular elements. On triangular grids, this results in a highly overdetermined problem, and therefore we consider its solution that minimizes the residuals in the least‐squares norm. The standard second‐order least‐squares scheme is extended to third‐order by adding a high‐order correction term in the residual. The resulting high‐order method is shown to give sufficiently accurate solutions on relatively coarse grids. Combined with a multigrid technique, the method then becomes a highly accurate and efficient solver. We present some results to demonstrate its accuracy and efficiency, including both structured and unstructured triangular grids. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical modeling of multiphase first-contact miscible flows. Part 1. Analytical Riemann solver

In this series of two papers, we present a front-tracking method for the numerical simulation of first-contact miscible gas injection processes. The method is developed for constructing very accurate (or even exact) solutions to one-dimensional initial-boundary-value problems in the form of a set of evolving discontinuities. The evolution of the discontinuities is given by analytical solutions to Riemann problems. In this paper, we present the mathematical model of the problem and the complete Riemann solver, that is, the analytical solution to the one-dimensional problem with piecewise constant initial data separated by a single discontinuity, for any left and right states. The Riemann solver presented here is the building block for the front-tracking/streamline method described and applied in the second paper.     

黏弹性管道水柱分离弥合水锤有限体积法模型

基于有限体积法二阶Godunov求解格式对黏弹性输水管道中水柱分离弥合现象进行建模和模拟研究. 在传统的弹性管道模型基础上考虑管道黏弹性效应的影响. 在瞬变流控制方程中引入管道黏弹性项和动态摩阻项, 采用有限体积法进行求解, 考虑压力修正系数来模拟自由气体对计算单元的影响, 同时为避免数值模拟结果产生虚假震荡引入斜率限制器MINMOD函数; 通过虚拟单元法进行边界构建, 实现了计算区域的统一计算. 将所建模型计算结果与已有模型结果、试验结果进行对比, 并对影响模型的各参数进行敏感性分析. 结果表明, 本文模型能够准确模拟出纯水锤、水柱分离弥合水锤两种情况下的瞬态压力变化, 均能与试验数据高度吻合; 与传统的特征线方法相比, 当库朗数C_r小于1时, 有限体积法二阶Godunov格式计算结果更准确、稳定; 在压力波动的衰减过程中, 黏弹性效应相比于管道摩阻起主导作用; 与弹性管道模型相比, 考虑管道黏弹性效应后可显著提高模拟结果的准确性, 尤其是压力波峰值的相对误差明显降低.      

A generalized Rusanov method for the Baer‐Nunziato equations with application to DDT processes in condensed porous explosives

The paper addresses a numerical approach for solving the Baer‐Nunziato equations describing compressible 2‐phase flows. We are developing a finite‐volume method where the numerical flux is approximated with the Godunov scheme based on the Riemann problem solution. The analytical solution to this problem is discussed, and approximate solvers are considered. The obtained theoretical results are applied to develop the discrete model that can be treated as an extension of the Rusanov numerical scheme to the Baer‐Nunziato equations. Numerical results are presented that concern the method verification and also application to the deflagration‐to‐detonation transition (DDT) in porous reactive materials.     

Finite‐volume multi‐stage schemes for shallow‐water flow simulations

A finite‐volume multi‐stage (FMUSTA) scheme is proposed for simulating the free‐surface shallow‐water flows with the hydraulic shocks. On the basis of the multi‐stage (MUSTA) method, the original Riemann problem is transformed to an independent MUSTA mesh. The local Lax–Friedrichs scheme is then adopted for solving the solution of the Riemann problem at the cell interface on the MUSTA mesh. The resulting first‐order monotonic FMUSTA scheme, which does not require the use of the eigenstructure and the special treatment of entropy fixes, has the generality as well as simplicity. In order to achieve the high‐resolution property, the monotonic upstream schemes for conservation laws (MUSCL) method are used. For modeling shallow‐water flows with source terms, the surface gradient method (SGM) is adopted. The proposed schemes are verified using the simulations of six shallow‐water problems, including the 1D idealized dam breaking, the steady transcritical flow over a hump, the 2D oblique hydraulic jump, the circular dam breaking and two dam‐break experiments. The simulated results by the proposed schemes are in satisfactory agreement with the exact solutions and experimental data. It is demonstrated that the proposed FMUSTA schemes have superior overall numerical accuracy among the schemes tested such as the commonly adopted Roe and HLL schemes. Copyright © 2007 John Wiley & Sons, Ltd.     

An implicit numerical algorithm for solving non‐hydrostatic free‐surface flow problems

Details are given of the development of a two‐dimensional vertical numerical model for simulating unsteady free‐surface flows, using a non‐hydrostatic pressure distribution. In this model, the Reynolds equations and the kinematic free‐surface boundary condition are solved simultaneously, so that the water surface elevation can be integrated into the solution and solved for, together with the velocity and pressure fields. An efficient numerical algorithm has been developed, deploying implicit parameters similar to those used in the Crank–Nicholson method, and generating a block tri‐diagonal algebraic system of equations. The model has been applied to simulate a range of unsteady flow problems involving relatively strong vertical accelerations. The results show that the numerical algorithm described is able to produce accurate predictions and is also easy to apply. Copyright © 2001 John Wiley & Sons, Ltd.     

An upstream flux‐splitting finite‐volume scheme for 2D shallow water equations

An upstream flux‐splitting finite‐volume (UFF) scheme is proposed for the solutions of the 2D shallow water equations. In the framework of the finite‐volume method, the artificially upstream flux vector splitting method is employed to establish the numerical flux function for the local Riemann problem. Based on this algorithm, an UFF scheme without Jacobian matrix operation is developed. The proposed scheme satisfying entropy condition is extended to be second‐order‐accurate using the MUSCL approach. The proposed UFF scheme and its second‐order extension are verified through the simulations of four shallow water problems, including the 1D idealized dam breaking, the oblique hydraulic jump, the circular dam breaking, and the dam‐break experiment with 45° bend channel. Meanwhile, the numerical performance of the UFF scheme is compared with those of three well‐known upwind schemes, namely the Osher, Roe, and HLL schemes. It is demonstrated that the proposed scheme performs remarkably well for shallow water flows. The simulated results also show that the UFF scheme has superior overall numerical performances among the schemes tested. Copyright © 2005 John Wiley & Sons, Ltd.     

A parallel unstructured dynamic mesh adaptation algorithm for 3‐D unsteady flows

An unstructured dynamic mesh adaptation and load balancing algorithm has been developed for the efficient simulation of three‐dimensional unsteady inviscid flows on parallel machines. The numerical scheme was based on a cell‐centred finite‐volume method and the Roe's flux‐difference splitting. Second‐order accuracy was achieved in time by using an implicit Jacobi/Gauss–Seidel iteration. The resolution of time‐dependent solutions was enhanced by adopting an h‐refinement/coarsening algorithm. Parallelization and load balancing were concurrently achieved on the adaptive dynamic meshes for computational speed‐up and efficient memory redistribution. A new tree data structure for boundary faces was developed for the continuous transfer of the communication data across the parallel subdomain boundary. The parallel efficiency was validated by applying the present method to an unsteady shock‐tube problem. The flows around oscillating NACA0012 wing and F‐5 wing were also calculated for the numerical verification of the present dynamic mesh adaptation and load balancing algorithm. Copyright © 2005 John Wiley & Sons, Ltd.     

Modelling the Damage Induced by Pressure Transients in Elasto-Plastic Pipes

This paper is concerned with the analysis of pressure transients in damageable elasto-plastic piping systems. The fluid dynamics and pipewall deformation are modelled by the classical water hammer theory, whereas pipewall mechanical behavior is described by an internal variable constitutive theory. The constitutive model coupling plasticity and damage used herein gives rise to a nonlinear hyperbolic problem in which the wavespeeds are altered by damage evolution. The problem is numerically approximated by means of a technique based on an additive decomposition together with the Glimm's method and a special Euler-type time integration scheme. Examples concerning the structural integrity analysis of a reservoir-pipe-valve installation, where hydraulic transients are generated by valve slam, are presented to illustrate the applicability of both theory and numerical method.     

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.     

The reactive Riemann problem for thermally perfect gases at all combustion regimes

In this work we analyze the reactive Riemann problem for thermally perfect gases in the deflagration or detonation regimes. We restrict our attention to the case of one irreversible infinitely fast chemical reaction; we also suppose that, in the initial condition, one state (for instance the left one) is burnt and the other one is unburnt. The indeterminacy of the deflagration regime is removed by imposing a (constant) value for the fundamental flame speed of the reactive shock. An iterative algorithm is proposed for the solution of the reactive Riemann problem. Then the reactive Riemann problem and the proposed algorithm are investigated from a numerical point of view in the case in which the unburnt state consists of a stoichiometric mixture of hydrogen and air at almost atmospheric condition. In particular, we revisit the problem of 1D plane‐symmetric steady flames in a semi‐infinite domain and we verify that the transition from one combustion regime to another occurs continuously with respect to the fundamental flame speed and the so‐called piston velocity. Finally, we use the ‘all shock’ solution of the reactive Riemann problem to design an approximate (‘all shock’) Riemann solver. 1D and 2D flows at different combustion regimes are computed, which shows that the approximate Riemann solver, and thus the algorithm we use for the solution of the reactive Riemann problem, is robust in both the deflagration and detonation regimes. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical flux functions for Reynolds‐averaged Navier–Stokes and kω turbulence model computations with a line‐preconditioned p‐multigrid discontinuous Galerkin solver

We present an eigen‐decomposition of the quasi‐linear convective flux formulation of the completely coupled Reynolds‐averaged Navier–Stokes and kω turbulence model equations. Based on these results, we formulate different approximate Riemann solvers that can be used as numerical flux functions in a DG discretization. The effect of the different strategies on the solution accuracy is investigated with numerical examples. The actual computations are performed using a p‐multigrid algorithm. To this end, we formulate a framework with a backward‐Euler smoother in which the linear systems are solved with a general preconditioned Krylov method. We present matrix‐free implementations and memory‐lean line‐Jacobi preconditioners and compare the effects of some parameter choices. In particular, p‐multigrid is found to be less efficient than might be expected from recent findings by other authors. This might be due to the consideration of turbulent flow. Copyright © 2012 John Wiley & Sons, Ltd.     

Lattice Boltzmann methods for shallow water flow applications

We apply the lattice Boltzmann (LB) method for solving the shallow water equations with source terms such as the bed slope and bed friction. Our aim is to use a simple and accurate representation of the source terms in order to simulate practical shallow water flows without relying on upwind discretization or Riemann problem solvers. We validate the algorithm on problems where analytical solutions are available. The numerical results are in good agreement with analytical solutions. Furthermore, we test the method on a practical problem by simulating mean flow in the Strait of Gibraltar. The main focus is to examine the performance of the LB method for complex geometries with irregular bathymetry. The results demonstrate its ability to capture the main flow features. Copyright © 2007 John Wiley & Sons, Ltd.     

A new numerical model for simulations of wave transformation,breaking and long‐shore currents in complex coastal regions

In this paper, we propose a model based on a new contravariant integral form of the fully nonlinear Boussinesq equations in order to simulate wave transformation phenomena, wave breaking, and nearshore currents in computational domains representing the complex morphology of real coastal regions. The aforementioned contravariant integral form, in which Christoffel symbols are absent, is characterized by the fact that the continuity equation does not include any dispersive term. A procedure developed in order to correct errors related to the difficulties of numerically satisfying the metric identities in the numerical integration of fully nonlinear Boussinesq equation on generalized boundary‐conforming grids is presented. The Boussinesq equation system is numerically solved by a hybrid finite volume–finite difference scheme. The proposed high‐order upwind weighted essentially non‐oscillatory finite volume scheme involves an exact Riemann solver and is based on a genuinely two‐dimensional reconstruction procedure, which uses a convex combination of biquadratic polynomials. The wave breaking is represented by discontinuities of the weak solution of the integral form of the nonlinear shallow water equations. The capacity of the proposed model to correctly represent wave propagation, wave breaking, and wave‐induced currents is verified against test cases present in the literature. The results obtained are compared with experimental measures, analytical solutions, or alternative numerical solutions. Copyright © 2015 John Wiley & Sons, Ltd.     

二维洪水演进数值模拟

利用非结构化的有限体积方法,建立了二维浅水方程高精度、高分辨率模型。以Roe类型的近似Rie-mann解计算界面通量,通过MUSCL和两步TVD Runge-Kutta法获得了空间和时间都具有二级精度的TVD格式。采用特征分解的方法处理底坡源项和采用半隐式方法处理摩擦源项均能保证了格式的稳定性与和谐性。通过水滴算例对模型进行验证,并应用此模型对98年胖头泡分滞洪区分洪过程进行模拟,获得滞洪区不同时段的淹没范围和淹没水深,为防洪救灾提供了依据。