Mathematical modeling for extraction of essential oil from Aquilaria crassna by hydrodistillation and quality of agarwood oil

Agarwood oil is widely used in medicine, perfume and cosmetic industries due to its unique properties and high medicinal value. However, knowledge on distillation process and standard for quality of the oil are far from sufficient. In this work, agarwood oil was hydrodistilled at various temperatures in order to study the distillation kinetics and quality of the oil. Higher oil yield and better quality of the oil were obtained at higher operating temperature. The mechanism of oil hydrodistillation consisted of two stages: washing stage and diffusion stage. The diffusion stage was the controlling step of oil hydrodistillation. Both mass transfer kinetic coefficients increased with increasing distillation temperatures. The oil gain at a higher temperature had a darker color and seemed to stick longer to the human skin. The different physical properties of the oil were a result of different chemical constituents exhibited in GC-MS profiles. Article Info Received: 29 June 2011 Accepted: 30 June 2011 Available Online: 3 July 2011

improve distillation process without the need to run too many experiments (Cassel et al., 2009). By modeling the kinetics of distillation, one can contribute not only to the fundamental understanding of the process but also to its better control and higher efficiency. Cassel et al. (2009) studied steam distillation modelling for essential oil extraction from rosemary (Rosmarinus officinalis L.), basil (Ocimum basilicum L.) and lavender (Lavandula dentate L.). Their kinetic model fitted very well with their experimental data and the model parameters could be useful during the scale-up of the extraction process. Koul et al. (2004) developed a simple first order kinetic model for steam distillation of lemongrass. The behaviour of oil production inside the distillation unit is not instantaneous, but it requires time for grass wetting, diffusion and osmosis of oil.
The models for steam distillation and hydrodistillation consider essential oil as one pseudo-component and most of them assume that diffusion of the solute to the particle surface is the controlling step. Milojevič et al. (2008) investigated the hydrodistillation kinetics of essential oil from comminuted ripe juniper berries (Juniperus communis L.) using a two-parameter model of unsteady-state diffusion through the plant material. The model was well correlated with their experimental results.
One of the major limiting factors of increasing trade and market penetration of agarwood oil has been a lack of quality standards. Hidayat et al. (2010) reported that agarwood oil grading has been traditionally performed by trained human graders (sensory panels). The disadvantages of this technique include subjectivity, poor reproducibility, time consumption and large labor expense (Keller, 1999). A human nose cannot tolerate a high number of samples and work in continuous production because it fatigues rapidly with increasing number of samples (Hidayat et al., 2010;Keller, 1999). Identification of agarwood oil would seem to be impossible without the use of laboratory techniques (Barden et al., 2000). However, in the trading market agarwood oil is graded according to the basic specification, such as the oleoresin content, the color and the aroma/smell. The objective of this study is to investigate the effect of distillation temperature on extraction of agarwood oil from A. crassna in terms of oil yield and oil quality (physical properties and chemical composition). Then a developed mass transfer mathematical model is validated with the experimental results to describe the process behavior.

