Next Article in Journal
Impact Identification of Carbon-Containing Carboniferous Clays on Surfaces of Friction Nodes
Next Article in Special Issue
Fischer–Tropsch Synthesis as the Key for Decentralized Sustainable Kerosene Production
Previous Article in Journal
Transient Thermal Analysis of NH000 gG 100A Fuse Link Employing Finite Element Method
Previous Article in Special Issue
Solar Prosumers in the German Energy Transition: A Multi-Level Perspective Analysis of the German ‘Mieterstrom’ Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Objective Life Cycle Optimization Model of an Integrated Algal Biorefinery toward a Sustainable Circular Bioeconomy Considering Resource Recirculation

by
Celine Marie A. Solis
1,*,
Jayne Lois G. San Juan
1,2,
Andres Philip Mayol
3,
Charlle L. Sy
1,2,
Aristotle T. Ubando
2,3 and
Alvin B. Culaba
2,3
1
Industrial Engineering Department, De La Salle University, Manila 0922, Philippines
2
Center for Engineering and Sustainable Development Research, De La Salle University, Manila 0922, Philippines
3
Mechanical Engineering Department, De La Salle University, Manila 0922, Philippines
*
Author to whom correspondence should be addressed.
Energies 2021, 14(5), 1416; https://doi.org/10.3390/en14051416
Submission received: 26 January 2021 / Revised: 25 February 2021 / Accepted: 26 February 2021 / Published: 4 March 2021

Abstract

:
Biofuel production from microalgae biomass has been considered a viable alternative to harmful fossil fuels; however, challenges are faced regarding its economic sustainability. Process integration to yield various high-value bioproducts is implemented to raise profitability and sustainability. By incorporating a circular economy outlook, recirculation of resource flows is maximized to yield economic and environmental benefits through waste minimization. However, previous modeling studies have not looked into the opportunity of integrating productivity reduction related to the continuous recirculation and reuse of resources until it reaches its end of life. In this work, a novel multi-objective optimization model is developed centered on an algal biorefinery that simultaneously optimizes cost and environmental impact, adopts the principle of resource recovery and recirculation, and incorporates the life cycle assessment methodology to properly account for the environmental impacts of the system. An algal biorefinery involving end-products such as biodiesel, glycerol, biochar, and fertilizer was used for a case study to validate the optimization model. The generated optimal results are assessed and further analyzed through scenario analysis. It was seen that demand fluctuations and process unit efficiencies have significant effect on the optimal results.

Graphical Abstract

1. Introduction

The overwhelming environmental concerns associated with fossil fuel production and consumption are only expected to rise in the following years as global energy demand grows by around 1.3% every year [1]. This proves that further research into alternative renewable sources of energy is necessary [2]. Biodiesel, which is considered as the main substitute for fossil fuel, has numerous biomass feedstock alternatives, with microalgae being considered the most valuable due to its high growth productivity and photosynthetic efficiency [3,4,5]. Microalgae are aquatic microorganisms that grow in the process of photosynthesis, converting nutrients (carbon dioxide, water, etc.) into energy [6,7]. In addition, growing microalgae uses a lower area of arable land compared to other biofuel candidates such as corn, soybean, canola, and jatropha [8]. One of the promising attributes that makes microalgae a viable option for alternative biofuel is the ability to harness carbon dioxide from the atmosphere [9].
However, economic limitations are faced when looking into microalgal biofuel production due to its higher overall costs as compared to the production of fossil fuels. Proper investment strategies were found to be an integral factor to ensure long-term sustainability and growth in the diffusion of algal biofuel [10]. Aside from that, the processes involved in biofuel production typically generate a lot of waste including unused biomass waste. To maximize the algal biomass used in biofuel production, integrated process methods are incorporated into the production system through biorefineries. An algal biorefinery is built to produce energy and along with other products derived from biomass. Examples of these are biogas, biodiesel, bioethanol, fertilizers, and cosmetics. The incorporation of additional high-value products leads to increased profitability and biomass utilization for the biorefinery.
With that, numerous researchers have investigated the integration of processes with the aim of a zero-waste biofuel production system. Mitra and Mishra (2019) and Mohan et al., (2019) investigated self-sustaining algal biorefinery design models with a closed-loop approach regarding the byproduct utilization in the system [11,12]. In particular, Mitra and Mishra (2019) discussed the current state of the art and potential advancements of a zero-waste biorefinery from Arthrospira spp. through the valorization of multiple products, such as biofertilizers, jet fuels, nutraceuticals, and c-phycocyanin, while Mohan et al., (2019) looked into the co-production of monomeric sugars, bioalcohols, and biohydrogen [11,12]. De Bhowmick et al., (2019) proposed biofuel and biochar co-production to minimize and possibly achieve net zero waste discharge [13]. Hemalatha et al., (2019) studied the potential of a microalgae biorefinery with a cascading resource recovery design through the integration of dairy wastewater treatment processes [14]. Wu and Chang (2019) evaluated an integrated microalgal biorefinery from a process systems engineering standpoint using life cycle assessment and techno-economic assessment methodologies to appraise the biorefinery’s environmental and economic benefits [15]. However, the previously mentioned studies simply review the zero-waste concept applied to an algal biorefinery without mathematically supported decision making presented. The inherent complexity of designing biorefineries with several significantly interdependent processes requires appropriate systematic management tools, such as optimization modeling, to ensure that the potential benefits of closed-loop biorefineries are achieved, such as reduced resource consumption and waste or pollution generation, while minimizing any unintended negative consequences [16,17]. Looking specifically into optimization models on algal biorefineries, Garcia Prieto et al., (2017) developed a mixed integer nonlinear programming model for the design of an integrated algal biorefinery aiming to maximize its net present value [18]. However, the study only focused on a single objective involving economic impact, which indicates that their solution may not be environmentally sustainable. This has been addressed through a superstructure optimization model proposed by Gong and You (2014) covering algae cultivation, harvesting, drying, and lipid extraction to minimize costs and greenhouse gas emissions, achieving zero emissions throughout the process through reutilization [19]. The same reutilization of all wastes and products may be achieved with proper design. A multi-objective target-oriented robust optimization model was developed by Sy et al., (2018) for an integrated algal biorefinery polygeneration system focusing on profit maximization and environmental footprint minimization, which was later extended by Culaba et al., (2019) into a multi-period model to capture fluctuations in demand through planning periods and allow for material storage to buffer disparities between available supply and spikes in demand [16,17]. However, these studies did not take into account the specific inputs and outputs of each process in the integrated algal biorefinery by treating the algal biorefinery as a black box process unit. Each process within an algal biorefinery entails several potentially complicated and conflicting decisions that should also be taken into account. To address these, Ching et al., (2021) make use of artificial intelligence-based methodologies to model and possibly optimize the vacuum drying process of Chlorococcum infusionum for algal biofuel production [20]. Focusing instead on the cultivation and harvesting processes, San Juan et al., (2020) proposed a scheduling and planning optimization algorithm that maximizes profit and minimizes environmental impact through decisions such as species selection and cultivation and harvesting technology selection and scheduling [21]. Nonetheless, these processes within an algal biorefinery cannot only be modeled as stand-alone systems, under the assumption that no interaction exists between them. Accounting for interactions between these processes and other process systems that can operate symbiotically with algal biorefineries may be helpful in improving economic and environmental performances. A superstructure optimization model was formulated by Galanopoulos et al., (2019) integrating an algal biorefinery with a straw wheat biorefinery to minimize a single objective (i.e., biofuel production costs) [22]. Caligan et al., (2020) modeled the integration of a wastewater treatment plant and an algal biorefinery connected through the potential of cultivating algae on wastewater [23]. Extending this, San Juan et al., (2020) proposed a multi-objective optimization model to design an integrated system consisting of a wastewater treatment park, an algal biorefinery, and a sludge-based bioenergy conversion system [24]. Čuček et al., (2014) formulated a mixed integer linear programming model for the synthesis of a bioenergy supply network utilizing a combination of first-, second-, and third-generation biofuels (e.g., bioethanol, biohydrogen, and biodiesel) to maximize economic performance [25]. However, it is important to note that these systems do not feature a self-sustaining algal biorefinery because of their dependence on other supplementary systems. Integrating these facilities may not necessarily be practical when the facilities have already been established as stand-alone systems geographically far away from each other, as installing connections between them may be costly, inefficient, or impossible.
To be able to accurately calculate the environmental impact generated by a certain product, the life cycle assessment (LCA) methodology is used. According to the ISO 14040 standard, a full product life cycle is measured starting from the extraction of raw materials until the product end-of life. The LCA methodology has been applied to existing studies relating to algal biorefineries. Chowdhury et al., (2018) made use of the LCA methodology to simulate the effect of residence time on the production of algal bioenergy from dairy manure [26]. The LCA model developed was able to predict that with diluted waste and shorter residence time, energy requirement and greenhouse gas (GHG) emissions for the bioenergy production can be reduced. Moreover, Barlow et al., (2016) used the LCA methodology along with techno-economic assessment (TEA) to identify the sustainability of incorporating hydrothermal liquefaction and wastewater treatment into the algal biofuel production system [27]. It was determined that incorporating wastewater treatment into the system was able to reduce the environmental impact of the biorefinery. However, the LCA methodology simply provides an evaluation for a pre-determined system and cannot alone be used in system design.
The life cycle optimization (LCO) methodology was created to be able to integrate the advantages of the LCA methodology into an optimization model for application in system design. With the use of TEA and LCA methodologies, Wu et al., (2018) were able to develop an optimization model to obtain the optimal combination of cultivation and pretreatment process chains for an algal biorefinery [28]. The optimization model aimed to minimize environmental impact and maximize revenue. While having resource recirculation in the form of heat and CO2 recovery, the study did not account for waste recovery and reuse for wastewater and biomass residue, which are possible sources of revenue and can lead to worsened environmental impact in the long run.
The LCA methodology is incorporated into this optimization model through the environmental optimization objective. In the previously mentioned study by Wu et al., (2018), the impact assessment method used was life cycle GHG emissions [28]. Bussa et al., (2018) stated that the ReCiPe impact assessment method is most recommended for microalgal biorefinery systems as it covers numerous environmental factors, not just concerning energy and greenhouse gas emissions [29]. Therefore, it was decided that the ReCiPe impact assessment method would be used for the life cycle optimization model.
Therefore, there is a need to develop a mathematical optimization model on an integrated microalgal biorefinery that optimizes profit and environmental impact, integrates the life cycle assessment methodology to account for the environmental impacts associated with a microalgal biofuel process flow, and incorporates resource recirculation among materials recirculated within the biorefinery centered on a closed-loop production system.

