SciELO - Scientific Electronic Library Online

vol.38 número4Improved dynamic window approach by using Lyapunov stability criteriaEvaluation and selection of discrete-event simulation software for the oil industry índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados




  • No hay articulos citadosCitado por SciELO

Links relacionados

  • No hay articulos similaresSimilares en SciELO


Latin American applied research

versión impresa ISSN 0327-0793

Lat. Am. appl. res. v.38 n.4 Bahía Blanca oct. 2008


Simulation of deacidification process by molecular distillation of deodorizer distillate

M. A. Marttinello1, I. Leone1 and M. Pramparo2

1 Facultad de Ingeniería, UNRC, 5800 Río Cuarto, Córdoba, Argentina.
2 Facultad de Ingeniería, UNRC, 5800 Río Cuarto, Córdoba, Argentina.

Abstract — A computer program has been developed to simulate the deacidification of sunflower deodorizer distillate. The developed mathematical model is based on momentum, mass and energy balances. The Langmuir-Knudsen constitutive equation was used in order to represent the kinetics of evaporation and condensation. Physical and transport properties have been evaluated using correlations, which are functions of temperature and composition. The obtained partial differential equation system was solved by means of the Crank-Nicholson technique, and the calculation programs employed were developed and debugged using Matlab 7.1.
The mathematical model was used to analyze the phenomena that take place during the molecular distillation process, and it also allowed studying the influence of operating variables on the performance of the process. In this investigation, we analyzed the influence of the operating temperature on purity and yield of the components. Trials at three evaporating temperatures (110 ºC, 130 ºC, and 140 ºC) were carried out in a KDL4 (UIC) molecular distillation apparatus to verify the model.

Keywords — Molecular Distillation. Simulation. Free Fatty Acids. Tocopherols. Sterols.


Molecular distillation operation, which is also known as short-path distillation, is a useful separation technique in the purification of compounds which are thermally unstable and have low volatility. Molecular distillation is based on the evaporation of compounds in a mixture. These phenomena take place in two types of distillers: falling film or centrifugal. In the falling film method, one film is in contact with a heated surface and the other film is in contact with a cold surface. The cold surface is near the heated surface, at a distance shorter than the mean free path of the residue gas. The main operating characteristic is the working pressure, and it ranges from 10-2 to 10-4 KPa. Under these conditions, the volatility of the components increases and the operation temperature decreases, allowing the compounds to separate at lower temperatures. The molecules that leave the evaporation surface suffer virtually no collisions at lower operating pressures. As a result, the molecules take a short path before their condensation; that is, they arrive at the condenser surface in a brief period of time. It is for this reason that very high evaporation velocities are reached. Under these two conditions, short residence time and low temperature, the thermal decomposition of the components is reduced and separation occurs at acceptable velocities (Weissberger, 1951; Perry et al., 1984).

The separation level that is obtained in a molecular distillation process is not only a function of the volatility of the components, but it is also a function of the mass and heat transfers in the liquid phase and its molecular kinetics play an important role in the performance of the operation. When the liquid is evaporated, the vapor-liquid interface becomes cold and, in mixtures, it decreases the composition of the more volatiles components. This leads to the formation of driving forces for the diffusive transfer of mass and heat. All these resistance affect not only the evaporation velocity, but also the product purity (Bose and Palmer, 1984).

Advances in the theoretical modelling of this operation have been reported. Bose and Palmer (1984) showed that the existence of concentration and temperature gradients in the liquid phase decreases the separation efficiency. Bandarkar and Ferron (1991) worked in the mass and heating transference in a liquid film on a centrifugal still. Batistella and Maciel (1996) have shown results of efficiency and performance with two types of distillers. The pressure effect in the molecular distillation process has been analyzed by Cvengros and Lutisan (1995) using both Boltzmann proposals and the Maxwell distribution function for the molecular velocity in the vapor phase. Cvengros et al. (1997) used a model to theoretically analyze the operation of molecular distillation for binary mixtures.

Tocopherols and sterols are compounds with important economical value. They decompose at the high temperatures that are required in the conventional distillation operation (Yankah and Jones, 2001; Dunford, 2001). These compounds can be obtained from a by-product of the deodorization stage of the vegetable oil refining process: the 'deodorizer distillate'. Besides sterols and tocopherols, the distillate generally contains lighter compounds such as free fatty acids, and heavier components such as mono, di and triglycerides, and esters from sterols and pigments (Winters, 1994). This mixture must be treated with separation techniques in order to obtain the purified compounds. The separation of free fatty acids from the remaining components of the mixture is generally the first stage of distillation in the separation and purification process, because the fatty acids are more volatile than the rest of the components. However, the conventional distillation requires high temperatures, which tocopherol decomposition. In this context, the molecular distillation method becomes a promising process.

