连续分层海洋中内波传播的一种数值模式

基于Euler方程,使用有限体积法建立了一种密度为连续分层情况下、适应水深变化的水域中内波传播的数值模式.为了使计算格式能够达到二阶精度,对流项的处理使用了TVD (total variation diminishing)格式.将SIMPLE算法引入连续分层海洋中内波的数值计算,为了简化计算并方便地适应多种TVD格式,在计算预估速度场时采用了显式格式,而没有采用传统的隐式格式;鉴于在原始的SIMPLE算法中没有涉及到由于密度扰动而引起的静水压力场的改变问题,给出了该问题的计算方法.因此改进了SIMPLE算法.出流边界的处理采用阻尼消波和Sommerfeld辐射条件相结合的方式,以使内波得到有效的衰减和释放.将等水深水域的数值解和理论解进行了比较,两者吻合较好;并对存在潜堤时数值计算的不同时刻密度变化的空间分布进行了详细的定性分析.计算结果表明,所建立的数值模式能有效地模拟内波的传播和变形.      

High-order divergence-free velocity reconstruction for free surface flows on unstructured Voronoi meshes

In this paper, we present an efficient semi-implicit scheme for the solution of the Reynolds-averaged Navier-Stokes equations for the simulation of hydrostatic and nonhydrostatic free surface flow problems. A staggered unstructured mesh composed by Voronoi polygons is used to pave the horizontal domain, whereas parallel layers are adopted along the vertical direction. Pressure, velocity, and vertical viscosity terms are taken implicitly, whereas the nonlinear convective terms as well as the horizontal viscous terms are discretized explicitly by using a semi-Lagrangian approach, which requires an interpolation of the three-dimensional velocity field to integrate the flow trajectories backward in time. To this purpose, a high-order reconstruction technique is proposed, which is based on a constrained least squares operator that guarantees a globally and pointwise divergence-free velocity field. A comparison with an analogous reconstruction, which is not divergence-free preserving, is also presented to give evidence of the new strategy. This allows the continuity equation to be satisfied up to machine precision even for high-order spatial discretizations. The reconstructed velocity field is then used for evaluating high-order terms of a Taylor method that is here adopted as ODE integrator for the flow trajectories. The proposed semi-implicit scheme is validated against a set of academic test problems, and proof of convergence up to fourth-order of accuracy in space is shown.     

Semi-implicit finite difference methods for three-dimensional shallow water flow

A semi-implicit finite difference method for the numerical solution of three-dimensional shallow water flows is presented and discussed. The governing equations are the primitive three-dimensional turbulent mean flow equations where the pressure distribution in the vertical has been assumed to be hydrostatic. In the method of solution a minimal degree of implicitness has been adopted in such a fashion that the resulting algorithm is stable and gives a maximal computational efficiency at a minimal computational cost. At each time step the numerical method requires the solution of one large linear system which can be formally decomposed into a set of small three-diagonal systems coupled with one five-diagonal system. All these linear systems are symmetric and positive definite. Thus the existence and uniquencess of the numerical solution are assured. When only one vertical layer is specified, this method reduces as a special case to a semi-implicit scheme for solving the corresponding two-dimensional shallow water equations. The resulting two- and three-dimensional algorithm has been shown to be fast, accurate and mass-conservative and can also be applied to simulate flooding and drying of tidal mud-flats in conjunction with three-dimensional flows. Furthermore, the resulting algorithm is fully vectorizable for an efficient implementation on modern vector computers.     

An efficient operator splitting scheme for three-dimensional hydrodynamic computations

A new efficient numerical method for three-dimensional hydrodynamic computations is presented and discussed in this paper. The method is based on the operator splitting method and combined with Eulerian–Lagrangian method, finite element method and finite difference method. To increase the efficiency and stability of the numerical solutions, the operator splitting method is employed to partition the momentum equations into three parts, according to physical phenomena. A time step is divided into three time substeps. In the first substep, advection and Coriolis force are solved using the explicit Eulerian–Lagrangian method. In the second substep, horizontal diffusion is approximated by implicit FEM in each horizontal layer. In the last substep, the continuity equation is solved by implicit FEM, and vertical diffusion and pressure gradient are discretized by implicit FDM in each nodal column. The stability analysis shows that this method is unconditionally stable. A number of numerical experiments have been performed. The results simulated by the present scheme agree well with analytical solutions and the other documented model results. The method is efficient for 3D shallow water flow computations and fully fits complicated configurations. © 1998 John Wiley & Sons, Ltd.     

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

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

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 divergence-free semi-implicit finite volume scheme for ideal,viscous, and resistive magnetohydrodynamics

