<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AJCM</journal-id><journal-title-group><journal-title>American Journal of Computational Mathematics</journal-title></journal-title-group><issn pub-type="epub">2161-1203</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ajcm.2016.64033</article-id><article-id pub-id-type="publisher-id">AJCM-72910</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Comparative Analysis of Hydrodynamics Behavior of Microalgae Suspension Flow in Circular, Square and Hexagonal Shape Photo Bioreactors
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mushfique</surname><given-names>Shahriar</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mohammad</surname><given-names>Iftekhar Monir</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ujjwal</surname><given-names>Kumar Deb</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Mathematics, Premier University, Chittagong, Bangladesh</addr-line></aff><aff id="aff1"><addr-line>Department of Mechanical Engineering, Chittagong University of Engineering &amp;amp; Technology, Chittagong, Bangladesh</addr-line></aff><aff id="aff3"><addr-line>Department of Mathematics, Chittagong University of Engineering &amp;amp; Technology, Chittagong, Bangladesh</addr-line></aff><pub-date pub-type="epub"><day>11</day><month>11</month><year>2016</year></pub-date><volume>06</volume><issue>04</issue><fpage>320</fpage><lpage>335</lpage><history><date date-type="received"><day>October</day>	<month>17,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>December</month>	<year>18,</year>	</date><date date-type="accepted"><day>December</day>	<month>21,</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  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.
 