Since these processes are carried out on multi-components mixtures, the application of mathematical models for the previous study of the operation is very complicated since it requires the use of many component properties. Consequently, the study of a simplification in the treatment is justified. The aim of this work is to analyze the deacidification process of sunflower deodorizer distillates through molecular distillation using a mathematical model and its respective experimental verification considering the multi-components mixture as a set of two identities.


A. Mathematical Modelling

Figure 1: Molecular distillation equipment scheme.

In Fig. 1 the performance of a falling film molecular distiller is shown in a schematic way. The transport equations, which describe the evaporation process in the film, are obtained using momentum, mass and energy balances.

To develop the mathematical model, the following assumptions have been made:

1. Condition of steady state throughout the process.
2. Temperature of condensation set low to avoid re-evaporation from the condensation surface.
3. A rectangular coordinates system to describe variables over the film, because thickness film is very small in comparison with the radius of the cylindrical evaporation surface.
4. Negligible collisions between molecules in the vapor phase.

A binary mixture (A+B) is considered. The main free fatty acids are considered as entity A (lighter compounds of the mixture), and the remaining less volatile components (the heavier ones) as the other entity (B). Their properties have been calculated considering the relative quantity of each component, in order to have averaged properties of each entity.

For movement description in the system considering laminar flow, the momentum balance is formulated as follows.


with the following boundary conditions:


From the continuity equation for component A of a binary mixture (A+B), the following is obtained:


with the following conditions:


where xAf is the molar fraction of A in the feed.

For boundary condition in z=h, Fick's first law for diffusion is used:


where n represents normal direction to the evaporation surface and s represents the liquid surface.

The molar flow density for the components is expressed as a function of the evaporation rate:


where xAs is the molar fraction of component A in the evaporation surface.

The evaporation rate of each component is evaluated by Langmuir, which is based on the kinetic gas theory:


where Ts is the evaporation surface temperature.

As a result of formulating a global mass balance, a variation equation of molar total flow is obtained:


where R1 is the average radio of the film. The boundary condition related to this equation is formulated as:


where I f is the feeding molar flow.

The energy balance equation for this system has the following form:


The boundary conditions for this equation are:


where Tf and Tw are feed and wall temperatures, respectively.

For boundary condition in z=h, the equation has the following shape:


The thickness value of the film, h, is related with molar flow I by means of:


Figure 2: Scheme of the net employed to solve the numerical method.

In Fig. 2, the scheme of the net employed to solve the numerical method is shown. The partial differential equations were formulated in terms of the implicit finite difference method, considering the curved borders present in the falling films. The set of linear equations from finite difference scheme is solved by the Gauss-Seidel method.

The α value is determined by:


The equations used to formulate the condensation mathematical model are the conservation equations of mass, movement and energy already formulated for the evaporation film with the following changes:

1- The film initial velocity is zero.


2- The initial temperature is equal to the temperature of the cooling fluid:


3- Initial compositions of A and B can be determined from:


The complete mathematical model has been numerically solved using MatLab 7.1. In order to determinate the physical properties, correlations available in bibliographical data were used. An average molecular weight was used taking into account the relative mass from the principal fatty acids in the deodorizer distillate (Madhavi et al., 1996) and to represent the less volatile compounds, an average molecular weight was calculated taking into account the compositions of the more abundant ones (Lide, 1992). The average vapor pressure has been considered as a function of the temperature (Madhavi et al., 1996). The compound's critical properties were calculated with Ambrose's proposal (functional group contribution) (Reid and Sherwood, 1968). The transport properties thermal conductivity, viscosity and diffusivity have been obtained by means of the Sheffy-Johnson, Thomas and Sitaraman-Ibrahim-Kuloor correlations, respectively (Reid et al., 1987). Considering the thermodynamic properties, the heat capacity was obtained by means of the Chow-Bright correlation (Perry, 1984), the density by means of the Benson-Fishtine (Madhavi et al., 1996) equation and the latent heat of evaporation by means of Vetere equation (Reid and Sherwood, 1968). All properties were averaged pondering with mass fractions. The unknown compounds were not considered for the determination of the properties average; instead, they were assigned the calculated average properties.

B. Experimentals

