Numerical approximation of the generalized regularized long wave equation using Petrov–Galerkin finite element method

The generalized regularized long wave (GRLW) equation has been developed to model a variety of physical phenomena such as ion‐acoustic and magnetohydrodynamic waves in plasma, nonlinear transverse waves in shallow water and phonon packets in nonlinear crystals. This paper aims to develop and analyze a powerful numerical scheme for the nonlinear GRLW equation by Petrov–Galerkin method in which the element shape functions are cubic and weight functions are quadratic B‐splines. The proposed method is implemented to three reference problems involving propagation of the single solitary wave, interaction of two solitary waves and evolution of solitons with the Maxwellian initial condition. The variational formulation and semi‐discrete Galerkin scheme of the equation are firstly constituted. We estimate rate of convergence of such an approximation. Using Fourier stability analysis of the linearized scheme we show that the scheme is unconditionally stable. To verify practicality and robustness of the new scheme error norms L₂, L_∞ and three invariants I₁, I₂, and I₃ are calculated. The computed numerical results are compared with other published results and confirmed to be precise and effective.     

Application of the Sinc method to a dynamic elasto-plastic problem

This paper presents the application of Sinc bases to simulate numerically the dynamic behavior of a one-dimensional elastoplastic problem. The numerical methods that are traditionally employed to solve elastoplastic problems include finite difference, finite element and spectral methods. However, more recently, biorthogonal wavelet bases have been used to study the dynamic response of a uniaxial elasto-plastic rod [Giovanni F. Naldi, Karsten Urban, Paolo Venini, A wavelet-Galerkin method for elastoplasticity problems, Report 181, RWTH Aachen IGPM, and Math. Modelling and Scient. Computing, vol. 10, 2000]. In this paper the Sinc–Galerkin method is used to solve the straight elasto-plastic rod problem. Due to their exponential convergence rates and their need for a relatively fewer nodal points, Sinc based methods can significantly outperform traditional numerical methods [J. Lund, K.L. Bowers, Sinc Methods for Quadrature and Differential Equations, SIAM, Philadelphia, 1992]. However, the potential of Sinc-based methods for solving elastoplasticity problems has not yet been explored. The aim of this paper is to demonstrate the possible application of Sinc methods through the numerical investigation of the unsteady one dimensional elastic-plastic rod problem.     

A well‐balanced explicit/semi‐implicit finite element scheme for shallow water equations in drying–wetting areas

Numerical solutions of the shallow water equations can be used to reproduce flow hydrodynamics occurring in a wide range of regions. In hydraulic engineering, the objectives include the prediction of dam break wave propagation, fluvial floods and other catastrophic flooding phenomena, the modeling of estuarine and coastal circulations, and the design and optimization of hydraulic structures. In this paper, a well‐balanced explicit and semi‐implicit finite element scheme for shallow water equations over complex domains involving wetting and drying is proposed. The governing equations are discretized by a fractional finite element method using a two‐step Taylor–Galerkin scheme. First, the intermediate increment of conserved variable is obtained explicitly neglecting the pressure gradient term. This is then corrected for the effects of pressure once the pressure increment has been obtained from the Poisson equation. In order to maintain the ‘well‐balanced’ property, the pressure gradient term and bed slope terms are incorporated into the Poisson equation. Moreover, a local bed slope modification technique is employed in drying–wetting interface treatments. The proposed model is well validated against several theoretical benchmark tests. Copyright © 2014 John Wiley & Sons, Ltd.     

A greedy meshless local Petrov–Galerkin methodbased on radial basis functions

The meshless local Petrov–Galerkin (MLPG) method with global radial basis functions (RBF) as trial approximation leads to a full final linear system and a large condition number. This makes MLPG less efficient when the number of data points is increased. We can overcome this drawback if we avoid using more points from the data site than absolutely necessary. In this article, we equip the MLPG method with the greedy sparse approximation technique of (Schaback, Numercail Algorithms 67 (2014), 531–547) and use it for numerical solution of partial differential equations. This scheme uses as few neighbor nodal values as possible and allows to control the consistency error by explicit calculation. Whatever the given RBF is, the final system is sparse and the algorithm is well‐conditioned. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 847–861, 2016     

On Galerkin models for transport in groundwater

The Galerkin finite element model (GFEM) may provide oscillatory results when employed to predict contaminant transport in groundwater unless a very fine mesh is used. Adaptation of a very fine mesh may make the application of the GFEM impractical to field problems. The Petrov—Galerkin finite element models (PGFEMs) can provide oscillation free results for relatively coarser mesh. However, the PGFEM violates the Galerkin principle and introduces large “numerical” dispersion. The objective of this paper has been to develop accurate criteria to improve the applicability of the GFEM to obtain oscillation free accurate results for coarser mesh and compare its performance with that of the PGFEM. It has been shown that the GFEM provides oscillation free accurate results for coarser mesh with Peclet number Pe 20. Further, the GFEM prediction has always been more accurate than the PGFEM for a variety of source configurations and flow fields.     

Well‐posedness of the equations governing the motions of a one‐dimensional hybrid thermo‐elastic structure

In this paper a model for the vibrations of a one‐dimensional hybrid thermo‐elastic structure consisting of an extensible thermo‐elastic beam which is hinged at one end, with a rigid body attached to its free end, is studied with a view to establishing the existence of a unique solution in a weak sense. The model takes account of the effect of stretching on bending and rotational inertia. By treating eigenvalue problems with the spectral parameter also in the boundary conditions, we are able to employ the method of Faedo–Galerkin approximations. Copyright © 2004 John Wiley & Sons, Ltd.     

Space–time integrated least squares: a time‐marching approach

A time‐marching formulation is derived from the space–time integrated least squares (STILS) method for solving a pure hyperbolic convection equation and is numerically compared to various known methods. Copyright © 2004 John Wiley & Sons, Ltd.     

An invariant manifold approach to nonlinear normal modes of oscillation

Summary A method for determining the amplitude-dependent mode shapes and the corresponding modal dynamics of weakly nonlinear vibratory systems is described. The method is a combination of a Galerkin projection and invariant manifold techniques and is applied to a class of distributed parameter vibratory systems. In this paper the general theory for a class of conservative systems is outlined and applied to determine the nonlinear mode shapes and modal dynamics of a linear Euler-Bernoulli team attached to a nonlinear elastic foundation.     

Galerkin eigenvector approximations

How close are Galerkin eigenvectors to the best approximation available out of the trial subspace? Under a variety of conditions the Galerkin method gives an approximate eigenvector that approaches asymptotically the projection of the exact eigenvector onto the trial subspace--and this occurs more rapidly than the underlying rate of convergence of the approximate eigenvectors. Both orthogonal-Galerkin and Petrov-Galerkin methods are considered here with a special emphasis on nonselfadjoint problems, thus extending earlier studies by Chatelin, Babuska and Osborn, and Knyazev. Consequences for the numerical treatment of elliptic PDEs discretized either with finite element methods or with spectral methods are discussed. New lower bounds to the of a pair of operators are developed as well.