The present paper deals with the numerical analysis to the initial-boundary-value problem (ibvp) of pseudo-parabolic type in
where are continuously differentiable functions in and bounded above and below by positive constants. The right hand side of (1.1) involves linear and nonlinear terms ; they are assumed to be continuously differentiable functions of and .
In , we perform an extensive study, by computational means, of the use of spectral discretizations, of Galerkin and collocation type, based on Legendre and Chebyshev polynomials for models of type (1.1)-(1.3). The resulting semidiscrete systems are fully discretized there by suitable time integrators, with the aim at overcoming the possible midly stiff character of the ordinary differential problems and the difficulties to simulate nonsmooth data. This computational study is complemented by the present paper with a numerical analysis of the spectral discretization. More specifically, for Galerkin and collocation methods based on a family of Jacobi polynomials (which includes the Legendre and Chebyshev cases described in ) existence of solution of the semidiscrete systems and error estimates in suitable Sobolev norms are derived. As usual for this kind of approaches, the results proved here establish the rate of convergence of the spectral approximation in terms of the regularity of the data of the problem. In particular, they justify those experiments in  concerning spectral convergence in the smooth case.
Pseudo-parabolic equations of the form (1.1), in one or more dimensions, are used for modelling in different areas of Physics and Engineering. Relevant examples are the BBM-Burgers equation, a dissipative modification of the BBM equation for water waves, , and the pseudo-parabolic Buckley-Leverett equation describing two-phase flow in porous media, . We refer to the rich bibliography on it commented in .
The mathematical theory of pseudo-parabolic equations can be covered by [40, 41, 52, 53, 27, 13, 17, 20, 50]. Existence and uniqueness of weak solutions to nonlinear pseudo-parabolic equations are proved in , whereas the existence of weak solutions for degenerate cases is studied in [44, 43]. A homogenization of a closely related pseudo-parabolic system is considered in . Traveling wave solutions and their relation to non-standard shock solutions to hyperbolic conservation laws are investigated in [21, 23] for linear higher order terms. Uniqueness of weak solutions for a pseudo-parabolic equations modelling flow in porous media can be found in [18, 36, 17]. In , the authors study existence and uniqueness of weak solutions of the initial and boundary value problem for a fourth-order pseudo-parabolic equation with variable exponents of non-linearity, along with a long-time behaviour of weak solutions. Finally, existence of weak solutions for a nonlocal pseudo-parabolic model for Brinkman two-phase flow model in porous media has been recently established, .
Concerning the numerical approximation of equations of the form (1.1), the literature contains many references involving finite differences, [4, 5, 54, 6, 22, 28], as well as finite elements and finite volumes, [38, 45, 2, 56]
. We also mention some convergence results. First, stability and convergence of difference approximations to pseudo-parabolic partial differential equations is discussed in[29, 30] and the time stepping Crank-Nicolson Galerkin method to approximate several nonlinear Sobolev-type problems is analyzed in [26, 25]. Of particular relevance for the present study is the finite element approach for the nonlinear periodic-initial-boundary-value problem, , where Arnold and collaborators obtain optimal error estimates, in and norm, of a standard Galerkin method with continuous piecewise polynomials, and a nodal superconvergence. Moreover, Fourier spectral methods of Galerkin and collocation type for quasilinear pseudo-parabolic equations are analysed in . This is, to our knowledge, the main reference about the use of spectral methods to approximate Sobolev equations. More recent convergence results can be found in , where an analysis of a linearization scheme for an interior penalty discontinuous Galerkin for a pseudo-parabolic model in porous media applications is considered. High-order finite differences are employed in 
and B-spline quasi-interpolation methods in. In addition, an adaptive mesh approach for pseudo-parabolic-type problems is introduced in  and a Meshless RBFs method is considered in 
. Finally, unconditionally stable vector splitting schemes for pseudo-parabolic equations are constructed and analyzed in. It is worthwhile to mention that standard operator splitting may fail to capture the correct behavior of the solutions for pseudo-parabolic type differential models. In , the authors presented a non-splitting numerical method which is based on a fully coupled space-time mixed hybrid finite element/volume discretization approach to account for the delicate nonlinear balance between the hyperbolic flux and the pseudo-parabolic term linked to the full pseudo-parabolic differential model.
The structure of the paper is as follows. Section 2 is devoted to some theoretical aspects of (1.1)-(1.3) as the weak formulation (already mentioned in  for Legendre and Chebyshev cases) and assumptions on well-posedness. These preliminaries also include a summary on inverse inequalities, as well as projection and interpolation error estimates for the family of Jacobi polynomials under consideration. The contents of Section 2 will be used to the numerical analysis of the spectral Galerkin approximation in Section 3, and the collocation approximation in Section 4. Both contain, under suitable hypotheses on the data of the problem, results on the existence of numerical solution and convergence to the solution of (1.1)-(1.3). Concluding remarks and perspectives for future work are outlined in Section 5.
We now describe the main notation used throughout the paper. For positive integer , denotes the normed space of -functions on with as associated norm, while for nonnegative integer , is the space of -th order continuously differentiable functions on .
Let and define the Jacobi weight function
( corresponds to the Legendre case and to the Chebyshev case.) Then will denote the space of squared integrable functions with respect to the weighted inner product
and associated norm . For the Sobolev spaces integer (where ) the corresponding norm will be denoted by
We will also consider the spaces of functions such that . For , (and ) are defined by interpolation theory, . Note that in the case of the Legendre approximation () the spaces are the standard Sobolev spaces .
For an integer , will stand for the space of polynomials of degree at most on and
If and , stands for the space of functions on with norm . For an integer , the space of -th order continuously differentiable functions , where or , will be denoted by . Additionally, if , will stand for the normed space of functions with associated norm
We also denote by the space of functions with finite norm
where stands for the essential spectrum. Furthermore, (resp. ) will stand for the space of continuously differentiable (resp. uniformly bounded, continuously differentiable) functions in .
The analysis of the collocation methods requires the introduction of discrete norms. Let be the nodes and weights of the Gauss-Lobatto quadrature related to , [42, 16, 12]. For continuous on , the discrete inner product based on the Gauss-Lobatto data is denoted by
with associated norm . We recall that, 
Let be an integer. Then there exist positive constants , independent of , such that for any
Finally will be used to denote a generic, positive constant, independent of and , but that may depend on (this will be denoted by ).
2.1. Weak formulation
The analysis of the spectral discretizations that will be made below requires some hypotheses, properties and technical results concerning (1.1)-(1.3) and the approximation in weighted norms. From now on we will fix and consider the weight (1.4). The first property to be mentioned is the weak formulation of (1.1)-(1.3), cf. 
Since is bounded above and below by positive constants, then is equivalent to
Let and assume that , . Given , then there is a unique solution of (2.1) with bounded by a constant depending only on and the data of the problem. Furthermore, if with integer, , then for all and
where is a constant depending only on and the data of the problem.
2.2. Projection and interpolation errors with Jacobi polynomials
Here we collect several results concerning projection and interpolation errors with respect to the weighted inner product (1.5) and that will be used below. We refer to, e. g., [32, 42, 39, 11, 16, 15, 12, 51] for details and additional properties.
The estimates in the weighted Sobolev spaces , concern the use of the family of Jacobi polynomials , which are orthogonal to each other in . Particular cases such as Legendre and Chebyshev families correspond to and , respectively. Most properties of this Jacobi family (a particular one of the more general Jacobi polynomials , orthogonal in with ) are extension of the corresponding properties of the Legendre family, [11, 12].
We start with projection errors. Let be an integer, and let be the orthogonal projection of with respect to the inner product (1.5), and be the orthogonal projection of with respect to the inner product in
A third projection operator used below concerns the bilinear form given by (2.2). If then the orthogonal projection of with respect to is defined as such that
For this projection, we have, 
for . Furthermore, a generalized estimate can be obtained as follows. If , let be a polynomial such that, 
Finally, an additional estimate comparing the continuous and discrete inner products will be necessary: if and , then
3. Spectral Galerkin approximation
Let be an integer, . The semidiscrete Galekin approximation is defined as the function satisfying, cf. 
) are ensured by standard theory of ordinary differential equations (ode) when (3.1) is considered as a finite system for the coefficients of in some basis of , by using the property of ellipticity of and continuity of . Concerning this last point, some estimates on will be required in order to prove a global existence result. This is discussed in the following remark.
for some constant . Then
As far as the error estimates are concerned, let be the projection defined in (2.8) and
The right hand side of (3.10) is written as
Since (3.2) implies that and therefore , then writing yields
Actually ( see ) and
Similarly, we write