OpenFOAM® (OF) is gaining ever-increasing momentum for simulating reactive flows like combustion. However, the accuracy of capturing the complex interaction between chemistry and transport is hardly studied in the literature. A first theoretical investigation of OF’s operator splitting to dived chemistry and transport has been studied in a recent work by the author and revealed that OF cannot capture near-ignition and near-extinction events and that steady-state results are influenced by time step size. This work investigates strategies to improve OF’s operator splitting algorithm. A well-stirred scalar test case and combustion test case known to be challenging for operator splitting schemes are used to evaluate the accuracy and steady-state preserving properties of the OF and variations of the OF splitting scheme. Additionally, the stability properties and steady-state preservence are evaluated analytically. The analytical analysis revealed that OF employs a source term linearization approach instead of a proper operator splitting scheme. The investigated OF splitting variations show similar properties to the original one. The investigations reveal that the OF splitting scheme is inferior to well-established schemes like the Strang or Simpler splitting scheme. However, implementing these classical splitting schemes is troublesome due to OF’s code structure. At the current code base, variations of the OF linearization approach must be used to improve OF’s reactive flow modeling capabilities.