This research explores the complex dynamics of a Novel Four-Dimensional Fractional Supply Chain System (NFDFSCS) that integrates a quadratic interaction term involving the actual demand of customers and the inventory level of distributors. The introduction of the quadratic term results in significantly larger maximal Lyapunov exponents (MLE) compared to the original model, indicating increased system complexity. The existence, uniqueness, and Ulam–Hyers stability of the proposed system are verified. Additionally, we establish the global Mittag-Leffler attractive set (MLAS) and Mittag-Leffler positive invariant set (MLPIS) for the system. Numerical simulations and MATLAB phase portraits demonstrate the chaotic nature of the proposed system. Furthermore, a dynamical analysis achieves verification via the Lyapunov exponents, a bifurcation diagram, a 0–1 test, and a complexity analysis. A new numerical approximation method is proposed to solve non-linear fractional differential equations, utilizing fractional differentiation with a non-singular and non-local kernel. These numerical simulations illustrate the primary findings, showing that both external and internal factors can accelerate the process. Furthermore, a robust control scheme is designed to stabilize the system in finite time, effectively suppressing chaotic behaviors. The theoretical findings are supported by the numerical results, highlighting the effectiveness of the control strategy and its potential application in real-world supply chain management (SCM).