Analysis of the Nonlinear Shallow Water Equations Over Nonplanar Topography

The role played by the beach bottom profile on coastal inundation phenomena is analyzed here by means of approximate analytical solutions of the nonlinear shallow water equations (NSWEs) over uneven bottoms. These are obtained by only using the assumptions of small waves at the seaward boundary and small topographic forcing. Our work, built on the Carrier and Greenspan [ 1 ] hodographic transformation and on the solution of the boundary value problem (BVP) for the NSWEs proposed by Antuono and Brocchini [ 2 ], focuses on the propagation of nonlinear non-breaking waves over quasi-planar beaches. Since the terms associated with the perturbed bottom only appear in the second-order perturbed solutions, the breaking conditions for the planar-beach bathymetry also predict well the breaking occurring on the nonplanar beaches analyzed here. The most important results, concerning the shoreline position and the near-shoreline velocity, are given for both pulse-like and periodic input waves propagating over two types of nonplanar bathymetries. The solution proposed here is a fundamental benchmark for any numerical and theoretical analyzes concerned with estimates of wave run-up on beaches of complex shape.     

An elastodynamic Galerkin Boundary Element Formulation for semi-infinite domains in time-domain

The present work focuses on the problem of modelling wave propagation phenomena within a 3–d elastodynamic halfspace by use of a symmetric Galerkin Boundary Element formulation. Unfortunately, this formulation requires the evaluation of hypersingular integral kernels which are regularized by integration by parts. In Boundary Element Methods semi–infinite domains are commonly approximated in space by considering just a sufficiently large enough region. Applying this simple discretization to the symmetric formulation implies the evaluation of the hypersingular bilinear form on a truncated mesh which will fail due to the regularization approach. To overcome this drawback a methodology based on infinite elements is presented. The numerical tests show that this approach is promising for treating semi–infinite domains with a symmetric Galerkin scheme. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Legendre-Laguerre coupled spectral element methods for second- and fourth-order equations on the half line

Some Legendre spectral element/Laguerre spectral coupled methods are proposed to numerically solve second- and fourth-order equations on the half line. The proposed methods are based on splitting the infinite domain into two parts, then using the Legendre spectral element method in the finite subdomain and Laguerre method in the infinite subdomain. C⁰ or C¹-continuity, according to the problem under consideration, is imposed to couple the two methods. Rigorous error analysis is carried out to establish the convergence of the method. More importantly, an efficient computational process is introduced to solve the discrete system. Several numerical examples are provided to confirm the theoretical results and the efficiency of the method.     

Three‐Dimensional Surface Water Waves Governed by the Forced Benney–Luke Equation

This paper studies the propagation of three‐dimensional surface waves in water with an ambient current over a varying bathymetry. When the ambient flow is near the critical speed, under the shallow water assumptions, a forced Benney–Luke (fBL) equation is derived from the Euler equations. An asymptotic approximation of the water's reaction force over the varying bathymetry is derived in terms of topographic stress. Numerical simulations of the fBL equation over a trough are compared to those using a forced Kadomtsev–Petviashvilli equation. For larger variations in the bathymetry that upstream‐radiating three‐dimensional solitons are observed, which are different from the upstream‐radiating solitons simulated by the forced Kadomtsev–Petviashvilli equation. In this case, we show the fBL equation is a singular perturbation of the forced Kadomtsev–Petviashvilli equation which explains the significant differences between the two flows.     

Interfacial long traveling waves in a two‐layer fluid with variable depth

Long wave propagation in a two‐layer fluid with variable depth is studied for specific bottom configurations, which allow waves to propagate over large distances. Such configurations are found within the linear shallow‐water theory and determined by a family of solutions of the second‐order ordinary differential equation (ODE) with three arbitrary constants. These solutions can be used to approximate the true bottom bathymetry. All such solutions represent smooth bottom profiles between two different singular points. The first singular point corresponds to the point where the two‐layer flow transforms into a uniform one. In the vicinity of this point nonlinear shallow‐water theory is used and the wave breaking criterion, which corresponds to the gradient catastrophe is found. The second bifurcation point corresponds to an infinite increase in water depth, which contradicts the shallow‐water assumption. This point is eliminated by matching the “nonreflecting” bottom profile with a flat bottom. The wave transformation at the matching point is described by the second‐order Fredholm equation and its approximated solution is then obtained. The results extend the theory of internal waves in inhomogeneous stratified fluids actively developed by Prof. Roger Grimshaw, to the new solutions types.     

Application of the least-squares spectral element method using Chebyshev polynomials to solve the incompressible Navier-Stokes equations

