The aim of the present study is to obtain approximate solutions of the fractional counterpart of a boundary value problem that appears in electrohydrodynamic flows by using generalized differential transform method (in short, GDTM). Convergence of the solution problem is proved. The efficiency and applicability of the present method are demonstrated by solving the equation for selected values of the parameters that appear in the model. In addition, residual error computation is adopted to verify the validity of the results. © 2019, Società Italiana di Fisica and Springer-Verlag GmbH Germany, part of Springer Nature.