Simulation of mixing intensity profile for bioethanol production via two-step fermentation in an unbaffled agitator reactor

: Bioethanol synthesis techniques have been studied intensively due to the energy crisis and various environmental concerns. A two-step bioethanol production process was carried out multiple times in an unbaffled agitator tank. The parameters varied, including the fermentation temperature, the pH level, the amount of yeast, and the impeller type. Then, a simulation was used to obtain an image of the agitation behavior inside the agitator tank to compare the velocity profile of each type of impeller design. The impeller with eight blades was found to produce the highest flow velocity: 0.28 m/s. The highest concentration of bioethanol generated from the fermentation was 34 g/L, which was produced by using an eight-blade impeller at 30 °C, a pH level of 5, an agitation speed of 70 rpm, and 2 wt % yeast. The two-blade impeller produced the lowest bioethanol concentration, 18 g/L, under the same conditions. Ethanol concentration was found to peak at 40 °C and a pH level of 5. The geometry of the impeller, the fermentation temperature, and the pH level were each found to have a significant effect on the resulting bioethanol concentration according to the results of an ANOVA test. The amount of yeast had no effect on the fermentation reaction. Finally, the results demonstrated the possibility of using computational fluid dynamic modeling to determine the impeller’s behavior for the development of the bioethanol fermentation process. The simulation and experimental results from this research support the scaling up of a bioethanol production facility.


Introduction
The efficient production method for bio-fuel, such as bioethanol and biodiesel, from sustainable resources is now mandatory for humanity today and in the near future [1]. A shortage of petroleum reserves, increasing oil prices, and the degradation of the environment have driven countries around the world to seek scientific knowledge of how to produce more and better-quality biofuels [2]. In Thailand, the Ministry of Energy recognized the potential for bioethanol production from various biomass, such as sugar cane and agricultural wastes. Many financial subsidies from the government and Kyoto protocols have been issued as an incentive to establish a production base in Thailand that is capable of generating at least 11.9 million liters of bioethanol daily by the year 2036 according to the Alternative Energy Development Plan [3]. Bioethanol is produced by fermentation of reduced sugar (glucose, xylose, and fructose) in aqueous solution inhabited by saccharomyces cerevisiae yeast, which is also the kind of yeast that is used to make bread. However, in the case of agricultural waste or biomass fermentation it is important for the carbohydrate in these materials to first be broken down to fermentable sugar through a hydrolysis process in acidic solution.
In addition to biological activity, another important parameter for the fermentation process is the mechanical motion of the agitator. Mixing intensity in typical operating conditions depends on the formation of a vortex as it interacts with the liquid solution. A graphic representation of these interactions in the form of a vector is known as vorticity. Vortex depth is an important factor for a reaction that is affected by the air-to-liquid interface. This is especially relevant in the case of ethanol fermentation, which was found to improve significantly with the amount of oxygen in the reactor [4].
Additionally, in order to achieve a deep vortex, the water height-to-agitator tank diameter (h:d) should be greater than 1. Bimlesh Kumar et al. [5] demonstrated that the critical speed for the impeller's rotation is just enough for the vortex to reach the impeller. For small and medium-sized tanks, the critical speed is 60 rpm and 72 rpm, respectively. Under these conditions, and with a Reynolds number higher than 20,000, the turbulence flow system is usually adopted when simulating agitator and vortex behavior [6,7]. Basudeb Munshi et al. [8] reported the simulation of a mixing tank using the Reynolds Averaged Navier-Stokes (RANS) equations with k-ε turbulence to predict the surface pressure of flow inside the mixing tank. Since their main concern is to observe the gas-toliquid interaction, they proposed the volume of fluid (VOF) approach, which resulted in accurate simulation data. Recent studies have characterized molasses slurry as a non-Newtonian fluid that demonstrates shear thinning behavior during simulation based on computational fluid dynamics using a combined Eulerian-LES (Large Eddy Simulation) model [9].
One of the main obstacles in the fermentation of biomass is the viscosity of the biomass slurry. Energy intensive and time-consuming technologies have been adapted to help prepare feedstock in order to enhance fermentation. Consequently, bioreactors have to be designed to achieve efficient mass transfer, low shear stress force, and inexpensive material cost. For this reason, a continuousmode reactor with a bubble generator as its agitation method has been designed and referred to as a gas-lift bioreactor [9]. Plasma-assisted detoxification has been applied to treat sugarcane biomass [10]. Bae et al. [11] suggested a bamboo treatment with a hydrogen peroxide-acetic acid solution. A more convenient approach would be to improve bioethanol by understanding the contribution of agitation inside the bioreactor. Agitation speed and impeller geometry are important parameters that influence flow patterns, which prevent dead-zone areas (no-flow areas) in the case of a high viscosity liquid medium. The Reynolds number of a given flow was greatly affected by the impeller geometry. Marine blade and pitched paddle impellers were observed to provide flows with the highest Reynolds numbers (6000 and 8000) compared with both the Rushton turbine and Pitch Blade turbine [12]. In this experiment, the marine blade impeller offered superior biological activity even though its Reynolds number is smaller than that of the pitched paddle impeller. This is because the marine blade impeller implemented both radial and axial flow behavior. Ruihong Zhang et al. [13] reported a successful utilization of a hybrid helical ribbon and screw impeller for an aerobic compost bioreactor [13]. Currently, a new type of impeller with a curved blade (Scaba 6SRGT) was reported to significantly increase the well-mixed region but also increase power consumption [14]. Apart from measuring the improvement in biological activity due to the type of impeller, this research also monitors and takes into consideration the power consumption.
This research aimed to simulate the agitation profiles of impellers with different geometries, including pitched two-blade, pitched four-blade, Rushton four-blade, and modified propellers. An unbaffled agitator tank (constructed in the laboratory) was used to carry out the fermentation of treated Napier grass to bioethanol. An unbaffled tank was chosen because it was found to consume less power compared with a baffle tank [15]. Simulation of each type of impeller was controlled and executed by Simflow 3.1 software (SIMFLOW Technologies, Warsaw, Poland). The size of fresh Napier grass was mechanically reduced and treated using a solution of low-concentration sulfuric acid to liberate reduced sugar. The resulting slurry containing reduced sugar was then transferred and fermented in a separate reactor. Operating conditions such as temperature (30 °C, 40 °C, and 50 °C), amount of yeast (2.0 wt %, 3.5 wt %, and 5.0 wt %), and pH level (4, 5, and 6), were adjusted. This research also investigated the effect of additional acetic acid during the fermentation process. Finally, a statistical analysis of the data was obtained from the experiments to confirm the significance of the impeller's geometry and the operating conditions at a confidence level of 95%.

