The second iteration of the optimal homotopy asymptotic technique (OHAM-2) has been protracted to fractional order partial differential equations in this work for the first time (FPDEs). Without any transformation, the suggested approach can be used to solve fractional-order nonlinear Zakharov–Kuznetsov equations. The Caputo notion of the fractional-order derivative, whose values fall within the closed interval [0, 1], has been taken into consideration. The method's appeal is that it provides an approximate solution after just one iteration. The suggested method's numerical findings have been contrasted with those of the variational iteration method, residual power series method, and perturbation iteration method. Through tables and graphs, the proposed method's effectiveness and dependability are demonstrated.