2. System Definition

The algal biorefinery featured in the study is divided into two main sections. The process units included in the first section are microalgae cultivation, harvest and dewatering, and lipid extraction, which are the upstream processes [30]. In this study, different process unit alternatives are selected for the upstream processes, and a comparison of each is presented in Table 1.
The second section includes the downstream processes of the biorefinery, the conversion of algal biomass into its intended end-products such as biodiesel, glycerol, biochar, and fertilizer. The main output of the oil extraction process is algal lipids turned into biodiesel and glycerol through transesterification [34]. Aside from that, waste is generated from the lipid extraction process in the form of solid and liquid algal biomass residues. These residues are converted into bioproducts such as biochar and fertilizer through pyrolysis and anaerobic digestion, respectively [16]. With the methane generated from the anaerobic digestion process, a link to a combined heat and power plant is established for conversion to power and heat with the aim of recirculation back to the biorefinery. The biorefinery process flow is presented in Figure 1. Aside from the main processes involving microalgal biofuel production, recovery process units are also incorporated into the system. These are added to incorporate resource recirculation into the biorefinery through the zero-waste concept. Process units such as microalgae cultivation, harvesting, and lipid extraction each have recovery processing alternatives that are to be decided by the optimization model. Moreover, model decisions are also made regarding the usage of resource recovery process units.
For the optimization model, the following cost components will be considered namely: investment costs, operating costs, material purchase costs, utility costs, and inventory costs. Investment costs refer to the initial capital expenses incurred to build the algal biorefinery process units expressed on a cost per process unit basis. Operating costs include labor, overhead, and material costs associated with each process unit expressed as a cost per time period. Material purchase costs are set for each material input that is used up in the production system, expressed as a cost per kg unit. Given that the optimization model considers multiple time periods, inventory costs are included for each material expressed as a periodic cost per kg material. Since all process units involved in the microalgal biorefinery are all assumed to be in a single facility, transportation costs are not considered in the optimization model. These cost components are deducted from the revenue streams which would be derived from the sale of products. Demand and selling price parameters are determined for the final end-products of the biorefinery and presented in Table A1, while cost parameters are presented in Table A2 and Table A3 of the Appendix A.
The material inputs that are handled in resource recovery process units include water, solvent, flocculant, and catalyst. The optimization model involves multiple time periods to be able to incorporate the cyclability of resource recirculation. Furthermore, capital and operating costs are properly accounted for in the algal biorefinery with the consideration of multiple periods. For this study, a base period of 10 years will be used.

3. Life Cycle Assessment

For the environmental objective of the model, measures for environmental impact were obtained with the life cycle assessment (LCA) methodology. The goal of the study was to evaluate the total environmental impacts of each process path in the algal biorefinery. The system boundary of the study is presented in Figure 1. The study looked into each biorefinery product from its raw material extraction until its usage and end-of-life disposal scenarios. For LCA normalization, the functional unit to be utilized in the model was 1 kg biodiesel produced in the biorefinery. Various research articles such as those written by Gnansounou and Raman (2017), Dasan et al., (2019), Biller et al., (2013), Sy et al., (2018), Garcia Prieto et al., (2017), and Wu et al., (2018) were used as the basis for the material inputs and outputs as well as the energy consumption of each process unit [16,18,28,35,36,37]. After assessing the life cycle inventory, the researchers determined the overall impact of the system using the SimaPro software. For the impact assessment portion of the study, the ReCiPe 2016 Endpoint model was used in the simulation to be able to yield a single score impact value for each configuration in the microalgal production system. The environmental impact values generated from the LCA methodology were used as inputs for the objective function of the model to represent the environmental impact of each process path. Once a single score impact was determined for each production process configuration in the algal biorefinery, the one with the lowest single score was considered as the most sustainable configuration. However, this does not automatically mean that a certain configuration is optimal since a balance must be met between the economic and environmental objectives.
Data for the life cycle inventory for the different process paths are presented in Table A4 and Table A5 of the Appendix A. Each process path was evaluated with a single score that represents its overall environmental impact. Impact categories were subject to normalization and weighting factors (see Table A6 of Appendix A) based on the ReCiPe 2016 method, aggregated to obtain the single score environmental impact. Looking firstly into upstream processes only, namely cultivation, harvesting, and oven drying, the resulting environmental impacts are presented in Figure 2. It is evident that the path with the greatest impact includes the photobioreactor cultivation and centrifugation for harvesting. The greatest contributor to this impact is the land use category due to the use of compost for cultivation and the larger input requirement for centrifugation. Moreover, it is seen that the upstream process with the highest contribution to environmental impact is centrifugation, attributable to its large biomass requirement for biodiesel production. On the other hand, the configuration with the least environmental impact is open pond cultivation paired with flocculation.
Consequently, the environmental impacts of the downstream process paths are displayed in Figure 3. The path with microwave extraction is shown to have the biggest environmental impact since it is more energy-intensive. The categories with the most contributions to the environmental impact of both scenarios are human non-carcinogenic toxicity and human health global warming attributable to the large energy requirements of anaerobic digestion.
Finally, looking into all process paths, the impacts are all displayed in Figure 4. The breakdown of environmental impacts for each process configuration and impact category is also numerically presented in Table A7 of the Appendix A. The paths with the largest environmental impact are those involving photobioreactor cultivation and centrifugation, since it is more energy-intensive. Similar to the earlier result, land use is the greatest contributor to the environmental impact values because of the use of compost for the nutrients used to cultivate the algal biomass. Among all the process paths, the one with the least impact is open pond cultivation with flocculation harvesting and solvent extraction, which makes it the most likely process path for the algal biorefinery.

4. Model Formulation

A mixed integer linear programming (MILP) model was developed for the study with the objectives of optimizing the profit and environmental impacts while satisfying demand and capacity limitations. The indices, parameters, and variables are displayed in Table A8, Table A9 and Table A10 of the Appendix A.

