3.1. Computational fluid dynamics of reactors
The CFD models of CGTR, STR and CAR have been developed using the commercial software Ansys Student software 2022 (Ansys Inc.). The reactor models of respective dimensions were defined (Kaviya and Kanmani 2022) and modelled using the AnsysSpaceclaim tool (Fig. 3).
The inlet and outlet tubes were modelled for the three reactor models as they have a great impact on the flow field. The quartz tubes in CAR were simulated as an empty cylinder. The meshed reactors are displayed in Fig. 4.
The Element size was defined as 2mm and high quality smoothening was utilized for building mesh. The reactor volume was discretized in approximately 940000 structured and unstructured volume cells using Ansys Meshing tool. This number of cells has been found to be high enough to give mesh-independent results, corresponding to a mean cell volume of 0.01 mm3. A denser mesh was defined at the inlet and outlet as there may be greater velocity gradient at both inlet and outlet areas.
The Fluent analysis was used to predict the hydrodynamic behaviour of the three photocatalytic reactors. Figure 5 shows the streamlines of velocity magnitude calculated for the CGTR, STR and CAR assuming a steady state model. In case of CGTR, the streamlines of flow are disturbed at the inlet due to the sudden change in flow direction from vertical inlet to circular pathway and further through the flow the streamlines become smooth and even. It is observed that in CGTR, there is minimal amount of transverse mixing. The Reynolds number (Re) of CGTR is calculated to be 300, which makes it clear that the flow is laminar. The low Re is due to the less characteristic length of the micro-depth reactor of 2mm (0.002m). The maximum velocity (0.15m/s) was taken for the calculation of Re for all reactors. The theoretical narrow laminar profile (in RTD curve) was not reached due to the axial dispersion.
It was observed that at CGTR inlet (Fig. 6), there was small eddies causing some extent of mixing and further the flow velocity became laminar throughout the reactor. In case of STR, initially the flow is even and later the flow becomes transient. There is dead zone found at the junction between glass tube and inlet/outlet, where the velocity is lowered due to the sudden increase in the cross-section of flow. The characteristic length dimension of STR is calculated to be 0.0175m. Reynolds number of STR is calculated to be 2625, which makes it clear that the flow is a transition between laminar and turbulent flow. In the case of CAR, it was observed that velocity streamlines is non-uniform through the annular region. There are dead zones found at the walls of reactor and at the walls of quartz tube where velocity magnitude is very low. The Re was calculated to be 9300, which is very high implying that CAR follows turbulent flow. The characteristic length of CAR was calculated to be 0.062cm.Similarly Cintia et al. (2016) work was related to the hydrodynamic modelling of annular photocatalytic reactor, which demonstrated the existence of dead and recirculation zones along with non-uniform flow through the reactor domain. The longitudinal dimensions of CGTR, STR and CAR are 0.6m, 0.24m and 0.14m respectively. The longer longitudinal dimension enhances the laminar flow and shorter longitudinal dimension enhance turbulent flow (Xiao 2019).
3.2. Residence time distribution
The study of hydrodynamics of reactors is very important for the reactions involving fluids. The residence time distribution indicates the type of chemical reactor and Reynolds number indicates the type of flow.
To confirm the hydrodynamic behaviour and flow pattern of the three reactors, Residence time distribution (RTD) approach was performed. One of the most important attributes that need to be considered while examining the hydrodynamic behaviour of the reactor is residence time distribution (RTD) using Methylene blue tracer. The Fig. 7 shows the concentration curve of RTD data, where there is a short delay during which little or no tracer was detected at the outlet till 4 seconds. At 5th second, there was a sharp rise in tracer concentration followed by a steady drop in tracer concentration. The flow behaviour of CGTR indicated ideal plug flow model.
Further Fig. 8 shows the Exit age distribution E(t) for three types of reactor from experimental analysis. The RTD of CGTR denoted that all the tracer fluid elements systematically passed the reactor without transverse mixing. The RTD of CAR indicated the profile of ideal continuously stirred tank reactor which is characterised by perfect back mixing. At t = 0 second, the tracer concentration was identified at the outlet. This showed that CAR is perfectly mixed flow and each portion of the fluid elements has the same chance to be discharged at the outlet, regardless of how long it has already been inside the continuously stirred tank reactor. There is a slow decrease in the tracer concentration at the outlet till 4 seconds.
The RTD of STR showed that there was sudden rise in tracer concentration at outlet at 1 second, later there was a decline in the concentration of tracer till 4 seconds. There was some degree of plug flow due to the initial sharp peak and there was also some degree of mixed flow due to the slow decrease in concentration. This showed that STR was operated in between plug flow reactorand continuously stirred tank reactor. This type of flow is characterised as arbitrary flow (Dong and Zongshu 2012). Residence time and variance values were calculated as reported in Sahleet al. (2003).The mean residence time (T) were calculated to be 5.2 sec (σ2 = 0.45s), 2.9 sec (σ2 = 2s) and 1.5 sec (σ2 = 2.5s). The theoretical residence time was kept constant as 7.98 sec (0.133 min). The difference in the theoretical and mean residence time is due to the limited back mixing and presence of dead zones in case of STR and CAR (Peter et al. 2019). Similar RTD study was done by Ina et al. (2021) for fixed bed reactor, where plug flow behaviour was obtained with least dispersion.
3.3. Comparison of reactors
In our previous work (Kaviya and Kanmani 2022), it was shown that Photocatalysiswith simultaneous adsorption and hydrodynamic caviation (PCSA + HC) exhibited the maximum colour removal of 100% in CGTR. Similarlyunderoptimised conditions, experiments (1. photocatalysis with initial adsorption - PCIA, 2. photocatalysis with initial adsorption and hydrodynamic cavitation - PCIA + HC, 3. Photocatalysis under simultaneous adsorption – PCSA, 4. Photocatalysis under simultaneous adsorption and hydrodynamic cavitation - PCSA + HC) including control experiments (Photolysis and hydrodynamic cavitation- HC)were performed for STR and CAR.The % colour removal achieved in the three reactors is illustrated in Fig. 9. The experiments optimised conditions were initial pH 11 and graphene dosage 0.3g/L, AO dye concentration 10 ppm and irradiation time 35 min (Kaviya and Kanmani 2022). The experiments were carried out in triplicate to obtain the standard deviation error of < 5%.
The STR exhibited the second highest removal (85.8%) under PCSA + HC. In CAR, 57% of colour was removed under PCSA + HC. These results implied that CGTR performed better than STR and CAR. Under PCSA + HC, there is synergetic effect of photocatalysis, adsorption and hydrodynamic cavitation for increased degradation rate (Esrafil et al. 2021).
To get a better understanding of the performance and energy efficiency among Micro-reactor (MR) (Aura et al. 2014), Packed bed reactor (PBR) [1], CAR, STR and CGTR, the modified Space time yield (STYmodified) and modified Photocatalytic space time yield (PSTYmodified)is plotted in Fig. 10 and represented in Table 1.
Table 1
Calculated modified performance parameters
Reactors | Characteristic length (m) | Velocity (m/s) | Volume throughput (L) | STYmodified (s− 1) | PSTYmodified (s− 1/kW/m3) |
CAR | 0.062 | 0.15 | 1 | 0.04 | 0.0017 |
STR | 0.017 | 0.15 | 1 | 85.7 | 0.357 |
CGTR | 0.002 | 0.15 | 1 | 225 | 9.375 |
MR | 0.00005 | 0.0067 | 1.46x10− 6 | 133 | 16.625x10− 7 |
PBR | 0.001 | 0.0036 | 0.15 | 3.6 | 0.0132 |
Space time yield and Photocatalytic space time yield is a benchmark proposed by Enis et al. (2015), which is dependent on the hydrodynamic behaviour (either continuously stirred tank reactor or plug flow reactor) and reaction rate constant. Space time yield is dependent on reaction rate constant and hydrodynamic flow behavior whereas Photocatalytic space time yield is dependent on Space time yield and illumination efficiency. The modified parameters STYmodified and PSTYmodified are calculated using the hydrodynamic flow behavior influencial parameter, which is called characteristic length of photocatalytic reactor. The reactor which performs best in pollutant removal (productivity) is present on the top side of graph whereas a reactor which performs best in energy utilisation is present on the right side of graph.
It was observed from Fig. 10that CAR and PBR (Thomas et al. 2019) are positioned near the origin. This implies that these reactors have less productivity (less STYmodified) and less energy efficiency (less PSTYmodified). The less values of STYmodified of CAR and PBR is attributed to fluid velocity and characteristic length, where the higher the characteristic length and lower the fluid velocity resulted in the lower values ofSTYmodified. On the other hand, PSTYmodified depends on the STYmodified and lamp power required to illuminate the volume of reactor.Also the reactors which are operated in recirculation has higher modified parameters compared to continuously operated reactors. STR has lowercharacteristic length and higher fluid velocity compared to CAR impling higher modified parameters. Further the position of MR (Aura et al. 2014) in the right bottom in the graph (high STYmodified and less PSTYmodified), reveals that the fluid velocity also has significant influence on modified parameters eventhough the characteristic length is least in MR (Aura et al. 2014). In the case of CGTR, both high productivity and high energy efficiency is achieved due to its highfluid velocityand less characteristic lengthand operated in recirculation process. Therefore this study showed that CGTR was the most efficient in term of both performance and energy efficiency.
The interpretations of modified parameters revealed that the maximum modified Space Time Yield (STYmodified) 225 s− 1and maximum modified PhotocatalyticSpace Time Yield (PSTYmodified) 9.375 s− 1/kW/m3 was obtained for plug flow reactor with lesser characteristic length when compared to continuously stirred tank reactor (CAR) and arbitrary flow reactor (STR) which has higher characteristic length. The treated throughput (1L) and fluid velocity (0.15m/s) is kept constant in CAR, STR and CGTR.