After the discovery of fractal structures of financial markets, fractional partial differential equations (fPDEs) became very popular in studying dynamics of financial markets. Available research results involves two key modelling aspects; firstly, derivation of tractable asset pricing models, those that closely reflects the actual dynamics of financial markets. Secondly, the development of robust numerical solution methods. Often times, most the effective models are of a nonlinear nature, and as such, reliable analytical solution methods are seldomly available. On the other hand, the accurate value of American options strongly lies on the unknown free boundaries associated with these types of derivative contracts. The free boundaries emanates from the flexibility of the early exercise features with American options. To the best of our knowledge, the approach of pricing American options under the fractional calculus framework has not been extensively explored in literature, and an obvious wider research gap still exist on the design of robust solution methods for pricing American option problems formulated under the fractional calculus framework. Therefore, this paper serve to propose a robust numerical scheme for solving time-fractional Black–Scholes PDEs for pricing American put option problems. The proposed scheme is based on the front-fixing algorithm, under which the early exercise boundaries are transformed into fixed boundaries, allowing for a simultaneous computation of optimal exercise boundaries and their corresponding fair premiums. Results herein indicate that, the proposed numerical scheme is consistent, stable, convergent with order $${\mathcal {O}}(h^2,k)$$
O
(
h
2
,
k
)
, and also does guarantee positivity of solutions under all possible market conditions.