# Impact of Dropwise Condensation on the Biomass Production Rate in Covered Raceway Ponds

Jack Hoeniges, Keyong Zhu, Jeremy Pruvost, Jack Legrand, Dominique

Grizeau, El-Khider Si-Ahmed, Laurent Pilon

► To cite this version:

Jack Hoeniges, Keyong Zhu, Jeremy Pruvost, Jack Legrand, Dominique Grizeau, et al. Impact of Dropwise Condensation on the Biomass Production Rate in Covered Raceway Ponds. Energies, 2021, 14, 10.3390/en14020268. [https://doi.org/10.3390/en14020268](https://doi.org/10.3390/en14020268)

HAL Id: hal-03609173

https://hal.science/hal-03609173

Submitted on 15 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Impact of Dropwise Condensation on the Biomass Production Rate in Covered Raceway Ponds

Jack Hoeniges<fnref n="1" />,<fnref n="2" /> ID, Keyong Zhu<fnref n="3" />, Jeremy Pruvost<fnref n="2" />, Jack Legrand<fnref n="2" /> ID, El-khider Si-Ahmed<fnref n="2" /> and Laurent Pilon<fnref n="1" />,* ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Hoeniges, J.; Zhu, K.; Pruvost, J.; Legrand, J.; Si-Ahmed, E.-k.; Pilon, L. Impact of Dropwise Condensation on the Biomass Production Rate in Covered Raceway Ponds. *Energies* **2021**, *14*, 268. https://doi.org/10.3390/en14020268

Received: 2 December 2020

Accepted: 22 December 2020

Published: 6 January 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** This study investigates the effect of condensed water droplets on the areal biomass productivity of outdoor culture systems with a free surface, protected by a transparent window or cover to prevent contamination and to control the growth conditions. Under solar radiation, evaporation from the culture causes droplets to condense on the interior surface of the cover. To quantify the effect of droplets on the system's performance, the bidirectional transmittance of a droplet-covered window was predicted using the Monte Carlo ray-tracing method. It was combined with a growth kinetics model of *Chlorella vulgaris* to predict the temporal evolution of the biomass concentration on 21 June and 23 September in Los Angeles, CA. A droplet contact angle of 30° or 90° and a surface area coverage of 50% or 90% were considered. Light scattering by the condensed droplets changed the direction of the incident sunlight while reducing the amount of light reaching the culture by up to 37%. The combined effect decreased the daily areal biomass productivity with increasing droplet contact angle and surface area coverage by as much as 18%. Furthermore, the areal biomass productivity of the system was found to scale with the ratio $X_0/a$ of the initial biomass concentration $X_0$ and the specific illuminated area $a$, as previously established for different photobioreactor geometries, but even in the presence of droplets. Finally, for a given day of the year, the optical thickness of the culture that yielded the maximum productivity was independent of the window condition. Thus, the design and operation of such a system should focus on maintaining a small droplet contact angle and surface area coverage and an optimum optical thickness to maximize productivity.

**Keywords:** light transfer; photobioreactor; biofuel; bioprocess engineering; design; modeling

## 1. Introduction

Biofuels derived from microalgae have the potential to serve as a renewable and carbon-neutral alternative to conventional transportation fuels [1]. Furthermore, microalgae have garnered interest as a source of food and dietary supplements that use significantly less resources than conventional agriculture due to their rapid growth and high nutrient and protein content [2]. Cost-effective production of algal biomass is essential to produce economically viable algae-based products. While microalgae grow readily in a variety of conditions, optimizing their growth poses various technical challenges [3]. Cultivation systems must be optimized to deliver light and nutrients to the microalgae while maintaining optimum growth conditions such as temperature, pH, and dissolved O₂ and balancing costs associated with construction, operation, land use, water use, and nutrient sourcing [4].

Microalgae cultivation systems are generally classified as open or closed. Open systems such as raceway ponds are inexpensive to build and operate and have been widely

<footnote n="1">¹ Mechanical and Aerospace Engineering Department, Henry Samueli School of Engineering and Applied Science, University of California, Los Angeles, CA 90095, USA; jhoeniges@g.ucla.edu</footnote>

<footnote n="2">² Université de Nantes, Oniris, CNRS, GEPEA, UMR 6144, F-44600 Saint-Nazaire, France; jeremy.pruvost@univ-nantes.fr (J.P.); jack.legrand@univ-nantes.fr (J.L.); el-khider.si-ahmed@univ-nantes.fr (E.-k.S.-A.)</footnote>

<footnote n="3">³ School of Aeronautical Science and Engineering, Beihang University, Beijing 100191, China; zhukeyong@buaa.edu.cn</footnote>

<footnote>* Correspondence: pilon@seas.ucla.edu</footnote>
used for large-scale microalgae cultivation [4]. However, they can be easily contaminated; they lack effective temperature control; they suffer from water loss by evaporation; and they have relatively high harvesting costs due to their low cell concentration compared to closed systems [3]. Their operation also results in CO₂ degassing [5]. By contrast, closed cultivation systems, also referred to as photobioreactors (PBRs), typically feature a lower risk of contamination and higher biomass productivity than open systems [3,6]. However, PBRs are typically more costly, and their scale up remains challenging [3]. Specifically, in flat horizontal or inclined culture systems, a transparent window covering the culture and containing an over-pressurized gas space can be used to (i) limit culture contamination, (ii) reduce CO₂ degassing [5], and (iii) improve the control of the growth conditions [3]. Such an outdoor PBR consisting of a race pond with a paddle wheel sealed by a transparent plastic cover is illustrated in Figure 1a. The use of a transparent window enables the cultivation of a wider variety of species and thus greater potential for the production of high value-added products without significantly increasing construction and maintenance costs [3,6]. Regardless of the specifics of the cultivation system, water evaporation occurs when an outdoor microalgae culture is exposed to sunlight. When the culture is covered by a transparent window, evaporation leads to the formation of water droplets on the interior surface of the PBR cover, as illustrated in Figure 1b. However, the presence of droplets on the back side of a window has been shown to reduce its transmittance and change the direction of the transmitted light [7,8]. Both of these phenomena could impact the PBR biomass productivity [4,9]. However, to the best of our knowledge, the impact of droplet-covered windows on PBR performance has not been investigated to date.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Photographs of (a) an outdoor photobioreactor at the R&D facility AlgoSolis (Saint-Nazaire, France) with a clear window cover and (b) condensed droplets formed on the back side of a photobioreactor (PBR) window during a summer day.</figcaption></figure>

The objective of this study is to quantify the effect of pendant droplets on the cover of outdoor microalgae cultivation systems on biomass productivity. To do so, light transfer was simulated through the PBR cover supporting pendant condensate droplets. The predicted bidirectional transmittance served as an input into combined light transfer and growth kinetics models in the culture of *Chlorella vulgaris*. Particular attention was paid to the effect of droplet contact angle, surface area coverage, culture depth, initial biomass concentration, and day of the year on the daily biomass production rate.

## 2. Background

### 2.1. *Light Transfer through Droplet-Covered PBR Windows*

The interaction between light and droplet-covered windows has been studied both analytically [7,8,10–12] and experimentally [11,13–17]. Zhu et al. [8] simulated visible light transfer through a transparent window covered with non-absorbing cap-shaped water droplets on its back side. Light transfer was modeled using the Monte Carlo ray-tracing method for a wide range of droplet contact angle $\theta_c$ and surface area coverage $f_A$. Interestingly, both the normal-hemispherical transmittance $T_{nh}$ and the directional-hemispherical transmittance $T_{dh}$ of the droplet-covered window were found to be independent of droplet
size distribution and droplet spatial arrangement. Similar results were observed for the bidirectional transmittance $T_{bd}$ [7]. The normal-hemispherical transmittance $T_{nh}$ was nearly independent of contact angle $\theta_c$ until it reached the critical angle $\theta_{cr}$ for total internal reflection at the droplet/air interface, i.e., for $\theta_c < \theta_{cr}$. However, for $\theta_c > \theta_{cr}$, the normal-hemispherical $T_{nh}$ and directional-hemispherical transmittance $T_{dh}$ decreased monotonically with increasing droplet surface area coverage $f_A$ and/or increasing contact angle $\theta_c$ until reaching a minimum at $\theta_c = 90^{\circ}$. Similar results were found experimentally by previous studies [11,17].

Zhu and Pilon [18] investigated the effect of absorption by the window and/or droplets on the transmittance of a window covered with cap-shaped droplets on its back side. In the case of a weakly absorbing window or droplets with absorption indices of $k_w = 10^{-6}$ and $k_d = 10^{-4}$, respectively, the normal-hemispherical transmittance $T_{nh}$ was observed to decrease only slightly. Furthermore, the dependence on droplet contact angle $\theta_c$ was the same as that for a non-absorbing window and droplets. Thus, in the case of water on plexiglass, whose absorption indices are less than $1.6 \times 10^{-7}$ and $1.8 \times 10^{-6}$, respectively, over much of the photosynthetically active radiation (PAR) region from 400 to 700 nm, absorption by the window and droplets can be neglected [19,20].

## 2.2. Light Transfer in Microalgae Culture

Light transfer through a well-mixed suspension of microalgae is governed by the radiative transfer equation (RTE) for a homogeneous absorbing, scattering, and non-emitting medium. The two-flux approximation can be used as an analytical solution to the one-dimensional RTE for flat plate PBRs or raceway ponds [21,22]. This approximation has been used extensively [21–25] and has been successfully validated against results from a 3D RTE solver for ponds and flat plate PBRs [22].

Souliès et al. [9] theoretically and experimentally investigated the effect of incidence angle on the growth of *Chlorella vulgaris* in thin flat panel PBRs artificially illuminated by collimated visible light with an incidence angle $\theta_i$ of 0° or 60°. The intensity of the light source was adjusted such that the PBRs received the same radiative flux for both incidence angles. Nonetheless, the experimentally measured volumetric biomass productivity was observed to decrease from $0.029 \text{ kg m}^{-3} \text{ h}^{-1}$ for the PBR under normal incidence to $0.026 \text{ kg m}^{-3} \text{ h}^{-1}$ for the PBR under oblique incidence at $\theta_i = 60^{\circ}$. Compared to normal incidence, light propagating through the culture at an oblique angle had a longer pathlength to reach a given depth. This caused light to be fully absorbed at shallower depths and prevented it from penetrating deeper into the PBR. Thus, a smaller fraction of the culture volume was illuminated, and the volumetric PBR biomass productivity decreased.

Pruvost et al. [25] theoretically investigated the maximum achievable productivity for outdoor PBRs cultivating *Arthrospira platensis*. Two scenarios were considered. The first was an idealized tracking PBR wherein the solar flux was normally incident throughout the day. The second was a fixed horizontal PBR that experienced oblique incidence during the day. Here, the oblique incidence was observed to decrease (i) the photon flux density incident on the PBR and (ii) the penetration depth of sunlight into the PBR culture, as demonstrated by Souliès et al. [9]. Together, these two effects resulted in a 31% decrease in the areal biomass productivity from 55 tons ha$^{-1}$ year$^{-1}$ to 38 tons ha$^{-1}$ year$^{-1}$ for the “ideal” tracking PBR and the PBR under oblique incidence, respectively.

## 2.3. Microalgae Growth Kinetics

The time rate of change of the biomass concentration $X(t)$ in a PBR operated in batch mode can be modeled as [26]:

$$ \frac{dX}{dt} = \bar{\mu}(t)X(t) \qquad (1) $$

where $\bar{\mu}(t)$ is the average specific growth rate (in h$^{-1}$) at time t. Various models have been proposed to predict the average specific growth rate $\bar{\mu}$ of microalgae as a function of the local spectral fluence rate $G_{\lambda}(z)$. Takache et al. [27] developed a growth kinetics
model based on an energetic analysis of the excitation energy transfer in the photosynthetic apparatus. This model accounted for light limitation, photoinhibition, and cell respiration activity as a function of the local fluence rate [27]. Souliès et al. [9] further developed this model to account for the cell respiration activity based on the local rate of photon absorption (LRPA).

Lee et al. [22] studied the effect of PBR geometry on its daily biomass productivity. Microalgae growth in outdoor open ponds, vertical flat plate PBRs, and tubular PBRs cultivating *Chlamydomonas reinhardtii* was modeled accounting for light saturation and photoinhibition using the Haldane model [28]. The biomass loss rate from cell respiration was assumed to be constant [22,29]. Remarkably, for all PBR configurations, the biomass productivity was found to scale with the ratio of the initial biomass concentration $X_0$ and the specific illuminated area $a$ (in m$^{-1}$) given by $a = S/V$ where S is the illuminated surface area and V is the culture volume [22,30]. In addition, the maximum PBR biomass productivity occurred at the same value of $X_0/a$ for all PBR geometries. Similar results were obtained from both experimental and theoretical data reported in previous studies [31,32] despite the use of different microalgae species and growth models.

The present study aims to quantify the negative effect of condensed droplets on the window of outdoor covered raceway ponds on the time-dependent LRPA, biomass concentration, growth rate, and daily biomass productivity. This was achieved by coupling simulations of light transfer through droplet-covered windows with a simplified model for light transfer through the microalgae culture and a growth kinetics model for various solar conditions, droplet contact angle, and surface area coverage.

# 3. Methods

## 3.1. Problem Statement

Let us consider an outdoor raceway pond located in Los Angeles, CA, USA (34.07° N, 118.44° W) exposed to collimated solar radiation $G_{S,\lambda}$ on the summer solstice (21 June) and on the autumn equinox (23 September). The PBR was covered by a transparent horizontal window of thickness $H$ and refractive index $n_w$, as illustrated in Figure 2. The solar zenith angle $\theta_z$ was defined with respect to the outward pointing normal vector of the window surface, and the solar azimuth angle $\gamma_s$ was defined with respect to the due south direction where $\gamma_s = 90^\circ$ corresponds to due west and $\gamma_s = -90^\circ$ corresponds to due east. Figure 2a shows the situation when the PBR window was dry and droplet-free. Figure 2b schematically illustrates the PBR window partially covered by condensed droplets on its back side with surface area coverage $f_A$. The droplets were non-absorbing and cap-shaped with refractive index $n_d$, diameter $d$, projected diameter $d_p$, and contact angle $\theta_c$. The incident radiation was either back-scattered or transmitted through the dry or through the droplet-covered window into the PBR culture undergoing reflection or refraction at each interface. The PBR contained a well-mixed culture of depth $L$ and effective refractive index $n_m$ growing *Chlorella vulgaris* with time-dependent biomass concentration $X(t)$ (in $g L^{-1}$) and spectral average mass absorption $\bar{A}_{abs,\lambda}$ and scattering $\bar{S}_{sca,\lambda}$ cross-sections (in $m^2 kg^{-1}$), and spectral scattering phase function $\Phi_{\lambda}(\Theta)$. The culture depth $L$ varied from 0.1 m to 0.3 m, and the initial biomass concentration $X_0$ varied from 0.01 $g L^{-1}$ to 1 $g L^{-1}$.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Side view (not to scale) of the covered outdoor PBRs studied with (a) a dry window and (b) a droplet-covered window. In both cases, a portion of the solar radiation was reflected at each interface before the radiative flux $q_{m,\lambda,j}$ was transmitted into the microalgae culture at angle $\theta_m$. Scattering by cap-shaped droplets on the window's back side was described by the bidirectional transmittance $T_{bd}$.</figcaption></figure>

## 3.2. Assumptions

Light transfer through the droplet-covered windows was simulated using the Monte Carlo ray-tracing method [7,8,12]. To do so, the following assumptions were made [8]: (1) all interfaces were considered optically smooth. In other words, surface roughness was considered to be much less than the wavelength of light in the photosynthetically active region ($\lambda$ = 400–700 nm) such that reflection and transmission were specular and Snell's law and Fresnel's equations were valid. (2) The dimensions of the droplets and windows were much larger than the wavelength of light such that geometric optics prevailed. (3) Both the window and the droplets were considered to be non-absorbing over the PAR region. (4) The droplets were considered to be cap-shaped with a constant curvature. (5) The droplets were monodisperse with an arbitrary diameter $d = 100$ µm since the bidirectional transmittance of transparent windows with non-absorbing droplets was found to be independent of the droplet size [7,8]. (6) The windows simulated had a surface area coverage $f_A$ of 50% or 90%. Although local conditions may impact the amount of condensate present on the window, the range of $f_A$ was selected to represent a moderate and a worst-case-scenario droplet coverage. (7) Droplet contact angle $\theta_c$ was 30° or 90°. These two angles were chosen to approximately represent water on glass and water on a hydrophobic cover, respectively. (8) The refractive index of the air, window, water droplets, and microalgae culture were taken as constant across the PAR region and equal to $n_a = 1.0$, $n_w = 1.5$, $n_d = 1.33$, and $n_m = 1.33$, respectively. Thus, light scattering by the droplets was independent of wavelength. (9) The diffuse part of the incident sunlight was neglected such that all solar radiation incident upon the PBR window was collimated, and its direction ($\theta_s, \gamma_s$) varied during the course of the day.

Predictions of light transfer and microalgae growth in the outdoor PBR were based on the assumptions that (10) the PBR was operated in the light-limited regime wherein growth was only a function of the amount of light available to the microalgae [9,22,25]. Thus, the PBR was not limited by any other operating parameters such as temperature, pH, and nutrient availability, and changes in PBR biomass productivity could be attributed to light scattering by the droplets. (11) The microalgae culture was well-mixed with uniform biomass concentration $X(t)$ throughout the PBR. (12) The cells were randomly oriented. (13) The liquid medium was non-scattering and non-absorbing over the PAR region. (14) The radiation characteristics of *Chlorella vulgaris* were constant throughout the day and taken from the literature [9]. (15) All walls of the PBR were considered non-reflective and perfectly absorbing. (16) Edge effects and shading from the side walls were considered negligible so that light transfer in the PBR culture could be treated as one-dimensional and depended only on the culture depth $z$.
### 3.3. Light Transfer through Droplet-Covered PBR Windows

For PBRs covered by dry windows, Figure 2a illustrates that the incoming collimated spectral solar irradiance $G_{S,λ}$ (in µmol$_{h/v}$ m$^{-2}$ s$^{-1}$) was reflected and refracted at the air/window interface, at the window/air interface, and at the air/microalgae culture interface before being transmitted into the PBR culture. The transmissivity at the air/window and the air/microalgae culture interface, denoted by $τ_{a/w}$ and $τ_{a/m}$, respectively, was determined from Fresnel’s relation given by [33]:

$$τ_{1/2} = 1 - \frac{1}{2} \left[ \frac{\tan^2 (\theta_1 - \theta_2)}{\tan^2 (\theta_1 + \theta_2)} + \frac{\sin^2 (\theta_1 - \theta_2)}{\sin^2 (\theta_1 + \theta_2)} \right] \quad (2)$$

where $θ_1$ and $θ_2$ are the incident and refracted angles across the Medium 1/Medium 2 interface given by Snell’s law $n_1 \sin θ_1 = n_2 \sin θ_2$. Then, the spectral radiative flux at the upper boundary of the microalgae culture $q_{m,λ}$ (Figure 2a) was expressed as a function of the solar position ($θ_z, γ_s$) and the transmissivity at each interface as:

$$q_{m,\lambda}(\theta_z, \gamma_s) = \tau_w \tau_{a/m} G_{S,\lambda}(\theta_z, \gamma_s) \cos \theta_z. \quad (3)$$

where $τ_w$ is the transmissivity of the window given by $τ_w = τ_{a/w} / (2 - τ_{a/w})$ [34].

For droplet-covered windows, the incident sunlight was scattered into the transmitted direction described by the polar $θ_t$ and azimuthal $φ_t$ transmission angles (Figure 2b). The resulting angular distribution of the transmitted light was described by the bidirectional transmittance $T_{bd}$ in the transmitted direction ($θ_t, φ_t$) defined as [33,34]:

$$T_{bd}(\theta_z, \gamma_s, \theta_t, \varphi_t) = \frac{I_{t,\lambda}(\theta_z, \gamma_s, \theta_t, \varphi_t)}{G_{S,\lambda}(\theta_z, \gamma_s) \cos \theta_z} \quad (4)$$

where $I_{t,λ}$ is the spectral radiative intensity transmitted in the direction ($θ_t, φ_t$). Note that $T_{bd}$ is wavelength independent since the optical properties of the window and droplets were assumed to be constant over the PAR region.

In order to be coupled with the one-dimensional light transfer analysis in the microalgae culture, the bidirectional transmittance $T_{bd}$ ($θ_z, γ_s, θ_t, φ_t$) was averaged over the transmitted azimuthal angle $φ_t$ to yield the one-dimensional bidirectional transmittance $\bar{T}_{bd}$ ($θ_z, γ_s, θ_t$) according to [7]:

$$\bar{T}_{bd}(\theta_z, \gamma_s, \theta_t) = \frac{1}{2\pi} \int_0^{2\pi} T_{bd}(\theta_z, \gamma_s, \theta_t, \varphi_t) d\varphi_t. \quad (5)$$

Then, $\bar{T}_{bd}$ ($θ_z, γ_s, θ_t$) (in sr$^{-1}$) gave the average intensity transmitted in direction $θ_t$ due to solar radiation from the incident direction ($θ_z, γ_s$). The one-dimensional bidirectional transmittance $\bar{T}_{bd}$ ($θ_z, γ_s, θ_t$) was then used to express the transmitted intensity $I_{t,λ}$ through the droplet-covered window into direction $θ_t$ as:

$$I_{t,\lambda}(\theta_z, \gamma_s, \theta_t) = G_{S,\lambda}(\theta_z, \gamma_s) \cos \theta_z \bar{T}_{bd}(\theta_z, \gamma_s, \theta_t). \quad (6)$$

### 3.4. Light Transfer in Microalgae Culture

For wet windows, droplets scattered the incident sunlight into different directions in the transmission hemisphere. The latter can be discretized into M directions ($θ_{t,j}$) $1 \le j \le M$. The spectral radiative flux $q_{m,λ,j}$ reaching the upper boundary of the microalgae culture in the $θ_{t,j}$ direction can be written as:

$$q_{m,\lambda,j}(\theta_z, \gamma_s, \theta_{t,j}) = \tau_{a/m} I_{t,\lambda,j} \cos \theta_{t,j} \Delta\Omega_{t,j} \quad (7)$$

where the solid angle is given by $\Delta\Omega_{t,j} = 2\pi \sin \theta_{t,j} \Delta\theta_{t,j}$. Due to refraction, radiation transmitted through the window at angle $θ_{t,j}$ was then transmitted through the air/microalgae
culture interface at an angle $\theta_{m,j} = \sin^{-1}(n_a \sin \theta_{t,j}/n_m)$. Then, according to the two-flux approximation, the resulting local fluence rate $G_{\lambda,j}(z)$ at depth $z$ within the culture due to the spectral radiative flux $q_{m,\lambda,j}$ incident at angle $\theta_{t,j}$ on the upper boundary of the culture was expressed as [21]:

$$ \frac{G_{\lambda,j}(z, \theta_{t,j})}{q_{m,\lambda,j}(\theta_{t,j})} = \frac{2}{\cos \theta_{m,j}} \frac{(1 + \alpha_{\lambda})e^{\delta_{\lambda,j}(L-z)} - (1 - \alpha_{\lambda})e^{-\delta_{\lambda,j}(L-z)}}{(1 + \alpha_{\lambda})^2 e^{\delta_{\lambda,j}L} - (1 - \alpha_{\lambda})^2 e^{-\delta_{\lambda,j}L}} \quad (8) $$

where $\alpha_{\lambda}$ and $\delta_{\lambda,j}$ are given by [21]:

$$ \alpha_{\lambda} = \sqrt{\frac{\bar{A}_{abs,\lambda}}{\bar{A}_{abs,\lambda} + 2b_{\lambda}\bar{S}_{sca,\lambda}}} \quad \text{and} \quad \delta_{\lambda,j} = \frac{\alpha_{\lambda}X}{\cos \theta_{m,j}} (\bar{A}_{abs,\lambda} + 2b_{\lambda}\bar{S}_{sca,\lambda}). \quad (9) $$

here, the average spectral mass absorption $\bar{A}_{abs,\lambda}$ and scattering $\bar{S}_{sca,\lambda}$ cross-sections in $m^2 kg^{-1}$ (shown in Figure A1) are the properties of the microalgae species being cultivated and were taken from [9] for a cell mean radius of 2 µm and chlorophyll a, chlorophyll b, and photoprotective carotenoid (PPC) concentrations of 3.3 wt.%, 0.85 wt.%, and 0.85 wt.%, respectively. Assuming axisymmetric scattering by the culture, the backward scattering ratio $b_{\lambda}$ is defined as [21]:

$$ b_{\lambda} = \frac{1}{2} \int_{\pi}^{\pi/2} \Phi_{\lambda}(\Theta) \sin \Theta d\Theta. \quad (10) $$

here, $\Phi_{\lambda}(\Theta)$ is the so-called scattering phase function representing the probability that photons in the solid angle $d\Omega_i$ about direction $\hat{s}_i$ are scattered into the solid angle $d\Omega$ about direction $\hat{s}$ with $\Theta$ defined as the angle between the scattered direction $\hat{s}$ and the incident direction $\hat{s}_i$. The scattering phase function was predicted on a spectral basis for *Chlorella vulgaris* using the Lorenz–Mie theory for a coated sphere according to the method described in [35]. The resulting spectral backward scattering ratio $b_{\lambda}$ is shown in Figure A1 of Appendix A.

For PBRs with droplet-covered windows, the solar radiation was transmitted through the window into many directions as depicted in Figure 2b. Thus, for a given solar position ($\theta_z, \gamma_s$), the spectral local fluence rate $G_{\lambda}(z)$ inside the culture covered by a droplet-covered window was found by summing the contribution from all transmission directions $\theta_{t,j}$ according to:

$$ G_{\lambda}(z) = \sum_{j=1}^{M+1} G_{\lambda,j}(z, \theta_{t,j}). \quad (11) $$

For the dry window PBR, the light reaching the upper boundary of the microalgae culture was collimated, and $G_{\lambda}(z)$ is given by Equation (8) for that specific direction.

The ability of the microalgae to absorb the local spectral fluence rate $G_{\lambda}(z)$ was then described by the local rate of photon absorption (LRPA) $A(z)$ (in $\mu mol_{hv} kg^{-1} s^{-1}$) defined as [9]:

$$ A(z) = \int_{PAR} \bar{A}_{abs,\lambda} G_{\lambda}(z) d\lambda. \quad (12) $$

### 3.5. Microalgae Growth Kinetics

The microalgae growth kinetic model proposed by Takache et al. [27] and further developed by Souliès et al. [9] for *Chlorella vulgaris* was employed to predict the average growth rate $\bar{\mu}$ of the microalgae culture. As discussed previously, this model accounted for light limitation, photoinhibition, and cell respiration activity as a function of the LRPA $A(z)$ [27]. First, the local specific rate of oxygen production or consumption $J_{O_2}(z)$ (in $mol_{O_2} kg^{-1} s^{-1}$) was calculated according to [9]:

$$ J_{O_2}(z) = \left[ \rho_M \frac{K}{K + A(z)} \Phi'_{O_2} A(z) - \frac{J_{NADH_2}}{\nu_{NADH_2-O_2}} \frac{K_r}{K_r + A(z)} \right]. \quad (13) $$
here, $\rho_M$ is the maximum energy yield for photon conversion, $\bar{\phi}'_{O_2}$ (in mol$_{O_2}$ µmol$^{-1}$) is the mole quantum yield of $O_2$ for the Z-scheme of photosynthesis, $K$ (in µmol$_{hv}$ kg$^{-1}$ s$^{-1}$) is the half-saturation constant for photosynthesis, $J_{NADH_2}$ (in mol$_{NADH_2}$ kg$^{-1}$ s$^{-1}$) is the specific rate of cofactor regeneration on the respiratory chain, related to the oxygen consumption by the stoichiometric coefficient of cofactor regeneration on the respiratory chain $v_{NADH_2-O_2}$, and $K_r$ (in µmol$_{hv}$ kg$^{-1}$ s$^{-1}$) is the half saturation constant describing the inhibition of respiration in light. These parameters were taken from [9] and are summarized in Table A1. These parameters and the absorption $\bar{A}_{abs,\lambda}$ and scattering $\bar{S}_{sca,\lambda}$ cross-sections corresponded to the same culture [9]. The average specific rate of oxygen production $\bar{J}_{O_2}$ over the depth $L$ of the PBR was defined as [9]:

$$ \bar{J}_{O_2} = \frac{1}{L} \int_{0}^{L} J_{O_2}(z) dz. \qquad (14) $$

Then, $\bar{J}_{O_2}$ was used to predict the average specific growth rate $\bar{\mu}$ (in s$^{-1}$) via the stoichiometric relationship between the production of oxygen and the production of biomass [9]:

$$ \bar{\mu} = \frac{\bar{J}_{O_2} M_x}{v_{O_2-X}} \qquad (15) $$

where $M_x$ (in kg$_x$ mol$^{-1}$) is the C-molar mass in the biomass given by $C_m H_p O_n$ and $v_{O_2-X}$ is the stoichiometric coefficient of the oxygen production.

Finally, the temporal evolution of the biomass concentration $X(t)$ was found by integrating Equation (1). Then, to assess the impact of droplets on PBR performance, the areal biomass productivity for a given day $P$ (in kg m$^{-2}$ day$^{-1}$) was expressed as:

$$ P = (X_f - X_0)L \qquad (16) $$

where $X_0$ and $X_f$ are the initial and final biomass concentrations, respectively, given by $X_0 = X(t_0 = 7 \text{ am})$ and $X_f = X(t_f = 7 \text{ pm})$ on 21 June and $X_0 = X(t_0 = 9 \text{ am})$ and $X_f = X(t_f = 5 \text{ pm})$ on 23 September.

Another important PBR performance metric that can be used to assess the effect of droplets is the so-called illuminated fraction $\gamma$ given by [25]:

$$ \gamma = \frac{z_c}{L} \qquad (17) $$

where $z_c$ is the maximum culture depth at which the local LRPA $A(z)$ still yields positive photosynthetic growth, i.e., $A(z_c) = A_c$ where $A_c$ is the minimum LRPA required to prevent biomass loss due to respiration [9]. When $\gamma = 1$, positive photosynthetic growth occurs throughout the entire culture, and all incident light available for photosynthesis is absorbed by the culture. Under these conditions, the maximum average growth rate $\bar{\mu}_{max}$ is reached for given light conditions [25]. When $\gamma < 1$, dark zones are present in the culture volume, and biomass loss occurs due to respiration. By contrast, $\gamma > 1$ indicates that excess light is available for photosynthesis.

### 3.6. Initial and Boundary Conditions

The solar spectrum incident upon the PBR was determined by the PBR location and the time of year. Table 1 shows the value of the solar zenith $\theta_z$ and azimuth $\gamma_s$ angles associated with each simulated time on 21 June and 23 September. Figure 3a,b shows the spectral solar direct normal irradiance $G_{S,\lambda}$ (in µmol$_{hv}$ m$^{-2}$ s$^{-1}$) in Los Angeles predicted by the Simple Model of the Atmospheric Radiative Transfer of Sunshine (SMARTS) [36] for wavelengths in the PAR region from 400 to 700 nm between 7 a.m. and 7 p.m. on 21 June and 9 a.m. and 5 p.m. on 23 September, respectively. As mentioned previously, edge effects were neglected, and the bottom of the PBR was considered to be black. The initial biomass concentration $X_0 = X(t_0)$ was varied between 0.01 and 1.0 g L$^{-1}$.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Incident collimated solar irradiance $G_{S,\lambda}$ over the PAR region at different times of day simulated on (a) 21 June and (b) 23 September in Los Angeles, CA.</figcaption></figure>

<table><caption>**Table 1.** Solar zenith $\theta_z$ and azimuth $\gamma_s$ angles (in degrees) for the different times simulated on 21 June 21 and 23 September 23 in Los Angeles, CA.</caption><thead><tr><th rowspan="2">Time</th><th colspan="2">21 June</th><th colspan="2">23 September</th></tr><tr><th>θ<sub>z</sub></th><th>γ<sub>s</sub></th><th>θ<sub>z</sub></th><th>γ<sub>s</sub></th></tr></thead><tbody><tr><td>7 a.m.</td><td>76.2</td><td>-109.2</td><td>-</td><td>-</td></tr><tr><td>9 a.m.</td><td>52.0</td><td>-94.7</td><td>62.9</td><td>-69.6</td></tr><tr><td>11 a.m.</td><td>27.3</td><td>-74.8</td><td>42.3</td><td>-41.6</td></tr><tr><td>1 p.m.</td><td>10.7</td><td>5.5</td><td>34.4</td><td>6.2</td></tr><tr><td>3 p.m.</td><td>29.1</td><td>76.9</td><td>46.5</td><td>49.5</td></tr><tr><td>5 p.m.</td><td>53.8</td><td>95.8</td><td>68.4</td><td>74.2</td></tr><tr><td>7 p.m.</td><td>77.9</td><td>110.3</td><td>-</td><td>-</td></tr></tbody></table>

### 3.7. Method of Solution

Figure 4 shows a block diagram of the procedure used to predict the final biomass concentration $X_f$ and the PBR daily areal biomass productivity $P$. The droplet parameters ($\theta_c, f_A$), the spectral solar irradiance $G_{S,\lambda}$, and the solar incidence angles ($\theta_z, \gamma_s$) at initial time $t = t_0$ were the input parameters for modeling light transfer through the PBR window. The appropriate light transfer model was then applied according to the window condition. Then, the resulting spectral radiative flux at the upper boundary of the microalgae culture $q_{m,\lambda}(t)$ and the PBR operating parameters ($X_0, L$) were used to model light transfer and growth kinetics in the culture to predict the biomass concentration $X(t)$. This process was then repeated for subsequent time steps until the final biomass concentration $X_f = X(t_f)$ and the daily areal biomass productivity $P$ were computed.
<figure><img src="image_7.png" /><figcaption>Figure 4. Block diagram of the procedure used to predict the final daily biomass $X_f$ and daily areal biomass productivity $P$ for PBRs with dry and droplet-covered windows.</figcaption></figure>

Light transfer through the dry window and into the microalgae culture was modeled using Equations (2) and (3). For a droplet-covered window, the Monte Carlo ray-tracing method was used to predict the bidirectional transmittance $T_{bd}(\theta_z, \gamma_s, \theta_t, \gamma_t)$. Here, a large number of discrete collimated photon bundles or “rays” were launched in direction $(\theta_z, \gamma_s)$ towards the window covered by computationally-generated droplets. The method and algorithm were described in detail in [7,8,12] and need not be repeated. In brief, the polar transmission angle $\theta_t$ ranged from 0° to 90° and was discretized into $M+1$ angular intervals where $M = 90$, with $\Delta\theta_{t,j=1} = \Delta\theta_{t,j=M+1} = \pi/4M$ (in rad) for the first and last intervals at angles $\theta_{t,1} = \pi/8M$ (in rad) and $\theta_{t,M+1} = \pi/2 - \pi/8M$ (in rad) while the interval $(\Delta\theta_{t,j})_{2\le j\le M} = \pi/2M$ (in rad) was used for all other transmission angles. The azimuthal transmission angle $\varphi_t$ ranged from 0° to 360° in uniform intervals of $\Delta\varphi_t = \pi/180$ (in rad). At each interface a ray encountered, the interface reflectivity was calculated using Fresnel’s equations. Then, a random number selected from the uniform distribution was compared to the interface reflectivity to determine if the ray was reflected or refracted. The new ray direction was governed by specular reflection or Snell’s law for refracted rays. The next interface the ray encountered was then determined using the geometry of the 3D simulation domain. For a non-absorbing window and droplets, this process continued until the ray was either (i) transmitted through the window and droplets or (ii) reflected away from the window. Then, the fraction of solar energy transmitted through the droplet-covered window for a given solar position $(\theta_z, \gamma_s)$ was represented by the directional-hemispherical transmittance $T_{dh}$ expressed as [8]:

$$T_{dh}(\theta_z, \gamma_s) = \frac{N_t}{N_i} \qquad (18)$$

where $N_i$ is the total number of rays incident on the window and $N_t$ is the total number of transmitted rays. Furthermore, Equation (4) was recast to express the bidirectional transmittance for a given solar position $(\theta_z, \gamma_s)$ and transmitted direction $(\theta_{t,j}, \varphi_{t,j})$ from the output of the Monte Carlo ray-tracing simulation according to [7]:

$$T_{bd,jk}(\theta_z, \gamma_s, \theta_{t,j}, \varphi_{t,j}) = \frac{N_{t,jk}}{N_i \cos \theta_t \Delta\Omega_{t,jk}} \qquad (19)$$
where $N_{t,jk}$ is the number of rays transmitted into the solid angle $\Delta\Omega_{t,jk} = \sin\theta_{t,j}\Delta\theta_{t,j}\Delta\phi_{t,k}$. In order to achieve numerical convergence for $T_{bd,jk}$, $N_i = 10^7$ rays were simulated [7]. Then, the resulting spectral radiative flux at the upper boundary of microalgae culture $q_{m,\lambda,j}$ was calculated using the bidirectional transmittance $T_{bd,jk}$ and Equations (5)–(7).

Next, the LRPA $\mathcal{A}(z, t)$ accounting for all transmitted directions was predicted via the two-flux approximation according to Equations (8)–(12). Then, the growth kinetics model was used to predict the average specific growth rate $\bar{\mu}(t)$ based on Equations (13)–(15). The biomass concentration at time $t + \Delta t$ was predicted from the biomass concentration $X(t)$ at time $t$ by integrating Equation (1) according to [22]:

$$X(t + \Delta t) = X(t)[1 + \bar{\mu}(t)\Delta t]. \qquad (20)$$

Light transfer through the droplet-covered window was simulated in time increments of $\Delta t = 2$ h from 7 a.m. to 7 p.m. on 21 June and 9 a.m. to 5 p.m. on 23 September. To check for numerical convergence, $X(t)$ and $\bar{\mu}(t)$ were predicted by simulating light transfer through a dry window PBR using a time increment of $\Delta t = 0.5$ h and $\Delta t = 2$ h (see Figure A2). The maximum relative errors throughout the day for the biomass $X(t)$ were less than 0.1% and 0.3% on 21 June and 23 September, respectively. Thus, a time increment of 2 h was considered to be an appropriate compromise between reasonable computation time and sufficient accuracy.

# 4. Results and Discussion

## 4.1. Light Transfer through Droplet Covered PBR Windows

Figure 5a,b show the one-dimensional bidirectional transmittance $\bar{T}_{bd}(\theta_z, \gamma_s, \theta_t)$ as a function of the transmission angle $\theta_t$ at 1 p.m. and 5 p.m. for various window conditions on 21 June and 23 September, respectively. On 21 June at 1 p.m., all window conditions showed a strong peak in $\bar{T}_{bd}$ at the transmission angle equal to the solar zenith angle (i.e., $\theta_t = \theta_z$) with a large surrounding lobe. At 5 pm, as the solar zenith angle $\theta_z$ increased, the corresponding peak in $\bar{T}_{bd}$ decreased, and the surrounding lobes became wider for a given droplet configuration. This was indicative of increased scattering at larger solar zenith angles. Indeed, for a window with droplets of a contact angle of $\theta_c = 90^\circ$ and surface area coverage $f_A = 90\%$ on 23 September at 5 p.m., so much scattering occurred that $\bar{T}_{bd}$ peaked around $\theta_t = 45^\circ$ despite a solar zenith angle $\theta_z$ of 68$^\circ$ at this time.

Figure 5c,d show the temporal evolution of the average transmission angle $\bar{\theta}_t$ defined as the transmission angle $\theta_t$ weighted by the one-dimensional bidirectional transmittance $\bar{T}_{bd}(\theta_z, \gamma_s, \theta_t)$ for a given solar position $(\theta_z, \gamma_s)$, i.e.,

$$\bar{\theta}_t = 2\pi \int_0^{\pi/2} \bar{T}_{bd}(\theta_z, \gamma_s, \theta_t) \theta_t \sin \theta_t d\theta_t. \qquad (21)$$

In the case of a dry window, no scattering occurred, and the transmission angle was equal to the incident angle, i.e., $\theta_t = \theta_z$. However, for condensate-covered windows, at all times except between 11 a.m. and 3 p.m. on 21 June, light scattering by the droplets reduced the average transmission angle by as much as 27$^\circ$ and 20$^\circ$ on 21 June and 23 September, respectively. In other words, for large solar zenith angles $\theta_z$, the droplets reduced the transmittance, but scattered sunlight in directions closer to normal incidence, causing the transmitted light to penetrate deeper into the culture compared to a culture covered by a dry window.
<figure><img src="image_8.png" /><figcaption>(b) September 23rd</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>(c) June 21st</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>(d) September 23rd</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>(e) June 21st</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>(f) September 23rd</figcaption></figure>

<figure><img src="image_13.png" /></figure>

**Figure 5.** One-dimensional bidirectional transmittance $\bar{T}_{bd}$ at 1 p.m. and 5 p.m. (a,b), average transmission angle $\bar{\theta}_t$ (c,d), and directional-hemispherical transmittance $T_{dh}(\theta_z, \gamma_s)$ (e,f) as functions of time for a variety of window conditions on (a,c,e) 21 June and (b,d,f) 23 September.
Finally, Figure 5e,f show the temporal evolution of the directional-hemispherical transmittance $T_{dh}$ for various window conditions on 21 June and 23 September, respectively. On both days, Figure 5e,f indicate that the directional-hemispherical transmittance $T_{dh}$ for the dry window increased to reach a maximum around midday, before decreasing in the afternoon. This was due to the larger reflectance of glass at larger incidence angles [33]. Furthermore, Figure 5e,f show that the presence of droplets reduced $T_{dh}$ for nearly all solar positions investigated on both days. Indeed, droplets with contact angle $\theta_c = 90^\circ$ and surface area coverage $f_A = 90\%$ reduced $T_{dh}$ by up to 37\% on 21 June. However, on the same day, $T_{dh}$ increased slightly compared with a dry window for droplet contact angle $\theta_c = 30^\circ$ and surface area coverage $f_A = 90\%$ at 11 a.m. and 1 p.m. and for $\theta_c = 90^\circ$ and $f_A = 90\%$ at 7 a.m. and 7 p.m. This was due to the smaller refractive index mismatch between the window ($n_w = 1.5$) and the water droplets ($n_d = 1.33$) for droplet-covered windows than between the window ($n_w = 1.5$) and the air ($n_a = 1.0$) for dry windows. Then, the presence of droplets reduced reflection at the back side of the window, resulting in slightly larger $T_{dh}$ compared to the dry window case. In addition, Figure 5 indicates that for droplet contact angle $\theta_c = 90^\circ$, increasing the droplet surface area coverage $f_A$ resulted in smaller transmittance $T_{dh}$ of the PBR window between 9 a.m. and 5 p.m. on 21 June and at all times except 5 p.m. on 23 September. Moreover, under normal incidence, previous studies have shown that for a given surface area coverage $f_A$, the normal-hemispherical transmittance $T_{nh}$ decreased significantly with increasing droplet contact angle $\theta_c$ greater than the critical angle $\theta_{cr}$ for total internal reflection at the droplet/air interface such that $\theta_{cr} = \sin^{-1}(1/n_d) = 48.8^\circ$ [8]. This phenomenon was also observed in Figure 5e,f at near-normal incidence angles (i.e., small solar zenith angles $\theta_z$).

## 4.2. Light Transfer in Microalgae Culture

Figure 6a,b show the normalized LRPA $A/A_c$ as a function of the normalized culture depth $z/L$ for various window conditions at 1 p.m. and 5 p.m. in Los Angeles, CA, on 21 June and 23 September, respectively. Here, the initial biomass concentration of the PBR was $X_0 = 0.07 \text{ g L}^{-1}$, and the culture depth was $L = 0.3 \text{ m}$. For both days, Figure 6 indicates that the normalized LRPA $A/A_c$ exceeded 1.0 at 1 p.m. throughout the culture while at 5 p.m. up to 58\% of the culture volume features $A/A_c < 1.0$, i.e., an illuminated fraction $\gamma$ as low as 0.42. This was caused by (i) the smaller solar irradiance later in the day, (ii) the increase in the pathlength of light through the culture for larger solar zenith angles $\theta_z$, and (iii) the photosynthetic response that occurred between 1 p.m. and 5 p.m. Thus, more of the incident radiation was absorbed at shallower culture depths at 5 p.m. than at 1 p.m. These observations were also reported in previous studies [9,25]. On both days at 1 p.m., the droplets were responsible for a decrease in the local rate of photon absorption throughout the PBR depth, with the exception of a window with $f_A = 90\%$ and $\theta_c = 30^\circ$ on 21 June at 1 p.m. due to the increase in $T_{dh}$ observed in Figure 5e. On the other hand, on both days at 5 p.m., Figure 6 shows that the presence of droplets with contact angle $\theta_c = 90^\circ$ resulted in an increase of up to a 14\% in the illuminated fraction $\gamma$ of the PBR, compared to a dry window. This was due to light scattering by the droplets observed in Figure 5a–d, which redistributed the incoming light in many directions and caused light to penetrate deeper into the culture.
<figure><img src="image_14.png" /><figcaption>**Figure 6.** Normalized local rate of photon absorption (LRPA) $A/A_c$ as a function of the normalized culture depth $z/L$ at 1 p.m. and 5 p.m. for a variety of window conditions on (a) 21 June and (b) 23 September in Los Angeles, CA. Here, $L = 0.3$ and $A_c = 2800 \ \mu\text{mol}_{\text{h}\nu} \text{kg}^{-1} \text{s}^{-1}$.</figcaption></figure>

### 4.3. Microalgae Growth Kinetics

Figure 7a,b plot the temporal evolution of the average growth rate $\bar{\mu}(t)$ (in h$^{-1}$) for PBRs with various window conditions on 21 June and 23 September, respectively. Here also, the initial biomass concentration $X_0$ was 0.07 g L$^{-1}$, and the culture depth $L$ was 0.3 m. For most of the day on 21 June, Figure 7a indicates that droplets reduced the average growth rate $\bar{\mu}(t)$. Furthermore, the average growth rate $\bar{\mu}(t)$ decreased with increasing droplet contact angle $\theta_c$ and surface area coverage $f_A$. A similar trend was observed on 23 September. However, at 5 p.m. on 23 September and 7 a.m. and 7 p.m. on 21 June, droplets with $\theta_c = 90^\circ$ and $f_A = 90\%$ slightly increased $\bar{\mu}(t)$. This was attributed to the positive effect of light scattering by droplets towards the culture observed previously, which increased the illuminated fraction $\gamma$ and light penetration in the culture late in the day (see Figure 6). Interestingly, on 23 September, this was the case despite a relative decrease in the directional-hemispherical transmittance $T_{dh}$ of up to 24% compared to a dry window. Nonetheless, this phenomena occurred in the morning and evening, when the average growth rate $\bar{\mu}(t)$ was a fraction of its daily peak value and the impact on the final biomass concentration $X_f$ was minimal.

Indeed, Figure 7c,d plot the temporal evolution of the biomass concentration $X(t)$ as a function of time for various window conditions on 21 June and 23 September, respectively. The figures indicate that the final biomass concentration $X_f$ was smaller on 23 September compared to 21 June for all configurations considered. This was due to the decrease in spectral solar irradiance $G_{S,\lambda}$ on 23 September (see Figure 3a,b) and the larger solar zenith angles $\theta_z$ (see Table 1), which led to a smaller window transmittance $T_{dh}$ in the presence of droplets. The presence of droplets decreased the final biomass concentration $X_f$ by up to 10.0% and 8.2% on 21 June and 23 September, respectively.
<figure><img src="image_15.png" /><figcaption>**Figure 7.** Average PBR growth rate $\bar{\mu}$ as a function of time for various window conditions on (a) 21 June and (b) 23 September and corresponding biomass concentration $X(t)$ on (c) 21 June and (d) 23 September.</figcaption></figure>

4.4. PBR Biomass Productivity

Figure 8a,b show the daily areal biomass productivity $P$ (in kg m$^{-2}$ day$^{-1}$) as a function of the initial biomass concentration $X_0$ for various culture depths $L$ and window conditions for PBR operation between 7 a.m. and 7 p.m. on 21 June and 9 a.m. and 5 p.m. on 23 September, respectively. The trends in the daily biomass productivity $P$ were similar for both days although $P$ was smaller on 23 September than on 21 June. On both days, for a given culture depth $L$, the culture did not fully absorb the incoming radiation for small initial biomass concentrations $X_0$. However, as $X_0$ increased, more of the incoming light was absorbed by the culture, and $P$ increased to reach a maximum $P_{max}$. As previously discussed, the optimal biomass productivity occurs for a light-limited PBR when the incoming radiation is fully absorbed in the culture, i.e., the microalgae suspension is fully illuminated with no dark zones. As $X_0$ increased further, dark zones appeared in the culture causing the areal biomass productivity $P$ to decrease. Figure 8a,b indicate that as the culture depth $L$ increased, the maximum biomass productivity occurred for smaller values of $X_0$.
Lee et al. [22] demonstrated that PBR productivity *P* scaled with *X*₀/*a* where *a* is the specific illuminated area. In the case of the simulated raceway pond, *a* is given by *a* = *S*/*V* = 1/*L* where *S* is the illuminated surface area and *V* is the culture volume [30]. Furthermore, we note that the initial optical thickness of the PBR is defined as *β*<sub>λ,0</sub>*L* = (Â<sub>abs,λ</sub> + õ<sub>sca,λ</sub>)*X*₀*L* where *β*<sub>λ,0</sub> is the initial culture extinction coefficient (in m<sup>-1</sup>). Thus, the quantity *X*₀*L* = *X*₀/*a* is representative of the culture initial optical thickness [22]. In order to test this scaling relationship for PBRs covered by windows supporting pendant droplets, Figure 8c,d plots the biomass productivity *P* with respect to the initial optical thickness represented by *X*₀/*a* for various window conditions. Both figures indicate that, for a given day and window condition, all data collapse on the same line. Thus, the biomass productivity *P* depended solely on *X*₀/*a*, the PBR window conditions, and the day of the year. Furthermore, for a given day of the year, the optimum initial optical thickness (*X*₀/*a*)<sub>opt</sub> corresponding to the maximum biomass productivity *P*<sub>max</sub> was essentially independent of the window condition. Note that the same results were obtained by Lee et al. [22], albeit using a different growth kinetics model for open ponds, as well as vertical and tubular PBRs without droplets cultivating *Chlamydomonas reinhardtii*. These results confirm the importance of the optical thickness represented by *X*₀/*a* for maximizing PBR biomass productivity through the PBR design (via *a* = 1/*L*) and its operation (via *X*₀) regardless of the window condition and/or growth kinetics model. Figure 8d shows that on 23 September, all droplet configurations decreased the maximum biomass productivity *P*<sub>max</sub>. By contrast, Figure 8c shows that, on 21 June, the effect of droplets with *θ*<sub>c</sub> = 30° and *f*<sub>A</sub> = 90% was essentially negligible, while droplets with *θ*<sub>c</sub> = 90° decreased *P*<sub>max</sub>.

Table 2 summarizes the maximum biomass productivity $P_{max}$ and the corresponding optimum initial optical thickness ($X_0/a)_{opt}$ for each droplet configuration for both 21 June and 23 September, as well as the change in $P_{max}$ relative to the reference case of a PBR covered by a dry window. It is evident that the droplet contact angle strongly affected the maximum achievable PBR biomass productivity $P_{max}$. Indeed, the maximum biomass productivity $P_{max}$ decreased by just 2.0% and 5.1% for droplet contact angle $\theta_c = 30^\circ$ and surface area coverage $f_A = 90\%$ on 21 June and 23 September, respectively. On the other hand, $P_{max}$ decreased by 17.5% and 18.1% for $\theta_c = 90^\circ$ and $f_A = 90\%$ on 21 June and 23 September, respectively. Although hydrophobic surfaces may inhibit condensation [37], these results suggest that a hydrophilic cover featuring small droplet contact angle $\theta_c$ should be selected whenever possible to mitigate the effect of droplets on the light reaching the culture and on the PBR biomass productivity. Indeed, on a given day, the daily biomass productivity $P$ was larger for droplet contact angle $\theta_c = 30^\circ$ and surface area coverage $f_A = 90\%$ than for $\theta_c = 90^\circ$ and $f_A = 50\%$ on both days and for all values of optical thickness $X_0/a$.
<figure><img src="image_16.png" /><figcaption>**Figure 8.** Daily areal PBR biomass productivity *P* as a function of (a,b) starting biomass concentration *X*₀ for various PBR depths *L* and window conditions and (c,d) the same daily areal PBR biomass productivity *P* as a function of the microalgae culture initial optical thickness represented by *X*₀/*a*, on (a,c) 21 June and (b,d) 23 September</figcaption></figure>
<table><caption>**Table 2.** Maximum biomass productivity $P_{max}$ and the corresponding optimum optical thickness $(X_0/a)_{opt}$ for each window condition on 21 June and 23 September. The percent change in $P_{max}$ compared to that of a dry window ($f_A = 0\%$) is also shown.</caption><thead><tr><th colspan="5">21 June</th></tr><tr><th>f<sub>A</sub></th><th>θ<sub>c</sub></th><th>(X<sub>0</sub>/a)<sub>opt</sub><br>(kg m<sup>-2</sup>)</th><th>P<sub>max</sub><br>(kg m<sup>-2</sup> day<sup>-1</sup>)</th><th>P<sub>max</sub><br>Change</th></tr></thead><tbody><tr><td>0%</td><td>N/A</td><td>0.018</td><td>0.0245</td><td>0%</td></tr><tr><td>90%</td><td>30°</td><td>0.018</td><td>0.0240</td><td>-2.0%</td></tr><tr><td>50%</td><td>90°</td><td>0.018</td><td>0.0219</td><td>-10.6%</td></tr><tr><td>90%</td><td>90°</td><td>0.018</td><td>0.0202</td><td>-17.5%</td></tr></tbody></table><table><thead><tr><th colspan="5">23 September</th></tr><tr><th>f<sub>A</sub></th><th>θ<sub>c</sub></th><th>(X<sub>0</sub>/a)<sub>opt</sub><br>(kg m<sup>-2</sup>)</th><th>P<sub>max</sub><br>(kg m<sup>-2</sup> day<sup>-1</sup>)</th><th>P<sub>max</sub><br>Change</th></tr></thead><tbody><tr><td>0%</td><td>N/A</td><td>0.021</td><td>0.0177</td><td>0%</td></tr><tr><td>90%</td><td>30°</td><td>0.021</td><td>0.0168</td><td>-5.1%</td></tr><tr><td>50%</td><td>90°</td><td>0.021</td><td>0.0156</td><td>-11.9%</td></tr><tr><td>90%</td><td>90°</td><td>0.021</td><td>0.0145</td><td>-18.1%</td></tr></tbody></table>

## 5. Conclusions

This study demonstrates the effect of condensate droplets on the biomass productivity of covered outdoor raceway ponds. Numerical simulations of light transfer were coupled to a microalgae growth kinetics to predict the performance of PBRs with and without condensate droplet-covered windows. The bidirectional transmittance of droplet-covered windows in Los Angeles, CA, at different times throughout the day on 21 June and 23 September was predicted by the Monte Carlo ray-tracing method. The local fluence rate within the horizontal culture of *Chlorella vulgaris* was predicted using the two-flux approximation for oblique incidence. A growth kinetics model accounting for photolimitation, photoinhibition, and respiration was applied. Droplets with contact angle $\theta_c = 90°$ decreased the maximum areal PBR biomass productivity $P_{max}$ by up to 17.5% and 18.1% on 21 June and 23 September, respectively. On both days, the maximum daily areal PBR biomass productivity decreased with increasing droplet contact angle and surface area coverage compared with a dry window. Thus, hydrophilic PBR covers are preferable to minimize the adverse impacts of droplets on PBR biomass productivity. Additionally, the areal PBR biomass productivity scaled with the ratio $X_0/a$ of the initial biomass concentration $X_0$ and the specific illuminated area $a = 1/L$ such that $X_0/a = X_0L$ was proportional to the initial optical thickness of the PBR. Finally, for a given day of the year, the maximum biomass productivity occurred at the same initial optical thickness for all window conditions simulated. This optimal initial optical thickness can be used to maximize PBR biomass productivity through PBR design (via $L$) and operation (via $X_0$).

**Author Contributions:** Conceptualization, J.P., J.L., E.-k.S.-A., and L.P.; data curation, J.H., K.Z. and J.P.; formal analysis, J.H. and K.Z.; funding acquisition, L.P.; investigation, J.H. and K.Z.; methodology, J.H., K.Z., J.P., and L.P.; software, J.H., K.Z., and J.P.; supervision, J.P., J.L., and L.P.; validation, J.H. and K.Z.; visualization, J.H. and L.P.; writing, original draft, J.H.; writing, review and editing, J.H., J.P., J.L., E.-k.S.-A., and L.P. All authors read and agreed to the published version of the manuscript.

**Funding:** This research was supported in part by the National Science Foundation NRT-INFEWS: Integrated Urban Solutions for Food, Energy, and Water Management (Grant No. DGE-1735325) and the DISCUS project of the NExT Initiative (Nantes Excellence Trajectory) International Research Partnership.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.
**Data Availability Statement:** Data presented in this manuscript are available upon request from the
corresponding author.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; nor in the decision to publish the results.

**Nomenclature**

<table><tr><td>a</td><td>specific illuminated area, m<sup>-1</sup></td></tr><tr><td>A&#773;<sub>abs,&lambda;</sub></td><td>average spectral mass absorption cross-section, m<sup>2</sup> kg<sup>-1</sup></td></tr><tr><td>A(z)</td><td>local specific rate of photon absorption, &mu;mol<sub>h&nu;</sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>A<sub>c</sub></td><td>specific rate of photon absorption at the compensation point, &mu;mol<sub>h&nu;</sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>b<sub>&lambda;</sub></td><td>spectral back-scattering ratio</td></tr><tr><td>d</td><td>droplet diameter, &mu;m</td></tr><tr><td>d<sub>p</sub></td><td>droplet projected diameter, &mu;m</td></tr><tr><td>f<sub>A</sub></td><td>droplet surface area coverage, %</td></tr><tr><td>G</td><td>fluence rate, &mu;mol<sub>h&nu;</sub> m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>H</td><td>window thickness, m</td></tr><tr><td>I</td><td>radiative intensity, &mu;mol<sub>h&nu;</sub> m<sup>-2</sup> s<sup>-1</sup> sr<sup>-1</sup></td></tr><tr><td>J<sub>NADH<sub>2</sub></sub></td><td>specific rate of cofactor regeneration, mol<sub>NADH<sub>2</sub></sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>J<sub>O<sub>2</sub></sub>(z)</td><td>local specific rate of oxygen production/consumption, in mol<sub>O<sub>2</sub></sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>J&#773;<sub>O<sub>2</sub></sub></td><td>average specific rate of oxygen production/consumption, mol<sub>O<sub>2</sub></sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>k</td><td>absorption index</td></tr><tr><td>K</td><td>half saturation constant for photosynthesis, &mu;mol<sub>h&nu;</sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>K<sub>r</sub></td><td>half saturation constant for respiration, &mu;mol<sub>h&nu;</sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>L</td><td>microalgae culture depth, m</td></tr><tr><td>M</td><td>interval number of transmission angles</td></tr><tr><td>M<sub>x</sub></td><td>carbon molar mass in the biomass, kg<sub>x</sub> mol<sup>-1</sup></td></tr><tr><td>N</td><td>number of rays</td></tr><tr><td>n</td><td>refractive index</td></tr><tr><td>P</td><td>daily PBR biomass productivity, kg m<sup>-2</sup> day<sup>-1</sup></td></tr><tr><td>q</td><td>radiative flux, &mu;mol<sub>h&nu;</sub> m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>S</td><td>culture illuminated surface area, m<sup>2</sup></td></tr><tr><td>S&#773;<sub>sca,&lambda;</sub></td><td>average spectral mass scattering cross-section, m<sup>2</sup> kg<sup>-1</sup></td></tr><tr><td>t</td><td>time, hour</td></tr><tr><td>T<sub>bd</sub></td><td>bidirectional transmittance, sr<sup>-1</sup></td></tr><tr><td>T&#773;<sub>bd</sub></td><td>one-dimensional bidirectional transmittance, sr<sup>-1</sup></td></tr><tr><td>T<sub>dh</sub></td><td>directional-hemispherical transmittance</td></tr><tr><td>T<sub>nh</sub></td><td>normal-hemispherical transmittance</td></tr><tr><td>V</td><td>culture volume, m<sup>3</sup></td></tr><tr><td>X</td><td>biomass concentration, g L<sup>-1</sup></td></tr><tr><td>z</td><td>culture depth, m</td></tr></table>
## Greak Symbols

<table><tr><td>α<sub>λ</sub></td><td>spectral linear scattering modulus for the two-flux approximation</td></tr><tr><td>β<sub>λ</sub></td><td>spectral extinction coefficient m<sup>-1</sup></td></tr><tr><td>γ</td><td>PBR illuminated fraction</td></tr><tr><td>γ<sub>s</sub></td><td>solar azimuth angle, °</td></tr><tr><td>δ<sub>λ</sub></td><td>spectral extinction coefficient for the two flux approximation, m<sup>-1</sup></td></tr><tr><td>Θ</td><td>scattering angle, °</td></tr><tr><td>θ</td><td>polar angle, °</td></tr><tr><td>θ<sub>c</sub></td><td>droplet contact angle, °</td></tr><tr><td>θ<sub>cr</sub></td><td>critical angle for total internal reflection, °</td></tr><tr><td>θ&#773;<sub>t</sub></td><td>average transmission angle, °</td></tr><tr><td>θ<sub>z</sub></td><td>solar zenith angle, °</td></tr><tr><td>λ</td><td>light wavelength, nm</td></tr><tr><td>μ&#773;</td><td>average specific growth rate, h<sup>-1</sup></td></tr><tr><td>ν<sub>NADH<sub>2</sub>-O<sub>2</sub></sub></td><td>stoichiometric coefficient of cofactor regeneration</td></tr><tr><td>ν<sub>O<sub>2</sub>-X</sub></td><td>stoichiometric coefficient of oxygen production</td></tr><tr><td>ρ<sub>M</sub></td><td>maximum energy yield for photon conversion</td></tr><tr><td>τ</td><td>transmissivity</td></tr><tr><td>Φ</td><td>scattering phase function</td></tr><tr><td>φ&#773;<sub>O<sub>2</sub></sub></td><td>mean mass quantum yield for the Z-scheme, mol<sub>O<sub>2</sub></sub> µmol<sub>hν</sub><sup>-1</sup></td></tr><tr><td>φ</td><td>azimuth angle, °</td></tr><tr><td>Ω</td><td>solid angle, sr</td></tr></table>

## Subscripts

<table><tr><td>0</td><td>refers to initial conditions</td></tr><tr><td>a</td><td>refers to air</td></tr><tr><td>d</td><td>refers to droplet</td></tr><tr><td>f</td><td>refers to final conditions</td></tr><tr><td>i</td><td>refers to incidence</td></tr><tr><td>j</td><td>index for transmitted polar angle θ<sub>t</sub></td></tr><tr><td>k</td><td>index for transmitted azimuthal angle φ<sub>t</sub></td></tr><tr><td>m</td><td>refers to microalgae culture</td></tr><tr><td>s</td><td>refers to solar variable</td></tr><tr><td>t</td><td>refers to transmission</td></tr><tr><td>w</td><td>refers to window</td></tr><tr><td>X</td><td>refers to biomass</td></tr><tr><td>λ</td><td>refers to a spectral variable</td></tr></table>

# Appendix A

<table><caption>Table A1. Kinetic growth parameters of *Chlorella vulgaris* from Ref. [9].</caption><thead><tr><th>Parameter</th><th>Value</th><th>Units</th></tr></thead><tbody><tr><td>ρ<sub>M</sub></td><td>0.8</td><td>-</td></tr><tr><td>J<sub>NADH<sub>2</sub></sub></td><td>2.8×10<sup>-3</sup></td><td>mol<sub>NADH<sub>2</sub></sub> kg<sub>X</sub><sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>ν<sub>O<sub>2</sub>-X</sub></td><td>1.13</td><td>-</td></tr><tr><td>φ&#773;<sup>′</sup><sub>O<sub>2</sub></sub></td><td>1.1×10<sup>-7</sup></td><td>mol<sub>O<sub>2</sub></sub> µmol<sub>hν</sub><sup>-1</sup></td></tr><tr><td>M<sub>X</sub></td><td>0.024</td><td>kg<sub>X</sub> mol<sub>C</sub><sup>-1</sup></td></tr><tr><td>ν<sub>NADH<sub>2</sub>-O<sub>2</sub></sub></td><td>2</td><td>-</td></tr><tr><td>K</td><td>40,000</td><td>µmol<sub>hν</sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>K<sub>r</sub></td><td>556.5</td><td>µmol<sub>hν</sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr><tr><td>A<sub>c</sub></td><td>2800</td><td>µmol<sub>hν</sub> kg<sup>-1</sup> s<sup>-1</sup></td></tr></tbody></table>
<figure><img src="image_17.png" /></figure>

**Figure A1.** Average spectral mass (a) absorption $\bar{A}_{abs,\lambda}$ and (b) scattering $\bar{S}_{sca,\lambda}$ cross-sections taken from [9] for *Chlorella vulgaris* with a mean cell radius of 2 µm and chlorophyll a, chlorophyll b, and photoprotective carotenoid (PPC) concentrations of 3.3 wt.%, 0.85 wt.%, and 0.85 wt.%, respectively. (c) Spectral backward scattering ratio $b_{\lambda}$ as a function of wavelength $\lambda$ calculated using the method presented in [35].
<figure><img src="image_18.png" /><figcaption>**Figure A2.** Comparison of the average specific growth rate $\bar{\mu}(t)$ on (a) 21 June and (b) 23 September and the biomass concentration $X(t)$ on (c) 21 June and (d) 23 September for droplet-free PBRs simulated using a time increment $\Delta t$ of 0.5 h and 2 h. The PBRs had an initial biomass concentration $X_0$ of 0.07 $gL^{-1}$ and a culture depth $L$ of 0.3 m. The maximum relative error in the biomass concentration $X(t)$ was 0.1% and 0.3% on 21 June and 23 September, respectively.</figcaption></figure>

References

1. Chisti, Y. Biodiesel from microalgae. *Biotechnol. Adv.* **2007**, *25*, 294–306. [CrossRef] [PubMed]

2. Chew, K.W.; Yap, J.Y.; Show, P.L.; Suan, N.H.; Juan, J.C.; Ling, T.C.; Lee, D.J.; Chang, J.S. Microalgae biorefinery: High value products perspectives. *Bioresour. Technol.* **2017**, *229*, 53–62. [CrossRef] [PubMed]

3. Tredici, M.R. Mass Production of Microalgae: Photobioreactors. In *Handbook of Microalgal Culture: Biotechnology and Applied Phyiology*; Richmond, A., Ed.; Blackwell Publishing Ltd.: Cornwall, UK, 2004; Chapter 9; pp. 178–214.

4. Pruvost, J.; Cornet, J.-F.; Pilon, L. Large-scale production of algal biomass: Photobioreactors. In *Algae Biotechnology: Products and Processes*; Bux, F., Chisti, Y., Eds.; Springer International Publishing: Cham, Switzerland, 2016; pp. 41–66.

5. Pruvost, J.; Borgne, F.; Artu, A.; Legrand, J. Development of a thin-film solar photobioreactor with high biomass volumetric productivity (AlgoFilm©) based on process intensification principles. *Algal Res.* **2017**, *21*, 120–137. [CrossRef]

6. Grobbelaar, J.U. Factors governing algal growth in photobioreactors: The “open” versus “closed” debate. *J. Appl. Phycol.* **2009**, *21*, 489–492. [CrossRef]

7. Huang, Y.; Feng, C.; Hoeniges, J.; Zhu, K.; Pilon, L. Bidirectional transmittance of transparent windows with external or backside condensation of nonabsorbing cap-shaped droplets. *J. Quant. Spectrosc. Radiat. Transf.* **2020**, *251*, 107039. [CrossRef]

8. Zhu, K.; Pilon, L. Transmittance of semitransparent windows with non-absorbing cap-shaped droplets condensed on their backside. *J. Quant. Spectrosc. Radiat. Transf.* **2017**, *194*, 98–107. [CrossRef]
9. Souliès, A.; Legrand, J.; Marec, H.; Pruvost, J.; Castelain, C.; Burghelea, T.; Cornet, J.F. Investigation and modeling of the effects of light spectrum and incident angle on the growth of *Chlorella Vulgaris* in photobioreactors. *Biotechnol. Prog.* **2016**, *32*, 247–261. [CrossRef]

10. Hsieh, C.K.; Rajvanshi, A.K. The effect of dropwise condensation on glass solar properties. *Sol. Energy* **1997**, *19*, 389–393. [CrossRef]

11. Briscoe, B.J.; Galvin, K.P. The effect of surface fog on the transmittance of light. *Sol. Energy* **1991**, *46*, 191–197. [CrossRef]

12. Zhu, K.; Li, S.; Pilon, L. Light transfer through windows with external condensation. *J. Quant. Spectrosc. Radiat. Transf.* **2018**, *208*, 164–171. [CrossRef]

13. Pollet, I.V.; Pieters, J.G. Condensation and radiation transmittance of greenhouse cladding materials, Part 3: Results for glass plates and plastic films. *J. Agric. Eng. Res.* **2000**, *77*, 419–428. [CrossRef]

14. Pollet, I.V.; Pieters, J.G. Condensation and radiation transmittance of greenhouse cladding materials, Part 2: Results for a complete condensation cycle. *J. Agric. Eng. Res.* **2000**, *75*, 65–72. [CrossRef]

15. Pollet, I.V.; Pieters, J.G. PAR transmittances of dry and condensate covered glass and plastic greenhouse cladding. *Agric. For. Meteorol.* **2002**, *110*, 285–298. [CrossRef]

16. Bhardwaj, R.; Kortenaar, M.V.t.; Mudde, R.F. Influence of condensation surface on solar distillation. *Desalination* **2013**, *326*, 37–45. [CrossRef]

17. Simsek, E.; Zhu, K.; Kashanchi, G.N.; Williams, M.J.; Galy, T.; Marszewski, M.; Tolbert, S.H.; Pilon, L. Light transfer through semi-transparent glass panes supporting pendant droplets. *J. Quant. Spectrosc. Radiat. Transf.* **2020**, under review. [CrossRef]

18. Zhu, K.; Pilon, L. Transmittance of semitransparent windows with absorbing cap-shaped droplets condensed on their backside. *J. Quant. Spectrosc. Radiat. Transf.* **2017**, *201*, 53–63. [CrossRef]

19. Rubin, M. Optical properties of soda lime silica glasses. *Sol. Energy Mater.* **1985**, *12*, 275–288. [CrossRef]

20. Hale, G.M.; Query, M.R. Optical constants of water in the 200-nm to 200-µm wavelength region. *Appl. Opt.* **1973**, *12*, 555–563. [CrossRef]

21. Pottier, L.; Pruvost, J.; Deremetz, J.; Cornet, J.-F.; Legrand, J.; Dussap, C.-G. A fully predictive model for one-dimensional light attenuation by *Chlamydomonas Reinhardtii* A Torus Photobioreactor. *Biotechnol. Bioeng.* **2005**, *91*, 569–582. [CrossRef]

22. Lee, E.; Pruvost, J.; He, X.; Munipalli, R.; Pilon, L. Design tool and guidelines for outdoor photobioreactors. *Chem. Eng. Sci.* **2014**, *106*, 18–29. [CrossRef]

23. Cornet, J.-F.; Dussap, C.G.; Dubertret, G. A structured model for simulation of cultures of the cyanobacterium *Spirulina Platensis* Photobioreactors: I. Coupling between Light Transf. Growth Kinet. *Biotechnol. Bioeng.* **1992**, *40*, 817–825. [CrossRef] [PubMed]

24. Cornet, J.-F.; Dussap, C.G.; Gros, J.B.; Binois, C.; Lasseur, C. A simplified monodimensional approach for modeling coupling between radiant light transfer and growth kinetics in photobioreactors. *Chem. Eng. Sci.* **1995**, *50*, 1489–1500. [CrossRef]

25. Pruvost, J.; Cornet, J.-F.; Goetz, V.; Legrand, J. Theoretical investigation of biomass productivities achievable in solar rectangular photobioreactors for the cyanobacterium *Arthrospira platensis*. *Biotechnol. Prog.* **2012**, *28*, 699–714. [CrossRef] [PubMed]

26. Dunn, I.; Heinzle, E.; Ingham, J.; Prenosil, J. *Biological Reaction Engineering: Dynamic Modelling Fundamentals with Simulation Examples*, 2nd ed.; Wiley-VCH: Weinheim, Germany, 2003.

27. Takache, H.; Pruvost, J.; Cornet, J.-F. Kinetic modeling of the photosynthetic growth of *Chlamydomonas Reinhardtii* A Photobioreactor. *Biotechnol. Prog.* **2012**, *28*, 681–692. [CrossRef] [PubMed]

28. Versyck, K.J.; Claes, J.E.; van Impe, J.-F. Practical identification of unstructured growth kinetics by application of optimal experimental design. *Biotechnol. Prog.* **1997**, *13*, 524–531. [CrossRef]

29. Fouchard, S.; Pruvost, J.; Degrenne, B.; Titica, M.; Legrand, J. Kinetic modeling of light limitation and sulfur deprivation effects in the induction of hydrogen production with *Chlamydomonas reinhardtii*: Part I. Model Dev. Parameter Identification. *Biotechnol. Bioeng.* **2009**, *102*, 232–245. [CrossRef]

30. Cornet, J.-F.; Dussap, C.-G. A simple and reliable formula for assessment of maximum volumetric productivities in photobioreactors. *Biotechnol. Prog.* **2009**, *25*, 424–435. [CrossRef]

31. Pruvost, J.; Vooren, G.; Gouic, B.; Couzinet-Mossion, A.; Legrand, J. Systematic investigation of biomass and lipid productivity by microalgae in photobioreactors for biodiesel application. *Bioresour. Technol.* **2011**, *2*, 150–158. [CrossRef]

32. Slegers, P.M.; Wijffels, R.H.; van Straten, G.; van Boxtel, A.J.B. Design scenarios for flat panel photobioreactors. *Appl. Energy* **2011**, *88*, 3342–3353. [CrossRef]

33. Modest, M.F. *Radiative Heat Transfer*, 3rd ed.; Elsevier Inc.: Oxford, UK, 2013.

34. Howell, J.R.; Mengüç, M.P.; Siegel, R. *Thermal Radiation Heat Transfer*, 6th ed.; CRC Press: Taylor and Francis: Boca Raton, FL, USA, 2016.

35. Kandilian, R.; Pruvost, J.; Artu, A.; Lemasson, C.; Legrand, J.; Pilon, L. Comparison of experimentally and theoretically determined radiation characteristics of photosynthetic microorganisms. *J. Quant. Spectrosc. Radiat. Transf.* **2016**, *175*, 30–45. [CrossRef]

36. Gueymard, C. Simple Model of the Atmospheric Radiative Transfer of Sunshine (SMARTS). Version 2.9.5. 2005. Available online: https://www.nrel.gov/grid/solar-resource/smarts.html (accessed on 30 December 2020).

37. Figgis, B.; Nouviaire, A.; Wubulikasimu, Y.; Javed, W.; Guo, B.; Ait-Mokhtar, A.; Belarbi, R.; Ahzi, S.; Rémond, Y.; Ennaoui, A. Investigation of factors affecting condensation on soiled PV modules. *Sol. Energy* **2018**, *159*, 488–500. [CrossRef]