-
A hypocoercivity-exploiting stabilised finite element method for Kolmogorov equation
Authors:
Zhaonan Dong,
Emmanuil H. Georgoulis,
Philip J. Herbert
Abstract:
We propose a new stabilised finite element method for the classical Kolmogorov equation. The latter serves as a basic model problem for large classes of kinetic-type equations and, crucially, is characterised by degenerate diffusion. The stabilisation is constructed so that the resulting method admits a \emph{numerical hypocoercivity} property, analogous to the corresponding property of the PDE pr…
▽ More
We propose a new stabilised finite element method for the classical Kolmogorov equation. The latter serves as a basic model problem for large classes of kinetic-type equations and, crucially, is characterised by degenerate diffusion. The stabilisation is constructed so that the resulting method admits a \emph{numerical hypocoercivity} property, analogous to the corresponding property of the PDE problem. More specifically, the stabilisation is constructed so that spectral gap is possible in the resulting ``stronger-than-energy'' stabilisation norm, despite the degenerate nature of the diffusion in Kolmogorov, thereby the method has a provably robust behaviour as the ``time'' variable goes to infinity. We consider both a spatially discrete version of the stabilised finite element method and a fully discrete version, with the time discretisation realised by discontinuous Galerkin timestepping. Both stability and a priori error bounds are proven in all cases. Numerical experiments verify the theoretical findings.
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Convergence of a steepest descent algorithm in shape optimisation using $W^{1,\infty}$ functions
Authors:
Klaus Deckelnick,
Philip J. Herbert,
Michael Hinze
Abstract:
Built upon previous work of the authors in (Deckelnick, Herbert, and Hinze, ESAIM: COCV 28 (2022)), we present a general shape optimisation framework based on the method of mappings in the $W^{1,\infty}$ topology together with a suitable finite element discretisation. For the numerical solution of the respective discrete shape optimisation problems we propose a steepest descent minimisation algori…
▽ More
Built upon previous work of the authors in (Deckelnick, Herbert, and Hinze, ESAIM: COCV 28 (2022)), we present a general shape optimisation framework based on the method of mappings in the $W^{1,\infty}$ topology together with a suitable finite element discretisation. For the numerical solution of the respective discrete shape optimisation problems we propose a steepest descent minimisation algorithm with Armijo-Goldstein stepsize rule. We show that the sequence generated by this descent method globally converges, and under mild assumptions also, that every accumulation point of this sequence is a stationary point of the shape functional. Moreover, for the mesh discretisation parameter tending to zero we under mild assumptions prove convergence of the discrete stationary shapes in the Hausdorff complementary metric. To illustrate our approach we present a selection of numerical examples for PDE constrained shape optimisation problems, where we include numerical convergence studies which support our analytical findings.
△ Less
Submitted 23 October, 2023;
originally announced October 2023.
-
Shape optimization in $W^{1,\infty}$ with geometric constraints: a study in distributed-memory systems
Authors:
Philip J. Herbert,
Jose A. Pinzon Escobar,
Martin Siebenborn
Abstract:
In this paper we present a shape optimization scheme which utilizes the alternating direction method of multipliers (ADMM) to approximate a direction of steepest descent in $W^{1,\infty}$. The followed strategy is a combination of the approaches presented in Deckelnick, Herbert, and Hinze, ESAIM: COCV 28 (2022) and Müller et al. SIAM SISC 45 (2023). This has appeared previously for relatively simp…
▽ More
In this paper we present a shape optimization scheme which utilizes the alternating direction method of multipliers (ADMM) to approximate a direction of steepest descent in $W^{1,\infty}$. The followed strategy is a combination of the approaches presented in Deckelnick, Herbert, and Hinze, ESAIM: COCV 28 (2022) and Müller et al. SIAM SISC 45 (2023). This has appeared previously for relatively simple elliptic PDEs with geometric constraints which were handled using an ad-hoc projection. Here, however, the optimization problem is expanded to include geometric constraints, which are systematically fulfilled. Moreover, this results in a nonlinear system of equations, which is challenging from a computational perspective. Simulations of a fluid dynamics case study are carried out to benchmark the novel method. Results are given to show that, compared to other methods, the proposed methodology allows for larger deformations without affecting the convergence of the used numerical methods. The mesh quality is studied across the surface of the optimized obstacle, and is further compared to previous approaches which used descents in $W^{p,\infty}$. The parallel scalability is tested on a distributed-memory system to illustrate the potential of the proposed techniques in a more complex, industrial setting.
△ Less
Submitted 27 September, 2023;
originally announced September 2023.
-
The Stochastic Steepest Descent Method for Robust Optimization in Banach Spaces
Authors:
Neil K. Chada,
Philip J. Herbert
Abstract:
Stochastic gradient methods have been a popular and powerful choice of optimization methods, aimed at minimizing functions. Their advantage lies in the fact that that one approximates the gradient as opposed to using the full Jacobian matrix. One research direction, related to this, has been on the application to infinite-dimensional problems, where one may naturally have a Hilbert space framework…
▽ More
Stochastic gradient methods have been a popular and powerful choice of optimization methods, aimed at minimizing functions. Their advantage lies in the fact that that one approximates the gradient as opposed to using the full Jacobian matrix. One research direction, related to this, has been on the application to infinite-dimensional problems, where one may naturally have a Hilbert space framework. However, there has been limited work done on considering this in a more general setup, such as where the natural framework is that of a Banach space. This article aims to address this by the introduction of a novel stochastic method, the stochastic steepest descent method (SSD). The SSD will follow the spirit of stochastic gradient descent, which utilizes Riesz representation to identify gradients and derivatives. Our choice for using such a method is that it naturally allows one to adopt a Banach space setting, for which recent applications have exploited the benefit of this, such as in PDE-constrained shape optimization. We provide a convergence theory related to this under mild assumptions. Furthermore, we demonstrate the performance of this method on a couple of numerical applications, namely a $p$-Laplacian and an optimal control problem. Our assumptions are verified in these applications.
△ Less
Submitted 11 August, 2023;
originally announced August 2023.
-
PDE constrained shape optimisation with first-order and Newton-type methods in the $W^{1,\infty}$ topology
Authors:
Klaus Deckelnick,
Philip J. Herbert,
Michael Hinze
Abstract:
We present a general shape optimisation framework based on the method of mappings in the $W^{1,\infty}$ topology. We propose steepest descent and Newton-like minimisation algorithms for the numerical solution of the respective shape optimisation problems. Our work is built upon previous work of the authors in (Deckelnick, Herbert, and Hinze, ESAIM: COCV 28 (2022)), where a $W^{1,\infty}$ framework…
▽ More
We present a general shape optimisation framework based on the method of mappings in the $W^{1,\infty}$ topology. We propose steepest descent and Newton-like minimisation algorithms for the numerical solution of the respective shape optimisation problems. Our work is built upon previous work of the authors in (Deckelnick, Herbert, and Hinze, ESAIM: COCV 28 (2022)), where a $W^{1,\infty}$ framework for star-shaped domains is proposed. To illustrate our approach we present a selection of PDE constrained shape optimisation problems and compare our findings to results from so far classical Hilbert space methods and recent $p$-approximations.
△ Less
Submitted 21 December, 2023; v1 submitted 20 January, 2023;
originally announced January 2023.
-
Shape Optimisation with $W^{1,\infty}$: A connection between the steepest descent and Optimal Transport
Authors:
Philip J. Herbert
Abstract:
In this work, we discuss the task of finding a direction of optimal descent for problems in Shape Optimisation and its relation to the dual problem in Optimal Transport. This link was first observed in a previous work which sought minimisers of a shape derivative over the space of Lipschitz functions which may be closely related to the $\infty$-Laplacian. We provide some results of shape optimisat…
▽ More
In this work, we discuss the task of finding a direction of optimal descent for problems in Shape Optimisation and its relation to the dual problem in Optimal Transport. This link was first observed in a previous work which sought minimisers of a shape derivative over the space of Lipschitz functions which may be closely related to the $\infty$-Laplacian. We provide some results of shape optimisation using this novel Lipschitz approach, highlighting the difference between the Lipschitz and $W^{1,\infty}$ semi-norms. After this, we provide an overview of the necessary results from Optimal transport in order to make a direct link to the Shape optimisation of star-shaped domains. Demonstrative numerical experiments are provided.
△ Less
Submitted 19 January, 2023;
originally announced January 2023.
-
An evolving space framework for Oseen equations on a moving domain
Authors:
Ana Djurdjevac,
Carsten Gräser,
Philip J. Herbert
Abstract:
This article considers non-stationary incompressible linear fluid equations in a moving domain. We demonstrate the existence and uniqueness of an appropriate weak formulation of the problem by making use of the theory of time-dependent Bochner spaces. It is not possible to directly apply established evolving Hilbert space theory due to the incompressibility constraint. After we have established th…
▽ More
This article considers non-stationary incompressible linear fluid equations in a moving domain. We demonstrate the existence and uniqueness of an appropriate weak formulation of the problem by making use of the theory of time-dependent Bochner spaces. It is not possible to directly apply established evolving Hilbert space theory due to the incompressibility constraint. After we have established the well-posedness, we derive and analyse a first order time discretisation of the system.
△ Less
Submitted 8 September, 2023; v1 submitted 15 June, 2021;
originally announced June 2021.
-
A novel $W^{1,\infty}$ approach to shape optimisation with Lipschitz domains
Authors:
Klaus Deckelnick,
Philip J. Herbert,
Michael Hinze
Abstract:
This article introduces a novel method for the implementation of shape optimisation with Lipschitz domains. We propose to use the shape derivative to determine deformation fields which represent steepest descent directions of the shape functional in the $W^{1,\infty}-$ topology. The idea of our approach is demonstrated for shape optimisation of $n$-dimensional star-shaped domains, which we represe…
▽ More
This article introduces a novel method for the implementation of shape optimisation with Lipschitz domains. We propose to use the shape derivative to determine deformation fields which represent steepest descent directions of the shape functional in the $W^{1,\infty}-$ topology. The idea of our approach is demonstrated for shape optimisation of $n$-dimensional star-shaped domains, which we represent as functions defined on the unit $(n-1)$-sphere. In this setting we provide the specific form of the shape derivative and prove the existence of solutions to the underlying shape optimisation problem. Moreover, we show the existence of a direction of steepest descent in the $W^{1,\infty}-$ topology. We also note that shape optimisation in this context is closely related to the $\infty-$Laplacian, and to optimal transport, where we highlight the latter in the numerics section. We present several numerical experiments in two dimensions illustrating that our approach seems to be superior over a widely used Hilbert space method in the considered examples, in particular in developing optimised shapes with corners.
△ Less
Submitted 10 December, 2021; v1 submitted 25 March, 2021;
originally announced March 2021.
-
A formula for membrane mediated point particle interactions on near spherical biomembranes
Authors:
Charles M. Elliott,
Philip J. Herbert
Abstract:
We consider a model of a biomembrane with attached proteins. The membrane is represented by a near spherical continuous surface and attached proteins are described as discrete rigid structures which attach to the membrane at a finite number of points. The resulting surface minimises a quadratic elastic energy (obtained by a perturbation of the Canham-Helfrich energy) subject to the point constrain…
▽ More
We consider a model of a biomembrane with attached proteins. The membrane is represented by a near spherical continuous surface and attached proteins are described as discrete rigid structures which attach to the membrane at a finite number of points. The resulting surface minimises a quadratic elastic energy (obtained by a perturbation of the Canham-Helfrich energy) subject to the point constraints which are imposed by the attachment of the proteins. We calculate the derivative of the energy with respect to protein configurations. The proteins are constrained to move tangentially by translation and by rotation in the axis normal to a reference point. Previous studies have typically restricted themselves to a nearly flat membrane and circular inclusions. A numerically accessible representation of this derivative is derived and employed in some numerical experiments.
△ Less
Submitted 8 December, 2020;
originally announced December 2020.
-
Small deformations of spherical biomembranes
Authors:
Charles M. Elliott,
Luke Hatcher,
Philip J. Herbert
Abstract:
In this contribution to the proceedings of the 11th Mathematical Society of Japan (MSJ) Seasonal Institute (July 2018) we give an overview of some recent work on a mathematical model for small deformations of a spherical membrane. The idea is to consider perturbations to minimisers of a surface geometric energy. The model is obtained from consideration of second order approximations to a perturbed…
▽ More
In this contribution to the proceedings of the 11th Mathematical Society of Japan (MSJ) Seasonal Institute (July 2018) we give an overview of some recent work on a mathematical model for small deformations of a spherical membrane. The idea is to consider perturbations to minimisers of a surface geometric energy. The model is obtained from consideration of second order approximations to a perturbed energy. In particular, the considered problems involve particle constraints and surface phase field energies.
△ Less
Submitted 7 November, 2019;
originally announced November 2019.
-
Second order splitting of a class of fourth order PDEs with point constraints
Authors:
Charles M. Elliott,
Philip J. Herbert
Abstract:
We formulate a well-posedness and approximation theory for a class of generalised saddle point problems with a specific form of constraints. In this way we develop an approach to a class of fourth order elliptic partial differential equations with point constraints using the idea of splitting into coupled second order equations. An approach is formulated using a penalty method to impose the constr…
▽ More
We formulate a well-posedness and approximation theory for a class of generalised saddle point problems with a specific form of constraints. In this way we develop an approach to a class of fourth order elliptic partial differential equations with point constraints using the idea of splitting into coupled second order equations. An approach is formulated using a penalty method to impose the constraints. Our main motivation is to treat certain fourth order equations involving the biharmonic operator and point Dirichlet constraints for example arising in the modelling of biomembranes on curved and flat surfaces but the approach may be applied more generally. The theory for well-posedness and approximation is presented in an abstract setting. Several examples are described together with some numerical experiments.
△ Less
Submitted 28 April, 2020; v1 submitted 7 November, 2019;
originally announced November 2019.