Next Article in Journal
Ultra-High-Performance Micellar Liquid Chromatography Comparing Tween 20 and Tween 40 for the Determination of Hydroxycinnamic Acids
Previous Article in Journal
Research and Application of In Situ Sample-Processing Methods for Rapid Simultaneous Detection of Pyrethroid Pesticides in Vegetables
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Screening for New Efficient and Sustainable-by-Design Solvents to Assist the Extractive Fermentation of Glucose to Bioethanol Fuels

by
Effie Marcoulaki
* and
Pantelis Baxevanidis
Institute of Nuclear and Radiological Sciences and Technology, Energy and Safety, National Centre for Scientific Research “Demokritos”, 15341 Aghia Paraskevi, Greece
*
Author to whom correspondence should be addressed.
Separations 2022, 9(3), 60; https://doi.org/10.3390/separations9030060
Submission received: 20 December 2021 / Revised: 1 February 2022 / Accepted: 15 February 2022 / Published: 25 February 2022

Abstract

:
The production of bioethanol fuels using extractive fermentation increases the efficiency of the bioconversion reaction by reducing the toxic product inhibition. The choice of appropriate solvents to remove the bioethanol product without inhibiting the fermentation is important to enable industrial scale application. This work applies computer-aided molecular design technologies to systematically screen a wide variety of candidate solvents to enhance the separation, also considering the microorganisms that perform the fermentation. The performance of the candidates was evaluated using a rigorous process simulator for extractive fermentation, assisted by functional group-contribution (QSPR/QSAR) models for the prediction of various solvent properties, including toxicity and life cycle impacts. The solvent designs generated through this approach can provide powerful insights on the kind of molecular structures and functionalities that satisfy the process objectives and constraints, as well the desired sustainability features.

1. Introduction

The European Green Deal (EUGD) (COM(2019) 640 final) aims to make Europe the world’s first climate-neutral continent by 2050, and deliver benefits from sustainable green transition to European citizens and businesses [1]. EUGD raised the target for the renewable energy share to 40% of the European energy production by 2030, compared to the current regulations [2]. Despite the expected market penetration of electricity/hydrogen-powered vehicles and possible energy savings, biofuels and bioliquids will have a significant share in the transportation sector, due to the large long-term demand for dense liquid fuels for aviation, ocean shipping, and long-haul trucking [3]. The expected increase in biofuel production to meet these demands requires assessment of their life cycle impacts, and the implementation of more sustainable biofuel production systems [4].
Bioethanol is a worldwide leader in the transition to more sustainable liquid fuels, used since the 1970s to partially replace gasoline [5]. Its advantages compared to other biofuel alternatives include its cost, global warming potential, accessibility, and socioeconomic impact [6]. An additional advantage is that variable feedstocks can be used for bioethanol production, mainly sugars, starch, and lignocellulosic biomass (e.g., crops, aquatic plants, forest materials, or agricultural residues). The process of converting biomass to ethanol depends on the biochemical composition of this feedstock. The conversion parameters affect the total production costs (including costs for raw material, operation, and capital depreciation) and the cradle-to-grave environmental impacts [7,8].
The extractive fermentation (EF) process combines fermentation with a primary liquid–liquid (LL) extractive separation stage. Removing the produced ethanol and other inhibitory compounds in situ eliminates inhibitions caused by ethanol and increase the ethanol productivity. Recent research acknowledges the potential of this fermentation strategy to reduce operation costs, increase conversion efficiency, and reduce the inhibition of products and by-products. Lemos et al. [9] reported significant productivity increase compared to conventional extraction techniques, demonstrating the potential of EF for use at larger scales. More work is still needed, however, to make EF suitable for industrial application, as shown in recent experimental studies [10,11] on fed batch EF schemes.
Similar to conventional LL extraction processes, EF uses a solvent to selectively extract a desired fermentation product from a mixture. In addition to their extraction efficiency, EF solvents must also be biocompatible with the microorganisms performing the fermentation. When the known solvents (or other tabulated materials) fail to meet the EF process requirements, alternative solvent designs must be sought.
Borhani and Wang [12] classified solvent screening approaches into four different categories according to their reliance on experiments, process and equilibrium models, prediction models, and finally systematic solvent design. Considering the enormous number of candidate materials, trial-and-error using experiments or calculations would be very time-consuming and expensive to apply. Therefore, systematic methods for the optimal selection and molecular design of solvents by means of computer-aided molecular design (CAMD) are becoming increasingly popular for a variety of solvent design applications [13].
Even until recently, contributions in the field of solvent design considered only the solvent properties related to process efficiency. In CAMD, the solvent properties are predicted using quantitative structure-property/activity models (QSPRs/QSARs), to reverse-engineer solvents with desired/optimal properties. For organic solvents, these QSPR/QSAR models usually take the form of functional group-contribution (GC) correlations [14], which provide a good approximation of the property values of unknown materials, using a set of tabulated group-contribution coefficients. Linear GC models are usually applied for [15] the prediction of the solvent physical properties, while UNIFAC models are the dominant choice for accurate prediction of activity coefficients and phase equilibria [12,16].
Given the rising interest in sustainable chemical manufacturing, recent research on CAMD of separation solvents additionally considers health, safety, and environmental aspects [16,17]. For instance, Khor et al. [18] applied a fuzzy optimization approach to the design of palm oil extraction solvents to substitute hexane, considering a set of target properties along with health and safety aspects. Papadopoulos et al. [19] proposed a design framework for phase-change solvents used in CO2 capture, linked to databases for the evaluation of various sustainability-related indicators. Ten et al. [20] developed a CAMD framework for optimal liquid–liquid extraction solvents considering solvent properties, as well as health and safety criteria using Analytical Hierarchy Process (AHP). Baxevanidis et al. [21] proposed new GC models for lifecycle assessment of environmental impacts and applied them on the CAMD of LL extraction solvents. Solvent optimization studies usually consider a set of solvent properties to improve the process efficiency, e.g., solvent selectivity, solvent capacity, and various primary properties. Detailed process models (if/when applied) are usually considered after CAMD, as a further assessment of the CAMD results [22].
Despite the large progress in the design of separation solvents, challenges remain for the industrial use of the CAMD approach, due to the variety of process-specific properties and the need to capture accurately the desirable solvent profile [17]. This work combines state-of-the art CAMD optimization tools with a rigorous EF process simulator and recently developed lifecycle assessment (LCA) models, to design new, efficient, and sustainable solvents for EF of glucose to bioethanol fuels. The CAMD tool employed here is based on the original implementation of Marcoulaki and Kokossis [23], demonstrated on the design of optimal separation-enhancing organic materials [24]. Since then, the tool has been used extensively on a variety of solvent design applications. These include solvents for reactive separation processes [25], organic Rankine cycle solvents [26], CO2 capture solvents [27], also considering sustainability aspects [19,21].
Section 2 outlines the CAMD optimizer of Marcoulaki and Kokossis [23]. Section 3 gives the mathematical formulation of the solvent design problem considered here, including the EF process simulator and the GC models required for the prediction of the desirable solvent properties. Section 4 presents solvent design results for different sets of design constraints and objectives. Section 5 discusses the results and Section 6 concludes this work.

