Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2015, Article ID 810843, 8 pages http://dx.doi.org/10.1155/2015/810843
Research Article A Finite-Difference Solution of Solute Transport through a Membrane Bioreactor B. Godongwana,1 D. Solomons,2 and M. S. Sheldon1 1
Department of Chemical Engineering, Cape Peninsula University of Technology, P.O. Box 652, Cape Town 8000, South Africa Department of Mathematics and Applied Mathematics, University of Cape Town, Private Bag X3, Rondebosch 7700, South Africa
2
Correspondence should be addressed to B. Godongwana;
[email protected] 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.
1. Introduction Membrane bioreactors (MBRs) are finding increasing use in the production of primary and secondary metabolites such as amino acids, antibiotics, anticancer drugs, and tissue cells [1β3]. This technology is favoured by recent trends towards environmentally-friendly technologies, particularly because MBRs do not require additives, function at moderate operating conditions, and reduce by-product formation [1]. The efficiency of MBRs is dependent mainly on the transport of solutes through the bioreactor, and this is influenced by biochemical, geometric, and hydrodynamic parameters [2, 4]. This paper considers the numerical solution of the convection-diffusion equation for solute transport through a fixedfilm MBR. This analysis is important for simulation of the performance (i.e., efficiency and effectiveness) of the bioreactor. The governing equation for mass transport of solutes through the bioreactor is the convection-diffusion equation with Monod kinetics [5] as follows: π’
2
π π 1 π ππ ππ ππ ππ (π ) + 2 ] β π , + V = π·π΄π΅ [ ππ§ ππ π ππ ππ ππ§ πΎπ + π
(1)
where π is the local substrate concentration, π’ and V are the axial and radial velocity components, respectively, π·π΄π΅ is the substrate diffusion coefficient, ππ is the maximum rate of reaction, and πΎπ is the saturation (or Michaelis) constant. Equation (1) is made dimensionless by introducing the following variables: π=
π’ ; π’0
π=
πΎ β πΎπ = π; π0
V ; V0
πΆ=
π§ π= ; πΏ
π ; π0
π ; π
= π
1
π=β
πππ
12 , π0 π·π΄π΅
(2)
π
π = 1. πΏ
Equation (1) then becomes πPeπ’ π
ππΆ ππΆ + PeV π ππ ππ
π2 πΆ π2 πΆ 1 π ππΆ =[ (π
) + π2 2 ] β β , π
ππ
ππ
ππ πΎπ + πΆ
(3)
2
Mathematical Problems in Engineering
where the axial and radial Peclet numbers (Peπ’,V ) are, respectively, defined as Peπ’ =
π’0 π
1 ; π·π΄π΅
PeV =
V0 π
1 . π·π΄π΅
(4)
The velocity profiles, π and π, in (3) are solutions of the π§ and π-components of the Navier-Stokes equations, respectively [6]: 1 ππ Re π = β (1 β π
2 ) ( β ), 4 ππ Fr π’ π
π
2 π2 π , π = π ( 0 ) [ (1 β )] V0 8 2 ππ2
(5)
where π is the dimensionless hydrostatic pressure which is a function of the membrane hydraulic permeability π
and Re and Fr are the Reynolds and Froude numbers, respectively. When the membrane hydraulic permeability π
is much smaller than unity, (3) reduces to πβ
ππΆ ππΆ + Peπ’ π
π½ (2π
β π
3 ) ππ ππ
π2 πΆ π2 πΆ 1 π ππΆ =[ (π
) + π2 2 ] β β , π
ππ
ππ
ππ πΎπ + πΆ
(6)
where πβ = β4Peπ’ π
π½ [
1 + π] (1 β π
2 ) , (π β 1)
π =ΜΈ 1.
(7)
The fraction retentate, π, is defined as the ratio of the outlet to the inlet axial velocity (π = 0 for the dead-end mode and π βΌ 1 for the closed-shell mode), and π½ is the dimensionless transmembrane pressure. The corresponding boundary conditions are B.C.1
at π = 0 βπ
πΆ = 1,
B.C.2
at π
= 0 βπ
B.C.3
ππΆ 2 ππΆ at π
= 1 βπ = . ππ πPeπ’ ππ
ππΆ = 0, ππ
and externally unskinned region of microvoids. The nutrient solution is supplied by a peristaltic pump and permeates from the lumen-side to the shell-side of the MGR. The microorganism is immobilised on either the lumen-side or the shell-side of the MGR. Humidified air is supplied on the shell-side, and two pressure transducers are fitted at the inlet and outlet of the MGR.
2. Analytical Models The Graetz problem [8] is one of the oldest forced-convection problems describing the steady-temperature distribution and rate of heat transfer in tube flow. The evaluation of (6) for concentration profiles in a tubular reactor is mathematically analogous to the Graetz problem [5, 9]. In the original Graetz problem, however, there is no reaction (or source) term and axial diffusion, and radial convection is ignored. The assumption of negligible axial diffusion and radial convection is common in the majority of analytical models currently in use [7, 10β12]. Radial convective flows have been shown to significantly improve MBR efficiency [4, 13β15]. In the dead-end ultrafiltration mode, particularly, the assumption of negligible radial convective flow is not justifiable. Nagy [15] investigated the effect of radial convective flows (PeV ) on the mass transfer rates of solutes through a biocatalytic membrane layer. Analytical solutions of (6) for the zero-order and first-order limits of the Monod equation were provided. This analysis, however, was restricted to the matrix/fiber region of the membrane and hence the radial velocity was assumed constant, and axial convective and diffusive flows were ignored. At axial Peclet numbers smaller than one (Peπ’ βͺ 1) large concentration gradients exist in the membrane lumen, and ignoring axial diffusion is also not justified [6, 7]. The model proposed by Godongwana et al. [6] follows the approach suggested by Davis [9], that is, writing the solution of (6) in terms of known functions. The model accounts for radial convective flow and axial diffusion for the limiting case of first-order kinetics. In that model, (6) was solved by separation of variables and regular perturbation, resulting in the asymptotic expansion:
(8)
Boundary condition 1 (B.C.1) corresponds to a uniform inlet substrate concentration; B.C.2 corresponds to cylindrical symmetry at the center of the membrane lumen; B.C.3 corresponds to continuity of the substrate flux at the lumen-matrix interface [7]. The solution of (6) is based on the following general assumptions: (i) the system is isothermal; (ii) the flow regime is laminar and fully developed; (iii) the fluid is Newtonian and homogenous and has constant physical and transport properties; (iv) the membrane hydraulic permeability is constant. A schematic of the MBR is shown in Figure 1. The MBR consists of a single hollow-fibre, made of surface modified polysulphone, encased in a glass bioreactor. The membranes are asymmetric and characterized by an internally skinned
β
π
πΆ (π, π₯) = β β π΅π πΉπ (π) ππ (π₯) π
π .
(9)
π=1 π=0
Making use of the following changes of variables: π=β
2Peπ’ π
π½ 1 [ + π] , π2 (π β 1)
π = β(
π2 ) π2 , 4Peπ’ π
π½
(10)
π₯ = π π π
, where πΉ(π) in (9) is represented by the Kummer function [16] as follows: β ) 1 β (π2π + π2 /πΎπ πΉπ (π) = π [ , , π] . 4Peπ’ π
π½ 2
(11)
Mathematical Problems in Engineering
3
Outlet pressure transducer Humidified air
Glass bioreactor Capillary membrane
Permeate collection and air
Inlet pressure transducer
Nutrient solution
Peristaltic pump
Figure 1: A schematic diagram of the single capillary membrane bioreactor.
The zero-order and first-order approximations of π(π₯) in (9) are, respectively, π0 (π₯) = π½0 (π₯) , π1 (π₯) = π1 [
(π₯)3 π½3 (π₯) (π₯)2 π½2 (π₯) (π₯)4 π½4 (π₯) + π2 + π3 ], 3!! 5!! 7!! (12)
where π π are the eigenvalues and π½π is the Bessel function of the first kind of order π [16]: π1 = β
3PeV π½ (π’0 /V0 ) , 2π2π
π2 = β
20 , 3π2π
π3 =
35 . 4π2π (13)
The eigenvalues are obtained from B.C.3 of (8) and are roots of the following equation: β β (π2π + π2 /πΎπ ) π2 ππ 3 2 (π π + β ) π [ + 1, , π] 2π
π½ πΎπ 4Peπ’ π
π½ 2 β β (π2π + π2 /πΎπ ) 1 π½ (π) = 2π π 1 π[ , , π] . π½0 (π) 4Peπ’ π
π½ 2
the inlet condition B.C.1 of (8) and employing Lommel integrals to give [6] π΅π =
2 π π π [β (π2π
+
β ) /4Pe π
π½, 1/2, π ] π2 /πΎπ π’ 0
π½1 (π π ) β
[ 2 ]. π½0 (π π ) + π½12 (π π )
(15)
β β« πΆ), used in The assumption of first-order kinetics (πΎπ the above analysis, allows for analytical evaluation of (6); however it limits the range of inlet substrate concentrations. The complete nonlinear form of (6) is not amenable to analytical evaluation and hence was solved using a finite-difference scheme described in Section 3.
3. Finite-Difference Scheme A finite-difference representation of (6) is obtained by employing first-order upwind difference quotients for the derivatives on the LHS:
(14)
The first ten eigenvalues, at different axial positions, are listed in the appendix. The coefficient π΅π is obtained by imposing
(
πΆπ,π β πΆπβ1,π ππΆ ) = + π (β) , ππ
π,π β
πΆπ,π β πΆπ,πβ1 ππΆ + π (π) . ( ) = ππ π,π π
(16)
4
Mathematical Problems in Engineering
And second-order central differences for the derivatives on the RHS of (6) are (
πΆπ+1,π β 2πΆπ,π + πΆπβ1,π π2 πΆ ) = + π (β)2 , 2 β2 ππ
π,π
πΆπ,π+1 β 2πΆπ,π + πΆπ,πβ1 π2 πΆ ( 2) = + π (π)2 , π2 ππ π,π
(17)
resulting in the general difference equation as follows: πΌ1 πΆπβ1,π + (πΌ2 +
β πΎπ
πΌ3 ) πΆπ,π + πΌ4 πΆπ+1,π + πΆπ.π
(18)
+ πΌ5 πΆπ,πβ1 + πΌ6 πΆπ,π+1 = 0,
2
2
β β β ) πβ + βPeπ’ π
π½ (2π
β π
3 ) + 2π2 ( 2 ) + 2 β , π π π
2 2
πΌ3 = β π , πΌ4 = β1, πβ π2 + 2), π π
πΌ6 = βπ2 (
β2 ). π2 (19)
The corresponding difference quotients of the boundary conditions in (8) are βPeV ππ
B.C.1
at π = 1 βπ πΆπ,π = π
B.C.2
at π = 1 βπ
B.C.3
at π = π βπ
πΆπ,π β πΆπβ1,π β πΆπ+1,π β πΆπ,π β
,
= 0, =
πPeπ’ πΆπ,π β πΆπ,πβ1 , 2 π
(20)
πΆπ,π β πΆπ,πβ1
= 0. π Boundary condition 1 is an approximation of the inlet condition in (8) and is obtained by assuming the net diffusive and convective flux in the radial direction are negligible at the membrane entrance. Boundary condition 4 is applicable in the dead-end mode (one end of the lumen-side is closed). The solution domain is a regular 2-dimensional grid and is subdivided into π-intervals (of size β) in the π-dimension and π-intervals (of size π) in the π§-dimension. The difference equation (18), including the boundary conditions of (20), is solved by making use of the multivariate Newton-Rapshon iteration scheme: B.C.4
at π = π βπ
C(π+1) = Cπ β
π΅
0
π·
π΅
πΈ .. .
π·
0
0
0
0
0
0
β
β
β
0
0
0
0
0
0
0
0 .. .
d β
β
β
π·
π΅
πΈ
π·
0
πΈ
] ] ] 0 ] ] ] ] ]. ] ] 0 ] ] π΅ ] ] π· ]
(22)
π΅ = πΌ4 + πΌ5 ,
β πΌ1 = β [βPeπ’ π
π½ (2π
β π
3 ) + 1 β ] , π
πΌ5 = ββ2 (
π· [ [ πΈ [ [ 0 [ [ [ J=[ [ [ [ 0 [ [ 0 [ [ 0
The matrix elements π΅, π·, and πΈ are, respectively,
where
πΌ2 = (
where F(C) is the residual equation (18) and J is the tridiagonal Jacobian matrix as follows:
F (Cπ ) , Jπ
(21)
πΆπ.π 1 ], + π· = πΌ2 + πΌ3 [ β πΎπ + πΆπ.π (πΎβ + πΆ )2 π.π ] π [
(23)
πΈ = πΌ1 + πΌ6 . The Newton-Raphson iteration scheme was implemented on MATLAB R2014a and the procedure is shown in Figure 2. The algorithm begins with an initial guess of the solute concentration at each grid point; an initial guess of zero was used. Based on this guess, the residual column vector and the Jacobian matrix can be evaluated. The magnitude (Euclidean norm) of the quotient of the residual vector and Jacobian matrix, πC, is evaluated. The iteration is repeated with new solute concentration guess values until the Euclidean norm is less than the prescribed tolerance.
4. Results 4.1. Numerical Solution. The implicit finite-difference scheme was shown to be unconditionally stable for the different values of β and π listed in Table 1. The results are shown in Figure 3 for the parameter values listed in Table 2. Figure 3 illustrates the effect of the fraction retentate π on the solute concentration profiles. In the dead-end mode (π = 0) there is increased radial convective flow as shown by the streamlines in Figure 3(a). This increased radial flow allows for more solute contact with the biofilm and hence improved conversion, resulting in higher MBR efficiency. In this mode however the solute is limited to only the entrance half of the MBR as shown in Figure 3(a). Increasing the fraction retentate to π = 0.8 allows for a uniform distribution of the solute (Figure 3(b)); however radial convective flow is significantly reduced. This result implies that an optimum π value should be sought for enhanced MBR efficiency. The developed finitedifference scheme also allows for the evaluation of the effect on the concentration profiles of the radial and axial convective velocity, the convective mass transfer rates, the reaction rates, and the aspect ratio. The sensitivity analysis of these parameters however has been omitted in the current paper. Figure 4 is a plot of the solute concentration profile, from the FDM scheme, for different dimensionless saturation
Mathematical Problems in Engineering
5
Input physical parameters
Guess values of C for the m Γ n nodes
Final results
Yes
Loop through the nodes
Build the residual column vector F(C) from Eq. (18)
No
dπ β€ 0.0001
B.C., Eq. (20)
Build the Jacobian matrix J from Eq. (22)-(23)
Solve for C from Eq. (21)
Find the Euclidean norm (or l2 norm) of dπ, where dπ = π
(π)/π
Figure 2: The Newton-Raphson algorithm for solving (18). Table 1: Computation times for different sizes of β and π. Simulation 1 2 3
Number of radial nodes (π) 32 64 128
Number of axial nodes (π) 32 64 128
Total number of solution nodes 1,024 4,096 16,384
Total iteration time (s) 13.302 55.130 505.108
Table 2: Parameter values used to determine the concentration profile [6]. Model parameter Membrane hydraulic permeability Membrane inner radius Effective membrane length Lumen-side entrance axial velocity Permeation velocity Lumen-side inlet hydrostatic pressure Shell-side hydrostatic pressure Glucose diffusivity Solution density Solution viscosity Glucose inlet concentration Kinetic constants
Symbol ππ π
1 πΏ π’0 V0 π0 ππ π·π΄π΅ π π π0 πΎπ /ππ
Unit m/Pas m m m sβ1 m sβ1 Pa Pa m2 s kg mβ3 Pas g dmβ3 sβ1
Basic measured value 3.82 Γ 10β11 1.30 Γ 10β4 5.7 Γ 10β2 1.67 Γ 10β3 1.91 Γ 10β7 106 325 101 325 1.0 Γ 10β10 998.0 9.7 Γ 10β4 2.00 0.10
Mathematical Problems in Engineering 1
1
0.9
0.9
0.8
0.8
0.7
0.7
0.6
0.6
Z = z/L
Z = z/L
6
0.5 0.4
0.5 0.4
0.3
0.3
0.2
0.2
0.1
0.1
0 0
0.2
0.4
0.6
0.8
1
0 0
0.2
R = r/R1
0.4
0.6
0.8
1
R = r/R1
Concentration Streamlines
Concentration Streamlines (a)
(b)
Figure 3: Solute concentration profiles for π = π = 64 when (a) π = 0 and (b) π = 0.8 (for the parameter values listed in Table 2).
1.000
C = c/c0
0.800
0.600
0.400
0.200
0.000 0.0
0.2
0.4 0.6 Z = z/L
0.8
1.0
FDM (Km = 0.009) FDM (Km = 0.09) FDM (Km = 0.9)
Figure 4: Solute concentration profiles from the finite-difference method (FDM) for different saturation constants πΎπ , when π = 0.8. β constants when π = 0.8. At a low πΎπ value of 0.009 the dimensionless outlet solute concentration is 0.472, whereas β to 0.9 reduces the outlet concentration to increasing πΎπ β value only 0.853. This result is expected since a lower πΎπ is consistent with a higher biofilm/enzyme affinity for the substrate [17]. This result also illustrates the significance of the appropriate choice of substrate on bioreactor design.
4.2. Comparison with Analytical Solution. In Figure 5(a) the FDM scheme is compared with the analytical model presented in Section 2 for the open-shell mode (π = 0.8).
The analytical model predicts a linear decrease in the solute concentration inside the membrane lumen to 47% of the original concentration. This result is consistent, qualitatively, with the result of Heath and Belfort [7] for the parameter values listed in Table 2. The FDM scheme predicts the same outlet concentration; however the decrease is gradually close to the entrance and rises with increasing length. The discrepancies between the two profiles arise from the assumption of firstorder kinetics, assumed in developing the analytical solution. These two profiles suggest that the open-shell mode is suitable for microbial growth since the substrate is not depleted inside the lumen. The rapid decline in the solute concentration in Figure 5(b) is due to increased radial convective flow in the dead-end mode. This results in nonuniform microbial growth/tapering as observed by Godongwana et al. [18] for the bacterium Streptomyces coelicolor on a ceramic membrane. This phenomenon can be reduced by either increasing the solute flow rate or increasing the fraction retentate π. The numerical scheme matches the analytical model approximately on a small interval close to the origin. The divergence again is attributed to the assumption of first-order kinetics.
5. Conclusion A numerical solution of the dimensionless convectiondiffusion equation, with nonlinear kinetics, was developed. The numerical scheme was performed using the NewtonRaphson method and was shown to be unconditionally stable for different step-sizes (β and π). The analysis provides for evaluation of concentration profiles of solutes through a membrane bioreactor. The numerical solution was compared to a regular perturbation solution for two modes of operation, that is, the dead-end mode and open-shell mode. In
7
1.00
1.000
0.80
0.800
0.60
0.600
C = c/c0
C = c/c0
Mathematical Problems in Engineering
0.40
0.400 0.200
0.20
0.000
0.00 0
0.2
0.4
0.6
0.8
0
1
0.2
0.4
0.6
0.8
1
Z = z/L
Z = z/L Exact
Exact FDM
FDM (a)
(b)
Figure 5: A comparison of the analytical versus the FDM solution for solute concentration profiles: (a) π = 0.8 and (b) π = 0 (for the parameter values listed in Table 2). Table 3: Eigenvalues. (a) Positive roots of (14), π = 0
π1 π2 π3 π4 π5 π6 π7 π8 π9 π 10
π=0 2.142 4.968 7.891 10.885 13.929 17.004 20.101 23.211 26.331 29.458
π = 0.2 2.191 5.060 8.000 10.997 14.037 17.107 20.198 23.305 26.422 29.547
π = 0.4 2.242 5.162 8.131 11.140 14.181 17.249 20.337 23.440 26.554 29.677
π = 0.6 2.294 5.273 8.285 11.320 14.376 17.450 20.541 23.644 26.757 29.88
π = 0.8 2.349 5.394 8.462 11.541 14.630 17.728 20.835 23.949 27.070 30.197
π = 1.0 3.832 7.016 10.173 13.324 16.471 19.616 22.760 25.904 29.047 32.190
(b) Positive roots of (14), π = 0.8
π1 π2 π3 π4 π5 π6 π7 π8 π9 π 10
π=0 1.518 4.234 7.256 10.347 13.463 16.589 19.721 22.857 25.995 29.134
π = 0.2 1.538 4.248 7.265 10.354 13.468 16.594 19.726 22.861 25.999 29.138
π = 0.4 1.559 4.264 7.275 10.362 13.475 16.599 19.730 22.865 26.003 29.142
the dead-end mode the numerical results closely matched the perturbation solution. The assumption of linear kinetics, commonly used in literature models, was shown to result in inaccuracies in the open-shell mode. The numerical solution allows for the evaluation of the influence of the general operating parameters of a MBR on the concentration profiles.
π = 0.6 1.581 4.280 7.286 10.370 13.481 16.605 19.736 22.870 26.007 29.146
π = 0.8 1.604 4.298 7.299 10.379 13.489 16.611 19.741 22.875 26.012 29.151
π = 1.0 1.628 4.318 7.312 10.389 13.497 16.618 19.748 22.881 26.017 29.156
The fraction retentate (π) was shown to be an important optimisation parameter for improved MBR efficiency.
Appendix See Table 3.
8
Mathematical Problems in Engineering
Nomenclature
References
π΅π :
[1] L. Giorno and E. Drioli, βBiocatalytic membrane reactors: applications and perspectives,β Trends in Biotechnology, vol. 18, no. 8, pp. 339β349, 2000. [2] C. Charcosset, βMembrane processes in biotechnology: an overview,β Biotechnology Advances, vol. 24, no. 5, pp. 482β492, 2006. [3] D. F. Stamatialis, B. J. Papenburg, M. GironΒ΄es et al., βMedical applications of membranes: drug delivery, artificial organs and tissue engineering,β Journal of Membrane Science, vol. 308, no. 1-2, pp. 1β34, 2008. [4] S. Curcio, V. Calabr`o, and G. Iorio, βA theoretical and experimental analysis of a membrane bioreactor performance in recycle configuration,β Journal of Membrane Science, vol. 273, no. 1-2, pp. 129β142, 2006. [5] B. R. Bird, W. E. Stewart, and E. N. Lightfoot, Transport Phenomena, John Wiley & Sons, New York, NY, USA, 2nd edition, 2002. [6] B. Godongwana, D. Solomons, and M. S. Sheldon, βA solution of the convective-diffusion equation for solute mass transfer inside a capillary membrane bioreactor,β International Journal of Chemical Engineering, vol. 2010, Article ID 738482, 12 pages, 2010. [7] C. Heath and G. Belfort, βImmobilization of suspended mammalian cells: analysis of hollow fiber and microcapsule bioreactors,β Advances in Biochemical Engineering/Biotechnology, vol. 34, pp. 1β31, 1987. Β¨ [8] L. Graetz, βUber die wΒ¨armeleitungsfΒ¨ahigkeit von flΒ¨ussigkeiten,β Annalen der Physik und Chemie, vol. 18, pp. 79β94, 1883. [9] E. J. Davis, βExact solutions for a class of heat and mass transfer problems,β The Canadian Journal of Chemical Engineering, vol. 51, no. 5, pp. 562β572, 1973. [10] L. R. Waterland, A. S. Michaels, and C. R. Robertson, βA theoretical model for enzymatic catalysis using asymmetric hollow fiber membranes,β AIChE Journal, vol. 20, no. 1, pp. 50β59, 1974. [11] V. K. Jayaraman, βThe solution of hollow-fiber bioreactor design equations,β Biotechnology Progress, vol. 8, no. 5, pp. 462β464, 1992. [12] E. Nagy, Basic Equations of the Mass Transport Through a Membrane Layer, Elsevier, Boston, Mass, USA, 2012. [13] L. J. Kelsey, M. R. Pillarella, and A. L. Zydney, βTheoretical analysis of convective flow profiels in a hollow-fiber membrane bioreactor,β Chemical Engineering Science, vol. 45, no. 11, pp. 3211β3220, 1990. [14] V. Calabr`o, S. Curcio, and G. Iorio, βA theoretical analysis of transport phenomena in a hollow fiber membrane bioreactor with immobilized biocatalyst,β Journal of Membrane Science, vol. 206, no. 1-2, pp. 217β241, 2002. [15] E. Nagy, βBasic equations of mass transfer through biocatalytic membrane layer,β Asia-Pacific Journal of Chemical Engineering, vol. 4, no. 3, pp. 270β278, 2009. [16] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, 1965. [17] M. L. Shuler and F. K. Kargi, Bioprocess Engineering: Basic Concepts, Pearson, Upper Saddle River, NJ, USA, 2nd edition, 2014. [18] B. Godongwana, D. De Jager, M. S. Sheldon, and W. Edwards, βThe effect of Streptomyces coelicolor development on the hydrodynamics of a vertically orientated capillary membrane gradostat reactor,β Journal of Membrane Science, vol. 333, no. 1-2, pp. 79β87, 2009.
Coefficient of series solution, defined in text π: Substrate concentration (g dmβ3 ) Substrate feed concentration (g dmβ3 ) π0 : Dimensionless substrate concentration πΆ = π/π0 : Substrate diffusivity (m2 sβ1 ) π·π΄π΅ : Fraction retentate π = π’1 /π’0 : β: Step-size in the π-dimension (m) π: Grid point index in the π-dimension π: Grid point index in the π§-dimension Bessel function of order π of the first π½π (π): kind π: Step-size in the π§-dimension (m) Saturation (or Michaelis) constant πΎπ : (g dmβ3 ) β Dimensionless Michaelis constant πΎπ : πΏ: Membrane effective length (m) π(π, π, π): Kummer function of the first kind Peπ’ = π’0 π
πΏ /π·π΄π΅ : Axial Peclet number PeV = V0 π
πΏ /π·π΄π΅ : Radial Peclet number π: Radial spatial coordinate (m) Dimensionless radial spatial coordinate π
= π/π
1 : Membrane lumen radius (m) π
πΏ : π’: Axial velocity (m sβ1 ) Feed axial velocity (m sβ1 ) π’0 : Dimensionless axial velocity π = π’/π’0 : V: Radial velocity (m sβ1 ) Dimensionless radial velocity π = V/V0 : Maximum rate of reaction (g dmβ3 sβ1 ) ππ: π§: Axial spatial coordinate (m) π = π§/πΏ: Dimensionless axial spatial coordinate. Greek Letters πΌ:
Coefficients of finite-difference scheme, defined in text π½ = π0 βπ2 : Dimensionless transmembrane pressure π
: Dimensionless membrane hydraulic permeability π: Thiele modulus π = π
1 /πΏ: Aspect ratio Eigenvalues, π = 1, 2, . . .. ππ:
Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments The authors would like to thank the National Research Foundation (RSA) and the Fulbright Program (U.S. Department of State) for supporting this work. The MATLAB code for the Newton-Raphson algorithm was developed with the assistance of Professor Jeff Heys of Montana State University.
Advances in
Operations Research Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Advances in
Decision Sciences Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Applied Mathematics
Algebra
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Probability and Statistics Volume 2014
The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of
Differential Equations Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Submit your manuscripts at http://www.hindawi.com International Journal of
Advances in
Combinatorics Hindawi Publishing Corporation http://www.hindawi.com
Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Complex Analysis Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of Mathematics and Mathematical Sciences
Mathematical Problems in Engineering
Journal of
Mathematics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Discrete Mathematics
Journal of
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Discrete Dynamics in Nature and Society
Journal of
Function Spaces Hindawi Publishing Corporation http://www.hindawi.com
Abstract and Applied Analysis
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of
Journal of
Stochastic Analysis
Optimization
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014