Biolubricant Production through Double Transesteriﬁcation: Reactor Design for the Implementation of a Bioreﬁnery Based on Rapeseed

: The production and use of biolubricants as replacements for mineral lubricants align with the promotion of sustainable development goals, contributing to the sustainable economic growth of developing countries as well as the preservation of the environment. The implementation of bioreﬁneries (where the production of biolubricants through transesteriﬁcation could play an important role) is becoming important for these purposes, using natural feedstocks such as rapeseed, which is an interesting crop that can adapt to different climates under extreme weather conditions. The aim of this research work was to study the design of a reactor to produce rapeseed biolubricant through double transesteriﬁcation of the corresponding vegetable oil. Thus, the kinetics to design a reactor was studied, sizing it according to the properties of the biolubricant and the demand in Spain. As a result, a SAE 10W30 biolubricant was obtained, which was suitable for Diesel engines. A batch reactor was selected for the production of this biolubricant at a national level, with a pseudo-ﬁrst reaction order and a reactor volume of 9.66 m 3 .


Introduction
A biolubricant is a lubricant that is neither detrimental to health nor harmful to the environment. Their main advantages, compared to the use of conventional lubricants, are their adaptation to green chemistry, biodegradability, sustainability and compliance with current legislation and customers' needs, focused on environmentally friendly products [1]. Apart from that, other physical properties that are improved compared to mineral lubricants are: better lubricity, higher flash and combustion points, or higher viscosity index, among others. However, there are some disadvantages such as their lower oxidative stability and cold flow properties [2]. Nevertheless, the use of additives like antioxidants can improve the performance of this product [3].
The main raw materials for biolubricant production are vegetable oils, such as rapeseed, sunflower, soybean, safflower, etc., which are mainly composed of triglycerides. In that sense, rapeseed is an interesting oil crop, as it is resistant to cold and hot climates in drylands, having short vegetation periods and a fatty acid composition that is suitable for food use, among others. Moreover, one interesting characteristic of this crop is its long tap root, which makes it suitable to take part in crop rotations [4]. Finally, its high oleic acid content (with one single double bond, which makes it less prone to oxidation compared to linoleic or linolenic acid) could contribute to the stability of the final products obtained through transesterification, such as biodiesel or biolubricants. Indeed, the influence of the raw material on these products has been widely studied in the literature, presenting different viscosity or oxidative stability values depending on the fatty acid profile [5,6].
There are several processes to obtain products (from vegetable oils) that can be used as biolubricants, such as selective hydrogenation, estolide formation, epoxidation or double transesterification [3,7]. Concerning epoxidation, consisting of the interconnection of unsaturated C-C bonds by an oxygen atom (at low reaction temperatures), the biolubricants obtained presented better lubricity and thermo-oxidative stability. However, pour point increased, whereas viscosity index was reduced. Regarding transesterification, it is the reaction of an ester and an alcohol to produce another ester through alkyl group exchange. Normally, the biolubricants obtained have good low temperature properties and high thermo-oxidative stability, whereas high reaction temperatures are required [3]. Figure 1 shows this process. As it can be seen, the first transesterification implies the use of triglycerides and methanol to produce fatty acid methyl esters (FAMEs, biodiesel) and glycerol. For the second transesterification, FAMEs react with a superior alcohol to produce the corresponding ester, releasing methanol. As a consequence, double transesterification could provide plenty of products and by-products that can be valuable during biolubricant production. The whole process shows a high atom economy, as most by-products can be re-used in this same process or valued for market purposes, which is a considerable advantage. In that manner, there are many positive aspects in this process that can be suitable for the implementation of biorefineries [8]: • Use of natural and renewable sources. • Production of valuable products (like biodiesel and biolubricant). • Generation of exploitable by-products (such as glycerol, with a wide range of uses depending on its purity [9]). • Reutilization of by-products (like methanol).
This way, biorefineries can contribute to the economic growth of many developing countries or less-industrialized areas. In that sense, the Extremadura region is a traditionally agricultural area in Spain, which can profit from circular economy and current environmental policies, where the installation of biorefineries can perfectly fit, as pointed out in previous studies [10].
However, the implementation of biorefineries might be a challenge, as many factors have to be taken into account, such as raw materials, technologies, processing routes, products, and technical, economic or environmental aspects [8]. This way, the design of different components for biorefineries can help assess the feasibility of the implementation of these facilities. As far as we know, few specific studies have been carried out in order to prove the feasibility of biorefinery implementation based on vegetable oils through the design of its components (such as reactors), determining the kinetics and production for a certain area or region depending on the consumption of the products obtained.
The aim of this study was to assess the role of rapeseed oil as a starting point for a biorefinery, focusing on biolubricant production through double transesterification. Thus, the kinetics of the second transesterification as well as the reactor design and sizing (according to a specific production of the final biolubricant) were carried out.

