Stability of a Cartesian grid projection method for zero Froude number shallow water flows

In this paper a Godunov-type projection method for computing approximate solutions of the zero Froude number (incompressible) shallow water equations is presented. It is second-order accurate and locally conserves height (mass) and momentum. To enforce the underlying divergence constraint on the velocity field, the predicted numerical fluxes, computed with a standard second order method for hyperbolic conservation laws and applied to an auxiliary system, are corrected in two steps. First, a MAC-type projection adjusts the advective velocity divergence. In a second projection step, additional momentum flux corrections are computed to obtain new time level cell-centered velocities, which satisfy another discrete version of the divergence constraint. The scheme features an exact and stable second projection. It is obtained by a Petrov-Galerkin finite element ansatz with piecewise bilinear trial functions for the unknown height and piecewise constant test\ud functions. The key innovation compared to existing finite volume projection methods is a correction of the in-cell slopes of the momentum by the second projection. The stability of the projection is proved using a generalized theory for mixed finite elements. In order to do so, the validity of three different inf-sup conditions has to be shown. The results of preliminary numerical test cases demonstrate the method's applicability. On fixed grids the accuracy is improved by a factor four compared to a previous version of the scheme

Stability of a Cartesian grid projection method for zero Froude number shallow water flows

In this paper a Godunov-type projection method for computing approximate solutions of the zero Froude number (incompressible) shallow water equations is presented. It is second-order accurate and locally conserves height (mass) and momentum. To enforce the underlying divergence constraint on the velocity field, the predicted numerical fluxes, computed with a standard second order method for hyperbolic conservation laws, are corrected in two steps. First, a MAC-type projection adjusts the advective velocity divergence. In a second projection step, additional momentum flux corrections are computed to obtain new time level cell-centered velocities, which satisfy another discrete version of the divergence constraint. The scheme features an exact and stable second projection. It is obtained by a Petrov-Galerkin finite element ansatz with piecewise bilinear trial functions for the unknown incompressible height and piecewise constant test functions. The stability of the projection is proved using the theory of generalized mixed finite elements, which goes back to Nicolaïdes (1982). In order to do so, the validity of three different inf-sup conditions has to be shown. Since the zero Froude number shallow water equations have the same mathematical structure as the incompressible Euler equations of isentropic gas dynamics, the method can be easily transfered to the computation of incompressible variable density flow problems

A New Projection Method for the Zero Froude Number Shallow Water Equations

For non-zero Froude numbers the shallow water equations are a hyperbolic system of partial differential equations. In the zero Froude number limit, they are of mixed hyperbolic-elliptic type, and the velocity field is subject to a divergence constraint. A new semi-implicit projection method for the zero Froude number shallow water equations is presented. This method enforces the divergence constraint on the velocity field, in two steps. First, the numerical fluxes of an auxiliary hyperbolic system are computed with a standard second order method. Then, these fluxes are corrected by solving two Poisson-type equations. These corrections guarantee that the new velocity field satisfies a discrete form of the above-mentioned divergence constraint. The main feature of the new method is a unified discretization of the two Poisson-type equations, which rests on a Petrov-Galerkin finite element formulation with piecewise bilinear ansatz functions for the unknown variable. This discretization naturally leads to piecewise linear ansatz functions for the momentum components. The projection method is derived from a semi-implicit finite volume method for the zero Mach number Euler equations, which uses standard discretizations for the solution of the Poisson-type equations. The new scheme can be formulated as an approximate as well as an exact projection method. In the former case, the divergence constraint is not exactly satisfied. The "approximateness" of the method can be estimated with an asymptotic upper bound of the velocity divergence at the new time level, which is consistent with the method's second-order accuracy. In the exact projection method, the piecewise linear components of the momentum are employed for the computation of the numerical fluxes of the auxiliary system at the new time level. In order to show the stability of the new projection step, a primal-dual mixed finite element formulation is derived, which is equivalent to the Poisson-type equations of the new scheme. Using the abstract theory of Nicola"ides (1982) for generalized saddle point problems, existence and uniqueness of the continuous problem are proven. Furthermore, preliminary results regarding the stability of the discrete method are presented. The numerical results obtained with the new exact method show significant accuracy improvements over the version that uses standard discretizations for the solution of the Poisson-type equations. In the L-two as well as the L-infinity norm, the global error is about four times smaller for smooth solutions. Simulating the advection of a vortex with discontinuous vorticity field, the new method yields a more accurate position of the center of the vortex

A Semi-Implicit Multiscale Scheme for Shallow Water Flows at Low Froude Number

A new large time step semi-implicit multiscale method is presented for the solution of low Froude-number shallow water flows. While on small scales which are under-resolved in time the impact of source terms on the divergence of the flow is essentially balanced, on large resolved scales the scheme propagates free gravity waves with minimized diffusion. The scheme features a scale decomposition based on multigrid ideas. Two different time integrators are blended at each scale depending on the scale-dependent Courant number for gravity wave propagation. The finite-volume discretization is based on a Cartesian grid and is second order accurate. The basic properties of the method are validated by numerical tests. This development is a further step in the development of asymptotically adaptive numerical methods for the computation of large scale atmospheric flows

