A Petrov–Galerkin finite element model for one‐dimensional fully non‐linear and weakly dispersive wave propagation

A new finite element method is presented to solve one‐dimensional depth‐integrated equations for fully non‐linear and weakly dispersive waves. For spatial integration, the Petrov–Galerkin weighted residual method is used. The weak forms of the governing equations are arranged in such a way that the shape functions can be piecewise linear, while the weighting functions are piecewise cubic with C²‐continuity. For the time integration an implicit predictor–corrector iterative scheme is employed. Within the framework of linear theory, the accuracy of the scheme is discussed by considering the truncation error at a node. The leading truncation error is fourth‐order in terms of element size. Numerical stability of the scheme is also investigated. If the Courant number is less than 0.5, the scheme is unconditionally stable. By increasing the number of iterations and/or decreasing the element size, the stability characteristics are improved significantly. Both Dirichlet boundary condition (for incident waves) and Neumann boundary condition (for a reflecting wall) are implemented. Several examples are presented to demonstrate the range of applicabilities and the accuracy of the model. Copyright © 2001 John Wiley & Sons, Ltd.     

A fully non‐linear model for three‐dimensional overturning waves over an arbitrary bottom

An accurate three‐dimensional numerical model, applicable to strongly non‐linear waves, is proposed. The model solves fully non‐linear potential flow equations with a free surface using a higher‐order three‐dimensional boundary element method (BEM) and a mixed Eulerian–Lagrangian time updating, based on second‐order explicit Taylor series expansions with adaptive time steps. The model is applicable to non‐linear wave transformations from deep to shallow water over complex bottom topography up to overturning and breaking. Arbitrary waves can be generated in the model, and reflective or absorbing boundary conditions specified on lateral boundaries. In the BEM, boundary geometry and field variables are represented by 16‐node cubic ‘sliding’ quadrilateral elements, providing local inter‐element continuity of the first and second derivatives. Accurate and efficient numerical integrations are developed for these elements. Discretized boundary conditions at intersections (corner/edges) between the free surface or the bottom and lateral boundaries are well‐posed in all cases of mixed boundary conditions. Higher‐order tangential derivatives, required for the time updating, are calculated in a local curvilinear co‐ordinate system, using 25‐node ‘sliding’ fourth‐order quadrilateral elements. Very high accuracy is achieved in the model for mass and energy conservation. No smoothing of the solution is required, but regridding to a higher resolution can be specified at any time over selected areas of the free surface. Applications are presented for the propagation of numerically exact solitary waves. Model properties of accuracy and convergence with a refined spatio‐temporal discretization are assessed by propagating such a wave over constant depth. The shoaling of solitary waves up to overturning is then calculated over a 1:15 plane slope, and results show good agreement with a two‐dimensional solution proposed earlier. Finally, three‐dimensional overturning waves are generated over a 1:15 sloping bottom having a ridge in the middle, thus focusing wave energy. The node regridding method is used to refine the discretization around the overturning wave. Convergence of the solution with grid size is also verified for this case. Copyright © 2001 John Wiley & Sons, Ltd.     

A stabilization for three‐dimensional discontinuous Galerkin discretizations applied to nonhydrostatic atmospheric simulations

A discontinuous Galerkin nonhydrostatic atmospheric model is used for two‐dimensional and three‐dimensional simulations. There is a wide range of timescales to be dealt with. To do so, two different implicit/explicit time discretizations are implemented. A stabilization, based upon a reduced‐order discretization of the gravity term, is introduced to ensure the balance between pressure and gravity effects. While not affecting significantly the convergence properties of the scheme, this approach allows the simulation of anisotropic flows without generating spurious oscillations, as it happens for a classical discontinuous Galerkin discretization. This approach is shown to be less diffusive than usual spatial filters. A stability analysis demonstrates that the use of this modified scheme discards the instability associated with the usual discretization. Validation against analytical solutions is performed, confirming the good convergence and stability properties of the scheme. Numerical results demonstrate the attractivity of the discontinuous Galerkin method with implicit/explicit time integration for large‐scale atmospheric flows. Copyright © 2015 John Wiley & Sons, Ltd.     

