Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 810843, 8 pages
Research Article

A Finite-Difference Solution of Solute Transport through a Membrane Bioreactor

1Department of Chemical Engineering, Cape Peninsula University of Technology, P.O. Box 652, Cape Town 8000, South Africa
2Department of Mathematics and Applied Mathematics, University of Cape Town, Private Bag X3, Rondebosch 7700, South Africa

Received 26 December 2014; Revised 10 March 2015; Accepted 12 March 2015

Academic Editor: Sergio Preidikman

Copyright © 2015 B. Godongwana et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The current paper presents a theoretical analysis of the transport of solutes through a fixed-film membrane bioreactor (MBR), immobilised with an active biocatalyst. The dimensionless convection-diffusion equation with variable coefficients was solved analytically and numerically for concentration profiles of the solutes through the MBR. The analytical solution makes use of regular perturbation and accounts for radial convective flow as well as axial diffusion of the substrate species. The Michaelis-Menten (or Monod) rate equation was assumed for the sink term, and the perturbation was extended up to second-order. In the analytical solution only the first-order limit of the Michaelis-Menten equation was considered; hence the linearized equation was solved. In the numerical solution, however, this restriction was lifted. The solution of the nonlinear, elliptic, partial differential equation was based on an implicit finite-difference method (FDM). An upwind scheme was employed for numerical stability. The resulting algebraic equations were solved simultaneously using the multivariate Newton-Raphson iteration method. The solution allows for the evaluation of the effect on the concentration profiles of (i) the radial and axial convective velocity, (ii) the convective mass transfer rates, (iii) the reaction rates, (iv) the fraction retentate, and (v) the aspect ratio.