4.1. Model Assumptions

  • All parameters considered in this model are deterministic and known with certainty.
  • The outputs produced by the facility are transported to customers at the same time.
  • Processing of algal biomass in all facilities is instantaneous.

4.2. Optimization Model

4.2.1. Objective Functions

The profit is calculated by working out the difference between the revenue generated from all products and the overall costs throughout the biorefinery, presented in Equation (1). Equation (2) defines the revenue for each period, which is computed by multiplying the selling price of each product by the total final outputs.
Profit = t Revenue t Total   Cost
Revenue t = p t SP pt TO pt
The breakdown of the total costs is presented in Equation (3), which are investment costs, operating costs, material purchase costs, and inventory costs.
Total   Cost = Investment + t Operating t + t Purchase t + t Inventory t
The total investment is the sum product of the fixed costs per process unit and the binary variable for selecting the alternative, presented in Equation (4). The operating costs for each period in the biorefinery are defined as the sum product of the operating costs per process and the corresponding product output, as shown in Equation (5).
Investment = i FC i BC i + j FH j BH j + k FE k BE k + r FR r BR r + u FP u BP u
Operating t = i OC it PC it + j OH jt PH jt + k OE kt PE kt + r OR rt PR rt + u OP ut PP ut
Equations (6) and (7) indicate the calculations for purchase and inventory costs, respectively. The material purchase costs are calculated by multiplying the purchase cost per input material with the total amount purchased each period. Meanwhile, the inventory costs are the unit storage costs for each input and output material and the total inventory level in each process unit.
Purchase t = m MC mt MQ mt
Inventory t = m IC mt EI mt
The environmental impact minimization objective is defined in Equation (8). It is the sum product of the environmental impact per output of each process unit and their corresponding total output.
Impact = i EC i PC i + j EH j PH j + k EE k PE k + r ER r PR r + u EP u PP u
Since the model has dual objectives, there must be a balance between the economic and environmental objectives to generate the optimal solution. The objective function is defined as the maximization of the least desired value to balance the two objectives as seen in Equation (9) [38]. Equations (10) and (11) define the efficiencies for each objective, obtained by calculating the ratio of the attained improvement, which is the actual value subtracted from the worst possible one, and the potential improvement, which is the best possible value subtracted from the worst possible one. The best possible values for the two objectives are acquired through the optimization of each corresponding objective using a single objective linear programming model. The assumption is that the worst possible value for the environmental impact objective is its value when the profit objective is optimized, and vice versa.
Max   Z = min Eff Profit , Eff Impact
Eff Profit = Profit worst Profit Profit worst Profit best
Eff Impact = Impact worst Impact Impact worst Impact best
With the nonlinear nature of the objective function defined above, it is necessary to include linearizing constraints to the model to make sure that the optimization model generates the optimal solution. Equations (13) and (14) illustrate that the final value for efficiency is equal to the minimum of the efficiencies of the two objectives [38], while Equation (12) shows the final objective function.
Max Z = Efficiency
Efficiency ≤ EffProfit
Efficiency ≤ EffImpact

4.2.2. Constraints

The model constraints regarding the demand for the various products of the biorefinery are presented in Equation (15). The total amount of each product to be transported must be greater than or equal to the customer demand for each product, namely biodiesel, glycerol, biochar, and fertilizer.
TO pt D pt p   t
The capacity constraints define the production capability of each process as defined in Equations (16)–(20). The overall production output of each process is set to be less than or equal to the production capacity of that process unit multiplied by the binary variable assigned to each process.
PC it BC i * CC i i   t
PH jt BH j * CH j j   t
PE kt BE k * CE k k   t
PR rt BR r * CR r r   t
PP ut BP u * CP u u   t
The product input-to-output conversion is defined in process constraints shown in Equations (21)–(26). For the conversion of product inputs to their corresponding outputs, the output of the previous process is set to be greater than or equal to the overall production output of the process multiplied by the conversion output yield.
j PH jt YH j i PC it t
PP Dry , t * YP Dry j PH jt t
k PE kt YE k PP Drying , t t
PP Trans , t * YP Trans 0.6 * k PE kt t
PP AD , t * YP AD 0.3 * k PE kt t
PP Pyro , t * YP Pyro 0.1 * k PE kt t
PP CHP , t * YP CHP 0.6 * PP AD , t t
The equations regarding the relationships between the total final product outputs and the production outputs of the conversion processes are displayed in Equations (28)–(31). The final product outputs are limited by the conversion efficiencies of their respective processes.
TO Biodiesel , t 0.9 * PP Trans , t t
TO Glycerol , t 0.1 * PP Trans , t t
TO Biochar , t PP Pyro , t t
TO Fertilizer , t 0.4 * PP AD , t t
As for the equations presenting how the production outputs of select processes branch out into inputs and outputs of corresponding next processes, they are displayed in Equations (32)–(34). Equation (32) details the flow of materials from the extraction processes into the anaerobic digestion, pyrolysis, and transesterification process units. Similarly, Equation (33) describes biodiesel and glycerol to be the product outputs of the transesterification process, while the outputs of the anaerobic digestion unit include the final product fertilizer and inputs into combined heat and power process unit.
k PE kt PP Trans , t * YP Trans + PP AD , t * YP AD + PP Pyro , t * YP Pyro t
PP Trans , t TO Biodiesel , t + TO Glycerol , t t
PP AD , t TO Fertilizer , t + PP CHP , t * YP CHP t
The relationship between the recovered material inputs in each facility and the respective material input usage in production for each facility is presented in Equations (35)–(38). Equations (35) and (36) model the recovery of water and flocculant from harvesting, respectively, while Equation (37) represents the recovery of solvent from the extraction process alternatives. The same is done for the recovery of catalyst from the transesterification process in Equation (38).
j PH jt YH j 0.9 * PR Water , t t
j PH jt YH j 0.05 * PR Floc , t t
k PE kt YE k 0.1 * PR Solv , t t
PP Trans , t * YP Trans 0.2 * PR Cat , t t
Equations (39)–(42) describe the flow of the material inventory in the system understudy. The beginning inventory is defined in Equation (39) as a function of material purchases and the ending inventory carried over from the previous period, while the ending inventory is expressed as a function of the beginning inventory and recovered material inputs in Equation (40). Equations (41) and (42) ensure that the amount of material purchased and kept inventory are limited by the available inventory and purchase capacities.
BI mt = MQ mt + EI mt 1 t
EI mt = BI mt + PR mt t
EI mt IK mt t
MQ mt MK mt t
The available material inputs (i.e., water, flocculant, solvent, and catalyst) at the beginning of each period should be able to satisfy the requirements of their respective conversion processes (i.e., cultivation, harvesting, extraction, and transesterification). These constraints are modeled in Equations (43)–(46).
BI Water , t 0.9 * i PC it YC i t
BI Floc , t 0.05 * j PH jt YH j t
BI Solv , t 0.1 * k PE kt YE k t
BI Cat , t 0.2 * PP Trans , t * YP Trans t
Equations (47)–(49) require that the system activates one of the available alternatives in each of the cultivation, harvesting, and extraction processes, while Equation (50) ensures that other processes without alternative units must be installed through binary constraints.
i BC i = 1
j BH j = 1
k BE k = 1
BP u = 1 u

5. Model Validation

The mathematical optimization model was validated with the use of MATLAB R2020a software along with the Cplex optimization solver. A base case is considered as an initial scenario to be used for comparison in the scenario analysis. In order to attain the best and worst possible values involving the two objectives which were used in calculating the efficiency, optimization with single objectives was executed.

5.1. Profit Maximization

With the objective of profit maximization, the algal biorefinery chose to operate at its maximum yielding annual biodiesel output of 49,473.14 kg. The optimal selection of process alternatives with only profit consideration for cultivation, harvesting, and extraction is shown in Figure 1. Based on the results, the process units chosen were those that require the least fixed cost among other alternatives, which were open pond cultivation, centrifugation harvesting, and solvent extraction. However, the maximized production of biodiesel results in a trade-off for the environmental impact objective since this process configuration is one of those with the biggest impact as calculated from the LCA methodology.

5.2. Impact Minimization