Analytical methods. With the aim of knowing the characteristics of raw material and the quality of the obtained product, the following analytical determinations were used, by means of AOCS methods: Free fatty acids (AOCS Cd 3d-63), FFA composition (AOCS Ce 1-62-gas chromatography method) and Sterols, Tocopherols and Squalene (AOCS Ce 7-87- gas chromatography method) (AOCS, 1994).

Table 1. Composition of sunflower oil deodorizer distillate. (Deodorizer conditions: 3-4 torr, 240°C, 45 min)

Table 2. Free fatty acid composition for sunflower oil deodorizer distillate.

Raw material. A sample of deodorizer distillate was obtained from the sunflower oil refining process. The sample was liquid at ambient temperature and had a dark brown color, almost black. Table 1 shows the composition of the sample. Fatty acid distribution is shown in Table 2.

Figure 3: Scheme of molecular distiller

Molecular distillation. Molecular distillation was carried out in KDL4 (UIC) equipment (falling film, evaporation surface: 0.04 m2, condensation surface: 0.02 m2). The distiller is equipped with wiper rollers of variable velocity. The vacuum system consists of a diffusion pump and a mechanical pump. These are capable of obtaining a maximum vacuum of 10-4 KPa, with a maximum feeding of 0.5 Kg/h. A scheme of this apparatus is shown in Fig. 3.

Table 3. Operating conditions for the process

The experimental tests were planned in order to analyze the influence of the evaporator temperature. The feeding and condensation temperatures, the feeding flow rate, the pressure and the velocity of wiper rollers were established according to those typically used in these kinds of applications. Operative conditions used in the trials are shown in Table 3.

C. Results and Discussion

Table 4. Results of the molecular distillation process at different evaporator temperatures (from 100ºC to 140ºC). Mass balance and compositions.
* F: Feed, D: distillate, R: residue.

Table 5. Set of physical, thermodynamic and transport parameters of A and B at 120ºC.

Table 6. Vapor pressure for A and B.

Results of molecular distillation process at different evaporator temperatures from 100ºC to 140ºC are shown in Table 4. Mass balance and compositions of main components of these trials allow as comparing experimental results with theoretical results from the model developed. The physical parameters obtained in the application of the model for the sunflower deodorizer distillate used in the experimental trials at 120ºC can be observed in Table 5. Vapor pressure variation with the temperature is reported in Table 6, which shows the value of vapor pressure for A and B at two operating temperature.

Figure 4: Variation of relative quantity of heavier compounds in the residue (RB) with respect to the heavier compounds in the feeding (FB) for three evaporating temperature.

The change in the relative quantity of the heavier components in the residue at different non-dimensional vertical positions is shown in Fig. 4. This quantity shows the heavier remaining components in the residue at the vertical y-position. The results of three simulations for three different evaporation temperatures are shown. At higher evaporation temperature, the evaporated quantity of all components was higher, so the amount of residue was lower.

Figure 5: Variation of the total evaporated quantity with the evaporating temperature.

Results of the corroboration of the model are shown in Fig. 5 through 8. Figure 5 shows the variation of the total evaporated quantity, which is relative to the feeding, with the evaporation temperature. Based on the simulation results, the total evaporated quantity is approximately 35 % at 107 °C. Trials at higher temperatures were later carried out in order to evaporate, at least, the desired quantity of distillate (similar to the quantity of the lighter compounds that was intended to separate in the distillate).

Figure 6: Final compositions in residue (XA) and distillate (YA) at different evaporation temperatures.

The results of the final composition for residue and distillate at different temperatures are shown in Fig. 6. At higher temperatures, lighter compounds composition in the residue and in the distillate is lower. The more adequate temperature will depend on the particular objective of purification (to maximize the profitability of the process).

Figure 7: Mass flows of the lighter compounds (A) in the distillate and mass flows of heavier compounds (B) in the residue for different evaporating temperatures.

Figure 7 shows modeling and experimental results of mass flows of lighter compounds in the distillate and heavier compounds in the residue for different evaporation temperatures, relative to the respective feeding. Both flows should be as high as possible in order to get larger separation efficiency, but it does not occur simultaneously. There will be a technical trade-off solution.

The separation factor, Fs, is calculated as:


Figure 8: Separation factor as a function of the evaporation temperature.

Fs is shown in Fig. 8, in relation to the temperature. This factor suffers a larger decrease due to the resistance to the heat and mass transfer in the liquid. At higher temperatures, the factor decreases abruptly, so it is convenient to set this temperature in a value near to 110°C. This value allows for the evaporation of a distillate percentage to equal the percentage of the lighter compounds in feeding.