The anti‐dissipative,non‐monotone behavior of Petrov–Galerkin upwinding

The Petrov–Galerkin method has been developed with the primary goal of damping spurious oscillations near discontinuities in advection dominated flows. For time‐dependent problems, the typical Petrov–Galerkin method is based on the minimization of the dispersion error and the simultaneous selective addition of dissipation. This optimal design helps to dampen the oscillations prevalent near discontinuities in standard Bubnov–Galerkin solutions. However, it is demonstrated that when the Courant number is less than 1, the Petrov–Galerkin method actually amplifies undershoots at the base of discontinuities. This is shown in an heuristic manner, and is demonstrated with numerical experiments with the scalar advection and Richards' equations. A discussion of monotonicity preservation as a design criterion, as opposed to phase or amplitude error minimization, is also presented. The Petrov–Galerkin method is further linked to the high‐resolution, total variation diminishing (TVD) finite volume method in order to obtain a monotonicity preserving Petrov–Galerkin method.     

Simulation of non‐linear free surface motions in a cylindrical domain using a Chebyshev–Fourier spectral collocation method

When a liquid is perturbed, its free surface may experience highly non‐linear motions in response. This paper presents a numerical model of the three‐dimensional hydrodynamics of an inviscid liquid with a free surface. The mathematical model is based on potential theory in cylindrical co‐ordinates with a σ‐transformation applied between the bed and free surface in the vertical direction. Chebyshev spectral elements discretize space in the vertical and radial directions; Fourier spectral elements are used in the angular direction. Higher derivatives are approximated using a collocation (or pseudo‐spectral) matrix method. The numerical scheme is validated for non‐linear transient sloshing waves in a cylindrical tank containing a circular surface‐piercing cylinder at its centre. Excellent agreement is obtained with Ma and Wu's [Second order transient waves around a vertical cylinder in a tank. Journal of Hydrodynamics 1995; Ser. B4 : 72–81] second‐order potential theory. Further evidence for the capability of the scheme to predict complicated three‐dimensional, and highly non‐linear, free surface motions is given by the evolution of an impulse wave in a cylindrical tank and in an open domain. Copyright © 2001 John Wiley & Sons, Ltd.     

A discontinuous Galerkin‐based sharp‐interface method to simulate three‐dimensional compressible two‐phase flow

A numerical method for the simulation of compressible two‐phase flows is presented in this paper. The sharp‐interface approach consists of several components: a discontinuous Galerkin solver for compressible fluid flow, a level‐set tracking algorithm to follow the movement of the interface and a coupling of both by a ghost‐fluid approach with use of a local Riemann solver at the interface. There are several novel techniques used: the discontinuous Galerkin scheme allows locally a subcell resolution to enhance the interface resolution and an interior finite volume Total Variation Diminishing (TVD) approximation at the interface. The level‐set equation is solved by the same discontinuous Galerkin scheme. To obtain a very good approximation of the interface curvature, the accuracy of the level‐set field is improved and smoothed by an additional P_NP_M‐reconstruction. The capabilities of the method for the simulation of compressible two‐phase flow are demonstrated for a droplet at equilibrium, an oscillating ellipsoidal droplet, and a shock‐droplet interaction problem at Mach 3. Copyright © 2015 John Wiley & Sons, Ltd.     

On use of the Galerkin method to solve the fracture mechanics problem for a linear crack under normal loading

The problem of contact interaction of the opposite faces of a linear crack under a normally incident harmonic tension-compression wave is numerically solved by the Galerkin method with piecewise-linear continuous elements. The dependence of the stress intensity factor (opening mode) on the wave number is investigated Published in Prikladnaya Mekhanika, Vol. 41, No. 11, pp. 137–142, November 2005.     

An implicit three‐dimensional fully non‐hydrostatic model for free‐surface flows

