A Dynamic Optimization Model for Designing Open-Channel Raceway Ponds for Batch Production of Algal Biomass

This work focuses on designing the optimum raceway pond by considering the effects of sunlight availability, temperature fluctuations, and harvest time on algae growth, and introduces a dynamic programing model to do so. Culture properties such as biomass productivity, growth rate, and concentration, and physical properties, such as average velocity, pond temperature, and rate of evaporation, were estimated daily depending on the dynamic behavior of solar zenith angle, diurnal pattern of solar irradiance, and temperature fluctuations at the location. Case studies consider two algae species (Phaeodactylum. tricornutum and Isochrysis. galbana) and four locations (Tulsa, USA; Hyderabad, India; Cape Town, South Africa; and Rio de Janeiro, Brazil). They investigate the influences of the type of algae strain and geographical location on algae biomass production costs. From our case studies, the combination of I. galbana species grown in Hyderabad, India, with a raceway pond geometry of 30 cm channel depth, about a meter channel width, and 300 m in length, and a harvest interval of every six days yielded the minimum algal biomass production costs. The results of the sensitivity analysis reveal that smaller channel depths and longer ponds (within the ranges considered) are recommended to minimize the net present cost of algae biomass production.


Introduction
Large-scale production of micro-algal biomass is achieved by open pond cultivation technologies [1].Open pond cultivation technologies have been known since 1950s.They can be classified as circular ponds, unmixed ponds, inclined ponds, and raceway ponds, depending on their design and operation [2].The most preferable among them for algae cultivation are raceway ponds [3].They are single or multiple closed-loop oval-shaped recirculation channels with paddle wheel that enable circulation and prevent sedimentation.The principal advantages of raceway pond cultivation technology are small capital investment, use of free solar energy [2], easy maintenance, and lower energy requirements [4].The idea of cultivating algae for the purpose of biofuels production was first investigated by Oswald W. J. in 1970s.Later, in 1978, the U.S. Department of Energy's Aquatic Species Program attempted to identify the algal species and conditions that would generate a cost-effective algae biofuels production [5].This eighteen-year program explored open ponds and yielded a collection of high oil producing species.
There is a significant body of literature that studies the cultural properties of algal biomass for a given outdoor raceway pond or focuses on designing new outdoor pond configurations.In 1979, a theoretical framework was developed by Goldman for continuous culturing of algal biomass to predict the biomass productivity as a function of sunlight and temperature in open raceways [6].The model also took into account the effects of culture depth and nutrients on light intensity.The upper limit on biomass yield was determined to be 30-40 dry¨g¨m ´2¨day ´1.Hill and Lincoln developed a mathematical model for continuous flow mass culture of microalgae in multi rectangular open channels.The model predicted algal production, dissolved oxygen, pH, and concentrations of nutrients [7].Grobbelaar et al. developed a deterministic mathematical model to describe the production of the green microalgae Scenedesmus.obliquus (S. obliquus) and Coelastrum.sphaericum (C.sphaericum) in outdoor raceways [8].Assuming that CO 2 and other nutrients were supplied in excess during cultivation, the dynamic model incorporated sixteen months of irradiance and temperature measurements to evaluate the productivity.The productivity varied between 1.7 and 16.92 g¨m ´2¨day ´1.Sukenik et al. developed a deterministic simulation model to estimate productivity of the marine algae Isochrysis.Galbana (I.galbana) in an outdoor raceway pond [9].The effects of pond depth and chlorophyll concentration on production rate were evaluated in various seasons.The model predicted an annual average productivity of 9.7 g of carbon uptake per square meter per day (gC¨m ´2¨day ´1).A new regulatory model was developed by Geider et al. [10] to describe the acclimation of phytoplankton growth rate to irradiance, temperature and nutrient availability.The model used Poisson function to model photosynthesis, Arrhenius function to model temperature and Michaelis-Menten kinetics to model nutrient uptake.James et al. [11] simulated hydrodynamics coupled with growth kinetics of Phaeodactylum.tricornutum (P.tricornutum) in open-channel raceways using modified versions of Environmental Fluid Dynamics Codes (EFDC) [12] and U.S. Army Corp of Engineers' water quality code (CE-QUAL) [13].System parameters such as flow velocity and raceway depth were varied to improve growth rate.Their results show that the depth of the pond should vary in response to the atmospheric temperature to maximize algae yield.They also reported that increasing the flow rate above 6.25 L/s does not further improve the algae growth rate.
Ketheesan et al. developed a new airlift-driven raceway configuration, that is paddle wheel free, for energy-efficient algal cultivation [14].Harvesting was performed after 14 days when culture reached the end of linear growth phase.Sompech et al. used computational fluid dynamics to model the energy consumption, the extent and location of dead zones in various geometric configurations of raceway ponds at various flow velocities [15].They concluded that the standard configuration with baffles on each end are relatively inexpensive and prevented the formation of dead zones.Chiaramonti et al. designed an innovative raceway pond and compared its performance with conventional raceway pond [16].The innovative raceway pond was shallower and paddle wheel was replaced by centrifugal pump.They investigated the energy consumption and concluded that the new design saved 50% energy over an area of 500 m 2 .Talent et al. [17] developed a protocol for harvesting microalgae massive cultures from open ponds.According to them, a minimum of 10% of the pond volume per day should be harvested and replaced with new feedstock water.
Many studies focus on predicting the culture properties such as productivity, growth rate, and biomass concentration depending on the sunlight and temperature fluctuations at a given facility under known pond configurations.Other studies include varying the raceway pond geometry to understand the physical parameters such as average velocity, evaporation rate, energy consumption and their influence on operating costs.The optimum pond geometry depends on the growth characteristics such as biomass concentration, growth rate, and productivity of the algae species which are affected by sunlight availability and temperature fluctuations at a location.There is a gap in literature of a model that relates all three: pond geometry, location and type of algae species.The key contribution of this paper is the development of a dynamic programming formulation that considers algae strain characteristics, diurnal pattern of solar irradiance, temperature fluctuations, and harvesting intervals to yield the optimum design configuration of an outdoor raceway pond.The model defines optimum design as the one that minimizes total production cost for a plant life of 10 years, including both capital and operating costs.Our model also captures the dynamics of biomass concentration, average light inside the pond, growth rate, production, and radiation entering and leaving the pond.