Biodiesel and Biolubricant Production
Rapeseed oil was produced in CICYTEX (Centro de Investigaciones Científicas y Tecnológicas de Extremadura) by the Department of Energy and Extensive Crops. The oil had high-quality parameters such as low free fatty acid content (<1%), only requiring filtration after mechanical extraction. Afterwards, the oil was stored in 25-L opaque containers for further treatments and analysis.
Rapeseed biolubricant was obtained through double transesterification with different alcohols. Thus, the steps were the following (see Table 1 for the main chemical conditions):

•
First transesterification with methanol. For fatty acid methyl ester production (biodiesel), the chemical conditions were based on previous studies [11]. These experiments were also used for the kinetic determination and reactor design. Once the reaction took place, the sample was purified by vacuum filtration, separating the surplus alcohol.  1 The alcohol used for the first transesterification was methanol, whereas for the second transesterification 2-ethyl-2-hydroxymethyl-1,3-propanediol was used. 2 The catalyst used in both cases was sodium methylate.

Biodiesel and Biolubricant Characterization
Both biodiesel and biolubricant were characterized according to the UNE-EN 14,214 standard [12], as explained in previous studies [2,11,13]. Density was determined by using a pycnometer, correcting the value at 15 • C. Viscosity was measured by using Ostwald's viscometer at 40 and 100 • C, obtaining with these data the viscosity index. Flash and combustion points were obtained through the open-cup Cleveland method. Moisture was determined by using a Metrohm 870 KF Trinitro Plus (Metrohm, Herisau, Switzerland). Acid and iodine numbers were determined following the corresponding volumetric determination, according to the standard and cold filter plugging point (CFPP) was obtained (when necessary) for the samples obtained [13][14][15]. Regarding the yield of the reactions and fatty acid methyl ester profile of biodiesel, gas chromatography was used, following the conditions explained elsewhere [16]. Thus, the yield of the biolubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows:

