Article
Article
Article
Chemistry &
Engineering
View Article Online
REVIEW View Journal
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
A microorganism culture process is a complex system in which physical operating parameters and
biological responses strongly interact. Mathematical formulation and modelling of the different phenomena
involved in the process enable a better understanding of the behaviour of the process, and therefore
enable the process parameters to be defined accordingly. The contribution of a model, even a simple one,
is highly beneficial to the understanding of the process. The definition of a model for the particular case of
photobioreactors is not easy, however, and requires the integration of multiple and often complex
knowledge. This article reports a review on the biological aspects of the photosynthetic microorganisms
Received 1st December 2020, culture necessary to model kinetic growths, the designs of photobioreactor used for deep analysis of the
Accepted 14th April 2021
physiologic aspects of the microalgae culture and for the industrial culture. The different approaches to
DOI: 10.1039/d0re00450b
model the kinetic growth are described together with the modelling of the radiative field and its coupling
to a simple biological model in order to illustrate the particular influence of the light factor, which is the
rsc.li/reaction-engineering main specific feature of photobioreactors.
This journal is © The Royal Society of Chemistry 2021 React. Chem. Eng.
View Article Online
carbon source. Mixotrophy is obtained when the two sources In the first part, a reminder of photosynthesis and the
are mixed. In this article, only autotrophic cultures are associated production of microalgae is given, as well as the
considered. Two main families of culture systems can be factors limiting their growth. Then, some examples of culture
distinguished: open reactors and closed reactors (PBR). system design are described, with focus on laboratory culture
Although expensive due to their high level of technology, the systems, which are necessary to study the behaviour of
latter allow total control of the reactor and growing conditions microalgae and to get information on growth kinetics, on
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
and give access to performances far superior to those of the open, closed and developing systems. The last part concerns
obtained in extensive cultivation systems (raceway type). This is the different modelling approaches.
mainly due to that PBRs allow a certain control level of
operating conditions obtaining reproducible conditions for the Photosynthesis and production of
culture9 and all nutritional limitations can be deleted with the microalgae
exception of light. The growth rate of the microorganisms is
then controlled solely by the thermodynamic efficiency of The term “microalgae” includes, strictly speaking, any
photon use,10 corresponding, to photolimited growth. This microscopic algae. This includes microalgae (eukaryotes) as
mode of cultivation makes it possible to obtain performances well as prokaryotic organisms (cyanobacteria) and
close to maximum productivity on the surface (mass of photosynthetic bacteria. Three modes of culture are to be
biomass produced per unit area). For a reduction of the distinguished: autotrophy, which consists in the feeding of
impacts of biomass production, several ways are being studied, an inorganic carbon source, heterotrophy, an organic carbon
such the recycling of industrial flue gases as a source of carbon, source, and mixotrophy, when inorganic and organic sources
the recycling of the culture media to limit water consumption, are mixing. In this article only autotrophic cultures, then
the use of waste water as a basis for culture media with in cultivation by photosynthesis, are considered. Photosynthesis
addition environmental benefits. In fact, the CO2 from is described by the following equation:
industrial plants can be fixed for mitigation of the greenhouse
light
gas emission as well as nitrogen or phosphorus removal in CO2 þ H2 O → ½CH2 O þ O2 (1)
wastewater treatment processes.11 The majority of the culture
systems are built in a small scale medium (maximum a few ha) This process converts light energy into chemical energy (ATP
for a total worldwide biomass production around 30 000–50 000 and NADPH) which is then used by cells to synthesize
tonnes of dry matter per year. This production is obviously not organic carbon from inorganic carbon (CO2 or other
enough to meet the demand for biofuel, an insufficient dissolved inorganic carbon). This conversion is made
quantity for the real worldwide demand.12,13 Moreover, the possible by the succession of two types of reactions: (i) light
production costs of biofuels from microalgae are too high to (or photochemical) reactions and dark (or biochemical)
compete in energy markets. To further reduce the production reactions. The light reactions take place in the thylakoids.
cost of the system, it is necessary to employ the use of Thylakoid membranes have photosynthetic pigments
advanced control strategies to ensure an optimization of the grouped into photosystems. Two photosystems called PSI and
system is a way to reduce the production cost.14 To control the PSII absorb light at slightly different wavelengths. Incident
pH of the culture and reducing significantly the cost and the photons having a wavelength within a domain of radiation,
CO2 losses of the system, techniques based on model predictive called photosynthetically active radiation (PAR, 400 nm < λ <
control have obtained successful results in this target.15–17 It is 700 nm) are picked up by a photon collecting antenna,
necessary to take into account all the system variables to have a named here light harvesting complex II (LHCII), composed of
cost reduction by a good use of the resources. For example, pigments (chlorophyll a and b and protective carotenoids
Ifrim et al. (2013)14 developed a nonlinear multivariable (PPC, photo-protective carotenoids)). The energy of photons
controller based on dynamical model with an exact feedback having an energy equivalent higher than that of a photon
linearization to control biomass concentration and pH by emitted at 680 nm or 700 nm (E = hc/λ) is degraded to the
acting on the dilution rate and the injected carbon dioxide gas same energy level accepted by the corresponding
flow rate. Depending on the control complexity, different types photosystem. Excess energy is dissipated in the form of heat
of system models are developed. Fernández et al. (2014)18 has or fluorescence by the collector antenna. The energy
developed a dynamic model for microalgae production in equivalent to one photon emitted at 680 nm (PSII) or 700 nm
tubular photobioreactor for the prediction of the main (PSI) is then transmitted to the reaction centers (primary
parameters influencing the microalgae growth rate: electron donors), called P680 and P700) by proteins D1 and
temperature, pH, dissolved oxygen and biomass concentration. D2. This contribution of energy at the reaction center will
Calibration and validation tests have been made in an outdoor change it from a P680 state to an excited state noted P680*,
tubular photobioreactor. A review on the use of Computational thus releasing an electron. The latter will be transferred by
Fluid Dynamics (CFD) for the optimization of bioreactor design the plastoquinone (PQ) which will take two protons available
and for the study of the interaction of hydrodynamics, light in the stroma (inside the chloroplast) to switch to its oxidized
supply, heat and mass transfer with biological kinetics has state PQH2. The electrons recovered by the PQ are removed
been done by Pires et al. (2017).19 from the water by oxidation (photolysis of water). This
React. Chem. Eng. This journal is © The Royal Society of Chemistry 2021
View Article Online
reaction takes place at the level of the tyrosine Z complex group. The G3P molecule is a three-carbon sugar. For every
and allows the release of 4 protons, 4 electrons and 1 oxygen twelve moles of G3P synthesized, only two are transported to
molecule per molecule of oxidized water. The protons are the metabolic pathway of sugars, the remainder being
released in the lumen (inside the thylacoid). The PQ directed to the third phase of the cycle, the regeneration of
transmits its two electrons to the cytochrome b6f and rejects RuBP. This last phase completes the Calvin cycle by
two protons into the lumen. The cytochrome b6f also pumps synthesizing RuBP molecules from G3P under the action of
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
two protons from the stroma into the lumen, thus ATP to regenerate the first CO2 acceptor.
contributing to accentuate the proton gradient between the The efficiency of photosynthesis is directly related to the
inside (highly charged proton) and the outside of the amount at the photon rate received at a point. This rate is
thylacoid (weakly charged). The two electrons are then called irradiance (G) and its unit is the micromole of photons
transferred to a second electron carrier, the plastocyanine absorbed per square meter of surface area per second (μmolh
(PC). The PC will transmit one electron at a time to the P700 m−2 s−1). Three behaviours are distinguished: (i) a
(PSI). The P700 is similar to P680, except that it operates at photolimited regime, in which photosynthesis increases with
an energy equivalent of one photon at 700 nm. The P700 will irradiance, (ii) a second, so-called photosaturated regime, in
accumulate two electrons, which will take it to its excited which photosynthesis is independent of the light received,
state P700*. He will then give up his electrons to ferredoxin and (iii) a photoinhibition regime in which the efficiency of
(Fd). The ferredoxin transfers its electrons to the enzyme photosynthesis decreases with the increase in irradiance
ferredoxin NADP reductase (FNR) which will allow the received by the microorganism. It should be added to this
reduction of NADP+ to NADPH, H+, which induces the that below a certain irradiance value, called compensatory
pumping of an additional proton from the stroma to the irradiance, and characterized by a zero oxygen balance at the
lumen. The proton gradient generated by the protolysis of cell level, the phenomenon of respiration is predominant
the water is used as a proton-motor force for the synthesis of over photosynthesis. This is characterised by negative
ATP by the enzyme ATP synthase at from ADP and inorganic biomass production (consumption of carbon reserves).
phosphate. There are two electron transfer paths: an acyclic
pathway, leading to the synthesis of ATP and NADPH, H+ by Factors limiting growth of microalgae
photophosphorylation (Z-schema), and a second, called cyclic
electron transfer. This is a direct transfer of a electron from Many parameters can affect the productivity of a
ferredoxin to PQ. This results in the pumping of protons by photobioreactor. Biological parameters such as bacterial and
the cytochrome, which increases the proton gradient and fungal contamination, predation by protozoa, or even
promotes the production of ATP, without producing of sometimes contamination by another microalgae than the
NADPH, H+. The specific production rates of ATP and desired one.22 Several physico-chemical parameters are also
NADPH, H+ directly affect the metabolism and growth of a influential: light energy supply, temperature, pH, salinity and
photosynthetic microorganism. The ratio of these two terms the supply of nutrients necessary for growth (including
(defined as the P/2e− by Cornet et al. (1998)20 will reflect the inorganic carbon). The homogeneity of temperature, pH and
adjustment of the cells energy metabolism. Similarly, the salinity, as well as the nutrient accessibility to microalgae and
phthosynthetic quotient, QP, is denoted, as the ratio between light access in the system are controlled by the hydrodynamics
the production rate of one mole of oxygen and the of the system, making it a key parameter for the optimization
consumption rate of one mole of carbon dioxide. Similarly, of PBR production. The dynamics of operation under solar
the ratio QRO2 (quantum requirement for oxygen production) conditions brings an additional particularity, a notion of
denotes the number of photons necessary to produce one instability over time, which does not allow the system to
mole of oxygen. Theoretically this value is eight photons per operate at the optimum throughout the day. This applies in
oxygen molecule.21 The ATP and NADPH, H+ molecules will particular to light, or even temperature and pH depending on
then be used for Calvin cycle (dark reactions). The dark the control conditions of the culture system.
reactions are related to carbon fixation take place through a
cycle called Calvin cycle. This cycle takes place in the stroma Light energy in photobioreactors
and uses the energy stored during the light phase in the form The efficiency of photosynthesis is directly related to the
of ATP and NADPH, H+ to fix inorganic carbon and amount of light absorbed by the microalgae. At the cell level, it
synthesize organic carbon. The first phase is the carbon is possible to know the quantity of light necessary to obtain the
fixation by the synthesis of 3-phosphoglycerate (3PGA) from maximum performance of the microalgae. When considering
CO2 and ribulose-1,5-biphosphate (RuBP), a reaction the photobioreactor system, it is easy to understand that if
triggered by an enzyme, RuBisCO. The second phase is the microalgae absorb part of the light transmitted to it, the rest of
reduction phase, during which each 3PGA molecule receives the microalgae contained in this system will share the
an additional phosphate from ATP (synthesized in clear remaining flux. As a consequence, light energy is
phase), then an electron pair released by the NADPH, H+ heterogeneous in the volume of the photobioreactor. It will be
molecule reduces the 1,3-biphosphoglycerate molecule to then impossible to maintain maximum performance for all the
D-glycéraldéhyde-3-phosphate (G3P), releasing a phosphate microorganisms in the PBR. When considering a flat PBR, the
This journal is © The Royal Society of Chemistry 2021 React. Chem. Eng.
View Article Online
attenuation of light, when transmitted homogeneously to the considering the light received by the microalgae in the reactor,
PBR, occurs in only one dimension, in the thickness of the characterized by irradiance G. This parameter γ is the ratio of
culture. The light attenuation profile depends on several the illuminated volume to the total volume of the PBR. The
parameters: (i) the amount of light transmitted at the surface boundary between these two zones is then given by the
of the PBR also called PFD (photon flux density in μmol m−2 compensating irradiance (Gc). When the local irradiance is
s−1), (ii) the concentration of microalgal biomass (CX in gX l−1) below this value, the zone is considered to be a dark volume, in
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
and (iii) the radiative properties of the microalgae, depending which the phenomenon of respiration is predominant, and
in particular on the pigment content of the microalgae (wpig in thus negatively impacts biomass productivity. The three
%) and the shape of the cultured microalgae.23 Pruvost and operating regimes are as follows (Fig. 1):
Cornet (2012)24 identified three light attenuation regimes - Photo-limited regime (case a): in this case γ < 1, the
corresponding to distinct performances. These are entire photon flux transmitted to the crop is absorbed. A dark
differentiated by the parameter γ, representing the illuminated zone is present due to the high biomass concentration. This
volume fraction in a PBR. In the literature, this is defined by is achieved by a high residence time of the microalgae and
Fig. 1 Light attenuation profiles for a planar PBR and associated productivities:24 (a) case for γ < 1, (b) case for γ = 1, (c) case for γ > 1 (this figure
was published in J. Pruvost and J.-F. Cornet, Microalgal Biotechnology: Potential and Production, 2012, De Gruyter, Berlin, Germany, 181–224,
Copyright De Gruyter).
React. Chem. Eng. This journal is © The Royal Society of Chemistry 2021
View Article Online
does not allow maximum biomass productivity due to the consumption of the cell's carbohydrates (energy source) in the
negative impact of respiration. metabolism of cell maintenance and protein synthesis.32
- Strict physical limitation or luminostat regime (case b): Photosynthetic microorganisms are unequal to the phenomenon
here γ = 1, i.e. all the light transmitted to the system is of respiration. Indeed, cyanobacteria will lose on average less
absorbed by the microalgae without the appearance of a biomass at night than the eukaryotic cells.33 Moreover, within
respiration zone (G(L) = Gc). It is this regime that leads to the eukaryotic microorganisms, there is a great disparity in
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
maximum performance of a PBR for a given light flux. respiration rates at night depending on the cultured.31 The most
- Kinetic regime (case c): here γ > 1, i.e. part of the light is influential parameter on the loss of biomass in the dark phase
not absorbed in the culture volume and is therefore of a given microalgae is temperature. Edmundson and
transmitted. This translates in energy terms into biomass Huesemann (2015)31 showed that when the temperature at night
productivity below the maximum achievable productivity was reduced from 25 °C to 10 °C, it was possible to reduce
(because not all the transmitted energy is converted). This biomass decrease. The most impressive case is that of the
regime is very particular because the residence time is low, microalga Nanochloropsis salina, which loses 20% of its dry
which leads to a decrease in biomass concentration and a matter concentration after ten hours in the dark at 25 °C,
risk of culture leaching. Moreover, the system performance whereas reducing the temperature to 10 °C over the same period
can be reduced by poor absorption of the luminous flux. only results in a decrease in concentration of around 2%. The
Indeed, due to the low light attenuation, the amount of slowing down of metabolic activity, and therefore of the decrease
energy absorbed by the microalgae is important, which can in biomass by temperature reduction has been confirmed by
induce a phenomenon of photosynthesis oversaturation or numerous studies on numerous microorganisms.32–35
even photoinhibition, damaging the photosynthetic In order to guarantee a good thermal management in a
apparatus of the microalgae and significantly reducing the PBR (to maximise production) and to be able to estimate the
performance of the PBR. energy needs linked to the control of a PBR, it is essential to
be able to represent the different exchanges taking place
within it. Knowing the exchanges taking place in a growing
Temperature and microalgal growth system allows to model the temperature evolution in a given
A non-optimal culture temperature affects the growth of culture system throughout the year. It also serves as a tool for
microalgae. Although not directly affecting photochemical dimensioning the thermal exchanger necessary for the
reactions, the temperature at which the microorganism grows regulation of a solar PBR (and for estimating the associated
plays a major role in enzyme activity.25 This has the costs), and thus as a tool for comparing temperature-
consequence of reducing the growth capacity of a given regulated closed PBRs. Solar radiation is responsible for the
microorganism when growth deviates from the optimum, heating of microalgae cultures, but other thermal exchanges
and in some cases even leading to cell death. In addition, a take place between a PBR and its environment.
change in the temperature of the medium can force the
microalgae to change its composition. This has been shown
for example on lipids.25 For the culture of the study Nutrients and carbon supply
microalga, C. vulgaris, the optimal growth temperature is For a photosynthetic microorganism, essential nutrients are
estimated to be between 20 and 30 °C.26–28 Most microalgae needed, at least a source of nitrate (or ammonium), sulphates
are so-called mesophilic, i.e. they have an optimal growth and phosphates ions. These nutrients are provided under
temperature between 15 and 40 °C. For the latter, a chemical salts dissolved in aqueous solution. The
temperature of the environment surrounding the microalgae composition of the culture medium is adjusted according to
above 45 °C for more than 24 hours can have an irreversible the specific needs of the cultured microorganisms. Nutrients
effect on the culture.29 There are other strains which, due to can be supplied in three different ways. A high dose is only
their natural environment, have totally different optimal given at the beginning of the cultivation, this is called a
growth temperatures, such as the so-called psychrophilic batch culture. The solution is brought continuously to the
(Topt <15 °C) and thermophilic (Topt >50 °C) microalgae.28 culture, a part of the culture is then continuously renewed),
No photosynthetic microorganisms with an optimal growth it is called continuous culture mode. The fed-batch mode
temperature above 75 °C have so far been identified, corresponds to the addition of new medium is done
probably due to the instability of chlorophylls at this punctually during the cultivation process. For the continuous
temperature.28 Generally, the temperature control must be mode, two cases can be distinguished: the chemostat mode
addressed for the evaluation of the technical feasibility of for which the supply of the medium is fixed by the flow rate
large-scale algae production.30 of the feed pump, and the turbidostat mode for which the
When producing microalgal biomass under solar conditions, supply of the medium is regulated according to the quantity
the night period is often a period when no production takes of biomass in the culture.
place. However, the loss of biomass due to respiration can be Carbon is the majority element in the biomass, it
major, sometimes up to 20% biomass loss for a ten-hour represents about half of the dry mass of the microalgae. In
night.31 The night period is generally accompanied by a autotrophic culture, it is brought into the medium in the
This journal is © The Royal Society of Chemistry 2021 React. Chem. Eng.
View Article Online
form of inorganic salts dissolved in the culture medium agitation speed should be around 15 cm s−1.38 In a raceway-
(HCO3− or CO32−) or in gaseous form (CO2) transferred to the type basin, the preferred agitation system is the paddle wheel
culture. When grown in an open system, the liquid is in for several reasons.36 The presence of vertical velocities
equilibrium with the atmosphere containing a very low allows good agitation and distribution of biomass in all
quantity of CO2(g), i.e. about 400 ppm. Consequently, if the dimensions. Moreover, this system causes very low shear and
medium contains dissolved carbon in a quantity greater than is suitable for all microalgae strains, including the most
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
that given by the equilibrium with the atmosphere, a fragile. Finally, the electrical consumption linked to agitation
desorption of CO2 in its aqueous form towards its gaseous with this type of system remains low (around 6 W m−2).
form will be generated by this concentration gradient. In Nevertheless, recent studies highlight the interest of stirring
addition, the form in which inorganic carbon is present is raceway systems by airlift, which would be up to three times
pH dependent. With a pH above 9, desorption is low because less energy consuming than the paddle wheel, while
only the aqueous form of CO2 is in equilibrium with the maintaining the same quality of agitation.39
gaseous phase. Therefore, at high pH, it is not necessary to A wide variety of PBRs technologies exist, such as tubular,
cover a culture system to avoid carbon limitation, this cylindrical or flat panel systems with specific mixing
explains why Spirulina platensis culture is generally grown in problems. However, the criteria for a “good” mixing can be
an open system. As is the case with temperature, each stated: avoid sedimentation of the biomass, give the same
microalga has an optimal pH for growth. For example, story about the light received to the whole algal population,
Chlorella vulgaris grows at pH values between 5 and 9 with an satisfy heat transfer to maintain an optimal temperature and
optimal pH between 7 and 8.36 For the cyanobacterium gas–liquid mass transfer for CO2 supply and O2 removal,
Anabaena, this optimum is between 9 and 10. During avoid biofilm formation and shear rates too high in relation
cultivation, when growth occurs, the total carbon to cell fragility. From the conception point of view, different
concentration decreases. Since CO2 is an acid, the pH then criteria have to be taken into account: cost, limitation of the
increases.37 Increasing this pH may reduce growth. One dark zones, limitation of the thickness of the culture media
method of ensuring the presence of carbon (avoiding to increase the volumic productivity, easiness to scale up and
limitation) and maintaining the pH at the optimal level is energy efficiency. The various PBR technologies available
therefore to occasionally add gaseous (acidic) CO2 to the have advantages and limitations in terms of hydrodynamics
culture when the pH exceeds the set value. conditions and biomass productivity relatively to the
construction cost.50,69,76,81
React. Chem. Eng. This journal is © The Royal Society of Chemistry 2021
View Article Online
This journal is © The Royal Society of Chemistry 2021 React. Chem. Eng.
View Article Online
the term “farming” appeared in the 1960s.22 The first systems to reduce the thickness of the culture as much as possible for
dedicated to the cultivation of microalgae were artificial the same ground surface area. For this purpose, a slightly
ponds, also called lagoons. This is the simplest production inclined plane system was developed (inclination of 1.7%,67
system. It consists of a body of water between 20 and 40 on which the microalgae culture flows down and is then
centimeters deep without agitation. This system was used by reinjected at the reactor head via a pump (Fig. 6).
the Cognis company in Hutt Lagoon (Australia) to cultivate This cultivation system reduces the costs of biomass
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
the microalga Dunaliella salina (Fig. 4). These microalgae are production and processing. Indeed, since production is only
cultivated for their high carotene content. Today, the global dependent on the lighting area, a reduction in thickness
company BASF farms these 700 hectares of land, making it will have under identical conditions no impact on the
the largest production site in the world.60 quantity of biomass produced. Agitation and separation of
Subsequently, agitated basins were developed in order to biomass from water will therefore be cheaper due to the
avoid sedimentation and thus optimize light distribution in small amount of water to be treated. An industrial version
the culture volume. Moreover, this avoids the presence of of this concept has been implemented in Portugal by the
concentration gradients or thermal stratification, reducing the company A4F in 2014. This system comprises a succession
overall performance of the process.42 Thus, in order to of four inclined planes for a total production area of 3000
guarantee homogeneity in the culture basin, Andersen (2005)61 m2 (Fig. 7).
recommends an average velocity between 20 and 30 cm s−1 to Even today, open systems are still the most widely used on
avoid the risk of sedimentation. This generally guarantees an industrial scale because of their low manufacturing cost
minimum local velocities in the basin greater than 10 cm s−1. and ease of extrapolation. However, they are not suitable for
This velocity is defined as the velocity at which the risks of all types of production. Indeed, these systems are reserved
biomass sedimentation appear (for most strains). The most for the production of extremophilic strains, because of the
popular open agitated systems on a large scale are circular and risk of contamination by airborne microorganisms. There is
raceway-type systems.114 Circular systems are agitated via an no problem for strains such as Dunaliella salina or Arthrospira
arm with a size of the radius of the basin set in motion by the platensis because these two photosynthetic microorganisms
axis. This type of basin has two major disadvantages, a non- grow in conditions where most contaminants cannot
optimal surface area (due to the circular shape) and difficulty proliferate (hypersalinity and highly alkaline medium,
in extrapolating beyond 1000 m2.62 respectively). The second point limiting production in this
Notoriety of circular systems remained relatively low as they type of system is the input of inorganic carbon. Since these
were supplanted by basins called raceways because of their systems are in contact with air (containing very little CO2),
shape. Initially developed in the field of wastewater treatment the dissolved inorganic carbon in the basins will tend to
by Oswald and Goleuke (1967)63 and Benemann and Oswald equilibrate with the atmosphere (except in the case of highly
(1996),64 this technology consists in a loop (or sometimes alkaline pH). This difference in inorganic carbon
several interconnected loops) consisting of two straight zones concentration will induce desorption of the latter from the
and two turning zones (Fig. 5) in which a thin layer of basin into the ambient air, which may lead to a limitation in
microalgae culture (about 20 to 40 centimeters for industrial nutrients, which will tend to significantly reduce the
systems) is set in motion, usually by means of a paddle wheel.
A final type of open system was developed in the 1960s in
the Czech Republic, in Trebon66,67 The aim of this system is
Fig. 5 Raceway of 3000 m2 in NBT, Israël (https://www.israel21c.org/ Fig. 6 Inclined plane system (200 litres) developed in Trebon, Czech
wp-content/uploads/2016/02/algae-NBT.jpg)65 (this figure was Republic68 (original figure was published in J. R. F. Malapascua, C. G.
published in P. J. Harvey and A. Ben-Amotz, Algal Research, 2020, 50, Jerez, M. Sergejevová, F. L. Figueroa and J. Masojídek, Aquat Biol.,
102002, Copyright Elsevier). 2014, 22, 123–140, under the Creative Commons CC-BY License).
React. Chem. Eng. This journal is © The Royal Society of Chemistry 2021
View Article Online
This journal is © The Royal Society of Chemistry 2021 React. Chem. Eng.
View Article Online
React. Chem. Eng. This journal is © The Royal Society of Chemistry 2021
View Article Online
Fig. 11 SunOleo photobioreactor (with the kind permission of the company SunOleo).
simply increasing the number of panels (Fig. 12). This observed. The parameters of the model are adjusted for
conception brings many advantages on key points in the each case, but are not necessarily representative of the
culture of photosynthetic microorganisms, namely (i) an physical, chemical or biological reactions that take place.
optimised contribution of light within the culture itself, The most widely used are those of Monod (1942)90 and
with thus an optimal exploitation of the photons emitted, Andrews–Haldane.91 Monod's model is commonly used to
(ii) a design limiting the risks of adhesion on the lighting describe bacterial growth. This model is based on the
structures, easily dismantled and facilitating cleaning, and assumptions that the yield of growth (in g of biomass per
(iii) a modular production that can be simply extrapolated g of substrate) of conversion of a given substrate
in volume while maintaining constant surface and volume (concentration of this substrate noted S, in g l−1) into
productivity. biomass remains constant over time, and that the
New types of photobioreactor with a coupling between a evolution the specific growth rate (μ, in h−1) with the
photobioreactor and a photovoltaic cell are currently in evolution of the substrate concentration considered follows
development.115,116 A pilot-scale spectrally-selective, a hyperbolic trend:
insulated-glazed photovoltaic flat plate photobioreactor was
developed with an infrared reflecting system embedded in μ = μmaxS/(KS + S) (2)
the illumination surface for the thermal regulation of
outdoor photobioreactors. The interest is to produce both where μmax is the maximum specific growth rate and KS
microalgal biomass and electricity and to increase the net represents the half-saturation constant of the microorganism
energy ratio.117 with respect to the substrate under consideration. This
constant is dependent on the microorganism and its culture
Modelling of the PBR conditions, which express the affinity of the microorganism
to the substrate. A limitation of the Monod model is the non-
Growth models inclusion of the inhibitor effect of a substrate present in
The overall approach generally consists of a mathematical excess. The Andrews–Haldane model91 introduced a term for
law that correctly translates the experimental behaviour the substrate inhibition:
This journal is © The Royal Society of Chemistry 2021 React. Chem. Eng.
View Article Online
value of the growth rate. The type of model is often used for metabolism in response to changes in the illumination
the good approximation of the experimental behaviours,92,93 conditions. Current state of the art of constraint-based
although its formulation is problematic in the case of modeling and computational method development are
photosynthetic microorganisms. There is thus no discussed by Tibocha-Bonilla et al. (2018).100
proportionality between respiration and luminous intensity.94 Light attenuation modelling. Light attenuation can be
Some approaches are based on artificial intelligence with represented by the Beer–Lambert equation using an
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
deep learning technology and data-driven surrogate modeling extinction coefficient representing the biomass absorption
framework.95 A stochastic optimization algorithm was used coefficient.101 Beer–Lambert's law is a relatively simple law
in order to develop models from hydrodynamic and allowing to model correctly the exponential attenuation of
biochemical kinetic results to optimise biosystems and monochromatic radiation through a homogeneous medium,
obtain decision-making for choosing the best parameters for absorbent and non-diffusive. It is a law which is very widely
microalgae production. used in spectrophotometry, especially for the determination
of chemical species in solution. This law does not take into
account the phenomena of radiation diffusion. Its simplicity
Predictive models of use means that this law, even if it neglects an important
A dynamic semi-predictive model for microalgal culture in phenomenon, namely diffusion, is still widely used for
tubular photobioreactor have been developed by Fernández modeling, in the first instance, the approximation, light
et al. (2014).18 The model is based on a photosynthesis rate transfer to PBR.101–104
equation taking into account light intensity and the others However, this model does not take into account the effects
most important variables (temperature, pH, and dissolved of light scattering by microalgae. For this reason, Cornet
oxygen) that influence the growth and performance of the et al. (1992, 1995)97,105 used a different approach to
culture in any microalgal production system. The accurately describe light attenuation in a photobioreactor. By
thermodynamic approach of the irreversible processes was adapting Schuster's (1905)106 model describing the behaviour
applied by Stucki (1978)96 to describe the biological of light in a foggy atmosphere, the proposed radiative
reactions observed in microorganism cultures. This transfer model describes light attenuation by taking into
approach consists in considering that, in a complex account light absorption by pigments and light scattering
biological process such as photosynthesis, the average state throughout the cell. The assumptions are as follows: the
of the measured reactions (not thermodynamically stable in medium is assumed to be absorbing, scattering and non-
instantaneous value, due to the autocatalytic steps of fluorescent. Therefore, it is sufficient to determine three
enzyme activation and inhibition) is stable on a parameters to characterize the light path: the mass
appropriated time scale.97 Cornet's analysis of the Z-scheme coefficient of light absorption (Eaλ in m2 kg−1), the mass
of photosynthesis has led to the establishment of a growth coefficient of scattering (Esλ in m2 kg−1) and the coefficient
model, which was compared with experiments on cultures of backscattering (b2λ, without unit). When this model
of two microorganisms: Arthrospira platensis and considers the light attenuation as monodirectional (on the
Rhodospirillaceae sp.10,98 This predictive method requires a z-axis), it is called a two-flux model, and takes into account
high degree of theoretical knowledge, since it is based on the propagation of light, for a given wavelength (in nm) along
the precise analysis of the major phenomena governing the z-axis and in two opposite directions, Iλ+ and Iλ−.97 The
growth, as the energetics of photosynthetic conversion, sum of these two specific intensities gives the local irradiance
anabolism and catabolism. G (eqn (4)):
Flux metabolic modeling. Metabolic flux analysis allows ðð
the modeling of intracellular metabolism in response to Gλ ¼ I λ dω ¼ I λ þ þ I λ − (4)
4π
genetic and/or environmental variations in a system
biological. All modeling is based on the reconstruction of a
where dω represents the solid angle defining a radiation
metabolic network associating each reaction involved for the
beam. The system of differential equations to be solved as a
conversion of a given substrate into products of interest.
function of culture depth is as follows:
Cogne et al. (2003)99 established the metabolic network of
autotrophic growth of Arthrospira platensis; 121 reactions dI λ þ
¼ −Eaλ ·CX ·I λ þ − Esλ ·b̄ 2λ ·CX ·ðI λ þ þ I λ − Þ (5)
associated with 134 metabolites made it possible to dz
reconstruct a metabolic network linked to the production of
C-phycocyanin. Cogne et al. (2011),48 based on two dI λ −
computational approaches of metabolic flux, modeled of the ¼ −Eaλ ·CX ·I λ − − Esλ ·b̄ 2λ ·CX ·ðI λ− þ I λ þ Þ (6)
dz
behavior of the eukaryotic microalgae Chlamydomonas
reinhardtii grown under photoautotrophic conditions. The The boundary conditions for a planar system, illuminated by
reconstruction of a metabolic network comprising 280 a collimated light source at normal incidence, having a rear
metabolic reactions linked to 278 metabolites was carried out. face of reflectivity ρ are as follows (qλ,0 is the surface incident
The study was able to demonstrate a reorientation of flux of the PBR):
React. Chem. Eng. This journal is © The Royal Society of Chemistry 2021
View Article Online
(
z ¼ 0; I λ þ ¼ qλ;0 and 10 times the wavelength of the received radiation. To
− þ
(7) solve these equations, the following parameters are
z ¼ L; I λ ¼ ρI λ
required:107,108 (i) the complex refractive index of the
Finally, the expression of irradiance as a function of the particle mλ, (ii) its size, and (iii) the refractive index of the
thickness of a planar PBR is given by Pottier et al. 2005:107 medium surrounding the particles nm,λ. The complex
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
This journal is © The Royal Society of Chemistry 2021 React. Chem. Eng.
View Article Online
ðL ðL J O2 C X M X
1 1 — rX ¼
hi ¼ ¼ EaGðzÞdz (17) (22)
L 0 L 0
vO2 –X
The MRPA can also be determined by carrying out a balance where MX and vO2–X represent, respectively, the C-molar mass of
on the photonic phase in the reactor volume.109 the biomass (kgX mol−1) and the stoichiometric coefficient of
S as oxygen production. In the reactor volume (VR), the MRPA value,
hi ¼ q0 − qL ¼ q0 − qL (18)
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
React. Chem. Eng. This journal is © The Royal Society of Chemistry 2021
View Article Online
Fig. 14 Comparison of the volumic productivity in PBR of different sizes and comparison the model of eqn (27)–(29) with for Neochloris
oleoabundans: ρM = 0.8, φM = 1.83 × 10−9 kg per μmole, K = 90 μmole per m2 s−1.
This journal is © The Royal Society of Chemistry 2021 React. Chem. Eng.
View Article Online
PFD. The positive effect of increasing the PFD on biomass 12 L. Brennan and P. Owende, Renewable Sustainable Energy
productivity and the negative effect of increasing the depth of Rev., 2010, 14, 557–577.
culture can be seen. 13 S. Taras and A. Woinaroschy, Comput. Chem. Eng., 2012, 43,
10–22.
Conclusion 14 G. A. Ifrim, M. Titica, M. Barbu, L. Boillereaux, G. Cogne, S.
Caraman and J. Legrand, Chem. Eng. J., 2013, 218, 191–203.
The development of microalgae culture processes requires
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
React. Chem. Eng. This journal is © The Royal Society of Chemistry 2021
View Article Online
37 D. Kaplan, A. E. Richmond, Z. Dubinsky and S. Aaronson, Vitamins, Pigments and Growth Factors, ed. E. J. Vandamme,
Algal Nutrition, in Handbook of Microalgal Mass Culture, Elsevier Science Publishers LTD, England, 1989, Elsevier
CRC Press, Inc., Boca Raton, Florida, 1986, pp. 147–98. Applied Biotechnology Series, pp. 15–26.
38 J. C. Weissman, R. P. Goebel and J. R. Benemann, 61 R. A. Andersen, Algal Culturing Techniques, Academic press,
Biotechnol. Bioeng., 1988, 31(4), 336–344. 2005.
39 S. Nacir, F. Zouhir, N. Ouazzani, H. Jupsin, J.-C. Vasel and 62 D. Chaumont, J. Appl. Phycol., 1993, 5(6), 593–604.
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
L. Mandi, Rev. Sci. Eau, 2012, 25(3), 287–299. 63 W. J. Oswald and C. G. Goleuke, Technical Report, 1967,
40 M. A. Borowitzka, J. Biotechnol., 1999, 70(1–3), 313–321. Large-Scale Production of Algae, PB-296411, California
41 C.-Y. Chen, K.-L. Yeh, R. Aisyah, D.-J. Lee and J.-S. Chang, Univ., Berkeley (USA), http://www.osti.gov/scitech/biblio/
Bioresour. Technol., 2011, 102(1), 71–81. 5398149.
42 O. Pulz and K. Scheibenbogen, Apoptosis, Advances in 64 J. R. Benemann and W. J. Oswald, Final Report, 1994,
Biochemical Engineering Biotechnology, Springer, Berlin Report to DOE-NETL http://www.osti.gov/bridge/purl.cover.
Heidelberg, 1998, vol. 59, pp. 123–152. jsp?purl=/137315-0uSjuX/webviewable/.
43 Z. Csögör, M. Herrenbauer, I. Perner, K. Schmidt and C. 65 P. J. Harvey and A. Ben-Amotz, Algal Res., 2020, 50, 102002.
Posten, Chem. Eng. Process.: Process Intesif., 1999, 38(4–6), 66 I. Šetlík, V. Šust and I. Málek, Algol. Stud., 1970, 111–164.
517–523. 67 J. Doucha and K. Lívanský, Algol. Stud., 1995, 129–147.
44 J. Pruvost, L. Pottier and J. Legrand, Chem. Eng. Sci., 68 J. R. F. Malapascua, C. G. Jerez, M. Sergejevová, F. L.
2006, 61(14), 4476–4489. Figueroa and J. Masojídek, Aquat. Biol., 2014, 22, 123–140.
45 H. Takache, G. Christophe, J.-F. Cornet and J. Pruvost, 69 O. Pulz, Appl. Microbiol. Biotechnol., 2001, 57(3), 287–293.
Biotechnol. Prog., 2010, 26(2), 431–440. 70 J. Doucha and K. Lívanský, J. Appl. Phycol., 2006, 18(6),
46 H. Takache, J. Pruvost and H. Marec, Algal Res., 2015, 8, 811–826.
192–204. 71 P. Chen, Int. J. Agric. Biol. Eng., 2009, 2(4), 1–30.
47 A. Martzolff, E. Cahoreau, G. Cogne, L. Peyriga, J.-C. 72 M. Tredici, Handbook of Microalgal Culture, ed. A.
Portais, E. Dechandol, F. Le Grand, S. Massou, O. Richmond, Blackwell Publishing Ltd., 2003, pp. 178–214.
Gonçalves, J. Pruvost and J. Legrand, Biotechnol. Bioeng., 73 C. Jiménez, B. R. Cossio and F. Xavier Niell, Aquaculture,
2012, 109(12), 3030–3040. 2003, 221(1–4), 331–345.
48 G. Cogne, M. Rugen, A. Bockmayr, M. Titica, C. G. Dussap, 74 M. R. Tredici and R. Materassi, J. Appl. Phycol., 1992, 4(3),
J.-F. Cornet and J. Legrand, Biotechnol. Prog., 2011, 27(3), 221–231.
631–640. 75 Y.-K. Lee, S.-Y. Ding, C.-S. Low, Y.-C. Chang, W. L. Forday
49 A. Taleb, J. Pruvost, J. Legrand, H. Marec, B. Le Gouic, B. and P.-C. Chew, J. Appl. Phycol., 1995, 7(1), 47–51.
Mirabella, B. Legeret, S. Bouvet, G. Peltier, Y. Li-Beisson, S. 76 E. Molina Grima, F. G. Acién Fernández, F. García Camacho
Taha and H. Takache, Bioresour. Technol., 2015, 177, and Y. Chisti, J. Biotechnol., 1999, 70(1–3), 231–247.
224–232. 77 Y. Chisti, Biotechnol. Adv., 2007, 25(3), 294–306.
50 A. P. Carvalho, L. A. Meireles and F. Xavier Malcata, 78 M. R. Tredici and G. Chini Zittelli, Biotechnol. Bioeng.,
Biotechnol. Prog., 2006, 22(6), 1490–1506. 1998, 57(2), 187–197.
51 J. C. Ogbonna, T. Soejima and H. Tanaka, J. Biotechnol., 79 A. Prokop, R. K. Bajpai and M. E. Zappi, Algal
1999, 70, 289–297. Biorefineries: Volume 2: Products and Refinery Design,
52 A. Ramos de Ortega and J. C. Roux, Biomass, 1986, 10(2), Springer, 2015.
141–156. 80 M. R. Tredici, N. Bassi, M. Prussi, N. Biondi, L. Rodolfi, G.
53 B. Tamburic, F. W. Zemichael, P. Crudge, G. C. Maitland Chini Zittelli and G. Sampietro, Appl. Energy, 2015, 154,
and K. Hellgardt, Int. J. Hydrogen Energy, 2011, 36(11), 1103–1111.
6578–6591. 81 M. Morweiser, O. Kruse, B. Hankamer and C. Posten, Appl.
54 J. Degen, A. Uebele, A. Retze, U. Schmidt-Staigar and W. A. Microbiol. Biotechnol., 2010, 87(4), 1291–1301.
Trosch, J. Biotechnol., 2001, 92, 89–94. 82 C. Posten, Eng. Life Sci., 2009, 9(3), 165–177.
55 G. Chini Zittelli, L. Rodolfi and M. R. Tredici, J. Appl. 83 J. P. Diaz, C. Inistrosa and F. G. Acien Fernandez, Process
Phycol., 2003, 15, 107–114. Biochem., 2019, 86, 1–8.
56 M. Harker, A. J. Tsavalos and A. J. Young, J. Ferment. 84 V. Goetz, F. Le Borgne, J. Pruvost, G. Plantard and J.
Bioeng., 1996, 82, 113–118. Legrand, Chem. Eng. J., 2011, 175, 443–449.
57 K. Kaewpintong, A. Shotipruk, S. Powtongsook and P. 85 W. Blanken, M. Janssen, M. Cuaresma, Z. Libor, T. Bhaiji
Pavasant, Bioresour. Technol., 2007, 98, 288–295. and R. H. Wijffels, Biotechnol. Bioeng., 2014, 111(12),
58 D. O. Hall, F. G. A. Fernandez, E. C. Guerrero, K. K. Rao 2436–2445.
and E. M. Grima, Biotechnol. Bioeng., 2003, 82, 62–73. 86 B. Ji, W. Zhang, N. Zhang, J. Wang, G. A. Lutzu and T. Liu,
59 Y. Watanabe and H. Saiki, Energy Convers. Manage., Bioprocess Biosyst. Eng., 2013, 37(7), 1369–1375.
1997, 38, 499–503. 87 A. Ozkan, K. Kinney, L. Katz and H. Berberoglu, Bioresour.
60 L. J. Borowitzka and M. A. Borowitzka, β-Carotene Technol., 2012, 114, 542–548.
(Provitamin A) Production with Algae, in Biotechnology of 88 J.-F. Cornet, Chem. Eng. Sci., 2010, 65(2), 985–998.
This journal is © The Royal Society of Chemistry 2021 React. Chem. Eng.
View Article Online
89 J. Pruvost, J.-F. Cornet and J. Legrand, World Pat., 105 J.-F. Cornet, C. G. Dussap, J.-B. Gros, C. Binois and C.
WO2018096107, 2018. Lasseur, Chem. Eng. Sci., 1995, 50(9), 1489–1500.
90 J. Monod, Recherches sur la croissance des cultures 106 A. Schuster, Astrophys. J., 1905, 21, 1.
bactériennes, Hermann & cie, Paris, 1942. 107 L. Pottier, J. Pruvost, J. Deremetz, J.-F. Cornet, J. Legrand
91 J. F. Andrews, Biotechnol. Bioeng., 1968, 10(6), 707–723. and C. G. Dussap, Biotechnol. Bioeng., 2005, 91(5),
92 C.-Y. Chen, P.-C. Kao, C.-J. Tsai, D.-J. Lee and J.-S. Chang, 569–582.
Published on 20 April 2021. Downloaded by University of California - Santa Barbara on 5/16/2021 9:07:10 PM.
Bioresour. Technol., 2013, 145, 307–312. 108 R. Kandilian, J. Pruvost, A. Artu, C. Lemasson, J. Legrand
93 O. Bernard and B. Rémond, Bioresour. Technol., and L. Pilon, J. Quant. Spectrosc. Radiat. Transfer, 2016, 175,
2012, 123(0), 520–527. 30–45.
94 P. Falkowski, Z. Dubinsky and K. Wyman, Limnol. 109 J. Pruvost, J.-F. Cornet and J. Legrand, Chem. Eng. Sci.,
Oceanogr., 1985, 30(2), 311–321. 2008, 63, 3679–3694.
95 E. A. del Rio-Chanona, J. L. Wagner, A. Haider, F. Fiorelli, 110 R. Kandilian, J. Pruvost, J. Legrand and L. Pilon, Bioresour.
D. Zhang and K. Hellgardt, AIChE J., 2019, 65, 915–923. Technol., 2014, 163, 308–319.
96 J. Stucki, Prog. Biophys. Mol. Biol., 1978, 33, 99–187. 111 J. Pruvost, J.-F. Cornet, F. Le Borgne, V. Goetz and J.
97 J.-F. Cornet, C. G. Dussap and G. Dubertret, Biotechnol. Legrand, Algal Res., 2015, 10, 87–99.
Bioeng., 1992, 40(7), 817–825. 112 J. Pruvost, F. Le Borgne, A. Artu, J.-F. Cornet and J. Legrand,
98 J.-F. Cornet, L. Favier and C. G. Dussap, Biotechnol. Prog., Advances in Chemical Engineering, Elsevier, 2016, vol. 48, pp.
2003, 19, 1216–1227. 257–310.
99 G. Cogne, J.-B. Gros and C. G. Dussap, Biotechnol. Bioeng., 113 A. Souliès, C. Castellain, T. Burghelea, J. Legrand, J.-F.
2003, 84(6), 667–676. Cornet, H. Marec and J. Pruvost, Biotechnol. Prog., 2016, 32,
100 J. D. Tibocha-Bonilla, C. Zuñiga, R. D. Godoy-Silva and K. 247–261.
Zengler, Biotechnol. Biofuels, 2018, 11, 241, DOI: 10.1186/ 114 I. Hamed, Compr. Rev. Food Sci. Food Saf., 2016, 15,
s13068-018-1244-3. 1104–1123.
101 C.-G. Lee, Biotechnol. Bioprocess Eng., 1999, 4(1), 78–81. 115 E. G. Nwoba, D. A. Parlevliet, D. W. Laird, K. Alameh and
102 E. Molina Grima, A. F. G. Fernandez, G. F. Camacho and Y. N. R. Moheimani, Algal Res., 2019, 39, 101433.
Chisti, J. Biotechnol., 1999, 70, 231–247. 116 E. G. Nwoba, D. A. Parlevliet, D. W. Laird, K. Alameh and
103 C. H. Hsieh and W. T. Wu, Biochem. Eng. J., 2009, 46(3), N. R. Moheimani, Algal Res., 2020, 45, 101731.
300–305. 117 E. G. Nwoba, D. A. Parlevliet, D. W. Laird, K. Alameh, J.
104 E. Molina Grima, J. Fernandez, F. G. Acién and Y. Chisti, Louveau, J. Pruvost and N. R. Moheimani, Appl. Energy,
J. Biotechnol., 2001, 92(2), 113–131. 2020, 271, 115403.
React. Chem. Eng. This journal is © The Royal Society of Chemistry 2021