In this paper, we present a novel pressure-based semi-implicit finite volume solver for the equations of compressible ideal, viscous, and resistive magnetohydrodynamics (MHD). The new method is conservative for mass, momentum, and total energy, and in multiple space dimensions, it is constructed in such a way as to respect the divergence-free condition of the magnetic field exactly, also in the presence of resistive effects. This is possible via the use of multidimensional Riemann solvers on an appropriately staggered grid for the time evolution of the magnetic field and a double curl formulation of the resistive terms. The new semi-implicit method for the MHD equations proposed here discretizes the nonlinear convective terms as well as the time evolution of the magnetic field explicitly, whereas all terms related to the pressure in the momentum equation and the total energy equation are discretized implicitly, making again use of a properly staggered grid for pressure and velocity. Inserting the discrete momentum equation into the discrete energy equation then yields a mildly nonlinear symmetric and positive definite algebraic system for the pressure as the only unknown, which can be efficiently solved with the (nested) Newton method of Casulli et al. The pressure system becomes linear when the specific internal energy is a linear function of the pressure. The time step of the scheme is restricted by a CFL condition based only on the fluid velocity and the Alfvén wave speed and is not based on the speed of the magnetosonic waves. Being a semi-implicit pressure-based scheme, our new method is therefore particularly well suited for low Mach number flows and for the incompressible limit of the MHD equations, for which it is well known that explicit density-based Godunov-type finite volume solvers become increasingly inefficient and inaccurate because of the more and more stringent CFL condition and the wrong scaling of the numerical viscosity in the incompressible limit. We show a relevant MHD test problem in the low Mach number regime where the new semi-implicit algorithm is a factor of 50 faster than a traditional explicit finite volume method, which is a very significant gain in terms of computational efficiency. However, our numerical results confirm that our new method performs well also for classical MHD test cases with strong shocks. In this sense, our new scheme is a true all Mach number flow solver.     

A composite multigrid method for calculating unsteady incompressible flows in geometrically complex domains

A time-accurate, finite volume method for solving the three-dimensional, incompressible Navier-Stokes equations on a composite grid with arbitrary subgrid overlapping is presented. The governing equations are written in a non-orthogonal curvilinear co-ordinate system and are discretized on a non-staggered grid. A semi-implicit, fractional step method with approximate factorization is employed for time advancement. Multigrid combined with intergrid iteration is used to solve the pressure Poisson equation. Inter-grid communication is facilitated by an iterative boundary velocity scheme which ensures that the governing equations are well-posed on each subdomain. Mass conservation on each subdomain is preserved by using a mass imbalance correction scheme which is secondorder-accurate. Three test cases are used to demonstrate the method's consistency, accuracy and efficiency.     

A semi-implicit lagrangian scheme for 3D shallow water flow with a two-layer turbulence model

A semi-implicit Lagrangian finite difference scheme for 3D shallow water flow has been developed to include an eddy viscosity model for turbulent mixing in the vertical direction. The α-co-ordinate system for the vertical direction has been introduced to give accurate definition of bed and surface boundary conditions. The simple two-layer mixing length model for rough surfaces is used with the standard assumption that the shear stress across the wall region at a given horizontal location is constant. The bed condition is thus defined only by its roughness height (avoiding the need for a friction formula relating to depth-averaged flow, e.g. Chezy, used previously). The method is shown to be efficient and stable with an explicit Lagrangian formulation for convective terms and terms for surface elevation and vertical mixing handled implicitly. The method is applied to current flow around a circular island with gently sloping sides which produce periodic recirculation zones (vortex shedding). Comparisons are made with experimental measurements of velocity using laser Doppler anemometry (time histories at specific points) and surface particle-tracking velocimetry.     

流固耦合分析的一种改进CBS有限元算法

针对流固耦合问题, 发展了一种基于任意拉格朗日-欧拉(ALE)描述有限元法的弱耦合分区算法. 运用半隐式特征线分裂算法求解Navier-Stokes方程, 在压力Poisson 方程中引入质量源项以满足几何守恒律; 运用子块移动技术更新动态网格, 并配以光滑处理防止网格质量下降; 采用Newmark-β 法求解结构运动方程. 为保持流体-结构界面处速度和动量守恒, 利用修正结合界面边界条件方法求解界面处速度通量和动量通量. 运用本方法分别模拟了不同雷诺数下单圆柱横向和两向流致振动、串列双圆柱两向流致振动. 计算表明, 本文方法计算效率高, 计算结果与已有实验和数值计算数据吻合.     

