Abstract: The use of microfluidic devices in the field of plasma-liquid interaction can unlock unique possibilities to investigate the effects of plasma-generated reactive species for environmental and biomedical applications. So far, very little simulation work has been performed on microfluidic devices in contact with a plasma source. We report on the modelling and computational simulation of physical and chemical processes taking place in a novel plasma-microfluidic platform. The main production and transport pathways of reactive species both in plasma and liquid are modelled by a novel modelling approach that combines 0D chemical kinetics and 2D transport mechanisms. This combined approach, applicable to systems where the transport of chemical species occurs in unidirectional flows at high Péclet numbers, decreases calculation times considerably compared to regular 2D simulations. It takes advantage of the low computational time of the 0D reaction models while providing spatial information through multiple plug-flow simulations to yield a quasi-2D model. The gas and liquid flow profiles are simulated entirely in 2D, together with the chemical reactions and transport of key chemical species. The model correctly predicts increased transport of hydrogen peroxide into the liquid when the microfluidic opening is placed inside the plasma effluent region, as opposed to inside the plasma region itself. Furthermore, the modelled hydrogen peroxide production and transport in the microfluidic liquid differs by less than 50% compared with experimental results. To explain this discrepancy, the limits of the 0D–2D combined approach are discussed.

Abstract: Plasma-treated liquids have great potential for biomedical applications. However, insight into the underlying mechanisms and the exact chemistry is still scarce. In this study, we present the combination of a 0D chemical kinetics and a 2D fluid dynamics model to investigate the plasma treatment of a buffered water solution with the kINPen (R) plasma jet. Using this model, we calculated the gas and liquid flow profiles and the transport and chemistry of all species in the gas and the liquid phase. Moreover, we evaluated the stability of the reactive oxygen and nitrogen species after plasma treatment. We found that of all species, only H2O2, HNO2/NO2-, and HNO3/NO3- are stable in the buffered solution after plasma treatment. This is because both their production and loss processes in the liquid phase are dependent on short-lived radicals (e.g. OH, NO, and NO2). Apart from some discrepancy in the absolute values of the concentrations, which can be explained by the model, all general trends and observations in our model are in qualitative agreement with experimental data and literature.