•
Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 • C and viscosity index and comparing with traditional lubricants. • Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level. • Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. lubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. yield of the reactions and fatty acid methyl ester profile of biodiesel, gas chromatography was used, following the conditions explained elsewhere [16]. Thus, the yield of the biolubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. (CFPP) was obtained (when necessary) for the samples obtained [13][14][15]. Regarding the yield of the reactions and fatty acid methyl ester profile of biodiesel, gas chromatography was used, following the conditions explained elsewhere [16]. Thus, the yield of the biolubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. Switzerland). Acid and iodine numbers were determined following the corresponding volumetric determination, according to the standard and cold filter plugging point (CFPP) was obtained (when necessary) for the samples obtained [13][14][15]. Regarding the yield of the reactions and fatty acid methyl ester profile of biodiesel, gas chromatography was used, following the conditions explained elsewhere [16]. Thus, the yield of the biolubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. volumetric determination, according to the standard and cold filter plugging point (CFPP) was obtained (when necessary) for the samples obtained [13][14][15]. Regarding the yield of the reactions and fatty acid methyl ester profile of biodiesel, gas chromatography was used, following the conditions explained elsewhere [16]. Thus, the yield of the biolubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. . Acid and iodine numbers were determined following the corresponding volumetric determination, according to the standard and cold filter plugging point (CFPP) was obtained (when necessary) for the samples obtained [13][14][15]. Regarding the yield of the reactions and fatty acid methyl ester profile of biodiesel, gas chromatography was used, following the conditions explained elsewhere [16]. Thus, the yield of the biolubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. ture was determined by using a Metrohm 870 KF Trinitro Plus (Metrohm, Herisau, Switzerland). Acid and iodine numbers were determined following the corresponding volumetric determination, according to the standard and cold filter plugging point (CFPP) was obtained (when necessary) for the samples obtained [13][14][15]. Regarding the yield of the reactions and fatty acid methyl ester profile of biodiesel, gas chromatography was used, following the conditions explained elsewhere [16]. Thus, the yield of the biolubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. . Acid and iodine numbers were determined following the corresponding volumetric determination, according to the standard and cold filter plugging point (CFPP) was obtained (when necessary) for the samples obtained [13][14][15]. Regarding the yield of the reactions and fatty acid methyl ester profile of biodiesel, gas chromatography was used, following the conditions explained elsewhere [16]. Thus, the yield of the biolubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. . Acid and iodine numbers were determined following the corresponding volumetric determination, according to the standard and cold filter plugging point (CFPP) was obtained (when necessary) for the samples obtained [13][14][15]. Regarding the yield of the reactions and fatty acid methyl ester profile of biodiesel, gas chromatography was used, following the conditions explained elsewhere [16]. Thus, the yield of the biolubricant was obtained by the decrease in FAMEs. Equally, the reaction evolution was also assessed by FAME decrease. Finally, oxidative stability was calculated through the Rancimat method, obtaining the induction point (expressed in hours) [17].

Kinetic Study and Volume Reactor Determination
In order to design a reactor for biolubricant production from biodiesel, different factors were taken into account, as follows: • Characterization of the biolubricant, which determines its use in industry, paying special attention to viscosity at 40 and 100 °C and viscosity index and comparing with traditional lubricants.

•
Determination of the order of reaction of the second transesterification process, obtaining the main kinetic parameters. For that purpose, the decrease in FAMEs from biodiesel to biolubricant production at different reaction times (0, 5, 15, 25 and 90 min) was measured to determine the reaction evolution. Different parameters like temperature and catalyst addition were studied. • Calculation of biolubricant production at a national level (in this case, Spain), comparing with the equivalent lubricant production at the same level.

•
Selection of the suitable reactor depending on production, obtaining the design equation by a mass balance. • Reactor volume determination.

Biodiesel Characterization
First, a characterization of the biodiesel used as a starting point for biolubricant production was carried out. As it can be observed in Table 2, this biofuel complied with almost all the parameters established in the UNE-EN 14,214 standard, except for the oxidative stability. Thus, the chemical conditions were suitable to obtain a high-yield product, showing a high flash point, which implies safety during its storage. Regarding FAME content of biodiesel, the main components were methyl oleate (66.4%), linoleate (16.6%), linolenate (9.1%) and palmitate (3.8%). The higher content in methyl oleate, which is a mono-unsaturated compound, could explain the higher oxidative stability compared to other biodiesel samples found in the literature, needing lower amounts of antioxidant additives to comply with the standard [18][19][20].

