Werner Bauer

Werner Bauer, Colin J Cotter

We describe an energy-enstrophy conserving discretisation for the rotating shallow water equations with slip boundary conditions. This relaxes the assumption of boundary-free domains (periodic solutions or the surface of a sphere, for example) in the energy-enstrophy conserving formulation of McRae and Cotter (2014). This discretisation requires extra prognostic vorticity variables on the boundary in addition to the prognostic velocity and layer depth variables. The energy-enstrophy conservation properties hold for any appropriate set of compatible finite element spaces defined on arbitrary meshes with arbitrary boundaries. We demonstrate the conservation properties of the scheme with numerical solutions on a rotating hemisphere.

Rüdiger Brecht, Werner Bauer, Alexander Bihlo et al.

We develop a variational integrator for the shallow-water equations on a rotating sphere. The variational integrator is built around a discretization of the continuous Euler-Poincaré reduction framework for Eulerian hydrodynamics. We describe the discretization of the continuous Euler-Poincaré equations on arbitrary simplicial meshes. Standard numerical tests are carried out to verify the accuracy and the excellent conservational properties of the discrete variational integrator.

Werner Bauer, François Gay-Balmaz

This paper presents a geometric variational discretization of compressible fluid dynamics. The numerical scheme is obtained by discretizing, in a structure preserving way, the Lie group formulation of fluid dynamics on diffeomorphism groups and the associated variational principles. Our framework applies to irregular mesh discretizations in 2D and 3D. It systematically extends work previously made for incompressible fluids to the compressible case. We consider in detail the numerical scheme on 2D irregular simplicial meshes and evaluate the scheme numerically for the rotating shallow water equations. In particular, we investigate whether the scheme conserves stationary solutions, represents well the nonlinear dynamics, and approximates well the frequency relations of the continuous equations, while preserving conservation laws such as mass and total energy.

Werner Bauer, François Gay-Balmaz

The anelastic and pseudo-incompressible equations are two well-known soundproof approximations of compressible flows useful for both theoretical and numerical analysis in meteorology, atmospheric science, and ocean studies. In this paper, we derive and test structure-preserving numerical schemes for these two systems. The derivations are based on a discrete version of the Euler-Poincaré variational method. This approach relies on a finite dimensional approximation of the (Lie) group of diffeomorphisms that preserve weighted-volume forms. These weights describe the background stratification of the fluid and correspond to the weighed velocity fields for anelastic and pseudo-incompressible approximations. In particular, we identify to these discrete Lie group configurations the associated Lie algebras such that elements of the latter correspond to weighted velocity fields that satisfy the divergence-free conditions for both systems. Defining discrete Lagrangians in terms of these Lie algebras, the discrete equations follow by means of variational principles. Descending from variational principles, the schemes exhibit further a discrete version of Kelvin circulation theorem, are applicable to irregular meshes, and show excellent long term energy behavior. We illustrate the properties of the schemes by performing preliminary test cases.

Robert Piel, Werner Bauer

This paper introduces discrete differential form spaces over two-dimensional manifold meshes that feature enhanced subdivision-induced inter-element regularity compared to conventional finite element (FE) spaces. This increase in smoothness is achieved by pulling back refined subdivision basis functions along a hierarchy of increasingly fine meshes that are generated by a subdivision algorithm. We introduce a framework that casts several known instances of $k$-form subdivision schemes in the language of FE and derive conditions under which the resulting subdivision-induced hierarchy of FE function spaces satisfies a discrete de Rham complex. The paper further illustrates the enforcing of zero boundary conditions by discarding basis functions close to the mesh boundary and shows that this does not compromise the de Rham complex. To analyse our novel subdivision $k$-form spaces we solve the Maxwell eigenvalue problem to confirm the absence of spurious modes and to study the accuracy of the computed eigenvalues. Recovering accurately the expected analytic eigenvalue spectrum shows that our novel subdivision $k$-form spaces indeed preserve the de Rham complex, since this test case is known to be challenging for methods not preserving this structure. Further, we numerically investigate the approximation errors of these subdivision spaces for given analytic functions. The presented study shows that our method can be employed in two ways. Upon a suitable choice of parameters, the subdivision $k$-form spaces are up to $1.5$ orders of magnitude more accurate in the $L^2$ norm than conventional lowest-order FE spaces with the same number of degrees of freedom. Alternatively, for a given target accuracy, the number of required degrees of freedom can be significantly reduced, resulting in a speed-up by a factor of up to 6 for the discussed test cases.

Werner Bauer, Jörn Behrens

We introduce a new finite element (FE) discretization framework applicable for covariant split equations. The introduction of additional differential forms (DF) that form pairs with the original ones permits the splitting of the equations into topological momentum and continuity equations and metric-dependent closure equations that apply the Hodge-star operator. Our discretization framework conserves this geometrical structure and provides for all DFs proper FE spaces such that the differential operators hold in strong form. We introduce lowest possible order discretizations of the split 1D wave equations, in which the discrete momentum and continuity equations follow by trivial projections onto piecewise constant FE spaces, omitting partial integrations. Approximating the Hodge-star by nontrivial Galerkin projections (GP), the two discrete metric equations follow by projections onto either the piecewise constant (GP0) or piecewise linear (GP1) space. Our framework gives us three schemes with significantly different behavior. The split scheme using twice GP1 is unstable and shares the dispersion relation with the P1-P1 FE scheme that approximates both variables by piecewise linear spaces (P1). The split schemes that apply a mixture of GP1 and GP0 share the dispersion relation with the stable P1-P0 FE scheme that applies piecewise linear and piecewise constant (P0) spaces. However, the split schemes exhibit second order convergence for both quantities of interest. For the split scheme applying twice GP0, we are not aware of a corresponding standard formulation to compare with. Though it does not provide a satisfactory approximation of the dispersion relation as short waves are propagated much too fast, the discovery of the new scheme illustrates the potential of our discretization framework as a toolbox to study and find FE schemes by new combinations of FE spaces.