I develop and mantain Algoim, a collection of high-order accurate numerical methods and C++ algorithms for working with implicitly defined geometry and level set methods, available on GitHub.

Motivated by multi-phase multi-physics applications, particularly those with evolving dynamic interfaces, the algorithms in Algoim target core, fundamental techniques in level set methods. They have been designed with a view to standard finite difference implementations as well as more advanced finite element and discontinuous Galerkin implementations, multi-threading and massively parallel MPI computation. The collection includes:

Algoim-based algorithms have been used across much of my research in high-order accurate computational physics, including complex flow in non-trivial geometry, free surface flow driven by intricate surface tension dynamics, multi-scale models of thin-film foam dynamics, multi-phase fluid flow, and petascale simulation of rotary bell atomisation dynamics. These algorithms have also found use in a number of other research projects around the world, including the DOE AMReX exacale initiative as well as work on extended finite element methods and cut cell finite volume methods in computational physics, chemistry, and materials.


Visit the Algoim page on GitHub.

Comments or suggestions

Feel free to contact me if you have any comments or suggestions: rsaye {at} lbl {dot} gov