Materials and Methods
Plant material: The wood of A. crassna was obtained from Kritsana Cooperative Club, Trat, Thailand. The plants were grown for 7 years before resin activated for 2 years and then harvested. The wood was dried until the moisture content was 3% (dry basis) before being milled and sieved to a size of 0.7 cm.
Extraction of essential oil: The comminuted wood (3,000 g) was soaked in distilled water (27,000 g) for up to 7 d (Liu et al., 2008;Bhuiyan et al., 2009) in order to break down the parenchymatous and oil glands. This ratio of agarwood to distilled water (1:9) is usually applied in the industrial production of agarwood oil. The soaked agarwood was then placed with water into a distillation apparatus. Hydrodistillation was operated at various temperatures: 80, 100 and 120°C. The distillation temperature was not higher than 120°C to prevent burning of plant material. During distillation the vapor, composed of steam and volatile compounds, rises from the extractor to the condenser, where the two immiscible liquid phases, namely aromatic water and agarwood oil, are formed. In the separator, the lighter agarwood oil is separated from the aromatic water, which returns to the extractor via a recycle pipe. Agarwood oil was collected every 8 h and left in the separator funnel for essential oil separation. The distillation was conducted in triplicates and the average values are presented.
The distillation unit consists of four major components ( Figure 1): a stainless steel extractor with an inner diameter of 0.4 m and a height of 0.5 m, a stainless steel condenser with an inner diameter of 0.4 m and a height of 1.2 m, a stainless steel separator with an inner diameter of 0.3 m and a height of 1.2 m, and a cooling tower driven by a 1.5kW motor. The extractor was heated by using a gas stove which was connected to a controller with an accuracy of ±0.1 °C to control the gas usage and extraction temperature.
Study of distillation kinetic model: The mass transport of essential oil through plant particles during hydrodistillation occurs as unsteady-state diffusion. The diffusion of the solute to the particle surface is assumed to be the controlling step. For batch hydrodistillation and no chemical reaction, Fick's second law (Eq. 1) of one-dimensional unsteady-state diffusion can be applied. (1) Where qp is the average concentration of essential oil in the plant particles at time t, x is the distance along the direction of diffusion, and Deff is the effective diffusion coefficient. The following assumptions at the level of plant particles in the suspension should be made to solve Eq. 1: 2) the essential oil is considered as a pseudocomponent; 3) a part of essential oil is located at the external surfaces of the broken agarwood and the rest is uniformly distributed in the plant particles; 4) the effective diffusion coefficient is constant; 5) the concentration of essential oil on the external surfaces of the plant particles at any moment during the distillation is zero because of its instantaneous "washing" from the surface; and 6) the resistance to the mass transport of essential oil from the external surface of the plant particles is assumed to be zero.
Integrating Eq. (1) with the above suitable initial and boundary conditions (Treybal, 1975) leads to Eq. (2): ( 2) where q is the average yield of agarwood oil at time t, qo is the initial average concentration of agarwood oil (found by extraction with hexane in a Soxhlet apparatus), A is a constant, and k is a kinetic constant including the effective diffusion coefficient. Based on this model, the extraction of agarwood oil by hydrodistillation consists of two successive stages, a washing stage and diffusion stage, as proposed by Milojevič et al. (2008): In washing stage, the essential oil is washed from and near the external surface of the plant particles during an initial, short period of distillation (theoretically at t = 0). It can be characterized by a rapid increase in the oil yield at the very beginning of the process (Eq. 3): (3) Where, at t = 0; qw is the essential oil yield and b is the relative amount of essential oil extracted.
In diffusion stage, the essential oil diffuses from the interior of the plant particles towards their external surfaces followed by its distillation. It is characterized by a slow exponential increase in the oil yield with the progress of distillation. The basic kinetic equation for agarwood oil hydrodistillation was derived by combining Eqs. (2)  (5) Study of physical properties: Relative density and refractive index were analyzed using two standard procedures: ISO 279, a method for determining the relative density of essential oils (ISO, 1998); and ISO 280, a method for determining the refractive index of essential oils (ISO, 1998). Surface tension was determined by a tensiometer (Model OS, Torsion, UK).

Study of chemical compounds by GC-MS:
Agarwood oil composition was studied by GC-MS analysis using a Hewlett Packard gas chromatograph (GC 5890) coupled with a mass selective detector (5972) (Hewlett Packard, Palo Alto, USA). Separation of the analytes by gas chromatography was carried out using a silica capillary column (30 m length, 0.25 mm diameter, 0.25 mm film thickness) of HP-5MS (Hewlett Packard). Separation of the compounds involved injection of 1.0 mL of the oil into the front inlet of the gas chromatograph operating at 250°C in the splitless mode. The flow rate of the carrier gas, helium, was 2.0 mL/min with a 1:50 split ratio. The oven program commenced at 80°C, where it was held for 2 min and then increased at a rate of 10°C/ min to 250°C, where it was held for 10 min. The interface temperature was 250°C. Ionization of the analytes by electron impact (EI) was obtained using an emission current of 70 eV. The ion source temperature was set at 250°C and the scan scope was set from 32 to 500 amu. The compounds were characterized by database matching and comparison of their MS spectra with existing data in the Wiley and Adams library search data (Adam, 2001).