In this paper the extension of the Legendre least-squares spectral element formulation to Chebyshev polynomials will be explained. The new method will be applied to the incompressible Navier-Stokes equations and numerical results, obtained for the lid-driven cavity flow at Reynolds numbers varying between 1000 and 7500, will be compared with the commonly used benchmark results. The new results reveal that the least-squares spectral element formulations based on the Legendre and Chebyshev Gauss-Lobatto Lagrange interpolating polynomials are equally accurate.     

Infinite elements in a poroelastodynamic FEM

Wave propagation phenomena in unbounded domains occur in many engineering applications, e.g., soil structure interactions. When simulating unbounded domains, infinite elements are a possible choice to describe the far field behavior, whereas the near field is described through conventional finite elements. Finite element formulations for porous materials in terms of Biot's theory [1] have been published, e.g., by Zienkiewicz [2]. For infinite elements, several approaches are described in [3,4]. Infinite elements are based on special shape functions to approximate the semi-infinite geometry as well as the Sommerfeld radiation condition, i.e., the waves decay with distance and are not reflected at infinity. If there is only one wave traveling in the media, a formulation in time-domain can be established. But in poroelastodynamics, there are three body waves and eventually also a Rayleigh wave. Unfortunately, the extension to more than one wave is not straight forward. Here, an infinite element is presented which can handle all wave types, as it is needed in poroelasticity. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A boundary parametric approximation to the linearized scalar potential magnetostatic field problem

We consider the linearized scalar potential formulation of the magnetostatic field problem in this paper. Our approach involves a reformulation of the continuous problem as a parametric boundary problem. By the introduction of a spherical interface and the use of spherical harmonics, the infinite boundary conditions can also be satisfied in the parametric framework. That is the field in the exterior of a sphere is expanded in a ‘harmonic series’ of eigenfunctions for the exterior harmonic problem. The approach is essentially a finite element method coupled with a spectral method via a boundary parametric procedure. The reformulated problem is discretized by finite element techniques which leads to a discrete parametric problem which can be solved by well conditioned iteration involving only the solution of decoupled Neumann type elliptic finite element systems and L² projection onto subspaces of spherical harmonics. Error and stability estimates given show exponential convergence in the degree of the spherical harmonics and optimal order convergence with respect to the finite element approximation for the resulting fields in L².     

The ultra weak variational formulation for the modified mild-slope equation

This paper is devoted to evaluate a new alternative numerical method to capture accurately the diffraction-refraction process of waves in coastal areas. The modified mild-slope equation (MMSE) has been used to predict the water wave transformation when waves approach the shoreline. We perform numerical simulations in order to illustrate the efficiency of the ultra weak variational formulation (UWVF) method in comparison with the finite elements method (FEM). The UWVF method uses plane wave solutions on each element and has been shown to reduce the computational complexity at high wave numbers. We also present an alternative method to seek the angle of attack of the wave front on the domain boundary and show that the UWVF method reproduces effectively the numerical experimental data. We compare the FEM and the UWVF method for the MMSE and show that the UWVF method solves some of the difficulties that arise when the FEM is applied.     

一个二流体系统中非线性水波的Hamilton描述

讨论了一个二流体系统中非线性水波的Hamilton描述,该系统由水平固壁之上的两层常密度不可压无粘流体组成,上表面为自由面.文中将速度势函数展开成垂向坐标的幂级数,在浅水长波的假定下,取下层流体的“动厚度”与上层流体的“折合动厚度”为广义位移、界面上和自由面上的速度势为广义动量,根据Hamilton原理并运用Legendre变换导出该系统的Hamilton正则方程,从而将单层流体情形的结果推广到分层流体的情形.     

Mixed Laguerre–Legendre Spectral Method for Incompressible Flow in an Infinite Strip

This paper concerns the mixed Laguerre–Legendre spectral approximation and its application to numerical simulation of incompressible flow in an infinite strip. Some approximation results in weighted Sobolev spaces are given. A Laguerre–Legendre spectral scheme for the stream function form of Navier–Stokes equations is constructed. The stability and the convergence of the proposed scheme are proved. The numerical experiments show the high accuracy of this method. The main techniques used in this paper are also applicable to other nonlinear partial differential equations in an infinite strip.     

A Legendre spectral element method on a large spatial domain to solve the predator–prey system modeling interacting populations

A Legendre spectral element method is developed for solving a one-dimensional predator–prey system on a large spatial domain. The predator–prey system is numerically solved where the prey population growth is described by a cubic polynomial and the predator’s functional response is Holling type I. The discretization error generated from this method is compared with the error obtained from the Legendre pseudospectral and finite element methods. The Legendre spectral element method is also presented where the predator response is Holling type II and the initial data are discontinuous.     

