2. Biodiesel feedstock and the reactor device

While corn, sunflower, and palm oil are readily available in the market, their utilization may raise a strong debate on land for food vs. land for energy. To rule out this debate, waste cooking oil (WCO) is used instead, because it is abundantly available feedstock, inexpensive and most of the time its disposal into sewerage systems is associated with environmental concerns. The supply chain of collecting this abundant source is beyond the scope of this chapter. However, it is important to state that residential communities are in general in favor of trapping this problematic sewerage source that is responsible for clogging the plumbing systems. Collection also can be facilitated by using drum type or smaller plastic containers provided to the local restaurant, school/university canteens, residential communities as a privately owned small business, or through municipality.

The collected WCO generally requires pretreatment that can be facilitated at a moderate temperature to maintain its liquid and less viscous form. One can use the heat of the summer (45–50C) at which unsaturated and saturated steric fatty acid stays in liquid form, and 10– 20 μm filtration can be used to eliminate any suspended oil solid residuals. Dehydrating of the WCO is also required, during which any water content brought by the processed food into the waste cooking oil is liberated through evaporation. At the laboratory scale, stirring and heating pad at near 110C for several hours can perfect this task. Process methanol and catalyst NaOH or KOH can also be substituted with commercially available grade instead of high purity pharmaceutical grade/Sigma-Aldrich that can also leverage process economically. Once pretreatment of the feedstock is done, the NaOH solid catalyst in the form of small ballets is dissolved into methanol at the stipulated ratio, i.e. 0.5–1% by mass of WCO. This process can be facilitated under moderate heating and a temperature below 60C and stirring forming the meth oxide reactant solution. In the lab, multiple transesterification reaction experiments can be conducted simultaneously under the same temperature and stirring rate to reduce experimental sequence and human error. This can be carried out using a multiple dissolution apparatus such as those provided by Agilent Technologies, featuring 6–12 reactor vessels as depicted in Figure 1 [6]. The caps are tightly fitted and are equipped with direct access ports for sampling without process interruption. They all are also set on thermally controlled wells. The process in these individual small-scale batch reactors resemble those carried out in larger

promotes the desired forward reaction [2]. Contrary to well-known hydrocarbon fuels that are characterized by fixed thermodynamic and physical properties, the TG, DG, MG have no fixed chemical formula and neither their thermodynamic properties, such as standard enthalpy or specific heats, nor physical ones like density or viscosity, are consistent throughout the literature [3, 4]. Therefore, material characterization is an essential step in the modeling of the transesterification process. The extent of these properties depends on the complexity and comprehensiveness of the simulation, from a simple incompressible flow that requires only viscosity and density, to a complex non-isothermal flow that requires heat of formation, specific heat, thermal conductivity and their associated diffusions. However, as these properties can be derived following the American Society for Testing and Materials (ASTM) standards, their reactions are more complex. Table 1 summarizes some of the utilized properties for the waste oil, TG, DG, and MG used in the work of Noureddine and Zhu who were

Molecular weight

Methanol CH4O 32 3.96E-4 1.470E3 791.8

Diglyceride C37H72O5 596 — — 880 Monoglyceride C20H40O4 344 — — 875 Biodiesel C18H36O6 284 1.12E-3 1.187E3 870 Glycerol C3H9O3 93 1.412E0 0238.6 1261

Viscosity (kg/m.s)

C54H105O6 848 1.61E-2 2.2E3 883.3

Cp (J/kg.C) Density (kg/m<sup>3</sup> )

(g/mol)

Setting up a reaction mechanism of numerous species or elements and hundreds of reactions, while accounting for reaction radicals, is rather impractical for engineers. The overall reaction can be captured through well controlled conditions and yield assessment procedures that can save the pain of the development or use these reaction mechanisms. This chapter undertakes the conventional transesterification at different process temperatures, highlighting their influ-

While corn, sunflower, and palm oil are readily available in the market, their utilization may raise a strong debate on land for food vs. land for energy. To rule out this debate, waste cooking oil (WCO) is used instead, because it is abundantly available feedstock, inexpensive and most of the time its disposal into sewerage systems is associated with environmental concerns. The supply chain of collecting this abundant source is beyond the scope of this

amongst the pioneers of quantifying transesterification reaction kinetics [5].

ence on the yield and their distribution inside the reactor.

Species Chemical

48 Biofuels - Challenges and opportunities

Waste oil or Triglyceride formula

Table 1. Summary of species properties and molecular weight (MW) [5].

2. Biodiesel feedstock and the reactor device

Figure 1. Dissolution apparatus representing eight multiple 1-L batch reactors and HomeBiodeselKits 500-liter batch reactor.

3. Transesterification method

Figure 3. Transesterification reaction mechanism.

Figure 3 shows the three-step details of the three equilibrium reactions process of transester-

This equilibrium given in Figure 1 is represented by six coupled first order differential equations including the overall shunt reaction as per Eqs. (1)–(6) [10] and these are written as:

d DG ½ �=dt ¼ �k3½ � DG ½ �þ AL k4½ � E ½ �þ MG k1½ � TG ½ � AL –k2½ � E ½ � DG (2)

d MG ½ �=dt ¼ �k5½ � MG ½ �þ AL k6½ � E ½ �þ GL k3½ � DG ½ � AL –k4½ � E ½ � MG (3)

–k8½ � <sup>E</sup> <sup>3</sup>

Chemical Kinetic and High Fidelity Modeling of Transesterification

http://dx.doi.org/10.5772/intechopen.80008

d AL ½ �=dt ¼ �d E½ �=dt (6)

ð1Þ

51

ð4Þ

½ � GL (5)

ification reaction as indicated by Noureddine et al. [5] and Janajreh et al. [9].

d GL ½ �=dt <sup>¼</sup> <sup>k</sup>5½ � MG ½ � AL –k6½ � <sup>E</sup> ½ �þ GL <sup>k</sup>7½ � TG ½ � AL <sup>3</sup>

In these equations, k is the rate constant of each of the forward and backward reaction of Figure 3 and (x) is the concentration of species x, the TG, DG, MG, AL, E, and GL are

Figure 2. Patented tubular continuous biodiesel reactor 1–8L/hr [19].

scale ones provided by HomeBiodeselKits, as depicted in Figure 1. Beside batch method, a patented continuous reactor is also depicted in Figure 2 featuring compactness, continuity, and scalability. This reactor reported better efficacy observed by the higher throughput and better quality of product than the batch reactor.

The continuous reactor consists of a tubular reactor of two or more upright concentric cylinders. Each is equipped with inflow at the bottom and outflow at the top port, connected by 5-mm ϕ chemically resistive hoses. The inflow and outflow are configured circumferentially in each cylinder, rendering the flow more residence time due to their helical trajectory. The details of the reactor's dimensions and geometry can be found in the work of Janajreh et al. [7]. A peristaltic or a diaphragm pump type is used to inject the two reactants, i.e. pretreated oil and methoxide into the reactor at the stipulated molar ratio. The mass flow rate is typically being set beyond laminar limits for turbulence mixing that help the transesterification reaction. The mixture from the batch dissolution is periodically recovered in 10 mm vials, and is refrigerated to halt their progressive reactions for downstream species analysis. The analysis of the compositions of TG, DG, MG, FAME, glycerol (GL) and alcohol (AL) is carried out using a standard Gas Chromatography Mass Spectrometry (GC/MS) equipment as those provided by Thermo Scientific DSQ II equipped with a flame ionization detector.

