Energy-efficient bioethanol recovery process using deep eutectic solvent as entrainer

Bioethanol is commonly recovered from fermentation broth via distillation because it is the most economical and reliable process for large-scale industrial operations. Because extraction is a highly energy-efficient process applicable to low-composition bio-product separation, high-performance solvents are necessary for efficient bioethanol recovery. 2-Methyl pentanol, a branched long-chain alcohol, which was used as the extraction solvent, and a deep eutectic solvent, choline chloride and ethylene glycol (1:2), was employed as an entrainer for product refinement. Thermodynamic models of vapor–liquid equilibrium (VLE) and liquid–liquid equilibrium (LLE) systems were derived using molecular simulations and experimental results to develop the bioethanol extraction and refinement processes. The heat duty of the designed process was reduced by a quarter compared to that of previous recovery processes.


Introduction
Bioethanol is a sustainable fuel used extensively in automobiles.Bioethanol recovery is generally achieved via distillation because it is the most economical and reliable solution for large-scale industrial operations.Distillation is suitable and has long been used for ethanol recovery because ethanol has a significantly lower boiling point than water, unlike other bio-alcohols [1].However, as the distillation process consumes a large amount of energy, improvements are necessary to enhance sustainability [2][3][4][5][6].The recovery of ethanol from bioethanol involves removing approximately 90% of the water from the feedstock.Various techniques can be implemented for ethanol-only recovery from low-ethanolcontent water mixtures, including membrane distillation [7,8], supercritical extraction [9,10], hybrid extraction/distillation [11][12][13], pressure swing adsorption [14][15][16], direct gasohol production [17], and extractive distillation [18][19][20].
Extraction is a favorable recovery process owing to its low energy demand.Although the solvent recycling step in this process requires energy, an efficient solvent with a low solvent-to-feedstock ratio and a sufficiently large boilingpoint difference with the product can reduce the energy demand.Extraction can be performed with high efficiency at room temperature without heating, and recycling a lowtemperature solvent provides heat recovered from the hightemperature recycled solvent drawn from the solvent recovery column.Extraction is an established technique that has undergone substantial design and operation advancements.Therefore, the practical application of a new solvent can lead to improvements when it is developed using a known liquid-liquid equilibrium (LLE) model [21].
The high affinity between ethanol and water that gives rise to its azeotrope makes it difficult to produce ethanol with a purity of over 95%.Either extractive or azeotropic distillation [22][23][24][25] can be applied to produce pure ethanol, but azeotropic distillation consumes more energy because of the high vapor flow in the rectifying section of the distillation column.Benzene and cyclohexane are used for azeotropic distillation, while ethylene glycol, glycerin, and deep eutectic solvents (DESs), such as choline/ethylene glycol [26,27], choline/urea [28][29][30], and choline/glycerol [31], are suitable for extractive distillation.A comparison of the energy demands of various ethanol recovery processes indicates that extractive distillation is preferable for ethanol concentration.
In this study, 2-methyl pentanol was utilized for ethanol extraction from a 10% ethanol fermentation broth, and ethanol concentration was achieved via extractive distillation using a DES, namely, choline/ethylene glycol (1:2).The energy efficiency of the developed system was compared to those of existing processes applied to the same ethanol feedstock to demonstrate superior performance of the former.A similar procedure was applied for the production of 2,3-butanediol [32], ethanol [13], and butanol [33] from fermented feedstock.

Solvent selection
The proposed bioethanol recovery process encompasses one extractor and three distillation columns.Most of the water is removed from the feedstock via the extraction step and first distillation, and the two subsequent extractive distillations purify the pre-processed ethanol.Two solvents, 2-methyl pentanol and a DES, were used during the recovery process.Generally, distillation is employed for the pre-processing of low-content ethanol, whereas in this study, extraction was applied for this purpose [34].The processed ethanol was purified in the next extractive distillation step using choline chloride/ethylene glycol at a molar ratio of 1:2.The suitability of choline chloride as an entrainer for the ethanol/ water azeotrope has been proposed and experimentally demonstrated in a previous study [35].
There are three major requirements an extraction solvent must satisfy: (1) high selectivity toward ethanol, (2) low solubility in water, and (3) a sufficient difference between its boiling point and that of ethanol [36].Overall, 2-methyl pentanol is proposed as an optimal solvent for ethanol extraction considering the toxicity of higher alcohols and the boiling point requirement.The wastewater from the extraction process contains a small amount of the extraction solvent, which is recycled in the ethanol fermentation process [37], as in most bio-processes; hence, its toxicity to microorganisms must be considered.