Space‐time spectral method for two‐dimensional semilinear parabolic equations

In this paper, a high‐order accurate numerical method for two‐dimensional semilinear parabolic equations is presented. We apply a Galerkin–Legendre spectral method for discretizing spatial derivatives and a spectral collocation method for the time integration of the resulting nonlinear system of ordinary differential equations. Our formulation can be made arbitrarily high‐order accurate in both space and time. Optimal a priori error bound is derived in the L²‐norm for the semidiscrete formulation. Extensive numerical results are presented to demonstrate the convergence property of the method, show our formulation have spectrally accurate in both space and time. John Wiley & Sons, Ltd.     

A Posteriori Error Estimation of Spectral and Spectral Element Methods for the Stokes/Darcy Coupled Problem

In this paper, we carry out an a posteriori error analysis of Legendre spectral approximations to the Stokes/Darcy coupled equations. The spectral approximations are based on a weak formulation of the coupled equations by using the Beavers-Joseph-Saffman interface condition. The main contribution of the paper consists of deriving a number of posteriori error indicators and their upper and lower bounds for the single domain case. An extension of the upper bounds to the multi-domain case in the spectral element framework is also given.     

Wave Packet Defocusing Due to a Highly Disordered Bathymetry

Slowly modulated water waves are considered in the presence of a strongly disordered bathymetry. Previous work is extended to the case where the random bottom irregularities are not smooth and are allowed to be of large amplitude. Through the combination of a conformal mapping and a multiple‐scales asymptotic analysis it is shown that large variations of a disordered bathymetry can affect the nonlinearity coefficient of the resulting damped nonlinear Schrödinger equations. In particular it is shown that as the bathymetry fluctuation level increases the critical point (separating the focusing from the defocusing region) moves to the right, hence enlarging the region where the dynamics is of a defocusing character.     

A spectral element method using the modal basis and its application in solving second‐order nonlinear partial differential equations

We present a high‐order spectral element method (SEM) using modal (or hierarchical) basis for modeling of some nonlinear second‐order partial differential equations in two‐dimensional spatial space. The discretization is based on the conforming spectral element technique in space and the semi‐implicit or the explicit finite difference formula in time. Unlike the nodal SEM, which is based on the Lagrange polynomials associated with the Gauss–Lobatto–Legendre or Chebyshev quadrature nodes, the Lobatto polynomials are used in this paper as modal basis. Using modal bases due to their orthogonal properties enables us to exactly obtain the elemental matrices provided that the element‐wise mapping has the constant Jacobian. The difficulty of implementation of modal approximations for nonlinear problems is treated in this paper by expanding the nonlinear terms in the weak form of differential equations in terms of the Lobatto polynomials on each element using the fast Fourier transform (FFT). Utilization of the Fourier interpolation on equidistant points in the FFT algorithm and the enough polynomial order of approximation of the nonlinear terms can lead to minimize the aliasing error. Also, this approach leads to finding numerical solution of a nonlinear differential equation through solving a system of linear algebraic equations. Numerical results for some famous nonlinear equations illustrate efficiency, stability and convergence properties of the approximation scheme, which is exponential in space and up to third‐order in time. Copyright © 2014 John Wiley & Sons, Ltd.     

The boundary integral method for the linearized rotating Navier-Stokes equations in exterior domain

In this paper, we apply the boundary integral method to the linearized rotating Navier-Stokes equations in exterior domain. Introducing some open ball which decomposes the exterior domain into a finite domain and an infinite domain, we obtain a coupled problem by the linearized rotating Navier-Stokes equations in finite domain and a boundary integral equation without using the artificial boundary condition. For the coupled problem, we show the existence and uniqueness of solution. Finally, we study the finite element approximation for the coupled problem and obtain the error estimate between the solution of the coupled problem and its approximation solution.     

Axisymmetric stokes exterior problem and its numerical computation

We establish variational formulation and prove the existence and uniqueness of the three dimensional axisymmetric Stokes exterior problem in weighted spaces.Error estimates and convergence for P2-P0 elements with infinite element methods are also obtained.Numerical experiments are presented to verify the theoretical analysis.     

Studies of an infinite element method for acoustical radiation

Infinite element computations are very efficient for predicting the vibro-acoustic response and sensitivities of a vibrating structure for an exterior acoustic domain. In addition, domain decomposition methods are very powerful algorithms for solving large linear systems in parallel. In this paper, an infinite element method is proposed and analyzed for parallel computations purpose. An original formulation of this method with Lagrange multipliers defined on (semi-)infinite space is presented. The implementation aspects of this method in an industrial acoustic software (SYSNOISE) are discussed. New numerical results illustrate the efficiency of the proposed method for realistic acoustical radiation problems.