In the GC/MS analytical equipment, the FAME column is initially calibrated using standard biodiesel and glycerol samples to ensure precise qualitative and quantitative analysis. The collected and refrigerated vials are obtained at numerous and progressive reaction intervals of 15-time steps over 2 hr of reaction time [8]. The breakdown of sample species composition is acquired using methods similar to those carried out by Noureddini and Zhu [5] and Janajreh et al. [7]. Transesterification is performed at 6:1 alcohol to oil molar ratio and 0.5% NaOH by mass of WCO, and both at different temperatures of 50 and 60C.

#### 3. Transesterification method

scale ones provided by HomeBiodeselKits, as depicted in Figure 1. Beside batch method, a patented continuous reactor is also depicted in Figure 2 featuring compactness, continuity, and scalability. This reactor reported better efficacy observed by the higher throughput and

The continuous reactor consists of a tubular reactor of two or more upright concentric cylinders. Each is equipped with inflow at the bottom and outflow at the top port, connected by 5-mm ϕ chemically resistive hoses. The inflow and outflow are configured circumferentially in each cylinder, rendering the flow more residence time due to their helical trajectory. The details of the reactor's dimensions and geometry can be found in the work of Janajreh et al. [7]. A peristaltic or a diaphragm pump type is used to inject the two reactants, i.e. pretreated oil and methoxide into the reactor at the stipulated molar ratio. The mass flow rate is typically being set beyond laminar limits for turbulence mixing that help the transesterification reaction. The mixture from the batch dissolution is periodically recovered in 10 mm vials, and is refrigerated to halt their progressive reactions for downstream species analysis. The analysis of the compositions of TG, DG, MG, FAME, glycerol (GL) and alcohol (AL) is carried out using a standard Gas Chromatography Mass Spectrometry (GC/MS) equipment as those provided by Thermo

In the GC/MS analytical equipment, the FAME column is initially calibrated using standard biodiesel and glycerol samples to ensure precise qualitative and quantitative analysis. The collected and refrigerated vials are obtained at numerous and progressive reaction intervals of 15-time steps over 2 hr of reaction time [8]. The breakdown of sample species composition is acquired using methods similar to those carried out by Noureddini and Zhu [5] and Janajreh et al. [7]. Transesterification is performed at 6:1 alcohol to oil molar ratio and 0.5% NaOH by

better quality of product than the batch reactor.

Figure 2. Patented tubular continuous biodiesel reactor 1–8L/hr [19].

50 Biofuels - Challenges and opportunities

Scientific DSQ II equipped with a flame ionization detector.

mass of WCO, and both at different temperatures of 50 and 60C.

Figure 3 shows the three-step details of the three equilibrium reactions process of transesterification reaction as indicated by Noureddine et al. [5] and Janajreh et al. [9].

This equilibrium given in Figure 1 is represented by six coupled first order differential equations including the overall shunt reaction as per Eqs. (1)–(6) [10] and these are written as:

$$d[TG]/dt = -k\_1[TG][AL] + k\_2[E][DG] - k\_7[TG][AL]^3 + k\_8[E]^3[GL] \tag{1}$$

$$d[\text{DG}]/dt = -k\_3[\text{DG}][AL] + k\_4[E][MG] + k\_1[TG][AL] - k\_2[E][DG] \tag{2}$$

$$d[\text{MG}]/dt = -k\_5[\text{MG}][AL] + k\_6[E][GL] + k\_3[\text{DG}][AL] - k\_4[E][MG] \tag{3}$$

ð4Þ

$$d[\text{GL}]/dt = k\_5[\text{MG}][AL] - k\_6[E][\text{GL}] + k\_7[\text{TG}][AL]^3 - k\_8[E]^3[\text{GL}] \tag{5}$$

$$d[AL]/dt = -d[E]/dt\tag{6}$$

In these equations, k is the rate constant of each of the forward and backward reaction of Figure 3 and (x) is the concentration of species x, the TG, DG, MG, AL, E, and GL are


Figure 3. Transesterification reaction mechanism.

respectively the Triglyceride, Diglyceride, Monoglyceride, Alcohol, Fatty Acid Methyl Ester and Glycerol. This system of PDEs Eqs. (1)–(6) is solved for k1 through k8 at minimum root mean square following the temporal measurements of each specie and using Eq. (7) as:

$$A\mathbf{x} = b \text{ or } \mathbf{x} = A^{-1}b \tag{7}$$

where A is the coefficient matrix of the species concentration at each time step obtained through GC/MS analysis, x is the evaluated rate constant vector (k1 through k8) and b is obtained by evaluating the time derivative of the measured concentrations. The rate constant (k) of each reaction is expressed by the Arrhenius Eq. (8) [11] as follows:

$$k = A\_0 e^{-E/RT} \text{ or } \ln k = \ln A\_0 - \frac{E}{R} \frac{1}{T} \tag{8}$$

role to kinetically propel these reactions. In general, an increase in temperature provides a higher energy state for the molecules leading to more collisions and also improves the solubility of the reactants. As these reactions are reversible, the consumption of any of the reactants

Figure 4. Transesterification of WCO conversion results at 50C (top left) and 60C (top right). Corresponding chemical

Chemical Kinetic and High Fidelity Modeling of Transesterification

http://dx.doi.org/10.5772/intechopen.80008

53

As can be shown in Figure 4 of the GC/MS results obtained at the two temperatures, DG and MG are intermediate species. They appear following mass transfer limitation for a relatively shorter time, before they asymptotically vanish. The amount of unreacted alcohol at T = 60C is shown to be lower than those obtained at T = 50C, as per Figure 4. On the contrary, the FAME production and its byproduct GL is noticeably higher at higher temperature (T = 60C) than the lower process temperature (T = 50C). The rate constants for the two process temperatures are depicted side by side in the bottom left 'staked bar plot' of Figure 4. These values are relatively comparable, except for k2 and k4. Considering the lower temperature stacked bar (left side), it shows the low value of k1 and k2 reaction constants compared to the rest of reaction constants placing them as the rate limiting reactions. As soon as the product of these reactions becomes available, higher rates for their intermediate yield takeover. Comparing the lower and higher temperature ks, the higher temperature results in higher k1 and k2, which are the rate limiting reactions, while the rest are comparable. This explains the higher initial yield beyond the twophase mass transfer limitation that lasts within the first minutes of the reaction. These evaluated activation energies are also depicted in the bar plot of Figure 4 on the bottom right. These

must be avoided to steer and ensure these reactions of Figure 3 forward.

kinetic reaction rate constants (bottom left) and activation energy (bottom right).

