Source author record

Marc Gerritsma

Marc Gerritsma appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2021arXiv

A mass-, kinetic energy- and helicity-conserving mimetic dual-field discretization for three-dimensional incompressible Navier-Stokes equations, part I: Periodic domains

We introduce a mimetic dual-field discretization which conserves mass, kinetic energy and helicity for three-dimensional incompressible Navier-Stokes equations. The discretization makes use of a conservative dual-field mixed weak formulation where two evolution equations of velocity are employed and dual representations of the solution are sought for each variable. A temporal discretization, which staggers the evolution equations and handles the nonlinearity such that the resulting discrete algebraic systems are linear and decoupled, is constructed. The spatial discretization is mimetic in the sense that the finite dimensional function spaces form a discrete de Rham complex. Conservation of mass, kinetic energy and helicity in the absence of dissipative terms is proven at the discrete level. Proper dissipation rates of kinetic energy and helicity in the viscous case is also proven. Numerical tests supporting the method are provided.

preprint2016arXiv

A mass, energy, enstrophy and vorticity conserving (MEEVC) mimetic spectral element discretization for the 2D incompressible Navier-Stokes equations

In this work we present a mimetic spectral element discretization for the 2D incompressible Navier-Stokes equations that in the limit of vanishing dissipation exactly preserves mass, kinetic energy, enstrophy and total vorticity on unstructured grids. The essential ingredients to achieve this are: (i) a velocity-vorticity formulation in rotational form, (ii) a sequence of function spaces capable of exactly satisfying the divergence free nature of the velocity field, and (iii) a conserving time integrator. Proofs for the exact discrete conservation properties are presented together with numerical test cases on highly irregular grids.

preprint2015arXiv

Mimetic spectral element method for Hamiltonian systems

There is a growing interest in the conservation of invariants when numerically solving a system of ordinary differential equations. Methods that exactly preserve these quantities in time are known as geometric integrators. In this paper we apply the recently developed mimetic framework to the solution of a system of first order ordinary differential equations. Depending on the discrete Hodge-* employed, two classes of arbitrary order time integrators are derived. It is shown that the one based on a canonical Hodge-* results in a symplectic integrator, whereas the one based on a Galerkin Hodge-* results in an energy preserving integrator. A set of numerical tests confirms these theoretical results.

preprint2013arXiv

Higher-order compatible discretization on hexahedrals

We derive a compatible discretization method that relies heavily on the underlying geometric structure, and obeys the topological sequences and commuting properties that are constructed. As a sample problem we consider the vorticity-velocity-pressure formulation of the Stokes problem. We motivate the choice for a mixed variational formulation based on both geometric as well as physical arguments. Numerical tests confirm the theoretical results that we obtain a pointwise divergence-free solution for the Stokes problem and that the method obtains optimal convergence rates.

preprint2013arXiv

Mimetic Spectral Element advection

We present a discretization of the linear advection of differential forms on bounded domains. The framework previously established is extended to incorporate the Lie derivative, $\mathcal L$, by means of Cartan's homotopy formula. The method is based on a physics-compatible discretization with spectral accuracy . It will be shown that the derived scheme has spectral convergence with local mass conservation. Artificial dispersion depends on the order of time integration.

preprint2013arXiv

Mixed Mimetic Spectral Element method applied to Darcy's problem

We present a discretization for Darcy's problem using the recently developed Mimetic Spectral Element Method. The gist lies in the exact discrete representation of integral relations. In this paper, an anisotropic flow through a porous medium is considered and a discretization of a full permeability tensor is presented. The performance of the method is evaluated on standard test problems, converging at the same rate as the best possible approximation.

preprint2013arXiv

Physics-compatible discretization techniques on single and dual grids, with application to the Poisson equation of volume forms

This paper introduces the basic concepts for physics-compatible discretization techniques. The paper gives a clear distinction between vectors and forms. Based on the difference between forms and pseudo-forms and the $\star$-operator which switches between the two, a dual grid description and a single grid description are presented. The dual grid method resembles a staggered finite volume method, whereas the single grid approach shows a strong resemblance with a finite element method. Both approaches are compared for the Poisson equation for volume forms.

preprint2012arXiv

Mixed Mimetic Spectral Element Method for Stokes Flow: A Pointwise Divergence-Free Solution

In this paper we apply the recently developed mimetic discretization method to the mixed formulation of the Stokes problem in terms of vorticity, velocity and pressure. The mimetic discretization presented in this paper and in [50] is a higher-order method for curvilinear quadrilaterals and hexahedrals. Fundamental is the underlying structure of oriented geometric objects, the relation between these objects through the boundary operator and how this defines the exterior derivative, representing the grad, curl and div, through the generalized Stokes theorem. The mimetic method presented here uses the language of differential $k$-forms with $k$-cochains as their discrete counterpart, and the relations between them in terms of the mimetic operators: reduction, reconstruction and projection. The reconstruction consists of the recently developed mimetic spectral interpolation functions. The most important result of the mimetic framework is the commutation between differentiation at the continuous level with that on the finite dimensional and discrete level. As a result operators like gradient, curl and divergence are discretized exactly. For Stokes flow, this implies a pointwise divergence-free solution. This is confirmed using a set of test cases on both Cartesian and curvilinear meshes. It will be shown that the method converges optimally for all admissible boundary conditions.

preprint2012arXiv

On Core and Curvature Corrections used in Straight-Line Vortex Filament Methods

The convergence characteristics of two viscous core corrections as used in straight-line segmentation methods are rigorously analysed. These are \emph{curvature corrections} that account for the induced velocity contribution at a point on a vortex filament due to local curvature and \emph{core corrections} that remove unrealistically high induced velocities near vortex segments. Two alternative versions of the latter correction are studied: the original as introduced by Scully and a recently introduced improved correction by the authors. The problem is analysed using a viscous vortex ring. A high-order numerical filament method is presented that uses rational B-spline curves to model the vortex ring geometry exactly. First, the separate influence of the two corrections on induced velocity values are studied using the numerical filament method. Afterwards, the corrections are combined in four separate cases that are studied using the segmentation method. It is found that in order to get accurate results at coarse discretisations, curvature corrections cannot be neglected. When used together with the original core correction model, results start to diverge from the reference values for discretisations smaller than approximately ten degrees. Combining the curvature correction with the improved core correction extends the discretisation region for which accurate results are found down by approximately one order. Finally, it is shown that segmentation methods that use the original core correction model without curvature correction completely fail to converge to the reference values.

preprint2011arXiv

Mimetic framework on curvilinear quadrilaterals of arbitrary order

In this paper higher order mimetic discretizations are introduced which are firmly rooted in the geometry in which the variables are defined. The paper shows how basic constructs in differential geometry have a discrete counterpart in algebraic topology. Generic maps which switch between the continuous differential forms and discrete cochains will be discussed and finally a realization of these ideas in terms of mimetic spectral elements is presented, based on projections for which operations at the finite dimensional level commute with operations at the continuous level. The two types of orientation (inner- and outer-orientation) will be introduced at the continuous level, the discrete level and the preservation of orientation will be demonstrated for the new mimetic operators. The one-to-one correspondence between the continuous formulation and the discrete algebraic topological setting, provides a characterization of the oriented discrete boundary of the domain. The Hodge decomposition at the continuous, discrete and finite dimensional level will be presented. It appears to be a main ingredient of the structure in this framework.