In this paper, we first introduce the 2-variables Konhauser matrix polynomials; then, we investigate some properties of these matrix polynomials such as generating matrix relations, integral representations, and finite sum formulae. Finally, we obtain the fractional integrals of the 2-variables Konhauser matrix polynomials.

In this paper, we first introduce the incomplete extended Gamma and Beta functions with matrix parameters; then, we establish some different properties for these new extensions. Furthermore, we give a specific application for the incomplete Bessel matrix function by using incomplete extended Gamma and Beta functions; at last, we construct the relation between the incomplete confluent hypergeometric matrix functions and incomplete Bessel matrix function.

The principal object of this paper is to introduce two variable Shivley’s matrix polynomials and derive their special properties. Generating matrix functions, matrix recurrence relations, summation formula and operational representations for these polynomials are deduced. Finally, Some special cases and consequences of our main results are also considered.

In recent years, much attention has been paid to the role of degenerate versions of special functions and polynomials in mathematical physics and engineering. In the present paper, we introduce a degenerate Euler zeta function, a degenerate digamma function, and a degenerate polygamma function. We present several properties, recurrence relations, infinite series, and integral representations for these functions. Furthermore, we establish identities involving hypergeometric functions in terms of degenerate digamma function.

We propose a hybrid numerical framework for solving time-fractional Navier–Stokes equations with nonlinear damping. The method combines the finite difference L1 scheme for time discretization of the Caputo derivative ($0<\alpha <1$) with mixed finite element methods (P1b–P1 and $P_2$–$P_1$) for spatial discretization of velocity and pressure. This approach addresses the key challenges of fractional models, including nonlocality and memory effects, while maintaining stability in the presence of the nonlinear damping term ${\gamma \left|\mathbf{u}\right|}^{r-2}\mathbf{u}$, for $r\ge 2$. We prove unconditional stability for both semi-discrete and fully discrete schemes and derive optimal error estimates for the velocity and pressure components. Numerical experiments validate the theoretical results. Convergence tests using exact solutions, along with benchmark problems such as backward-facing channel flow and lid-driven cavity flow, confirm the accuracy and reliability of the method. The computed velocity contours and streamlines show close agreement with analytical expectations. This scheme is particularly effective for capturing anomalous diffusion in Newtonian and turbulent flows, and it offers a strong foundation for future extensions to viscoelastic and biological fluid models.

We present in this paper a generalization of the fractional kinetic equation using the generalized incomplete Wright hypergeometric function. The pathway-type transform technique is then used to investigate the solutions to a fractional kinetic equation with specific fractional transforms. Furthermore, exceptional cases of our outcomes are discussed and graphically illustrated using MATLAB software. This work provides a thorough overview for further investigation into these topics in order to gain a better understanding of their implications and applications.