where Ao (s�<sup>1</sup> ) is the pre-exponential constant, E (J/mol) is the activation energy corresponding to each of the eight reactions, R (J/mol K) is the molar universal gas constant and T is the process temperature expressed in degree Kelvin. Therefore, from the slope and intercept of the natural log line of K and 1/T one can infer the activation energy (which is normalized by the gas constant) and the natural log of the pre-constant, respectively. Following previous work on chemical kinetics [5, 12], the reaction rate constants Eqs. (1)–(8) are solved simultaneously using MATLAB which enables direct inference of E and Ao [13].

### 4. Transesterification kinetics evaluation

Transesterification is a slowly reversible reaction, where only the forward reactions are desired for the production of biodiesel. It is easy to say that lower activation energies are favored for forward reactions (E1, E3, E5, E7) while higher activation energies are preferred for backward reactions (E2, E4, E6, E8). Unfortunately, this is not as simple as will be stated below. The obtained GC/MS results for this conventional transesterification at 0.5% NaOH, 6:1 molar ratio and temperature of 50 and 60�C, in terms of mole fraction, are depicted in Figure 4. It shows the rates of consumption of TG and Alcohol into the formation of methylesters (E) and GL as well as the intermediate species, i.e. DG and MG.

The E production and its byproduct GL demonstrate a sigmoidal trend: a slow delay in the reaction rate at the start is then followed by a quick surge in the production, and a final slow rate near the reaction completion. This trend characterizes the changing mechanisms reaction that swings from the slow mass transfer controlled to the fast kinetically controlled and back to slow reactivity as equilibrium is attained [5]. This can be explained by the poor diffusion between the two phase reactants system, i.e. TG and AL. As soon as the methyl ester is formed the solubility of the two phases into one single phase is increased drastically which causes the production surge. This is stipulated by any sample taken at the first minute of the reaction, which tends to form two separate phases. Samples taken thereafter appear to be a more homogenized single phase as sampling time progress. At the onset of the formation of the single phase, mixing becomes less important and hence temperature can play more important Chemical Kinetic and High Fidelity Modeling of Transesterification http://dx.doi.org/10.5772/intechopen.80008 53

respectively the Triglyceride, Diglyceride, Monoglyceride, Alcohol, Fatty Acid Methyl Ester and Glycerol. This system of PDEs Eqs. (1)–(6) is solved for k1 through k8 at minimum root mean square following the temporal measurements of each specie and using Eq. (7) as:

Ax <sup>¼</sup> b or x <sup>¼</sup> <sup>A</sup>�<sup>1</sup>

where A is the coefficient matrix of the species concentration at each time step obtained through GC/MS analysis, x is the evaluated rate constant vector (k1 through k8) and b is obtained by evaluating the time derivative of the measured concentrations. The rate constant

to each of the eight reactions, R (J/mol K) is the molar universal gas constant and T is the process temperature expressed in degree Kelvin. Therefore, from the slope and intercept of the natural log line of K and 1/T one can infer the activation energy (which is normalized by the gas constant) and the natural log of the pre-constant, respectively. Following previous work on chemical kinetics [5, 12], the reaction rate constants Eqs. (1)–(8) are solved simultaneously

Transesterification is a slowly reversible reaction, where only the forward reactions are desired for the production of biodiesel. It is easy to say that lower activation energies are favored for forward reactions (E1, E3, E5, E7) while higher activation energies are preferred for backward reactions (E2, E4, E6, E8). Unfortunately, this is not as simple as will be stated below. The obtained GC/MS results for this conventional transesterification at 0.5% NaOH, 6:1 molar ratio and temperature of 50 and 60�C, in terms of mole fraction, are depicted in Figure 4. It shows the rates of consumption of TG and Alcohol into the formation of methylesters (E) and GL as

The E production and its byproduct GL demonstrate a sigmoidal trend: a slow delay in the reaction rate at the start is then followed by a quick surge in the production, and a final slow rate near the reaction completion. This trend characterizes the changing mechanisms reaction that swings from the slow mass transfer controlled to the fast kinetically controlled and back to slow reactivity as equilibrium is attained [5]. This can be explained by the poor diffusion between the two phase reactants system, i.e. TG and AL. As soon as the methyl ester is formed the solubility of the two phases into one single phase is increased drastically which causes the production surge. This is stipulated by any sample taken at the first minute of the reaction, which tends to form two separate phases. Samples taken thereafter appear to be a more homogenized single phase as sampling time progress. At the onset of the formation of the single phase, mixing becomes less important and hence temperature can play more important

) is the pre-exponential constant, E (J/mol) is the activation energy corresponding

(k) of each reaction is expressed by the Arrhenius Eq. (8) [11] as follows:

using MATLAB which enables direct inference of E and Ao [13].

4. Transesterification kinetics evaluation

well as the intermediate species, i.e. DG and MG.

