Abstract
In this article, the occurrence of dead core in catalytic particles containing immobilized enzymes is analyzed for the Michaelis–Menten kinetics. An assessment of numerical methods is performed to solve the boundary value problem generated by the mathematical modeling of diffusion and reaction processes under steady state and isothermal conditions. Two classes of numerical methods were employed: shooting and collocation. The shooting method used the ode function from Scilab software. The collocation methods included: that implemented by the bvode function of Scilab, the orthogonal collocation, and the orthogonal collocation on finite elements. The methods were validated for simplified forms of the Michaelis–Menten equation (zero-order and first-order kinetics), for which analytical solutions are available. Among the methods covered in this article, the orthogonal collocation on finite elements proved to be the most robust and efficient method to solve the boundary value problem concerning Michaelis–Menten kinetics. For this enzyme kinetics, it was found that the dead core can occur when verified certain conditions of diffusion–reaction within the catalytic particle. The application of the concepts and methods presented in this study will allow for a more generalized analysis and more accurate designs of heterogeneous enzymatic reactors.
Similar content being viewed by others
References
Šekuljica NŽ, Prlainović NŽ, Jovanović JR, Stefanović AB, Djokić VR, Dušan ŽM, Knežević-Jugović ZD (2016) Immobilization of horseradish peroxidase onto kaolin. Bioproc Biosyst Eng. doi:10.1007/s00449-015-1529-x
Edet E, Ntekpe M, Omereji S (2013) Current trend in enzyme immobilization: a review. Int J Mod Biochem 2:31–49
Illanes A (2008) Enzyme biocatalysis: principles and applications. Springer, Netherlands
Bailey JE, Ollis DF (1986) Biochemical engineering fundamentals. McGraw-Hill, New York
Doran PM (2013) Bioprocess engineering principles. Academic Press, Waltham
Fogler HS (2005) Elements of chemical reaction engineering. Prentice Hall PTR, Upper Saddle River
Oliveira SC (1999) Evaluation of effectiveness factor for immobilized enzymes using Runge–Kutta–Gill method: how to solve mathematical undetermination at particle center point? Bioproc Eng. doi:10.1007/s004490050579
Granato MA, Queiroz LC (2003) Dead core in porous catalysts: modeling and simulation of a case problem using Mathematica. In: Proceedings of the workshop modelling and simulation in chemical engineering, CIM–Centro Internacional de Matemática, Coimbra, 30 June–4 July 2003
Bird RB, Stewart WE, Lightfoot EN (2002) Transport phenomena. Wiley, New York
Froment GF, Bischoff KB, De Wilde J (2011) Chemical reactor analysis and design. Wiley, Hoboken
Mahalakshmi M, Hariharan G (2016) An efficient Chebyshev wavelet based analytical algorithm to steady state reaction–diffusion models arising in mathematical chemistry. J Math Chem. doi:10.1007/s10910-015-0560-0
Devi MR, Sevukaperumal S, Rajendran L (2015) Non-linear reaction diffusion equation with Michaelis–Menten kinetics and Adomian decomposition method. Appl Math. doi:10.5923/j.am.20150501.04
Shanthi D, Ananthaswamy V, Rajendran L (2013) Analysis of non-linear reaction-diffusion processes with Michaelis–Menten kinetics by a new Homotopy perturbation method. Nat Sci. doi:10.4236/ns.2013.59128
Joy RA, Meena A, Loghambal S, Rajendran L (2011) A two-parameter mathematical model for immobilized enzymes and Homotopy analysis method. Nat Sci. doi:10.4236/ns.2011.37078
Hindmarsh AC (1983) In: Stepleman RS et al (eds) Scientific computing. IMACS Transactions on Scientific Computation, Amsterdam
Brown PN, Hindmarsh AC (1989) Reduced storage matrix methods in stiff ODE systems. Appl Math Comput. doi:10.1016/0096-3003(89)90110-0
Description and use of LSODE, the Livermore solver for ordinary differential equations (1993) NASA, Livermore. https://computation.llnl.gov/casc/nsde/pubs/u113855.pdf. Accessed 16 Feb 2016
Shampine LF (1996) Some practical Runge–Kutta formulas. Math Comput. doi:10.2307/2008219
Shampine LF, Watts HA (1980) DEPAC–design of a user oriented package of ODE solvers. Sandia National Laboratories, Albuquerque
Fehlberg E (1970) Klassische Runge–Kutta-Formeln vierter und niedrigerer Ordnung mit Schrittweiten-Kontrolle und ihre Anwendung auf Wärmeleitungs probleme. Computing. doi:10.1007/BF02241732
Ascher U, Christiansen J, Russell RD (1981) Collocation software for boundary-value ODEs. ACM T Math Softw. doi:10.1145/355945.355950
Bader G, Ascher U (1987) A new basis implementation for a mixed order boundary value ode solver. SIAM J Sci Stat Comp. doi:10.1137/0908047
Ascher U, Christiansen J, Russell RD (1979) A collocation solver for mixed order systems of boundary value problems. Math Comput. doi:10.1090/S0025-5718-1979-0521281-7
Deboor C, Weiss R (1980) SOLVEBLOK: a package for solving almost block diagonal linear systems. ACM T Math Softw. doi:10.1145/355873.355880
Villadsen J, Michelsen ML (1978) Solution of differential equation models by polynomial approximation. Prentice-Hall, Englewood Cliffs
Wylie CR, Barret LC (1982) In: Wylie CR, Barret LC (eds) Advanced engineering mathematics. McGraw-Hill, Singapore
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Pereira, F.M., Oliveira, S.C. Occurrence of dead core in catalytic particles containing immobilized enzymes: analysis for the Michaelis–Menten kinetics and assessment of numerical methods. Bioprocess Biosyst Eng 39, 1717–1727 (2016). https://doi.org/10.1007/s00449-016-1647-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00449-016-1647-0