Spectral discretization of the Stokes problem with mixed boundary conditions

The variational formulation of the Stokes problem with three independent unknowns, the vorticity, the velocity and the pressure, was born to handle non standard boundary conditions which involve the normal component of the velocity and the tangential components of the vorticity. We propose an extension of this formulation to the case of mixed boundary conditions in a three-dimensional domain. Next we consider a spectral discretization of this problem. A detailed numerical analysis leads to error estimates for the three unknowns and numerical experiments conrm the interest of the discretization

Analysis of Compatible Discrete Operator Schemes for the Stokes Equations on Polyhedral Meshes

Compatible Discrete Operator schemes preserve basic properties of the continuous model at the discrete level. They combine discrete differential operators that discretize exactly topological laws and discrete Hodge operators that approximate constitutive relations. We devise and analyze two families of such schemes for the Stokes equations in curl formulation, with the pressure degrees of freedom located at either mesh vertices or cells. The schemes ensure local mass and momentum conservation. We prove discrete stability by establishing novel discrete Poincar\'e inequalities. Using commutators related to the consistency error, we derive error estimates with first-order convergence rates for smooth solutions. We analyze two strategies for discretizing the external load, so as to deliver tight error estimates when the external load has a large irrotational or divergence-free part. Finally, numerical results are presented on three-dimensional polyhedral meshes

Discrete conservation properties for shallow water flows using mixed mimetic spectral elements

A mixed mimetic spectral element method is applied to solve the rotating shallow water equations. The mixed method uses the recently developed spectral element histopolation functions, which exactly satisfy the fundamental theorem of calculus with respect to the standard Lagrange basis functions in one dimension. These are used to construct tensor product solution spaces which satisfy the generalized Stokes theorem, as well as the annihilation of the gradient operator by the curl and the curl by the divergence. This allows for the exact conservation of first order moments (mass, vorticity), as well as quadratic moments (energy, potential enstrophy), subject to the truncation error of the time stepping scheme. The continuity equation is solved in the strong form, such that mass conservation holds point wise, while the momentum equation is solved in the weak form such that vorticity is globally conserved. While mass, vorticity and energy conservation hold for any quadrature rule, potential enstrophy conservation is dependent on exact spatial integration. The method possesses a weak form statement of geostrophic balance due to the compatible nature of the solution spaces and arbitrarily high order spatial error convergence

Incorporating variable viscosity in vorticity-based formulations for Brinkman equations

In this brief note, we introduce a non-symmetric mixed finite element formulation for Brinkman equations written in terms of velocity, vorticity and pressure with non-constant viscosity. The analysis is performed by the classical Babu\v{s}ka-Brezzi theory, and we state that any inf-sup stable finite element pair for Stokes approximating velocity and pressure can be coupled with a generic discrete space of arbitrary order for the vorticity. We establish optimal a priori error estimates which are further confirmed through computational example

Spectral method for the unsteady incompressible Navier-Stokes equations in gauge formulation

A spectral method which uses a gauge method, as opposed to a projection method, to decouple the computation of velocity and pressure in the unsteady incompressible Navier-Stokes equations, is presented. Gauge methods decompose velocity into the sum of an auxilary field and the gradient of a gauge variable, which may, in principle, be assigned arbitrary boundary conditions, thus overcoming the issue of artificial pressure boundary conditions in projection methods. A lid-driven cavity flow is used as a test problem. A subtraction method is used to reduce the pollution effect of singularities at the top corners of the cavity. A Chebyshev spectral collocation method is used to discretize spatially. An exponential time differencing method is used to discretize temporally. Matrix diagonalization procedures are used to compute solutions directly and efficiently. Numerical results for the flow at Reynolds number Re = 1000 are presented, and compared to benchmark results. It is shown that the method, called the spectral gauge method, is straightforward to implement, and yields accurate solutions if Neumann boundary conditions are imposed on the gauge variable, but suffers from reduced convergence rates if Dirichlet boundary conditions are imposed on the gauge variable