Consequently, looking into the environmental objective, biodiesel production was performed only to meet the minimum annual demand of 25,000 kg. A larger production amount would result in a higher impact given that the single score values are calculated based on a unit output for each process. Figure 1 presents the process path that yields the least environmental impact as calculated from the life cycle assessment. However, the resulting profit value with this process flow drastically changes, resulting in a net loss (see Table 2) given that the revenue generated is lower than the total production cost since less biofuel and other bioproducts are produced from the microalgal biorefinery.
The results presented above validate the hypothesis that a balance must be established between the two objectives to obtain the optimal result. With that, the developed model with the efficiency objectives highlighted in the model formulation section were added onto the final model.

5.3. Multi-Objective Model

Using the obtained results from the single-objective execution of the optimization model, the best and worst values for the profit and environmental impact were retrieved. With that, the efficiency could be calculated for both objectives which was then maximized for the optimization model.
The resulting optimal efficiency value obtained from the model was 92.69%, as presented in Table 3, which means that the optimal multi-objective values are very close to their respective single objective values. Annual biodiesel production equates to 49,951.43 kg, which is higher than that of the first two single-objective runs while having a lower profit than the maximum profit single objective. This shows the significance of incorporating certain process paths in the system. Looking simply into the economic objective, a bigger profit is achieved with less output due to the use of low-cost processes. Figure 1 shows the optimal process path chosen for the biorefinery, which is identical to that of the minimum impact single objective. This shows the balance that was established between the two objectives to yield the optimal result.
Looking into the behavior of product purchasing given the incorporation of the recovery processes, purchase amounts for each product are presented in Figure 5.
For the first period, the maximum purchase amounts are seen since there is no initial inventory count and no material recovered yet. By the second year, due to the recovered materials from the previous period, a decrease in the purchase amount is observed. However, since the recovery efficiency decreases over time, purchase amounts start to increase after the second year. This behavior is beneficial to the algal biorefinery due to the purchase cost savings and an overall decrease in use of new materials. To discuss further, the inventory level of each material over time is assessed, as presented in Figure 6.
It can be observed that the behaviors of the purchase amounts and inventory levels of each material follow an inverse relationship. The inventory rose to maximum levels during the second period due to the newly recovered materials from the first period. Moreover, the inventory started to decrease, attributable to the decreasing efficiency level of the recovery processes. Eventually, the inventory of materials would reach a minimum, which would then indicate that a full purchase of new materials is needed.

6. Scenario Analysis

To assess the capabilities of the optimization model further, a scenario analysis was conducted looking into the impacts of applying efficiencies to each process unit and fluctuations in demand for biodiesel.

6.1. Impact of Process Unit Efficiencies

For the base case, it was assumed that all process units are always 100% efficient. However, this is not true most of the time due to the occurrence of waste and differences in raw material input quality. To look into the effect of the application of efficiencies to each process unit, different efficiency levels were assigned to all process units to test the corresponding impact on the objective functions. An efficiency level lower than 100% indicates that a unit of input yields less output compared to its full capacity.
Decreasing the efficiency levels of each process unit yields less environmental impact and less profit for the algal biorefinery (see Figure 7a). Although the decrease in impact is favorable, the decrease in profit is not. This behavior occurs since less products are produced within the system, which yields less revenue and overall environmental impact. Aside from that, it is noted that the efficiency objective functions which are related to the worst and best values for each of the objectives have an inverse relationship with the process unit efficiencies, as seen in Figure 7b. This indicates that the resulting optimal values from the multi-objective model move nearer to the optimal results from the single-objective runs for each objective as efficiency decreases. However, this is not a favorable outcome if looking at the big picture since the profit has an opposite trend.

6.2. Impact of Demand Fluctuations

Product demand does not always assume a deterministic value. Most of the time, fluctuations in demand are caused by various industry factors affecting the operations of the production line. To determine the impact when accounting for demand fluctuations, a demand index, taking values from 1.5 to 0.5 at an interval of 0.1, was incorporated to be multiplied to the current demand parameter to represent the changes in demand. The results of the optimization model with the adapted fluctuations in demand are presented in Figure 8a,b.
A higher demand, as represented by the demand indices ranging 1.5 to 1.1, indicates higher production of biodiesel, which leads to a higher profit and environmental impact (see Figure 8a). A similar trend is seen with regards to having a lower demand. This indicates that the objective values are highly sensitive to the demand parameters that are set in the model. The varying demand parameters come with changes in the efficiency objective values as seen in Figure 8b. A higher demand setting yields a higher efficiency value since the denominator for the profit efficiency, which is the potential improvement for profit, takes on a lower value. This is because the worst possible profit is increased due to taking on a bigger demand. With lower demand comes the lower worst possible profit value, which then results in a decrease in the efficiency objective value. This is an indicator that the demand estimate for the base case may have been too low for the system. This shows the importance of appropriate demand forecasting alongside a system that is able to overcome demand fluctuations.

7. Conclusions and Recommendations

This study introduced a multi-objective optimization model for an algal biorefinery incorporating a circular bioeconomy outlook with objectives of maximizing profit and minimizing environmental impact. To determine the single score impact assessment for the process units involved in the system boundary, a life cycle analysis (LCA) methodology was utilized, looking into the life cycle of the products from cultivation until usage and disposal. Incorporating this technique into the optimization model created a new approach in determining the environmental impact associated with each product to be applied in mathematical modeling. The model aimed to design a system for an algal biorefinery to determine which processes to take in to consideration.
By looking into the profit and environmental impact objectives separately through single-objective optimization, the conflict between the two was observed through their inverse relationship. When profit maximization is prioritized, the resulting environmental impact increases, which is unfavorable due to the harm it will cause for the environment. On the other hand, when environmental impact is minimized, the algal biorefinery obtains a net loss since the generated revenue is not enough to cover the overall production expenses. Considering both objectives simultaneously through the optimization model would strike a balance between profit and environmental impact for the algal biorefinery.
Consideration of process unit efficiencies as well as demand fluctuations causes significant changes in the objective function values. A lower process unit efficiency results in a decrease in profit and impact due to less product outputs. Similarly, having a lower demand parameter yields less profit and environmental impact since biodiesel production is decreased. This shows the impact of the changes in demand and process unit efficiencies that are usually caused by numerous other factors outside the system.
A possible extension of this research may be through the incorporation of uncertain parameters into the optimization model such as demand, microalgae nutrient content, and process unit product yield. The LCA section of the study currently lacks further investigation on the breakdown of environmental impacts in each process. A more in-depth analysis on the results of the LCA methodology can also be added to the research, expounding on the impacts of specific material inputs in the algal biorefinery. Moreover, using real-life data for the estimating system parameters would lead to more accurate findings for applied industry projects.

Author Contributions

