Abstract The homotopy analysis method (HAM), a general analytic technique for non-linear problems, is applied to analyze the solute dispersion process in non-Newtonian Carreau-Yasuda and Carreau fluids flow in a straight tube with the effect of wall absorption/reaction. Unlike the other analytical methods such as perturbation method, the HAM provides a simple way to get the convergent series solution. Any assumptions of small or high physical quantities are not required for the HAM. It provides us a great freedom to choose the so-called convergence control parameter which is used to guarantee the convergence of series solution. The convergent series solution is obtained by choosing the optimal value of convergence control parameter for which the series converges fastest. The optimal value of convergence control parameter is obtained by minimizing the square residual which also provides the convergence region for the series solution. The previous analytical studies on solute dispersion fail to justify the convergence of the series solution, whereas in this investigation, our results are convergent and valid for all physical parameters. In addition, present results are validated by the numerical and some existing results. This study explains elaborately about the advantage of the homotopy analysis method over perturbation and eigenfunction expansion methods for nonlinear problems. Owing to the great potential and flexibility of the homotopy analysis method, it is a more suitable analytical approach to solve the different types of non-linear problems in science and engineering. Besides, this study helps to gain some knowledge on the transportation process of drugs in the blood flow.