-
A Test Problem for Flow Codes
Authors:
Henry von Wahl,
L. Ridgway Scott
Abstract:
We propose a test problem for Navier-Stokes solvers based on the flow around a cylinder. We choose a range of Reynolds numbers for which the flow is time-dependent but can be characterized as essentially two-dimensional. The test problem requires accurate resolution of chaotic dynamics over a long time interval. It also requires the use of a relatively large computational domain, part of which is…
▽ More
We propose a test problem for Navier-Stokes solvers based on the flow around a cylinder. We choose a range of Reynolds numbers for which the flow is time-dependent but can be characterized as essentially two-dimensional. The test problem requires accurate resolution of chaotic dynamics over a long time interval. It also requires the use of a relatively large computational domain, part of which is curved, and it requires evaluation of derivatives of the solution and pressure on the curved boundary. We review the performance of different finite element methods for the proposed range of Reynolds numbers. These tests indicate that some of the most established methods do not capture the correct behavior.
△ Less
Submitted 26 April, 2024; v1 submitted 25 April, 2024;
originally announced April 2024.
-
A conservative Eulerian finite element method for transport and diffusion in moving domains
Authors:
Maxim Olshanskii,
Henry von Wahl
Abstract:
The paper introduces a finite element method for an Eulerian formulation of partial differential equations governing the transport and diffusion of a scalar quantity in a time-dependent domain. The method follows the idea from Lehrenfeld & Olshanskii [ESAIM: M2AN, 53(2): 585-614, 2019] of a solution extension to realise the Eulearian time-stepping scheme. However, a reformulation of the partial di…
▽ More
The paper introduces a finite element method for an Eulerian formulation of partial differential equations governing the transport and diffusion of a scalar quantity in a time-dependent domain. The method follows the idea from Lehrenfeld & Olshanskii [ESAIM: M2AN, 53(2): 585-614, 2019] of a solution extension to realise the Eulearian time-stepping scheme. However, a reformulation of the partial differential equation is suggested to derive a scheme which conserves the quantity under consideration exactly on the discrete level. For the spatial discretisation, the paper considers an unfitted finite element method. Ghost-penalty stabilisation is used to release the discrete solution extension and gives a scheme robust against arbitrary intersections between the mesh and geometry interface. The stability is analysed for both first- and second-order backward differentiation formula versions of the scheme. Several numerical examples in two and three spatial dimensions are included to illustrate the potential of this method.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
A coupled high-accuracy phase-field fluid-structure interaction framework for Stokes fluid-filled fracture surrounded by an elastic medium
Authors:
Henry von Wahl,
Thomas Wick
Abstract:
In this work, we couple a high-accuracy phase-field fracture reconstruction approach iteratively to fluid-structure interaction. The key motivation is to utilize phase-field modelling to compute the fracture path. A mesh reconstruction allows a switch from interface-capturing to interface-tracking in which the coupling conditions can be realized in a highly accurate fashion. Consequently, inside t…
▽ More
In this work, we couple a high-accuracy phase-field fracture reconstruction approach iteratively to fluid-structure interaction. The key motivation is to utilize phase-field modelling to compute the fracture path. A mesh reconstruction allows a switch from interface-capturing to interface-tracking in which the coupling conditions can be realized in a highly accurate fashion. Consequently, inside the fracture, a Stokes flow can be modelled that is coupled to the surrounding elastic medium. A fully coupled approach is obtained by iterating between the phase-field and the fluid-structure interaction model. The resulting algorithm is demonstrated for several numerical examples of quasi-static brittle fractures. We consider both stationary and quasi-stationary problems. In the latter, the dynamics arise through an incrementally increasing given pressure.
△ Less
Submitted 12 December, 2023; v1 submitted 29 August, 2023;
originally announced August 2023.
-
A discontinuous Galerkin approach for atmospheric flows with implicit condensation
Authors:
Sabine Doppler,
Philip L. Lederer,
Joachim Schöberl,
Henry von Wahl
Abstract:
We present a discontinuous Galerkin method for moist atmospheric dynamics, with and without warm rain. By considering a combined density for water vapour and cloud water, we avoid the need to model and compute a source term for condensation. We recover the vapour and cloud densities by solving a pointwise non-linear problem each time step. Consequently, we enforce the requirement for the water vap…
▽ More
We present a discontinuous Galerkin method for moist atmospheric dynamics, with and without warm rain. By considering a combined density for water vapour and cloud water, we avoid the need to model and compute a source term for condensation. We recover the vapour and cloud densities by solving a pointwise non-linear problem each time step. Consequently, we enforce the requirement for the water vapour not to be supersaturated implicitly. Together with an explicit time-stepping scheme, the method is highly parallelisable and can utilise high-performance computing hardware. Furthermore, the discretisation works on structured and unstructured meshes in two and three spatial dimensions. We illustrate the performance of our approach using several test cases in two and three spatial dimensions. In the case of a smooth, exact solution, we illustrate the optimal higher-order convergence rates of the method.
△ Less
Submitted 12 December, 2023; v1 submitted 23 May, 2023;
originally announced May 2023.
-
Unfitted Trefftz discontinuous Galerkin methods for elliptic boundary value problems
Authors:
Fabian Heimann,
Christoph Lehrenfeld,
Paul Stocker,
Henry von Wahl
Abstract:
We propose a new geometrically unfitted finite element method based on discontinuous Trefftz ansatz spaces. Trefftz methods allow for a reduction in the number of degrees of freedom in discontinuous Galerkin methods, thereby, the costs for solving arising linear systems significantly. This work shows that they are also an excellent way to reduce the number of degrees of freedom in an unfitted sett…
▽ More
We propose a new geometrically unfitted finite element method based on discontinuous Trefftz ansatz spaces. Trefftz methods allow for a reduction in the number of degrees of freedom in discontinuous Galerkin methods, thereby, the costs for solving arising linear systems significantly. This work shows that they are also an excellent way to reduce the number of degrees of freedom in an unfitted setting. We present a unified analysis of a class of geometrically unfitted discontinuous Galerkin methods with different stabilisation mechanisms to deal with small cuts between the geometry and the mesh. We cover stability and derive a-priori error bounds, including errors arising from geometry approximation for the class of discretisations for a model Poisson problem in a unified manner. The analysis covers Trefftz and full polynomial ansatz spaces, alike. Numerical examples validate the theoretical findings and demonstrate the potential of the approach.
△ Less
Submitted 9 June, 2023; v1 submitted 23 December, 2022;
originally announced December 2022.
-
A high-accuracy framework for phase-field fracture interface reconstructions with application to Stokes fluid-filled fracture surrounded by an elastic medium
Authors:
Henry von Wahl,
Thomas Wick
Abstract:
This work considers a Stokes flow in a deformable fracture interacting with a linear elastic medium. To this end, we employ a phase-field model to approximate the crack dynamics. Phase-field methods belong to interface-capturing approaches in which the interface is only given by a smeared zone. For multi-domain problems, the accuracy of the coupling conditions is, however, of utmost importance. He…
▽ More
This work considers a Stokes flow in a deformable fracture interacting with a linear elastic medium. To this end, we employ a phase-field model to approximate the crack dynamics. Phase-field methods belong to interface-capturing approaches in which the interface is only given by a smeared zone. For multi-domain problems, the accuracy of the coupling conditions is, however, of utmost importance. Here, interface-tracking methods are preferred, since the interface is resolved on mesh edges up to discretization errors, but it does not depend on the length scale parameter of some smeared zone. The key objective of this work is to construct a robust framework that computes first a crack path via the phase-field method (interface-capturing) and then does an interface-tracking reconstruction. We then discuss several approaches to reconstruct the Eulerian description of the open crack domain. This includes unfitted approaches where a level-set of the crack interface is constructed and an approach where the geometry is re-meshed. Using this reconstructed domain, we can compute the fluid-structure interaction problem between the fluid in the crack and the interacting solid. With the explicit mesh reconstruction of the two domains, we can then use an interface-tracking Arbitrary-Lagrangian-Eulerian (ALE) discretisation approach for the resulting fluid-structure interaction (FSI) problem. Our algorithmic procedure is realised in one final numerical algorithm and one implementation. We substantiate our approach using several numerical examples based on Sneddon's benchmark and corresponding extensions to Stokes fluid-filled regimes.
△ Less
Submitted 19 May, 2023; v1 submitted 15 December, 2022;
originally announced December 2022.
-
Error analysis for a parabolic PDE model problem on a coupled moving domain in a fully Eulerian framework
Authors:
Henry von Wahl,
Thomas Richter
Abstract:
We introduce an unfitted finite element method with Lagrange-multipliers to study an Eulerian time stepping scheme for moving domain problems applied to a model problem where the domain motion is implicit to the problem. We consider a parabolic partial differential equation (PDE) in the bulk domain, and the domain motion is described by an ordinary differential equation (ODE), coupled to the bulk…
▽ More
We introduce an unfitted finite element method with Lagrange-multipliers to study an Eulerian time stepping scheme for moving domain problems applied to a model problem where the domain motion is implicit to the problem. We consider a parabolic partial differential equation (PDE) in the bulk domain, and the domain motion is described by an ordinary differential equation (ODE), coupled to the bulk partial differential equation through the transfer of forces at the moving interface. The discretisation is based on an unfitted finite element discretisation on a time-independent mesh. The method-of-lines time discretisation is enabled by an implicit extension of the bulk solution through additional stabilisation, as introduced by Lehrenfeld & Olshanskii (ESAIM: M2AN, 53:585-614, 2019). The analysis of the coupled problem relies on the Lagrange-multiplier formulation, the fact that the Lagrange-multiplier solution is equal to the normal stress at the interface and that the motion of the interface is given through rigid body motion. This paper covers the complete stability analysis of the method and an error estimate in the energy norm, under an assumption on the discrete interface velocity. This includes the dynamic error in the domain motion resulting from the discretised ODE and the forces from the discretised PDE. To the best of our knowledge this is the first error analysis of this type of coupled moving domain problem in a fully Eulerian framework. Numerical examples illustrate the theoretical results.
△ Less
Submitted 1 August, 2022; v1 submitted 10 November, 2021;
originally announced November 2021.
-
Using a deep neural network to predict the motion of under-resolved triangular rigid bodies in an incompressible flow
Authors:
Henry von Wahl,
Thomas Richter
Abstract:
We consider non-spherical rigid body particles in an incompressible fluid in the regime where the particles are too large to assume that they are simply transported with the fluid without back-coupling and where the particles are also too small to make fully resolved direct numerical simulations feasible. Unfitted finite element methods with ghost-penalty stabilisation are well suited to fluid-str…
▽ More
We consider non-spherical rigid body particles in an incompressible fluid in the regime where the particles are too large to assume that they are simply transported with the fluid without back-coupling and where the particles are also too small to make fully resolved direct numerical simulations feasible. Unfitted finite element methods with ghost-penalty stabilisation are well suited to fluid-structure-interaction problems as posed by this setting, due to the flexible and accurate geometry handling and for allowing topology changes in the geometry. In the computationally under resolved setting posed here, accurate computations of the forces by their boundary integral formulation are not viable. Furthermore, analytical laws are not available due to the shape of the particles. However, accurate values of the forces are essential for realistic motion of the particles. To obtain these forces accurately, we train an artificial deep neural network using data from prototypical resolved simulations. This network is then able to predict the force values based on information which can be obtained accurately in an under-resolved setting. As a result, we obtain forces on very coarse and under-resolved meshes which are on average an order of magnitude more accurate compared to the direct boundary-integral computation from the Navier-Stokes solution, leading to solid motion comparable to that obtained on highly resolved meshes that would substantially increase the simulation costs.
△ Less
Submitted 10 June, 2021; v1 submitted 23 February, 2021;
originally announced February 2021.
-
Falling balls in a viscous fluid with contact: Comparing numerical simulations with experimental data
Authors:
Henry von Wahl,
Thomas Richter,
Stefan Frei,
Thomas Hagemeier
Abstract:
We evaluate a number of different finite element approaches for fluid-structure (contact) interaction problems against data from physical experiments. For this we take the data from experiments by Hagemeier [Mendeley Data, doi: 10.17632/mf27c92nc3.1]. This consists of trajectories of single particles falling through a highly viscous fluid and rebounding off the bottom fluid tank wall. The resultin…
▽ More
We evaluate a number of different finite element approaches for fluid-structure (contact) interaction problems against data from physical experiments. For this we take the data from experiments by Hagemeier [Mendeley Data, doi: 10.17632/mf27c92nc3.1]. This consists of trajectories of single particles falling through a highly viscous fluid and rebounding off the bottom fluid tank wall. The resulting flow is in the transitional regime between creeping and turbulent flows. This type of configuration is particularly challenging for numerical methods due to the large change of the fluid domain and the contact between the wall and particle. In the numerical simulations we consider both rigid body and linear elasticity models for the falling particles. In the first case, we compare results obtained with the well established Arbitrary Lagrangian Eulerian (ALE) approach and a moving domain CutFEM method together with a simple and common approach for contact avoidance. For the full fluid-structure interaction (FSI) problem with contact, we use a fully Eulerian approach in combination with a unified FSI-contact treatment using Nitsche's method. For higher computational efficiency we use the geometrical symmetry of the experimental set up to reformulate the FSI system into two spatial dimensions. Finally, we show full three dimensional ALE computations to study the effects of small perturbations in the initial state of the particle to investigate deviations from a perfectly vertical fall observed in the experiment. The methods are implemented in open-source finite element libraries and the results are made freely available to aide reproducibility.
△ Less
Submitted 17 November, 2020;
originally announced November 2020.
-
An unfitted Eulerian finite element method for the time-dependent Stokes problem on moving domains
Authors:
Henry von Wahl,
Thomas Richter,
Christoph Lehrenfeld
Abstract:
We analyse a Eulerian Finite Element method, combining a Eulerian time-stepping scheme applied to the time-dependent Stokes equations using the CutFEM approach with inf-sup stable Taylor-Hood elements for the spatial discretisation. This is based on the method introduced by Lehrenfeld \& Olshanskii [ESAIM: M2AN 53(2):585--614] in the context of a scalar convection-diffusion problems on moving doma…
▽ More
We analyse a Eulerian Finite Element method, combining a Eulerian time-stepping scheme applied to the time-dependent Stokes equations using the CutFEM approach with inf-sup stable Taylor-Hood elements for the spatial discretisation. This is based on the method introduced by Lehrenfeld \& Olshanskii [ESAIM: M2AN 53(2):585--614] in the context of a scalar convection-diffusion problems on moving domains, and extended to the non-stationary Stokes problem on moving domains by Burman, Frei \& Massing [arXiv:1910.03054 [math.NA]] using stabilised equal-order elements. The analysis includes the geometrical error made by integrating over approximated levelset domains in the discrete CutFEM setting. The method is implemented and the theoretical results are illustrated using numerical examples.
△ Less
Submitted 6 February, 2020;
originally announced February 2020.
-
Numerical benchmarking of fluid-rigid body interactions
Authors:
Henry von Wahl,
Thomas Richter,
Christoph Lehrenfeld,
Jan Heiland,
Piotr Minakowski
Abstract:
We propose a fluid-rigid body interaction benchmark problem, consisting of a solid spherical obstacle in a Newtonian fluid, whose centre of mass is fixed but is free to rotate. A number of different problems are defined for both two and three spatial dimensions. The geometry is chosen specifically, such that the fluid-solid partition does not change over time and classical fluid solvers are able t…
▽ More
We propose a fluid-rigid body interaction benchmark problem, consisting of a solid spherical obstacle in a Newtonian fluid, whose centre of mass is fixed but is free to rotate. A number of different problems are defined for both two and three spatial dimensions. The geometry is chosen specifically, such that the fluid-solid partition does not change over time and classical fluid solvers are able to solve the fluid-structure interaction problem. We summarise the different approaches used to handle the fluid-solid coupling and numerical methods used to solve the arising problems. The results obtained by the described methods are presented and we give reference intervals for the relevant quantities of interest.
△ Less
Submitted 19 September, 2019; v1 submitted 13 August, 2019;
originally announced August 2019.