Thermodynamic models
The proposed ethanol recovery process was developed using Aspen Plus ver.8.8, which uses vapor-liquid equilibrium (VLE) and LLE models for thermodynamic computations.

NRTL model
Although the design of the extraction and distillation processes can be readily achieved using the commercial design software Aspen Plus [38], its database does not contain thermodynamic information for 2-methyl pentanol and DES.For VLE and LLE computations, the non-random two-liquid (NRTL) model [39] was used, and molecular simulations were conducted to obtain the equilibrium data because experimental measurements were not available.Because the extraction process is crucial in the proposed system, several LLE data were experimentally obtained to evaluate the validity of the simulated LLE results [13].
The parameters of the NRTL model were estimated from the experimental and molecular simulation results, and the activity coefficients were calculated using the following relationship: where N is the component number and x j is the liquid composition of the jth component.The binary interaction parameters (a ji and b ij ) and non-randomness factors (α ji ) were used as the parameters in the above model: The relation between vapor and liquid systems is given by: where P i is the pressure of pure component i at absolute temperature T, and P T is the system pressure.φ i is the fugacity coefficient and y i is the vapor composition.

Molecular simulations
The feedstock components were water and ethanol, and two solvents, 2-methyl pentanol and a DES, were utilized for extraction and extractive distillation, respectively.The binary interaction parameters of the VLE system among water, ethanol, and ethylene glycol are available in the Aspen Plus database, but those of 2-methyl pentanol and the DES are not.For the VLE system, DES-related binary parameters are given in [27].Therefore, 2-methyl pentanol-related VLE and LLE parameters were estimated using molecular simulation results [13].
The LLE binary parameters were estimated based on the molecular distribution in two liquid phases, computed using molecular simulations to determine the minimum potential among the constituent molecules in the respective liquid phases.The liquid phase was assigned as a cube in the molecular simulation.The force fields in molecular systems (1) consist of intramolecular and intermolecular potentials; the former includes bond stretching, angle bending, and torsion, and the latter includes van der Waals and Coulomb forces, and are thus calculated as follows [40]: where r is the spatial location of a molecule, r ij is the distance between atoms i and j in adjacent molecules, and k is the potential constant [41].The bend angle (θ) and torsion angle (φ) are with respect to the reference axis and plane, respectively [40].The constants ε, σ, and q in the last term are for specific atoms.Molecular simulations determine the molecular distribution in the two boxes, which in turn reveals the compositions of the two phases at equilibrium.There are numerous parameters in the simulation that significantly affect the simulation results.By using the experimental LLE data in a reference system consisting of various solvents, the tuning parameter was adjusted to its optimum value for the proposed LLE system.The molecular simulation was conducted using RASPA, an open software, the installation details of which are explained in the installation and simulation manual [42].RASPA performance is discussed in [43], and an example of the input file used in the simulation is listed in File S1 of the Supplementary Material.In the file, the TranslationProbability, RotationProbability, ReinsertionProbability, and GibbsSwapProbability parameters are adjustable; among them, GibbsSwapProbability was the most important for our application.Considering previous molecular simulation studies wherein LLE was applied in butanol and 2,3-butanediol extractions [32,33], a value of 0.05 was deemed optimal in the LLE simulation, while the other parameters were set to the values in the file.In this study, phase transition probability was selected as an adjustable tuning parameter, while all other parameters were kept constant during the simulation.The molecular simulation results were partially compared with the experimental results to examine their accuracy.Similarly, the molecular distributions in the vapor and liquid phases were obtained from the potential computation, and the VLE data were used to estimate the VLE binary interaction parameters [44][45][46].The liquid compositions at equilibrium determined using molecular simulation were compared with the experimental results obtained in [13].Although numerous studies [32,33,47] have implemented molecular simulations for establishing the thermodynamic and physical properties of various compounds, the reliability of the results is generally inadequate.In particular, process designs that necessitate large ( 5) investments toward plant construction require experimental verification of the computed results.

