Next Article in Journal
Enhanced Biomechanical Properties of Polyvinyl Alcohol-Based Hybrid Scaffolds for Cartilage Tissue Engineering
Previous Article in Journal
Impact of MiRNA-181a2 on the Clinical Course of IDH1 Wild Type Glioblastoma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Zero-Discharge Process for Recycling of Tetrahydrofuran–Water Mixtures

by
Karina Schuldt
,
Torsten Brinkmann
* and
Prokopios Georgopanos
*
Helmholtz-Zentrum Hereon, Institute of Membrane Research, Max-Planck-Str. 1, 21502 Geesthacht, Germany
*
Authors to whom correspondence should be addressed.
Processes 2021, 9(5), 729; https://doi.org/10.3390/pr9050729
Submission received: 25 March 2021 / Revised: 12 April 2021 / Accepted: 19 April 2021 / Published: 21 April 2021
(This article belongs to the Section Chemical Processes and Systems)

Abstract

:
The sustainable design of separation and polymer synthesis processes is of great importance. Therefore, an energy-efficient process for the purification of tetrahydrofuran (THF)–water (H2O) solvent mixtures from an upstream polymer synthesis process in pilot scale was developed with the aim to obtain high purity separation products. The advantages and limitations of a hybrid process in the pilot scale were studied utilizing an Aspen Plus Dynamics® simulation at different pressures to prove the feasibility and energy efficiency. For the rough separation of the two components, distillation was chosen as the first process step. In this way, a separation of a water stream of sufficient quality for further precipitations after polymer synthesis could be achieved. In order to overcome the limitations of the distillation process posed by the azeotropic point of the mixture, a vapor permeation is used, which takes advantage of the heat of evaporation already used in the distillation column. For the purpose of achieving the required low water contents, an adsorption column is installed downstream for final THF purification. This leads to a novel hybrid separation process that is energy efficient and thus allows also the use of the solvents again for upstream polymer synthesis achieving the high purity requirements in a closed-loop process.

1. Introduction

Nowadays, it is increasingly important to work more sustainably in chemistry. One area to be addressed is purification, e.g., for enabling the reuse of organic solvents. The increased usage of such applications is expressed by the worldwide development of the market value of solvent recycling equipment: while in 2016 the market value was approximately USD 900 million, a global value of more than USD 1 billion is expected in 2022 [1]. In the chemical industry are numerous combinations of organic solvent and water, which are waste products in various processes and, in many cases, form an azeotropic mixture. For the separation of those mixtures, several techniques are used, including membranes, pressure swing distillation, and extractive distillation, because conventional distillation is not sufficient for azeotropic mixtures [2,3,4,5,6]. However, the achievement of high purity at low cost and energy consumption is still a challenge. Therefore, intelligently designed hybrid processes combining the advantages of different separation techniques are attractive, especially when the focus is put on low energy consumption. Since separation processes are among the most energy-consuming processes in the chemical industry, the energy-saving design of such processes has a high priority.
Membrane separation processes are known for their energy efficiency and are straightforwardly integrated with distillation. Moreover, they convince with several advantages in their application, such as high selectivity, the moderate cost-to-performance ratio, avoidance of entrainers, steady-state operation, and compact and modular design [7]. However, the membrane has to be selected carefully with respect to the desired separation and the stability of the membrane materials at the anticipated operating pressure and temperature [8]. The use of membrane processes for the dehydration of process streams has already been investigated with different membrane materials, but no high-quality solvent purification, as aimed for in this study, has been achieved [9,10]. Due to the fact that the operating principle depends on the difference of fugacities across the membrane that only can be maintained at extensive costs for ultrapure high-pressure side products, a membrane separation process cannot provide an ultrapure separation of a mixture. Classical distillation processes, e.g., pressure swing distillation or azeotropic or extractive distillation, also suffer from distinct drawbacks, e.g., the energy requirement of operating at different pressure levels or the addition of an entrainer, resulting in narrow operating parameter windows and not allowing for flexibility in operation. The usage of a hybrid process consisting of distillation and a membrane stage can overcome the limitations of the individual processes. The membrane processes well suited for integration into processes aiming at the separation of close boiling or azeotropic mixtures are vapor permeation or pervaporation [11]. In general, the combination of a membrane process with a classical process is more economical than the use of a membrane process alone [12]. Hybrid processes combining distillation and vapor permeation are significantly more energy efficient, particularly in terms of heating and cooling power consumption, than conventional distillation processes [13,14]. Studies by Roth et al. on the dehydration of ethanol, for example, show an energy saving of up to 29% using a combined process with vapor permeation and adsorption compared to configurations with distillation [14]. Since dry organic solvents are important for many industrial sectors, the optional integration of an adsorption process as a final polishing stage can ensure ultrapure products [15]. Therefore, the intended process should be designed for high purity products and low energy consumption. It is known that optimization with regard to additional parameters is also successfully implemented with sensitivity analysis and multiobjective optimization [16]; however, this kind of optimization was not part of this work because the focus lies on the closed-loop character and the low energy consumption. The recycling of solvents and nonsolvents that are used in polymer production is of great importance since this implies high processing and high waste management costs [17]. The combination of modern methods for the polymer synthesis such as continuous flow chemistry [18] model-assisted polymer synthesis [19], as well as the use of the emulsion or dispersion polymerization techniques [20,21] in which water replaces some of the used organic solvents, do not solve the problem since, again, the large volume of solvent mixtures from the synthesis function through the material processing (i.e., precipitation of polymer after the synthesis) lead to the large volume of solvents that need to be processed.
In this publication, the purification of solvent waste from polymerization experiments, exemplary tetrahydrofuran (THF)–water (H2O), is investigated because THF, in particular, is widely used in polymer synthesis, for example, in the synthesis of block copolymers [22]. The recycling of THF from waste streams is generally of great interest because of the versatile and multiple uses of this solvent due to its ability to dissolve both protic and aprotic compounds, which leads to a large field of application [23]. Furthermore, THF has a particularly high energy consumption in both production and recovery by distillation, compared to other organic solvents, and is therefore not recommended from an environmental perspective [24]. The development of a new sustainable recycling process that makes THF available for reuse in high quality would contribute to the sustainable development of polymer synthesis. Based on the process of precipitating the synthesized polymer in water after a successful polymerization, a THF–H2O mixture with 20 mass percent THF is used as the feed for the process developed. The synthesis of polymers in pilot scale is considered, which is defined here as batch experiments in reactors with volumes of up to 20 L. The subsequent precipitation of the polymerization solution with H2O results in a mixture of organic solvent with water with a total volume of about 100 L. The components of the process are designed in a manner that the purification is possible within one day. Therefore, a pilot scale plant is developed.
Since the main goal is the reuse of the solvent THF, it is defined that THF is necessary to be obtained from the separation process in high purity. To create a closed-loop process, it is examined how the water of the mixture, previously considered as waste, can also be treated in order to allow for the reuse for precipitation of the polymer synthesized next. In this way, the potential of sustainable design approaching a zero-discharge process is to be fully exploited. Figure 1 shows how the solvents used in a polymerization process can be reused in the next polymerization process. The organic solvent–water mixture produced after polymer synthesis and subsequent precipitation of polymer product is to be processed in the hybrid process so that the organic solvent can be used again in a new polymer synthesis experiment and the water can be used again for the precipitation of the polymer. In this hybrid process, distillation, membrane separation and adsorption are combined in such a way that purification of the solvent and nonsolvent is achieved with low energy input. The individual advantages of the different processes are combined by using distillation for coarse separation, the membrane process to overcome the azeotropic point and adsorption for ultrafine purification. Due to the complexity of the targeted hybrid process (safety aspects, expenses, and demands for large rooms to host this pilot plant), the simulation-based development of the process before the construction of a plant is of great importance. Additionally, the models for the simulation can be used for the advancement of the process and for the monitoring of the process as well as the examination of the possibilities of using this plant for additional separation tasks. Therefore, this work proposes a novel closed loop process, supported by the simulation that can be further applied for the utilization of the organic waste mixtures from the upstream polymer synthesis, which differs from other processes in particular due to the high purity of the processed products.