Unstructured mesh methods for stratified turbulent flows

Developments are reported of unstructured-mesh methods for simulating stratified, turbulent and shear flows. The numerical model employs nonoscillatory forward in-time integrators for anelastic and incompressible flow PDEs, built on Multidimensional Positive Definite Advection Transport Algorithm (MPDATA) and a preconditioned conjugate residual elliptic solver. Finite-volume spatial discretisation adopts an edge-based data structure. Tetrahedral-based and hybrid-based median-dual options for unstructured meshes are developed, enabling flexible spatial resolution. Viscous laminar and detached eddy simulation (DES) flow solvers are developed based on the edge-based NFT MPDATA scheme. The built-in implicit large eddy simulation (ILES) capability of the NFT scheme is also employed and extended to fully unstructured tetrahedral and hybrid meshes. Challenging atmospheric and engineering problems are solved numerically to validate the model and to demonstrate its applications. The numerical problems include simulations of stratified, turbulent and shear flows past obstacles involving complex gravity-wave phenomena in the lee, critical-level laminar-turbulence transitioning and various vortex structures in the wake. Qualitative flow patterns and quantitative data analysis are both presented in the current study

Non-Linear Shallow Water Equations numerical integration on curvilinear boundary-conforming grids

An Upwind Weighted Essentially Non-Oscillatory scheme for the solution of the Shallow Water Equations on generalized curvilinear coordinate systems is proposed. The Shallow Water Equations are expressed in a contravariant formulation in which Christoffel symbols are avoided. The equations are solved by using a high-resolution finite-volume method incorporated with an exact Riemann Solver. A procedure developed in order to correct errors related to the difficulties of numerically satisfying the metric identities on generalized boundary-conforming grids is presented; this procedure allows the numerical scheme to satisfy the freestream preservation property on highly-distorted grids. The capacity of the proposed model is verified against test cases present in literature. The results obtained are compared with analytical solutions and alternative numerical solutions

Centered-potential regularization for the advection upstream splitting method

International audienceThis paper is devoted to a centered IMEX scheme in a multidimensional framework for a wide class of multicomponent and isentropic flows. The proposed strategy is based on a regularized model where the advection velocity is modified by the gradient of the potential of the conservative forces in both mass and momentum equations. The stability of the scheme is ensured by the dissipation of mechanic energy, which stands for a mathematical entropy, under an advective CFL condition. The main physical properties, such as positivity, conservation of the total momentum, and conservation of the steady state at rest, are satisfied. In addition, asymptotic preserving properties in the regimes (“incompressible” and “acoustic”) are analyzed. Finally, several simulations are presented to illustrate our results in a simplified context of oceanic flows in one dimension

On Simulation-based Ship Maneuvering Prediction in Deep and Shallow Water

A simulation-based framework for the prediction of ship maneuvering in deep and shallow water is presented. A mathematical model for maneuvering represented by coupled nonlinear differential equations stemming from Newtonian mechanics is derived. Hydrodynamic forces are modeled by multivariat polynomials, and therein included are coefficients representing ship-specific hydrodynamic properties which are determined by way of captive maneuvering tests using Computational Fluid Dynamics (CFD). The development and evaluation of efficacy of the proposed framework encompasses verification and validation studies on numerical methods for maneuvering and flows around ships in shallow water. The flow field information available from numerical simulations are used to discuss hydrodynamic phenomena related to viscous and free surface effects, as well as squat

Numerical Simulation of a Marine Current Turbine in Turbulent Flow

The copyright of this thesis rests with the author and no quotation from it or information derived from it may be published without the prior written consent of the authorThe marine current turbine (MCT) is an exciting proposition for the extraction of renewable tidal and marine current power. However, the numerical prediction of the performance of the MCT is difficult due to its complex geometry, the surrounding turbulent flow and the free surface. The main purpose of this research is to develop a computational tool for the simulation of a MCT in turbulent flow and in this thesis, the author has modified a 3D Large Eddy Simulation (LES) numerical code to simulate a three blade MCT under a variety of operating conditions based on the Immersed Boundary Method (IBM) and the Conservative Level Set Method (CLS). The interaction between the solid structure and surrounding fluid is modelled by the immersed boundary method, which the author modified to handle the complex geometrical conditions. The conservative free surface (CLS) scheme was implemented in the original Cgles code to capture the free surface effect. A series of simulations of turbulent flow in an open channel with different slope conditions were conducted using the modified free surface code. Supercritical flow with Froude number up to 1.94 was simulated and a decrease of the integral constant in the law of the wall has been noticed which matches well with the experimental data. Further simulations of the marine current turbine in turbulent flow have been carried out for different operating conditions and good match with experimental data was observed for all flow conditions. The effect of waves on the performance of the turbine was also investigated and it has been noticed that this existence will increase the power performance of the turbine due to the increase of free stream velocity