Conceptualization, C.M.A.S., J.L.G.S.J., and A.P.M.; model formulation, verification, and validation, C.M.A.S.; writing—original draft preparation, C.M.A.S.; writing—review and editing, J.L.G.S.J., A.P.M., C.L.S., A.T.U., and A.B.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Demand and selling price for each product [16].
Table A1. Demand and selling price for each product [16].
Demand (kg)Selling Price (USD/kg)
Biodiesel25,00012.50
Glycerol25000.78
Biochar10000.50
Fertilizer2000.25
Table A2. Purchase and inventory costs for each recoverable material [18].
Table A2. Purchase and inventory costs for each recoverable material [18].
Purchase Cost (USD/kg)Inventory Cost (USD/kg/pd)
Water0.3670.15
Potassium sulfate0.120.15
Hexane0.410.15
Sulfuric acid0.740.15
Table A3. Fixed and operating costs for each process unit [16,18,36].
Table A3. Fixed and operating costs for each process unit [16,18,36].
Fixed Cost (USD)Operating Cost (USD)
Open Pond734,713.760.50
Photobioreactor1,192,058.790.50
Flocculation923,823.860.50
Centrifugation285,617.510.50
Filtration807,400.000.50
Oven Drying632,120.000.50
Solvent Extraction720,000.000.50
Microwave Extraction936,000.000.50
Transesterification1,050,441.710.50
Pyrolysis357,7600.50
Anaerobic Digestion693,6000.50
Combined Heat and Power459,0000.25
Table A4. Life cycle inventory for 1 kg of biodiesel (open pond systems) [15,16,18,35,36,37].
Table A4. Life cycle inventory for 1 kg of biodiesel (open pond systems) [15,16,18,35,36,37].
StagesMaterialOP-FL-SOOP-CN-SOOP-FN-SOOP-FL-MIOP-CN-MIOP-FN-MI
Cultivation—Open PondWater (kg)6.4852221.76026.33851.754659.99931.7150
Urea (g)2.577988.14972.51960.697523.84970.6817
Diammonium phosphate (g)2.318579.27932.26600.627321.44980.6131
Electricity (MJ)0.359912.30770.35180.09743.33000.0952
* Algal broth (kg)6.6203226.38026.47061.791261.24931.7507
Harvesting—FlocculationPotassium sulfate (kg)0.0794 0.0215
Electricity (MJ)7.7385 2.0937
Algal broth (kg)6.6203 1.7912
* Wet biomass (kg)5.4353 1.4706
Harvesting—CentrifugationChitosan (g) 45.9690 12.4374
Electricity (MJ) 12.4740 3.3750
Algal broth (kg) 226.3802 61.2493
* Wet biomass (kg) 5.4353 1.4706
Harvesting—FiltrationElectricity (MJ) 190.2355 51.4700
Algal broth (kg) 6.4706 1.7507
* Wet biomass (kg) 5.4353 1.4706
Oven DryingWet biomass (kg)5.43535.43535.43531.47061.47061.4706
Heat (MJ)62.451662.451662.451616.896916.896916.8969
* Dry biomass (kg)4.62004.62004.62001.25001.25001.2500
Extraction—SolventHexane (g)2.95682.95682.9568
Electricity (MJ)0.33260.33260.3326
Heat (MJ)2.31002.31002.3100
Dry biomass (kg)4.62004.62004.6200
* Algal oil (g)1049.98841049.98841049.9884
* Liquid residue (g)170.0000170.0000170.0000
* Solid residue (g)40.000040.000040.0000
Extraction—MicrowaveElectricity (MJ) 36.749636.749636.7496
Dry biomass (kg) 1.25001.25001.2500
* Algal oil (kg) 1.05001.05001.0500
* Liquid residue (g) 45.995145.995145.9951
* Solid residue (g) 10.822410.822410.8224
TransesterificationAlgal oil (g)1049.98841049.98841049.98841049.98841049.98841049.9884
Methanol (g)124.8787124.8787124.8787124.8787124.8787124.8787
Sodium hydroxide (g)10.487410.487410.487410.487410.487410.4874
Sulfuric acid (g)15.800415.800415.800415.800415.800415.8004
Electricity (MJ)0.16630.16630.16630.16630.16630.1663
Heat (MJ)5.59025.59025.59025.59025.59025.5902
Water (kg)0.13860.13860.13860.13860.13860.1386
* Biodiesel (g)1000.00001000.00001000.00001000.00001000.00001000.0000
* Glycerol (g)113.2825113.2825113.2825113.2825113.2825113.2825
PyrolysisSolid residue (g)40.000040.000040.000010.822410.822410.8224
Heat (MJ)6.41686.41686.41681.73611.73611.7361
* Biochar (g)5.99705.99705.99701.62251.62251.6225
Anaerobic DigestionLiquid residue (g)170.0000170.0000170.000045.995145.995145.9951
Electricity (MJ)54.726054.726054.726014.806614.806614.8066
Heat (MJ)416.2671416.2671416.2671112.6250112.6250112.6250
* Methane (g)19.406419.406419.40645.25065.25065.2506
* Fertilizer (g)69.842869.842869.842818.896618.896618.8966
Combined Heat and PowerMethane (g)19.406419.406419.40645.25065.25065.2506
* Electricity (MJ)0.21560.21560.21560.05830.05830.0583
* Heat (MJ)0.40970.40970.40970.11080.11080.1108
* Output Material; OP—Open Pond; FL—Flocculation; CN—Centrifugation; FN—Filtration; SO—Solvent; MI—Microwave.
Table A5. Life cycle inventory for 1 kg of biodiesel (photobioreactor) [15,16,18,35,36,37].
Table A5. Life cycle inventory for 1 kg of biodiesel (photobioreactor) [15,16,18,35,36,37].
StagesInputPB-FL-SOPB-CN-SOPB-FN-SOPB-FL-MIPB-CN-MIPB-FN-MI
Cultivation—Open PondWater (kg)5.6105191.84765.48361.518051.90621.4836
Chicken compost (kg)0.22447.67390.21930.06072.07620.0593
Atmospheric air (kg)0.22447.67390.21930.06072.07620.0593
Inoculum (kg)0.561019.18480.54840.15185.19060.1484
Electricity (MJ)5.5890191.11525.46261.512251.70801.4780
* Algal broth (kg)6.6203226.38026.47061.791261.24931.7507
Harvesting—FlocculationPotassium sulfate (kg)0.0794 0.0215
Electricity (MJ)7.7385 2.0937
Algal broth (kg)6.6203 1.7912
* Wet biomass (kg)5.4353 1.4706
Harvesting—CentrifugationChitosan (g) 45.9690 12.4374
Electricity (MJ) 12.4740 3.3750
Algal broth (kg) 226.3802 61.2493
* Wet biomass (kg) 5.4353 1.4706
Harvesting—FiltrationElectricity (MJ) 190.2355 51.4700
Algal broth (kg) 6.4706 1.7507
* Wet biomass (kg) 5.4353 1.4706
Oven DryingWet biomass (kg)5.43535.43535.43531.47061.47061.4706
Heat (MJ)62.451662.451662.451616.896916.896916.8969
* Dry biomass (kg)4.62004.62004.62001.25001.25001.2500
Extraction—SolventHexane (g)2.95682.95682.9568
Electricity (MJ)0.33260.33260.3326
Heat (MJ)2.31002.31002.3100
Dry biomass (kg)4.62004.62004.6200
* Algal oil (g)1049.98841049.98841049.9884
* Liquid residue (g)170.0000170.0000170.0000
* Solid residue (g)40.000040.000040.0000
Extraction—MicrowaveElectricity (MJ) 36.749636.749636.7496
Dry biomass (kg) 1.25001.25001.2500
* Algal oil (kg) 1.05001.05001.0500
* Liquid residue (g) 45.995145.995145.9951
* Solid residue (g) 10.822410.822410.8224
TransesterificationAlgal oil (g)1049.98841049.98841049.98841049.98841049.98841049.9884
Methanol (g)124.8787124.8787124.8787124.8787124.8787124.8787
Sodium hydroxide (g)10.487410.487410.487410.487410.487410.4874
Sulfuric acid (g)15.800415.800415.800415.800415.800415.8004
Electricity (MJ)0.16630.16630.16630.16630.16630.1663
Heat (MJ)5.59025.59025.59025.59025.59025.5902
Water (kg)0.13860.13860.13860.13860.13860.1386
* Biodiesel (g)1000.00001000.00001000.00001000.00001000.00001000.0000
* Glycerol (g)113.2825113.2825113.2825113.2825113.2825113.2825
PyrolysisSolid residue (g)40.000040.000040.000010.822410.822410.8224
Heat (MJ)6.41686.41686.41681.73611.73611.7361
* Biochar (g)5.99705.99705.99701.62251.62251.6225
Anaerobic DigestionLiquid residue (g)170.0000170.0000170.000045.995145.995145.9951
Electricity (MJ)54.726054.726054.726014.806614.806614.8066
Heat (MJ)416.2671416.2671416.2671112.6250112.6250112.6250
* Methane (g)19.406419.406419.40645.25065.25065.2506
* Fertilizer (g)69.842869.842869.842818.896618.896618.8966
Combined Heat and PowerMethane (g)19.406419.406419.40645.25065.25065.2506
* Electricity (MJ)0.21560.21560.21560.05830.05830.0583
* Heat (MJ)0.40970.40970.40970.11080.11080.1108
* Output Material; PB—Photobioreactor; FL—Flocculation; CN—Centrifugation; FN—Filtration; SO—Solvent; MI—Microwave.
Table A6. Normalization and weighting factors for ReCiPe 2016 impact categories.
Table A6. Normalization and weighting factors for ReCiPe 2016 impact categories.
CategoryNormalization FactorWeighting Factor
Damage to human health11.2400
Damage to ecosystems1186400
Damage to resource availability 0.000357200
Retrieved from SimaPro ReCiPe 2016 Egalitarian impact assessment method.
Table A7. Life cycle assessment results for 1 kg of biodiesel.
Table A7. Life cycle assessment results for 1 kg of biodiesel.
Impact CategoryOP-FL-SOOP-CN-SOOP-FN-SOOP-FL-MIOP-CN-MIOP-FN-MIPB-FL-SOPB-CN-SOPB-FN-SOPB-FL-MIPB-CN-MIPB-FN-MI
Global Warming, Human health4.42273.900924.87447.03416.39195.88205.74785.226126.19958.35937.71717.2071
Stratospheric ozone depletion1.79 × 10−45.37 × 10−31.42 × 10−31.79 × 10−12.83 × 10−45.36 × 10−32.35 × 10−45.43 × 10−31.48 × 10−31.79 × 10−13.39 × 10−45.42 × 10−3
Ionizing radiation1.14 × 10−47.14 × 10−52.01 × 10−35.75 × 10−41.39 × 10−49.79 × 10−51.33 × 10−49.05 × 10−52.03 × 10−35.94 × 10−41.58 × 10−41.17 × 10−4
Fine particulate matter formation0.56000.49053.33940.96560.89820.83030.74410.67473.52361.14971.08241.0145
Ozone formation, Human health7.96 × 10−47.66 × 10−43.85 × 10−32.81 × 10−31.18 × 10−31.15 × 10−31.05 × 10−31.02 × 10−34.10 × 10−33.06 × 10−31.43 × 10−31.40 × 10−3
Human carcinogenic toxicity2.10921.874912.16354.14943.61773.38862.85352.619112.90784.89364.36194.1329
Human non-carcinogenic toxicity14.717011.9791122.803329.182924.746622.070619.543616.8058127.629934.009529.573226.8973
Water consumption, Human health7.01 × 10−35.32 × 10−42.28 × 10−16.33 × 10−36.93 × 10−35.98 × 10−47.12 × 10−36.44 × 10−42.28 × 10−16.44 × 10−37.04 × 10−37.10 × 10−4
Global Warming, Terrestrial ecosystems0.93810.82745.27631.48941.35361.24541.21881.10805.55691.77001.63431.5260
Ozone formation, Terrestrial ecosystems1.21 × 10−21.20 × 10−25.90 × 10−25.58 × 10−21.79 × 10−21.78 × 10−21.58 × 10−21.57 × 10−26.28 × 10−25.95 × 10−22.17 × 10−22.16 × 10−2
Terrestrial acidification5.06 × 10−24.34 × 10−23.28 × 10−18.29 × 10−27.14 × 10−26.44 × 10−26.51 × 10−25.79 × 10−23.43 × 10−19.74 × 10−28.60 × 10−27.89 × 10−2
Terrestrial ecotoxicity8.60 × 10−43.99 × 10−41.71 × 10−21.33 × 10−31.02 × 10−35.69 × 10−49.86 × 10−45.24 × 10−41.72 × 10−21.45 × 10−31.14 × 10−36.94 × 10−4
Water consumption, Terrestrial ecosystems4.70 × 10−37.04 × 10−41.47 × 10−11.08 × 10−25.02 × 10−31.12 × 10−34.89 × 10−38.93 × 10−41.48 × 10−11.10 × 10−25.21 × 10−31.30 × 10−3
Land use0.043528.53680.3170974.63660.054327.90310.054328.54760.3278974.64740.065127.9139
Global Warming, Freshwater ecosystems2.56 × 10−52.25 × 10−51.44 × 10−44.06 × 10−53.69 × 10−53.39 × 10−53.32 × 10−53.02 × 10−51.52 × 10−44.82 × 10−54.45 × 10−54.16 × 10−5
Freshwater eutrophication9.07 × 10−38.64 × 10−33.43 × 10−21.93 × 10−21.75 × 10−21.71 × 10−2 × 10−21.28 × 10−21.24 × 10−23.80 × 10−22.31 × 10−22.12 × 10−22.08 × 10−2
Freshwater ecotoxicity6.78 × 10−44.20 × 10−41.00 × 10−21.23 × 10−31.01 × 10−37.54 × 10−48.43 × 10−45.85 × 10−41.02 × 10−21.39 × 10−31.17 × 10−39.19 × 10−4
Water consumption, Aquatic ecosystems4.77 × 10−64.60 × 10−61.27 × 10−56.73 × 10−64.85 × 10−64.67 × 10−69.24 × 10−69.07 × 10−61.72 × 10−51.12 × 10−59.31 × 10−69.14 × 10−6
Marine ecotoxicity0.79120.62577.22391.56431.31511.15331.04330.87777.47601.81631.56711.4053
Marine eutrophication1.10 × 10−51.13 × 10−52.20 × 10−53.23 × 10−51.23 × 10−51.26 × 10−52.09 × 10−52.12 × 10−53.19 × 10−54.22 × 10−52.22 × 10−52.25 × 10−5
Mineral resource scarcity4.23 × 10−47.07 × 10−51.31 × 10−21.06 × 10−34.39 × 10−49.48 × 10−54.39 × 10−48.64 × 10−51.31 × 10−21.07 × 10−34.55 × 10−41.10 × 10−4
Fossil resource scarcity2.78 × 10−21.68 × 10−24.22 × 10−14.38 × 10−24.10 × 10−23.02 × 10−23.46 × 10−22.35 × 10−24.28 × 10−15.06 × 10−24.77 × 10−23.69 × 10−2
TOTAL23.696048.3245177.26411019.42738.541362.612731.349455.9779184.91751027.08146.194670.2661
* Output Material; OP—Open Pond; PB—Photobioreactor; FL—Flocculation; CN—Centrifugation; FN—Filtration; SO—Solvent; MI—Microwave.
Table A8. Model indices.
Table A8. Model indices.
IndicesDescription
iIndex of cultivation process alternatives (1…I)
jIndex of harvesting process alternatives (1…J)
kIndex of extraction process alternatives (1…K)
rIndex of recovery processes (1…R)
uIndex of single process units (1…U)
mIndex of recoverable material inputs (1…M)
pIndex of final products (1…P)
tIndex of time periods (1…T)
Table A9. Model parameters.
Table A9. Model parameters.
ParameterDescription
DptDemand of final product p for period t
SPptSelling price of final product p for period t
FCiFixed cost of cultivation process alternative i
FHjFixed cost of harvesting process alternative j
FEkFixed cost of extraction process alternative k
FRrFixed cost of recovery process r
FPuFixed cost of single process unit u
OCitOperating cost of cultivation process alternative i for period t
OHjtOperating cost of harvesting process alternative j for period t
OEktOperating cost of extraction process alternative k for period t
ORrtOperating cost of recovery process r for period t
OPutOperating cost of single process unit u for period t
MCmtPurchase cost of material input m for period t
MKmtPurchase capacity of material input m for period t
ICmtInventory cost of material input m for period t
IKmtInventory capacity of material input m for period t
CCiOutput capacity of cultivation process alternative i
CHjOutput capacity of harvesting process alternative j
CEkOutput capacity of extraction process alternative k
CRrOutput capacity of recovery process r
CPuOutput capacity of single process unit u
YCiOutput yield of cultivation process alternative i per input material
YHjOutput yield of harvesting process alternative j per input material
YEkOutput yield of extraction process alternative k per input material
YRrOutput yield of recovery process r per input material
YPuOutput yield of single process unit u per input material
ECiEnvironmental impact per output of cultivation process alternative i
EHjEnvironmental impact per output of harvesting process alternative j
EEkEnvironmental impact per output of extraction process alternative k
ERrEnvironmental impact per output of recovery process r
EPuEnvironmental impact per output of single process unit u
Table A10. Model variables.
Table A10. Model variables.
VariablesDescription
BCiBinary variable, 1 if cultivation process alternative i is used
BHjBinary variable, 1 if harvesting process alternative j is used
BEkBinary variable, 1 if extraction process alternative k is used
BRrBinary variable, 1 if recovery process r is used
BPuBinary variable, 1 if single process unit u is used
PCitProduction output of cultivation process alternative i for period t
PHjtProduction output of harvesting process alternative j for period t
PEktProduction output of extraction process alternative k for period t
PRrtProduction output of recovery process r for period t
PPutProduction output of single process unit u for period t
MQmPurchase quantity of material input m for period t
BImtBeginning inventory of material input m for period t
EImtEnding inventory of material input m for period t
TOptTotal output of final product p for period t