2. CAMD Optimization Algorithm

Figure 1 outlines the CAMD optimization algorithm used here. The optimization starts from an initial guess for the functional groups comprising the solvent, e.g., {2xCH3, 3xCH2, 1xCH, 1xOH}. The search proceeds via random perturbations, applied on the initial solvent to generate a new solvent. These perturbations consist of modifying the groups comprising the initial solvent (i.e., addition of new groups, removal of existing groups, replacement of existing groups with new ones). The perturbation strategies and the restoration procedures to ensure structural feasibility of the new solvent are according to Marcoulaki and Kokossis [23].
The new solvent is then compared to the initial one using a set of stochastic criteria according to the Simulated Annealing optimization algorithm [28] as adjusted by Marcoulaki and Kokossis [29]. The comparison is in terms of the performances of the two solvents, based on a set of user-defined objective functions and constraints. For the evaluation of objectives and constraints, the CAMD optimizer is linked to process simulators, GC models, databases, etc., depending on the process task (see Section 3). If the new solvent is accepted, it replaces the initial solvent, otherwise we continue with the initial solvent. This procedure runs iteratively until a set of termination criteria are met, and returns the final solvent, i.e., a final set of functional groups. Note that the typical solution space for this type of problems stretches along trillions of molecular configurations.
Since the mapping between molecules and functional groups is usually many to one, we may have more than one way of connecting the groups together. For example, {2xCH3, 3xCH2, 1xCH, 1xOH} can be any of the following molecules: CCCC(C)CO, CCC(C)CCO, CC(C)CCCO, CCC(CC)CO, CCCCC(O)C, or CCCC(O)CC. Even though the algorithm always generates structurally feasible sets of groups, certain design configurations may be uncommon and potentially hard to synthesize. If this is the case, it is possible to include a complexity assessment in the solvent performance calculations. Additionally or alternatively, one can control the generation of extremely large molecules by limiting the total number of groups, or avoid very complex (possibly unstable) aromatics by limiting the number of rings or the number of substituents per ring, etc. However, setting a very low number of allowable groups may trap the optimizer in a locally optimal domain with limited chance of escaping.
Due to the stochastic nature of the CAMD optimizer, each time we run the tool using the same setup, we may obtain a slightly different final design. The convergence quality of the optimizer is evidenced by the low deviations between the objective function values of the final designs. We also need to ensure that the final objective function value in each run is adequately close (if not the same) with the best objective function value observed during the run.

3. Formulation of the Solvent Design Problem

A good EF solvent has several desirable characteristics. These include biocompatibility, good extraction performance, and low solubility in the aqueous phase, and that the solvent density is adequately different to that of the aqueous phase and that the solvent is chemically stable. Additionally, the solvent selection should contribute to the sustainability of the process, as well as the economic and performance aspects of the process [9]. Based on these criteria, this section presents the mathematical models used here to evaluate the performance of the different solvents generated during the search. The models include a rigorous process simulator to calculate the process objectives, and GC models for the prediction of the various properties appearing in the solvent constraints.

3.1. EF Process Simulator

Several simulation models of variable complexity are available for the mathematical formulation of the EF process, to assist the assessment of candidate extraction solvents [30]. The efficiency of the CAMD tool used here is not affected by the complexity of the process simulator, so we are able to use the rigorous model proposed by Kollerup and Daugulis [31] for continuous stirred tank fermenters. This model is suitable for ethanol production via EF and calculates the extraction efficiency and the productivities in the nutrient and the extract phases.

3.1.1. Simulation Model Assumptions/Requirements

The simulation model assumes the following:
  • The broth density, ρ1, is a linear function of the ethanol and glucose concentration in the aqueous phase. Note that linear mixing rules are the default assumption.
  • The bioreaction system produces only cells, ethanol, and carbon dioxide, with yield coefficients YX/S, YP/S, and YC/S, respectively. There are no cells present in the extract phase.
  • The specific growth rate, μ, follows the modified Monod equation:
    μ = μ max S 1 K S + S 1 K P K P + P 1
    where μmax denotes the maximum specific growth rate, KP and KS are the product inhibition and the substrate saturation constants, respectively; S1 denotes the effluent substrate concentration in the aqueous phase; and P1 is the effluent product (ethanol) concentration in the aqueous phase.
The simulator model also requires that the solvent is completely biocompatible, and the only function of the solvent with respect to the fermentation system is the extraction of the bioreaction product from the nutrient. It is additionally required that there are negligible solvent losses to the nutrient phase and negligible substrate present in the extract phase. These are solvent-specific requirements enforced under the constraints of Section 3.2.

3.1.2. Simulation Model Formulation