Results and Discussion
From soxhlet extraction the initial average concentration of agarwood oil was 1000 mg/g. With hydrodistillation for 64 hours the average agarwood oil yields at temperatures of 120, 100 and 80 °C were 611.67, 515.33 and 455.33 mg/g dried wood, respectively. The oil yields at various distillation temperatures as a function of distillation time are shown in Figure 2. It was notable that the curves are typical for the hydrodistillation of essential oils from aromatic plant materials and verify the proposed mechanism of oil isolation (Milojevič et al., 2008). The mechanism consisted of two stages: fast oil distillation and slow oil distillation. The oil yield increased rapidly during the initial stage because the oil was rapidly washed from the external surfaces of the wood during short period of hydrodistillation. This is called fast oil distillation stage. In the latter stage (slow oil distillation), the extraction rate gradually slowed down until nearly constant as a result of diffusion of oil from interior of the wood towards it surface. Figure 3. It can be seen that Eq. (5) described the experimental data quite well. This result confirms that Figure 2: Yield of essential oil at various distillation temperatures Figure 3: Kinetics of essential oil from experimental data fitted to the mass transfer kinetic model at various distillation temperatures the hydrodistillation of agarwood oil was controlled by diffusion stage. The operation with a small particle size of wood and long soaking time might help to promote the distillation efficiency. The parameters of the kinetic model, b and k, were calculated from the experimental data by means of linear regression method using Eq. (5), as shown in Table I.

The dependence of ln [(q0-q)/q0] versus time is shown in
The relationship of washing coefficient (b) and diffusion coefficient (k) on the distillation temperature can be expressed as shown in Eq. (6) and Eq. (7), respectively.
where T is the distillation temperature (K); the coefficients of linear correlation of k and b were 0.92 and 1.00, respectively. Both kinetic parameters increased with increasing distillation temperatures. This result implies that the oil was washed and diffused from the wood faster at higher temperatures as the diffusivities and driving forces of solute (oil) and solvent increased (Mani et al., 2007).
Physical properties of agarwood oils distilled at various temperatures are shown in Table II. It is noticeable that higher values of refractive index, relative density and surface tension were observed at higher operating temperatures. The oil with higher refractive index presents darker color. The oil with greater relative density is assumed to contain higher amount of aroma resinous compounds that distinguishes agarwood oil from the others. The oil with greater surface tension shows higher strength and longevity of the oil fragrance. With higher surface tension, the oil seems to stick longer to human skin. With these physical 023 .     properties along with the smell agarwood oil can be simply graded. In addition, the agarwood oil quality can be basically tested according to some assumptions: pure agarwood oil will not give different color at the base and surface inside the bottle; when put the oil in front of the high power lamp the shining light from the bulb cannot see trough; the oil should remain on hand more than 5 hours; and the very finest agarwood oil will sink in water.
From GC-MS analyses the agarwood oil distilled at 120, 100 and 80°C gave the major compounds of 21, 13 and 8, respectively. This indicates that distillation temperature had an effect on chemical compounds of extracted oil. At higher distillation temperature, the oil components with higher boiling point can be extracted as a result of a higher number of chemical compounds. As shown in Table III, agarospirol, alloaromadendrene, valencene and valerenol were the major compounds (chromatogram area > 6%) of the oil hydrodistilled at 80 °C. Aristolene, 1,5-Diphenyl-3-pentanone and agarospirol were found to be the major compounds of the oil distilled at 100 °C (Table IV). From Table V, aromadendrenepoxide, agarospirol, alpha-Gurjunene and 2-naphthaleneethanol were the major compounds of the oil distilled at 120 °C. It is noticed that agarospirol, hexadecanoic acid and 1,5-Diphenyl-3pentanone were found in all oil samples. However, hexadecanoic acid or palmitic acid and 1,5-Diphenyl-3pentanone are not the fragrant agents. Agarospirol was the compound identified in agarwood oil by many reports Bhuiyan et al., 2009). Therefore, agarospirol should be a marker compound of the agarwood oil from A. crassna. Agarospirol is used in the pharmaceutical industry for sedative (Takemoto et al., 2009;Okugawa et al., 1996). Aristolene and benzylaccetone demonstrate positive effects on the central nervous system (Okugawa et al., 1996). 10 epi gamma eudesmol and junipene present sweet woody odor and are used in the fragrant industry (Takemoto et al., 2009).
In conclusion, hydrodistillation of essential oil from A. crassna at 120°C gave maximum agarwood oil yield of 611.67 mg/g dried wood in 64 hours. The agarwood oils obtained at higher temperature have better physical properties; darker colour, heavier, higher strength of oil fragrance and longer adhesion to the human skin. Grading of agarwood oil could be performed by testing the physical properties as well as analysis of the chemical composition by GC-MS.