where Ao (s�<sup>1</sup>

52 Biofuels - Challenges and opportunities

b (7)

ð8Þ

Figure 4. Transesterification of WCO conversion results at 50C (top left) and 60C (top right). Corresponding chemical kinetic reaction rate constants (bottom left) and activation energy (bottom right).

role to kinetically propel these reactions. In general, an increase in temperature provides a higher energy state for the molecules leading to more collisions and also improves the solubility of the reactants. As these reactions are reversible, the consumption of any of the reactants must be avoided to steer and ensure these reactions of Figure 3 forward.

As can be shown in Figure 4 of the GC/MS results obtained at the two temperatures, DG and MG are intermediate species. They appear following mass transfer limitation for a relatively shorter time, before they asymptotically vanish. The amount of unreacted alcohol at T = 60C is shown to be lower than those obtained at T = 50C, as per Figure 4. On the contrary, the FAME production and its byproduct GL is noticeably higher at higher temperature (T = 60C) than the lower process temperature (T = 50C). The rate constants for the two process temperatures are depicted side by side in the bottom left 'staked bar plot' of Figure 4. These values are relatively comparable, except for k2 and k4. Considering the lower temperature stacked bar (left side), it shows the low value of k1 and k2 reaction constants compared to the rest of reaction constants placing them as the rate limiting reactions. As soon as the product of these reactions becomes available, higher rates for their intermediate yield takeover. Comparing the lower and higher temperature ks, the higher temperature results in higher k1 and k2, which are the rate limiting reactions, while the rest are comparable. This explains the higher initial yield beyond the twophase mass transfer limitation that lasts within the first minutes of the reaction. These evaluated activation energies are also depicted in the bar plot of Figure 4 on the bottom right. These


5. High-fidelity transesterification model development

finite volume CFD based software of Ansys/Fluent 17.1 [18].

∂ ∂t ð Þþ ∅

(ε). These two equations in steady state are expressed as:

5.1. Governing equations

(rui

In order to sharpen our understanding of the transesterification reaction progression and species distribution, a high fidelity reactive flow model for transesterification is developed. The basis of the model is computational fluid dynamics (CFD). Several research software have adopted CFD environment to solve multiple physics, whereby some are open sources and others are more commercially tuned. The transesterification model is developed within the

We here focus on the modeling of the tubular reactor that depicted earlier in Figure 2 which has numerous advantages over the batch reactor, i.e. continuous, compactness, better yield [19]. Modeling involves the application of flow continuity, momentum, and energy equations. Furthermore, the flow is characterized as a mixture of multiple reacting species, incompressible, viscous, and turbulent. The reaction is assumed to start as soon as the reactant components are met inside the reactor. The flow is governed by the Navier-Stokes equation, which is

associated with temporal, advective, viscous, and any source term and is written as:

ð Þ¼ ui<sup>∅</sup> �<sup>∂</sup> ∂xi

where ui is the velocity and S<sup>∅</sup> is the source term due to the dispersed phase interaction. The ∅ is the flow dependent variable and when represents the density (r), velocity density multiple

The left terms govern the advective while the right hand terms govern respectively the generation, the diffusion, and destruction of the turbulent quantities. In these equations, μ<sup>t</sup> is the

where f <sup>μ</sup> and C<sup>μ</sup> are flow dependent constants and C1ε, C2ε,σ<sup>k</sup> and σ<sup>s</sup> are tuning empirical constants. The transport equations that govern the of species m<sup>i</sup> take the following form:

eddy viscosity parameter and it overwhelms the laminar viscosity and is written as:

) and the temperature (T) it yields the continuity, the momentum, and the energy equations, respectively. It takes also the turbulent scalars, i.e. kinetic energy (k) and dissipation rate

Γ<sup>∅</sup> ∂∅ ∂xi 

þ S<sup>∅</sup> (9)

Chemical Kinetic and High Fidelity Modeling of Transesterification

http://dx.doi.org/10.5772/intechopen.80008

55

ð10Þ

ð11Þ

ð12Þ

∂ ∂xi

Table 2. Summary of the rate constants values of the reactions at two temperatures and their activation energy.

are consistent with those reported experimentally by Feedman and coworkers [14, 15]. Based on the Arrhenius equations, they reported activation energies in the transesterification of the soybean oil in the range of 33,472–83,400 J/mol [16, 17]. These values are tabulated in Table 2, along with the rate constant. The activation energy is the highest for the reverse of the first reaction, DG + E—TG + AL suggesting this reaction is not as spontaneous as the others and reverse formation of TG is rather less expected. This is due to surplus of the Methanol alcohol that forces the forward reactivity of the TG + AL—DG + E reaction. This is also evident from all the forward reactions which are at lower activation energy than their associated reversal reaction.

The evaluated kinetics will be utilized in the development of a high fidelity reaction model. The model will be based on reactive flow of the six species, TG, AL, DG, MG, E and GL and their eight reactions. These are types of volume homogenous reactions. The considered reactor will be a tube type subjected to appropriate conditions and enabling the assessment of the yield and species distribution.

In summary, the required kinetics that governs the reactions and change of species has been evaluated. Results are intuitively correct, as a higher temperature results in higher rate constants. These data will be used in the high fidelity reactive flow analysis that will be detailed next. It is important to emphasize that importance of these simulations in capturing the overall species and their distribution in the reaction device. This is becoming the tool for the development of the new and innovative reactors. The observations procured from performing the chemical kinetic study enables the shift in conventional transesterification towards a modified sonication kinetic data for further reactive flow analysis. Comparison between regular mixing and ultra- sound assisted procedures proves the supremacy of the latter in terms of fastening the forward reactions by increasing reaction rates and decreasing activation energies of forward reactions. Note that the comparative study is performed to obtain an idea of the trend of sonication and not to get exact chemical kinetics.

#### 5. High-fidelity transesterification model development

In order to sharpen our understanding of the transesterification reaction progression and species distribution, a high fidelity reactive flow model for transesterification is developed. The basis of the model is computational fluid dynamics (CFD). Several research software have adopted CFD environment to solve multiple physics, whereby some are open sources and others are more commercially tuned. The transesterification model is developed within the finite volume CFD based software of Ansys/Fluent 17.1 [18].

#### 5.1. Governing equations

are consistent with those reported experimentally by Feedman and coworkers [14, 15]. Based on the Arrhenius equations, they reported activation energies in the transesterification of the soybean oil in the range of 33,472–83,400 J/mol [16, 17]. These values are tabulated in Table 2, along with the rate constant. The activation energy is the highest for the reverse of the first reaction, DG + E—TG + AL suggesting this reaction is not as spontaneous as the others and reverse formation of TG is rather less expected. This is due to surplus of the Methanol alcohol that forces the forward reactivity of the TG + AL—DG + E reaction. This is also evident from all the forward reactions which are at lower activation energy than their associated reversal

0.001 0.001 3089 77,983,755.8

Rate constant at T = 60C

TG + AL—E + DG 0.1213 0.1157 4103 77,983,755.8 E + DG—TG + AL 0.1167 0.3276 90,540 449,020.085 DG + AL—E + MG 0.8177 0.8947 7896 1.0976 E12 E + MG—DG + AL 3.5379 3.1354 10,597 4570,1701.3 MG + AL—E + GL 1.5943 1.7922 10,267 767.86 GL + E—MG + AL 0.1661 0.1322 19,980 32,870.58

GL + 3E—TG + 3A 0.8149 0.7568 6490 449,020.085

Table 2. Summary of the rate constants values of the reactions at two temperatures and their activation energy.

Activation energy

Pre-constant (sec<sup>1</sup> )

(J/mol)

The evaluated kinetics will be utilized in the development of a high fidelity reaction model. The model will be based on reactive flow of the six species, TG, AL, DG, MG, E and GL and their eight reactions. These are types of volume homogenous reactions. The considered reactor will be a tube type subjected to appropriate conditions and enabling the assessment of the

In summary, the required kinetics that governs the reactions and change of species has been evaluated. Results are intuitively correct, as a higher temperature results in higher rate constants. These data will be used in the high fidelity reactive flow analysis that will be detailed next. It is important to emphasize that importance of these simulations in capturing the overall species and their distribution in the reaction device. This is becoming the tool for the development of the new and innovative reactors. The observations procured from performing the chemical kinetic study enables the shift in conventional transesterification towards a modified sonication kinetic data for further reactive flow analysis. Comparison between regular mixing and ultra- sound assisted procedures proves the supremacy of the latter in terms of fastening the forward reactions by increasing reaction rates and decreasing activation energies of forward reactions. Note that the comparative study is performed to obtain an idea of the trend of

reaction.

TG + 3AL— 3E + GL

yield and species distribution.

Designated K or E Rate constant at

54 Biofuels - Challenges and opportunities

T = 50C

sonication and not to get exact chemical kinetics.

We here focus on the modeling of the tubular reactor that depicted earlier in Figure 2 which has numerous advantages over the batch reactor, i.e. continuous, compactness, better yield [19]. Modeling involves the application of flow continuity, momentum, and energy equations. Furthermore, the flow is characterized as a mixture of multiple reacting species, incompressible, viscous, and turbulent. The reaction is assumed to start as soon as the reactant components are met inside the reactor. The flow is governed by the Navier-Stokes equation, which is associated with temporal, advective, viscous, and any source term and is written as:

$$\frac{\partial}{\partial t}(\mathfrak{Q}) + \frac{\partial}{\partial \mathbf{x}\_i}(u\_i \mathfrak{Q}) = \frac{-\mathfrak{q}}{\mathfrak{d}\mathbf{x}\_i} \left(\Gamma\_{\mathcal{B}} \frac{\partial \mathfrak{Q}}{\partial \mathbf{x}\_i}\right) + \mathfrak{S}\_{\mathcal{Q}} \tag{9}$$

where ui is the velocity and S<sup>∅</sup> is the source term due to the dispersed phase interaction. The ∅ is the flow dependent variable and when represents the density (r), velocity density multiple (rui ) and the temperature (T) it yields the continuity, the momentum, and the energy equations, respectively. It takes also the turbulent scalars, i.e. kinetic energy (k) and dissipation rate (ε). These two equations in steady state are expressed as:

$$\begin{split} \rho u\_{\boldsymbol{i}} \frac{\partial \boldsymbol{k}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} &= \mu\_{\boldsymbol{i}} \Big( \frac{\partial \boldsymbol{u}\_{\boldsymbol{i}}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} + \frac{\partial \boldsymbol{u}\_{\boldsymbol{i}}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} \Big) \frac{\partial \boldsymbol{u}\_{\boldsymbol{i}}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} + - \frac{\partial}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} \Big( \frac{\mu\_{\boldsymbol{e}}}{\sigma\_{\boldsymbol{e}}} \frac{\partial \boldsymbol{k}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} \Big) - \rho \boldsymbol{e} \\ \rho \boldsymbol{u}\_{\boldsymbol{i}} \frac{\partial \boldsymbol{\mathcal{E}}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} &= \boldsymbol{C}\_{\boldsymbol{1}\boldsymbol{e}} \frac{\mu\_{\boldsymbol{e}} \boldsymbol{e}}{k} \Big( \frac{\partial \boldsymbol{u}\_{\boldsymbol{e}}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} + \frac{\partial \boldsymbol{u}\_{\boldsymbol{i}}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} \Big) \frac{\partial \boldsymbol{u}\_{\boldsymbol{i}}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} + \frac{\partial}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} \Big( \frac{\mu\_{\boldsymbol{e}}}{\sigma\_{\boldsymbol{e}}} \frac{\partial \boldsymbol{\mathcal{E}}}{\partial \boldsymbol{x}\_{\boldsymbol{i}}} \Big) - \boldsymbol{C}\_{\boldsymbol{1}\boldsymbol{e}} \frac{\rho \boldsymbol{e}}{k}^{\boldsymbol{1}} \end{split} \tag{10}$$