2. Process Description

The separation of a 20 mass percent THF–H2O mixture as an example for organic solvent recovery is theoretically investigated by using process simulation to create a design for a closed-loop process. The challenge is that THF and water form an azeotropic mixture. Therefore, special distillation processes or the use of different separation principles are required, especially to meet the demand of high purity. A single distillation column, for example, is not sufficient for azeotropic mixtures but can remove the largest share of water. Consequently, distillation is chosen as the first process step. In this first step, the feed stream is separated into two streams. The aim of the column design is that the bottom product provides water in a quality that makes it reusable. Thus, the bottom product of the distillation is water containing very small amounts of solvent. This by-product can be reused in the preceding polymer precipitation process from which the waste stream to be separated originates. Depending on the azeotropic point and the pressure in the distillation column, in the second product stream, the top product still contains a large amount of water. To overcome the limitations of the vapor–liquid equilibrium-based separation process, a membrane process is used for the further treatment of the top product, whose separation principle is a selective solution–diffusion mechanism. The main objective of this process step is the further separation of water from the top product of the distillation process. Since the first process step delivers a vapor stream, vapor permeation, instead of pervaporation, is chosen to avoid an additional phase separation and to reduce the energy demands. Apart from this advantage, vapor permeation is more efficient in comparison to pervaporation if the aim is to overcome the azeotropic composition [25]. Additionally, due to in general less swelling of the membrane material, the durability is longer, and there is a lower sensitivity to concentration polarization, compared to the pervaporation process [26]. In view of the fact that the efficiency of membrane processes depends on the driving force generated, increased pressure is used on the feed side of the membrane. In order to make the overall process as efficient as possible, this pressure is generated upstream of the distillation column. Even though this leads to a shift of the azeotropic point to a higher water content in the top product, this is not to the detriment of the overall process because the membrane process is a less energy-intensive process, compared to the classical distillation process and is particularly efficient for higher water concentrations in the feed stream.
For the vapor permeation, a membrane module containing envelopes of thin-film composite membranes was assumed. Specific information on the membrane type and membrane area is obtained and described in the vapor permeation simulation part. However, the membrane module is designed in such a way that the retentate stream is as pure as possible. In order to achieve this, a permeate stream of lower purity containing considerable amounts of THF besides the main component water was acceptable. The purity of this stream was of no large concern since it should be fed back to the inlet of the distillation column. Hence, the THF in the permeate is not lost but is returned to the purification process. For ultrafine purification of the main THF stream, an adsorption process is used downstream of the membrane module. For this purpose, liquid-phase adsorption was selected. The operation of an adsorption column in liquid instead of gas phase has two main advantages: there is no risk of undesirable condensation and also the process more closely approaches isothermal operation since the release of the contribution of the heat of condensation to the heat of adsorption is avoided. In combination with this part of the process, the only waste is produced when the loaded adsorption column is regenerated with nitrogen. Moreover, nitrogen is used for providing inert conditions inside the plant before the process starts in order to ensure safe operation. In sum, the design results in a three-stage hybrid process integrating the unit operations distillation, vapor permeation, and adsorption. The focus in the design of this novel hybrid process lies in the selection and combination of the various individual separation processes in which efficiency is based on their coordination with each other, enabling an energy-efficient process. The process is intended to distinguish itself from previous purification ways developed for the separation of organic solvent and water, in particular by producing high-purity product streams, both of water and THF, so that all components can be reused without losses. Figure 2 shows a flowsheet of the hybrid process for the separation of the described azeotropic mixture.

3. Methods

A mixture of THF and water is investigated as the solvent mixture to be purified. Since a typical mixture of a waste product from the precipitation step of polymerization on a pilot plant scale is considered, a solution with 20 mass percent THF in water is investigated. In order to process a mixture of organic solvent and water from a previous polymerization synthesis resulting from a 20 L reactor within one working day, the feed mass flow into the solvent preparation plant should be 50 kg·h−1. For the design of the plant in the pilot scale, models were created with Aspen Tech® process simulation software packages (Aspen Technology Inc., Bedford, MA, USA). The distillation column is designed by employing the equilibrium stage method for the creation of a rigorous equilibrium stage RadFrac model in Aspen Plus® v.11. For the vapor permeation module, an equation-oriented model developed and implemented in Aspen Custom Modeler® by Brinkmann [27] was employed. The permeation process herein is described by the free-volume model [28]. A model for the adsorption column was also implemented in Aspen Custom Modeler® V11, applying the Langmuir adsorption isotherm. The models created are described in detail below.

4. Hybrid Process Simulation

The software Aspen Plus Dynamics® v.11 was used to represent the entire hybrid process in one simulation. For a precise representation of the distillation column, it was first configured in Aspen Plus® and then imported into Aspen Plus Dynamics® as flow driven file. For the calculation of the membrane and adsorption part of the process, the models created in Aspen Custom Modeler® were imported into Aspen Plus Dynamics®. Details of the Aspen Custom Modeler® models are provided in the Appendix A. Thus, all required models for representing the interconnected hybrid process were available in Aspen Plus Dynamics®.
The simulation of the distillation column was first performed with the software Aspen Plus® v.11 based on a rigorous model, whereby the initial values for the design were determined using the McCabe–Thiele method employing the vapor–liquid equilibrium diagram of the binary mixture that is to be investigated. Since this varies with the pressure, the operating pressure of the column must be set as an assumption first. The selection of the operating pressure depends on the distillation process and the membrane process because with the intended design, further compression of the process stream between the distillation column and membrane module can be avoided.
The analysis of the vapor–liquid equilibrium plotted in terms of THF vapor phase mass fraction vs. THF liquid phase mass fraction reveals that an increase in pressure results in a shift of the azeotropic point. The higher the pressure is, the lower the THF mass fraction in the vapor phase (see Figure 3). Hence, too high pressure is undesirable because of the shift to more significant amounts of water in the top product. Furthermore, Figure 3 indicates that the vapor in equilibrium with the feed at a THF mass fraction of 20 wt.-% is in equilibrium with a vapor containing THF in excess of 85 wt.-%. This observation led to the design of the distillation section of the process as a stripping column only since no rectifying section was necessary to increase the THF content of the distillate further. For the subsequent membrane process, higher pressure is useful for an increase in the driving force. Additionally, the working range of the membrane module, including the thermal stability of the membrane material, and the energy costs had to be considered when fixing the operating pressure and hence minimum temperature for operation of the vapor permeation module. For this hybrid process, a pressure of 4 bar was selected to support good separation results in the vapor permeation stage, ensuring the stability of the membrane material and having a high enough THF content in the distillate.
Another important aspect in selecting the operating pressure is the existence of a miscibility gap in the THF–H2O mixture. Similar to the azeotropic point of this mixture, the miscibility gap depends on the pressure [29]. However, since the first step of the stripping process does not exceed the value of 86 mass percent THF in the vapor phase that lies before the miscibility gap, which starts at a value of 89 mass percent THF, and the separation mechanisms of the following process steps operate independently of the vapor–liquid equilibrium, choosing the hybrid process to purify this mixture, the miscibility gap does not pose a problem.
After setting the working pressure, the McCabe–Thiele method can be applied to the corresponding vapor–liquid equilibrium diagram to determine a first estimate for the number of stages of the stripping column. As initial values for the calculation of the stripping column, a number of five stages was selected.
In the following, Aspen Plus® v.11 software was used for creating a rigorous equilibrium stage model based on the RadFrac model. The process was designed with regard to the composition of the bottom product. For this purpose, a water mass fraction in the bottom product of 99.9 wt.-% was assumed, whereas the THF composition of the distillate was assumed to be in equilibrium with the liquid exiting the feed stage, resulting in a vapor mass fraction of 88.57 wt.-% (Figure 3). A material balance around the column results in a bottoms mass flow rate of 38.75 kg·h−1 that was employed as specification for the RadFrac model. A total number of five stages including the reboiler was simulated, although a McCabe–Thiele design indicated that a number of four stages was sufficient. The additional stage was employed in order to provide flexibility in the intended operation of the process in a research environment.
The temperature inside the stripping column decreases from bottom to top and varies between 146 and 109 °C. The main points of the column’s configuration are summarized in Table 1. With this design of the stripping column, it is possible to produce a recycled water stream with a purity of 99.8 mass percent from the feed stream with a water content of 80 mass percent. In addition to the recycled water, this process step results in a THF–H2O mixture containing 11.7 mass percent water. The further purification of the mixture takes place in the following process steps.
In addition, the size of the distillation plant was checked to see if it meets the requirements of an application on a pilot plant scale. Therefore, the diameter of the required column was determined based on the simulation results with Aspen Plus®. The following equation was used to calculate the column diameter d C [30]:
d C = 0.00542 · G ˙ · T p · w F
The pressure p and temperature T are given by the operating conditions and the flow rate G ˙ is derived from the results of the previous Aspen Plus® simulation, which specifies the configuration of the distillation column. In order to determine the flood point velocity w F , it is necessary to consider further mathematical coherences and parameters, which express the properties of the structured packing (porosity ε , hydraulic diameter of the structured packing d h , resistance coefficient of the dry packing c R ) and the droplet diameter d D r , see Equation (2) [30,31]. The calculation of the flood point velocity w F is based on the drop–float–bed model [32]:
w F = 0.565 · c R 1 6 · ( d h d D r ) 1 / 4 · ε 6 / 5 · [ d D r · Δ ϱ · g ϱ v ] 1 / 2 ,
where Δ ϱ is the difference between liquid and vapor density, g is the gravity, and ϱ v is the vapor density.
Choosing a structured packing Mellapak 752.Y from Sulzer Chemtec [33,34] results in a column diameter of about 7 cm. The height of the column is determined by the number of stages and the packing material, which specifies the number of stages per meter column height. For the selected column material and number of stages, determined in the Aspen Plus® simulation, the column height is approximately 1.6 m. With these dimensions, this column is well suited for operation in the intended hybrid process on a pilot plant scale.
In order to ensure that the feed to the stripping column is in saturated liquid state, a heat exchanger was included into the Aspen Plus® simulation. The heat exchanger, in combination with the stripping column, was exported into a file format readable by Aspen Plus Dynamics® in order to combine this detailed process description with the following process steps.
In the second process stage, the top product of the distillation process is further purified by vapor permeation. This kind of membrane process is characterized by the fact that the feed, retentate, and permeate are in the vapor phase. This is an advantage compared to the related pervaporation process in combination with distillation, in which there would be a phase transformation within the membrane modules, i.e., the evaporation of the permeate, which leads to a cooling effect. This cooling would have to be counteracted by intermediate heating, in case of using a pervaporation process with several modules connected in series. Therefore, by choosing vapor permeation as process, additional energy for reheating can be saved. A common characteristic of the two processes is that the driving force, as for every membrane process, can be described by the difference of chemical potential. It can be shown that this can be replaced by the difference in fugacities for vapor permeation. For the introduced process, a pressure of 4 bar on the feed side and a vacuum of 0.05 bar absolute on the permeate side are employed to generate the driving force. Mass flow, composition, temperature, and pressure of the feed stream correspond to the top product of the previous distillation process.
For the simulation, a model developed by Brinkmann [27] was implemented in the Aspen Plus Dynamics® file, which already contains the distillation part. This model was validated with experimental results collected with a pilot plant for ethanol dehydration. Changing variables along the membrane area, such as concentration and mass flow, required the solution of the system of differential and algebraic equations containing, e.g., material balance and permeation relationship. For the simulation of the membrane module, pressure drops on feed and permeate side were neglected, and no impact of concentration polarization and isothermal operation were assumed. The permeation through the membrane is described by
n . i P · R · T S p S = V ˙ S T P , i P A M = L i · ( f i , H M f i , L M ) ,
where n . i P is the permeate molar flux of component i , R is the universal gas constant, T S is the standard temperature, p S is the standard pressure, V ˙ S T P , i P is the permeate volumetric flow rate at standard temperature and pressure of component i , A M is the membrane area, L i is the permeance of component i , and f i is the fugacity of component i for the high pressure ( H M ) and low pressure ( L M ) side of the membrane.
The permeance is determined by the free-volume model (FVM), which was developed by Fang et al. [28] as a form of the solution–diffusion model for swelling polymers, defined as follows:
L i = L i 0 · e x p ( j = 1 n c ( σ i σ j ) 2 · m j · 0.5 · ( f i , H M + f i ,   LM ) ) ,
where L i 0 describes the permeance of a component i for the case of infinite dilution and σ is the Lennard–Jones coefficient.
In order to assess the feasibility of this process, the parameters for the permeance L i are based on investigations of the permeation of isopropanol–water mixtures with a polymer membrane having poly(vinyl alcohol) as the selective layer [35]. For this study, the permeation behavior of THF is assumed to be represented by that of isopropanol. The parameters used in the free-volume model are summarized in Table 2.
Condition for the design of the membrane module is a purification of the THF stream in the retentate stream to 99 mole percent. The result is a membrane area of 517 cm² under the design conditions summarized in Table 3. This size is well suited for implementation in a pilot plant. A membrane module type that can be applied for the flat sheet membranes is an envelope type membrane module specifically designed for pilot plant application and introduced, e.g., in [38]. For this membrane area, a feed stream with a THF concentration of 65.15 mole percent results in a THF concentration in the permeate of 6.5 mole percent. This value is similar to the composition of the feed stream into the stripping column so that it can be recycled to the column after condensation for further processing without significantly affecting the operating conditions of the column.
About 99 mass percent of the water stream fed to the pilot plant is already separated from the organic solvent in the first distillation step. Nevertheless, the results confirm that ultrafine purification is not possible with the membrane vapor permeation separation process alone since the decrease of water driving force toward the retentate side renders this impossible. The composition and mass flow of the retentate is used in the following as input flow for the adsorption column.
As last step of the hybrid process, a model for an adsorption column is added to the Aspen Plus Dynamics® simulation to examine the scale and the time frame for a high-quality purification of the THF product stream from the combination of distillation and vapor permeation process. To calculate the breakthrough curve, the dimension of the column must be specified. Based on the assumption that the stream resulting from the previous simulation of the combined process of stripping column with membrane module should be purified at a typical liquid adsorption superficial velocity value of 0.17 cm·s−1, the diameter of the adsorption column is 5 cm. For the dehydration of the stream in the adsorption column, molecular sieve 4Å (Mole Sieve 4A) was selected as the proper adsorption material in form of spheres. The quantity of water to be adsorbed, over a period of approximately 2 h, requires a calculated quantity of 400 g of this adsorption material. This results in a minimum height of the adsorber bed of 30 cm.
The adsorption behavior of a material is described by its adsorption isotherms. The adsorption isotherm for water on Mole Sieve 4A is following the Langmuir adsorption model. This type of isotherm is valid under the assumptions of energetic equality and accessibility of all adsorption sites, no interactions between the adsorbate molecules and a monomolecular coverage of the adsorbent surface. The Langmuir coefficients were taken from Pahl [39] and are 8287 m³·mole−1 for the Langmuir constant b and 9.11 mole·kg−1 for the maximum adsorption qmon. Equation (5) describes the loading of the adsorption material q H 2 O .
q H 2 O = q m o n · b · c H 2 O 1 + b · c H 2 O
The model for the adsorption column consists of differential material balances for both components. The velocity is assumed to be constant due to the low water content, while axial dispersion is accounted for by a fixed dispersion coefficient. Mathematically, the material balance is expressed by a partial differential equation, first order with respect to time and second order with respect to axial coordinate. For each spatial location, instantaneous equilibrium between liquid and adsorbate phase is assumed according to Equation (5) for water, whereas no adsorption of THF is considered. Since no extra- or intraparticle mass transfer is considered, liquid and solid phase are combined within one material balance for each component. As already pointed out, the adsorption process is assumed to be isothermal. The created model provides the calculation of the breakthrough curve to check the process conditions and the application in the pilot plant scale. The used adsorption bed specifications, which result from the requirement of the whole hybrid process operation during one working day and the selected adsorption material for efficient dehydration of organic solvent, are summarized in Table 4. The result of the simulation is discussed in the next paragraphs. At a temperature of 25 °C and a pressure of 4 bar, the adsorption material is fully loaded after a time of about 2 h and cannot adsorb any further water molecules. To ensure sufficient dehydration by the adsorption process, the dimension of the column could, therefore, be adjusted. This can easily be realized by increasing the column height, which, even if doubled, still fits nicely into the desired dimensions of a pilot plant. The model predicts a THF stream with 100% purity up to the onset of breakthrough. Since bypass flows within the adsorption column can occur during actual operation, the purity of the product has to be verified experimentally for an accurate statement; however, due to the complexity of the construction of such a pilot-scale plant, this will be investigated in future work. The suitability of adsorption processes for the removal of water to produce high-purity organic solvents has already been proven in experimental studies [40,41].
After the adsorption column has been fully loaded, the adsorption material must be regenerated before the column can be put into operation again. The inert gas nitrogen is particularly suitable for this purpose. This is the only waste product generated during the hybrid process described because the solvent mixture itself can be separated and purified in sufficient quality for reutilization.
Figure 4 shows the whole process as simulated in Aspen Plus Dynamics® with all steps in combination. The original feed stream is first fed into a vessel (V1) that allows mixing with the permeate stream of the membrane module. Due to the similar composition of the permeate, compared to the feed of the entire process, its addition to the inlet stream of the stripping column has no significant effect on the configuration of this column. Figure 5 shows the simulation results for the molar composition of the liquid and vapor phase inside the stripping column for each stage. In Figure 5a, the values for the vapor and liquid mole fraction for THF are shown, while in Figure 5b, the results for H2O are shown. Stage one has the same composition as the feed due to the fact that in this stage, the feed inlet is located. With increasing stage number, the THF mole fraction decreases in the stripping column. The difference between stages four and five is hardly visible in the plotted graph. The simulation with Aspen Plus®, therefore, supports the previous estimation according to the McCabe–Thiele method, i.e., a large part of the purification can already be achieved with four stages. The composition in stage five shows the purity in the bottom product. The bottom product consists of 99.99 mass percent H2O. Thus, one of the main requirements for the products of this hybrid process is fulfilled since this product stream can be reused, for example, for the precipitation of new polymers. The top product continues to be a mixture of THF and H2O with above 14 mass percent H2O. It is only in the next process step that a purification larger than 99 mass percent THF is achieved by using the membrane process. This further purification of the top product by linking it to a membrane module is shown in Figure 6a, which describes the change in composition of the retentate stream. The green line shows the increasing THF mole fraction in the retentate over the membrane area. After reaching 80% of the membrane area, the retentate stream already has a purity of 0.98 mole·mole−1. The purification of the organic solvent in high quality takes place in the last process step, the adsorption process. Figure 6b shows the breakthrough curve of the adsorption column generated by the simulated hybrid process. The breakthrough of the adsorption column starts after approximately 115 min and ends after approximately 125 min. These times correspond to the required feasibility of the process, which, in a realistic scenario, should be able to recycle a quantity of 100 l of mixture of organic solvent and water within one working day. As a result, in addition to the recycling of H2O, THF can also be recovered in sufficient quality for repeated use as a solvent in polymerization processes. In this way, a zero-discharge process is created that reprocesses both components of the original mixture.

Investigation of Energy Consumption

In view of the fact that a sustainable, resource-saving process is to be developed, its energy consumption is a key parameter. Therefore, the energy consumption of the process was investigated. Since this depends to a large extent on the operating pressure, the influence of this factor was analyzed by variation. For this purpose, seven different feed pressures between 2 and 8 bar were investigated. Figure 7a shows the energy consumption for the whole process depending on different feed pressures, which includes the individual energy consumptions of pumps, heat exchangers, and the reboiler for the operation of the stripping column. The main part of the energy consumption consists of the heating power of the hybrid process with approximately 80%. Just under 20% results from the necessary cooling power, while the remaining energy is used for electrical power. Besides the total energy consumption, Figure 7a shows the data for the energy consumption of the stripping column as a single data set. This includes the reboiler of the column itself and the heater for the feed stream. It can be observed that the first process step with stripping accounts for the largest share of the energy consumption of the entire process. The membrane process, on the other hand, consumes only 9 percent of the total energy consumption due to the additionally required vacuum pump and heat exchanger. The distillation process step shows a strong dependence on the pressure in the range investigated. To account for the effect of pressure change on the efficiency of the processes, Figure 7b shows the molar fraction of THF in the top product of the stripping column and that in retentate and permeate of the membrane process. The top product of the stripping column shows a decreasing THF content with increasing pressure. The permeate stream from the membrane process shows a sharp increase in THF content from a pressure of 4 bar, while the retentate stream shows only a small increase in THF content. This is expected since a molar content of 0.99 is already present at 4 bars. The increase of the THF content in the permeate stream above a value of 6 mole percent is not desired since this stream is to be added to the feed of the stripping column and would increase the load on the column. It is not surprising that the figures shown indicate favorable operation at 4 bars since the configuration of the column and the membrane area are designed for this operating condition. Nevertheless, a consideration in comparison with the dependence of the energy consumption on the pressure is useful, since ultimately, it must be decided whether the energy consumption of the recycling process represents a significant improvement over the production of new solvent, taking into account the achieved purification target. In order to save an additional compressor upstream of the membrane process and take into account the increasing strain on membrane material with increasing pressure, the choice of a feed pressure of 4 bar represents an acceptable compromise. The configuration of the process at a feed pressure of 4 bar results in an energy consumption of 191.8 W per kg·h−1 of feed mixture to be purified. Calculating the energy required to produce one kg of pure THF from the waste stream that would otherwise be discarded, the value is 3.5 kJ·kg−1. In addition to reducing the raw material requirement, the use of the process described here for recycling and reuse of the THF and water in the preceding process of polymer synthesis can also save the energy that would otherwise be required for the proper disposal of organic solvent.

5. Conclusions

A hybrid process to separate and purify a water–solvent mixture was developed based on simulation models. The individual process steps were simulated with the simulation software from Aspen Tech® with the aim to develop a pilot-scale plant. This plant should meet the requirements of producing solvent in reusable quality for polymer synthesis processes with the aim to carry out more environmentally friendly and sustainable polymer synthesis on a large—and small—scale and hence approach a closed cycle, zero-discharge process. A process consisting of three process steps containing distillation, vapor permeation, and adsorption was chosen. Hence, the separation principles of the individual process steps are based on different mechanisms. The combination of stripping column with downstream vapor permeation allows us to overcome the azeotropic point of the used sample mixture THF–H2O. The less highly purified permeate stream from the single-stage membrane process was returned to the stripping column by adding it to the column feed stream. The result of this process combination, which was simulated with Aspen Plus Dynamics®, is a THF stream with a purity of 99.7 mass percent. Moreover, the bottom product of the stripping column is a water stream with 99.6 mass percent purity. Therefore, this stream is of sufficient quality for reuse in the precipitation process of polymerization products. The retentate stream of the membrane module is further purified with adsorption in the last process step. This ensures that THF is available in high purity. The calculations show that a separation of a THF–H2O mixture of 112 L with 20 mass percent THF can be realized in a pilot plant within one working day. With the development of this process for the separation of an exemplary mixture of organic solvent and water, which accumulates as a waste product in many industrial processes, into output streams of high purity, a significant contribution can be made to a sustainable transformation of the chemical industry. In the future, the construction of the plant is planned so that the results of the simulation can be compared with real measured values and the models can be adjusted if necessary. Furthermore, the models created should be used to evaluate the solution from separation tasks of other organic solvent mixtures and prove in which quality a separation is possible with this hybrid process.

Author Contributions

Conceptualization, T.B., K.S. and P.G.; methodology, T.B. and K.S.; data curation, K.S.; writing—original draft preparation, K.S. and P.G.; writing—review and editing, K.S., T.B. and P.G.; visualization, K.S.; supervision, P.G. and T.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.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank Heiko Notzke for the 3D symbol image of the hybrid process used in the scheme of the closed-loop process and Holger Pingel for fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The equations used in the Aspen Custom Modeler® model, which is implemented in the Aspen Plus Dynamics® file for the simulation of the vapor permeation module part in the whole hybrid process are as follows.
The differential material balance for the membrane module was approximated using a first-order upstream finite differencing scheme, where the membrane area A is discretized into a number of elements i = 1 , , n z .
( n ˙ k , i R n ˙ k , i 1 R ) + n . k , i M · Δ A = 0
n . k , i M · R · T S p S = L k , i · ( f k , i R f k , i P )
n . k , i M = n . k , i P
where n ˙ k , i R is the molar flow rate of component k exiting element i on the retentate side R, n . k , i M is the molar flux of component i through the membrane at element i, Δ A = A n z is the discretized area element, R is the universal gas constant, T S is the temperature at standard conditions, p S is the pressure at standard conditions, L k , i is the permeance which is calculated using Equation (4) employing local composition variables, f k , i R and f k , i R   are the retentate and permeate side fugacities of component k at discretization element i, respectively, and n . k , i P is the local permeate flux perpendicular to the membrane. For discretization element i = 1   the retentate molar flow rates at i 1 = 0   are calculated from the feed definition, i.e., the stripping column’s vapor phase product stream. The conditions on the permeate side are specified by supplying a value for the permeate pressure. Required compositions on retentate and permeate sides are calculated from molar retentate side flow rates and transmembrane fluxes, respectively.
The equations used in the Aspen Custom Modeler® model for the description of a simple, isothermal adsorption column where extra- and intraparticle mass transfer are considered negligible, and the superficial velocity is considered constant are detailed as follows. The model is implemented in the Aspen Plus Dynamics® file for the simulation of the adsorption column module in the whole process:
The material balance is numerically solved by using a first-order upwind finite differencing scheme for the first-order spatial derivatives and the central differencing scheme for the second-order spatial derivatives. For this, the bed height is discretized into nz nodes. The derivatives with respect to time are solved using the implicit Euler, variable step-size method implemented in Aspen Custom Modeler® and Aspen Dynamics®.
n ˙ F e e d = Σ ( c k , 0 ) · v · π 4 · D C 2
n ˙ P r o d u c t = Σ ( c k ,   n z ) · v · π 4 · D C 2
c k , n z c k , n z 1 Δ z = 0
ε · d c k , i d t + v · c k , i c k , i 1 Δ z D a x · c k , i + 1 2 · c k , i + c k , i 1 ( Δ z ) 2 + ρ B u l k · d q k , i d t = 0
where n ˙ F e e d and n ˙ P r o d u c t are feed and product molar flow rates of the adsorption column, c k , i is the molar concentrations of component k at node i, v is the superficial velocity, D c is the adsorption column diameter, ε is the bed voidage, z is the spatial co-ordinate, D a x is the axial dispersion coefficient, ρ B u l k is the adsorbent bulk density, and q k , i the loading of component k at node i calculated according to Equation (5).