An implicit method is developed for solving the complete three‐dimensional (3D) Navier–Stokes equations. The algorithm is based upon a staggered finite difference Crank‐Nicholson scheme on a Cartesian grid. A new top‐layer pressure treatment and a partial cell bottom treatment are introduced so that the 3D model is fully non‐hydrostatic and is free of any hydrostatic assumption. A domain decomposition method is used to segregate the resulting 3D matrix system into a series of two‐dimensional vertical plane problems, for each of which a block tri‐diagonal system can be directly solved for the unknown horizontal velocity. Numerical tests including linear standing waves, nonlinear sloshing motions, and progressive wave interactions with uneven bottoms are performed. It is found that the model is capable to simulate accurately a range of free‐surface flow problems using a very small number of vertical layers (e.g. two–four layers). The developed model is second‐order accuracy in time and space and is unconditionally stable; and it can be effectively used to model 3D surface wave motions. Copyright © 2004 John Wiley & Sons, Ltd.     

A non‐linear coupled fluid–structure aeroelastic analysis of a slender bridge deck

This paper presents a new numerical methodology for the integral aeroelastic instability analysis of slender structures, based on the appropriate conjugation of an algorithm for dynamic and geometrically non‐linear analysis of structures based on the finite element method with another algorithm of computational fluid dynamics (Finite volume method). It is considered a viscous incompressible unsteady turbulent bidimensional air flow solved on a structured control volume mesh. The computer code developed on the basis of this methodology is applied to the aeroelastic study of a simply supported slender bridge deck in order to find out the critical wind velocity leading to instability. Some of the most significant results associated with the analysis of the corresponding aeroelastic behaviour are presented. Copyright 2004 John Wiley & Sons, Ltd.     

A novel implicit numerical treatment for non‐linear turbulence models using high and low Reynolds number formulations

Non‐linear turbulence models can be seen as an improvement of the classical eddy‐viscosity concept due to their better capacity to simulate characteristics of important flows. However, application of non‐linear models demand robustness of the numerical method applied, requiring a stable discretization scheme for convergence of all variables involved. Usually, non‐linear terms are handled in an explicit manner leading to possible numerical instabilities. Thus, the present work shows the steps taken to adapt a general non‐linear constitutive equation using a new semi‐implicit numerical treatment for the non‐linear diffusion terms. The objective is to increase the degree of implicitness of the solution algorithm to enhance convergence characteristics. Flow over a backward‐facing step was computed using the control volume method applied to a boundary‐fitted coordinate system. The SIMPLE algorithm was used to relax the algebraic equations. Classical wall function and a low Reynolds number model were employed to describe the flow near the wall. The results showed that for certain combination of relaxation parameters, the semi‐implicit treatment proposed here was the sole successful treatment in order to achieve solution convergence. Also, application of the implicit method described here shows that the stability of the solution either increases (high Reynolds with non‐orthogonal mesh) or preserves the same (low Reynolds number applications). Additional advantages of the procedure proposed here lie in the possibility of testing different non‐linear expressions if one considers the enhanced robustness and stability obtained for the entire numerical algorithm. Copyright © 2010 John Wiley & Sons, Ltd.     

An hp‐adaptive discontinuous Galerkin method for shallow water flows

An adaptive spectral/hp discontinuous Galerkin method for the two‐dimensional shallow water equations is presented. The model uses an orthogonal modal basis of arbitrary polynomial order p defined on unstructured, possibly non‐conforming, triangular elements for the spatial discretization. Based on a simple error indicator constructed by the solutions of approximation order p and p?1, we allow both for the mesh size, h, and polynomial approximation order to dynamically change during the simulation. For the h‐type refinement, the parent element is subdivided into four similar sibling elements. The time‐stepping is performed using a third‐order Runge–Kutta scheme. The performance of the hp‐adaptivity is illustrated for several test cases. It is found that for the case of smooth flows, p‐adaptivity is more efficient than h‐adaptivity with respect to degrees of freedom and computational time. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulation of three‐dimensional bubbles using desingularized boundary integral method