Process design
The proposed process consists of two stages: (1) water removal via extraction and (2) ethanol refinement via extractive distillation.Because no single solvent is suitable for both these processes, different solvents were used for each stage in this study.2-Methyl pentanol was used as the water removal solvent, considering its water solubility and boiling point.A DES consisting of choline chloride/ethylene glycol (1:2) was used for ethanol concentration, as previously reported [26,27].Figure 1 illustrates a simplified scheme of the proposed bioethanol recovery process.The process was designed using Aspen Plus ver.8.8.
The feed contained 10% ethanol, and the targets for ethanol purity and recovery were 99% and 98%, respectively.The tray number for the extractor was determined for bioethanol recovery.The distillation columns were optimized in terms of their minimum total annualized cost (TAC), which is a combination of investment and operating costs.The cost equations used are reported in the literature [48,49], and a cost index of 1638.2 (the Marshall and Swift equipment cost index in 2018) was used.The payback time for the investment was 10 years, and the number of annual operation days was 300.Optimization was conducted using the response surface method [50].

Results and discussion
The thermodynamic parameter estimation, process design, and performance evaluation results are presented below.

NRTL parameters
Figure 2 shows a ternary diagram for water/ethanol/2-methyl pentanol.The plus ( +) symbols and squares represent the molecular simulation and experimental results, respectively [13]; the circles represent the values computed using the NRTL model and the estimated parameters listed in Table 1; and the multiplication symbols are the Aspen Plus computed values.The connecting curves are the LLE compositions computed using the NRTL model.Agreement between the NRTL computation and Aspen Plus estimation ensures that Aspen Plus process design is supported by experimental and molecular simulation results.The LLE data are listed in Table S1 of the Supplementary Material, and the standard deviation was 0.0026.In the design of the LLE measurement experiment, the feed composition should provide a system of liquid-liquid separation.A preliminary ternary diagram constructed using the molecular simulation results was used to determine the feed compositions on the tie lines of the split liquid phases.The curve computed using the NRTL model agreed well with the molecular simulation and experimental results.The LLE results of the molecular simulation and NRTL prediction of the water/ethanol/2-methyl pentanol system are summarized in Table S2 of the Supplementary Material.The mean deviation between the molecular simulation results and the predicted NRTL values was 0.0127.
The VLE parameters of the system consisting of water/ ethanol/choline chloride/ethylene glycol listed in Table 1 were obtained from the Aspen Plus database and reported in [27].The 2-methyl pentanol parameters were estimated from the VLE distribution predicted by molecular simulations.As a demonstration of the VLE molecular simulation results, the Txy diagram for the ethanol/2-methyl pentanol system is illustrated in Fig. S1 of the Supplementary Material.The circles and plus ( +) symbols of the ethanol liquid and vapor compositions agree well with the curves computed using the NRTL parameters.The VLE values obtained from molecular simulations and NRTL calculations are summarized in Table S3 of the Supplementary Material, and the mean deviation in the ethanol vapor composition was 0.0121.

