The Voronoi implicit interface method
Research
My research interests are in mathematical modelling and scientific computation, particularly numerical methods for interfacial fluid dynamics, high-order accurate algorithms for implicitly defined geometry, and multi-scale and multi-physics simulation. Below are some examples of my work.
Multiscale modelling of foam dynamics
Interfacial gauge methods
High-order accurate closest point calculation for implicitly defined surfaces
Implicit mesh discontinuous Galerkin methods
High-order accurate quadrature for implicitly defined domainsHigh-order accurate quadrature for implicitly defined domains
An implicitly defined domain is either a volumetric region or codimension-one surface whose shape is characterised implicitly by an isosurface of a continuous scalar function. A variety of applications involving implicitly defined geometry require the evaluation of integrals over such domains, including level set methods for propagating interfaces in computational physics, embedded boundary methods for solving partial differential equations on curved domains, and in treating jump conditions and singular source terms in weak formulations. A specific example is that of implicit mesh discontinuous Galerkin methods, which has been developed to facilitate high-order accurate modelling of interfacial fluid dynamics.
In practice, to calculate integrals over implicitly defined domains, a quadrature scheme must be computed. The high-order accurate numerical quadrature algorithms I have developed are based on the idea of converting the implicitly defined geometry into the graph of an implicitly defined height function, and leads to a recursive algorithm on the number of spatial dimensions, requiring only one-dimensional root finding and simple one-dimensional Gaussian quadrature schemes. The resulting quadrature schemes have strictly positive quadrature weights and inherits the high-order accuracy of Gaussian quadrature, e.g., 20th order accuracy can be demonstrated. The resulting algorithms have been a key component in implicit mesh discontinuous Galerkin methods.
Some examples illustrating the quadrature schemes constructed by the algorithm (based on a one-dimensional Gaussian quadrature scheme of order 4) for surface integrals (left column) and volume integrals (middle/right columns) [1]. The weights are coloured according to a scale that is normalised for each particular case: pale indicates a low-valued weight and dark blue indicates a high-valued weight.