A semi-implicit finite volume scheme for a simplified hydrostatic model for fluid-structure interaction

Simulating fluid-structure interaction problems usually requires a considerable computational effort. In this article, a novel semi-implicit finite volume scheme is developed for the coupled solution of free surface shallow water flow and the movement of one or more floating rigid structures. The model is well-suited for geophysical flows, as it is based on the hydrostatic pressure assumption and the shallow water equations. The coupling is achieved via a nonlinear volume function in the mass conservation equation that depends on the coordinates of the floating structures. Furthermore, the nonlinear volume function allows for the simultaneous existence of wet, dry and pressurized cells in the computational domain. The resulting mildly nonlinear pressure system is solved using a nested Newton method. The accuracy of the volume computation is improved by using a subgrid, and time accuracy is increased via the application of the theta method. Additionally, mass is always conserved to machine precision. At each time step, the volume function is updated in each cell according to the position of the floating objects, whose dynamics is computed by solving a set of ordinary differential equations for their six degrees of freedom. The simulated moving objects may for example represent ships, and the forces considered here are simply gravity and the hydrostatic pressure on the hull. For a set of test cases, the model has been applied and compared with available exact solutions to verify the correctness and accuracy of the proposed algorithm. The model is able to treat fluid-structure interaction in the context of hydrostatic geophysical free surface flows in an efficient and flexible way, and the employed nested Newton method rapidly converges to a solution. The proposed algorithm may be useful for hydraulic engineering, such as for the simulation of ships moving in inland waterways and coastal regions.     

On numerical solutions of the time-dependent Euler equations for incompressible flow

This paper presents finite element methods for the non-stationary Euler equations of a two dimensional inviscid and incompressible flow. For the time discretization, we compare numerical results obtained by the use of a leap-frog scheme and a semi-implicit scheme of order two.     

A well‐balanced upstream flux‐splitting finite‐volume scheme for shallow‐water flow simulations with irregular bed topography

This study extends the upstream flux‐splitting finite‐volume (UFF) scheme to shallow water equations with source terms. Coupling the hydrostatic reconstruction method (HRM) with the UFF scheme achieves a resultant numerical scheme that adequately balances flux gradients and source terms. The proposed scheme is validated in three benchmark problems and applied to flood flows in the natural/irregular river with bridge pier obstructions. The results of the simulations are in satisfactory agreement with the available analytical solutions, experimental data and field measurements. Comparisons of the present results with those obtained by the surface gradient method (SGM) demonstrate the superior stability and higher accuracy of the HRM. The stability test results also show that the HRM requires less CPU time (up to 60%) than the SGM. The proposed well‐balanced UFF scheme is accurate, stable and efficient to solve flow problems involving irregular bed topography. Copyright © 2009 John Wiley & Sons, Ltd.     

A mixed‐interpolation finite element method for incompressible thermal flows of electrically conducting fluids

A new mixed‐interpolation finite element method is presented for the two‐dimensional numerical simulation of incompressible magnetohydrodynamic (MHD) flows which involve convective heat transfer. The proposed method applies the nodal shape functions, which are locally defined in nine‐node elements, for the discretization of the Navier–Stokes and energy equations, and the vector shape functions, which are locally defined in four‐node elements, for the discretization of the electromagnetic field equations. The use of the vector shape functions allows the solenoidal condition on the magnetic field to be automatically satisfied in each four‐node element. In addition, efficient approximation procedures for the calculation of the integrals in the discretized equations are adopted to achieve high‐speed computation. With the use of the proposed numerical scheme, MHD channel flow and MHD natural convection under a constant applied magnetic field are simulated at different Hartmann numbers. The accuracy and robustness of the method are verified through these numerical tests in which both undistorted and distorted meshes are employed for comparison of numerical solutions. Furthermore, it is shown that the calculation speed for the proposed scheme is much higher compared with that for a conventional numerical integration scheme under the condition of almost the same memory consumption. Copyright © 2016 John Wiley & Sons, Ltd.     

A finite element/volume method model of the depth‐averaged horizontally 2D shallow water equations