The left terms govern the advective while the right hand terms govern respectively the generation, the diffusion, and destruction of the turbulent quantities. In these equations, μ<sup>t</sup> is the eddy viscosity parameter and it overwhelms the laminar viscosity and is written as:

$$
\mu\_i = f\_{\mu} \mathbf{C}\_{\mu} \rho k^{\mathbb{Z}} / \mathfrak{e} \tag{11}
$$

where f <sup>μ</sup> and C<sup>μ</sup> are flow dependent constants and C1ε, C2ε,σ<sup>k</sup> and σ<sup>s</sup> are tuning empirical constants. The transport equations that govern the of species m<sup>i</sup> take the following form:

$$\frac{\partial}{\partial t} \left( \rho \boldsymbol{m}\_{\boldsymbol{\cdot}} \right) + \frac{\partial}{\partial \boldsymbol{x}\_{\boldsymbol{\cdot}}} \left( \rho \boldsymbol{u}\_{\boldsymbol{\cdot}} \boldsymbol{m}\_{\boldsymbol{\cdot}} \right) = \frac{\partial}{\partial \boldsymbol{x}\_{\boldsymbol{\cdot}}} \left( \rho \boldsymbol{D}\_{\boldsymbol{x}, \boldsymbol{\cdot}} + \boldsymbol{\mu}\_{\boldsymbol{\cdot}} \right) \boldsymbol{S} \boldsymbol{c}\_{\boldsymbol{\cdot}} \right) \frac{\partial \boldsymbol{m}\_{\boldsymbol{\cdot}}}{\partial \boldsymbol{x}\_{\boldsymbol{\cdot}}} + \boldsymbol{R}\_{\boldsymbol{\cdot}} + \boldsymbol{S}\_{\boldsymbol{\cdot}} \tag{12}$$

Here Di,m is the diffusion coefficient of i specie within the m bulk species and Sct is the turbulent Schmidt number defined as the ratio of the eddy viscosity t to the eddy diffusivity Di,m. Eq. 10 incorporate an additional source term Ri which accounts for species reaction. This term is governed by the reaction stoichiometry of the specie mi and is written as:

$$\sum\_{\ell=1}^{N} \mathbf{v}\_{\iota,r}^{\prime} \mathbf{S}\_{\iota} \quad \bigoplus\_{k\_{b,r}}^{k\_{f,r}} \quad \sum\_{\iota=1}^{N} \mathbf{v}\_{\iota,r}^{\prime} \mathbf{S}\_{\iota} \tag{13}$$

The reaction rate is proportional to the concentration of the reaction species (both reactants and products) to an ordered of specified power coefficients that can be written as:

