Historically, CFD codes implementing arbitrary Lagrangian Eulerian interface tracking (ALEIT) for the simulation of bubbly flows were highly customised for this specific application. However, taking a broader look at multiphase flows in general, it is noticeable that they resemble other multiphysics systems in their basic structure. Analogous to Fluid-Structure-Interaction (FSI) or conjugate heat transfer problems (CHT), multiphase systems also involve several regions with special physical properties that interact with each other across interfaces. Exploiting these structural similarities, a novel unified multiphysics framework for OpenFOAM named multiRegionFoam is introduced, which incorporates the ALE-IT and is tested for its application towards reactive bubbly flows. In this context, the parallelisability of the new framework is addressed in particular, since interface coupled multiphysics simulations in general and simulations of reactive bubbly flows in particular are very computationally intensive. Furthermore, a Dirichlet-Neumann algorithm convergence control based on interface residuals is implemented and tested.