Design results
The design parameters of the column structure and operating conditions for the proposed process are listed in Table 2. Figure 3 illustrates the entire design process, showing the material and energy flows among the extractor and distillation columns, including their structural information.In Aspen Plus process design, the reflux flow was adjusted to establish the specified product specification with a fixed tray number and column operating pressure.The tray numbers were TAC optimized for achieving an ethanol recovery of 98.3% and purity of 99.3%.The composition profiles of the organic phases in the extractor are shown in Fig. S2 of the Supplementary Material.The extract composition indicates high ethanol content at the top of the extractor, as shown in the figure, and the extract was composed of mostly ethanol and solvent.The low ethanol content at the bottom of the extractor is related to the LLE composition of the aqueous phase, as shown in Fig. 2, and is indicative of minimal ethanol loss into the wastewater from the extractor.The ethanol recovery set by the process design objectives is related to the loss, and the number of extractor stages was determined as the design target.To improve ethanol recovery, the solvent flow rate must be increased, which in turn increases the equipment size and energy demand.In this study, a recovery of 98.3% was selected, avoiding significant solvent flow increments and loss of ethanol.Recently reviewed bioethanol dehydration processes reported purities of 99.2-99.9%[51].Although an initial solvent charge is necessary as holdup in the extractor, a small amount of solvent make-up is required, resulting in minimized process costs.The diameter of the designed extractor was 1.13 m, as determined using design guidelines [52], and the height was 7.3 m, based on two-foot tray spacing.
Bioethanol is fermented from grains, which necessitates the use of copious quantities of water.Wastewater can only be recycled in the fermentation process if the extraction solvent is non-toxic to fermentation microorganisms.The toxicity of 2-methyl pentanol is classed as Category 4 [53].A toxicity study on yeast found that low-molecular-weight alcohols inhibited yeast growth at high concentrations [54].Although these solvents inhibited yeast growth at saturation, they were partially or non-inhibitory at concentrations lower than 10% saturation [55].While the solubility of 2-methyl pentanol is 0.31%, its content in the wastewater was less than 0.01%; thus, the toxicity was not a concern for the recycling step.Recycling the wastewater minimizes product loss and maximizes recovery.The extraction process removes 97% of the water from the feed, leading to lower energy consumption in the subsequent distillation process.For energy recovery, preheaters were used for feedstock heating in the three distillation columns.Because the bottom temperature in the first solvent recovery column is high and the extractor is operated at near room temperature, 58% of its reboiler duty is recovered from its cooler.The recovered heat accounts for 55% of the preheater heat demand.The availability of recovered heat for reuse in the preheaters was evaluated using a heat exchanger network design with pinch technology [56], as shown in Fig. 4. With a pinch temperature of 10 °C, sufficient utility was available for feedstock preheating in the distillation columns.

DES entrainer
The second stage of the process entails distillation to increase the ethanol purity from 79 to 99.3% after the solvent removal step.The purification process comprises a number of steps, including azeotropic and extractive distillation using various solvents.Although a detailed comparison of the energy demand is given in the following section, the design result for extractive distillation using choline chloride/ethylene glycol as the entrainer is presented in Fig. 5.The extractive distillation column has two feed streams, including the entrainer stream.The entrainer was introduced near the top, and the feed from the first solvent recovery column was near the bottom to obtain the highest ethanol composition in the overhead product, the target product of this project.Most of the water in the feed from the previous column was confined to the bottom product, which consisted of water and the entrainer.The entrainer was recovered in solvent recovery column II, and its overhead product, a mixture of water and ethanol, was recycled to the extractor.
Figure 6 depicts the extractive distillation column profile, demonstrating the efficiency of the entrainer in overcoming the ethanol/water azeotrope.The plus symbols in the x-y Fig. 4 A diagram of the heat exchanger network for heat recovery.The numbers below the circles represent the heat transfer rate in kW, and those on the lines indicate the temperature in centigrade Fig. 5 Vapor composition profiles in the extractive distillation column with and without a deep eutectic solvent entrainer, demonstrating the differences in azeotrope treatment diagram of the ethanol/water mixture without the entrainer show the changes in azeotropic composition at high ethanol contents.The circles indicate higher ethanol vapor compositions compared to liquid compositions.Although ethylene glycol has been utilized previously as an entrainer to concentrate ethanol [18,57,58], the amount used was > fourfold that used in this study for equal ethanol production.Choline chloride/urea (1:2) is another DES that has been used for flash-tank ethanol concentration [30,59], and in the distillation column [29] for solvent recovery.