Experiment
The Napier grass was obtained from a plantation in Chaing Rai, Thailand. The biomass was stored in a low-temperature refrigerator (4 °C) to preserve the integrity of the biological structure. Prior to being used, the long Napier grass leaves were mechanically transformed to smaller pieces with sizes of 30-50 mm, as shown in Figure 1. A 15 wt % concentrated sulfuric acid was prepared and mixed with the biomass in a batch reactor at an agitation rate of less than 70 rpm for 5 h. The liquid slurry was then removed and passed through a vacuum filtration system in order to remove the liquid from the biomass, and it was washed with deionized water three times. The treated biomass was then transferred to another batch reactor for fermentation. Fermentation was carried out at different temperatures, amounts of yeast, and pH levels for 10 h. Each experimental trial was repeated three times to certify reproducibility. Calcium hydroxide was added in order to adjust the pH level of the fermentation solution. Acetic acid was added after fermentation to stop the bioethanol process. The concentration of bioethanol was sampled and analyzed every hour using a highperformance liquid chromatography system (Series 1100; Hewlett-Packard, Palo Alto, CA, USA) containing a model 168 UV detector. The injection temperature was set at 70 °C, the injection volume was 10 µL, and the mobile phase's flow rate was 0.9 mL/min sulfuric acid. The agitation reactor's dimensions are illustrated in Figure 2. Each operating condition was tested using a different type of impeller. The dimensions of each type of impeller used in this experiment are shown in Figure 3.