Problem Statement
Given an algae species, a geographical location to grow algae, and algal biomass demand, the goal is to obtain the optimum raceway pond geometry considering the effect of physical parameters such as temperature fluctuations in the growth media, diurnal pattern of sunlight, harvest schedule, light limitation, and evaporation losses.Here, the optimum design is the one that yields the minimum net present cost for biomass production with a plant life of 10 years.The net present cost is a combination of both capital costs and operating costs.Capital costs include site preparation, grading, compaction, pond levees, paddle wheel, CO 2 diffusers, harvesting, flocculation, centrifugation, water and nutrient supply system, waste treatment system, buildings and structures, electrical supply and distribution, instrumentation and machinery, land, and engineering and contingency costs.Operating costs include electricity, water, nutrients, flocculants, waste disposal, maintenance and repairs, operating supplies, depreciation, taxes and insurance, financing, labor and overhead costs, and general expenses.

Algal Production Model for Raceway Ponds
This section explains the changes in culture properties of the algae species and their dependencies on the geometry of the raceway pond and weather conditions at the location.
Algae growth: Figure 1 shows the aerial view of the raceway pond and its geometry.The pond is considered to be perfectly mixed and of uniform concentration.However, the concentration changes with time due to algae growth.Equation (1) explains the dynamic changes occurring in biomass concentration with the production of dry algae in a given pond volume V (assuming constant density).In the pond, with time as algae cells begin to multiply, biomass concentration rises.In our model, we considered seeding the pond just before sunrise, i.e., biomass concentration, has been initialized to certain value just before sunrise after harvest.
here, V (m 3 ) is the volume of the raceway pond and BC (g¨m ´3) and X DA (g) represent the biomass concentration and the amount of dry algae in the pond at any given time.
The amount of sunlight reaching the algae cells inside the pond varies depending on the accumulation of biomass, the culture's light absorption coefficient, and the depth of the pond.During the cultivation of algae in outdoor raceway ponds, light intensity is attenuated because of gradual increase in the concentration of biomass [18] and due to uneven absorption of light by algae cells.The algae cells closer to light source receive higher irradiance than cells further away.Grima et al. proposed an empirical relationship between light attenuation and biomass concentration, and verified its accuracy for several types of freshwater microalga, P. tricornutum [19].The model was later validated for different algae species and photobioreactor orientations in literature [20].The model first, estimates the average solar irradiance, I avg (t) (µE¨m ´2¨s ´1), experienced by a single cell moving inside the culture medium (via Equation ( 2)) based on the position of the sun relative to the pond and Beer-Lambert relationship [21].
I avg ptq " I o ptq Ka ¨BCptq ¨CD eq ptq " 1 ´e´pKa¨BCptq¨CD eq ptqq ı In Equation ( 2), Ka is light absorption coefficient and varies with algae species.I o (t) (µE¨m ´2¨s ´1) is the incident light on the surface which is a function of daily maximum and minimum temperatures as explained by Hargreaves [22].It is shown via Equation (24) of the Supplementary Materials provided with this manuscript.CD eq (t) is the length of light path from the surface to any point inside the pond.The length of light path depends on the depth of the pond, CD (m), and the position of the sun relative to the earth's surface called solar zenith angle as shown in Equation (3).The greater the pond depth, the longer is the length of the light path.The lower the sun's position from the Earth's surface, during sunrise and sunset, the longer is the length of the light path in the pond.However, when the sun's position is high in the sky, the length of the light path is shorter.
Growth rate of algae cells depend on their light extinction coefficient, I k (µE¨m ´2¨s ´1), and the availability of light inside the pond, I avg (t) as shown in Equation ( 4).The equation does not take into account photoinhibition.Although it is well documented, it has often been disregarded.Studies suggest that growth models that express µ in terms of the average irradiance raised to some power greater than unity better fit experimental observations [21,23].Equation ( 4) is a light-limited growth equation which was developed and tested by Grima, Camacho, Pérez, Sevilla, Fernández, and Gómez [19].Based on the study, material balance performed in a chemostat where specific growth rate is expressed as a function of irradiance, it was observed that the algae growth rate predictions of this model exhibit a good agreement with the measured ones.It models the effect of light attenuation on the observed algae growth rate in the cultivation unit.It defines a hyperbolic relationship between the specific growth rate of algae and average solar irradiance inside the culture.
µptq " µ max ptq here, n is an exponent that describes the abruptness of the transition from weakly-illuminated to strongly-illuminated regions and is obtained from non-linear regression analysis on light intensity as a function of biomass concentration [19].
The maximum specific growth rates, µ max , (day ´1) of algae species depend strongly on the temperature of the pond and the species itself [24].Several models describe the influence of temperature on growth rate of algae, we have chosen the model presented in [25] where eight species of marine phytoplankton were grown at different temperatures ranging from 10 to 25 ˝C.They experimentally modeled the general response of growth rate to temperature as µ max " α.exp pβ.Tq where α and β are species dependent constants and T is temperature of medium.In our work, we have used their data to model the growth response of the selected species.We used the same relationship in our model via Equation (5).
where T pond is temperature of the pond.The species specific constants α and β are taken from our previous findings [26].The experimental data modeled the response of growth rates of eight species of marine phytoplankton to changes in temperature under continuous light [27].Growth rate, biomass concentration, and algae cells in the pond are interrelated.The biomass concentration and the specific growth rate determine the rate of change of algal biomass in the pond.This relation is shown via Equation (6).As the biomass concentration and growth rate increase from sunrise to further during the day, the algae cells multiply.However, as concentration further increases, the growth rate begins to decline due to mutual shading, thus effecting the production of algae cells.
Temperature of the pond is essential for proper design and optimal operation of algal ponds.Based on the pond temperature, algae could either exhibit optimal growth characteristics or could perish above a certain threshold temperature.In this work, pond temperature is calculated using the temperature model that is "universally applicable" for shallow ponds containing opaque water body [28].According to Béchet, Q., et al., the model was tested against experimental data that has been collected for a period of one year from a wastewater treatment high-rate algal pond and was validated with a "high degree of accuracy" [28].The model is shown explicitly via Equations ( 7) through (17).Here, the pond temperature, T pond (t) (K), is calculated via an overall energy balance around the pond (Equation ( 17)).According to [28], conduction was of low importance in the total heat balance for shallow algae ponds.Thus conductive heat flux is neglected for this model since we assume no heat transfer between the pond bottom and the soil beneath.We also ignored the changes in pond temperature caused due to rainfall, and hence, removed the rain heat flux term from the model.Therefore, temperature of the outdoor pond is determined by six heat fluxes: (1) Heat flux due to pond radiation; (2) heat flux due to solar radiation; (3) heat flux due to air radiation; (4) Evaporation; (5) Convection; and (6) Inflow heat flux.
(1) Heat flux due to pond radiation: Radiation from the pond's surface to the atmosphere Q p (t) is estimated by Stefan-Boltzmann's fourth power law as shown in Equation (7) [28].Since heat transfer takes place from pond surface to the atmosphere, this heat flux is negative.The higher the radiation from the pond, the higher heat transfer from pond to the atmosphere becomes, which lead to lower pond temperature.Q p ptq " ´εw σT 4  pond ptqSA where ε w = 0.97 represents the emissivity of water [29], σ = 5.67 ˆ10 ´8 represents Stefan-Boltzmann constant (W¨m ´2¨K ´4) and SA is the surface area of the pond (m 2 ).(2) Heat flux due to solar radiation: Radiation received by the pond from the sun Q s (t) is calculated from incident irradiance on the surface in Equation ( 8) [28].The higher the radiation from the sun to the pond, the higher would be the temperature of the pond.
where f a is the fraction of sunlight converted by algae into chemical energy during photosynthesis.
It is assumed to be constant and equal to 2.5%.(3) Heat flux due to air radiation: Radiation from air, Q a (t), is given by Stefan-Boltzmann's fourth power law as shown in Equation ( 9) [28].T surr (t) (K) is calculated by the method outlined by Woodhead [30].The higher the radiation from the air to the pond, the higher would be the pond temperature.
where L w = 2.45 ˆ10 6 J¨kg ´1 represents the latent heat of water.Theoretically, applying Buckingham theorem to evaporation at the pond surface yields three dimensionless numbers namely, the Sherwood number Sh L , the Schmidt number Sch L , and the Reynolds number Re L .The rate of evaporation is shown to be dependent on the three dimensionless numbers [31].For mass transfer in a horizontal surface, the three dimensionless groups are correlated as follows: Sh L " 0.035 pRe L q 0.8 pSch L q 1{3 (11) where Sherwood number: Schmidt number: Reynolds number: where K is the mass transfer coefficient (m¨s ´1) determined via Equation ( 11), D h is the hydraulic diameter of the pond (m), D w,a = 2.4 ˆ10 5 m 2 ¨s´1 represents the mass diffusion coefficient of water vapor in air, ϑ a = 1.5 ˆ10 5 m 2 ¨s´1 denotes the air kinematic viscosity, and WV(t) represents the wind velocity (m¨s ´1).Equation ( 11) has been derived with an assumption that wind is constant with respect to height.The rate of evaporation is calculated by Equation ( 12).
meptq " Kptq where MW H 2 O is the molecular weight of water (kg¨mole ´1), R = 8.314 Pa¨m 3 mol ´1¨K ´1 is the ideal gas constant, and RH(t) represents the relative humidity of the air over the pond surface.Relative humidity is important with regard to both evaporation losses and pond cooling.When humidity is low, high rates of evaporation occur, particularly during windy periods.When humidity is high and there is little or no wind, and when sunlight is abundant, the water in shallow cultures may heat up.
In Equation (12), P a and P w are saturated vapor pressures (Pa) at T surr and T pond , respectively, and can be determined using Equation ( 13) [32].
P i ptq " 3385.5e´r´8.0929`0.97608pTi ptq`42.607´273.15q 0.5 s ( (5) Convection: Convective heat transfer, Q cv (t), is expressed similar to evaporative heat transfer as heat transfer and mass transfer obey the same laws.Applying Buckingham theorem to convection at the pond surface yields three dimensionless numbers namely, the Nusselt number, the Prandtl number, and the Reynolds number.Convection coefficient is calculated using the correlation between the three dimensionless numbers (Equation ( 14)) [31].
(6) Inflow heat flux: The heat transfer associated with the water inflow Q i (t) can be expressed as a steady flow thermal energy equation shown in Equation (16).Water is added to compensate for evaporation from the pond surface, and it is assumed to be at the surrounding temperature T surr (t) when added to the pond.Q i ptq " pmeptq ¨SAq Cp ´Tsurr ptq ´Tpond ptq ¯ (16) where the first set of parenthesis indicates mass flow rate of incoming water.It is the product of rate of evaporation and surface area of the pond.Cp represents the specific heat capacity of the pond water (J¨g ´1¨K ´1).The overall energy balance for the pond is given in Equation ( 17). .