A computer mathematical modeling system to simulate the molecular distillation process was developed. It was used to study the deacidification operation of sunflower deodorizer distillate. The physical model allowed us to analyze the influence of evaporation surface temperature.

The model was experimentally verified, obtaining good accordance with the simulation model results and the data coming from the trials.

This methodology of analysis can be carried out to investigate improvements in the separation efficiency with other influential parameters such as flow, condensation temperature and pressure.


c = molar density, [mol/m3]
h = film thickness, [m]
ki = rate of component evaporation, [mol/m2,s]
k = thermal conductivity, [Joule/m,ºC]
L = evaporator high, [m]
nz = intervals on z direction
v = velocity, [m/s]
xA= molar fraction of A component
y = vertical coordinate, [m]
z = horizontal coordinate, [m]
cv = heat capacity, [Joule/mol,ºC]
DAB = diffusion coefficient, [m2/s]
I = molar flow, [mol/s]
M = molecular weight
N = density of molar flow, [mol/m2,s]
Pv = vapor pressure, [Pa]
R = universal constant of gases
T = temperature, [ºC]
ρ = density, [kg/m3]
µ = viscosity, [kg/m,s]
ν = cinematic viscosity, [m2/s]
ΔHV = latent heat of vaporization, [Joule/mol]
Δz = horizontal increment, [m]

1. AOCS, Officials Methods & Recommended Practices of the American Oil Chemists'Society, 4th edition. AOCS Press, Champaign, IL (1994).         [ Links ]
2. Batistella, C. and M. Maciel, "Modeling, Simulation and Analysis of Molecular Distillators: Centrifugal and Falling Film", Proc. of European Symposium, Computers Aided Process Engineering, 20, 19-24 (1996).         [ Links ]
3. Bhandarkar, M. and J.R. Ferron, "Simulation of Rarefied Vapor Flows", Ind. Eng. Chem. Res., 30, 998-1007 (1991).         [ Links ]
4. Bose, A. and H.J. Palmer, "Influence of Heat and Mass Transfer Resistances on the Separation Efficiency in Molecular Distillations", Ind. Eng. Chem. Fundam., 23, 459-465 (1984).         [ Links ]
5. Cvengros, J. and J. Lutisan, "Effect of Inert Gas Pressure on the Molecular Distillation Process", Separ. Sci. Technol., 30, 3375-3389 (1995).         [ Links ]
6. Cvengros, J., J. Lutisan and M. Micov, "Balance Equations for Molecular Distillation", Separ. Sci. Technol., 32, 3051-3066 (1997).         [ Links ]
7. Dunford, N.T., "Health Benefits and processing of Lipid-Based Nutritionals" Food Technology, 55, 38-44 (2001).         [ Links ]
8. Lide, D., Handbook of Chemistry and Physics, 72nd edition, CRC Press, Boston, USA (1992).         [ Links ]
9. Madhavi, D., S. Deshpande and K. Salunkhe, Food Antioxidants, Technological, Toxicol. Health Perspect., M. Dekker (1996).         [ Links ]
10. Perry, R., D. Green and J. Maloney, Perry´s Chemical Engineers´ Handbook, Sixth Edition, McGraw Hill (1984).         [ Links ]
11. Reid, R. and T. Sherwood, Propiedades de Gases y Líquidos, 2da edición, THEA, México (1968).         [ Links ]
12. Reid, R., J. Prausnitz and B. Poling, The Properties of Gases & Liquid, McGraw-Hill, Inc. 4th edition (1987).         [ Links ]
13. Weissberger, A. "Destillation under High Vacuum in Technique of Organic Chemistry", Destillation. Interscience Publishers, IV, 495-602 (1951).         [ Links ]
14. Winters, R.L. "Valor y Utilización de los destilados del Proceso de Desodorización en el Procesamiento Aceites Comestibles", Aceites y Grasas, 4, 54-61 (1994)-         [ Links ]
15. Yankah, V.V. and P.J. Jones, Phytosterols and Health Implications: Efficacy and Nutritional Aspects, Inform, 12, 899-903 (2001).
        [ Links ]

Received: August 2, 2007.
Accepted: December 17, 2007.
Recommended by Subject Editor: Orlando Alfano

Creative Commons License Todo el contenido de esta revista, excepto dónde está identificado, está bajo una Licencia Creative Commons