A very simple model based on the three‐dimensional desingularized boundary integral method is applied to study the evolution of bubble(s) with or without the presence of solid structures. The choice of the desingularization parameters, which is crucial to the success of the method, is studied in the context of bubble dynamics. With the proper choice of parameters, the new model is far more efficient than previous models with virtually the same level of accuracy being achieved. This is largely attributed to the simplicity of the desingularization method. Furthermore, the new model offers a simple and attractive way for mesh refinement. Although it has limitations in the sense that, with this model the time stepping tends to slow down as two surfaces approach each other, this can be easily rectified by switching over to a direct method so that the two surfaces can be drawn closer as required in the context of jet impact. After this the new model can be reinstated to treat the complicated doubly connected geometry involving toroidal bubbles that would otherwise be very difficult to deal with. Copyright © 1999 John Wiley & Sons, Ltd.     

A local domain‐free discretization method to simulate three‐dimensional compressible inviscid flows

In this paper, the domain‐free discretization method (DFD) is extended to simulate the three‐dimensional compressible inviscid flows governed by Euler equations. The discretization strategy of DFD is that the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior‐dependent points are updated at each time step by extrapolation along the wall normal direction in conjunction with the wall boundary conditions and the simplified momentum equation in the vicinity of the wall. Spatial discretization is achieved with the help of the finite element Galerkin approximation. The concept of ‘osculating plane’ is adopted, with which the local DFD can be easily implemented for the three‐dimensional case. Geometry‐adaptive tetrahedral mesh is employed for three‐dimensional calculations. Finally, we validate the DFD method for three‐dimensional compressible inviscid flow simulations by computing transonic flows over the ONERA M6 wing. Comparison with the reference experimental data and numerical results on boundary‐conforming grid was displayed and the results show that the present DFD results compare very well with the reference data. Copyright © 2009 John Wiley & Sons, Ltd.     

Slope limiting for discontinuous Galerkin approximations with a possibly non‐orthogonal Taylor basis

The use of high‐order polynomials in discontinuous Galerkin (DG) approximations to convection‐dominated transport problems tends to cause a violation of the maximum principle in regions where the derivatives of the solution are large. In this paper, we express the DG solution in terms of Taylor basis functions associated with the cell average and derivatives at the center of the cell. To control the (derivatives of the) discontinuous solution, the values at the vertices of each element are required to be bounded by the means. This constraint is enforced using a hierarchical vertex‐based slope limiter to constrain the coefficients of the Taylor polynomial in a conservative manner starting with the highest‐order terms. The loss of accuracy at smooth extrema is avoided by taking the maximum of the correction factors for derivatives of order p and higher. No free parameters, oscillation detectors, or troubled cell markers are involved. In the case of a non‐orthogonal Taylor basis, the same limiter is applied to the vector of discretized time derivatives before the multiplication by the off‐diagonal part of the consistent mass matrix. This strategy leads to a remarkable gain of accuracy, especially in the case of simplex meshes. A numerical study is performed for a 2D convection equation discretized with linear and quadratic finite elements. Copyright © 2012 John Wiley & Sons, Ltd.     

A high‐order Petrov–Galerkin finite element method for the classical Boussinesq wave model

A high‐order Petrov–Galerkin finite element scheme is presented to solve the one‐dimensional depth‐integrated classical Boussinesq equations for weakly non‐linear and weakly dispersive waves. Finite elements are used both in the space and the time domains. The shape functions are bilinear in space–time, whereas the weighting functions are linear in space and quadratic in time, with C⁰‐continuity. Dispersion correction and a highly selective dissipation mechanism are introduced through additional streamline upwind terms in the weighting functions. An implicit, conditionally stable, one‐step predictor–corrector time integration scheme results. The accuracy and stability of the non‐linear discrete equations are investigated by means of a local Taylor series expansion. A linear spectral analysis is used for the full characterization of the predictor–corrector inner iterations. Based on the order of the analytical terms of the Boussinesq model and on the order of the numerical discretization, it is concluded that the scheme is fourth‐order accurate in terms of phase velocity. The dissipation term is third order only affecting the shortest wavelengths. A numerical convergence analysis showed a second‐order convergence rate in terms of both element size and time step. Four numerical experiments are addressed and their results are compared with analytical solutions or experimental data available in the literature: the propagation of a solitary wave, the oscillation of a flat bottom closed basin, the oscillation of a non‐flat bottom closed basin, and the propagation of a periodic wave over a submerged bar. Copyright © 2008 John Wiley & Sons, Ltd.     

