Mathematical Modeling of the Transesterification Reaction by Finite Elements: Optimization of Kinetic Parameters Using the Simplex Sequential Method

The use of biodiesel as fuel has been presented as a viable alternative in the search for renewable energies. It can be produced from the transesterification reaction of vegetable oils with methanol in the presence of sodium or potassium hydroxide as catalysts. In the present research, the transesterification reaction of soybean oil was modeled considering the three steps of the direct and reverse reactions following a second order general kinetics by the finite element method using the COMSOL Multiphysics software. The values of rate constants were determined using the simplex optimization method coupled with the desirability functions. The optimized rate constants for the forward reactions were 0.250 × 10 (k1f), 1.137 × 10 (k2f) and 3.134 × 10 (k3f); and for the reverse reactions were 0.202 × 10 (k1r), 0.884 × 10 (k2r) and 0.219 × 10 (k3r) all in m mol s. The kinetic model proposed for the reaction can be simulated by the finite element method (FEM) under realistic conditions.