helmholtz-x is a python library built upon open-source finite element framework, DOLFINx, to solve a non-homogeneous Helmholtz equation. We specifically address thermoacoustic Helmholtz equation.

The mesh is generated with Gmsh and the solver uses DOLFINx and UFL from FEniCSx. The performance, validity, stability and extensibility of the solver are demonstrated through several examples of thermoacoustic instability, from the one-dimensional Rijke tube to the three-dimensional MICCA combustor. The implementation of Bloch-type boundary conditions is explained and tested. The adjoint capability of the solver is also shown, and used to obtain derivatives of the eigenvalue with respect to shape parameters. This is exploited to find shape changes that reduce the thermoacoustic growth rate.

See README files in the /numerical_examples folder for detailed description of the thermoacoustic eigenmode and shape derivative calculations.