Analysis of surface water flows is of central importance in understanding and predicting a wide range of water engineering issues. Dynamics of surface water is reasonably well described using the shallow water equations (SWEs) with the hydrostatic pressure assumption. The SWEs are nonlinear hyperbolic partial differential equations that are in general required to be solved numerically. Application of a simple and efficient numerical model is desirable for solving the SWEs in practical problems. This study develops a new numerical model of the depth‐averaged horizontally 2D SWEs referred to as 2D finite element/volume method (2D FEVM) model. The continuity equation is solved with the conforming, standard Galerkin FEM scheme and momentum equations with an upwind, cell‐centered finite volume method scheme, utilizing the water surface elevation and the line discharges as unknowns aligned in a staggered manner. The 2D FEVM model relies on neither Riemann solvers nor high‐resolution algorithms in order to serve as a simple numerical model. Water at a rest state is exactly preserved in the model. A fully explicit temporal integration is achieved in the model using an efficient approximate matrix inversion method. A series of test problems, containing three benchmark problems and three experiments of transcritical flows, are carried out to assess accuracy and versatility of the model. Copyright © 2014 John Wiley & Sons, Ltd.     

二维洪水演进数值模拟

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

Multi-symplectic Runge-Kutta methods for Landau-Ginzburg-Higgs equation

Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic theory in the Hamilton space. The multi-symplectic Runge-Kutta method is reviewed, and a semi-implicit scheme with certain discrete conservation laws is constructed to solve the first-order partial differential equations （PDEs） derived from the Landau-Ginzburg-Higgs equation. The numerical re- sults for the soliton solution of the Landau-Ginzburg-Higgs equation are reported, showing that the multi-symplectic Runge-Kutta method is an efficient algorithm with excellent long-time numerical behaviors.     

Simulation of two- and three-dimensional internal subsonic flows using a finite element method

In this paper a finite element method is presented to predict internal subsonic flows. Using a low-Mach-number approximation, the pressure is decomposed into a mean thermodynamic contribution and a dynamic fluctuation to deal with the complex role of the pressure in internal aerodynamics. A semi-implicit time integration and a finite element method with a moving mesh are described to take into account complex geometries and moving boundaries. An Uzawa algorithm accelerated by a preconditioned residual method is introduced to solve the coupled non-symmetric linear system for the velocity components and the pressure. An efficient conjugate gradient method combined with an incomplete LU preconditioning is used to solve the non-symmetric linear systems arising from the discretization. The implementation of the numerical scheme on parallel supercomputers is also discussed. Efficient algorithms for the finite element assembly phase and for the solution of linear systems are described which take advantage of the parallel architecture of the new generation of supercomputers. With this technique a global speed-up of 10 is achieved on a supercomputer with eight processors. To illustrate the capabilities of the numerical method, 2D and 3D simulations of flows in the combustion chamber of a reciprocating engine and around the combustor dome of a gas turbine engine are presented.     

A well‐balanced weighted essentially non‐oscillatory scheme for pollutant transport in shallow water

In this paper, a well‐balanced finite difference weighted essentially non‐oscillatory scheme is presented for modeling transport and diffusion of pollutant in shallow water flows. The scheme balances exactly the flux gradients and the source terms. Extensive one‐dimensional and two‐dimensional numerical experiments on uniform and curvilinear meshes strongly suggest that high resolution results are achieved for both water depth and pollutant concentration. The scheme is efficient and robust and can be applied to practical numerical simulation of pollutant transport phenomena in shallow water flows. Copyright © 2012 John Wiley & Sons, Ltd.     

An efficient parallel high-order compact scheme for the 3D incompressible Navier–Stokes equations

This article provides a strategy for solving incompressible turbulent flows, which combines compact finite difference schemes and parallel computing. The numerical features of this solver are the semi-implicit time advancement, the staggered arrangement of the variables and the fourth-order compact scheme discretisation. This is the usual way for solving accurately turbulent incompressible flows. We propose a new strategy for solving the Helmholtz/Poisson equations based on a parallel 2d-pencil decomposition of the diagonalisation method. The compact scheme derivatives are computed with the parallel diagonal dominant (PDD) algorithm, which achieves good parallel performances by introducing a bounded numerical error. We provide a new analysis of its effect on the numerical accuracy and conservation features. Several numerical experiments, including two simulations of turbulent flows, demonstrate that the PDD algorithm maintains the accuracy and conservation features, while conserving a good parallel performance, up to 4096 cores.