A C 0 conforming dg finite element method for biharmonic equations on triangle/tetrahedron

A C 0 conforming discontinuous Galerkin (CDG) finite element method is introduced for solving the biharmonic equation. The first strong gradient of C 0 finite element functions is a vector of discontinuous piecewise polynomials. The second gradient is the weak gradient of discontinuous piecewise polynomials. This method, by its name, uses nonconforming (non C 1) approximations and keeps simple formulation of conforming finite element methods without any stabilizers. Optimal order error estimates in both a discrete H 2 norm and the L 2 norm are established for the corresponding finite element solutions. Numerical results are presented to confirm the theory of convergence.

On Rational Krylov and Reduced Basis Methods for Fractional Diffusion

We establish an equivalence between two classes of methods for solving fractional diffusion problems, namely, Reduced Basis Methods (RBM) and Rational Krylov Methods (RKM). In particular, we demonstrate that several recently proposed RBMs for fractional diffusion can be interpreted as RKMs. This changed point of view allows us to give convergence proofs for some methods where none were previously available. We also propose a new RKM for fractional diffusion problems with poles chosen using the best rational approximation of the function 𝑧 −𝑠 with 𝑧 ranging over the spectral interval of the spatial discretization matrix. We prove convergence rates for this method and demonstrate numerically that it is competitive with or superior to many methods from the reduced basis, rational Krylov, and direct rational approximation classes. We provide numerical tests for some elliptic fractional diffusion model problems.

POD-Galerkin Model Order Reduction for Parametrized Nonlinear Time Dependent Optimal Flow Control: an Application to Shallow Water Equations

In this work we propose reduced order methods as a reliable strategy to efficiently solve parametrized optimal control problems governed by Shallow Waters Equations in a solution tracking setting. The physical parametrized model we deal with is nonlinear and time dependent: this leads to very time consuming simulations which can be unbearable e.g. in a marine environmental monitoring plan application. Our aim is to show how reduced order modelling could help in studying different configurations and phenomena in a fast way. After building the optimality system, we rely on a POD-Galerkin reduction in order to solve the optimal control problem in a low dimensional reduced space. The presented theoretical framework is actually suited to general nonlinear time dependent optimal control problems. The proposed methodology is finally tested with a numerical experiment: the reduced optimal control problem governed by Shallow Waters Equations reproduces the desired velocity and height profiles faster than the standard model, still remaining accurate.

A posteriori error analysis of an enriched Galerkin method of order one for the Stokes problem

We derive optimal reliability and efficiency of a posteriori error estimates for the steady Stokes problem, with a nonhomogeneous Dirichlet boundary condition, solved by a stable enriched Galerkin scheme (EG) of order one on triangular or quadrilateral meshes in ℝ2, and tetrahedral or hexahedral meshes in ℝ3.

Guaranteed Upper Bounds For The Velocity Error Of Pressure-Robust Stokes Discretisations

This paper aims to improve guaranteed error control for the Stokes problem with a focus on pressure-robustness, i.e. for discretisations that compute a discrete velocity that is independent of the exact pressure. A Prager-Synge type result relates the velocity errors of divergence-free primal and perfectly equilibrated dual mixed methods for the velocity stress. The first main result of the paper is a framework with relaxed constraints on the primal and dual method. This enables to use a recently developed mass conserving mixed stress discretisation for the design of equilibrated fluxes and to obtain pressure-independent guaranteed upper bounds for any pressure-robust (not necessarily divergence-free) primal discretisation. The second main result is a provably efficient local design of the equilibrated fluxes with comparably low numerical costs. Numerical examples verify the theoretical findings and show that efficiency indices of our novel guaranteed upper bounds are close to one.

An energy, momentum and angular momentum conserving scheme for a regularization model of incompressible flow

We introduce a new regularization model for incompressible fluid flow, which is a regularization of the EMAC (energy, momentum, and angular momentum conserving) formulation of the Navier-Stokes equations (NSE) that we call EMAC-Reg. The EMAC formulation has proved to be a useful formulation because it conserves energy, momentum and angular momentum even when the divergence constraint is only weakly enforced. However, it is still a NSE formulation and so cannot resolve higher Reynolds number flows without very fine meshes. By carefully introducing regularization into the EMAC formulation, we create a model more suitable for coarser mesh computations but that still conserves the same quantities as EMAC, i.e., energy, momentum, and angular momentum. We show that EMAC-Reg, when semi-discretized with a finite element spatial discretization is well-posed and optimally accurate. Numerical results are provided that show EMAC-Reg is a robust coarse mesh model.

Adaptive space-time finite element methods for parabolic optimal control problems

We present, analyze, and test locally stabilized space-time finite element methods on fully unstructured simplicial space-time meshes for the numerical solution of space-time tracking parabolic optimal control problems with the standard L 2-regularization. We derive a priori discretization error estimates in terms of the local mesh-sizes for shape-regular meshes. The adaptive version is driven by local residual error indicators, or, alternatively, by local error indicators derived from a new functional a posteriori error estimator. The latter provides a guaranteed upper bound of the error, but is more costly than the residual error indicators. We perform numerical tests for benchmark examples having different features. In particular, we consider a discontinuous target in form of a first expanding and then contracting ball in 3d that is fixed in the 4d space-time cylinder.

Error analysis for a vorticity/Bernoulli pressure formulation for the Oseen equations

A variational formulation is analysed for the Oseen equations written in terms of vorticity and Bernoulli pressure. The velocity is fully decoupled using the momentum balance equation, and it is later recovered by a post-process. A finite element method is also proposed, consisting in equal-order Nédélec finite elements and piecewise continuous polynomials for the vorticity and the Bernoulli pressure, respectively. The a priori error analysis is carried out in the L2-norm for vorticity, pressure, and velocity; under a smallness assumption either on the convecting velocity, or on the mesh parameter. Furthermore, an a posteriori error estimator is designed and its robustness and efficiency are studied using weighted norms. Finally, a set of numerical examples in 2D and 3D is given, where the error indicator serves to guide adaptive mesh refinement. These tests illustrate the behaviour of the new formulation in typical flow conditions, and also confirm the theoretical findings.