The model inputs are the product distribution coefficient (MP), the influent substrate concentration in the aqueous phase ( S 1 0 ), and the dilution rate in the solvent phase (DE). The calculation of MP is discussed in Section 3.2. The S 1 0 , S1, and DE can be fixed, to allow comparison between the candidate solvents under typical fermentation conditions (see Section 4).
The analytical expressions used here in order to calculate P1 are the following:
A P 1 2 + B P 1 + C = 0
where the A, B, and C terms are formulated as:
A = M P { W S K P ρ P β D E + W Y ( D E W S K P ρ P ) }
B = K P { W S [ M P ρ P ( ρ W +   S 1   ( α Y C / S ) ) β ] M P D E + + W Y [ M P D E W S ( M P Y P / S S 1 ρ P 1 ) ] }
C = W S K P { S 1 S 1 0 [ ρ W + S 1 0 ( α Y C / S ) ] [ ρ W + S 1 ( α Y C / S ) ] }
where
W S = μ max S 1 / ( K S + S 1 )
W Y = ρ W + S 1 0   ( α Y C / S ) Y P / S S 1 0
Note that, Equations (2)–(7) are solved for each solvent generated during the search, using available solvers for non-linear equation systems.

3.1.3. Calculation of EF Process Objectives

Having estimated P1, it is possible to calculate the various process objectives typically used in EF solvent selection, i.e., the ethanol productivities, the substrate conversion, and the extraction efficiency.
The total ethanol productivity is the sum of productivities in the nutrient (PDW) and the extract (PDE) phases:
P D T = P D W + P D E
The ethanol productivity in the nutrient, PDW, depends on P1 as follows:
P D W = P 1 D
where the effluent infinite dilution rate in the aqueous phase, D’, is equal to μ of Equation (1).
The ethanol productivity in the extract, PDE, is given as:
P D E = P 2 D E
where P2 is the product concentration in the extract phase, i.e., P 2 = M P P 1 ; and DE’ denotes the effluent dilution rate in the extract phase calculated as:
D E = D E [ 1 + ( ρ P / ( M P P 1 ) 1 ) 1 ]
Using the above equations, the total productivity finally becomes:
P D T = P D W + P D E = P 1 ( D + M P D E )
The substrate conversion is calculated as follows:
C S = 1 D S 1 D S 1 0
where D is the influent infinite dilution rate in the aqueous phase, calculated as:
D = D ( S 1 + P 1 Y P / S ) + M P D E P 1 Y P / S S 1 0 ( 1 M P P 1 / ρ P )
The extraction efficiency (EE) is defined as the fraction of PDE over PDT, therefore:
E E = D D + M P D E

3.2. Constraints for EF Solvents

This section discusses the desirable properties for an appropriate EF solvent, and the models used here for property prediction. Particular emphasis is on the biocompatibility, the product distribution coefficient, and the solvent losses, to satisfy the requirements of the process simulator of Section 3.1.

3.2.1. Distribution Coefficient of the Product to Be Removed (High)

The MP is the main factor to affect the separation driving force, hence the main criterion in selecting extraction solvents [32]. In the simulator of Section 3.1, MP participates in the calculations of productivities and EE. According to Pretel et al. [33], the MP for LL extraction is predicted at infinite dilution conditions and at the expected process operating temperature as follows:
M P = γ A , B γ A , E M W B M W E  
where γ i , j denote the activity coefficients of mixture component i in component j; MWi is the molecular weight of component i; and A, B, and E are the solute, the aqueous feed, and the solvent/extract, respectively. Note that the actual MP value is taken here as a fixed fraction, η, of M P . Higher accuracy in MP could be achieved, but this would not alter the relative ranking of the obtained solvent designs.
The most common GC models for the prediction of γ’s are the UNIFAC model and its modifications [34]. Following the usual practice for solvent design, this work applies the original UNIFAC model using the VLE coefficients. The coefficient values for different pairs of functional groups can be found in Reid et al. [35].

3.2.2. Solvent Losses to the Aqueous Phase (Low)

Low solvent solubility to the aqueous phase leads to low solvent losses, and therefore reduces the required solvent makeup and the cost for downstream treatment of the aqueous effluent. In addition, the EF process simulator of Section 3.1 assumes negligible solvent losses to the aqueous phase. Solvent losses, S L , are predicted at infinite dilution conditions and at the expected process operating conditions as follows [33]:
S L = 1 γ E , B M W E M W B
Activity coefficients are predicted using the UNIFAC model as explained above.

3.2.3. Solvent Selectivity (High)

High selectivity of the solvent towards the solute is very important in LL extraction processes, and usually completes and/or consolidates the requirement for high MP. The solvent selectivity, S E , is predicted at infinite dilution conditions and at the expected process operating conditions as follows [33]:
S E = γ B , E γ A , E M W A M W B  
Activity coefficients are predicted using the UNIFAC model as explained above.

3.2.4. Distribution Coefficient for Essential Nutrients/Byproducts (Low)

The loss of substrate increases the investment/processing costs, and reduces the process efficiency, while the presence of substrate in the extract stream may complicate the recovery/purification stages following the EF stage. In addition, the EF process simulator (see Section 3.1) assumes negligible substrate concentration in the extract phase. Additionally, we may want to limit the mass transfer of reaction byproducts (e.g., acetic acid) from the nutrient to the extract phase. The distribution coefficient of component S (denoting glucose or a byproduct) is given as:
M S = γ S , B γ S , E M W B M W E  
Activity coefficients are again predicted using UNIFAC VLE as explained above.

3.2.5. Ease of Solvent—Solute Separation (High)

Solvent and solute recovery from the extract mixture is preferably performed by simple distillation. A substantial difference in the boiling point temperatures of the two components can guarantee the absence of azeotropic formations and the ease of separation. Boiling point temperature predictions are according to the well-established GC model of Constantinou and Gani [36], who provide the model formulation and the table of group-contribution coefficients.

3.2.6. Solvent Biocompatibility (High)

Solvent biocompatibility is the most important criterion in extractive fermentation [9]. Biocompatibility is in terms of the toxicity effects of the solvent on the microorganisms that perform the fermentation, so that their population size remains intact. It can be accounted by means of lethal concentration standards [37], or using the solvent octanol/water partition coefficient [31]. Trevizo and Nirmalakhandan [38] investigated the relation between solubilities and toxicity measures. The model used here is the GC model proposed by Gao et al. [39] for lethal concentrations, who provide the model formulation and the table of group-contribution coefficients.