Governing Equation for the Simulation
Vorticity is a quantity representing the local rotational motion of the fluid under investigation or simply expressed as the curl of the velocity profile as shown in Equation (1). This value indicates the tendency of a specific microscopic point in the fluid to form a vortex. If the value of vorticity is not equal to 0, then the fluid is referred to as rotational. However, if the value of vorticity is equal to 0, then the fluid is considered as irrotational.
where = vorticity and v = velocity flow vector.
The velocity vector in both the x and y-direction can be derived from the Navier-Stokes equation as shown in Equation (2). According to the equation of motion, the two-dimensional unsteady form of the equation illustrated the relationship between dynamic parameters such as velocity, pressure, temperature, and density. However, for the derivation of vorticity, the pressure value at a specific location inside the dynamic flow system can be calculated using Equation (2).
where x, y, and z are the directions inside the Cartesian coordinate system; , , and are the vectors corresponding to each direction; and the f value is related to density of the flowing liquid.
The equation for conservation of momentum was employed to calculate the pressure of the water-air interface during agitation. Assuming transient behavior and non-slip conditions, the pressure profile in the radial direction of a cylinder shape can be developed as shown in Equation (3).
The continuity equation or the Navier-Stokes equation was employed to calculate the pressure along the height of the simulated cylinder, as shown in Equation (4).

Simulating Conditions
Simulation of the agitating behavior inside a cylindrical batch reactor was performed using SimFlow software, which was developed to simplify OpenFOAM software (OpenFOAM v8, The OpenFOAM Foundation Ltd., London, England). Three-dimensional geometries of the impellers were drawn using AutoCAD software (AutoCAD 2017, Autodesk Inc., Bangkok, Thailand) and imported to SimFlow. Meshing of the agitating system was conducted using a radial division of 15, radial grading of 1.02, and circumferential division of 15. Agitation was assumed to create a turbulence flow inside the cylindrical tank. For this reason, a RANS (Reynolds-averaged Navier-Stokes) turbulence equation with the realizable k-ε for a transient system was used. The agitation speed for the simulation was fixed at 6.28 rad/s (60 rpm). The two isosurfaces created 1 minute after agitation at the water-air interface and bottom of the cylindrical reactors were rendered using ParaView software (Microsoft, Redmond, WA, USA).

Simulated Agitation Pattern
The simulated agitation system inside an unbaffled reactor was rendered, as shown in Figure 4. The isosurface layers of the bottom and water-air interface were also illustrated. It was observed that the bottom isosurface of a two-blade impeller indicated a low-velocity flow profile with very little change in velocity. In contrast, the isosurface of an eight-blade impeller indicated a high-velocity flow profile. The flow gradient changed dramatically from a low-velocity flow profile to a high-velocity flow profile and then back to a lower-velocity flow profile. Figure 5 demonstrates the simulated combined velocity at specific location from the center of the impeller. The impeller's center corresponds to the origin of the graph (0, 0) while the other end (0, 18) represents the wall of the reactor tank according to the boundary conditions. The maximum flow velocity from the four types of impellers was achieved very close to the impeller's tip. The highest simulated flow velocity from the two-blade, four-blade, Rushton, and eight-blade turbines were 0.132, 0.208, 0.154, and 0.281 m/s, respectively. It is clear that the two-blade and Rushton turbines both produce flows with lower velocity compared with the four-blade and eight-blade turbines. The velocity profile of the Rushton blade was significantly higher than that of the two-blade turbine. These data correspond well with another experiment conducted for non-Newtonian fluids [16]. It was also observed that the eightblade turbine was the only type of impeller where the velocity of the fluid increased beyond the impeller's tip. This is due to the momentum created by the collision of vortex rings generated by the shear force of the pitched blade [17,18].

