Comparative Analysis of Hydrodynamics Behavior of Microalgae Suspension Flow in Circular , Square and Hexagonal Shape Photo Bioreactors

Microalgae based biofuel is an emerging natural source of energy alternative to the fossil fuel. As microalgae are photosynthetic microorganisms, light is one of the limiting factors for its culture. Though many researches have been carried out for findings behind suitable culture system for the proper growth of microalgae, those are confined only to tubular Photo-bioreactor (PBR). This paper aims to make comparison among the horizontal loop photo bioreactors with different cross sections based on the analysis of hydrodynamics behavior. Three different geometrical shapes having vertical cross sections of circular, square and hexagonal PBR, have been proposed taking into account light intensity for microalgae culture. In this study, we simulate the flow dynamics of three types of PBRs and discuss the velocity, pressure and shear stress properties as microalgae endurance capacity depends on them. For the dimension of the three PBRs we considered here, each of them have radius of about 0.05 m while the length together with bending portion is approximately 20.5 m for a single loop. From the study, the hydrodynamic behaviors are observed to be quite dissimilar in case of three PBR’s. In the straight portion the velocity profile is quite parabolic in tubular but distorted minimally in case of square and hexagonal PBRs. In the middle of the U-loop, a haphazard fluid distribution is noticed. The velocity magnitude and agitation of microalgae cells are higher in hexagonal than in square and tubular. The shear rate is less in case of tubular compared to square and hexagonal. A linear pressure drop is found from the inlet to the outlet for three PBR’s. From this comparison, it can be said that the tubular one would be the best option for microalgae culture in case of industrial purposes.


Introduction
Sustainable development and efficient use of energy goes on hand which results in environment pollution minimisation and better socio-economic conditions.The Green energy is now a prime talk all around the globe to ensure zero toxic gas emission as global warming and other hazardous pollution caused by burning of fossil fuels are taking the world to an unstable condition.Also increasing consumption of petro-fuel with rapid growth of transportation and population compared to total deposited amount is leading it on the verge of gradual extinction as it is a depleted source of energy [1] [2].So in this context, oil extracted from organic matter (biofuel) has attracted the scientist and researchers as it is eco-friendly and sustainable.Biofuel has a great potential to mitigate the continuous demand in every aspects.The mission of negative carbon emission would be fulfilled to a great extent by large scale production of biofuel.
Biofuel production may become much more economical than fossil fuel.Extensive research has been carried out for the last five decades regarding the development of the biofuel technology.The developed countries of the world have emphasized and investing money to increase the productivity level of biofuel.So it can be foreseen to be the fuel of future [3] [4] [5].
Biofuel is referred to as solid, liquid or gaseous compound obtained from organic matter.The main advantage of it over fossil fuel is that it is non-toxic and biodegradable [2].Biofuel production is categorized into three generations.The first generation biofuel such as soybean, rapeseed, sunflower, palm oil has some environmental limitations as they compete with cultivable land for food production.The second generation biofuel such as agricultural residues, wood residues, non-edible oil or sugarcane has issues of environmental hazard [6] [7].So to overcome those drawbacks microalgae based biofuel has been appeared as a potential source of biodiesel which ensures the clean energy for better environment.
Microalgae are eukaryotic and prokaryotic photosynthetic microorganisms.They convert water and CO 2 into sugar i.e. lipid and carbohydrate in photosynthesis process by means of sunlight.The sugar contents are subsequently used to extract oil.In the very beginning, microalgae have harvested for some pharmaceutical purposes, waste water treatment, cosmetics and poultry food.As the time passing by, it has got importance in fuel sector due to its almost double productivity level of biomass and it grows 100 times faster than any terrestrial plants.So the challenging task for the researchers is to develop congenial culture system for high productivity of microalgae [8] [9] [10] [11].
Microalgae culture system requires supply of light, CO 2 , nutrients but availability of light is the first and foremost concern.Adequate exposure to light of every cell in the culture system ensures the uniform growth of microalgae, yet it is a challenging task due to proper design consideration of culture system.Generally, two types of culture systems are available such as traditional open pond system and closed photo bioreactor system.The open pond includes shallow big ponds, circular ponds, tanks and raceway pond but raceway pond has been commonly used for culture [12] [13].The raceway pond method facilitates economic and easier process of culture but has a lot of disadvantages include poor stirring mechanism so mass transfer rate is very low.Hence less mass productivity, insufficient utilization of sunlight close to the bottom surface of pond and contamination risk due to fast growing heterotrophs and predators have hindered the commercial production.On the contrary, the closed photo bioreactor technology is free from those drawbacks to a great extent.But, photobioreactor which is promising to culture of microalgae is still under consideration as better geometric shape and proper growth model are associated with it [14].Maximum sunlight capturing capacity and minimum space requirement are first and foremost condition in designing a photo bioreactor.Both the outdoor and the indoor cultivation process are in vogue but for large scale production outdoor cultivation is required as to utilize natural source of light.The photo bioreactors such as bubble column, torus, helical and stirred tank have low illuminated surface area.
But for outdoor cultivation large illuminated surface area is needed.In this context, U-loop horizontal tubular photo bioreactor is more efficient than other considering their large illumination area and easily scaling up capacity [15] [16] [17] [18].From this point of view, investigations were limited only to analyse the growth model in case of horizontal u-loop tubular photo bioreactor.In this work we have conducted CFD comparison analysis among different shape horizontal photo bioreactors as simulation plays a vital role to predict the result before going to experiment on a real test bed.We hereby proposed vertically square and hexagonal shape horizontal PBRs besides tubular PBR.As hydrodynamics, light regime, shear stress and pressure distribution are most important parameters to decide which one is the best suited to microalgae growth, so we consider these parameters for result analysis.
In our simulation only the microalgae suspension is considered to obtain fully developed single phase laminar flow with no slip condition at the wall.Here, we have focused on the comparison of fluid behaviour in U-loop portion other than straight part as shear rate and rate of cell damage is high there.
The rest of the paper is organized as follows.The theoretical framework is described in the Section 2; Section 3 describes the methodology of the mathematical model development; results and conclusions are presented in the Section 4 and the Section 5 respectively.

