Exploring different possibilities for second‐order well‐balanced Lagrange‐projection numerical schemes applied to shallow water Exner equations

Author(s):  
C. Chalons ◽  
A. Del Grosso
2009 ◽  
Vol 137 (4) ◽  
pp. 1422-1437 ◽  
Author(s):  
Jin-Luen Lee ◽  
Alexander E. MacDonald

Abstract An icosahedral-hexagonal shallow-water model (SWM) on the sphere is formulated on a local Cartesian coordinate based on the general stereographic projection plane. It is discretized with the third-order Adam–Bashforth time-differencing scheme and the second-order finite-volume operators for spatial derivative terms. The finite-volume operators are applied to the model variables defined on the nonstaggered grid with the edge variables interpolated using polynomial interpolation. The projected local coordinate reduces the solution space from the three-dimensional, curved, spherical surface to the two-dimensional plane and thus reduces the number of complete sets of basis functions in the Vandermonde matrix, which is the essential component of the interpolation. The use of a local Cartesian coordinate also greatly simplifies the mathematic formulation of the finite-volume operators and leads to the finite-volume integration along straight lines on the plane, rather than along curved lines on the spherical surface. The SWM is evaluated with the standard test cases of Williamson et al. Numerical results show that the icosahedral SWM is free from Pole problems. The SWM is a second-order finite-volume model as shown by the truncation error convergence test. The lee-wave numerical solutions are compared and found to be very similar to the solutions shown in other SWMs. The SWM is stably integrated for several weeks without numerical dissipation using the wavenumber 4 Rossby–Haurwitz solution as an initial condition. It is also shown that the icosahedral SWM achieves mass conservation within round-off errors as one would expect from a finite-volume model.


Author(s):  
Carl Trygve Stansberg ◽  
Trygve Kristiansen

Slowly varying motions and drift forces of a large moored ship in random waves at 35m water depth are investigated by an experimental wave basin study in scale 1:50. A simple horizontal mooring set-up is used. A second-order wave correction is applied to minimize “parasitic” long waves. The effect on the ship motion from the correction is clearly seen, although less in random wave spectra than in pure bi-chromatic waves. Empirical quadratic transfer functions (QTFs) of the surge drift force are found by use of cross-bi-spectral analysis, in two different spectra have been obtained. The QTF levels increase significantly with lower wave frequencies (except at the diagonal), which is special for finite and shallow water. Furthermore, the QTF levels frequencies at low frequencies increase significantly out from the QTF diagonal. Thus Newman’s approximation should preferrably not be used in these cases. Using the LF waves as a direct excitation in a “linear” ship force analysis gives random records that compare reasonably well with those from the cross-bi-spectral analysis. This confirms the idea that the drift forces in shallow water are closely correlated to the second-order potential, and thereby by the second-order LF waves.


2017 ◽  
Vol 145 (3) ◽  
pp. 751-772 ◽  
Author(s):  
Michael D. Toy ◽  
Ramachandran D. Nair

An energy and potential enstrophy conserving finite-difference scheme for the shallow-water equations is derived in generalized curvilinear coordinates. This is an extension of a scheme formulated by Arakawa and Lamb for orthogonal coordinate systems. The starting point for the present scheme is the shallow-water equations cast in generalized curvilinear coordinates, and tensor analysis is used to derive the invariant conservation properties. Preliminary tests on a flat plane with doubly periodic boundary conditions are presented. The scheme is shown to possess similar order-of-convergence error characteristics using a nonorthogonal coordinate compared to Cartesian coordinates for a nonlinear test of flow over an isolated mountain. A linear normal mode analysis shows that the discrete form of the Coriolis term provides stationary geostrophically balanced modes for the nonorthogonal coordinate and no unphysical computational modes are introduced. The scheme uses centered differences and averages, which are formally second-order accurate. An empirical test with a steady geostrophically balanced flow shows that the convergence rate of the truncation errors of the discrete operators is second order. The next step will be to adapt the scheme for use on the cubed sphere, which will involve modification at the lateral boundaries of the cube faces.


Author(s):  
Jo̸rgen Juncher Jensen

For bottom-supported offshore structures like oil drilling rigs and oil production platforms, a deterministic design wave approach is often applied using a regular non-linear Stokes’ wave. Thereby, the procedure accounts for non-linear effects in the wave loading but the randomness of the ocean waves is poorly represented, as the shape of the wave spectrum does not enter the wave kinematics. To overcome this problem and still keep the simplicity of a deterministic approach, Tromans, Anaturk and Hagemeijer (1991) suggested the use of a deterministic wave, defined as the expected linear Airy wave, given the value of the wave crest at a specific point in time or space. In the present paper a derivation of the expected second order short-crested wave riding on a uniform current is given. The analysis is based on the second order Sharma and Dean shallow water wave theory and the direction of the main wind direction can make any direction with the current. Numerical results showing the importance of the water depth, the directional spreading and the current on the conditional mean wave profile and the associated wave kinematics are presented. A discussion of the use of the conditional wave approach as design waves is given.


2014 ◽  
Vol 2014 ◽  
pp. 1-20 ◽  
Author(s):  
Sheng Bi ◽  
Jianzhong Zhou ◽  
Yi Liu ◽  
Lixiang Song

A second-order accurate, Godunov-type upwind finite volume method on dynamic refinement grids is developed in this paper for solving shallow-water equations. The advantage of this grid system is that no data structure is needed to store the neighbor information, since neighbors are directly specified by simple algebraic relationships. The key ingredient of the scheme is the use of the prebalanced shallow-water equations together with a simple but effective method to track the wet/dry fronts. In addition, a second-order spatial accuracy in space and time is achieved using a two-step unsplit MUSCL-Hancock method and a weighted surface-depth gradient method (WSDM) which considers the local Froude number is proposed for water depths reconstruction. The friction terms are solved by a semi-implicit scheme that can effectively prevent computational instability from small depths and does not invert the direction of velocity components. Several benchmark tests and a dam-break flooding simulation over real topography cases are used for model testing and validation. Results show that the proposed model is accurate and robust and has advantages when it is applied to simulate flow with local complex topographic features or flow conditions and thus has bright prospects of field-scale application.


Sign in / Sign up

Export Citation Format

Share Document