Boussinesq cut‐cell model for non‐linear wave interaction with coastal structures

Boussinesq models describe the phase‐resolved hydrodynamics of unbroken waves and wave‐induced currents in shallow coastal waters. Many enhanced versions of the Boussinesq equations are available in the literature, aiming to improve the representation of linear dispersion and non‐linearity. This paper describes the numerical solution of the extended Boussinesq equations derived by Madsen and Sørensen (Coastal Eng. 1992; 15 :371–388) on Cartesian cut‐cell grids, the aim being to model non‐linear wave interaction with coastal structures. An explicit second‐order MUSCL‐Hancock Godunov‐type finite volume scheme is used to solve the non‐linear and weakly dispersive Boussinesq‐type equations. Interface fluxes are evaluated using an HLLC approximate Riemann solver. A ghost‐cell immersed boundary method is used to update flow information in the smallest cut cells and overcome the time step restriction that would otherwise apply. The model is validated for solitary wave reflection from a vertical wall, diffraction of a solitary wave by a truncated barrier, and solitary wave scattering and diffraction from a vertical circular cylinder. In all cases, the model gives satisfactory predictions in comparison with the published analytical solutions and experimental measurements. Copyright © 2007 John Wiley & Sons, Ltd.     

A finite element method for two‐dimensional water‐wave problems

In this paper, the authors treat the free‐surface waves generated by a moving disturbance with a constant speed in water of finite and constant depth. Specifically, the case when the disturbance is moving with the critical speed is investigated. The water is assumed inviscid and its motion irrotational. The surface tension is neglected. It is well‐known that the linear theory breaks down when a disturbance is moving with the critical speed. As a remedy to overcome the invalid linear theory, approximate non‐linear theories have been applied with success in the past, i.e. Boussinesq and Korteweg de Vries equations, for example. In the present paper, the authors describe a finite element method applied to the non‐linear water‐wave problems in two dimensions. The present numerical method solves the exact non‐linear formulation in the scope of potential theory without any additional assumptions on the magnitude of the disturbances. The present numerical results are compared with those obtained by other approximate non‐linear theories. Also presented are the discussions on the validity of the existing approximate theories applied to two types of the disturbances, i.e. the bottom bump and the pressure patch on the free‐surface at the critical speed. Copyright © 1999 John Wiley & Sons, Ltd.     

A numerical method for the determination of weakly non‐hydrostatic non‐linear free surface wave propagation

A numerical method is described that may be used to determine the propagation characteristics of weakly non‐hydrostatic non‐linear free surface waves over a general, bottom topography. In shallow water of constant undisturbed depth, such waves are equivalent to the familiar cnoidal waves characterized by sharp crests and relatively flat troughs. For a certain range of parameters, these propagate without change of form by virtue of the weakly non‐hydrostatic balance in the vertical momentum equation. Effectively, this counters the tendency for the non‐linearity in a purely hydrostatic theory to lead to a continuously deforming surface wave profile. The realistic representation furnished by cnoidal wave theory of free surface waves in the shallow near‐shore zone has led to its utilization in evaluating their propagation characteristics. Nonetheless, the classic analytical theory is inapplicable to the case of wave propagation over a sloping beach or off‐shore sand bar topography. Under these conditions, a local change in form of the surface wave profile is anticipated before the waves break and knowing this is required in order to evaluate fully the propagation process. The efficacy of the numerical method is first demonstrated by comparing the solution for water of constant depth with the evaluation of the analytical solution expressed in terms of the Jacobian elliptic function cn. The general method described in the paper is then illustrated by experiments to determine the change in profile of weakly non‐hydrostatic non‐linear surface waves propagating over bed forms representative of those found in shallow coastal seas. Copyright © 2001 John Wiley & Sons, Ltd.