Reaction Optimization
Concerning the optimization of the second transesterification reaction, Figures 2-4 show the different conversion evolution depending on the different parameter changes (FAME/2-ethyl-2-hydroxymethyl-1,3-propanediol ratio; temperature and catalyst concentration, respectively). In the case of FAME/alcohol ratio, Figure 2 shows the main results, obtaining similar yields for 1:1 and 1:2 ratios. As the amount of alcohol increased, the reaction rate increased for the initial stages. The 1:2 ratio was discarded, as it implied higher costs and operational difficulties (due to the filtration process, which was more difficult as the amount of surplus alcohol was higher), selecting 1:1 ratio for further studies, as in the case of previous studies with similar alcohols (2,2-dimethyl-1,3-propanediol) [10]. Regarding temperature, Figure 3 shows the differences between the reaction temperatures selected for this experience. Thus, as temperature increased, the reaction took place at a faster rate, achieving high conversions (exceeding 80%) in 15 min, reaching a balance at 25 min for all cases. This behavior was also observed for other biolubricants (obtained from rapeseed biodiesel and 2-ethyl-1-hexanol), where temperature considerably increased the reaction rate at initial stages [2]. For the final production of biolubricant, the temperature at which the highest yield was obtained (exceeding 94%) was selected, that is, 120 • C, which is a relatively low reaction temperature compared to other biolubricant production processes observed in the literature, reaching 150 • C [18]. Finally, Figure 4 shows the effect of catalyst concentration. As expected, the rate of the reaction was higher as the catalyst addition was higher (and the equilibrium was achieved at shorter reaction times), reaching similar final conversions for 1.5 and 2.0% w/w. Consequently, an intermediate catalyst addition (1.5%) was selected, following economic criteria. Again, similar trends were obtained according to the literature, where the use of catalysts promoted the generation of products at faster reaction rates [2,10]. This way, it was proved that the use of catalysts is vital to make these processes competitive compared to traditional productions based on mineral products.

Biolubricant Characterization
In order to produce the final biolubricant, the best chemical conditions observed in the previous section (included in Table 1 for the second transesterification) were used, and the conversion evolution was included in Figure 5. Thus, conversion exceded 99%, with a considerable increase during the first 20 min of reaction (with around 97% at 15 min), obtaining a high-yield biolubricant under these chemical conditions. The main characteristics of the purified biolubricant are included in Table 3. According to the viscosity values obtained for 40 and 100 • C, this biolubricant could be comparable to a SAE10W30 lubricant (some characteristics of commercial lubricants are included in this table), which will be taken into account for the reactor design in following sections. Viscosity index, the resistance to viscosity changes over temperature, was relatively high compared to other products [21,22]. Flash and combustion points were high, within the range of other values found in the literature, which implies increased safety during storage [23,24]. According to data obtained for a SAE10W30 commercial lubricant (see Table 3), viscosity index and flash point were lower compared to these samples, although in the same order of magnitude.

Reactor Design
The choice of a continuous or discontinuous reactor mainly depends on the demand of the product obtained. In this case, a continuous reactor (steady state stirred tank reactor) was chosen, as it can be seen in Figure 6. The equation for reaction design is obtained through a mass balance. Figure 7 shows the basis for the design of a continuous reactor. Where F A0 and F A represent the mole flow of the limiting reagent at the reactor inlet and outlet, respectively; T E and T S are the temperatures of the stream at the reactor inlet and outlet, respectively; C A0 and C A are the concentrations of the limiting reactive at the reactor inlet and outlet, respectively.
If the balance of the limiting reactive is done, Equation (1) is obtained: were −r A is the reaction rate and dN A /dt is the accumulation of the limiting reagent. As the process is isothermal (there is not any density change) and it takes place in a steady state operation (dN A /dt = 0), Equation (1) can be written in the form of Equation (2), where X A represents the conversion of reactive A.
In order to determine the volume of the reactor, it is necessary to know the flow of reactive A at the reactor inlet or outlet, the reaction kinetics and the final conversion. Regarding conversion, it was considered to be 99.41%, that is, X A = 0.9941.
To determine F Ao , some data about biolubricants were considered. According to ASELUBE, 111,000 tonnes of lubricants for automotive industry were sold [25]. On the other hand, according to other sources, the percentages of sales of SAE10W30 lubricants was 5% [26]. As a consequence, and considering that the obtained biolubricant is similar to SAE10W30 characteristics (according to data provided in Table 3), a production of 5550 tm/year was established, in order to cover the market demand in Spain.
To determine the mole flow of reagent A (biodiesel or fatty acid methyl ester), it is necessary to know the average molecular weight of biodiesel through biodiesel composition. Table 4 shows these data, obtaining an average molecular weight for rapeseed biodiesel of 284.29 g/mol. Taking into account that the conversion of the process can be considered as 1 in practical terms, Table 5. shows the mass, volumetric and mole flows of biodiesel, 2-ethyl-2hydroxymethyl-1,3-propanediol (alcohol) and sodium methoxide (catalyst). These data were obtained according to a yearly typical production of SAE10W30 lubricant in a chemical plant located in Extremadura region [25]. Accordingly, biodiesel, alcohol and catalyst flows were calculated according to the optimum reaction conditions found in this experience (Biodiesel/alcohol mole ratio, 1/1 and catalyst concentration, 1.5%).