Effect of Operating Conditions on Bioethanol Fermentation
A two-step bioethanol production strategy was adopted for this experiment, including the primary hydrolysis reaction under acidic conditions to liberate reduced sugar inside of the biomass's structure. A fermentation reaction was promoted by the addition of acetic acid inside the solution in order to help digest the lignin present in most biomass materials [19]. Figure 6 demonstrates the effect of the type of impeller on the ethanol concentration of the fermented biofuel. After 40 h of fermentation, the eight-blade impeller was found to give the highest ethanol concentration of almost 34 g/L, followed by the four-blade impeller (30 g/L), the Rushton blade (20 g/L), and then the twoblade impeller (18 g/L). This is because of the increase in radial and axial mixing intensity as the number of pitched blades increases [20]. The ethanol concentration of the biofuel obtained in this experiment was slightly smaller than that of the fermentation of waste coke [21] and corn stover [22]. However, when an appropriate type of impeller was used, the ethanol concentration obtained from Napier grass for this experiment was higher than the concentration obtained from the fermentation of switchgrass [23]. This is because an ionic liquid was employed to break down the carbohydrate content inside the switchgrass instead of sulfuric acid, which was used in this experiment. Even though ionic liquids are currently less effective in releasing sugar from biomass compared with the chemical pretreatment method, they are environmentally friendly and can endure severe reaction conditions. Figure 7 demonstrates the effect of temperature on bioethanol fermentation. An increase in temperature from 30 °C to 40 °C resulted in an increase in ethanol concentration from 28.9 g/L to 34.3 g/L according to the samples collected after 40 h of fermentation. This is because an increase in temperature helps to promote lignin dissolution in the fermented solution. With lower lignin content, yeast habitants can easily penetrate to sugar molecules that are intertwined among larger macroscopic molecules, such as cellulose and hemicellulose. An increase in temperature also increases the kinetic energy, causing the fermentation rate to increase. However, when the system was exposed to higher temperatures, the ethanol concentration decreased gradually starting at 20 h after the beginning of fermentation. This result corresponds well with other research papers [24,25].
Techaparin et al. reported five thermotolerant yeasts that are capable of producing biofuel consisting of 48.51 g/L ethanol at a fermentation temperature of 40 °C [26].  In addition to the impeller type and fermentation temperature, the pH level was also varied to investigate its effect on bioethanol production from Napier grass. The pH level during fermentation is an important parameter because the pH level indicates the availability of protons in the system, which plays an important role in microbial activities [27]. An increase in the pH level from 4 to 5 resulted in a significant increase in bioethanol generation, as shown in Figure 8. When the pH level was increased further to 6, fermentation activity declined dramatically, producing biofuel with an ethanol concentration of only 18.9 g/L. These data align well with many other studies [28]. The amount of yeast added during fermentation does not affect the outcome of the fermentation process, as shown in Figure 9. However, an increase in yeast may affect power consumption because yeast increases the viscosity of the fermentation solution [29,30]. Additionally, yeast can created local dead spot inside the reactor with minimum heat transfer [31]. According to the results obtained from this experiment, it can be concluded that Napier grass can be used as second-generation feedstock for bioethanol production [32,33]. Comparisons with other studies on two-step ethanol production are shown in Table 1. According to Table 1, the bioethanol concentration of the fermentation product from Napier grass was found to be significantly higher than that of switchgrass.  A one-way analysis of variance (ANOVA) was conducted using Minitab 14.0 software (Minitab, Bangkok, Thailand) and data collected from the four types of impellers to statistically determine the significance of impeller geometry on bioethanol concentration. As illustrated in Table 2, the p-value is 0.0001, which is smaller than 0.05. Therefore, it was found that the type of impeller has a significant effect on the fermentation reaction. Additionally, the operating conditions were compared statistically using an ANOVA analysis of the regression model obtained by the Box-Behnken experimental design method. The temperature and pH level during fermentation were found to be statistically significant parameters. However, the amount of yeast added during fermentation was insignificant, as shown in Table 3.

Conclusions
The simulation of flow patterns produced from different types of impellers was performed in order to evaluate the microscopic velocity at a specific location in the radial direction of an agitation tank. Simulated results revealed the eight-blade impeller to provide the highest velocity (0.28 m/s) and the two-blade impeller to produce the lowest flow velocity (0.13 m/s). The flow velocity in the agitated area has a significant effect on the ethanol concentration of the produced biofuel. The twostep bioethanol production experiment revealed the effect of impeller type, fermentation temperature, pH level, and the amount of yeast on the ethanol concentration of the obtained biofuel. An ANOVA analysis revealed that the type of impeller, fermentation temperature, and pH level each have a significant effect on the bioethanol concentration in the product.