Tim van Beeck

Patrick E. Farrell, Tim van Beeck, Umberto Zerbinati

We analyse the nematic Helmholtz-Korteweg equation, a variant of the classical Helmholtz equation that describes time-harmonic wave propagation in calamitic fluids in the presence of nematic order. A prominent example is given by nematic liquid crystals, which can be modeled as nematic Korteweg fluids - that is, fluids whose stress tensor depends on density gradients and on a nematic director describing the orientation of the anisotropic molecules. These materials exhibit anisotropic acoustic properties that can be tuned by external electromagnetic fields, making them attractive for potential applications such as tunable acoustic resonators. We prove the existence and uniqueness of solutions to this equation in two and three dimensions for suitable (nonresonant) wave numbers and propose a convergent discretisation for its numerical solution. The discretisation of this problem is nontrivial as it demands high regularity and involves unfamiliar boundary conditions. We address these challenges by using high-order conforming finite elements and enforcing the boundary conditions with Nitsche's method. We illustrate our analysis with numerical simulations in two dimensions.

Tim Brüers, Christoph Lehrenfeld, Tim van Beeck et al.

We present a discrete Helmholtz--Hodge decomposition for H(div)-conforming Brezzi--Douglas--Marini (BDM) finite elements on triangulated surfaces of arbitrary topology. The divergence-free BDM subspace is split L2-orthogonally into rotated gradients of a continuous streamfunction space and a finite-dimensional space of discrete harmonic fields whose dimension equals the first Betti number of the surface. Consequently, any incompressible flow discretized on this subspace can be reformulated with a scalar streamfunction and finitely many harmonic coefficients as the only unknowns. This eliminates the pressure and the saddle-point structure while ensuring exact tangentiality, pointwise divergence-freeness, and pressure-robustness. We present a randomized algorithm for constructing the harmonic basis and discuss implementation aspects including hybridization, efficient treatment of the harmonic unknowns, and pressure reconstruction. Numerical experiments for unsteady surface Navier--Stokes equations on a trefoil knot and a multiply-connected sculpture surface demonstrate the method and illustrate the physical role of the harmonic velocity component.

Philip L. Lederer, Christoph Lehrenfeld, Christian Merdon et al.

This paper studies pressure-robustness for the axisymmetric Stokes problem. The transformation to cylindrical coordinates requires that the radially weighted velocity is divergence-free in the classical sense. Consequently, traditional divergence-free finite element methods from the Cartesian setting -- even if inf-sup stable -- are in general not divergence-free in the axisymmetric formulation. We therefore explore the approach that restores pressure-robustness via reconstruction operators for a low-order Bernardi--Raugel discretization. We show that an application of standard interpolation operators from the Cartesian setting to radially weighted test functions works in principle, but it lacks properties needed to derive optimal consistency error estimates. To address this, we introduce a reconstruction operator into a finite element space spanned by Raviart--Thomas functions that are modified such that they vanish on the rotation axis. This vanishing-on-axis property is the key to obtain optimal consistency error estimates. Numerical examples demonstrate the overall feasibility of the approach and include cases where the vanishing-on-axis property yields significantly better results.