$$\mathcal{R}\_{i,r} \quad = \mathcal{M}\_{i,r} \left( \mathbf{v}\_{i,r}^{\boldsymbol{\sigma}} - \mathbf{v}\_{i,r}' \right) \left( \mathbf{k}\_f \prod\_{f=1}^{N} \mathbf{C}\_{j,r}^{\eta\_{i,r}^\*} - \mathbf{k}\_b \prod\_{f=1}^{N} \mathbf{C}\_{j,r}^{\eta\_{i,r}^\*} \right) \tag{14}$$

where kf (k1, k3, k5, k7) and kb (k2, k4, k6, k8) are the reaction rate constants evaluated according to Eq. (8) for the forward and backward of the above Arrhenius equation Eq. (13) that summarizes equations of Figure 3. The Cj are the molar concentration of jth specie to the order of its stoichiometric coefficient ν while η is the reaction order. The Mi is the molecular weight of species i [20]. It is a good practice to pursue CFD as a non-reactive flow first before accounting the influence of reactions that numerically may induce earlier instability for the solution. After achieving near steady and consistent solution and influence of reaction can be included and analyzed.

#### 5.2. CFD setup and mesh sensitivity analysis

The model of the considered tubular reactor is shown below in Figure 5 along with its dimensions as summarized in Table 3.

A hybrid mesh of hexagonal and pyramid type is used to maintain the size within the processing capacity of the current laptop for an engineer. The mesh is established in ANSYS Design Modeler, consisting of two concentric tubular chambers along with the connecting tubing as shown in Figure 6. A baseline and another two refined meshes were created to study the dependency and the goodness of the mesh influence on the solution. The properties of the reactants and products are summarized in Table 4 as per the work of Narvaez et al. [21]. The ideal mixture is used to determine the mixture properties from the known properties of the individual specie as per Eq. (15).

ð15Þ

inlet boundary condition. The material initially is assumed as an ideal mixture that has the viscosity and density property proportion of 6:1 methanol-to-waste cooking oil mole fraction. Results of the unreacted isothermal flow distribution within the reactor, or cold flow as referred to in some literature, is depicted in Figure 7. Velocity vector colored by the residence time is used for the three level meshes. Results of residence time parameter is used as the most pronounced parameter for the intended reactive fluid to assess the mesh accuracy. These results are summarized in Table 5 for the three levels of the mesh for comparison. The baseline residence time is within 2% deviation from the next refined mesh, hence it is used for the rest of the analysis. However, the coarse mesh showed a sustainable deviation of nearly 10% that fails

Component Outer cylinder (glass) Middle cylinder (metal) Inner cylinder (metal)

Chemical Kinetic and High Fidelity Modeling of Transesterification

http://dx.doi.org/10.5772/intechopen.80008

57

Outer diameter (mm) 59.9 42.0 26.2 Inner diameter (mm) 51.9 36.7 21.8 Thickness (mm) 4.0 2.7 2.2 Length (mm) 272 313 353

to be qualified as an accurate and representative mesh.

Table 3. Dimensions of the reactor geometry components.

Figure 5. Tubular reactor model.

The boundary conditions are assigned as flow, constant velocity at the inlet and isothermal fluid. The boundary at the top outlet was defined as fixed pressure outlet at 0-Pa gauge pressure while the boundary walls are subjected to zero velocity, i.e. no slip and no penetration. The turbulence is accounted for via standard k-ɛ mode and model constants were kept at their common values [22–24]. The setup was initialized at the same values of velocity from the

Figure 5. Tubular reactor model.

Here Di,m is the diffusion coefficient of i specie within the m bulk species and Sct is the turbulent Schmidt number defined as the ratio of the eddy viscosity t to the eddy diffusivity Di,m. Eq. 10 incorporate an additional source term Ri which accounts for species reaction. This

The reaction rate is proportional to the concentration of the reaction species (both reactants and

where kf (k1, k3, k5, k7) and kb (k2, k4, k6, k8) are the reaction rate constants evaluated according to Eq. (8) for the forward and backward of the above Arrhenius equation Eq. (13) that summarizes equations of Figure 3. The Cj are the molar concentration of jth specie to the order of its stoichiometric coefficient ν while η is the reaction order. The Mi is the molecular weight of species i [20]. It is a good practice to pursue CFD as a non-reactive flow first before accounting the influence of reactions that numerically may induce earlier instability for the solution. After achieving near steady and consistent solution and influence of reaction can be included and

The model of the considered tubular reactor is shown below in Figure 5 along with its

A hybrid mesh of hexagonal and pyramid type is used to maintain the size within the processing capacity of the current laptop for an engineer. The mesh is established in ANSYS Design Modeler, consisting of two concentric tubular chambers along with the connecting tubing as shown in Figure 6. A baseline and another two refined meshes were created to study the dependency and the goodness of the mesh influence on the solution. The properties of the reactants and products are summarized in Table 4 as per the work of Narvaez et al. [21]. The ideal mixture is used to determine the mixture properties from the known properties of

The boundary conditions are assigned as flow, constant velocity at the inlet and isothermal fluid. The boundary at the top outlet was defined as fixed pressure outlet at 0-Pa gauge pressure while the boundary walls are subjected to zero velocity, i.e. no slip and no penetration. The turbulence is accounted for via standard k-ɛ mode and model constants were kept at their common values [22–24]. The setup was initialized at the same values of velocity from the

ð13Þ

ð14Þ

ð15Þ

term is governed by the reaction stoichiometry of the specie mi and is written as:

products) to an ordered of specified power coefficients that can be written as:

analyzed.

56 Biofuels - Challenges and opportunities

5.2. CFD setup and mesh sensitivity analysis

dimensions as summarized in Table 3.

the individual specie as per Eq. (15).


Table 3. Dimensions of the reactor geometry components.

inlet boundary condition. The material initially is assumed as an ideal mixture that has the viscosity and density property proportion of 6:1 methanol-to-waste cooking oil mole fraction.

Results of the unreacted isothermal flow distribution within the reactor, or cold flow as referred to in some literature, is depicted in Figure 7. Velocity vector colored by the residence time is used for the three level meshes. Results of residence time parameter is used as the most pronounced parameter for the intended reactive fluid to assess the mesh accuracy. These results are summarized in Table 5 for the three levels of the mesh for comparison. The baseline residence time is within 2% deviation from the next refined mesh, hence it is used for the rest of the analysis. However, the coarse mesh showed a sustainable deviation of nearly 10% that fails to be qualified as an accurate and representative mesh.

5.3. Reactive flow analysis

5.4. Simulation results

Table 6.

adjustable discharge capacity of up to 5 L/min.

described by Eq. (16) and is written as:

residence time to complete the reaction.

The reactor tube is setup vertically and is subjected to inlet flow rate at Re = 6000 at the bottom of the tube into the inner reactor chamber. There, the two fluids of TG and methoxide (mixture of the methanol and the catalyst) are injected circumferentially at the stipulated molar/mass ratio and inlet temperature. The outflow is subjected to atmospheric pressure which is admitted at exit boundary condition of the outer tube chamber located at the top. The no-lip no

Mesh type Coarse Baseline (Baseline � 1.5) Fine mesh (Baseline � 1.89)

Chemical Kinetic and High Fidelity Modeling of Transesterification

http://dx.doi.org/10.5772/intechopen.80008

59

To maintain single-phase flow avoiding the complication of other reactivation, the case of reactive flow simulation is carried out below the boiling point of methanol at a temperature of 60�C (333 K). A steady-state solution is sought for the flow. This is achieved by ignoring the temporal term of the governing equation. The flow is introduced to the reactor by means of an external peristaltic or diaphragm pump operating at relatively low head of nearly 2-m to overcome the viscous shear stresses and head losses and the reactor vertical head at an

Based on the stoichiometry of transesterification reactive 3 moles of methanol are consumed with 1 mole of triglycerides for the making of 1 mole of biodiesel. Using this ratio(3:1 methanol to WCO molar), which also corresponds to AL/WCO mass fraction of 0.102/0.892, results on the species contour plots depicted in Figure 8 and the outlet fractions are summarized in

At this velocity and molar ration, a low conversion of 28% is achieved. The conversion is

where C is the concentration and the subscripts o and f signify the initial and final state of the concentration of the waste cooking oil (WCO). When the model is subjected to higher molar ratio of 6:1 AL to WCO and at higher inlet mass flow results are depicted in Figure 9 of all the participating species. This figure illustrates the reaction takes place only within the inner tube ring. The conversion seems to be enhanced reaching as high as 89% based on the inlet and exit fraction evaluated in the inner loop of the tubular cylinder, however a considerable intermediate is still beside TG conversion and mandating the need for the flow to duel and additional

ð16Þ

penetration and insulated wall is also applied to all the bounded reactor walls.

No. of cells 281,382 427,742 534,228 Residence time (sec) 1.01 � 103 1.080 � <sup>10</sup><sup>3</sup> 1.100 � <sup>10</sup><sup>3</sup>

Error % 10% 2% —

Table 5. Cells for the two meshes for the mesh sensitivity analysis.

Figure 6. Discretized reactor domain using hybrid 3-D baseline mesh.


Table 4. Properties of the flow species.

Figure 7. Flow trajectory colored by the resident time for the (a) baseline, (b) fine and (c) very fine.

A total resident time of the order of 103 seconds is observed in each of the three meshes. If the flow was injected at relatively high velocity, it would ensure the required homogenize mixing of the reactants that also avoid any of the mass transfer limitation. The circumferential configuration of the entry and exit of the reactor enables it to maintain a long residence time, even at a higher inlet mass flow that forces the flow to move in a swirling trajectory.


Table 5. Cells for the two meshes for the mesh sensitivity analysis.

#### 5.3. Reactive flow analysis

The reactor tube is setup vertically and is subjected to inlet flow rate at Re = 6000 at the bottom of the tube into the inner reactor chamber. There, the two fluids of TG and methoxide (mixture of the methanol and the catalyst) are injected circumferentially at the stipulated molar/mass ratio and inlet temperature. The outflow is subjected to atmospheric pressure which is admitted at exit boundary condition of the outer tube chamber located at the top. The no-lip no penetration and insulated wall is also applied to all the bounded reactor walls.

To maintain single-phase flow avoiding the complication of other reactivation, the case of reactive flow simulation is carried out below the boiling point of methanol at a temperature of 60�C (333 K). A steady-state solution is sought for the flow. This is achieved by ignoring the temporal term of the governing equation. The flow is introduced to the reactor by means of an external peristaltic or diaphragm pump operating at relatively low head of nearly 2-m to overcome the viscous shear stresses and head losses and the reactor vertical head at an adjustable discharge capacity of up to 5 L/min.

#### 5.4. Simulation results

)