References

  1. International Energy Agency. Available online: www.iea:reports/world-energy-outlook-2019 (accessed on 2 March 2020).
  2. Morales, M.; Collet, P.; Lardon, L.; Hélias, A.; Steyer, J.; Bernard, O. Life-cycle assessment of microalgal-based biofuel. Biomass Biofuels Biochem. 2019, 507–550. [Google Scholar] [CrossRef]
  3. Khan, S.; Siddique, R.; Sajjad, W.; Nabi, G.; Hayat, K.M.; Duan, P.; Yao, L. Biodiesel production from algae to overcome the energy Crisis. Hayati J. Biosci. 2017, 24, 163–167. [Google Scholar] [CrossRef]
  4. Markou, G.; Vandamme, D.; Muylaert, K. Microalgal and cyanobacterial cultivation: The supply of nutrients. Water Res. 2014, 65, 186–202. [Google Scholar] [CrossRef] [Green Version]
  5. Karamerou, E.E.; Webb, C. Cultivation modes for microbial oil production using oleaginous yeasts—A review. Biochem. Eng. 2019, 151. [Google Scholar] [CrossRef]
  6. Brennan, L.; Owende, P. Biofuels from microalgae—a review of technologies for production, processing, and extractions of biofuels and co-products. Renew. Sustain. Energy Rev. 2010, 14, 557–577. [Google Scholar] [CrossRef]
  7. Cheirsilp, B.; Srinuanpan, S.; Mandik, Y.I. Chapter 6—Efficient harvesting of microalgal biomass and direct conversion of microalgal lipids into biodiesel. In Microalgae Cultivation for Biofuels Production; Yousuf, A., Ed.; Academic Press: Cambridge, MA, USA, 2020; pp. 83–96. [Google Scholar]
  8. Chisti, Y. Biodiesel from microalgae. Biotechnol. Adv. 2007, 25, 294–306. [Google Scholar] [CrossRef]
  9. De Assis, T.C.; Calijuri, M.L.; Assemany, P.P.; Pereira, A.S.A.D.P.; Martins, M.A. Using atmospheric emissions as CO2 source in the cultivation of microalgae: Productivity and economic viability. J. Clean. Prod. 2019, 215, 1160–1169. [Google Scholar] [CrossRef]
  10. Culaba, A.B.; Ching, P.L.; San Juan, J.L.; Mayol, A.P.; Sybingco, E.; Ubando, A.T. A dynamic sustainability assessment of algal biorefineries for Biofuel production. In Proceedings of the 2019 IEEE 11th International Conference on Humanoid, Nanotechnology, Information Technology, Communication and Control, Environment, and Management (HNICEM), Laoag, Philippines, 29 November–1 December 2019. [Google Scholar] [CrossRef]
  11. Mitra, M.; Mishra, S. Multiproduct biorefinery from Arthrospira spp. towards zero waste: Current status and future trends. Bioresour. Technol. 2019, 291. [Google Scholar] [CrossRef] [PubMed]
  12. Mohan, S.; Hemalatha, M.; Chakraborty, D.; Chatterjee, S.; Ranadheer, P.; Kona, R. Algal biorefinery models with self-sustainable closed loop approach: Trends and prospective for blue-bioeconomy. Bioresour. Technol. 2019. [Google Scholar] [CrossRef]
  13. De Bhowmick, G.; Sarmah, A.K.; Sen, R. Zero-waste algal biorefinery for bioenergy and biochar: A green leap towards achieving energy and environmental sustainability. Sci. Total Environ. 2019, 650. [Google Scholar] [CrossRef]
  14. Hemalatha, M.; Sravan, J.; Min, B.; Mohan, S. Microalgae-biorefinery with cascading resource recovery design associated to dairy wastewater treatment. Bioresour. Technol. 2019, 284, 424–429. [Google Scholar] [CrossRef]
  15. Wu, W.; Chang, J. Integrated algal biorefineries from process systems engineering aspects: A review. Bioresour. Technol. 2019, 291. [Google Scholar] [CrossRef]
  16. Sy, C.L.; Ubando, A.T.; Aviso, K.B.; Tan, R.R. Multi-objective target oriented robust optimization for the design of an integrated biorefinery. J. Clean. Prod. 2018, 170, 496–509. [Google Scholar] [CrossRef]
  17. Culaba, A.B.; San Juan, J.L.; Ching, P.L.; Mayol, A.P.; Sybingco, E.; Ubando, A.T. Optimal synthesis of algal biorefineries for Biofuel production based on techno-economic and environmental efficiency. In Proceedings of the 2019 IEEE 11th International Conference on Humanoid, Nanotechnology, Information Technology, Communication and Control, Environment, and Management (HNICEM), Laoag, Philippines, 29 November–1 December 2019. [Google Scholar] [CrossRef]
  18. García Prieto, C.V.; Ramos, F.D.; Estrada, V.; Villar, M.A.; Diaz, M.S. Optimization of an integrated algae-based biorefinery for the production of biodiesel, astaxanthin and PHB. Energy 2017, 139, 1159–1172. [Google Scholar] [CrossRef]
  19. Gong, J.; You, F. Optimal design and synthesis of algal Biorefinery processes for biological carbon sequestration and utilization with zero direct greenhouse gas emissions: MINLP model and global optimization algorithm. Ind. Eng. Chem. Res. 2014, 53, 1563–1579. [Google Scholar] [CrossRef]
  20. Ching, P.M.; Mayol, A.P.; San Juan, J.L.; Calapatia, A.M.; So, R.H.; Sy, C.L.; Ubando, A.T.; Culaba, A.B. AI methods for modeling the vacuum drying characteristics of Chlorococcum infusionum for algal Biofuel production. Process Integr. Optim. Sustain. 2021. [Google Scholar] [CrossRef]
  21. San Juan, J.L.; Mayol, A.P.; Sybingco, E.; Ubando, A.T.; Culaba, A.B.; Chen, W.H.; Chang, J.S. A scheduling and planning algorithm for microalgal cultivation and harvesting for biofuel production. IOP Conf. Ser. Earth Environ. Sci. 2020, 463. [Google Scholar] [CrossRef] [Green Version]
  22. Galanopoulos, C.; Kenkel, P.; Zondervan, E. Superstructure optimization of an integrated algae biorefinery. Comput. Chem. Eng. 2019, 130. [Google Scholar] [CrossRef]
  23. Caligan, C.J.; Garcia, M.M.; Mitra, J.L.; Mayol, A.P.; San Juan, J.L.; Culaba, A.B. Multi-objective optimization of water exchanges between a wastewater treatment facility and algal biofuel production plant. IOP Conf. Ser. Earth Environ. Sci. 2020, 463. [Google Scholar] [CrossRef]
  24. San Juan, J.L.; Caligan, C.J.; Garcia, M.M.; Mitra, J.L.; Mayol, A.P.; Sy, C.L.; Ubando, A.T.; Culaba, A.B. Multi-objective optimization of an integrated algal and sludge-based Bioenergy Park and wastewater treatment system. Sustainability 2020, 12, 7793. [Google Scholar] [CrossRef]
  25. Čuček, L.; Martín, M.; Grossman, I.E.; Kravanja, Z. Multi-period synthesis of optimally integrated biomass and bioenergy supply network. Comput. Chem. Eng. 2014, 66, 57–70. [Google Scholar] [CrossRef]
  26. Chowdhury, R.; Sadhukhan, J.; Traverso, M.; Keen, P.L. Effects of residence time on life cycle assessment of bioenergy production from dairy manure. Bioresour. Technol. Rep. 2018, 4, 57–65. [Google Scholar] [CrossRef]
  27. Barlow, J.; Sims, R.C.; Quinn, J.C. Techno-economic and life-cycle assessment of an attached growth algal biorefinery. Bioresour. Technol. 2016, 220, 360–368. [Google Scholar] [CrossRef]
  28. Wu, W.; Lin, K.; Chang, J. Economic and life-cycle greenhouse gas optimization of microalgae-tobiofuels chains. Bioresour. Technol. 2018, 267, 550–559. [Google Scholar] [CrossRef] [PubMed]
  29. Bussa, M.; Eisen, A.; Zollfrank, C.; Röder, H. Life cycle assessment of microalgae products: State of the art and their potential for the production of polylactid acid. J. Clean. Prod. 2018, 213, 1299–1312. [Google Scholar] [CrossRef]
  30. Giwa, A.; Adeyemi, I.; Dindi, A.; García-Bañoz Lopez, C.; Lopresto, C.; Curcio, S.; Chakraborty, S. Techno-economic assessment of the sustainability of an integrated biorefinery from microalgae and Jatropha: A review and case study. Renew. Sustain. Energy Rev. 2018, 88, 239–257. [Google Scholar] [CrossRef]
  31. Gifuni, I.; Pollio, A.; Safi, C.; Marzocchella, A.; Olivieri, G. Current Bottlenecks and Challenges of the Microalgal Biorefinery. Trends Biotechnol. 2019, 37, 242–252. [Google Scholar] [CrossRef] [PubMed]
  32. Branyikova, I.; Prochazkova, G.; Potocar, T.; Jezkova, Z.; Branyik, T. Harvesting of Microalgae by Flocculation. Fermentation 2018, 4, 93. [Google Scholar] [CrossRef] [Green Version]
  33. Su, Y.; Song, K.; Zhang, P.; Su, Y.; Cheng, J.; Chen, X. Progress of microalgae biofuel’s commercialization. Renew. Sustain. Energy Rev. 2017, 74, 402–411. [Google Scholar] [CrossRef]
  34. Yew, G.Y.; Lee, S.Y.; Show, P.L.; Tao, Y.; Law, C.L.; Nguyen, T.T.C.; Chang, J. Recent advances in algae biodiesel production: From upstream cultivation to downstream processing. Bioresour. Technol. Rep. 2019, 7. [Google Scholar] [CrossRef]
  35. Dasan, Y.K.; Lam, M.K.; Yusup, S.; Lim, J.W.; Lee, K.T. Life cycle evaluation of microalgae biofuels production: Effect of cultivation system on energy, carbon emission and cost balance analysis. Sci. Total Environ. 2019, 688, 112–128. [Google Scholar] [CrossRef] [PubMed]
  36. Biller, P.; Friedman, C.; Ross, A.B. Hydrothermal microwave processing of microalgae as a pre-treatment and extraction technique for bio-fuels and bio-products. Bioresour. Technol. 2013, 136, 188–195. [Google Scholar] [CrossRef] [PubMed]
  37. Gnansounou, E.; Raman, J.K. Chapter 7—Life Cycle Assessment of Algal Biorefinery. In Life-Cycle Assessment of Biorefineries; Elsevier: Amsterdam, The Netherlands, 2017; pp. 199–219. [Google Scholar]
  38. Solis, C.A.; Mayol, A.P.; San Juan, J.G.; Ubando, A.T.; Culaba, A.B. Multi-objective optimal synthesis of algal biorefineries toward a sustainable circular bioeconomy. IOP Conf. Ser. Earth Environ. Sci. 2020, 463. [Google Scholar] [CrossRef]
