50,406 research outputs found
A modified Galerkin/finite element method for the numerical solution of the Serre-Green-Naghdi system
A new modified Galerkin / Finite Element Method is proposed for the numerical
solution of the fully nonlinear shallow water wave equations. The new numerical
method allows the use of low-order Lagrange finite element spaces, despite the
fact that the system contains third order spatial partial derivatives for the
depth averaged velocity of the fluid. After studying the efficacy and the
conservation properties of the new numerical method, we proceed with the
validation of the new numerical model and boundary conditions by comparing the
numerical solutions with laboratory experiments and with available theoretical
asymptotic results
Dispersion and dissipation error in high-order Runge-Kutta discontinuous Galerkin discretisations of the Maxwell equations
Different time-stepping methods for a nodal high-order discontinuous Galerkin discretisation of the Maxwell equations are discussed. A comparison between the most popular choices of Runge-Kutta (RK) methods is made from the point of view of accuracy and computational work. By choosing the strong-stability-preserving Runge-Kutta (SSP-RK) time-integration method of order consistent with the polynomial order of the spatial discretisation, better accuracy can be attained compared with fixed-order schemes. Moreover, this comes without a significant increase in the computational work. A numerical Fourier analysis is performed for this Runge-Kutta discontinuous Galerkin (RKDG) discretisation to gain insight into the dispersion and dissipation properties of the fully discrete scheme. The analysis is carried out on both the one-dimensional and the two-dimensional fully discrete schemes and, in the latter case, on uniform as well as on non-uniform meshes. It also provides practical information on the convergence of the dissipation and dispersion error up to polynomial order 10 for the one-dimensional fully discrete scheme
Numerical wave propagation for the triangular - finite element pair
Inertia-gravity mode and Rossby mode dispersion properties are examined for
discretisations of the linearized rotating shallow-water equations using the
- finite element pair on arbitrary triangulations in planar
geometry. A discrete Helmholtz decomposition of the functions in the velocity
space based on potentials taken from the pressure space is used to provide a
complete description of the numerical wave propagation for the discretised
equations. In the -plane case, this decomposition is used to obtain
decoupled equations for the geostrophic modes, the inertia-gravity modes, and
the inertial oscillations. As has been noticed previously, the geostrophic
modes are steady. The Helmholtz decomposition is used to show that the
resulting inertia-gravity wave equation is third-order accurate in space. In
general the \pdgp finite element pair is second-order accurate, so this leads
to very accurate wave propagation. It is further shown that the only spurious
modes supported by this discretisation are spurious inertial oscillations which
have frequency , and which do not propagate. The Helmholtz decomposition
also allows a simple derivation of the quasi-geostrophic limit of the
discretised - equations in the -plane case, resulting in a
Rossby wave equation which is also third-order accurate.Comment: Revised version prior to final journal submissio
A Priori Error Estimates for Mixed Finite Element -Schemes for the Wave Equation
A family of implicit-in-time mixed finite element schemes is presented for
the numerical approximation of the acoustic wave equation. The mixed space
discretization is based on the displacement form of the wave equation and the
time-stepping method employs a three-level one-parameter scheme. A rigorous
stability analysis is presented based on energy estimation and sharp stability
results are obtained. A convergence analysis is carried out and optimal a
priori error estimates for both displacement and pressure are
derived
- …