The removal of the bulk liquor surrounding the pulp fibers using less concentrated liquor is known as pulp washing. In the present study, a pulp washing model involving diffusion-dispersion through packed beds of finite length is presented. Separation of variables is applied to solve system of governing partial differential equations and the resulting equations are solved using Mathematica. Results from the present case are compared with those of previous investigators. The present case is giving better results than the previous investigators. 1. Introduction The objective of pulp and paper industry to produce its target production with high efficiency and less environmental load can only be met by initiating a meticulously planned research on mathematical methods. Pulp washing plays an important role in reduction of black liquor solids in the pulp being carried forward for further processing. The efficiency of washing depends on the degree of mixing, rate of desorption, diffusion-dispersion of dissolved solids, and chemicals from the fibrous matrix. Modeling of pulp washing is done mainly using three approaches, namely, (a) process modeling (b) physical modeling, and (c) statistical modeling. A complete review of the various process models used to describe pulp washing has been presented by [1]. Initially researchers like in [2, 3] proposed the models based on axial dispersion. Pellett [4] introduced a mathematical model combining the effects of particle diffusion and axial dispersion. A detailed model related to mass transfer in fibrous particle was given by [5]; it was also restricted for axial dispersion only. Comprehensive models involving physical features of the fibers such as fiber porosity and fiber radius were presented by [6, 7]. Extensive study of axial dispersion model has been carried out by [8–26]. The model has been solved using analytic and numerical techniques like Laplace transform technique [2–4, 10, 15, 23, 26], finite difference technique [25], orthogonal collocation method [5, 7, 12], orthogonal collocation on finite elements [6, 20, 21], Galerkin/Petrov Galerkin method [8, 19], Hermite collocation method by [11, 17, 24] and Spline collocation method [13]. The accuracy of the analytic solution undoubtedly exceeds the limit of applicability of the theory to real situations. Moreover, it is highly desirable to have a simple and consistent model of the transport phenomenon based on essential features of real situation. Keeping this modest goal in mind axial dispersion model is solved along with linear adsorption isotherm. The method
References
[1]
M. Pekkanen and H. V. Norden, “Review of pulp washing models,” Paperi Ja Puu, vol. 67, pp. 689–696, 1985.
[2]
H. Brenner, “The diffusion model of longitudinal mixing in beds of finite length. Numerical values,” Chemical Engineering Science, vol. 17, no. 4, pp. 229–243, 1962.
[3]
W. R. Sherman, “The movement of a soluble material during the washing of a bed of packed solids,” AIChE Journal, vol. 10, pp. 855–860, 1964.
[4]
G. L. Pellett, “Longitudinal dispersion, intra particle diffusion and liquid-phase mass transfer during flow through multi particle systems,” TAPPI Journal, vol. 49, pp. 75–82, 1966.
[5]
L. E. Grahs, Washing of cellulose fibres: analysis of displacement washing operation [Ph.D. dissertation], Department of Chemical Engineering, Chalmers University of Technology, Gothenburg, Sweden, 1974.
[6]
S. Arora, S. S. Dhaliwal, and V. K. Kukreja, “Simulation of washing of packed bed of porous particles by orthogonal collocation on finite elements,” Computers & Chemical Engineering, vol. 30, no. 6-7, pp. 1054–1060, 2006.
[7]
N. S. Raghavan and D. M. Ruthven, “Numerical simulation of a fixed-bed adsorption column by the method of orthogonal collocation,” AIChE Journal, vol. 29, no. 6, pp. 922–925, 1983.
[8]
M. Al-Jabari, A. R. P. van Heiningen, and T. G. M. van de Ven, “Modeling the flow and the deposition of fillers in packed beds of pulp fibres,” Journal of Pulp and Paper Science, vol. 20, no. 9, pp. J249–J253, 1994.
[9]
S. Arora and F. Pot？？ek, “Modelling of displacement washing of packed bed of fibers,” Brazilian Journal of Chemical Engineering, vol. 26, no. 2, pp. 385–393, 2009.
[10]
B. V. Babu and A. S. Chaurasia, “Pyrolysis of biomass: Improved models for simultaneous kinetics and transport of heat, mass and momentum,” Energy Conversion and Management, vol. 45, no. 9-10, pp. 1297–1327, 2004.
[11]
K. D. Edoh, R. D. Russell, and W. Sun, “Computation of invariant tori by orthogonal collocation,” Applied Numerical Mathematics, vol. 32, no. 3, pp. 273–289, 2000.
[12]
L. T. Fan, G. K. C. Chen, and L. E. Erickson, “Efficiency and utility of collocation methods in solving the performance equations of flow chemical reactors with axial dispersion,” Chemical Engineering Science, vol. 26, no. 3, pp. 378–387, 1971.
[13]
B. Gupta and V. K. Kukreja, “Modelling & simulation of packed bed of porous particles by orthogonal spline collocation method,” Applied Mathematics & Computation, vol. 219, pp. 2087–2099, 2012.
[14]
V. K. Kukreja, Modeling of washing of brown stock on rotary vacuum washer [Ph.D. dissertation], Department of Pulp and Paper Technology, University of Roorkee, Roorkee, India, 1996.
[15]
V. K. Kukreja and A. K. Ray, “Mathematical modeling of a rotary vacuum washer used for pulp washing: a case study of a lab scale washer,” Cellulose Chemistry & Technology, vol. 43, no. 1–3, pp. 25–36, 2009.
[16]
A. Kumar, D. K. Jaiswal, and N. Kumar, “Analytical solutions of one-dimensional advection-diffusion equation with variable coefficients in a finite domain,” Journal of Earth System Science, vol. 118, no. 5, pp. 539–549, 2009.
[17]
C. P. Le？o and A. E. Rodrigues, “Transient and steady-state models for simulated moving bed processes: numerical solutions,” Computers & Chemical Engineering, vol. 28, no. 9, pp. 1725–1741, 2004.
[18]
H.-T. Liao and C.-Y. Shiau, “Analytical solution to an axial dispersion model for the fixed-bed adsorber,” AIChE Journal, vol. 46, no. 6, pp. 1168–1176, 2000.
[19]
F. Liu and S. K. Bhatia, “Application fo Petrov-Galerkin methods to transient boundary value problems in chemical engineering: adsorption with steep gradients in bidisperse solids,” Chemical Engineering Science, vol. 56, no. 12, pp. 3727–3735, 2001.
[20]
Y. Liu and E. W. Jacobsen, “On the use of reduced order models in bifurcation analysis of distributed parameter systems,” Computers & Chemical Engineering, vol. 28, no. 1-2, pp. 161–169, 2004.
[21]
Z. Ma and G. Guiochon, “Application of orthogonal collocation on finite elements in the simulation of non-linear chromatography,” Computers & Chemical Engineering, vol. 15, no. 6, pp. 415–426, 1991.
[22]
F. Pot？？ek, “Washing of pulp fibre bed,” Collection of Czechoslovak Chemical Communications, vol. 62, no. 4, pp. 626–644, 1997.
[23]
A. Rasmuson and I. Neretnieks, “Exact solution of a model for diffusion in particles and longitudinal dispersion in packed beds,” AIChE Journal, vol. 26, pp. 686–690, 1980.
[24]
A. La Rocca and H. Power, “A double boundary collocation Hermitian approach for the solution of steady state convection-diffusion problems,” Computers and Mathematics with Applications, vol. 55, no. 9, pp. 1950–1960, 2008.
[25]
L. M. Sun and F. Meunier, “An improved finite difference method for fixed-bed multicomponent sorption,” AIChE Journal, vol. 37, no. 2, pp. 244–254, 1991.
[26]
Y. Zheng and T. Gu, “Analytical solution to a model for the startup period of fixed-bed reactors,” Chemical Engineering Science, vol. 51, no. 15, pp. 3773–3779, 1996.
[27]
F. Shiraishi, “Highly accurate solution of the axial dispersion model expressed in S-system canonical form by Taylor series method,” Chemical Engineering Journal, vol. 83, no. 3, pp. 175–183, 2001.
[28]
M. K. Szukiewicz, “New approximate model for diffusion and reaction in a porous catalyst,” AIChE Journal, vol. 46, no. 3, pp. 661–665, 2000.