Figure 1. Algal biorefinery process flow.
Figure 1. Algal biorefinery process flow.
Energies 14 01416 g001
Figure 2. Life cycle environmental impacts of each upstream process path for 1 kg biodiesel (OP—open pond; PB—photobioreactor; FL—flocculation; CN—centrifugation; FN—filtration; OD—oven drying).
Figure 2. Life cycle environmental impacts of each upstream process path for 1 kg biodiesel (OP—open pond; PB—photobioreactor; FL—flocculation; CN—centrifugation; FN—filtration; OD—oven drying).
Energies 14 01416 g002
Figure 3. Life cycle environmental impacts of each downstream process path for 1 kg biodiesel.
Figure 3. Life cycle environmental impacts of each downstream process path for 1 kg biodiesel.
Energies 14 01416 g003
Figure 4. Life cycle environmental impacts of each process path for 1 kg biodiesel.
Figure 4. Life cycle environmental impacts of each process path for 1 kg biodiesel.
Energies 14 01416 g004
Figure 5. Behavior of purchase amounts for each input material over time.
Figure 5. Behavior of purchase amounts for each input material over time.
Energies 14 01416 g005
Figure 6. Behavior of inventory levels for each input material over time.
Figure 6. Behavior of inventory levels for each input material over time.
Energies 14 01416 g006
Figure 7. Results under varying process unit efficiencies: (a) Profit and environmental impact optimal values under varying process efficiencies; (b) objective function efficiency values under varying process efficiencies.
Figure 7. Results under varying process unit efficiencies: (a) Profit and environmental impact optimal values under varying process efficiencies; (b) objective function efficiency values under varying process efficiencies.
Energies 14 01416 g007
Figure 8. Results under varying demand parameters: (a) Profit and environmental impact optimal values under varying demand; (b) objective function efficiency values under varying demand.
Figure 8. Results under varying demand parameters: (a) Profit and environmental impact optimal values under varying demand; (b) objective function efficiency values under varying demand.
Energies 14 01416 g008
Table 1. Comparison of the different process unit alternatives.
Table 1. Comparison of the different process unit alternatives.
ProcessAlternativeAdvantagesDisadvantagesRef
CultivationOpen pondLow installation and operating costsVarying culture conditions[31]
PhotobioreactorHigh control over culture parametersHigh capital and operating costs
HarvestingFlocculationOver 90% recovery; possibly low costChemical contamination; longer settle time[32]
Filtration70–90% recoveryHigh capital and operating costs
CentrifugationOver 90% recovery; can utilize most algae speciesEnergy-intensive
ExtractionMicrowaveAround 90–95% yield with great qualityHigh operating cost; energy-intensive[33]
SolventLow-cost solventsLarge amount of solvent needed
Table 2. Summary of single objective optimization results.
Table 2. Summary of single objective optimization results.
Maximized ProfitMinimized Impact
Overall Profit—10 Periods (USD)5,923,239.93−36,518.68
Annual Impact (kPt)8,682.01592.40
Annual Biodiesel Output (kg)48,977.79 kg25,000 kg
Table 3. Summary of multi-objective optimal results.
Table 3. Summary of multi-objective optimal results.
Optimal ValueEfficiency
Overall Profit—10 Periods (USD)5,487,668.6992.69%
Annual Impact (kPt)1183.6592.69%
Annual Biodiesel Output (kg)49,951.43-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Solis, C.M.A.; San Juan, J.L.G.; Mayol, A.P.; Sy, C.L.; Ubando, A.T.; Culaba, A.B. A Multi-Objective Life Cycle Optimization Model of an Integrated Algal Biorefinery toward a Sustainable Circular Bioeconomy Considering Resource Recirculation. Energies 2021, 14, 1416. https://doi.org/10.3390/en14051416

AMA Style

Solis CMA, San Juan JLG, Mayol AP, Sy CL, Ubando AT, Culaba AB. A Multi-Objective Life Cycle Optimization Model of an Integrated Algal Biorefinery toward a Sustainable Circular Bioeconomy Considering Resource Recirculation. Energies. 2021; 14(5):1416. https://doi.org/10.3390/en14051416

Chicago/Turabian Style

Solis, Celine Marie A., Jayne Lois G. San Juan, Andres Philip Mayol, Charlle L. Sy, Aristotle T. Ubando, and Alvin B. Culaba. 2021. "A Multi-Objective Life Cycle Optimization Model of an Integrated Algal Biorefinery toward a Sustainable Circular Bioeconomy Considering Resource Recirculation" Energies 14, no. 5: 1416. https://doi.org/10.3390/en14051416

APA Style

Solis, C. M. A., San Juan, J. L. G., Mayol, A. P., Sy, C. L., Ubando, A. T., & Culaba, A. B. (2021). A Multi-Objective Life Cycle Optimization Model of an Integrated Algal Biorefinery toward a Sustainable Circular Bioeconomy Considering Resource Recirculation. Energies, 14(5), 1416. https://doi.org/10.3390/en14051416

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop