Energy Conservation in Ethanol-Water Distillation Column with Vapour Recompression Heat Pump

interest. However, the analysis of the bioethanol process shows that distillation is still the most widely used.


Introduction
Ethanol or ethyl alcohol CH 3 CH 2 OH, a colorless liquid with characteristic odor and taste; commonly called grain alcohol has been described as one of the most exotic synthetic oxygen-containing organic chemicals because of its unique combination of properties as a solvent, a germicide, a beverage, an antifreeze, a fuel, a depressant, and especially because of its versatility as a chemical intermediate for other organic chemicals. Ethanol could be derived from any material containing simple or complex sugars. The sugar-containing material is fermented after which the liquid mixture of ethanol and water is separated into their components using distillation.
Distillation is the most widely used separation operation in chemical and petrochemical industries accounting for around 25-40% of the energy usage. One disadvantage of distillation process is the large energy requirement. Distillation consumes a great deal of energy for providing heat to change liquid to vapour and condense the vapour back to liquid at the condenser. Distillation is carried out in distillation columns which are used for about 95% of liquid separations and the energy use from this process accounts for an estimated 3% of the world energy consumption (Hewitt et al, 1999). It has been estimated that the energy use in distillation is in excess. With rising energy awareness and growing environmental concerns there is a need to reduce the energy use in industry. The potential for energy savings therefore exists and design and operation of energy efficient distillation systems will have a substantial effect on the overall plant energy consumption and operating costs.
The economic competitiveness of ethanol has been heightened by concerns over prices and availability of crude oil as well as greenhouse gas emissions which have stimulated interest in alternatives to crude oil to provide for automotive power and also by the use of bioethanol in the production of hydrogen for fuel cells. Therefore, there is the need to explore ways of producing ethanol at competitive costs by the use of energy efficient processes. To cope with the high energy demand and improve the benefits from the process, the concept of polygeneration and hydrothermal treatment especially when dealing with small scale ethanol plants is fast gaining interest. However, the analysis of the bioethanol process shows that distillation is still the most widely used.
Over the years, there have been many searches for lower energy alternatives or improved efficiencies in distillation columns. One such search led to the use of heat pumps, the idea which was introduced in the 1950s. Also, Jorapur and Rajvanshi (1991) have used solar energy for alcohol distillation and concluded that it was not economically viable. Heat pumping, however, has been known as an economical energy integration technology for reduction in consumption of primary energy and to minimize negative impact of large cooling and heating demands to the environment. One of the heat pump cycles which have been widely studied is the recompression of the vapours where the reboiler is heated by adding a compressor to the column to recover some of the heat lost in the distillate.
Most studies have concluded that heat pumping is an effective means of saving energy and reducing column size without estimating the actual energy consumption and the parameters that are likely to have significant effect on energy consumption. Estimating the actual energy consumption is an important aspect towards the determination of the viability of the system in ethanol-water separation.
The purpose of this chapter was to study how previously neglected and/or assumed values of different parameters (the pressure increase across the compressor was ignored, column heat loss was assumed to be 10% of the reboiler heat transfer rate, and the overall heat transfer coefficient was determined without considering it as an explicit function of dimensionless numbers, and its dependence on fluid viscosity and thermal conductivity neglected) affect the process efficiency, energy consumption and the column size of a vapour recompression heat pump.
The conventional heat pumps are electrically driven vapour recompression types, which work on the principle that a liquid boils at a higher temperature if its pressure is increased. A low-pressure liquid passes into the evaporator, where it takes in heat causing the liquid to boil at low temperature. The low-pressure vapour is passed to the compressor where it is compressed by the application of work to a higher pressure. The resulting high pressure vapour flows to the condenser where it condenses, giving up its latent heat at a high temperature, before expanding back to a low pressure liquid.
The heat pump cycle may be connected to a distillation column in three ways (Fonyo and Benko, 1998) . The simplest alteration is to replace steam and cooling water with refrigerant (closed system). The other two types of heat pump system apply column fluids as refrigerant . When the distillate is a good refrigerant the vapour recompression can be used. If the bottom product is a good refrigerant the bottom flashing can be applied.
In this work, the direct vapour recompression system is studied due to its good economic figures ( Emtir et al, 2003). Also the vapour recompression is the most suitable as the boiling points of both key components (ethanol and water) are close to each other (Danziger, 1979) and the appropriate heat transfer medium (ethanol vapour) is available.

Use of vapour recompression in distillation columns
Vapour recompression system has been extensively studied since 1973, the year of drastic rise in energy (Null, 1976). The vapour recompression system is accomplished by using compressor to raise the energy level of vapour that is condensed in reboiler-condenser by exchange of heat with the bottoms. The condensate distillate is passed into reflux drum while the bottom product is vaporised into the column.
Vapour recompression consists of taking the overhead vapour of a column, condensing the vapour to liquid, and using the heat liberated by the condensation to reboil the bottoms liquid from the same column. The temperature driving force needed to force heat to flow from the cooler overhead vapours to the hotter bottoms product liquid is set up by either compressing the overhead vapour so that it condenses at a higher temperature, or lowering the pressure on the reboiler liquid so it boils at a lower temperature, then compressing the bottoms vapour back to the column pressure. While conventional column has a separate condenser and reboiler, each with its own heat transfer fluid such as cooling water and steam, the vapour recompression column has a combined condenser-reboiler, with external heat transfer fluids.
The advantage of vapour recompression lies in its ability to move large quantities of heat between the condenser and reboiler of the column with a small work input. This results from cases where there is only a small difference between the overhead and bottoms temperature. Also, the temperature, and therefore the pressure, at any point may be set where desired to achieve maximum separation. This effect is of particular importance where changing the pressure affects the relative volatility. By operating at more favourable conditions, the reflux requirement can be reduced and therefore the heat duties. These advantages can reduce a large amount of energy. ambient conditions, is preheated with bottom product and condensate in heat exchangers, preheaters PH1 and PH2, and fed to the column. An auxiliary reboiler (AR) is used to start the unit. This reboiler supplies the auxiliary heat duty, which is the heat of vaporization because the main reboiler can work only if there is some compressed vapour already available. The overhead vapours from the top are compressed in the compressor (CP) up to the necessary pressure in such a way that its condensing temperature is greater than the boiling temperature of the column bottom product. The vapour is then condensed by exchanging heat within the tubes of the reboiler-condenser (RC). In a condenser, the inlet temperature is equal to the outlet temperature. Ethanol vapour will only lose its latent heat of condensation. At the same time, the cold fluid (ethanol-water mixture) in the reboiler will absorb this latent heat and its temperature will increase to boil up the mixture to temperature T CEV. The liberated latent heat of condensation provides the boil-up rate to the column while the excess heat extracted from the condensate is exchanged with the feed in preheater PH2. The condensate, which is cooled in the cooler (CL) up to its bubble point at the column operating pressure, expands through the throttling valve (TV) at the same pressure and reaches the flash tank (FT). After expansion, the output phases are a vapour phase in equilibrium with a liquid phase. One part of the product in the liquid phase is removed as distillate and stored in the tank (DST), while the remainder is recycled into the column as reflux L1. The excess of vapour is recycled to the compressor.

Methodology
Like this work, nearly all publications in this field are based on modelling and simulation (Brousse et al., 1985;Ferre et al., 1985;Collura and Luyben, 1988;Muhrer et al, 1990;Oliveira et al. 2001). The mathematical modeling of the distillation system is derived by applying energy, composition and overall material balances together with vapour-liquid equilibrium under some assumptions (see Muhrer et al, 1990 andEnweremadu, 2007). These and other assumptions are aimed at simplifying the otherwise cumbersome heat-and mass-transfer, and the fluid flow equations Mori et al (2002).

Calculation of the distillation column
In this system, there is a direct coupling between the distillation column and the rest of the system, as the heat pump working fluid is the column's own fluid which, is a binary mixture of ethanol and water at composition X D . Therefore, the set of equations are not solved separately as in distillation column assisted by an external heat pump.
The detailed calculation of the overall material and component material balance such as the bottom flow rate, B and distillate flow rate, D; reflux ratio, R r ; the molar vapour flow rate which leaves the column top and feeds the condenser, V 1 ; feed vapour flow rate, V F ; feed vapour fraction, q; vapour molar flow rate remaining at the bottom of the column, L 2 are given (see Enweremadu, 2007).
The overall (global) energy balance equation applied to a control volume comprising the distillation column and the feed pre-heaters provides the total energy demand in the reboiler: where Q reb is the total heat load added to the reboiler, Q losses represents the heat losses in the column, which are to be determined; Q 1 and Q 2 are the heat loads of the pre-heaters; h LV,e is latent heat of vaporisation downstream of throttling valve; h D , is the enthalpy of the distillate; h B is the enthalpy of the bottom product; h F , is the enthalpy of the feed. The details of the mass balance variables are determined in Enweremadu (2007).
The first step in the design of a distillation column is the determination of the number of theoretical plates required for the given separation. The theoretical trays are numbered from the top down, and subscripts generally indicate the tray from which a stream originates with n and m standing for rectifying and stripping sections respectively. The design procedure for a tray distillation column consists of determining the liquid and vapour composition or fraction from top to bottom, along the column. In calculating the composition profile of the column two equations relating liquid mole fraction to temperature and vapour mole fraction to the liquid fraction are used. The compositions at the top (X D ) and bottom (X B ) of the column are previously pre-established data. In this work, the minimum number of theoretical stages (N min ) is calculated using Fenske's equation: where α is the relative volatility in the column. The actual number of plates is given by: where T  is the tray efficiency.

Heat losses from distillation column
The heat loss from the distillation column is the main factor that affects heat added and removed at the reboiler and condenser respectively. Most distillation columns operate above ambient temperature, and heat losses along the column are inevitable since insulating materials have a finite thermal conductivity. Heat loss along the distillation column increase condensation and reduces evaporation. Thus, the amount of vapour diminishes in the upward part of the column, where the flow of liquid is also less than at the bottom.
To prevent loss of heat, the distillation column should be well insulated. Insulation of columns using vapour recompression varies with the situation. Where the column is hot and extra reboiler duty is used, the column should be insulated (Sloley, 2001). The imperfect insulation of the column causes some heat output.
In determining the heat loss from the distillation column, it is assumed that the temperature is uniform in the space between two plates. The heat transfer between the column wall and the surrounding is then determined from the well-known relationship for overall heat transfer coefficient: where U p , the overall heat transfer coefficient is given by Gani, Ruiz and Cameron (1986), as where the temperature difference, p T  , is given as h o , the heat transfer coefficient between the surroundings and the column external surface, is given as h i is the heat transfer coefficient inside the column; K p is the thermal conductivity of the tray material; A o is the external area of heat exchange; A i is the internal area of heat exchange; A m is the logarithmic mean area; t ins is the thickness of insulation.
The heat output is calculated with the general expression for convection around cylindrical objects.
The column inner surface heat transfer resistance is neglected as the heat transfer coefficient for condensing vapor is large and therefore will have little effect on the overall heat transfer.
Based on the assumptions in Enweremadu (2007), the heat transfer due to free convection between the surroundings and the external column wall and due to conduction through the insulation materials is predicted.
Also, from geometry of the insulated cylinder ( Fig.2), the external diameter of insulation is given as Details of how the logarithmic mean diameter of the insulating layer (d ins,m ), external area of heat exchange (A o ) and the logarithmic mean area (A m ) can be found in Enweremadu (2007).
From dimensional analysis, where, t ins is the thickness of insulation; K ins -thermal conductivity of the insulation materials; Nu -Nusselt number; d o -external diameter of column; T amb -temperature of the surrounding; T p -plate temperature.

Fig. 2. Hypotethical Section of the Distillation Column with Insulation
For vertical cylinders, the commonly used correlations for free convection are adapted from Rajput (2002) where Gr is the Grashof number and Pr is Prandtl.
Based on the assumptions of neglecting h i , A i and the effect of thermal resistance, equation (5) reduces to: while equation (7) is given as The heat loss from the column trays is given by The total heat loss from the column is expressed as

Calculation of heat pump and compressor parameters
The heat pump is thermodynamically linked to the column through the heat load from the pump to the column Q HPC and from the column to the pump Q CHP , and reboiler-condenser temperature. These parameters provide the basis for the heat pump calculation.
The calculation of the heat pump parameters begins with the estimation of the working fluid condensation temperature obtained from the reboiler temperature and temperature drop across the heat exchangers.
where T CEV is the column vapourization (reboiler) temperature and ΔT CHP , a pre-established mean temperature difference across the heat exchangers (temperature drop in reboilercondenser). Next is the estimation of the relevant thermodynamic properties of the working fluid. These are obtained from thermodynamic correlations.
The thermodynamic properties are determined as functions of temperature. The relationships used for calculating the working fluid density, viscosity, thermal conductivity and heat capacity for input at various locations are presented in Enweremadu (2007). The condensation pressure, P CHP is expressed as a function of condensation temperature as while the condensation pressure is determined from ideal gas equation.
The latent heat of condensation from column to heat pump is numerically exactly equal to the latent heat of vaporisation, but has the opposite sign: latent heat of vaporisation is always positive (heat is absorbed by the substance), whereas latent heat of condensation is always negative (heat is released by the substance). Latent heat of condensation is expressed as a function of condensation temperature and is determined from the relationship (Ackland, 1990): is the heat of vapourisation at the boiling point of ethanol; T C is the critical temperature; A and B are constants.
The vapour specific volume at location "a" (entrance to the compressor) is expressed as a function of column condensation temperature, T CC and pressure at the top of the column, P TOP : Since compression is polytropic, at location "b" (compressor discharge), the vapour specific volume is determined by: where n is the polytropic index and ΔP is the pressure increase across the compressor.
The vapour specific enthalpy at "a" is a function of top pressure,P TOP and top temperature, T TOP .
h a =f(T TOP, P TOP ) www.intechopen.com The vapour specific enthalpy at "b" may be d e t e r m i n e d a s a f u n c t i o n o f c o m p r e s s o r discharge temperature T b and condensation pressure P CHP but in this study, it is determined by the development of numerical computation with calculations utilizing the Redlich -Kwong equation of state. The Redlich -Kwong equation of state is given as where Z is the compressibility factor, Cp o is the molar specific heat capacities of gases at zero pressure given as a polynomial in temperature.
Equations 24 The specific liquid enthalpies have been assumed to be simple functions of temperature.
The liquid specific enthalpy at location "c" is determined from EZChemDB Thermodynamic Properties while the liquid specific enthalpy at location e (at the exit of the throttling valve) is determined from h L,e =f(T CC ) The temperature at compressor discharge is determined from the knowledge of the compressor efficiency. The ideal discharge temperature (the temperature that gives an overall change in entropy equal to zero) is calculated before correcting with the compressor efficiency.
The dryness fraction after the isenthalpic expansion is given by where the molar latent heat of vaporization at location "e" is adapted from Ackland (1990): Since this is a throttling process, T d = T e and h d = h e The molar vapour flow rate which is recycled in the flash tank and conveyed to the compressor is calculated by Therefore, the molar flow rate across the compressor is expressed as While the dryness fraction at condenser exit is determined by where Q 23 is the distribution of excess heat between the pre-heater Q 2 and the cooler Q 3 and Cp L is the molar specific heat of the working fluid in the liquid phase.
The energy balance, applied to the heat pump working fluid, yields the available energy for exchange at the condenser, as follows: www.intechopen.com A comparison is made between this energy available at the condenser, Q cd , with the energy required by the column reboiler, Q reb . This brings about the following heat load control.
i. If the rate of energy available at the heat pump condenser, Q cd , is greater than the rate of energy required by the reboiler Q reb , then the condenser gives up Q reb to the reboiler and the remaining energy is conveyed to the preheaters (Q 2 ) and cooler (Q 3 ) ii. But if Q cd is smaller than or equal to Q reb , then all energy available is transferred to the reboiler and the auxiliary reboiler will provide the "extra" Q reb i.e.
if cd reb QQ  then HPC cd QQ  where Q HPC is the energy yield by the heat pump to the distillation column. The factor by which the heat pump contributes to the heat load of the reboiler is given as For a distillation column with vapour recompression, driving the compressor uses the most energy. Thus, the power consumption must be known so as to assess the feasibility of such a system. For a perfect gas, that is, a gas having a constant specific heat, Cp = Cp o , then the specific enthalpy rise between the compressor inlet and outlet is And if the change of state is isentropic, In reality, ideal gases do not exist and therefore improvements are made on equation (41).
Therefore, compression is polytropic and the isentropic index γ, is replaced by the polytropic index, n (see Enweremadu, 2007). The compressor polytropic efficiency pol  = 0.7 -0.8 is used.
Also, because a saturated vapour, especially at higher pressures, shows deviations from the ideal gas behaviour, the compressibility factor, Z is used. Hence equation (41) Equation (43) shows that the pressure ratio b a P P is crucial to the power requirement. This ratio or the pressure increase to be provided by the compressor of a column with vapour recompression is influenced by the following (Meili, 1990;Han et al, 2003):  Pressure drop in vapour ducts (pipes) and over valves and fittings, ∆P p .  Pressure drop across the column, ∆P cl .  The difference in boiling points between the top and bottom products, ∆P b .  Temperature difference in the reboiler, ∆P CHP .

Determination of the pressure increase over the compressor
Pressure drops in the vapour ducts may be caused by frictional loss, ∆P f ; static pressure difference, due to the density and elevation of the fluid, ∆P s ; and changes in the kinetic energy, ∆P k . Since, there are elbows, valves and other fittings along the pipes then the pressure drop is calculated with resistance coefficients specifically for the elements. Therefore, the pressure drop along a circular pipe with valves and fittings is given by and u is the fluid velocity; d p is the pipe diameter and ρ is the fluid density;  is the Fanning friction factor which is a function of Reynolds number; l p is the pipe length;  is the dynamic viscosity of the fluid and ξ is the resistance coefficient.
The pressure drop over the entire distillation column, ∆P cl is caused by losses due to vapour flowing through the connecting pipes and through pressure drop over the stages in rectifying and stripping section. This depends mainly on the column internals, number of stages, gas load and operating conditions. ΔP cl =0, if zero vapour boil up is assumed. But constant pressure drop is assumed in this work. The pressure drop over a stage consists of dry and wet pressure drop. The dry pressure is caused by vapour passing through the perforation of the sieve tray. The aerated liquid (static head) on the tray causes the wet pressure drop. Constant pressure drop per tray have been estimated from several authors to be equal to 5.3mmHg per tray (Muhrer, Collura and Luyben 1990 The temperature difference in the reboiler-condenser is expressed by means of the vapourpressure equation as a pressure difference, ∆P CHP . Temperature differences of 8 -17 o C are quite common for ethanol-water distillation (Gopichand et al, 1988;Canales and Marquez, 1992 where P TOP is inlet pressure (vapour pressure at top temperature).
Other compressor parameters are calculated by the following equations: i. Compressor power input is determined from equation (43) and (51) where Cv is empirical volumetric coefficient and r is the compressor clearance ratio.
iv. Compressor nominal capacity or compressor displacement rate where V c is compressor displacement volume (m 3 ) and ω is angular velocity (rad s -1 ).

Determination of the reboiler-condenser parameters
The overall heat transfer coefficient between condenser and reboiler is given by However, a careful analysis reveals that the overall heat transfer coefficient U is an explicit function of Prandtl, Reynolds and Nusselt numbers, and depends on other properties such as viscosity and thermal conductivity. The overall heat transfer coefficient referenced to inner surface is given by As thermal resistance of the wall is negligible, (K wall is large and ln(r o /r i )) ≈ 0, it is then compared with the inner tube diameter (r i /r o ≈1) where μ m is the mean bulk fluid viscosity and μ w all is the viscosity of the liquid at the wall.
The expression for condensation at low velocities inside tubes is adapted from ( From the condenser prescribed degree of sub-cooling, the temperature of the working fluid after cooling and before throttling is given by where T SC is the degree of sub-cooling (K) The corresponding latent heat (enthalpy) is given as

Analysis of distribution of excess heat rate
The distillation system uses the column's working fluid as refrigerant and does not execute a closed cycle. Therefore the excess heat which may occur is not assessed by an overall energy balance but by the method of Oliveira et al (2001). When the energy available at the condenser Q cd , is greater than the energy required by the reboiler Q reb , the column receives the amount Q cd and the energy left over corresponds to the excess. But if Q cd is smaller than or equal to Q reb , then all the energy available is transferred to the reboiler, i.e. there will be no excess. where Q 23 is the excess heat due to energy interactions between the heat pump and the reboiler.
The distribution of the excess heat rate, Q 23 , between the pre-heater (Q 2 ) and cooler (Q 3 ) is accomplished by controlling the feed condition pre-heated by Q 2 . In other words, the value of Q 2 should be such that the feed reaches a prescribed condition. The pre-heating of the feed is carried out by Q 1 (heat exchanged between the bottom product and the feed) and Q 2 (heat exchange between the heat pump working fluid and the feed), in the heat exchangers. The heat provided by the bottom product is determined as follows: where, T CEV is the column evaporation temperature. Cp B is the specific heat of the bottom product. T BE is the temperature at the bottom product flow exit.
The best feed condition is that of saturated liquid (Halvorsen, 2001). The energy required to pre-heat the feed to reach saturated condition is expressed as , .
( ) where, Tsat, F and T F are the saturation temperature of the feed source and the temperature of the feed source respectively. Cp F is the specific heat of the feed source.
To make the feed a saturated vapour, the energy required is given as It is important to verify whether Q 1 alone is capable of pre-heating the feed to reach the desired condition, otherwise the amount of heat that should be withdrawn from the second pre-heater Q 2 , will be determined as , , .
The value of heat at the second pre-heater Q 2 should be, at the most equal to Q withdrawn to prevent the feed reaching 50% dry. Therefore, a convenient heat load control could be made as follows:

Thermodynamic analysis
Since vapour recompression uses a refrigeration cycle rather than a Carnot cycle, the performance of the heat pump is defined according to the following relation; The thermodynamic efficiency of a separation process is the ratio of the minimum amount of thermodynamic work required for separation to the minimum energy required for the separation (Olujic et al, 2003). For a vapour recompression distillation column, the energy required for separation process is composed of the reboiler heat load, Q reb , and the compressor power input, For the separation of a binary mixture by distillation the minimum thermodynamic energy required to achieve complete separation is given by (Liu and Quian, 2000): Then the thermodynamic efficiency is expressed as:

Solution method and error analysis
The equations that model the system components were grouped together in one single system. The analyses of the status of the variables were carried out to identify those that were the input data and those which were the unknowns. The equations were then grouped together, resulting in a set of non-linear algebraic equations, which were solved iteratively based on the step by step use of the successive substitution method. Solution was obtained when convergence was attained. The convergence was checked by using the criterion: The model is coded in MATLAB environment and used to evaluate the unknowns. A control programme for column VRC was written to compare the actual column (column in which the parameters studied were considered, VRC ΔP ). Figure 3 shows how the compressor power input,

Effects of pressure increase over the compressor
. cp W varies with the pressure increase over the compressor, ΔP. It is obvious from the plots that an increase in pressure over the compressor increases the compression (pressure) ratio leading to increase in compressor power input. The curve in Figure 3 shows the effect of pressure increase over the compressor on coefficient of performance. As the pressure increase over the compressor, ΔP increases, the compression ratio increases and the coefficient of performance, COP decreases due to increase in compressor power input. In Figure 4, a negative non-linear relationship exists between the compressor volumetric efficiency and pressure increase over the compressor. For a given compressor nominal capacity, when the pressure increase across the compressor, ΔP increases, the pressure ratio also increases resulting in the reduction in volumetric efficiency.
The influence of pressure increase over the compressor, ΔP on the required compressor displacement rate (Vcω) is shown in Figure 4. The linear relationship that exists between the variables indicates that the compressor displacement rate is directly proportional to ΔP. A reduction in compressor volumetric efficiency caused by increase in ΔP, increases the compressor displacement rate. This implies an increase in the displacement volume required although this cannot be observed in the figures. The compressor displacement required for a given speed is related to compressor size. For the large specific volume obtained, 9.74m 3 kg -1 ethanol as the heat pump working fluid will require a compressor of greater capacity i.e large displacement rate.
It can be observed from Figure 5 that an increase in pressure increase across the compressor, ΔP increases the total energy consumption. For a given reboiler heat transfer rate, Q reb , it is obvious that as ΔP increases, the total energy consumption increases. The total energy consumption and hence the energy savings from the work of Oliveira et al (2001) in heat pump distillation gave lower values. This could be attributed to non-consideration of the effect of pressure increase over the compressor and the subsequent increase in compression ratio and in compression power input respectively.

Effects of column heat loss
The direct effect of column heat loss could be seen from equation (17). From the equation, it follows that, for a given reboiler heat load or heat expenditure, fewer trays are required for a given separation if heat losses are reduced. Where heat loss occurs, more vapour has to be produced in the reboiler, since the reboiler must provide not only the heat removed in the condenser but also the heat loss. The effect of this is a decrease in process and energy efficiency. Indirectly, heat loss affects the column size in terms of number of plates. In the control system, the reflux ratio was as high as 7.5 compared with 5.033 obtained for the actual system. Therefore, if heat losses are properly accounted for, there may not be any need for downward review of the number of plates in order to reduce the reflux ratio (Enweremadu and Rutto, 2010). Therefore, pressure drop across the column, ΔP cl and the difference in boiling points between the top and bottom products, ΔP b which have the most profound effects on the pressure increase across the compressor, ΔP will be properly predicted. The overall implication of this is that the column size would be determined properly.

Overall heat transfer coefficient
Analysis of the overall heat transfer coefficient, U of the heat pump reboiler-condenser revealed an increase in the value of U. This was expected as the value of U in boiling and condensation processes are high. Also, the value of the heat transfer coefficient of the condensing ethanol is dominated the relationship used in determining U. A low value of overall heat transfer coefficient U will result in an increase in the heat exchanger surface area which may be a disadvantage to ethanol-water system. But the results from this work showed an increase in the value of U with the implication of a reduction in the reboilercondenser heat transfer area.
The variation of the reboiler-condenser thermal conductance (UA) with the heat pump distribution factor, f, is shown in Figure 6. The plots show that the greater the heat load taken by the heat pump i.e. larger f's, the larger the thermal conductance, UA, and the larger the heat exchanger area. However, for better performance of any heat transfer system, the thermal resistance (R th ) which is the inverse of thermal conductance (R th = 1/UA), should be as low as possible. Therefore the value of the heat transfer area for the VRC ΔP system will be smaller compared to the VRC system. Hence, the reboiler-condenser studied has a better performance.
The relationship between the thermal conductance, UA and the reboiler-condenser temperature difference, ΔT CHP shows that the higher the reboiler-condenser temperature difference, the lower the thermal conductance. The implication of this is that a higher ΔT CHP causes a reduction of the necessary heat transfer area. However, beyond a certain limit of the thermal driving force, the heat transfer area and the performance of the heat pump reboiler-condenser reduces. This is expected as higher ΔT CHP leads to higher compression ratio, higher compressor power input and higher energy consumption.  Table 1 summarises the comparison of some parameters of the two vapour recompression columns studied, the control (VRC) and the actual (VRC ΔP ). Both systems operate at 101.2 kPa and the effect of pressure drop effect is considered to account properly for variations in heat duty. A pressure drop of 0.707kPa per tray is assumed here as a reasonable estimate for the purposes of this study. The pressure ratio as used throughout this work is the ratio of the condensation pressure, P CHP to the top pressure, P TOP , for the VRC system and sum of the condensation pressure, P CHP and pressure increase over the compressor, ΔP, to the top pressure, P TOP for the VRC ΔP system. Since the energy consumption changes linearly with the feed flow rate, and as the present work makes a comparison of the relative performances, the base case flow rate was taken to be 1.098x10 -4 kmol s -1 . Also the column heat loss and the effect of such parameters as pressure increase across the compressor, the overall heat transfer coefficient of reboiler-condenser as an explicit function of Prandtl, Reynolds and Nusselt numbers which in turn depend on fluid properties are considered to account properly for variations in heat duty (Enweremadu, 2007;Enweremadu et al, 2009). Table 1 shows that the VRC enables some energy savings when compared with VRC ΔP . The VRC has a slightly lower compression ratio and consumes less energy than VRC ΔP system. Although there was a marginal increase in ΔP, which increased the compressor power input slightly, the performances of the two systems differ greatly by 27.5%. However, in addition to the pressure increase, the energy consumption appears to depend more on the rate of heat transfer in the reboiler. Hence the total energy consumption is indirectly related to column heat loss and pressure increase across the compressor through reboiler heat transfer and compressor power input respectively.

Comparison of the vapour recompression distillation systems
Also from Table 1, the heat transfer duties of the reboiler-condenser in terms of the overall heat transfer coefficient for the two systems show that the VRC ΔP has a higher value when compared with the VRC system. This implies that the VRC ΔP will require a smaller heat  (Enweremadu, 2007;Enweremadu et al, 2008Enweremadu et al, & 2009 transfer area which is economical in terms of material conservation. Also, with higher U, the VRC ΔP will have a better performance. The column heat losses for the two vapour recompression distillation columns are shown in Table 1. The heat losses in distillation columns with heat pumps have been assumed to correspond to around 3% of the energy supplied to the reboiler (Danziger, 1979). Oliveira, Marques and Parise (2002) assumed it to be as high as 10%. However, in this study, the heat exchanged by the distillation column with the surroundings is considered and its effect included in the balance equation (1). The results obtained for the VRC ΔP showed a marked difference between the two systems.
A comparison of the main reboiler heat transfer rate for the two systems is presented in Table 1. It is evident that neglecting and /or assuming a value for column heat loss instead of determining it had a significant effect on the values of Q reb in both systems. The heat pump distribution factor, f, for the VRC ΔP system is 0.346 which is slightly less than that for the VRC system (0.451). Since low value of heat load taken by the heat pump, f implies lower thermal conductance, then the value of the heat transfer area for the VRC ΔP system is smaller when compared to the VRC system. However, lower value of thermal conductance, UA, for the VRC ΔP system indicates that the VRC system will have a better performance.
The simulation results also show that the coefficient of performance and the thermodynamic efficiency of the VRC ΔP system is lower when compared to the VRC system. Fonyo and Benko (1998) have shown that the electrically-driven compression heat pump should work with a COP not lower than 3-5. The value of 5.85 obtained from this study has shown that although there is a decrease in the effectiveness of the VRC ΔP system when compared with the VRC system with COP of 6.15, it is within the acceptable range. This may be due to the fact that the compressor power input and the total energy consumption in the VRC ΔP system is higher than in the VRC system.

Conclusions
From the outcome of the study, the following conclusions may be drawn: 1. Pressure increase across the compressor, ΔP increases the compression ratio, the compressor power input, temperature in the heat pump reboiler-condenser while the compressor volumetric efficiency decreases. The effect of these is the reduction in the heat pump coefficient of performance and the use of a compressor of greater capacity.