This study uses the optimal auxiliary function method to approximate solutions for fractional-order non-linear partial differential equations, utilizing Riemann–Liouville’s fractional integral and the Caputo derivative. This approach eliminates the need for assumptions about parameter magnitudes, offering a significant advantage. We validate our approach using the time-fractional Cahn–Hilliard, fractional Burgers–Poisson, and Benjamin–Bona–Mahony–Burger equations. Comparative testing shows that our method outperforms new iterative, homotopy perturbation, homotopy analysis, and residual power series methods. These examples highlight our method’s effectiveness in obtaining precise solutions for non-linear fractional differential equations, showcasing its superiority in accuracy and consistency. We underscore its potential for revealing elusive exact solutions by demonstrating success across various examples. Our methodology advances fractional differential equation research and equips practitioners with a tool for solving non-linear equations. A key feature is its ability to avoid parameter assumptions, enhancing its applicability to a broader range of problems and expanding the scope of problems addressable using fractional calculus techniques.

1. The ultrasonic-guided waves testing has become a common technique in the non-destructive testing (NDT) because of its capability of inspecting long distances to detect structural changes or defe...