References

  1. Statista. Available online: https://www.statista.com/statistics/1009366/market-value-solvent-recycling-equipment-globally/ (accessed on 27 October 2020).
  2. Gómez, P.; Gil, I. Simulation of the tetrahydrofuran dehydration process by extractive distillation in Aspen Plus. Latin Amer. Appl. Res. 2009, 39, 275–284. [Google Scholar]
  3. Mencarini, J., Jr.; Coppola, R.; Slater, C.S. Separation of tetrahydrofuran from aqueous mixtures by pervaporation. Sep. Sci. Technol. 1994, 29, 465–481. [Google Scholar] [CrossRef]
  4. Deorukhkar, O.A.; Deogharkar, B.S.; Mahajan, Y.S. Purification of tetrahydrofuran from its aqueous azeotrope by extractive distillation: Pilot plant studies. Chem. Eng. Process. Process Intensif. 2016, 105, 79–91. [Google Scholar] [CrossRef]
  5. Xu, S.; Wang, H. Separation of tetrahydrofuran–water azeotropic mixture by batch extractive distillation process. Chem. Eng. Res. Des. 2006, 84, 478–482. [Google Scholar] [CrossRef]
  6. Gu, J.; You, X.; Tao, C.; Li, J.; Shen, W.; Li, J. Improved design and optimization for separating tetrahydrofuran–water azeotrope through extractive distillation with and without heat integration by varying pressure. Chem. Eng. Res. Des. 2018, 133, 303–313. [Google Scholar] [CrossRef]
  7. Lipnizki, F.; Field, R.W.; Ten, P.-K. Pervaporation-based hybrid process: A review of process design, applications and economics. J. Membr. Sci. 1999, 153, 183–210. [Google Scholar] [CrossRef]
  8. Conde-Mejía, C.; Jiménez-Gutiérrez, A. Analysis of ethanol dehydration using membrane separation processes. Open Life Sci. 2020, 15, 122–132. [Google Scholar] [CrossRef]
  9. Chapman, P.D.; Oliveira, T.; Livingston, A.G.; Li, K. Membranes for the dehydration of solvents by pervaporation. J. Membr. Sci. 2008, 318, 5–37. [Google Scholar] [CrossRef]
  10. Morigami, Y.; Kondo, M.; Abe, J.; Kita, H.; Okamoto, K. The first large-scale pervaporation plant using tubular-type module with zeolite NaA membrane. Sep. Purif. Technol. 2001, 25, 251–260. [Google Scholar] [CrossRef]
  11. Melin, T.; Rautenbach, R. Membranverfahren: Grundlagen der Modul-und Anlagenauslegung; Springer: Berlin, Germany, 2007. [Google Scholar]
  12. Daviou, M.C.; Hoch, P.M.; Eliceche, A.M. Design of membrane modules used in hybrid distillation/pervaporation systems. Ind. Eng. Chem. Res. 2004, 43, 3403–3412. [Google Scholar] [CrossRef]
  13. Brinkmann, T.; Ebert, K.; Pingel, H.; Wenzlaff, A.; Ohlrogge, K. Prozessalternativen durch den Einsatz organisch-anorganischer Kompositmembranen für die Dampfpermeation. Chemie Ingenieur Technik 2004, 76, 1529–1533. [Google Scholar] [CrossRef]
  14. Roth, T.; Kreis, P.; Górak, A. Process analysis and optimisation of hybrid processes for the dehydration of ethanol. Chem. Eng. Res. Des. 2013, 91, 1171–1185. [Google Scholar] [CrossRef]
  15. Dijkstra, M.F.J.; Brinkmann, T.; Ebert, K.; Ohlrogge, K. Vorrichtung zur Trennung von flüssigen und/oder gasförmigen Gemischen. Patent DE10333049B3, 25 November 2004. Available online: https://patentimages.storage.googleapis.com/93/e6/a9/569ad857ddd5ad/DE10333049B3.pdf (accessed on 1 February 2021).
  16. Su, Y.; Jin, S.; Zhang, X.; Shen, W.; Eden, M.R.; Ren, J. Stakeholder-oriented multi-objective process optimization based on an improved genetic algorithm. Comput. Chem. Eng. 2020, 132, 106618. [Google Scholar] [CrossRef]
  17. Erdmenger, T.; Guerrero-Sanchez, C.; Vitz, J.; Hoogenboom, R.; Schubert, U.S. Recent developments in the utilization of green solvents in polymer chemistry. Chem. Soc. Rev. 2010, 39, 3317–3333. [Google Scholar] [CrossRef]
  18. Lauterbach, F.; Rubens, M.; Abetz, V.; Junkers, T. Ultrafast PhotoRAFT Block Copolymerization of Isoprene and Styrene Facilitated through Continuous-Flow Operation. Angew. Chem. Int. Ed. 2018, 57, 14260–14264. [Google Scholar] [CrossRef]
  19. Kandelhard, F.; Schuldt, K.; Schymura, J.; Georgopanos, P.; Abetz, V. Model-Assisted Optimization of RAFT Polymerization in Micro-Scale Reactors—A Fast Screening Approach. Macromol. React. Eng. 2021, 2000058. [Google Scholar] [CrossRef]
  20. Nieswandt, K.; Georgopanos, P.; Abetz, C.; Filiz, V.; Abetz, V. Synthesis of poly (3-vinylpyridine)-block-polystyrene diblock copolymers via surfactant-free RAFT emulsion polymerization. Materials 2019, 12, 3145. [Google Scholar] [CrossRef] [Green Version]
  21. Nieswandt, K.; Georgopanos, P.; Abetz, V. Well-defined polyvinylpyridine-block-polystyrene diblock copolymers via RAFT aqueous-alcoholic dispersion polymerization: Synthesis and isoporous thin film morphology. Polym. Chem. 2021, 12, 2210–2221. [Google Scholar] [CrossRef]
  22. Rangou, S.; Buhr, K.; Filiz, V.; Clodt, J.I.; Lademann, B.; Hahn, J.; Jung, A.; Abetz, V. Self-organized isoporous membranes with tailored pore sizes. J. Membr. Sci. 2014, 451, 266–275. [Google Scholar] [CrossRef] [Green Version]
  23. Chapman, P.D.; Tan, X.; Livingston, A.G.; Li, K.; Oliveira, T. Dehydration of tetrahydrofuran by pervaporation using a composite membrane. J. Membr. Sci. 2006, 268, 13–19. [Google Scholar] [CrossRef]
  24. Capello, C.; Fischer, U.; Hungerbühlera, K. What is a green solvent? A comprehensive framework for the environmental assessment of solvents. Green Chem. 2007, 9, 927–934. [Google Scholar] [CrossRef]
  25. Fontalvo, J.; Cuellar, P.; Timmer, J.; Vorstman, M.; Wijers, J.; Keurentjes, J. Comparing pervaporation and vapor permeation hybrid distillation processes. Ind. Eng. Chem. Res. 2005, 44, 5259–5266. [Google Scholar] [CrossRef]
  26. Pettersen, T.; Lien, K.M. Design of hybrid distillation and vapor permeation processes. J. Membr. Sci. 1995, 102, 21–30. [Google Scholar] [CrossRef]
  27. Brinkmann, T.; Dijkstra, M.; Ebert, K.; Ohlrogge, K. Improved simulation of a vapour permeation module. J. Chem. Technol. Biotechnol. Int. Res. Process Environ. Clean Technol. 2003, 78, 332–337. [Google Scholar] [CrossRef]
  28. Fang, S.M.; Stern, S.A.; Frisch, H.L. A “free volume” model of permeation of gas and liquid mixtures through polymeric membranes. Chem. Eng. Sci. 1975, 30, 773–780. [Google Scholar] [CrossRef]
  29. Schneider, G.M. Aqueous solutions at pressures up to 2 GPa: Gas–gas equilibria, closed loops, high-pressure immiscibility, salt effects and related phenomena. Phys. Chem. Chem. Phys. 2002, 4, 845–852. [Google Scholar] [CrossRef]
  30. Sattler, K. Thermische Trennverfahren: Grundlagen, Auslegung, Apparate, 3. Auflage; Wiley-VCH: Weinheim, Germany, 1995. [Google Scholar]
  31. Maćkowiak, J. Fluiddynamik von Füllkörpern und Packungen; Springer: Berlin, Germany, 2003. [Google Scholar]
  32. Maćkowiak, J. Fluiddynamik von Kolonnen mit modernen Füllkörpern und Packungen fur Gas./Flüssigkeitssysteme; Otto Salle: Berlin, Germany, 1991. [Google Scholar]
  33. Structured Packings for Distillation, Absorption and Reactive Distillation. Available online: https://www.sulzer.com/-/media/files/products/separation-technology/liquid_liquid_extraction/brochures/structured_packings.ashx (accessed on 1 February 2021).
  34. Olujić, Ž.; Behrens, M.; Spiegel, L. Experimental characterization and modeling of the performance of a large-specific-area high-capacity structured packing. Ind. Eng. Chem. Res. 2007, 46, 883–893. [Google Scholar] [CrossRef]
  35. Fröhlich, N. Charakterisierung einer Dampfpermeationsmembran für die Entwässerung organischer Flüssigkeiten. Diploma Thesis, HAW Hamburg, Library Institute of Membrane Research, Helmholtz-Zentrum Hereon, 2010. [Google Scholar]
  36. Drabowicz, W. Molecular dynamics study of the structural and dynamical properties of liquid tetrahydrofuran. Zeitschrift Für Naturforschung A 1990, 45, 1342–1344. [Google Scholar] [CrossRef]
  37. Poling, B.E.; Prausnitz, J.M.; O’Connell, J.P. The Properties of Gases and Liquids; Mc-Graw Hill: New York, NY, USA, 2001; Volume 5. [Google Scholar]
  38. Brinkmann, T.; Naderipour, C.; Pohlmann, J.; Wind, J.; Wolff, T.; Esche, E.; Müller, D.; Wozny, G.; Hoting, B. Pilot scale investigations of the removal of carbon dioxide from hydrocarbon gas streams using poly (ethylene oxide)–poly (butylene terephthalate) PolyActive™) thin film composite membranes. J. Membr. Sci. 2015, 489, 237–247. [Google Scholar] [CrossRef] [Green Version]
  39. Pahl, C. Adsorptive Wasserentfernung Im Unteren Ppm-Und Ppb-Bereich Zur Erzeugung Hochreiner Organischer Lösungsmittel; Shaker: Düren, Germany; Maastricht, The Netherlands, 2012. [Google Scholar]
  40. Gorbach, A.; Stegmaier, M.; Eigenberger, G. Measurement and modeling of water vapor adsorption on zeolite 4A—Equilibria and kinetics. Adsorption 2004, 10, 29–46. [Google Scholar] [CrossRef]
  41. Burrichter, B. Adsorptive Entfernung Von Wasser Im Unteren Ppm-Bereich Zur Erzeugung Hochreiner Organischer Lösungsmittel; Shaker: Düren, Germany; Maastricht, The Netherlands, 2015. [Google Scholar]
