We propose a new type of asymptotic expansion for the transition probability density function (or heat kernel) of certain parabolic partial differential equations (PDEs) that appear in option pricing. As other, related methods developed by Costanzino, Hagan, Gatheral, Lesniewski, Pascucci, and their collaborators, among others, our method is based on the computation of the truncated asymptotic expansion of the heat kernel with respect to a “small” parameter. What sets our method apart is that our small parameter is possibly different from the time to expiry and that the resulting commutator calculations go beyond the nilpotent Lie algebra case. In favorable situations, the terms of this asymptotic expansion can quickly be computed explicitly leading to a “closed-form” approximation of the solution, and hence of the option price. Our approximations tend to have much fewer terms than the ones obtained from short time asymptotics, and are thus easier to generalize. Another advantage is that the first term of our expansion corresponds to the classical Black-Scholes model. Our method also provides equally fast approximations of the derivatives of the solution, which is usually a challenge. A full theoretical justification of our method seems very difficult at this time, but we do provide some justification based on the results of (Siyan, Mazzucato, and Nistor, NWEJ 2018). We therefore mostly content ourselves to demonstrate numerically the efficiency of our method by applying it to the solution of the mean-reverting SABR stochastic volatility model PDE, commonly referred to as the [Formula: see text]SABR PDE, by taking the volatility of the volatility parameter [Formula: see text] (vol-of-vol) as a small parameter. For this PDE, we provide extensive numerical tests to gauge the performance of our method. In particular, we compare our approximation to the one obtained using Hagan’s formula and to the one obtained using a new, adaptive finite difference method. We provide an explicit asymptotic expansion for the implied volatility (generalizing Hagan’s formula), which is what is typically needed in concrete applications. We also calibrate our model to observed market option price data. The resulting values for the parameters [Formula: see text], [Formula: see text], and [Formula: see text] are realistic, which provides more evidence for the conjecture that the volatility is mean-reverting.