3.2.7. Solvent Environmental Impact (Low)

The European chemicals strategy for sustainability aims to consider global environmental impacts as early as possible in the material design process [40]. There are many aspects of sustainability that can be considered over the lifecycle of a new chemical. We herein consider three climate change indicators: the Global Warming Potential (GWP), the Cumulative Energy Demand (CED), and the EcoIndicator 99 (EI99). Baxevanidis et al. [21] provide the model formulations for these three LCA indices and the tables of group-contribution coefficients. Additional sustainability criteria can be considered, if they can be formulated using the same GC approach.

3.2.8. Other Physical Properties of the Solvent

Other property requirements include that the material should be liquid, thermally stable, and chemically inert at the process conditions. In addition, a large density difference between the nutrient and extract phases enhances immiscibility. All these are predicted using well-established GC models, which are generally available for a wide range of pure-component properties, such as critical properties, phase transition enthalpies, phase change temperatures, heat capacity, and viscosity. Constantinou and Gani [36] provide the model formulations and the tables of group-contribution coefficients for these predictions. These models could also be used to calculate flowrates, equipment size, conditions, etc., but these are beyond the scope of the present study.

3.2.9. Solvent Cost (Low)

Although the solvent cost is certainly extremely important in the selection of materials for large-scale industrial use, it may depend on many external factors making it difficult to estimate the market price of a new material, and this is outside the scope of this study.

4. EF Solvent Design Case Study

4.1. Design Case Parameters

Table 1 gives the parameters for the ethanol-water-glucose system used in the EF simulation model of Section 3.1
Table 2 gives the design constraints used here, according to the desired properties for the EF solvent discussed in Section 3.2. The bounds on S L and M S enforce the requirement for negligible solvent losses to the nutrient phase (<0.5%) and low substrate distribution in the extract phase (<10%), respectively. The bound on S E sets the selectivity above 200%. The M P lower bound is very relaxed, since M P is controlled via the objective for maximum ethanol productivity. The bounds on the boiling point temperature, Tb,E, allow recovery of the solvent via simple distillation, and also help to avoid solvents with subambient or extremely high Tb,E. We also require that the solvent is liquid at the process conditions, and that the solvent density, ρΕ, is significantly different from the broth density, ρ1. The two sets of bounds for the LCA indices are taken from Baxevanidis et al. [21]. Finally, we consider a range of tE,max values (from 1.5 to 10 mol/L) to study the effect of tE on the final designs.
As mentioned in Section 3.1, we use fixed values for the influent substrate concentration in the aqueous phase, S 1 0 , and the dilution rate for the solvent phase, DE, to run the simulator. Table 3 gives the data for the operation of the EF process, assuming 12 design cases, denoted as A1 to A12. The S 1 0 and DE values used here are within the range investigated by Kollerup and Daugulis [31].
According to the availability of tabulated group-contribution coefficients to predict the solvent properties, the functional groups used here are 57 listed in Baxevanidis et al. [21]. Each solvent design is allowed to have up to 40 functional groups and 5 aromatic rings. These numbers are adequately high to facilitate an extensive screening of solvent designs. The optimizer setup is according to Marcoulaki and Kokossis [23,29].

4.2. Solvent Design Results