Kinetic Study
The kinetic determination was carried out by applying the data analysis integral method. The reaction can be described according to Equation (3), where A is FAMEs, B is 2-ethyl-2-hydroxymethyl-1,3-propanediol, C is the biolubricant obtained and D is methanol.
The reaction rate is given by Equation (4), where k is the kinetic constant, C A is biodiesel concentration, C B is alcohol concentration and α and β are the reaction orders for A and B. − Taking into account the relationship between the reactant and its conversion and the stoichiometry of the reaction, θ B can be defined as follows (Equations (5)-(7)): where C A0 and C B0 are the initial concentration of A and B, and X A is FAME conversion. Replacing Equations (5)- (7) in Equation (4), the following expression is obtained (Equation (8)).
This equation is differential, whose integration after establishing the following limits (t = 0, X A = 0) and (t = t, X A = X A ) was (Equation (9)): The solution of this equation depends on α and β values. Although there are many possibilities, in the case of transesterification (as described in the literature) a first-order reaction regarding both reactants can be considered. Thus, α and β = 1. On the other hand, for most experiments, a θ B value of 1 was selected, which implies that the alcohol concentration was three times higher than the stoichiometric value. Under these circumstances, the concentration of this reagent can be considered constant and the kinetic model, originally of 1 + 1 order, can be considered a pseudo-first order reaction, that is, a model where the starting kinetic equation is shown in Equation (10), where k* is a complex constant that includes catalyst and alcohol concentration effects.
Equation (10) can be expressed as a function of conversion, obtaining Equation (11): As a result of the integration of Equation (11), Equation (12) is finally obtained, which is a linear equation, so the slope of the line is k*.
Equation (12) has been tested with the experimental results obtained. For instance, Figure 8 shows the representation for one of the experiments carried out. Thus, the linear trend was observed for the first stages of the reaction, as the transesterification balance was achieved quickly and, under these circumstances, the net reaction rate (that is, the difference between the direct and inverse reaction rate) is close to zero. As a consequence, the reaction model represented in Equation (3), which assumes an irreversible reaction, does not make any sense at longer times, where the reversible reaction could take place.  (12) for the following chemical conditions: FAME/alcohol mole ratio, 1:1; catalyst concentration, 1.5% w/w; temperature, 100 • C; reaction time, 90 min; stirring rate, 500 rpm.
Thus, through least squares adjustment, Table 6 shows the pseudo-kinetic constant, which was obtained from the slope of the linear representation observed in Figure 8, and the coefficient of determination, corresponding to experiments with different catalyst concentration. As it can be observed, the pseudo-first reaction order fits reasonably, as long as the reaction time considered is short (for the first 15-20 min).
As it was pointed out, k* is a complex constant, including catalyst and alcohol concentrations. If the order of reaction of alcohol and catalyst is assumed to be the unit, k* would be expressed by Equation (13), and according to this equation, a representation of k* versus catalyst concentration should imply a line passing though the origin.
This representation is included in Figure 9, which shows a linear trend (with a low R 2 value), passing through the origin as the kinetic model predicted. k * = k·C cat ·C B (13) Figure 9. Dependence of k* with catalyst concentration (pseudo-first reaction order).
An important aspect in every kinetic study is the dependence of kinetic constants on temperature. To determine this dependence, experiments at 100, 110 and 120 • C were carried out, adjusting the evolution of the reaction to a pseudo-first reaction order. Thus, Table 7 shows the kinetic constants (k*) and R 2 coefficients corresponding to this experimental set. The adjustment is reasonable, especially for the two first experiments, whereas in the case of the experiment carried out at the highest temperature the adjustment was worse, due to the fact that the balance was achieved earlier (at shorter reaction times) and, therefore, fewer experimental values were suitable to carry out the adjustment. Thus, points at longer reaction times (30-90 min) cannot be considered as the inverse reaction started to be considerable (which is not assumed in the kinetic model included in this study). Although k* is a complex constant, it is possible to assess its dependence on temperature. Thus, considering Equation (13), k constant would be the real kinetic constant. Assuming that k presents a dependence on temperature given by the Arrhenius equation, Equation (14) can be obtained from Equation (13), where all the terms that remain constant with temperature were regrouped, obtaining (Cte).
Processes 2021, 9, 1224 12 of 15 Calculating the Napierian logarithm, Equation (15) is obtained: where E is the activation energy, R is the ideal gas constant, T is absolute temperature and (Cte) is a constant where all the constant terms are included, including the pre-exponential factor k 0 . From Equation (15), through a least square adjustment of data corresponding to the experiments carried out at different temperatures, it is possible to determine the activation energy of the process. Table 8 shows the adjustment data, the activation energy and R 2 value. It should be pointed out the low activation energy, which is in accordance with the fact that it is a catalytic reaction carried out at relatively high temperatures (100-120 • C). This low activation energy value implies a very low dependence of the reaction rate on temperature. Thus, in the temperature range that has been studied, the reaction rate is not very sensitive to temperature. In any case, and even though according to previous results an optimization of reaction temperature would be justified, the reactor design, as explained in the following section, was carried out for a reaction temperature of 120 • C. Finally, considering the kinetic data obtained, the theoretical conversion corresponding to the final chemical conditions was determined. Thus, in Figure 10, a comparison between the experimental and theoretical data is included, pointing out a good convergence and ratifying the reliability of the kinetic data obtained. As it can be seen, the adjustment is especially good at the beginning of the reaction, for the first 10 min, where high conversions were observed. As a consequence, the kinetic model proposed could be suitable for the subsequent reactor design, as explained in further sections. Thus, the reactor design will be carried out with the abovementioned chemical conditions corresponding to the final biolubricant production, that is, 1:1 mole ratio, 1.5% catalyst concentration, 120 • C, a reaction time of 90 min and a stirring rate of 500 rpm.