Creeping Flow
where ρ is the density; u is the velocity vector; p is the pressure; τ is the viscous stress tensor and F is the body force vector.

Shear Stress
In case of all three kinds of PBR, hydrodynamic forces will obviously cause high shear stress which affects the culture of microalgae by causing constraint to its flow and growth which results in huge cell damage [5].Shear stress will also determine which PBR is conforming to microalgae culture.The shear stress is determined by the following equation.
where τ, u, y denote shear stress, velocity of flow and direction normal to the flow respectively.

Methodology
As in this paper, our focus is on the comparative analysis of flow phenomenon for three PBRs, thus we conduct the study in two steps: Mathematical model development and numerical simulation.

Geometry
Each of the three profiles has same hydraulic diameter, length and radius of the curvature.The radius is 0.025 m, length is approximately 20.4 m and radius of curvature at U-loop portion is 0.4 m.As all the parameters involved with geometric construction are constant so we can build up a common sketch for tubular, square and hexagonal shape PBR which is shown in Figure 1.In case of the tubular PBR the hydraulic diameter is equal to the diameter of profile.For calculating the hydraulic diameter for square and hexagonal shape the following equation is used where D h is the hydraulic diameter; A is the area of cross section; s is the wetted perimeter.

Computational Domain Development
For computation every domain is placed along x-y plane horizontally.Z axis is perpendicular to flow direction.Table 1 shows the faces, edges and intersecting points for the three PBR's.
The surface area and working volume are shown in Table 2.
The domains of three PBR's are presented in Figure 2.

Governing Equations
As temperature variation is low and the density is constant, microalgae suspension is considered to be incompressible fluid so Equation (1) reduces to and Equation ( 2) becomes ( ) where σ is the stress tensor and g is the gravity, σ can be expressed as ( ) where, η = viscosity of the fluid; D(v) = rate of deformation.The viscosity η(t) in Equation (7) is determined by The relative viscosity ( ) r t η relating to the concentration is then used and determined by ( ) ( ) where ε is the Einstein's coefficient [19].Based on the experiment conducted by Hon-Nami and Kunito [20], the concentration function C(t) in Equation ( 9) is given by the logistic Equation ( 10) where µ is the constant growth rate of microalgae cells; C 0 is the initial concentration of the suspension and A and B are constants.

Boundary and Initial Conditions
For simulation we considered the no slip condition on the wall of the tube and the zero normal stress at the outlet for all the three domains, as follows:

Mesh Generation
To implement the Navier-Stokes equation with incompressible flow in microalgae suspension, mesh generation is required for calculation.In our study we use normal mesh.
Table 3 shows the mesh elements, minimum and average quality for three PBR's.
The numbers of vertex, edge and boundary elements are given in Table 4.

Grid Sensitivity Analysis
A grid sensitivity test is performed in case of tubular photo bioreactor.Both normal and coarse mesh is formed to make comparative study of mesh quality for better result in case of time dependent study.Table 5 shows the comparison of elements between normal and coarse mesh.
From Table 5 we observed that mesh quality and elements both are better in normal mesh.

Simulation Parameters
The main goal of our study is to acquire in depth knowledge of flow behaviour for three different shape PBR's.For our simulation we use the COMSOL Multiphysics version 4.2a package.The parameters that are used as input data are given in Table 6.