</p></abstract><kwd-group><kwd>CFD</kwd><kwd> Microalgae</kwd><kwd> Biofuel</kwd><kwd> Photo Bioreactor</kwd><kwd> Simulation</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>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 [<xref ref-type="bibr" rid="scirp.72910-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref2">2</xref>] . 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 [<xref ref-type="bibr" rid="scirp.72910-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref5">5</xref>] .</p><p>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 [<xref ref-type="bibr" rid="scirp.72910-ref2">2</xref>] . 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 [<xref ref-type="bibr" rid="scirp.72910-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref7">7</xref>] . 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.</p><p>Microalgae are eukaryotic and prokaryotic photosynthetic microorganisms. They convert water and CO<sub>2</sub> 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 [<xref ref-type="bibr" rid="scirp.72910-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref10">10</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref11">11</xref>] .</p><p>Microalgae culture system requires supply of light, CO<sub>2</sub>, 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 [<xref ref-type="bibr" rid="scirp.72910-ref12">12</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref13">13</xref>] . 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 [<xref ref-type="bibr" rid="scirp.72910-ref14">14</xref>] . 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 [<xref ref-type="bibr" rid="scirp.72910-ref15">15</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref16">16</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.72910-ref18">18</xref>] . 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.</p><p>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.</p><p>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.</p></sec><sec id="s2"><title>2. Theoretical Framework</title><sec id="s2_1"><title>2.1. Creeping Flow</title><p>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.</p><p>The single phase creeping flow model follows the Navier-Stokes equation which is in general form as follows:</p><disp-formula id="scirp.72910-formula8"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x2.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.72910-formula9"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x3.png"  xlink:type="simple"/></disp-formula><p>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.</p></sec><sec id="s2_2"><title>2.2. Shear Stress</title><p>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 [<xref ref-type="bibr" rid="scirp.72910-ref5">5</xref>] . Shear stress will also determine which PBR is conforming to microalgae culture. The shear stress is determined by the following equation.</p><disp-formula id="scirp.72910-formula10"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x4.png"  xlink:type="simple"/></disp-formula><p>where τ, u, y denote shear stress, velocity of flow and direction normal to the flow respectively.</p></sec></sec><sec id="s3"><title>3. Methodology</title><p>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.</p><sec id="s3_1"><title>3.1. Mathematical Model</title><sec id="s3_1_1"><title>3.1.1. Geometry</title><p>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 <xref ref-type="fig" rid="fig1">Figure 1</xref>. 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</p><disp-formula id="scirp.72910-formula11"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x5.png"  xlink:type="simple"/></disp-formula><p>where D<sub>h</sub> is the hydraulic diameter; A is the area of cross section; s is the wetted perimeter.</p></sec><sec id="s3_1_2"><title>3.1.2. Computational Domain Development</title><p>For computation every domain is placed along x-y plane horizontally. Z axis is perpendicular to flow direction. <xref ref-type="table" rid="table1">Table 1</xref> shows the faces, edges and intersecting points for the three PBR’s.</p><p>The surface area and working volume are shown in <xref ref-type="table" rid="table2">Table 2</xref>.</p><p>The domains of three PBR’s are presented in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> (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</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x6.png"/></fig><fig-group id="fig2"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> (a) Tubular PBR showing inlet, outlet and U-loop; (b) Square PBR showing inlet, outlet and U-loop; (c) Hexagonal PBR showing inlet, outlet and U-loop.</title></caption><fig id ="fig2_1"><label>(b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x7.png"/></fig><fig id ="fig2_2"><label>(c)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x8.png"/></fig><fig id ="fig2_3"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x9.png"/></fig></fig-group><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Faces, edges and points of the PBR’s</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >PBR’s shape type</th><th align="center" valign="middle" >Faces</th><th align="center" valign="middle" >Edges</th><th align="center" valign="middle" >Points</th></tr></thead><tr><td align="center" valign="middle" >Tubular</td><td align="center" valign="middle" >8</td><td align="center" valign="middle" >18</td><td align="center" valign="middle" >12</td></tr><tr><td align="center" valign="middle" >Square</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >24</td><td align="center" valign="middle" >16</td></tr><tr><td align="center" valign="middle" >Hexagonal</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >42</td><td align="center" valign="middle" >24</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Surface area &amp; working volume of PBR’s</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >PBR’s shape type</th><th align="center" valign="middle" >Surface area (m<sup>2</sup>)</th><th align="center" valign="middle" >Working volume (m<sup>3</sup>)</th></tr></thead><tr><td align="center" valign="middle" >Tubular</td><td align="center" valign="middle" >6.351</td><td align="center" valign="middle" >0.1567</td></tr><tr><td align="center" valign="middle" >Square</td><td align="center" valign="middle" >8.124</td><td align="center" valign="middle" >0.2026</td></tr><tr><td align="center" valign="middle" >Hexagonal</td><td align="center" valign="middle" >7.036</td><td align="center" valign="middle" >0.1755</td></tr></tbody></table></table-wrap></sec><sec id="s3_1_3"><title>3.1.3. Governing Equations</title><p>As temperature variation is low and the density is constant, microalgae suspension is considered to be incompressible fluid so Equation (1) reduces to</p><disp-formula id="scirp.72910-formula12"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x10.png"  xlink:type="simple"/></disp-formula><p>and Equation (2) becomes</p><disp-formula id="scirp.72910-formula13"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x11.png"  xlink:type="simple"/></disp-formula><p>where σ is the stress tensor and g is the gravity, σ can be expressed as</p><disp-formula id="scirp.72910-formula14"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x12.png"  xlink:type="simple"/></disp-formula><p>where, η = viscosity of the fluid; D(v) = rate of deformation. The viscosity η(t) in Equation (7) is determined by</p><disp-formula id="scirp.72910-formula15"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x13.png"  xlink:type="simple"/></disp-formula><p>The relative viscosity <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-1100557x14.png" xlink:type="simple"/></inline-formula> relating to the concentration is then used and determined by</p><disp-formula id="scirp.72910-formula16"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x15.png"  xlink:type="simple"/></disp-formula><p>where ε is the Einstein’s coefficient [<xref ref-type="bibr" rid="scirp.72910-ref19">19</xref>] . Based on the experiment conducted by Hon- Nami and Kunito [<xref ref-type="bibr" rid="scirp.72910-ref20">20</xref>] , the concentration function C(t) in Equation (9) is given by the logistic Equation (10)</p><disp-formula id="scirp.72910-formula17"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x16.png"  xlink:type="simple"/></disp-formula><p>where &#181; is the constant growth rate of microalgae cells; C<sub>0</sub> is the initial concentration of the suspension and A and B are constants.</p></sec><sec id="s3_1_4"><title>3.1.4. Boundary and Initial Conditions</title><p>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:</p><disp-formula id="scirp.72910-formula18"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x17.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.72910-formula19"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-1100557x18.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s3_2"><title>3.2. Mesh Design</title><sec id="s3_2_1"><title>3.2.1. Mesh Generation</title><p>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. <xref ref-type="table" rid="table3">Table 3</xref> shows the mesh elements, minimum and average quality for three PBR’s.</p><p>The numbers of vertex, edge and boundary elements are given in <xref ref-type="table" rid="table4">Table 4</xref>.</p><p><xref ref-type="fig" rid="fig3">Figure 3</xref>(a) shows the mesh design of the U-loop of tubular, square and hexagonal shape PBR’s and <xref ref-type="fig" rid="fig3">Figure 3</xref>(b) shows the vertical cross sectional views of the corresponding PBR’s.</p></sec><sec id="s3_2_2"><title>3.2.2. Grid Sensitivity Analysis</title><p>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. <xref ref-type="table" rid="table5">Table 5</xref> shows the comparison of elements between normal and coarse mesh.</p><p>From <xref ref-type="table" rid="table5">Table 5</xref> we observed that mesh quality and elements both are better in normal mesh.</p></sec></sec><sec id="s3_3"><title>3.3. Simulation Parameters</title><p>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 <xref ref-type="table" rid="table6">Table 6</xref>.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Total elements &amp; quality of PBR’s</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >PBR’s shape type</th><th align="center" valign="middle" >Total elements</th><th align="center" valign="middle" >Minimum quality</th><th align="center" valign="middle" >Average quality</th></tr></thead><tr><td align="center" valign="middle" >Tubular</td><td align="center" valign="middle" >153,356</td><td align="center" valign="middle" >0.03063</td><td align="center" valign="middle" >0.6759</td></tr><tr><td align="center" valign="middle" >Square</td><td align="center" valign="middle" >194,771</td><td align="center" valign="middle" >0.07511</td><td align="center" valign="middle" >0.6691</td></tr><tr><td align="center" valign="middle" >Hexagonal</td><td align="center" valign="middle" >153,122</td><td align="center" valign="middle" >0.09728</td><td align="center" valign="middle" >0.6467</td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Mesh parameters of the PBR’s</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >PBR’s shape type</th><th align="center" valign="middle" >Vertex elements</th><th align="center" valign="middle" >Edge elements</th><th align="center" valign="middle" >Boundary elements</th></tr></thead><tr><td align="center" valign="middle" >Tubular</td><td align="center" valign="middle" >12</td><td align="center" valign="middle" >1825</td><td align="center" valign="middle" >28,472</td></tr><tr><td align="center" valign="middle" >Square</td><td align="center" valign="middle" >16</td><td align="center" valign="middle" >2416</td><td align="center" valign="middle" >19,158</td></tr><tr><td align="center" valign="middle" >Hexagonal</td><td align="center" valign="middle" >24</td><td align="center" valign="middle" >5145</td><td align="center" valign="middle" >30,746</td></tr></tbody></table></table-wrap><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Comparison between normal and coarse mesh</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Mesh Type</th><th align="center" valign="middle" >Total elements</th><th align="center" valign="middle" >Minimum quality</th><th align="center" valign="middle" >Average quality</th><th align="center" valign="middle" >Mesh volume (m<sup>3</sup>)</th></tr></thead><tr><td align="center" valign="middle" >Normal</td><td align="center" valign="middle" >153,526</td><td align="center" valign="middle" >0.1506</td><td align="center" valign="middle" >0.6948</td><td align="center" valign="middle" >0.1551</td></tr><tr><td align="center" valign="middle" >Coarse</td><td align="center" valign="middle" >53,189</td><td align="center" valign="middle" >0.06733</td><td align="center" valign="middle" >0.5986</td><td align="center" valign="middle" >0.1519</td></tr></tbody></table></table-wrap><table-wrap id="table6" ><label><xref ref-type="table" rid="table6">Table 6</xref></label><caption><title> Parameters used for simulation</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Name</th><th align="center" valign="middle" >Value</th><th align="center" valign="middle" >Description</th></tr></thead><tr><td align="center" valign="middle" >g</td><td align="center" valign="middle" >9.8 m/s^2</td><td align="center" valign="middle" >Gravity acceleration</td></tr><tr><td align="center" valign="middle" >η<sub>0</sub></td><td align="center" valign="middle" >0.001 [Pa*s]</td><td align="center" valign="middle" >Water viscosity</td></tr><tr><td align="center" valign="middle" >C<sub>0</sub></td><td align="center" valign="middle" >0.55</td><td align="center" valign="middle" >Constant parameter</td></tr><tr><td align="center" valign="middle" >B</td><td align="center" valign="middle" >200</td><td align="center" valign="middle" >Constant parameter</td></tr><tr><td align="center" valign="middle" >A</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >Constant parameter</td></tr><tr><td align="center" valign="middle" >μ</td><td align="center" valign="middle" >0.063 [1/h]</td><td align="center" valign="middle" >Maximum growth rate</td></tr></tbody></table></table-wrap><fig-group id="fig3"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> (a) Mesh design for three PBR’s showing planes C1, C2, C3 (C1: Entrance of U-loop, C2: Middle of U-loop, C3: Exit of U-loop); (b) Cross section of normal mesh for circular, square &amp; hexagonal PBR’s respectively.</title></caption><fig id ="fig3_1"><label>(b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x19.png"/></fig><fig id ="fig3_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x20.png"/></fig></fig-group></sec></sec><sec id="s4"><title>4. Results and Discussion</title><p>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.</p><p>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 [<xref ref-type="bibr" rid="scirp.72910-ref14">14</xref>] . As we previously discussed about the grid independency, some results regarding velocity magnitude are shown in <xref ref-type="fig" rid="fig4">Figure 4</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref> for normal and coarse mesh in case of tubular PBR.</p><p>In <xref ref-type="fig" rid="fig4">Figure 4</xref>(a), for coarse mesh the velocity magnitude is 0.9168 m/s and in <xref ref-type="fig" rid="fig4">Figure 4</xref>(b), for coarse mesh it is 0.9316 m/s. The both cross sections are taken at the middle of U-loop. From <xref ref-type="fig" rid="fig5">Figure 5</xref>(a) the maximum velocity at straight portion is 0.75 for coarse mesh and from <xref ref-type="fig" rid="fig5">Figure 5</xref>(b) it is 0.8 for normal mesh.</p><sec id="s4_1"><title>4.1. Velocity Profile</title><p>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 <xref ref-type="fig" rid="fig6">Figure 6</xref>, the velocity profiles along XY-plane are shown of three PBR’s at 50 second respectively.</p><p>From <xref ref-type="table" rid="table7">Table 7</xref> 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.</p><p>In <xref ref-type="fig" rid="fig7">Figure 7</xref> the vertical cross sectional views for three PBR’s at the straight portion, entrance of U-loop, and middle of U-loop and outlet of U-loop are shown for 50 s which also exhibits the same phenomenon as it is found along XY-plane. The cross sectional views are taken in ZX and YZ planes. From the cross sectional views, it is clear</p><fig-group id="fig4"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Velocity at the middle of U-loop of tubular PBR for (a) coarse mesh (b) normal mesh.</title></caption><fig id ="fig4_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x21.png"/></fig></fig-group><fig-group id="fig5"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Velocity profile at straight portion of tubular PBR for (a) coarse mesh (b) normal mesh.</title></caption><fig id ="fig5_1"><label>(b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x22.png"/></fig><fig id ="fig5_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x23.png"/></fig></fig-group><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Velocity profile along XY-plane for (d) tubular (e) square (f) hexagonal ducts at 50 s</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x24.png"/></fig><table-wrap id="table7" ><label><xref ref-type="table" rid="table7">Table 7</xref></label><caption><title> Maximum velocity magnitude at 50 s</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >PBR’s shape type</th><th align="center" valign="middle" >Max. Velocity (m/s)</th></tr></thead><tr><td align="center" valign="middle" >Tubular</td><td align="center" valign="middle" >0.9287</td></tr><tr><td align="center" valign="middle" >Square</td><td align="center" valign="middle" >0.95</td></tr><tr><td align="center" valign="middle" >Hexagonal</td><td align="center" valign="middle" >0.95</td></tr></tbody></table></table-wrap><fig-group id="fig7"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Velocity profile at arc length = 5 m; (a), (b), (c) at 10 m; (d), (e), (f) at 10.4 m; at middle portion of U-loop respectively at 50s [(a), (d), (b), (e) &amp; (c), (f) have same legend respectively].</title></caption><fig id ="fig7_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x25.png"/></fig><fig id ="fig7_2"><label> (c)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x26.png"/></fig><fig id ="fig7_3"><label>(d)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x27.png"/></fig><fig id ="fig7_4"><label> (e)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x28.png"/></fig></fig-group><p>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.</p><p>The velocity distribution graphs for the tubular, the square and the hexagonal shape PBR are shown in <xref ref-type="fig" rid="fig8">Figure 8</xref>. 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.</p><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> (a), (b), (c) at arc length = 5 m; (d), (e), (f) at 10 m (entrance of U-loop); (g), (h), (i) at 10.2 m (middle portion of U-loop); (j), (k), (l) at 10.4 m (outlet of U-loop) at 50 s. [In each row the graphs are for tubular, square and hexagonal PBR respectively from left to right]</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x29.png"/></fig></sec><sec id="s4_2"><title>4.2. Shear Stress</title><p>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 <xref ref-type="fig" rid="fig9">Figure 9</xref>.</p><p>By comparing three different geometrical shape of PBR, less shear stress is observed for the tubular PBR. <xref ref-type="table" rid="table8">Table 8</xref> shows the maximum shear stress for different cross sections of three PBR’s.</p><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> Shear stress distribution at three different cross sections: at the inlet of U-loop (Pc1: tubular, Pc2: square, Pc3: Hexagonal), middle of U-loop (Pc4: tubular, Pc5: square, Pc6: Hexagonal), outlet of U-loop (Pc7: tubular, Pc8: square, Pc9: Hexagonal)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x30.png"/></fig><table-wrap id="table8" ><label><xref ref-type="table" rid="table8">Table 8</xref></label><caption><title> Minimum shear stress for three PBR’s</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >PBR’s shape type</th><th align="center" valign="middle" >Inlet of U-loop</th><th align="center" valign="middle" >Middle of U-loop</th><th align="center" valign="middle" >Outlet of U-loop</th></tr></thead><tr><td align="center" valign="middle" >Tubular</td><td align="center" valign="middle" >66.645</td><td align="center" valign="middle" >89.281</td><td align="center" valign="middle" >66.645</td></tr><tr><td align="center" valign="middle" >Square</td><td align="center" valign="middle" >77.2</td><td align="center" valign="middle" >103</td><td align="center" valign="middle" >77.2</td></tr><tr><td align="center" valign="middle" >Hexagonal</td><td align="center" valign="middle" >75.1</td><td align="center" valign="middle" >106</td><td align="center" valign="middle" >75.1</td></tr></tbody></table></table-wrap></sec><sec id="s4_3"><title>4.3. Pressure Distribution</title><p>The pressure profiles are uniform from inlet to outlet for all three types of geometry. From <xref ref-type="fig" rid="fig1">Figure 1</xref>0, we can see pressure along the ducts in XY-plane of tubular PBR. <xref ref-type="table" rid="table9">Table 9</xref> 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 <xref ref-type="fig" rid="fig1">Figure 1</xref>1.</p></sec></sec><sec id="s5"><title>5. Conclusion</title><p>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</p><fig id="fig10"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> Pressure distribution along the ducts (pp1: Tubular, PP2: Square, PP3: Hexagonal)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x31.png"/></fig><table-wrap id="table9" ><label><xref ref-type="table" rid="table9">Table 9</xref></label><caption><title> Maximum pressure along the PBR’s</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Photobioreactor</th><th align="center" valign="middle" >Maximum pressure</th></tr></thead><tr><td align="center" valign="middle" >Tubular</td><td align="center" valign="middle" >83</td></tr><tr><td align="center" valign="middle" >Square</td><td align="center" valign="middle" >87.641</td></tr><tr><td align="center" valign="middle" >Hexagonal</td><td align="center" valign="middle" >84.8</td></tr></tbody></table></table-wrap><fig id="fig11"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>1</label><caption><title> Pressure profile along the PBR’s arc length (P1: Tubular, P2: Square, P3: Hexagonal)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-1100557x32.png"/></fig><p>of the PBR’s. So from our analysis it can be said that tubular PBR is the best choice for microalgae culture. Though pressure profile is always uniformly decreasing from inlet to outlet for all PBR’s, a little fluctuation is found in U-loop portion as for haphazard distribution of fluid. But the pressure is always less in tubular PBR. So the simulation result of fluid properties velocity, shear stress and pressure indicates that the tubular one shows better agreement for the culture of microalgae.</p></sec><sec id="s6"><title>Acknowledgements</title><p>The authors are gratefully acknowledged for the technical supports to the Centre of Excellence in Mathematics, Department of Mathematics, Mahidol University, Bangkok- 10400, Thailand, and the Simulation Lab, Department of Mathematics, Chittagong University of Engineering &amp; Technology, Chittagong, Bangladesh.</p></sec><sec id="s7"><title>Cite this paper</title><p>Shahriar, M., Monir, M.I. and Deb, U.K. (2016) Comparative Analysis of Hydrodynamics Behavior of Microalgae Suspension Flow in Circular, Square and Hexagonal Shape Photo Bioreactors. American Journal of Computational Mathematics, 6, 320-335. http://dx.doi.org/10.4236/ajcm.2016.64033</p></sec></body><back><ref-list><title>References</title><ref id="scirp.72910-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Brennan, L. and Owende, P. (2010) Biofuels from Microalgae—A Review of Technologies for Production, Processing, and Extractions of Biofuels and Co-Products. Renewable and Sustainable Energy Reviews, 14, 557-577. https://doi.org/10.1016/j.rser.2009.10.009</mixed-citation></ref><ref id="scirp.72910-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Lam, M.K. and Lee, K.T. (2012) Microalgae Biofuels: A Critical Review of Issues, Problems and the Way Forward. Biotechnology Advances, 30, 673-690. https://doi.org/10.1016/j.biotechadv.2011.11.008</mixed-citation></ref><ref id="scirp.72910-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Gavrilescu, M. and Chisti, Y. (2005) Biotechnology—A Sustainable Alternative for Chemical Industry. Biotechnology Advances, 23, 471-499. https://doi.org/10.1016/j.biotechadv.2005.03.004</mixed-citation></ref><ref id="scirp.72910-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Chisti, Y. (2007) Biodiesel from Microalgae. Biotechnology Advances, 25, 294-306.https://doi.org/10.1016/j.biotechadv.2007.02.001</mixed-citation></ref><ref id="scirp.72910-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Ramírez-Duque, J.L. and Ramos-Lucumi, M.A. (2011) Hydrodynamic Computational Evaluation in Solar Tubular Photobioreactors Bends with Different Cross Sections. CT&amp;F: Ciencia, Tecnología y Futuro, 4, 59-72.</mixed-citation></ref><ref id="scirp.72910-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Schenk, P.M., Thomas-Hall, S.R., Stephens, E., Marx, U.C., Mussgnug, J.H., Posten, C., Kruse, O. and Hankamer, B. (2008) Second Generation Biofuels: High-Efficiency Microalgae for Biodiesel Production. BioEnergy Research, 1, 20-43. https://doi.org/10.1007/s12155-008-9008-8</mixed-citation></ref><ref id="scirp.72910-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Nigam, P.S. and Singh, A. (2010) Production of Liquid Biofuels from Renewable Resources. Progress in Energy and Combustion Science, 37, 52-68. https://doi.org/10.1016/j.pecs.2010.01.003</mixed-citation></ref><ref id="scirp.72910-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Blanken, W., Postma, P.R., de Winter, L., Wijffels, R.H. and Janssen, M. (2016) Predicting Microalgae Growth. Algal Research, 14, 28-38.</mixed-citation></ref><ref id="scirp.72910-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Bitog, J.P., Lee, I.-B., Lee, C.-G., Kim, K.-S., Hwang, H.-S., Hong, S.-W., Seo, I.-H., Kwon, K.-S. and Mustafa, E. (2011) Application of Computational Fluid Dynamics for Modelling and Designing Photo Bioreactors for Microalgae Production: A Review. Computers and Electronics in Agriculture, 76, 131-147. https://doi.org/10.1016/j.compag.2011.01.015</mixed-citation></ref><ref id="scirp.72910-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Rengel, A. (2008) Promising Technologies for Biodiesel Production from Algae Growth Systems. Proceedings of the 8th European IFSA Symposium, Clermont-Ferrand, 6-10 July 2008, 683-692.</mixed-citation></ref><ref id="scirp.72910-ref11"><label>11</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Tredici</surname><given-names> M.R. </given-names></name>,<etal>et al</etal>. (<year>2010</year>)<article-title>Photobiology of Microalgae Mass Cultures: Understanding the Tools for the Next Green Revolution</article-title><source> Biofuels</source><volume> 1</volume>,<fpage> 143</fpage>-<lpage>162</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.72910-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Kunjapur, A.M. and Eldridge, R.B. (2010) Photobioreactor Design for Commercial Biofuel Production from Microalgae. Industrial &amp; Engineering Chemistry Research, 49, 3516-3526. https://doi.org/10.1021/ie901459u</mixed-citation></ref><ref id="scirp.72910-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Ugwu, C.U., Aoyagi, H. and Uchiyama, H. (2008) Photobioreactors for Mass Cultivation of Algae. Bioresource Technology, 99, 4021-4028. https://doi.org/10.1016/j.biortech.2007.01.046</mixed-citation></ref><ref id="scirp.72910-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Deb, U.K., Chayantrakon, K. and Lenbury, Y. (2012) Comparison of Single-Phase and Two-Phase Flow Dynamics in the HLTP for Microalgae Culture. International Journal of Mathematics and Computers in Simulation, 6, 496-503.</mixed-citation></ref><ref id="scirp.72910-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Pruvost, J., Pottier, L. and Legrand, J. (2006) Numerical Investigation of Hydrodynamic and Mixing Conditions in a Torus Photobioreactor. Chemical Engineering Science, 61, 4476-4489. https://doi.org/10.1016/j.ces.2006.02.027</mixed-citation></ref><ref id="scirp.72910-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Molina, E., Fernandez, J., Acien, F.G. and Chisti, Y. (2001) Tubular Photobioreactor Design for Algal Cultures. Journal of Biotechnology, 92, 113-131. https://doi.org/10.1016/S0168-1656(01)00353-4</mixed-citation></ref><ref id="scirp.72910-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Khanam, I.A. and Deb, U.K. (2016) Calculation of the Average Irradiance and the Microalgae Growth for a Year at CUET, Bangladesh. American Journal of Computational Mathematics, 6, 237-244. https://doi.org/10.4236/ajcm.2016.63024</mixed-citation></ref><ref id="scirp.72910-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Deb, U.K., Kittisak, C., Lenbury, Y. and Benchawan, W. (2012) Numerical Simulation of 2-p Laminar Flow for CO2 and Microalgae Suspension in the HLTP. In: Latest Advances in Systems Science and Computational Intelligence, Wseas LLC, 53-58.</mixed-citation></ref><ref id="scirp.72910-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Einstein, A. (1906) Eine neue Bestimmung der Moleküldimensionen. Annalen der Physik, 324, 289-306. https://doi.org/10.1002/andp.19063240204</mixed-citation></ref><ref id="scirp.72910-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Hon-Nami, K. and Kunito, S. (1998) Microalgae Cultivation in a Tubular Bioreactor and Utilization of Their Cells. Chinese Journal of Oceanology and Limnology, 16, 75-83. https://doi.org/10.1007/BF02849084</mixed-citation></ref></ref-list></back></article>