Reactor Volume
Reactor volume was determined by Equation (2). In order to apply this equation, F A0 , X A and -r A should be known. F A0 was previously determined (18,618.75 kmol/year). If expressed in international units, 0.5904 mol/s are obtained. For X A , the conversion obtained in the final experiment was considered (0.9941). In order to calculate -r A , Equation (10) was used, which implies the calculation of C A and k*. C A value, that is, the concentration of fatty acid methyl esters at the reactor outlet, can be obtained from Equation (5), previously obtaining C A0 with total volumetric (and mole) flow (included in Table 5). This way, C A0 = 2182.3 mol/m 3 and C A = 12.9 mol/m 3 . Finally, k* was obtained by using Equation (12) as previously explained. Table 9 shows the calculations made. Table 9. Reactor volume determination.

Conclusions
The main findings were the following:

•
Rapeseed biodiesel complied with almost all the standard requirements, highlighting the high flash and combustion points (exceeding 190 • C). However, oxidative stability did not comply with the standard (5.37 h), although it was relatively high compared to other equivalent samples found in the literature.

•
The characteristics of the biolubricant were calculated, being equivalent to a SAE 10W30 lubricant, suitable for Diesel engines.

•
For the reactor design, due to the operating conditions, the kind of reaction and desired production, a continuous reactor was selected.

•
The kinetics of the process could be considered a pseudo-first order reaction, obtaining an activation energy of 6.91 kJ·mol −1 . • According to the annual production (once the kinetics was known), a reactor volume of 9.66 m 3 was obtained (12 m 3 if oversizing is considered). Funding: This research was funded by the "Junta de Extremadura" and "FEDER" (IB18028 and GR18150).
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Acknowledgments: The authors would like to thank the "Junta de Extremadura" and "FEDER" for the financial support received to perform this study.

Conflicts of Interest:
The authors declare no conflict of interest.