Results and Discussion
During the simulation, the solver was configured as time dependent.For achieving better and comparative results, we ran the simulation for the seventh day of microalgae culture as microalgae growth can clearly be observed from this day.The time range was (540,000, 10, 540,050) seconds.From those analyses significant changes are noticed for three PBR's.For comparison, we analyse the results for 50 second.The topics of interest in our result analysis include grid independency, velocity distributions, pressure profiles through the ducts curvature and shear stress on the wall of ducts.The shear stress at the middle of the U-loop is the prime concern in this paper as it helps to find out which one is conforming to less cell damage.To identify the fluid behaviour for different profiles the three PBR's are fragmented in different three cross sections.
As mesh size plays a vital role in case of accuracy, so satisfactory computational accuracy can be achieved by continuously changing the meshes until the results from two trials lead to very close to each other [14].As we previously discussed about the grid independency, some results regarding velocity magnitude are shown in Figure 4 and Figure 5 for normal and coarse mesh in case of tubular PBR.
In Figure 4(a), for coarse mesh the velocity magnitude is 0.9168 m/s and in Figure 4(b), for coarse mesh it is 0.9316 m/s.The both cross sections are taken at the middle of U-loop.From Figure 5(a) the maximum velocity at straight portion is 0.75 for coarse mesh and from Figure 5(b) it is 0.8 for normal mesh.

Velocity Profile
The velocity profile helps to predict the fluid behaviour in the three PBR's.Higher velocity magnitude makes a haphazard distribution of fluid in U-loop in case of all ducts.
But the challenging task is to investigate which one has comparatively less velocity magnitude.In Figure 6, the velocity profiles along XY-plane are shown of three PBR's at 50 second respectively.
From Table 7 it is observed that at 50 s lower velocity magnitude is noticed in case of the tubular PBR whereas the square and the hexagonal show almost similar behaviour and higher value than the tubular with inlet velocity 0.5 m/s.
In Figure 7      that higher agitation and speed of particles region is adjacent to the wall of small radius of curvature for three ducts in U-loop portion.For 50 s in every case, tubular one shows less speed than square and hexagonal.
The velocity distribution graphs for the tubular, the square and the hexagonal shape

Shear Stress
As the shear stress distribution indicates to predict which one will show better performance to lessen cell damage, we have analysed the shear rate for three different cross sections of U-loop which are presented in Figure 9.
By comparing three different geometrical shape of PBR, less shear stress is observed for the tubular PBR.Table 8 shows the maximum shear stress for different cross sections of three PBR's.

Pressure Distribution
The pressure profiles are uniform from inlet to outlet for all three types of geometry.
From Figure 10, we can see pressure along the ducts in XY-plane of tubular PBR.Table 9 shows the value of maximum pressures.The pressure is slightly fluctuates at the middle of U-loop for all PBR's which are also observed from line graph shown in Fig-

Conclusion
In our study, we simulate the flow dynamics of three types of PBR's and discuss the velocity, pressure and shear stress properties as microalgae endurance capacity depends on them.For all the cases, in the U-loop portion, higher velocity exists than any other parts but always the speed is less and moderate in tubular PBR than others.The velocity distribution in tubular PBR is better suited to the culture of microalgae.As shear stress is mostly important factor for microalgae culture, we have analysed our results especially for U-loop portion and less shear stress is found in tubular shape than rest The creeping flow model hereby used in our simulation to explain the flow behaviour is a branch of single phase flow and works with fluid flow having very low Reynolds number.In case of creeping flow the inertia term of the Navier-Stokes equation has been neglected.As it occurs in fluid systems with high viscosity and micro scale geometry, so it can be used to analyse the flow behaviour of microalgae.The single phase creeping flow model follows the Navier-Stokes equation which is in general form as follows:

Figure 1 .Figure 2 .
Figure 1.(a) Dimension of length of straight portion of the PBR (tubular, square, hexagonal); (b) Inner and outer radius of the curvature; (c) Space between straight portions.

Figure 3 (
Figure 3(a) shows the mesh design of the U-loop of tubular, square and hexagonal shape PBR's and Figure 3(b) shows the vertical cross sectional views of the corresponding PBR's.

Figure 4 .
Figure 4. Velocity at the middle of U-loop of tubular PBR for (a) coarse mesh (b) normal mesh.

Figure 5 .
Figure 5. Velocity profile at straight portion of tubular PBR for (a) coarse mesh (b) normal mesh.

PBR are shown in Figure 8 .
The graphs are taken at time 50 s.The velocity distribution is completely parabolic in the straight portion of every PBR but at the entrance of U-loop, middle of U-loop and outlet of U-loop the shapes are distorted from its regular parabolic shape as the turbulence create in the U-loop portion.

Table 1 .
Faces, edges and points of the PBR's.

Table 2 .
Surface area & working volume of PBR's.

Table 4 .
Mesh parameters of the PBR's.

Table 5 .
Comparison between normal and coarse mesh.

Table 6 .
Parameters used for simulation.

Table 7 .
Maximum velocity magnitude at 50 s.

Table 8 .
Minimum shear stress for three PBR's.