The solvent design objective considered here is to maximize the ethanol productivity in the extract phase (Equation (10)). The objective is applied to each one of the 12 design cases of Table 3, using 6 different upper bounds for tE per design case (see Table 2) to study the effect of the toxicity constraint on the productivity of the EF solvents. Having a sample of 10 runs per combination of tE,max and design case, we end up with 720 computational experiments, which generated a set of 11 final solvent designs, denoted as D1 to D11. So, some of these solvents work well over a range of operating conditions.
Table 4 presents the solvents along with their property values predicted using the models of Section 3. The solvents are ranked according to tE, and they obey the constraints of Table 2 for S E , M P , Sl, Tb,E, Tf,E, and ρE. Regarding the LCA indices, different colors in the last three columns of the table indicate how the solutions conform to the GWP, CED, and EI99 bounds. We observe that tE ranges between 1.69 and 4.51, and cannot go below 1.69 without other constraints being violated. However, the system could tolerate higher toxicities, since the allowable solvent losses (see Table 2) ensure that the solvent amount present in the nutrient phase is extremely small. In terms of the bounds on LCA indices, only D1 satisfies all the low LCA bounds. All the designs satisfy the lower GWP bound. Most of the designs satisfy the lower CED bound and the rest satisfy the higher CED bound. EI99 is the harder to satisfy, with the majority of designs lying between the two EI99 bounds, and designs D7, D8, and D9 violating the upper EI99 bound.
As explained in Section 2, the output of the CAMD tool is actually a set of functional groups, e.g., D2 = {1xCH3, 1xCH2COO, 3xCH2, 1xCH2CN}, so the solvent SMILES shown in Table 4 are reasonable design interpretations using these sets. For instance, D2 could have the ester group closer to the nitrile group (CCCCCC(=O)OCCCCC#N), but usually the farther apart these groups are, the more stable the molecule might be. The proper consideration of stability or inertness is beyond the scope of solvent screening during early design stages, and should be considered at a later stage of validation through laboratory experiments.
Figure 2 shows the effect of the operational parameters of Table 3 on the optimization results when the objective is to maximize the productivity. Since we perform several runs per design case, we obtain a range of values for the ethanol productivity and the extraction efficiency. In the figure, the 12 design cases are first sorted by decreasing DE and secondly by decreasing S 1 0 . We observe that EE is mainly affected by the DE value, and we can increase the final EE values by increasing the DE. Looking at Equations (14) and (15) we see that there is a certain dependence between EE and DE which could justify this trend. Regarding the productivity, we observe a steady reduction of its maximal value (with the exception of case A8) as we first reduce DE and secondly S 1 0 .
Detailed information on the obtained productivities and extraction efficiencies for each design case, as well as the generated solvents, can be found at Supplementary Table S1. We can see that under design case A4 ( S 1 0 = 150 g/L and DE = 4/h) we obtain the maximum productivities (~100 gr/L/hr) with designs D5, D7, and D9 of Table 4. The associated efficiencies are also very high (~0.97).
Note that, additional results obtained for maximum EE (Equation (15)) indicate that the EE cannot be improved further compared to results using the productivity objective. In addition, the resulting designs are more diverse and complex in their structures, and the possibility of converging to common design features is extremely low. These designs also include the –CH2CN group, many of them include two nitrile groups; we also get (inferior) solutions featuring aromatic bonds.

5. Discussion of Results

The obtained solvent designs exhibit different extraction efficiencies, productivities, and conversions according to the initial substrate concentration ( S 1 0 ) and the dilution rate (DE). The solvent ranking is maintained regardless of the particular choice of operating conditions, as indicated in Table S1 of the supplementary material. For the productivity objective, the optimal designs (namely D9, D5, and D7) and their ranking remain the same regardless of the design case. The conversion shows strong dependence on the choice of operation, at least in the region of maximal PDE. We also observe that the most productive solvents are also the most efficient ones. It is, therefore, safe to say that high productivity in the extract phase is accompanied by high EE (at least near the optimum), without any indication that the reverse holds. In general, the productivity objective (Equation (10)) is found to be a very good choice to address the solvent design goals and generate promising solvents.
The solvent molecules give a clear picture of the structures and group compositions that seem to be promising for each design case. The nitrile group “−CH2CN” is dominant in the designed solvent molecules, regardless of the choice of the objective function and operation. The designs of Table 4 also include ester groups (D1, D2, D8), chloride groups (D1, D4, D6, D10), and double c-bonds (D3, D4, D6). As an indication of the molecules size, their molecular weights range between 100 and 170 g/mole. These design features differ from conventional solvents (e.g., dodecanol) and their types (e.g., large alcohols).
Some of the designed solvents have similar structure and extraction performance. The structural similarities usually refer to chain isomers, and this explains the proximity in their physical property values. Note that, for more elaborate forms of isomerism (i.e., location, cis–trans, and chiral isomerism) the group-wise representation fails to distinguish between the different configurations. However, the physical behavior differences between such isomeric forms are expected to be small.
Designs D9, D7, and D5 have the highest productivities among the results of the productivity objective, and very high extraction efficiencies. They are chain isomers with the general molecular formula C4H9−CH2CN. This opportunity to generalize, demonstrates the consistency of the results and robustness of the CAMD optimizer. In terms of productivity, design D9 maintains a clear edge, with D7 and D5 following with about 0.8% lower productivity. For instance, according to Table S1, the productivities for D9, D7, and D5 are, for set A4: 100.8, 100.1, and 100.0 g/L/h, and for sets A7: 40.8, 40.5, 40.5 g/L/h, respectively. Their efficiencies are around 0.97 and 0.93 for data sets A4 and A7, respectively. Among the three designs, D5 has a simple straight chain structure, it is the less toxic (higher biocompatibility), and it has the best LCA indications with GWP, CED, and EI99 at 3.7, 110, and 0.43, respectively. Therefore, D5 appears as the promising candidate given the current problem objectives.
On the other hand, design D1 is the only one that satisfies the strict set of constraints for all three LCA indices. This design features its best PDE and EE under design case A4, evaluated at 63.4 and 0.919, respectively. D1 could also be considered for further investigation, depending on the importance of LCA performance over the design decisions.
Note that, in practice, the analysis of the results may indicate the need for alternative, modified or additional constraints/objectives to guide the screening process. So, we may need to reconsider, reset, and rerun the design problem several times before obtaining a result worth for further investigation. In addition, these choices depend on the application, and may include additional system information once it becomes available. The CAMD tool used here can easily generate a new series of solvent designs, if the constraints, the process data, and/or the process simulation model are modified.

6. Conclusions

This work applies computer-aided molecular design (CAMD) tools to design new efficient and sustainable solvents for extractive fermentation (EF) of glucose to bioethanol fuels. The EF process combines fermentation with a primary liquid–liquid extraction stage, and has the potential to reduce operation costs, increase conversion efficiency, and reduce the inhibition of products and by-products. Finding appropriate solvents to enhance the extraction efficiency without inhibiting the fermentation is an important step towards industrial application of EF. CAMD technologies can perform the solvent screening task in a way that is systematic, general, and robust.
The solvent design procedure presented here uses a state-of-the-art CAMD optimizer to get useful insights on promising molecular structures, and incentives to perform more detailed assessments before a new material is considered for industrial use. The solvent assessment involves a rigorous simulation model for EF, as well as group-contribution methods for the prediction of activity coefficients and various primary properties, toxicity, and LCA indices.
The solvent design studies presented here use 57 different functional groups as building blocks for a great variety of smaller or larger molecules. Despite the opportunities for complexity, the proposed solvent designs are not particularly complex, and indicate a strong preference to medium-sized molecules and particular functional groups. They are also very efficient in meeting the EF process requirements and maximizing the production of ethanol.
The current set up of the problem seems to favor hydrocarbons with chloride and nitrile groups. In principle, the presence of these groups is a strong indication for environmental hazard. Note that, in the current model for solvent selection the considered lifecycle impacts are limited to carbon emissions. Consideration of additional or different environmental impact categories [41] will generate different results [21]. The problem set up is therefore very important to capture the precise design objectives, and the CAMD optimizer used here is not affected by the number or complexity of these models. This is a very important advantage of the approach proposed here, for fast, efficient, and extended screening of solvent alternatives before going to more detailed investigations, using laboratory experiments and detailed process optimization.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/separations9030060/s1, Table S1: Productivities (in gr/L/hr) and extraction efficiencies obtained for the design cases of Table 2. The solvents designs and their properties are given in Table 4 of the main article.

Author Contributions

Conceptualization, E.M.; data curation, P.B.; formal analysis, E.M. and P.B.; methodology, E.M.; software, E.M.; supervision, E.M.; validation, P.B.; writing—original draft, E.M.; writing—review & editing, E.M. and P.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A, B, Ccoefficients in Equation (2)
CSsubstrate conversion []
D, Ddilution rates based on influent and effluent infinite aqueous phase flowrate [h−1], respectively
DE, DEdilution rates based on influent and effluent infinite extract phase flowrate [h−1], respectively
EEextraction efficiency [g/L product in extract/g/L product in raffinate]
F1objective 1: maximum ethanol productivity in the extract phase (Equation (12))
F2objective 2: maximum extraction efficiency (Equation (14)).
KPproduct inhibition constant [g/L]
KSsubstrate saturation constant [g/L]
LC50concentration causing 50% mortality [mole/L]
MPdistribution coefficient of solute (produced ethanol) [g in raffinate/g in solvent]
P1, P2effluent product concentration in raffinate and extract phases [g/L], respectively
PDEproductivity in the extract phase [g/L/h]
PDTtotal product (ethanol) productivity [g/L/h]
PDWproductivity in the nutrient phase [g/L/h]
S 1 0 , S1influent and effluent substrate concentration in raffinate phase [g/L], respectively
SEsolvent selectivity towards solute [g solute/g raffinate]
SLsolvent loss to raffinate [g solvent/g raffinate]
titoxicity of compound i, defined as ti = −log(LC50i) [] (LC50 is in [mol/L])
Tb,iboiling point temperature of compound i [K]
Tf,ifreezing point temperature of compound i [K]
TFprocess temperature [K]
YC/Scarbon dioxide yield coefficient [g CO2 produced/g substrate consumed]
YP/Sethanol yield coefficient [g ethanol produced/g substrate consumed]
YX/Scells yield coefficient [g cells produced/g substrate consumed]
WS, WYauxiliary variables for Equations (3)–(5)
Greek letters
α, βcoefficients in Equations (3)–(5)
ηratio of actual MP over its value at infinite dilution conditions
μspecific growth rate [h−1]
ρidensity of pure component i [g/L]
Superscripts and subscripts
infinite dilution conditions
icompound i
Eextract, solvent phase
Psolute, product (ethanol)
Wraffinate, aqueous phase
min, maxminimum and maximum values, respectively

References

  1. [EUGD] European Green Deal: Commission Proposes Transformation of EU Economy and Society to Meet Climate Ambitions, Press Release 14 July 2021, Brussels. Available online: https://ec.europa.eu/commission/presscorner/detail/en/IP_21_3541 (accessed on 14 February 2022).
  2. Regulation (EU) 2018/1999 of the European Parliament and of the Council of 11 December 2018 on the Governance of the Energy Union and Climate Action, Amending Regulations (EC) No 663/2009 and (EC) No 715/2009 of the European Parliament and of the Council, Directives 94/22/EC, 98/70/EC, 2009/31/EC, 2009/73/EC, 2010/31/EU, 2012/27/EU and 2013/30/EU of the European Parliament and of the Council, Council Directives 2009/119/EC and (EU) 2015/652 and repealing Regulation (EU) No 525/2013 of the European Parliament and of the Council (Text with EEA Relevance). Available online: http://data.europa.eu/eli/reg/2018/1999/2021-07-29 (accessed on 14 February 2022).
  3. Fulton, L.M.; Lynd, L.R.; Körner, A.; Greene, N.; Tonachel, L.R. The need for biofuels as part of a low carbon energy future. Biofuels Bioprod. Biorefining 2015, 9, 476–483. [Google Scholar] [CrossRef]
  4. Correa, D.F.; Beyer, H.L.; Fargione, J.E.; Hill, J.D.; Possingham, H.P.; Thomas-Hall, S.R.; Schenk, P.M. Towards the implementation of sustainable biofuel production systems. Renew. Sustain. Energy Rev. 2019, 107, 250–263. [Google Scholar] [CrossRef]
  5. Vohra, M.; Manwar, J.; Manmode, R.; Padgilwar, S.; Patil, S. Bioethanol production: Feedstock and current technologies. J. Environ. Chem. Eng. 2014, 2, 573–584. [Google Scholar] [CrossRef]
  6. Ayas, N. Energy Materials. In Comprehensive Energy Systems; Dincer, I., Ed.; (Section 2.13.3.2.2 Advantages of ethanol fuel); Elsevier: Amsterdam, The Netherlands, 2018. [Google Scholar]
  7. Zabed, H.; Sahu, J.N.; Suely, A.; Boyce, A.N.; Faruq, G. Bioethanol production from renewable sources: Current perspectives and technological progress. Renew. Sustain. Energy Rev. 2017, 71, 475–501. [Google Scholar] [CrossRef]
  8. Soleymani Angili, T.; Grzesik, K.; Rödl, A.; Kaltschmitt, M. Life Cycle Assessment of Bioethanol Production: A Review of Feedstock, Technology and Methodology. Energies 2021, 14, 2939. [Google Scholar] [CrossRef]
  9. Lemos, D.A.; Sonego, J.L.; Boschiero, M.V.; Araujo, E.C.; Cruz, A.J.; Badino, A.C. Selection and application of nontoxic solvents in extractive ethanol fermentation. Biochem. Eng. J. 2017, 127, 128–135. [Google Scholar] [CrossRef]
  10. Farias, D.; Maugeri-Filho, F. Sequential fed batch extractive fermentation for enhanced bioethanol production using recycled Spathaspora passalidarum and mixed sugar composition. Fuel 2021, 288, 119673. [Google Scholar] [CrossRef]
  11. Lemos, D.A.; Sonego, J.L.; Cruz, A.J.; Badino, A.C. Improvement of ethanol production by extractive fed-batch fermentation in a drop column bioreactor. Bioprocess. Biosyst. Eng. 2020, 43, 2295–2303. [Google Scholar] [CrossRef]
  12. Borhani, T.N.; Wang, M. Role of solvents in CO2 capture processes: The review of selection and design methods. Renew. Sustain. Energy Rev. 2019, 114, 109299. [Google Scholar] [CrossRef]
  13. Papadopoulos, A.I.; Tsivintzelis, I.; Seferlis, P.; Linke, P. Computer-Aided Molecular Design: Fundamentals, Methods, and Applications. In Reference Module in Chemistry, Molecular Sciences and Chemical Engineering; Reedijk, J., Ed.; Elsevier: Waltham, MA, USA, 2018. [Google Scholar]
  14. Chemmangattuvalappil, N.G. Development of solvent design methodologies using computer-aided molecular design tools. Curr. Opin. Chem. Eng. 2020, 27, 51–59. [Google Scholar] [CrossRef]
  15. Gani, R. Group contribution-based property estimation methods: Advances and perspectives. Curr. Opin. Chem. Eng. 2019, 23, 184–196. [Google Scholar] [CrossRef]
  16. Zhou, T.; McBride, K.; Linke, S.; Song, Z.; Sundmacher, K. Computer-aided solvent selection and design for efficient chemical processes. Curr. Opin. Chem. Eng. 2020, 27, 35–44. [Google Scholar] [CrossRef]
  17. Chai, S.; Song, Z.; Zhou, T.; Zhang, L.; Qi, Z. Computer-aided molecular design of solvents for chemical separation processes. Curr. Opin. Chem. Eng. 2022, 35, 100732. [Google Scholar] [CrossRef]
  18. Khor, S.Y.; Liam, K.Y.; Loh, W.X.; Tan, C.Y.; Ng, L.Y.; Hassim, M.H.; Ng, D.K.; Chemmangattuvalappil, N.G. Computer aided molecular design for alternative sustainable solvent to extract oil from palm pressed fibre. Process. Saf. Environ. Prot. 2017, 106, 211–223. [Google Scholar] [CrossRef]
  19. Papadopoulos, A.I.; Shavalieva, G.; Papadokonstantakis, S.; Seferlis, P.; Perdomo, F.A.; Galindo, A.; Jackson, G.; Adjiman, C.S. An approach for simultaneous computer-aided molecular design with holistic sustainability assessment: Application to phase-change CO2 capture solvents. Comput. Chem. Eng. 2020, 135, 106769. [Google Scholar] [CrossRef]
  20. Ten, J.Y.; Liew, Z.H.; Oh, X.Y.; Hassim, M.H.; Chemmangattuvalappil, N. Computer-Aided Molecular Design of Optimal Sustainable Solvent for Liquid-Liquid Extraction. Process. Integr. Optim. Sustain. 2021, 5, 269–284. [Google Scholar] [CrossRef]
  21. Baxevanidis, P.; Papadokonstantakis, S.; Kokossis, A.; Marcoulaki, E. Group Contribution-based LCA models to enable screening for environmentally benign novel chemicals in CAMD applications. AIChE J. 2021, e17544. [Google Scholar] [CrossRef]
  22. Zhou, T.; Song, Z.; Zhang, X.; Gani, R.; Sundmacher, K. Optimal solvent design for extractive distillation processes: A multiobjective optimization-based hierarchical framework. Ind. Eng. Chem. Res. 2019, 58, 5777–5786. [Google Scholar] [CrossRef]
  23. Marcoulaki, E.C.; Kokossis, A.C. On the development of novel chemicals using a systematic synthesis approach. I–Optimisation framework. Chem. Eng. Sci. 2000, 55, 2529. [Google Scholar] [CrossRef]
  24. Marcoulaki, E.C.; Kokossis, A.C. On the development of novel chemicals using a systematic optimisation approach. II–Solvent design. Chem. Eng. Sci. 2000, 55, 2547. [Google Scholar] [CrossRef]
  25. Papadopoulos, A.I.; Linke, P. Integrated solvent and process selection for separation and reactive separation systems. Chem. Eng. Process. Process. Intensif. 2009, 48, 1047–1060. [Google Scholar] [CrossRef]
  26. Papadopoulos, A.I.; Stijepovic, M.; Linke, P.; Seferlis, P.; Voutetakis, S. Toward optimum working fluid mixtures for organic Rankine cycles using molecular design and sensitivity analysis. Ind. Eng. Chem. Res. 2013, 52, 12116–12133. [Google Scholar] [CrossRef]
  27. Papadopoulos, A.I.; Zarogiannis, T.; Seferlis, P. Computer-Aided Molecular Design of CO2 Capture Solvents and Mixtures. In Process Systems and Materials for CO2 Capture: Modelling, Design, Control and Integration; Papadopoulos, A.I., Seferlis, P., Eds.; John Wiley & Sons Ltd.: West Sussex, UK, 2017; p. 10173. [Google Scholar]
  28. Kirkpatrick, S.; Gelatt, C.D.; Vecchi, M.P. Optimization by simulated annealing. Science 1983, 220, 671–680. [Google Scholar] [CrossRef] [PubMed]
  29. Marcoulaki, E.C.; Kokossis, A.C. Scoping and screening complex reaction networks using stochastic optimization. AIChE J. 1999, 45, 1977–1991. [Google Scholar] [CrossRef]
  30. Gutieérrez, L.F.; Sánchez, O.J.; Cardona, C.A. Analysis and design of extractive fermentation processes using a novel short-cut method. Ind. Eng. Chem. Res. 2013, 52, 12915–12926. [Google Scholar] [CrossRef]
  31. Kollerup, F.; Daugulis, A.J. A mathematical model for ethanol production by extractive fermentation in a continuous stirred tank fermentor. Biotechnol. Bioeng. 1985, 27, 1335–1346. [Google Scholar] [CrossRef]
  32. Cockrem, M.; Flatt, J.; Lightfoot, E. Solvent selection for extraction from dilute solution. Sep. Sci. Technol. 1989, 24, 769. [Google Scholar] [CrossRef]
  33. Pretel, E.J.; Lopez, P.A.; Bottini, S.B.; Brignole, E.A. Computer-aided molecular design of solvents for separation processes. AIChE J. 1994, 40, 1349. [Google Scholar] [CrossRef]
  34. Fredenslund, A. Vapor-Liquid Equilibria Using UNIFAC: A Group-Contribution Method; Elsevier: Amsterdam, The Netherlands, 2012. [Google Scholar]
  35. Reid, R.C.; Prausnitz, J.M.; Poling, B.E. The Properties of Gases and Liquids, 4th ed.; McGraw-Hill: New York, NY, USA, 1988. [Google Scholar]
  36. Constantinou, L.; Gani, R. New group contribution method for estimating properties of pure compounds. AIChE J. 1994, 40, 1697. [Google Scholar] [CrossRef]
  37. Wang, W.; Achenie, L.E.K. Computer aided solvent design for extractive fermentation. Fluid Phase Equilibria 2002, 201, 1–18. [Google Scholar] [CrossRef]
  38. Trevizo, C.; Nirmalakhandan, N. Prediction of microbial toxicity of industrial organic chemicals. Water Sci. Technol. 1999, 39, 63. [Google Scholar] [CrossRef]
  39. Gao, C.; Govind, R.; Tabak, H.H. Application of the group contribution method for predicting the toxicity of organic chemicals. Environ. Technol. Chem. 1992, 11, 631. [Google Scholar] [CrossRef]
  40. [EC] Communication from the Commission to the European Parliament, the Council, the European Economic and Social Committee and the Committee of the Regions Chemicals Strategy for Sustainability towards a Toxic-Free Environment. COM(2020) 667 Final, Brussels. 14 October 2020. Available online: https://ec.europa.eu/environment/pdf/chemicals/2020/10/Strategy.pdf (accessed on 14 February 2022).
  41. Park, W.J.; Kim, R.; Roh, S.; Ban, H. Analysis of major environmental impact categories of road construction materials. Sustainability 2020, 12, 6951. [Google Scholar] [CrossRef]
Figure 1. Outline of the CAMD optimization algorithm used here.
Figure 1. Outline of the CAMD optimization algorithm used here.
Separations 09 00060 g001
Figure 2. Ranges of productivity and EE in the obtained results.
Figure 2. Ranges of productivity and EE in the obtained results.
Separations 09 00060 g002
Table 1. Process parameters for the solvent design case study.
Table 1. Process parameters for the solvent design case study.
ParameterValue
α0.411
β−0.163
ρw1000 g/L
YX/S0.10
YP/S0.46
YC/S0.44
μmax0.45 h−1
Kp23.0 g/L
Ks1.0 g/L
Tb,P351.4 K
TF298 K
S130 g/L
η80%
Table 2. Bounds on the designed solvent properties.
Table 2. Bounds on the designed solvent properties.
PropertyFeasible RangeBoundValueUnits
S E SE,minSE,min2.0-
M P MP,minMP,min1.0-
S L SlmaxSlmax0.005-
M S MS,maxMS,max0.1-
Tb,E∈[Tb,P + Tmin, 600] or [300, Tb,P − Tmin]ΔΤmin30. K
Tf,E≤TF − ΔΤminΔΤmin10. K
ρE≥ρ1⋅Δρ,maxΔρ,max125%-
≤ρ1⋅Δρ,minΔρ,min80%-
tE≤tE,maxtE,max1.5, 2.0, 2.5, 3.0, 3.5, 10.mol/L
GWP≤GWPE,maxGWPE,max8.0, 10.0g CO2-eq/g
CED≤CEDE,maxCEDE,max150, 200kJ-eq/g
EI99≤EI99E,maxEI99E,max0.4, 0.5-
Table 3. Operational parameters per solvent design case.
Table 3. Operational parameters per solvent design case.
ParameterParameter Value per Design Case A1 to A12
S 1 0 [g/L]A1, A2, A3, A4A5, A6, A7, A8A9, A10, A11, A12
150300600
DE [h−1]A1, A5, A9A2, A6, A10A3, A7, A11A4, A8, A12
0.51.02.04.0
Table 4. Solvent designs obtained for maximum ethanol productivity.
Table 4. Solvent designs obtained for maximum ethanol productivity.
NoSMILES *tEMESl × 103Tb,ErETf,EGWPCEDEI99
D1ClCCC(=O)OCC#N1.690.67110.94989562463.2386.20.389
D2CCC(=O)OCCCCC#N2.081.314.875107562453.7699.40.433
D3C=CCCCCC#N2.431.632.154626542214.981410.436
D4C=C(CCl)CC#N2.441.315.844658462293.361150.444
D5CCCCCC#N2.482.466.554426462053.681100.433
D6C=CC(Cl)CC#N2.751.275.164618162235.77151 0.442
D7C(C)(C)CCC#N2.912.456.534346541904.141320.503
D8C(C)(C)(C)C(=O)OCCC#N2.921.242.285137672525.16153 0.562
D9C(C)(C)(C)CC#N3.032.527.164256772034.94162 0.558
D10C(C)(C)C(C=O)CCl4.511.598.554458612286.71171 0.483
Notes: * The results of the CAMD tool may relate to more than one molecular structure, and each SMILES presents one of these alternatives; value between the two LCA bounds; value outside the LCA bounds of Table 2.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marcoulaki, E.; Baxevanidis, P. Screening for New Efficient and Sustainable-by-Design Solvents to Assist the Extractive Fermentation of Glucose to Bioethanol Fuels. Separations 2022, 9, 60. https://doi.org/10.3390/separations9030060

AMA Style

Marcoulaki E, Baxevanidis P. Screening for New Efficient and Sustainable-by-Design Solvents to Assist the Extractive Fermentation of Glucose to Bioethanol Fuels. Separations. 2022; 9(3):60. https://doi.org/10.3390/separations9030060

Chicago/Turabian Style

Marcoulaki, Effie, and Pantelis Baxevanidis. 2022. "Screening for New Efficient and Sustainable-by-Design Solvents to Assist the Extractive Fermentation of Glucose to Bioethanol Fuels" Separations 9, no. 3: 60. https://doi.org/10.3390/separations9030060

APA Style

Marcoulaki, E., & Baxevanidis, P. (2022). Screening for New Efficient and Sustainable-by-Design Solvents to Assist the Extractive Fermentation of Glucose to Bioethanol Fuels. Separations, 9(3), 60. https://doi.org/10.3390/separations9030060

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