mptqCp
BT pond ptq Bt " Q p ptq `Qs ptq `Qa ptq `Qev ptq `Qcv ptq `Qi ptq In Equation ( 17), Q p (t) is the rate of heat flow from the pond surface (W), Q s (t) is the rate of heat flow from sun to the pond (W), Q a (t) is the rate of heat flow from air to the pond (W), Q ev (t) is the rate of heat flow by evaporation (W), Q cv (t) is the rate of heat flow by convection at the surface of the pond (W), and Q i (t) is the rate of heat flow associated with water inflow (W).The mass flow rate, .m(t), can be calculated via Equation ( 18) using the concepts of fluid dynamics where mass flow rate is density times volumetric flow rate.
. mptq " ρ ¨Uavg ptq ¨CW ¨CD (18) where ρ is the density of water, U avg (t) is the average daily velocity of the pond, CW is the channel width, and CD is the channel depth.
The proposed model assumes no cloud coverage, and hence, it ignores the effects of cloud cover on biomass growth rate and productivity.The model also assumes the usage of atmospheric carbon dioxide at no additional cost.Nutrients, such as nitrogen and phosphorus, are assumed to be supplied, as needed, in their appropriate forms for algal growth in the growth medium.The pond is assumed to be completely mixed, and all physical properties are assumed to be that of water under standard conditions.

Dynamic Model
The objective of the optimization model is to find the raceway pond geometry that best suits the chosen location to grow the desired algae species.The objective function, to be minimized, considers the overall production costs involved.A complete listing of the subscripts, parameters, variables, and constraints of the program is given below.

Decision Variables
Depending on their function, variables used in designing the raceway pond can be categorized into three types: culture properties, physical properties and design properties.
Physical properties are defined as: (1) average velocity (m¨s ´1) of the pond (U avg (t)); (2) average irradiance (µE¨m ´2¨s ´1) inside the pond (I avg (t)); ( 3) length (m) of light path from the surface to any point in the pond (CD eq (t)); ( 4

The Model
The objective is to fulfill the biomass demand at a minimum net present cost, Z, for a raceway pond with a plant-life of 10 years.The objective function and constraints are as follows: where Z Water " WtCost l ¨Indwater (23) subject to: D algaebiomass,y ě demand @y (25) D algaebiomass,y " D DA,y `DWIA,y @y (26) X WW ptq " mptq " X DA ptq `XWIA ptq `XWW ptq q " 2CW (32) SA " V " SA ¨CD PrV " D DA,y V ¨360 PrA " PrV ¨V SA (37) PPptq " . mptq ¨g ¨hT ptq PWptq " 2 pmeptq ¨SAq ¨µ ¨length ¨p2CD `CWq 2 ¨Uavg ptq EPP " EPW "  20) assumes that the capital cost changes linearly with the pond surface area.Equations ( 21) and ( 22) calculate the cost of electricity required for mixing using paddle wheel and cost of electricity required for pumping additional water that is lost during evaporation.Equation ( 23) estimates the cost of water.Equation ( 24) calculates the total production cost.The estimation of cost coefficients for capital costs (TCI l ) and operating costs (TPI l ) are detailed in Section 6.1.
The first constraint, Equation ( 25), makes sure that the specified annual algae biomass demand is met.The second constraint, Equation ( 26), shows that the total annual production of wet algae biomass is the sum of annual production of dry algae biomass and water in algae.Equations ( 27) and (28) show that sum of water in algae and dry algae produced at different times in a year is equal to their annual productions.Equation ( 29) calculates the hourly production of WIA (water present in algae), given the percentage of dry algae in wet algae biomass for species s.The next constraint, Equation (30), shows the dependence of excess waste water on biomass concentration.If the biomass concentration is high, waste water production is low and vice versa.Also, the higher the biomass concentration, the higher is the dry algae present inside the pond.Total mass (DA, WIA, and WW) in the pond at any time is computed via Equation (31).
Equations ( 32) through (35) represent pond related geometric constraints.Ponds should have minimum number of bends to keep the energy required for recirculation to a minimum [33].Hence, this work focuses on the design of single channel raceway pond where the relationship between pond width and channel width is given in Equation (32).Pond length is the sum of channel length and pond width as shown in Equation (33).Surface area occupied by the pond is computed from Equation (34).To avoid the flow disturbance caused by the bends of the raceway pond, the ratio of pond length to width should be 10 or greater [33], and this requirement is enforced via Equation (47).Equation (35) calculates the volume of the raceway pond.
Equation (36) calculates the average daily volumetric productivity of microalgae on dry weight basis.Areal productivity is calculated as shown in Equation (37).Equations ( 45) and ( 46) are the bound constraints on pond depth and length.Pond depth should always be above an acceptable limit which may otherwise lead to pond evaporation and biomass solidifying issues.In our work, sensitivity analysis was performed on the lower bound of depth to see its impact on the dynamic changes of biomass concentration.Pond length is constrained to keep the head loss due to friction lower.The acceptable bounds on pond depth and length are taken from literature [34].Equation (48) specifies an upper bound on the average daily areal productivity of microalgae to be 60 g¨m ´2¨day ´1 [35].The cost coefficients in the objective function, TCI l and TPC l , were taken assuming the average daily areal productivity of 60 g¨m ´2¨day ´1.
Equations ( 38) (Manning's equation) and ( 39) calculate frictional and kinetic head losses [36].Total head loss is calculated in Equation ( 40) [36].Power required for maintaining flow in the raceway, i.e., power of paddlewheel is given by Equation ( 41) [36].Power required for pumping the inflow water that is added to compensate for annual evaporation losses is given by Equation (42).Annual energy requirements of paddle wheel and pumping the inflow water are calculated in Equations ( 43) and (44).
Mixing is required to expose all cells to light, distribute nutrients, and prevent settlement.The extent of thermal stratification in water bodies can be an indication of the degree of mixing.Velocities above 0.10 m¨s ´1 have been shown to prevent thermal stratification and sedimentation.Improper mixing may also result in poor CO 2 mass transfer rates causing low biomass productivity [37].Channel flow velocity is typically between 0.15 and 0.30 m¨s ´1, and Equation (49) enforces these bounds [38].Equation (50) constrains biomass concentration to be below 10,000 g¨m ´3 in order to inhibit light attenuation due to dense cultures and to maintain the Newtonian behavior of the growth medium while keeping the fluid properties similar to that of water [39].Equations ( 51) and (52) capture the dynamic changes in biomass concentration that occurs due to the production of algae biomass.These equations are analogues to the differential equations-Equations ( 1) and ( 6).

Solution Approach
The optimization model outlined in Sections 3 and 4 captures the dynamics of physical parameters such as solar irradiance and surrounding temperature and their impact on cultural properties such as biomass concentration and growth rate, and hence, on the outdoor open-channel raceway pond geometry.In order to capture the dynamics of the system with enough granularity, we discretized time into equal buckets of an hour, and assumed that the properties of algae biomass within each time bucket are constant.For an annual operational period (assuming operation of 360 days), this discretization yields 8640 time buckets.The resulting non-convex nonlinear program has 259,945 constraints and 268,587 variables.This nonlinear program is not tractable with global solvers, and it requires a feasible initial solution for yielding a local optimum solution with local solvers.
The following approach is developed to obtain approximate solutions of the problem quickly.Twelve representative days, each of which reflects the average hourly changes in one month of the year, are generated.It is, then, assumed that this representative day is repeated for 30 days, resulting in the annual production of 360 days.Considering only algal growth from sunrise to sunset of each day, the approximation reduces the number of equations and variables to 4219 and 4360, respectively.The details of this approximation is outlined in the next section.

Approach to Convert Annual Production to Hourly Production
One representative day is used to approximate the algal growth within a month, resulting in 12 representative days to simulate algal growth in a year.The first day of each calendar month is used as the marker for each representative day in the model.The harvest period (harvP), which is defined in terms of number of days algal biomass is allowed to accumulate in the pond before harvest, determines the corresponding harvest days for each month.For example, if the harvest period is equal to six days, for the month of January, 1 January is the marker for the representative day, and 6 January before sunset would be the harvest day.The data for daily maximum temperatures (T max l,d ), minimum temperatures (T min l,d ), hourly wind velocity (WV l,d (t)), and relative humidity (RH l,d (t)) for a year are collected from Wolfram Mathematica 8 Mathematica Weather data database [40].Monthly minimum and maximum of temperatures are calculated, and assigned to their respective representative day of the month.The averages for the wind velocity and relative humidity is calculated for each hour of the day within a month, and they are assigned to their corresponding representative day.Zenith angle is the position of the sun with respect to earth's surface.It not only helps in determining the incident irradiance at a location, I o (t), but also the sunrise and sunset times of the location.Hourly weighted averages of zenith angle were computed for all days in a month, and their monthly means were specified to their representative day of each month.By doing so, we have successfully created a day that represents the average changes in the relevant parameters within a day for each month.
With the hourly discretization of time and the representative day parameters, we can now calculate the algal growth rates, biomass concentrations, and the total algal accumulation at the end of the representative day.However, if the harvest period is longer than one day, this approximation method cannot be directly used to estimate the accumulation of algae biomass in the pond until the day of harvest.In order to understand the changes in culture properties such as biomass concentration and algal biomass accumulation for different harvest periods, we solved the original model for shorter periods, which were also set equal to the harvest period, than a year.We assumed that the representative day repeats itself.Here, we present the results from a 15-day period solution for the month of January.For this analysis, the data collected from Mathematica weather data database corresponds to year 2012.The rest of the analysis revealed similar conclusions to the one outlined below, and hence, the variables that change due to biomass accumulation are modelled using the same equation set described below.
Figure 2 shows changes in biomass concentration until harvest for the 15-day solution for the month of January.The x-axis shows day light time (sunrise through sunset) in a given day while y-axis shows the biomass concentration.Biomass concentration was constrained such that at sunrise of day-1, initial seeding took place.As the day progressed, the biomass concentration kept increasing until it reached the time of sunset.The concentration at sunrise of the next days was set equal to the concentration at sunset of the previous days until harvest (this assumption would yield a theoretical upper bound for the algal biomass production).The ratio of biomass concentration at sunrise to that at sunset is calculated for each day, and these ratios are given in Figure 2 at the bottom of x-axis.As can be seen from Figure 2, the ratios do not change significantly.The constant ratio allows the development of an expression that calculates the biomass concentration at harvest if the biomass concentrations at the time of seeding (which is a parameter for our original model and should be known) and at the time of sunset (which can be obtained if the dynamics of one representative day is included in the optimization model) are known.This expression is given in Equation (53).BCharv hd pt " sunsetq " BC pt " sunrise ´1q ¨" BC pt " sunsetq BC pt " sunrise ´1q In Equation ( 53), BCharv hd (t = sunset) is the biomass concentration at the time of harvest, BC(t = sunrise-1) and BC(t = sunset) are the biomass concentrations at sunrise and sunset of the representative day, and harvP is the harvest period.
The above relation can be extended to production variables as well.Equation (54) relates accumulation of algae biomass, Xharv c,hd (t), from the first day of harvest period to the last day of harvest period, to mass produced at sunrise and sunset, X c (t), and harvest period, harvP.Here, X c (t), (g) is the hourly production of the component c produced in a raceway pond on representative day.Xharv c,hd pt " sunsetq " X c pt " sunrise ´1q ¨" X c pt " sunsetq X c pt " sunrise ´1q  harvP @c P tDA, W I Au (54) now, the total annual production can be calculated by the sum of all the accumulated biomass, Xharv c,hd (t), from the sunrise of the first day to the sunset of the last day of harvest period times the number of harvest cycles, harvC, in a month as shown in Equation (55).This equation replaces Equations ( 27) and ( 28) in our model.
harvC ¨Xharv c,hd pt " sunsetq @y, @c P tDA, W I Au (55) Equations ( 43) and ( 44) can be rewritten as Equations ( 56) and ( 57) for calculating the energy requirements when considering only the times between sunrise and sunset for the whole year.In Equations ( 56) and ( 57), the multiplication with 30 converts the daily energy requirements to annual energy consumption.
EPP " 30 Annual industrial water requirements, Indwater (m 3 ), to run a raceway pond and to compensate for the evaporated water is calculated via Equation (58).
Equations ( 2) through ( 5), ( 7) through ( 26), ( 29) through (42), and ( 45 CONOPT option Rtmaxv = 1 ˆ10 30 was used as the upper limit on the variables.Maximum number of domain errors was set to 1000 and maximum number of solver iterations was set to 5000.Because CONOPT is a local optimization solver, we used a multiple-start approach, where the model variables were initialized to different values generated using Latin Hypercube Sampling technique, in an effort to find the global solution.

Case Study
For the case studies, the optimum raceway pond geometries are obtained for each combination of two species (P.tricornutum and I. galbana) and four locations (Tulsa, USA; Hyderabad, India; Cape Town, South Africa; and Rio de Janeiro, Brazil).The case study data for the species and model parameters are shown in Tables A1 and A2 of Appendix.The weather data-maximum temperatures, minimum temperatures, relative humidity and wind velocity for locations are presented in Tables A3 and A4 of Appendix.The breakdown of total capital investment coefficient (TCI l ) and total production cost coefficient (TPC l ) is shown in Section 6.1.In this work, all the cost coefficients have been inflated to the year 2012.The annual algae biomass demand in Equation ( 25) is five tons per annum.The harvest period and harvest cycle in each month were taken to be six days and five cycles, respectively.It was assumed that the raceway is seeded just before sunrise after each harvest to yield a biomass concentration of 0.6 g¨m ´3 [11].
Sensitivity analysis is performed on the best case results, i.e., species and location combination that yielded the lowest net present cost, to understand the influence of harvest time, demand, channel length to width ratio, and biomass initial seeding concentration on the optimum solution.

Economic Assumptions
The financial assumptions applied in the present study are summarized below.These costs are taken from literature for the year 1996 [35].
Estimation of capital investment cost (TCI l ) ˝Direct costs: They include site preparation, grading, compaction; pond levees, geotextiles; mixing (paddle wheels); CO 2 diffusers; harvesting; flocculation; centrifugation; water and nutrient supply system; waste treatment (blow down); building and structures; electrical supply and distribution; instrumentation and machinery; and land costs [35].The total direct costs add to $63,700/ha ˝Manufacturing costs: They include utilities-power required for harvesting, processing, lighting up buildings; nutrients; flocculants; waste disposal; maintenance and repairs = 5% of fixed-capital investment; operating supplies = 15% of maintenance and repairs; depreciation = fixed-capital investment/number of recovery years; local taxes and insurance = 4% of fixed-capital investment; financing = 5% of total capital investment; and labor and overhead costs.The total manufacturing costs add to $28,600.08/ha/year˝General expenses = 20% of total production cost = $7,150.02/ha/year˝Total production cost = manufacturing costs + general expenses = $35,750.1/ha/yearThese costs have been inflated to current values depending on chemical engineering cost index.Hence, the total capital investment (TCI l ) was calculated as $13.144 m ´2 and total production cost (TPC l ) was calculated as $5.475 m ´2¨year ´1.Here, cost index for 1996 and 2012 were taken as 382 and 585, respectively.Cost of electricity and water are taken from literature [35] as ElCost l : 0.07 $ (kW¨h)

´1
and WtCost l : 3.18 $ (1000 gal) ´1, respectively.It should be noted that these economic parameters are taken same for all locations to understand the sole influence of algae species, geographical location, and raceway pond geometry on one another.

Problem Statistics and Solution Approach
The resulting optimization problems (NLPs) have 4360 variables and 4219 equations for the case of Hyderabad location; 4484 variables and 4339 equations for the case of Tulsa location; 4391 variables and 4249 equations for the case of Cape Town location; and 4360 variables and 4219 equations for the case of Rio de Janeiro location.The length of the day (the period between sunrise and sunset) differs from location to location due to the difference in latitude, longitude, and time zone.This causes the variation in the total number of variables and equations.Based on these, the degrees of freedom is three: channel depth, channel width, and biomass concentrations.Note that biomass concentration is a time-dependent variable, and its size depends on the location.
The values of channel depth, channel width, and biomass concentrations were initialized using Latin Hypercube Sampling technique.The rest of the model variables were calculated using these initializations.This sampling method is chosen because it spreads the sample points more evenly across all possible values.It partitions each input distribution into N intervals of equal probability, and selects one sample from each interval.In our model, we generated 100 samples to initialize the decision variables.Each of the eight cases was solved 100 times using the Latin Hypercube initializations.As noted before, to obtain a local solution, the variables of the model should be initialized to a feasible or a near-feasible solution.An iterative initialization scheme is implemented in GAMS version 24.2.1 to draw model initialization to a feasible or a near-feasible solution.The model is first attempted to be solved with an initialization from the Latin Hypercube samples.If it is solved to optimality, the model is initialized with the next Latin Hypercube sample.If the model is not solved to optimality, the model is initialized using the values of channel depth and channel width of the feasible or infeasible solution obtained by the CONOPT version 3.15 M. The biomass concentrations are initialized to that at sunset of the first representative day.The initializations for the rest of the variables are calculated based on these three variable values.This process is repeated for eight iterations or until an optimum solution is obtained.If an optimum solution is not obtained at the end of eight iterations, the initialization scheme moves to the next Latin Hypercube initialization point.The solutions presented in the next section (Section 6) are the local optimum solutions that yielded the lowest costs for each algae species and location pair.

Results and Discussion
Table 1 summarizes the raceway pond geometries and averages of important culture properties for eight case studies.The table allows us to identify the best combination of species, location, and pond geometry, and helps in identifying the promising alternatives.The variables that are given in Table 1 are the net present cost (Z) which is the objective function for each case study; pond geometry such as the channel depth (CD), channel width (CW), pond length (length), Surface area occupied by the pond (SA) and volume occupied by the pond (V); average annual physical properties such as average irradiance inside the pond (I avg ), channel velocity (U avg ), and pond temperature (T pond ); average annual culture properties such as biomass concentration (BC) and specific growth rate (µ); and average daily productivities such as volumetric productivity (PrV) and areal productivity (PrA).Average annual culture and physical properties are highlighted here to give an idea of how much the average values would be for a given year.However, the culture properties and physical properties fluctuate on a daily basis based on the climatic conditions.These fluctuations are shown via Figure 3 through 8.
The results from Table 1 reveal that the geometry of the open-channel raceway pond is different for different algae species and geographical location because the culture properties such as biomass concentration, growth rate, and biomass productivities change with respect to the algae species and local climatic conditions.The optimal raceway pond designs presented in Table 1 suggest that amongst the two algae species, I. galbana species, and amongst the four locations, Hyderabad, India, (Case 6) meets the annual demand of 5 tons algae biomass at the minimum net present cost of $61,109 for a plant life of 10 years.The cost of dry algae was found to be $3.28 kg ´1.It can be observed that I. galbana species require shallower and wider ponds and Hyderabad, India location requires shallower and narrower ponds for growing algae.I. galbana species is selected over P. tricornutum species for algae biomass production because I. galbana species has higher light extinction coefficient, I k .This enables higher light absorption inside the pond which triggers higher productivity and lower costs.This behavior can be observed in Table 1.Cases 5 through 8 where I. galbana species is cultured requires smaller pond volumes when compared to Cases 1 through 4 where P. tricornutum species is cultured, except for Case 7. We hypothesize that the solution of Case 7 is a bad local minima.
Hyderabad, India is selected over the others for algae biomass production because the average irradiance inside the pond at this location triggers higher biomass concentration, growth rates, and productivities compared to that at other locations.As expected, at higher average pond temperatures, the average algae growth rates are higher, and higher average biomass concentrations in the pond lead to lower biomass production costs (Table 1).Cases 2 and 5 where I. galbana and P. tricornutum species are cultured in Hyderabad, India, require smaller pond volumes when compared to other cases.
Figure 3 shows the variation in sunlight at the best location to grow algae, i.e., Hyderabad, India.The x-axis defines the clock times in 24 h format for each representative day of the month, and the y-axis is the solar irradiance (µE¨m ´2¨s ´1).It is evident that as the day progresses, sunlight keeps increasing until mid-day and later decreases as it comes close to nightfall.It can also be observed that as summer approaches, solar irradiance tends to build up until the month of March and later declines as winter approaches.Figure 4 shows how average light intensity (I avg ) inside the pond varies with day light time of the representative day for I. galbana species cultured in Hyderabad, India.The x-axis defines the times in a 24 h format for each representative day of the month and y-axis is the average light intensity (µE¨m ´2¨day ´1).It follows a trend similar to that of solar irradiance (I o ).It can be observed that in a representative day, the average light intensity inside the pond gradually increases from the time of sunrise until mid-day and later decreases as the sun goes down.This trend continues for all representative days of the year.However, on an overall scale, the light intensity is maximum inside the pond during the months of February, March, April, and May. Figure 5 shows the biomass concentration profile (g¨m ´3) in the same format.On each representative day just before sunrise, biomass concentration is initiated to 0.6 g¨m ´3.During the course of the day with sunlight, the concentration in the pond builds up continuously until it reaches sunset.The concentration has rapidly raised from sunrise to sunset in the months of summer when compared to the winter months.Figure 6 shows the trend in biomass production (X f and X g ).The x-axis shows day light time (sunrise through sunset) in a given day while y-axis shows the biomass production.It follows the same trend as biomass concentration.It is clear that as production rises, concentration of biomass inside the pond rises.In Figure 7, x-axis is the clock times in 24 h format for each representative day of the month and y-axis is the growth rate (µ).For a given day, growth rate of algae biomass tends to increase in the first few hours until mid-day and later decreases.This behavior of deteriorating growth rate is mainly due to the decrease in the pond temperature in the afternoon hours.Algal growth rate is rapid in the summer months of March, April, and May because temperature of the pond oscillated around the optimal growth temperature of I. galbana (~27 ˝C).      Figure 8 shows the fluctuation of surrounding temperature (T surr ) and raceway pond temperature (T pond ) with time.Surrounding temperature depends on maximum and minimum temperature of the day, longitude, and time zone.It rises from sunrise until mid-day and then after goes down as the sun goes down.Temperature of the pond is driven by surrounding temperature and thermal gains/losses in the pond due to solar radiation, air radiation, pond radiation, evaporative heat transfer, convective heat transfer, and heat transfer due to inflow water.It is also dependent on the growth rate of algae as shown in Equations ( 4) and (5).The temperature of the pond rises as the day progresses and reaches a maximum after which it decreases.The decrease in pond temperature is mostly due to the cooling effects of evaporation.It can be observed that the surrounding temperature and pond temperature are high during the warmer months of March, April, and May.
Figure 9 shows the annual variation in heat flow in open-channel raceway ponds.The x-axis shows the representative days of the respective months in a year and y-axis shows the six heat fluxes considered for energy balance in the pond: radiation from the pond (Q p ), radiation from the air (Q a ), radiation from the sun (Q s ), evaporative heat flux (Q ev ), convective heat flux (Q cv ), and inflow heat flux (Q i ).Radiation from the pond (Q p ) depends on pond temperature (T pond ).Hence, rate of heat flow from the pond, (Figure 9A), follows a similar trend as pond temperature in Figure 8.It can be observed that at the beginning of the day, less radiation is emitted by the pond and as the day progresses, the pond temperature increases and so does the radiation from the pond.Radiation from the air (Q a ) (Figure 9B) follows a similar trend as surrounding temperature (T surr ) (Figure 8), and radiation from the sun (Q s ) (Figure 9C) follows a similar trend as irradiance at the location (I o ) (Figure 3) as they depend on one another.It can be observed that evaporative heat flux (Q ev ) (Figure 9D) is dominant in the month of May when the surrounding temperature is high.The increase in the heat transfer due to evaporation during this month has led to cooling in the pond, and hence, there is a decrease in the pond temperature.Evaporation increases from sunrise to mid-day and later decreases with the course of the day.However, in summer months, at sunrise when the surrounding temperature is higher than the pond, evaporation is high.Rate of heat flow due to convection (Q cv ) (Figure 9E) is higher in May and June similar to that of evaporation heat flux.Inflow heat flux (Q i ) (Figure 9F) is the rate of heat transfer between the water that is added to the pond, that is lost during evaporation, and the pond water itself.It can be observed that during the months of May and June, when the rate of evaporation is high, the pond water is relatively cooler.The water added to the pond is assumed to be at the surrounding temperature.Hence, heat transfer from inflow water to the pond water is high during those months.Among all the heat fluxes, rate of heat transfer due to evaporation plays a major role in determining the temperature of the pond.

Sensitivity Analysis
Sensitivity analysis was performed on the best case, i.e., Case 6, to understand the influence of model parameters on the overall design of raceway pond and on the cost of algae biomass production.The studied parameters include the harvest times (harvP), upper bound on pond length, lower bound on channel depth, and initial biomass concentration (BCini).
Table 2 summarizes the results obtained by changing the harvest periods (harvP).The harvest periods studied are one, two, three, six, and ten days.It is evident that for longer harvest periods, harvP, i.e., when harvesting frequency is low, biomass concentration increases until the time of harvest because of accumulation of algae biomass, and hence, growth rates are high.It can also be observed from Table 2 that higher biomass concentrations require lower channel depth, which enables higher average light intensity inside the pond.Longer harvest periods are also marked by higher volumetric and areal productivities.The productivities are high because longer harvest periods allows accumulation of the product for longer time inside the pond.For the harvest period of 10 days, the channel depth, channel width, and channel length are 4 cm, 1.45 m, and 48 m, respectively.However, at the channel depth of 4 cm, evaporation occurs rapidly leading to biomass concentrations that are higher than 10,000 g¨m ´3.This would violate the model assumptions used in this work.In the current model, we assumed biomass concentrations to be below 10,000 g¨m ´3 in order to inhibit light attenuation due to dense cultures and keeping the fluid properties similar to that of water.Thus, from our analysis, the next better option to harvest the algae biomass is six days.Since, the length of the pond for a six day harvest period case was found to be the upper bound on pond length, sensitivity analysis was performed by increasing the upper bound of pond length to see its influence on the net present cost.
Table 3 summarizes the results obtained by changing the upper bound on pond length to see its influence on the net present cost and the culture properties.To meet the demand, as the upper bound on the pond length increases, the length tends to follow the upper bound.This results in an increase in the channel depth, making the channel width smaller.In the case of 6-SA6, as the pond becomes longer, channel width became smaller than channel depth, making the pond narrower and deeper.From our analysis, it is evident that longer ponds are economical for algae biomass production.However, report from National Renewable Energy Laboratory suggests that pond length should range between 10 and 300 m.
Table 4 summarizes the results obtained by changing the lower bound on channel depth.The case with lower pond depth leads to higher biomass concentration, growth rate, and productivity.This is the result of higher average light intensity inside the pond due to lower pond depths.From Table 4, it is clear that lower pond depth results in lower net present cost.Our analysis suggests that ponds depth should be kept as low as the pumping and construction allows.
Table 5 summarizes the results obtained by changing the seeding, i.e., initial biomass concentration (BCini), in the pond at the time of sunrise.The case with higher initial biomass concentration proved to be economical, since higher initial seeding leads to better cultural properties such as average biomass concentration and productivities in the pond.As the initial biomass seeding concentration increases, the surface area required to meet the demand decreases.This leads to lower production costs of algae biomass.Therefore, the analysis suggests that higher initial seeding concentrations are required to obtain lower overall production costs.However, average light intensity inside the pond deteriorates with increasing biomass concentration because of mutual shading of the algae biomass.

Figure 1 .
Figure 1.Schematic of an open channel raceway pond.Here, length denotes the pond length (m), q the pond width (m), p the channel length (m), and CW the channel width (m).The schematic shows a paddle wheel used for mixing purpose and the direction of flow.

Figure 2 .
Figure 2. Schematic of the variables listed in Equation (53).The figure shows how biomass concentration changes from day-1 to day-15 (until harvest).The figure in the inset shows the accumulation over the period of day-1.
˝Indirect costs: Engineering and contingency = 15% of direct costs excluding land costs = $9,255/ha ˝Fixed-capital investment = direct costs + indirect costs = $72,955/ha ˝Working capital = 15% of total capital investment ˝Total capital investment = fixed-capital investment + working capital = $85,829/ha ˝10% rate of return ˝10 years plant life ˝Straight line depreciation ˝Cost index for 2012 Estimation of total production cost (TPC l )

Figure 4 .
Figure 4. Variation in average light inside the pond shown in a representative day of the year.

Figure 5 .
Figure 5. Concentration profile of I. galbana species cultured in Hyderabad, India.

Figure 6 .
Figure 6.Trend that shows the production of algae biomass with time on all representative days of the year.

Figure 7 .
Figure 7. Trend that shows growth rate with time on all representative days of the year.

Figure 8 .
Figure 8. Temperature profile of the pond during the course of the day.

Figure 9 .
Figure 9. (A) Rate of heat flow from the pond; (B) Rate of heat flow from air; (C) Rate of heat flow from the sun; (D) Rate of heat flow due to evaporation; (E) Rate of heat flow due to convection; (F) Inflow heat flux during the course of the day; (G) Variation in heat fluxes during the course of the day.
(19)objective function, Equation(19), is the net present cost.It includes capital costs: Z Cap-cost and operating costs: Z Electric-mixing , Z Electric-pumping , Z Water , and Z Product-cost .The Minimum Acceptable Rate of Return (MARR) is 15%.Equation (

Table 1 .
Summary of results for the eight Nonlinear Programming (NLP) Models.

Table A3 .
List of average temperature parameters for the locations.

Table A4 .
List of average relative humidity and wind velocity parameters for the locations.They are taken the same for all times of the day.