Figure 1. Scheme of a closed-loop process for the purification of a mixture of organic solvent and water resulting as waste product from the synthesis of polymers and their subsequent precipitation.
Figure 1. Scheme of a closed-loop process for the purification of a mixture of organic solvent and water resulting as waste product from the synthesis of polymers and their subsequent precipitation.
Processes 09 00729 g001
Figure 2. Flowsheet of the hybrid process for the separation of azeotropic mixtures.
Figure 2. Flowsheet of the hybrid process for the separation of azeotropic mixtures.
Processes 09 00729 g002
Figure 3. Vapor–liquid equilibrium graph for THF–H2O mixture in the pressure range of 1–10 bar, created with Aspen Plus® using the nonrandom two-liquid (NRTL) thermodynamic model. Exemplary pressures 1, 4, 10 bar depicted with three different colors.
Figure 3. Vapor–liquid equilibrium graph for THF–H2O mixture in the pressure range of 1–10 bar, created with Aspen Plus® using the nonrandom two-liquid (NRTL) thermodynamic model. Exemplary pressures 1, 4, 10 bar depicted with three different colors.
Processes 09 00729 g003
Figure 4. Process simulation in Aspen Plus Dynamics® with all subprocesses in combination.
Figure 4. Process simulation in Aspen Plus Dynamics® with all subprocesses in combination.
Processes 09 00729 g004
Figure 5. (a) THF mole fraction of the liquid and vapor phase inside the stripping column and (b) H2O mole fraction of the liquid and vapor phase inside the stripping column. As it can be inferred from the mole fractions, a large part of the purification can already be achieved with four stages.
Figure 5. (a) THF mole fraction of the liquid and vapor phase inside the stripping column and (b) H2O mole fraction of the liquid and vapor phase inside the stripping column. As it can be inferred from the mole fractions, a large part of the purification can already be achieved with four stages.
Processes 09 00729 g005
Figure 6. (a) Retentate mole fraction of THF and Water over the membrane area and (b) breakthrough curve of Water of adsorption column. It is observed from the data that the high purity of the THF can be achieved via the combination of the different sequential separation processes.
Figure 6. (a) Retentate mole fraction of THF and Water over the membrane area and (b) breakthrough curve of Water of adsorption column. It is observed from the data that the high purity of the THF can be achieved via the combination of the different sequential separation processes.
Processes 09 00729 g006
Figure 7. (a) Energy consumption of the whole process (Total) and the energy consumption of the stripping process part (Stripping column) and (b) THF mole fraction of the product streams in dependence of the used feed pressure for stripping column (top product) and membrane process (permeate and retentate stream).
Figure 7. (a) Energy consumption of the whole process (Total) and the energy consumption of the stripping process part (Stripping column) and (b) THF mole fraction of the product streams in dependence of the used feed pressure for stripping column (top product) and membrane process (permeate and retentate stream).
Processes 09 00729 g007
Table 1. Configuration stripping column with Aspen Plus® RadFrac model.
Table 1. Configuration stripping column with Aspen Plus® RadFrac model.
Value
Number of stages5
Feed inletStage 1 (stripping column)
Bottoms mass flow rate38.7 kg·h−1
Table 2. Parameters of the free-volume model (FVM) [35] for the 20 mass percent THF–H2O mixture separation.
Table 2. Parameters of the free-volume model (FVM) [35] for the 20 mass percent THF–H2O mixture separation.
FVM VariableUnitTHFH2O
L i 0 m³·m−2·h−1·bar−10.1678.12
m i bar−101.25
σ i [36,37]Å5.502.64
Table 3. Input data vapor permeation model.
Table 3. Input data vapor permeation model.
Process VariableUnitValue
pPermeatembar50
mFeedkg·hr−111.25
pFeedbar4
yFeed,THFmole percent65.14
Table 4. Adsorption bed specifications parameters.
Table 4. Adsorption bed specifications parameters.
UnitValue
Bed heightm0.30
Bed inner diameterm0.05
Bed void fraction-0.55
Bulk density molecular sievekg·m³670
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schuldt, K.; Brinkmann, T.; Georgopanos, P. Zero-Discharge Process for Recycling of Tetrahydrofuran–Water Mixtures. Processes 2021, 9, 729. https://doi.org/10.3390/pr9050729

AMA Style

Schuldt K, Brinkmann T, Georgopanos P. Zero-Discharge Process for Recycling of Tetrahydrofuran–Water Mixtures. Processes. 2021; 9(5):729. https://doi.org/10.3390/pr9050729

Chicago/Turabian Style

Schuldt, Karina, Torsten Brinkmann, and Prokopios Georgopanos. 2021. "Zero-Discharge Process for Recycling of Tetrahydrofuran–Water Mixtures" Processes 9, no. 5: 729. https://doi.org/10.3390/pr9050729

APA Style

Schuldt, K., Brinkmann, T., & Georgopanos, P. (2021). Zero-Discharge Process for Recycling of Tetrahydrofuran–Water Mixtures. Processes, 9(5), 729. https://doi.org/10.3390/pr9050729

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