A total resident time of the order of 103 seconds is observed in each of the three meshes. If the flow was injected at relatively high velocity, it would ensure the required homogenize mixing of the reactants that also avoid any of the mass transfer limitation. The circumferential configuration of the entry and exit of the reactor enables it to maintain a long residence time, even at

Species Chemical formula Molecular weight (g/mol) Viscosity (kg/m.s) Cp (J/kg.C) Density (kg/m<sup>3</sup>

Methanol CH4O 32 3.96e-4 1.470e3 791.8 Waste oil C54H104O6 849 1.61e-2 2.2e3 883.3 Biodiesel C18H36O2 284 1.12e-3 1.187e3 870 Glycerol C3H8O3 92 1.412e0 0238.6 1261

a higher inlet mass flow that forces the flow to move in a swirling trajectory.

Figure 7. Flow trajectory colored by the resident time for the (a) baseline, (b) fine and (c) very fine.

Figure 6. Discretized reactor domain using hybrid 3-D baseline mesh.

Table 4. Properties of the flow species.

58 Biofuels - Challenges and opportunities

Based on the stoichiometry of transesterification reactive 3 moles of methanol are consumed with 1 mole of triglycerides for the making of 1 mole of biodiesel. Using this ratio(3:1 methanol to WCO molar), which also corresponds to AL/WCO mass fraction of 0.102/0.892, results on the species contour plots depicted in Figure 8 and the outlet fractions are summarized in Table 6.

At this velocity and molar ration, a low conversion of 28% is achieved. The conversion is described by Eq. (16) and is written as:

$$Conv = \frac{c\_{0,WCO} - c\_{f,WCO}}{c\_{0,WCO}} \times 100\tag{16}$$

where C is the concentration and the subscripts o and f signify the initial and final state of the concentration of the waste cooking oil (WCO). When the model is subjected to higher molar ratio of 6:1 AL to WCO and at higher inlet mass flow results are depicted in Figure 9 of all the participating species. This figure illustrates the reaction takes place only within the inner tube ring. The conversion seems to be enhanced reaching as high as 89% based on the inlet and exit fraction evaluated in the inner loop of the tubular cylinder, however a considerable intermediate is still beside TG conversion and mandating the need for the flow to duel and additional residence time to complete the reaction.

Figure 8. Molar fraction of species across the reactor (global scale).


glycerol. It must also be noted that these reactions are reversible and hence they can go either way, such that high conversion in the 1st tubular loop may be counter reacted by the second

Port & ratio/Species C54H106O6 (WCO) CH4O C18H36O2 (FAME) C3H8O3 Conversion

Chemical Kinetic and High Fidelity Modeling of Transesterification

http://dx.doi.org/10.5772/intechopen.80008

61

Inlet 0.25 0.75 0 0 0.66

Inlet 0.143 0.857 0 0 0.66

Inlet 0.10 0.90 0 0 0.66

Outlet 0.085 0.254 0.496 0.165

Outlet 0.049 0.574 0.283 0.094

Outlet 0.034 0.702 0.198 0.066

Nevertheless, several sensitivity studies considering AL/WCO ratio, temperature, speed or mass flux of the feedstock etc. can be conducted. Results of the species concentration subjected to higher methanol ratio and larger influx of the feedstock of 600 mL/min (Re = 6000) are listed in Table 7. This result clearly marked a reasonable conversion of 66% which is higher than those obtained at lower influx %. The conversion seems also insensitive for the additional amount of the methanol which sounds counterintuitive initially and defies the experimental observation; this trend may be explained by the used kinetics which are inherited from the experimental work obtained at one ratio. These obtained values failed to capture the influence of the increase of the concentration of specific reactant to steer the reaction forward. Therefore, this may suggest another testing procedure at another molar ratio in parallel, mimicking those obtained and shown in Figure 4. Though, results emphasize the advantage of increasing residence time and efficient mixing of the flow which reduces the

Transesterification of WCO to biodiesel is presented in this work and experimental and numerical analysis to this green process is presented. Initially, the chemical kinetics for the reaction are evaluated to the multiple reactions of the transesterification and in both the forward and reversal pathways for these reactions. The obtained kinetics presented by the activation energy (E) and the reaction pre-constant are evaluated and used to establish a high fidelity and robust reactive flow model. This model is based on computational fluid dynamics that is governed by the Navies-Stokes equations for isothermal, multiple species reactive flow

tubular loop and marking lower conversion.

Table 7. Species molar fraction at practical conditions.

3:1 methanol-WCO molar ratio

6:1 methanol-WCO molar ratio

9:1 methanol-WCO molar ratio

additional methanol concentration.

6. Conclusions

Table 6. Species molar fraction at idealistic conditions.

Figure 9. Molar fraction of the participating species within the first inner tube of the reactor and at 6:1 AL to WCO ratio.

From Figure 9, the concentration of the TG is reduced as the flow climbs up in the reactor as does the AL due to the production of biodiesel (E) and the two-intermediate species (MG and DG). The figure also clearly shows how DG precedes the formation of the MG and GL observed in the fading blue color contours of the DG concentration near the bottom inlet. The MG and GL show progressive color contours from the nil value represented by dark blue contours and reaching to higher value of 0.07 and 0.02 for each of MG and GL, respectively. Therefore, this allows additional residence time for the reaction in the second reactor tube reduces these intermediates and lead to higher production of the biodiesel and its byproduct


Table 7. Species molar fraction at practical conditions.

glycerol. It must also be noted that these reactions are reversible and hence they can go either way, such that high conversion in the 1st tubular loop may be counter reacted by the second tubular loop and marking lower conversion.

Nevertheless, several sensitivity studies considering AL/WCO ratio, temperature, speed or mass flux of the feedstock etc. can be conducted. Results of the species concentration subjected to higher methanol ratio and larger influx of the feedstock of 600 mL/min (Re = 6000) are listed in Table 7. This result clearly marked a reasonable conversion of 66% which is higher than those obtained at lower influx %. The conversion seems also insensitive for the additional amount of the methanol which sounds counterintuitive initially and defies the experimental observation; this trend may be explained by the used kinetics which are inherited from the experimental work obtained at one ratio. These obtained values failed to capture the influence of the increase of the concentration of specific reactant to steer the reaction forward. Therefore, this may suggest another testing procedure at another molar ratio in parallel, mimicking those obtained and shown in Figure 4. Though, results emphasize the advantage of increasing residence time and efficient mixing of the flow which reduces the additional methanol concentration.

### 6. Conclusions

From Figure 9, the concentration of the TG is reduced as the flow climbs up in the reactor as does the AL due to the production of biodiesel (E) and the two-intermediate species (MG and DG). The figure also clearly shows how DG precedes the formation of the MG and GL observed in the fading blue color contours of the DG concentration near the bottom inlet. The MG and GL show progressive color contours from the nil value represented by dark blue contours and reaching to higher value of 0.07 and 0.02 for each of MG and GL, respectively. Therefore, this allows additional residence time for the reaction in the second reactor tube reduces these intermediates and lead to higher production of the biodiesel and its byproduct

Figure 9. Molar fraction of the participating species within the first inner tube of the reactor and at 6:1 AL to WCO ratio.

Port & ratio/species C54H104O6 (WCO) CH4O C18H36O2 (FAME) C3H8O3 Conversion Inlet 0.25 0.75 0 0 28%

Outlet 0.18\* 0.54 0.21 0.07

Figure 8. Molar fraction of species across the reactor (global scale).

Table 6. Species molar fraction at idealistic conditions.

\*

Conversion = (0.25�0.18)/0.25.

60 Biofuels - Challenges and opportunities

Transesterification of WCO to biodiesel is presented in this work and experimental and numerical analysis to this green process is presented. Initially, the chemical kinetics for the reaction are evaluated to the multiple reactions of the transesterification and in both the forward and reversal pathways for these reactions. The obtained kinetics presented by the activation energy (E) and the reaction pre-constant are evaluated and used to establish a high fidelity and robust reactive flow model. This model is based on computational fluid dynamics that is governed by the Navies-Stokes equations for isothermal, multiple species reactive flow in turbulent regime. The model is built around the newly patented multiple tube and continuous transesterification reactor. The flow enters the reactor laterally and induces swirling flow. In turn, this results in an order of magnitude residence time higher than the transfer time that is based on the tube length and inlet velocity. Reactive flow at low velocity results in low conversion of the WCO, while actual turbulent flow significantly increases the conversion rate. Excess methanol mass sensitivity was insignificant to the conversion parameter that suggest that modeling can be still limited and suggesting to integrate more kinetics data. This work emphasizes the efficient mixing of the flow, interplaying as a parameter to the additional methanol concentration, thereby avoiding their downstream separation. In closing, numerical simulation of transesterification undoubtedly demonstrates the effectiveness of this tool in analyzing complicated reactive flow as the student faces in transesterification. Still, the process fails to reach completion even at higher reactant concentration and beyond the process stoichiometry. This creates an opportunity in this area to go beyond the conventional methods, such as electrical stimulation or signification and to assess their needed kinetics accordingly. The results demonstrate the feasibility of reactive flow dynamics in capturing and numerically simulating the transesterification process. This work substantiates the practicability of using numerical methods to construct precise and insightful image of the distribution of reaction rates and associated species to design and develop more efficient reactors.

[2] Banga S, Varshney PK. Effect of impurities on performance of biodiesel: A review. Journal

Chemical Kinetic and High Fidelity Modeling of Transesterification

http://dx.doi.org/10.5772/intechopen.80008

63

[3] Li Z, Ramanathan S, Chen C-C. Predicting thermophysical properties of mono-and diglycerides with the chemical constituent fragment approach. Industrial & Engineering Chem-

[4] Ceriani R, Goncalves CB, Coutinho JAP. Prediction of viscosities of fatty compounds and

[5] Noureddini H, Zhu D. Kinetics of transesterification of soybean oil. Journal of the Amer-

[6] Janajreh I, ElSamad T, AlJaberi A, Diouri M. Transesterification of waste cooking oil:

[7] Janajreh I, Al Shrah M. Numerical simulation of multiple step transesterification of waste

[8] Janajreh I, Almusharekh M, Ghenai C. Transesterification process of waste cooking oil: Catalyst synthesis, kinetic study, and modeling sensitivity. In: International Conference on

[9] Janajreh I, ElSamad T. Mohammed Noorul Hussain, intensification of transesterification via sonication numerical simulation and sensitivity study. Applied Energy. 2017;185:2151-

[10] Janajreh I, Al Musharrekh M, Hussain M. Numerical simulation of transesterification and sensitivity study. International Journal of Modern Engineering. 2015;15(2):13-22

[11] Laidler KJ. The development of the Arrhenius equation. Journal of Chemical Education.

[12] Abd Rabu R, Janajreh I, Honnery D. Transesterification of waste cooking oil: Process optimization and conversion rate evaluation. Energy Conversion and Management.

[14] Freedman B, Butterfield RO, Pryde EH. Transeterfication kinetics of soybean oil. Journal of

[15] Freedman B, Pryde EH, Mounts TL. Variables affecting the yield of fatty esters from

[16] Farkas L, Schachter O, Vromen BH. On the rate of acetate. Journal of the American

[17] Dufeck EJ, Butterfield RO, Frankel EN. Esterification and transsesterification of 9(10) caboxystearic acid and its methyl esters, kinetic studies. Journal of the American Oil

[13] MATLAB. version 8.3.0 (R2014a). Natick, Massachusetts: The Math Works Inc.; 2014

the American Oil Chemists' Society. 1986;63:1375-1380

transesterfified vegetable oils. IBID. 1984;61:1638-1643

Chemical Society. 1949;71:1991-1994

Chemists' Society. 1972;49:302-306

biodiesel by group contribution. Energy & Fuels. 2011;25(8):3712-3717

Kinetic study and reactive flow analysis. Energy Procedia. 2015;75:547-553

oil in tubular reactor. Journal of Infrastructure Systems. 2016;22(4):A4014007

of Scientific and Industrial Research. 2010;69(8):575-579

ican Oil Chemists' Society. 1997;74(11):1457-1463

Sustainable Solid Waste Management. 2014

2159

1984;61(6):494

2013;65:764-769

istry Research. 2010;49(11):5479-5484

#### Acknowledgements

Khalifa University of Science and Technology, Masdar Institute Campus is highly acknowledged for their support.