Performance evaluation
The performance of the proposed process was evaluated by comparing its energy demands with those of previously reported processes.Because distillation is the most economical and reliable bioethanol recovery approach in large-scale industrial operations, many bioethanol processes apply conventional distillation for pre-processing fermented feedstock, as shown in Table 3. Biowaste generated from biomass pre-treatment during bioethanol fermentation is abundant and can be utilized as an energy source [60].Extraction has been used as an energy-efficient process for pre-processing fermented beer [1,36,61].At comparable feedstock compositions, the energy demands for pre-processing via distillation and extraction do not differ significantly owing to low solvent extraction efficiencies.Excluding azeotropic distillation with benzene, which has an excessive energy demand, the proposed process requires 25% less energy compared to the average demand of the previously reported operations listed in Table 3.The extraction solvent used in this study is more efficient than those used in other extraction processes.The solvent selection guidelines for bioethanol recovery [36] list branched long-chain alcohols as first -choice; the judicious selection of the solvent used in this study was based on its boiling point, which is a crucial factor in the solvent recovery process.Distillation continues to be widely implemented for bioethanol recovery [18][19][20], despite the high energy demand.The extraction efficiencies of previously used solvents were relatively low, resulting in high energy consumption.

Economics
The investment and operating costs for the extractor and the three distillation columns are listed in Table 4.While the investment includes column and tray costs and the costs of constructing the heat exchangers of the reboiler, condenser, preheater, and cooler, the operating cost represents the steam and water cooling costs.The steam cost was only calculated for reboiler duties because the preheaters use recovered heat.The steam cost was $24/t in 2017, as estimated by scaling the 2001 price given in [48].With a cash flow of 59% of utility cost and a return of investment of 23% at 7% interest rate, the payback time was 4.7 year at a discount rate of 21.4% [62].Considering operation for 300 days per year, the bioethanol production cost via recovery is US$0.10/Lethanol.This recovery cost is comparable to previously reported total production costs of US$0.66/L ethanol [63], US$0.84/L ethanol [64], and US$0.77/L ethanol [65].The downstream recovery cost was in the range of 30-40% of the total biofuel production costs [66].The recovery cost was compared to those of several reported processes, as listed in Table 5.The cost of the proposed process is 10% less than the average cost of five previously reported processes listed in the table, excluding the exceptionally high-cost example.The National Renewable Energy Laboratory (NREL) of the US Department of Energy examined bioethanol production from milled corn stover [67].The study entailed a practical technoeconomic evaluation that included cost estimation, and a process was developed with the guidance of Harris Group Inc., consultants with hands-on experience in commercial bioethanol plant engineering.The separation and purification costs reported in the above study were compared with those estimated herein.The reported 22.9 ton/h bioethanol production project entailed the use of two distillation columns and a molecular sieve adsorption tower for dehydration.The capital cost of the project was $26.98 m, equating to $2.73 m when a scaling-up exponent of 0.68 is applied [68].Meanwhile, the investment cost of the proposed process is $1.65 m, which is a reasonable improvement considering the higher solvent extraction efficiency.

Conclusions
Distillation is commonly applied for bioethanol recovery from fermentation broth as it is the most economical and reliable approach for large-scale industrial production.Preprocessing of low-content bioethanol was achieved via solvent extraction instead of conventional distillation to reduce the energy demand, and the product was concentrated via extractive distillation using a deep eutectic solvent to overcome the azeotrope of the ethanol/water mixture.Thus, a bioethanol recovery process entailing the combination of extraction pretreatment and extractive distillation using choline chloride/ethylene glycol (1:2) was designed and optimized for energy efficiency.The energy demand of the developed process was 25% less than the average demand of several previously reported processes.Furthermore, the proposed process is more cost efficient than the currently used processes.

Fig. 1 Fig. 2
Fig. 1 A simplified schematic diagram of the proposed bioethanol recovery process

Fig. 3
Fig. 3 Process flow diagram of the proposed bioethanol recovery process

Fig. 6
Fig. 6 Column profile for the extractive distillation column

Table 1
Binary interaction parameters in the LLE and VLE NRTL models

Table 4
Cost evaluation of the proposed recovery process.The units are in million US dollars and the utility cost is per annum coefficient (-) or torsion angle (°); γ: Activity coefficient (-); θ : Band angle (°S ubscripts i: Component i; j: Component j