Two different analytical solutions of the extended Graetz problem are analyzed in order to quantify the dynamics of particulate matter (PM) in a cylindrical element. Advection, diffusion and deposition of particles are calculated in a symmetrical circular tube with fully developed laminar flow applying the theory of confluent hypergeometric functions. The following in silico experiments are proposed: (i) the simulation of a thin tube deposition, acting as a filtering element for PM and (ii) the modeling of isokinetic sampling in the inner region of the tube. Two mathematical formulations are compared to obtain the general analytical solution of the advection-diffusion equation in cylindrical coordinates. After applying the variable separation method, a second order ordinary differential equation (ODE) is obtained, this ODE is solved with two distinct forms: (a) applying a transformation of variables to Whittaker’s function and (b) using a new variable transformation proposition. No studies were found that compare these analytical solutions, as well as the use of this new variable transformation, although both cases are essentially different ways of applying the Frobenius method. In this work, 75 eigenvalues and their constants are presented for the first time using five decimal places. The new series solution has high accuracy compared to the previous one, it can provide the concentration profile, the deposition rate and the determination of the PM profile in the isokinetic sampling region. Our solution proved to be more stable close to the tube wall, which may improve the techniques for measuring the flow of particles in tubes. Furthermore, using more eigenvalues improved the estimation of PM2.5 and PM10 deposition. The results presented here show that the proposed new analytic solution has the potential to contribute to numerical and experimental environmental analyses.
Read full abstract