Generalized semi-analytical solution of advection-diffusion-reaction in finite and semi-infinite cylindrical ducts

A semi-analytical solution for the transient advection-diffusion-reaction problem within finite and semi-infinite ducts is derived. The solution allows for general radial- and time-dependent inlet/outlet conditions, complex boundary conditions on the duct wall including adsorption and decay, and arbitrary velocity profiles of the transporting fluid. The only numerical step of the solution is the inverse Laplace transform in the time variable. Therefore, the approach also produces fully analytical steady-state solutions. The solution is verified against computational fluid dynamics (CFD) simulations under various boundary conditions and velocity profiles (Newtonian and power-law), and in all cases good agreement is obtained. Although theoretically applicable to all regimes, the solution is computationally difficult at very high Peclet numbers and very early times due to numerical instabilities as a result of finite precision arithmetic of computers. A convergence analysis is conducted to delineate the boundaries of this limit for two important cases. The solution was derived using a new approach for solving two-dimensional partial differential equations (PDEs) with non-constant coefficients which parallels the Frobenius and power series methods for solving ordinary differential equations (ODEs). The approach reduces the original PDE to a single infinite-order ODE with constant coefficients. The approach is suspected to provide solutions to a large class of PDEs of this type. The solution may find applications in a number of engineering and/or biomedical fields, it can be used to verify numerical simulators, and serve as a simple and easy-to-implement alternative where access to numerical simulators is not available.