Reactive bubbly flows are found in many chemical and biochemical processes. They are characterized by complex hydrodynamics that govern global mass transfer and reaction rates. Effects, which have to be taken into account when modelling and simulating the reaction progress, are enhanced bubble–bubble interaction and bubble-induced turbulence as well as swarm-induced macro-convection. This poses great challenges on both, experimental analysis and numerical simulation. Bubble swarms at higher gas fractions are opaque, which limits the use of optical flow measurement techniques. Moreover, simulations at industrial scale are only feasible with point-bubble-based Euler-Euler or Euler–Lagrange approaches. Such approaches require closure relations, which account for all relevant interfacial forces as well as bubble-induced turbulence and coalescence including swarm effects. This chapter describes the progress in the experimental analysis and CFD simulation of reactive bubbly flows achieved during the time of the DFG Priority Programme SPP 1740. Experimental studies in bubble columns reported in this chapter were carried out at TU Dresden and OVGU Magdeburg. CFD model development and simulation has been carried out at OVGU Magdeburg using the Euler-Euler and Euler–Lagrange method and at TU Kaiserslautern and Helmholtz-Zentrum Dresden-Rossendorf using the Euler-Euler approach. As reaction systems the chemisorption of CO₂ in NaOH solution and the reaction of Fe^II(edta) with NO were studied.