Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Analysis of newly developed fractal-fractional derivative with power law kernel for MHD couple stress fluid in channel embedded in a porous medium

## Abstract

Fractal-fractional derivative is a new class of fractional derivative with power Law kernel which has many applications in real world problems. This operator is used for the first time in such kind of fluid flow. The big advantage of this operator is that one can formulate models describing much better the systems with memory effects. Furthermore, in real world there are many problems where it is necessary to know that how much information the system carries. To explain the memory in a system fractal-fractional derivatives with power law kernel is analyzed in the present work. Keeping these motivation in mind in the present paper new concept of fractal-fractional derivative for the modeling of couple stress fluid (CSF) with the combined effect of heat and mass transfer have been used. The magnetohydrodynamics (MHD) flow of CSF is taken in channel with porous media in the presence of external pressure. The constant motion of the left plate generates the CSF motion while the right plate is kept stationary. The non-dimensional fractal-fractional model of couple stress fluid in Riemann–Liouville sense with power law is solved numerically by using the implicit finite difference method. The obtained solutions for the present problem have been shown through graphs. The effects of various parameters are shown through graphs on velocity, temperature and concentration fields. The velocity, temperature and concentration profiles of the MHD CSF in channel with porous media decreases for the greater values of both fractional parameter $$\alpha$$ and fractal parameter $$\beta$$ respectively. From the graphical results it can be noticed that the fractal-fractional solutions are more general as compared to classical and fractional solutions of CSF motion in channel. Furthermore, the fractal-fractional model of CSF explains good memory effect on the dynamics of couple stress fluid in channel as compared to fractional model of CSF. Finally, the skin friction, Nusselt number and Sherwood number are evaluated and presented in tabular form.

## Introduction

Fractional calculus is used to explain many real world problems with better memory effect. Fractional calculus FC is the generalization of integer order calculus which was not sufficient to explain some memory effect in some engineering and real world problems. Due to its enormous applications in many fields of sciences the FC is getting attention of the researchers from the last three decades. In other words FC is used to describe the historical states of different phenomenon which we call as memory. Many researchers developed different definitions of fractional derivatives for different physical situations. Some applications of FC are given in the Tenreiro et al.1. Dalir and Bashour2 investigated some applications of FC in applied sciences. Tavazoei et al.3 studied some uses of FC in modern sciences and in the processes of suppression of chaotic oscillations. Similarly, Sabatier et al.4 discussed some advanced applications of FC in modern sciences. Koca and Atangana5 investigated Cattaneo-Hristov model with the effect of heat diffusion using Caputo-Fabrizio and Atangana-Baleanu fractional derivatives. In another paper Atangana and Baleanu6 discussed new fractional derivatives to highlight some applications and theory of heat and mass transfer models. Arif et al.7 investigated CSF in channel. In this study they generalized the classical model of CSF into time fractional models of AB and CF and compare their results. Podlubny8 discussed fractional differential equations and their uses in different physical phenomenon. Abdeljawad and Baleanu9 studied fractional derivatives and their applications with generalized Mittag–Leffler function.

Recently, Atangana10 developed the new idea of fractal-fractional derivative in FC. This new idea is very suitable in many situations in dealing some complex problems. In the operator there are two orders, the first is known as fractional order and the second operator is called the fractal dimension. This new idea of fractal fractional derivative is better than the classical one and fractional derivatives as well. It is due to the fact that by dealing with fractal-fractional derivatives we can study the fractional operator as well as fractal dimension at the same time. Motivated from the advance and unique properties many researchers taking interest in fractal-fractional operator. Ali et al.11 investigated SIR model to study the dynamics of COVID-19 using fractal-fractional operator. Recently, Akgul and Siddique12 investigated the fluid flow in channel with MHD effect using the idea of fractal-fractional derivatives on the CSF. Esmonde13 studied fractal-fractional derivative modeling for the applications of material phase change. In another paper Akgul14 discussed some advance and new applications of fractal fractional differential equations using power law kernel.

Mainly, fluids are classified in two categories, Newtonian and non-Newtonian. Most of the real world problems are analyzed by non-Newtonian fluids. Many researchers have shown interest in non- Newtonian fluids due to enormous applications in daily life. Some applications of non-Newtonian fluids flow like laminar flow, micro channel flows and pipe flows are given in15. Some industrial and engineering applications of non-Newtonian fluids are given in the book of Chhabra and Richardson16. Pawar and Sunnapwar17 examined some experimental study of non-Newtonian fluids in helical coils with laminar and turbulent flow of the non-Newtonian fluids. The advanced applications of non-Newtonian fluids like, oil pipeline friction reduction, heating and cooling systems, scale-up and flow tracers are investigated by Hoyt18. Sohail et al.19 investigated the features of non-Newtonian fluid and thermophysical characteristics of yield exhibiting non-Newtonian fluid flow under gyrotactic microorganisms. Bao et al.20 discussed the applications and practical uses of non-Newtonian fluids and calculated numerical modeling for this fluid flow in fractures and porous media. Some other applications of non-Newtonian fluid are given in21,22,23.

Couple stress fluid CSF is in the class of non-Newtonian fluids and have many engineering applications in real world. The idea of couple stresses in fluid was first given by Stokes24. After stokes many researchers used the model of CSF in different situations. Arif et al.25 discussed CSF in channel using AB fractional derivatives with gold nanoparticles suspended in the blood. Researchers are taking interest to analyzed CSF fluid in different domain for different scientific reasons. Like Adesanya et al.26 discussed irreversibility process using CSF in inclined channel with isothermal boundaries. Krishna and Chamkha27 studied engineering applications using the concept of CSF through a porous media with slip effect. Khan et al.28 investigated some biomechanical problems of CSF with expanding or contracting porous channel. Akhtar et al.29 studied CSF in channel using fractional derivatives.

MHD play a vital role in different fluid flow problems and have many applications in biological sciences and engineering sciences. The applications of MHD is not limited to fluid flow problems but it can be used in many engineering problems like Kumar et al.30 discussed the impact of magnetic dipole on the fluid flow over a stretching cylinder. Gowda et al.31 investigated MHD effect in their study using magnetized flow with the theory of heat diffusion and Stefan blowing condition. In another paper Gowda et al.32 examined the contribution of MHD flow over a stretching sheet exploring the magnetic dipole on the fluid flow. Yusuf et al.33 discussed magneto-Bioconvection flow over an inclined plate with entropy generation. Kumar et al.34 discussed the impact of MHD on thermophoretic particle in the fluid flow over a stretching sheet.

Heat and mass transfer have useful applications in industries and different engineering problems. Researchers investigate the effect of heat and mass transfer in different circumstances like Naveen et al.35 examined the influence of heat and mass transfer on the fluid flow along a stretching cylinder. Similarly Punith et al.36 discussed the impact of binary chemical reaction with the combined effect of heat and mass transfer on the fluid flow problem. Similarly other researchers examined the applications of heat and mass transfer in different field of sciences which are given in37,38,39.

From the above mentioned literature no work is reported to investigate the MHD flow of CSF in channel with heat and mass transfer. The fluid is taken in channel with porous media. The left plate is moving with constant velocity which generates the fluid motion while the right plate is fixed. The main novelty of the present article is to apply the fractal-fractional derivatives on the classical model. Fractal-fractional derivative is a new class of fractional derivative with power Law kernel which has many applications in real world problems. This operator is used for the first time in such kind of fluid flow. The big advantage of this operator is that we can formulate models describing much better the systems with memory effects. Furthermore, in real world there are many problems where we need to know that how much information the system carries that is why need memory in a system which is explained by fractal-fractional derivatives with power law kernel. This distribution in memory effect in fractional derivatives follows a power-law distribution. Keeping these motivations in mind the classical model of CSF is transform to fractal-fractional model and the numerical solutions have been obtained via implicit finite difference method.

## Formulation of the problem

In the present study we have considered the incompressible unsteady MHD flow of CSF in channel. The couple stress fluid is taken in channel in the presence of external pressure with body couples. The fluid is taken between two plates and porous medium is considered. Initially, both plates and fluid were at rest after some time the left plate start moving with constant velocity due to which the fluid flow in channel and the right plate is fixed. The governing equations for the present flow regime are given below:

The continuity equation is given by25,29,40:

$$\nabla \cdot {\mathbf{\mathop{V}\limits^{\rightharpoonup} }}{ = 0,}$$
(1)

The momentum equation of the given problem is given by:

$$\rho \frac{{D{\mathbf{\mathop{V}\limits^{\rightharpoonup} }}}}{Dt} = - \nabla p - \mu \nabla \times \left( {\nabla \times {\mathbf{\mathop{V}\limits^{\rightharpoonup} }}} \right) - \eta \nabla \times \left( {\nabla \times \left( {\nabla \times \left( {\nabla \times {\mathbf{\mathop{V}\limits^{\rightharpoonup} }}} \right)} \right)} \right) + \rho \vec{b} + r,$$
(2)

The energy equation can be written as:

$$\rho c_{p} \frac{{\partial {\mathbf{\mathop{T}\limits^{\rightharpoonup} }}}}{\partial t} = k\nabla \times \nabla \times {\mathbf{\mathop{T}\limits^{\rightharpoonup} }},$$
(3)

The concentration equation can be written as:

$$\frac{{\partial {\mathbf{\mathop{C}\limits^{\rightharpoonup} }}}}{\partial t} = D\nabla \times \nabla \times {\mathbf{\mathop{C}\limits^{\rightharpoonup} }},$$
(4)

where in the momentum equation the term $$\rho \overrightarrow {b}$$ shows the body forces which can be expressed as:

$$\rho \mathop{b}\limits^{\rightharpoonup} = {\mathbf{\mathop{J}\limits^{\rightharpoonup} }} \times {\mathbf{\mathop{B}\limits^{\rightharpoonup} }} + \rho \mathop{g}\limits^{\rightharpoonup} .$$
(5)

Here $${\mathbf{\mathop{V}\limits^{\rightharpoonup} }},{\mathbf{\mathop{T}\limits^{\rightharpoonup} }}$$ and $${\mathbf{\mathop{C}\limits^{\rightharpoonup} }}$$ are the velocity, temperature and concentration vectors. $$p,\rho ,\mu ,\rho \vec{b},k,D,\, \, \vec{g}\,\,{\text{and}}\,\,{\text{r}}$$ is the pressure, density, dynamic viscosity, body forces, thermal conductivity, thermal diffusivity, gravitational acceleration and Darcy resistance respectively. $${\mathbf{\mathop{J}\limits^{\rightharpoonup} }}$$ and $${\mathbf{\mathop{B}\limits^{\rightharpoonup} }}$$ is the current density and total magnetic field. As we have considered unidirectional flow, therefore, velocity, temperature and concentration fields of the given flow are as:

$$\left. \begin{gathered} {\mathbf{\mathop{V}\limits^{\rightharpoonup} }} = \left( {u_{1} (y,t),0,0} \right), \hfill \\ {\mathbf{\mathop{T}\limits^{\rightharpoonup} }} = \left( {T(y,t),0,0} \right), \hfill \\ {\mathbf{\mathop{C}\limits^{\rightharpoonup} }} = \left( {C(y,t),0,0} \right). \hfill \\ \end{gathered} \right\}$$
(6)

In the given study we have consider Darcy resistance in CSF, therefor, Darcy’s law can be written in the following form:

$$- \frac{\mu \phi }{{k_{0} }}.\overrightarrow {V} = r,$$
(7)

where $$\phi$$ represent porous media and $$k$$ represent permeability of the porous medium. Maxwell equation can be defined as:

From Maxwell equation:

$$\overrightarrow {\nabla } \cdot {\varvec{B}} = 0,$$
(8)
$$\overrightarrow {\nabla } \times \overrightarrow {E} = - \frac{{\partial {\varvec{B}}}}{\partial E} = 0,{\text{Where }}E = 0,$$
(9)

here $$E$$ is the total electric field. By Ohm’s Law (generalized form):

$${\varvec{J}} = \sigma \left( {\overrightarrow {E} + \overrightarrow {V} \times {\varvec{B}}} \right) = \sigma \left( {\overrightarrow {V} \times {\varvec{B}}} \right).$$
(10)

Cross product with the magnetic field is:

$${\varvec{J}} \times {\varvec{B}} = \sigma \left( {\overrightarrow {V} \times {\varvec{B}}} \right) \times \user2{B, B} = B_{0} + b,$$
(11)

here $$B_{0}$$ is applied magnetic field and $$b$$ is induced magnetic field by polarization. Now applying the vector scalar triple product, Eq. (11) becomes:

$$\user2{J } \times {\varvec{B}} = - \sigma \left\{ {{\varvec{B}} \times \left( {\overrightarrow {V} \times {\varvec{B}}} \right)} \right\} = - \sigma \left\{ {\left( {{\varvec{B}} \cdot {\varvec{B}}} \right)\overrightarrow {V} - \left( {{\varvec{B}} \cdot \overrightarrow {V} } \right){\varvec{B}}} \right\}$$
(12)
$${\varvec{J}} \times {\varvec{B}} = - \sigma \left\{ {\left( {{\varvec{B}} \cdot {\varvec{B}}} \right)\overrightarrow {V} - 0} \right\} = - \sigma {\varvec{B}}_{0}^{2} u.$$
(13)

In this article we have considered the incompressible unsteady MHD flow of CSF in channel. The MHD CSF laminar fluid is considered to flow through an open channel of two parallel plates separated by a distance $$d$$. The medium is considered as a porous medium with porosity K in the presence of constant external pressure gradient G and the induced magnetic field B0 which is taken normal to the fluid flow. The motion of the fluid is considered in x-direction. Initially, for $$t \le 0$$ the fluid and both the plates are stationary with surrounding temperature $$T_{\infty }$$ and constant concentration $$C_{\infty }$$. At $$t = 0^{ + }$$, the temperature and concentration of the left plate raised to $$T$$ and $$C$$ respectively. As a result the left plate moving with constant velocity and the right plate is stationary. The geometrical sketch of the considered model is given in Fig. 1.

Using the assumptions which are considered in the problem the governing equations for the flow, energy and concentration equations are given by2,25:

$$\rho \frac{\partial u}{{\partial t}} = G + \mu \frac{{\partial^{2} u}}{{\partial y^{2} }} - \eta \frac{{\partial^{4} u}}{{\partial y^{4} }} - \sigma {\varvec{B}}_{0}^{2} u - \frac{\mu \phi }{{k_{0} }}u + g\rho \beta (T - T_{\infty } ) + g\rho \beta (C - C_{\infty } ),$$
(14)
$$\left( {\rho c_{p} } \right)\frac{{\partial T\left( {y,t} \right)}}{\partial t} = k\frac{{\partial^{2} T\left( {y,t} \right)}}{{\partial y^{2} }}$$
(15)
$$\frac{{\partial C\left( {y,t} \right)}}{\partial t} = D\frac{{\partial^{2} C\left( {y,t} \right)}}{{\partial y^{2} }}$$
(16)

with the physical initial and boundary conditions:

$$\left. \begin{gathered} u(y,0) = 0, \, T(y,0) = T_{\infty } , \, C(y,0) = C_{\infty } ,{\text{ for }}0 \le y \le d, \hfill \\ u(0,t) = H(t)U_{0} , \, T(0,t) = T_{w} , \, C(0,t) = C_{w} ,{\text{ for }}t > 0, \hfill \\ u(d,t) = 0, \, T(d,t) = T_{\infty } , \, C(d,t) = C_{\infty } ,{\text{ for }}t > 0, \hfill \\ \frac{{\partial^{2} u(0,t)}}{{\partial y^{2} }} = \frac{{\partial^{2} u(d,t)}}{{\partial y^{2} }} = 0,{\text{ for }}t > 0. \hfill \\ \end{gathered} \right\}$$
(17)

From the above initial and boundary conditions it can be observed that the fluid and plates were at rest initially. After some time the left plate $$u(0,t) = H(t)U_{0}$$ moving with constant velocity $$U_{0}$$, where $$H(t)$$ shows the Heaviside step function and the right plate is fixed. At the left plate the wall temperature and concentration and at the right plate there is ambient temeprature and concentration. Furthermore, $$\frac{{\partial^{2} u(0,t)}}{{\partial y^{2} }} = \frac{{\partial^{2} u(d,t)}}{{\partial y^{2} }} = 0$$, shows the couple shear stresses at the left and right plates.

For dimensional analysis, the following non-dimensional variables are introduced:

$$\left. \begin{gathered} \xi = \frac{y}{d},{\text{ w}} = \frac{u}{{U_{0} }}, \, \tau = \frac{{U_{0} t}}{d}, \, \eta_{A} = \frac{\eta }{{\mu d^{2} }}, \hfill \\ P = \frac{{d^{2} }}{{\mu U_{0} }}G, \, \Theta = \frac{{T - T_{\infty } }}{{T_{w} - T_{\infty } }}, \, \Phi = \frac{{C - C_{\infty } }}{{C_{w} - C_{\infty } }}, \hfill \\ \end{gathered} \right\}$$
(18)

After dimensionalization process we get the following dimensionless system of equations along with initial and boundary conditions.

$${\text{Re}} \frac{{\partial w\left( {\xi ,\tau } \right)}}{\partial \tau } = P + \frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} - Hw\left( {\xi ,\tau } \right) + Gr\Theta \left( {\xi ,\tau } \right) + Gm\Phi \left( {\xi ,\tau } \right),$$
(19)
$$A\frac{{\partial \Theta \left( {\xi ,\tau } \right)}}{\partial \tau } = \frac{{\partial^{2} \Theta \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}$$
(20)
$$B\frac{{\partial \Phi \left( {\xi ,\tau } \right)}}{\partial \tau } = \frac{{\partial^{2} \Phi \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}$$
(21)
$$\left. \begin{gathered} w\left( {\xi ,0} \right) = 0, \, \Theta \left( {\xi ,0} \right) = 0, \, \Phi \left( {\xi ,0} \right) = 0,{\text{ for }}0 \le \xi \le d \hfill \\ w\left( {0,\tau } \right) = 1, \, \Theta \left( {0,\tau } \right) = 1, \, \Phi \left( {0,\tau } \right) = 1,{\text{ for }}\tau > 0, \hfill \\ w\left( {d,\tau } \right) = 0, \, \Theta \left( {d,\tau } \right) = 0, \, \Phi \left( {d,\tau } \right) = 0,{\text{ for }}\tau > 0, \hfill \\ \frac{{\partial^{2} w\left( {0,\tau } \right)}}{{\partial \xi^{2} }} = \frac{{\partial^{2} w\left( {d,\tau } \right)}}{{\partial \xi^{2} }} = 0,{\text{ for }}\tau > 0, \hfill \\ \end{gathered} \right\}$$
(22)
$$Gr = \frac{{g\rho \beta_{T} d^{2} \left( {T_{w} - T_{\infty } } \right)}}{{U_{0} \mu }}, \, Gm = \frac{{g\rho \beta_{C} d^{2} \left( {C_{w} - C_{\infty } } \right)}}{{U_{0} \mu }},\,\,\Pr = \frac{{\mu c_{p} }}{k}, \, {\text{Re}} = \frac{{U_{0} d}}{\upsilon }, \, Sc = \frac{\upsilon }{D},$$
$$M = \frac{{\sigma {\varvec{B}}_{0}^{2} d^{2} }}{\mu },\,\frac{1}{K} = \frac{{d^{2} \phi }}{{k_{0} }},\,\,\,\,H = M + \frac{1}{K},\,\,A = \Pr .{\text{Re}} ,\,\,B = Sc.{\text{Re}} .$$

where $$Gr$$ represents Grashof number, $$Gm$$ mass Grashof number, $$\Pr$$ represents Prandtl number, $${\text{Re}}$$ represents Reynolds number, $$Sc$$ represents Schmidth number, $$M$$ represnets magnetic parameter,$$K$$ porosity parameter and $$H$$ Hartmann number.

## Definition of Fractal-Fractional Derivative with Power Law Kernal

Let assume that $$f(t)$$ is continuous in the interval (a, b) and let the function is fractal differentiable on (a, b) having order $$\beta$$, then the fractal-fractional derivative of $$f$$ having order $$\alpha$$ in Riemann–Liouville RL sense with power law kernel is given by10:

$${}_{a}^{FFP} D_{t}^{\alpha ,\beta } f(t) = \frac{1}{\Gamma [n - \alpha ]}\frac{d}{{dt^{\beta } }}\int\limits_{a}^{t} {f\left( y \right)\left( {t - y} \right)^{n - \alpha - 1} dy,\,\,} n - 1 < \alpha \le n,\,\,\,\,0 < n - 1 < \beta \le n,$$
(23)

where

$$\frac{df\left( y \right)}{{dy^{\beta } }} = \mathop {\lim }\limits_{t \to y} \frac{f\left( t \right) - f\left( y \right)}{{t^{\beta } - y^{\beta } }}.$$
(24)

## Solutions of CSF with fractal-fractional derivative

In order to transform the classical CSF model into fractal-fractional derivative Eqs. (1921) can be written in the following form:

$${}^{FFP}D_{\tau }^{\alpha ,\beta } w\left( {\xi ,\tau } \right) = \frac{1}{{\text{Re}}}\left( {P + \frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} - Hw\left( {\xi ,\tau } \right) + Gr\Theta \left( {\xi ,\tau } \right) + Gm\Phi \left( {\xi ,\tau } \right)} \right),$$
(25)
$${}^{FFP}D_{\tau }^{\alpha ,\beta } \Theta \left( {\xi ,\tau } \right) = \frac{1}{A}\left( {\frac{{\partial^{2} \Theta \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right),$$
(26)
$${}^{FFP}D_{\tau }^{\alpha ,\beta } \Phi \left( {\xi ,\tau } \right) = \frac{1}{B}\left( {\frac{{\partial^{2} \Phi \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right).$$
(27)

where $${}^{FFP}D_{t}^{\alpha ,\beta } (.,.)$$ shows the fractal-fractional derivative, $$0 < \alpha \le 1$$ is the fractional order and $$0 < \beta \le1$$ is the fractal dimension.

### Solution of Energy Equation

From Eq. (26), the following result is obtained:

$${}^{FFP}D_{\tau }^{\alpha ,\beta } \Theta \left( {\xi ,\tau } \right) = \frac{1}{A}\left( {\frac{{\partial^{2} \Theta \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right),$$
(28)

Equation (28) can be written as:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\frac{d}{d\tau }\int\limits_{0}^{\tau } {\Theta \left( {\xi ,t} \right)\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{A}} .\left( {\frac{{\partial^{2} \Theta \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right),$$
(29)

the above result can be written in the following form:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{\tau } {\frac{{\partial \Theta \left( {\xi ,t} \right)}}{\partial t}.\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{A}} .\left( {\frac{{\partial^{2} \Theta \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right) - \frac{{\tau^{ - \alpha } \Theta \left( {\xi ,0} \right)}}{{\Gamma \left( {1 - \alpha } \right)}},$$
(30)

using the given initial condition from Eq. (22), Eq. (30) reduces to the following form:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{\tau } {\frac{{\partial \Theta \left( {\xi ,t} \right)}}{\partial t}.\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{A}} .\left( {\frac{{\partial^{2} \Theta \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right),$$
(31)

by discretizing the above equation at $$\left( {\xi_{i} ,\,\tau = \tau_{n + 1} } \right),$$ the following form is obtained:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\frac{d}{d\tau }\int\limits_{0}^{{\tau_{n + 1} }} {\frac{{\partial \Theta \left( {\xi_{i} ,t} \right)}}{\partial t}\left( {\tau_{n + 1} - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau_{n + 1}^{\beta - 1} }}{A}} .\left( {\frac{{\Theta_{i + 1}^{n + 1} - 2\Theta_{i}^{n + 1} + \Theta_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }}} \right),$$
(32)

From the above equation the following results are obtained:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\sum\limits_{j = 0}^{n} {\frac{{\Theta_{i}^{j + 1} - 2\Theta_{i}^{j} }}{\Delta \tau }} .\int\limits_{{\tau_{j} }}^{{\tau_{j + 1} }} {\left( {\tau_{n + 1} - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau_{n + 1}^{\beta - 1} }}{A}} .\left( {\frac{{\Theta_{i + 1}^{n + 1} - 2\Theta_{i}^{n + 1} + \Theta_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }}} \right),$$
(33)
$$\frac{{\left( {\Delta \tau } \right)^{ - \alpha } }}{{\Gamma \left( {2 - \alpha } \right)}}.\sum\limits_{j = 0}^{n} {\frac{{\Theta_{i}^{j + 1} - 2\Theta_{i}^{j} }}{\Delta \tau }} .\left( {\left( {n - j + 1} \right)^{1 - \alpha } - \left( {n - j} \right)^{1 - \alpha } } \right) = \frac{\beta }{A}\left( {\left( {n + 1} \right)\Delta \tau } \right)^{\beta - 1} \left( {\frac{{\Theta_{i + 1}^{n + 1} - 2\Theta_{i}^{n + 1} + \Theta_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }}} \right).$$
(34)

### Solution of concentration equation

From equation (27), the following result is obtained:

$${}^{FFP}D_{\tau }^{\alpha ,\beta } \Phi \left( {\xi ,\tau } \right) = \frac{1}{B}\left( {\frac{{\partial^{2} \Phi \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right),$$
(35)

Equation (35) can be written as:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\frac{d}{d\tau }\int\limits_{0}^{\tau } {\Phi \left( {\xi ,t} \right)\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{B}} .\left( {\frac{{\partial^{2} \Phi \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right),$$
(36)

here

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{\tau } {\frac{{\partial \Phi \left( {\xi ,t} \right)}}{\partial t}.\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{B}} .\left( {\frac{{\partial^{2} \Phi \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right) - \frac{{\tau^{ - \alpha } \Phi \left( {\xi ,0} \right)}}{{\Gamma \left( {1 - \alpha } \right)}},$$
(37)

using the given initial condition from Eq. (22), Eq. (37) reduces to the following form:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{\tau } {\frac{{\partial \Phi \left( {\xi ,t} \right)}}{\partial t}.\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{B}} .\left( {\frac{{\partial^{2} \Phi \left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }}} \right),$$
(38)

by discretize the above equation at $$\left( {\xi_{i} ,\,\tau = \tau_{n + 1} } \right),$$ one can get the following result:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\frac{d}{d\tau }\int\limits_{0}^{{\tau_{n + 1} }} {\frac{{\partial \Phi \left( {\xi_{i} ,t} \right)}}{\partial t}\left( {\tau_{n + 1} - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau_{n + 1}^{\beta - 1} }}{B}} .\left( {\frac{{\Phi_{i + 1}^{n + 1} - 2\Phi_{i}^{n + 1} + \Phi_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }}} \right),$$
(39)

the above result can be written as:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\sum\limits_{j = 0}^{n} {\frac{{\Phi_{i}^{j + 1} - 2\Phi_{i}^{j} }}{\Delta \tau }} .\int\limits_{{\tau_{j} }}^{{\tau_{j + 1} }} {\left( {\tau_{n + 1} - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau_{n + 1}^{\beta - 1} }}{B}} .\left( {\frac{{\Phi_{i + 1}^{n + 1} - 2\Phi_{i}^{n + 1} + \Phi_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }}} \right),$$
(40)
$$\frac{{\left( {\Delta \tau } \right)^{ - \alpha } }}{{\Gamma \left( {2 - \alpha } \right)}}.\sum\limits_{j = 0}^{n} {\frac{{\Phi_{i}^{j + 1} - 2\Phi_{i}^{j} }}{\Delta \tau }} .\left( {\left( {n - j + 1} \right)^{1 - \alpha } - \left( {n - j} \right)^{1 - \alpha } } \right) = \frac{\beta }{B}\left( {\left( {n + 1} \right)\Delta \tau } \right)^{\beta - 1} \left( {\frac{{\Phi_{i + 1}^{n + 1} - 2\Phi_{i}^{n + 1} + \Phi_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }}} \right).$$
(41)

## Solution of momentum equation

From Eq. (25), one can get the following fractal-fractional form of momentum equation:

$${}^{FFP}D_{\tau }^{\alpha ,\beta } w\left( {\xi ,\tau } \right) = \frac{1}{{\text{Re}}}\left( {P + \frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} - Hw\left( {\xi ,\tau } \right) + Gr\Theta \left( {\xi ,\tau } \right) + Gm\Phi \left( {\xi ,\tau } \right)} \right),$$
(42)

Equation (42) can be written as:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\frac{d}{d\tau }\int\limits_{0}^{\tau } {w\left( {\xi ,\tau } \right)\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{{\text{Re}}}} .\left( \begin{gathered} P + \frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} \hfill \\ - Hw\left( {\xi ,\tau } \right) + Gr\Theta \left( {\xi ,\tau } \right) + Gm\Phi \left( {\xi ,\tau } \right) \hfill \\ \end{gathered} \right),$$
(43)

here

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{\tau } {\frac{{\partial w\left( {\xi ,t} \right)}}{\partial t}.\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{{\text{Re}}}} .\left( \begin{gathered} P + \frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} \hfill \\ - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} - Hw\left( {\xi ,\tau } \right) \hfill \\ + Gr\Theta \left( {\xi ,\tau } \right) + Gm\Phi \left( {\xi ,\tau } \right) \hfill \\ \end{gathered} \right) - \frac{{\tau^{ - \alpha } w\left( {\xi ,0} \right)}}{{\Gamma \left( {1 - \alpha } \right)}},$$
(44)

using the given initial condition from Eq. (22), Eq. (44) reduces to the following form:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{\tau } {\frac{{\partial w\left( {\xi ,t} \right)}}{\partial t}.\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{{\text{Re}}}} .\left( \begin{gathered} P + \frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} \hfill \\ - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} - Hw\left( {\xi ,\tau } \right) \hfill \\ + Gr\Theta \left( {\xi ,\tau } \right) + Gm\Phi \left( {\xi ,\tau } \right) \hfill \\ \end{gathered} \right),$$
(45)

by discretizing the above equation at $$\left( {\xi_{i} ,\,\tau = \tau_{n + 1} } \right),$$ the discretized form is given as under:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{{\tau_{n + 1} }} {\frac{{\partial w\left( {\xi_{i} ,t} \right)}}{\partial t}.\left( {\tau_{n + 1} - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau_{n + 1}^{\beta - 1} }}{{\text{Re}}}} .\left( \begin{gathered} P + \frac{{w_{i + 1}^{n + 1} - 2w_{i}^{n + 1} + w_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }} \hfill \\ - \eta_{A} \frac{{w_{i + 2}^{n + 1} - 4w_{i + 1}^{n + 1} + 6w_{i}^{n + 1} - 4w_{i - 1}^{n + 1} + w_{i - 2}^{n + 1} }}{{\left( {\Delta \xi } \right)^{4} }} \hfill \\ - Hw_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) + Gr\Theta_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) \hfill \\ + Gm\Phi_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) \hfill \\ \end{gathered} \right),$$
(46)

from the above step the following result is obtained:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\sum\limits_{j = 0}^{n} {\frac{{w_{i}^{j + 1} - 2w_{i}^{j} }}{\Delta \tau }} .\int\limits_{{\tau_{j} }}^{{\tau_{j + 1} }} {\left( {\tau_{n + 1} - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau_{n + 1}^{\beta - 1} }}{{\text{Re}}}} .\left( \begin{gathered} P + \frac{{w_{i + 1}^{n + 1} - 2w_{i}^{n + 1} + w_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }} \hfill \\ - \eta_{A} \frac{{w_{i + 2}^{n + 1} - 4w_{i + 1}^{n + 1} + 6w_{i}^{n + 1} - 4w_{i - 1}^{n + 1} + w_{i - 2}^{n + 1} }}{{\left( {\Delta \xi } \right)^{4} }} \hfill \\ - Hw_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) + Gr\Theta_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) \hfill \\ + Gm\Phi_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) \hfill \\ \end{gathered} \right),$$
(47)
$$\begin{gathered} \frac{{\left( {\Delta \tau } \right)^{ - \alpha } }}{{\Gamma \left( {2 - \alpha } \right)}}.\sum\limits_{j = 0}^{n} {\frac{{w_{i}^{j + 1} - 2w_{i}^{j} }}{\Delta \tau }} .\left( {\left( {n - j + 1} \right)^{1 - \alpha } - \left( {n - j} \right)^{1 - \alpha } } \right) \hfill \\ = \frac{\beta }{{\text{Re}}}\left( {\left( {n + 1} \right)\Delta \tau } \right)^{\beta - 1} \left( \begin{gathered} P + \frac{{w_{i + 1}^{n + 1} - 2w_{i}^{n + 1} + w_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }} - \eta_{A} \frac{{w_{i + 2}^{n + 1} - 4w_{i + 1}^{n + 1} + 6w_{i}^{n + 1} - 4w_{i - 1}^{n + 1} + w_{i - 2}^{n + 1} }}{{\left( {\Delta \xi } \right)^{4} }} \hfill \\ - Hw_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) + Gr\Theta_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) + Gm\Phi_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) \hfill \\ \end{gathered} \right). \hfill \\ \end{gathered}$$
(48)

## Limiting case

In this section the present obtained solutions are reduced to already published work in order to verify the obtained solutions. Therefore, by putting $$\left( {Gr = 0} \right),$$$$\left( {Gm = 0} \right),$$$$P = 0$$ and $$\frac{1}{K} \to 0$$ present solutions reduced to the solutions recently obtained by Akgül and Siddique12 which verify the present results.

Using the above assumptions Eq. (25) reduces to the following form:

$${}^{C}D_{\tau }^{\alpha ,\beta } w\left( {\xi ,\tau } \right) = \frac{1}{{\text{Re}}}\left( {\frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} - Mw\left( {\xi ,\tau } \right)} \right),$$
(49)

Equation (49) can be written as:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\frac{d}{d\tau }\int\limits_{0}^{\tau } {w\left( {\xi ,\tau } \right)\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{{\text{Re}}}} .\left( {\frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} - Mw\left( {\xi ,\tau } \right)} \right),$$
(50)

here

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{\tau } {\frac{{\partial w\left( {\xi ,t} \right)}}{\partial t}.\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{{\text{Re}}}} .\left( \begin{gathered} \frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} \hfill \\ - Mw\left( {\xi ,\tau } \right) \hfill \\ \end{gathered} \right) - \frac{{\tau^{ - \alpha } w\left( {\xi ,0} \right)}}{{\Gamma \left( {1 - \alpha } \right)}},$$
(51)

using the given initial condition from Eq. (22), Eq. (51) gives to the following form:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{\tau } {\frac{{\partial w\left( {\xi ,t} \right)}}{\partial t}.\left( {\tau - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau^{\beta - 1} }}{{\text{Re}}}} .\left( {\frac{{\partial^{2} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{2} }} - \eta_{A} \frac{{\partial^{4} w\left( {\xi ,\tau } \right)}}{{\partial \xi^{4} }} - Mw\left( {\xi ,\tau } \right)} \right),$$
(52)

by discretizing the above equation at $$\left( {\xi_{i} ,\,\tau = \tau_{n + 1} } \right),$$ one can get the following form:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\int\limits_{0}^{{\tau_{n + 1} }} {\frac{{\partial w\left( {\xi_{i} ,t} \right)}}{\partial t}.\left( {\tau_{n + 1} - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau_{n + 1}^{\beta - 1} }}{{\text{Re}}}} .\left( \begin{gathered} \frac{{w_{i + 1}^{n + 1} - 2w_{i}^{n + 1} + w_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }} \hfill \\ - \eta_{A} \frac{{w_{i + 2}^{n + 1} - 4w_{i + 1}^{n + 1} + 6w_{i}^{n + 1} - 4w_{i - 1}^{n + 1} + w_{i - 2}^{n + 1} }}{{\left( {\Delta \xi } \right)^{4} }} \hfill \\ - Mw_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) \hfill \\ \end{gathered} \right),$$
(53)

From the above equation the following result is obtaiend:

$$\frac{1}{{\Gamma \left( {1 - \alpha } \right)}}.\sum\limits_{j = 0}^{n} {\frac{{w_{i}^{j + 1} - 2w_{i}^{j} }}{\Delta \tau }} .\int\limits_{{\tau_{j} }}^{{\tau_{j + 1} }} {\left( {\tau_{n + 1} - t} \right)^{ - \alpha } dt = \frac{{\beta .\tau_{n + 1}^{\beta - 1} }}{{\text{Re}}}} .\left( \begin{gathered} \frac{{w_{i + 1}^{n + 1} - 2w_{i}^{n + 1} + w_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }} \hfill \\ - \eta_{A} \frac{{w_{i + 2}^{n + 1} - 4w_{i + 1}^{n + 1} + 6w_{i}^{n + 1} - 4w_{i - 1}^{n + 1} + w_{i - 2}^{n + 1} }}{{\left( {\Delta \xi } \right)^{4} }} \hfill \\ - Mw_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) \hfill \\ \end{gathered} \right),$$
(54)
$$\begin{gathered} \frac{{\left( {\Delta \tau } \right)^{ - \alpha } }}{{\Gamma \left( {2 - \alpha } \right)}}.\sum\limits_{j = 0}^{n} {\frac{{w_{i}^{j + 1} - 2w_{i}^{j} }}{\Delta \tau }} .\left( {\left( {n - j + 1} \right)^{1 - \alpha } - \left( {n - j} \right)^{1 - \alpha } } \right) \hfill \\ = \frac{\beta }{{\text{Re}}}\left( {\left( {n + 1} \right)\Delta \tau } \right)^{\beta - 1} \left( \begin{gathered} \frac{{w_{i + 1}^{n + 1} - 2w_{i}^{n + 1} + w_{i - 1}^{n + 1} }}{{\left( {\Delta \xi } \right)^{2} }} \hfill \\ - \eta_{A} \frac{{w_{i + 2}^{n + 1} - 4w_{i + 1}^{n + 1} + 6w_{i}^{n + 1} - 4w_{i - 1}^{n + 1} + w_{i - 2}^{n + 1} }}{{\left( {\Delta \xi } \right)^{4} }} \hfill \\ - Mw_{i}^{n + 1} \left( {\xi_{i} ,\tau_{n + 1} } \right) + \hfill \\ \end{gathered} \right). \hfill \\ \end{gathered}$$
(55)

## Nusselt number, Sherwood number and Skin friction:

### Nusselt number

Mathematically, Nusselt number for CSF can be written as:

$$Nu = - \left. {\frac{\partial \Theta }{{\partial \xi }}} \right|_{\xi = 0}$$
(56)

### Sherwood number:

Mathematically, Sherwood number for CSF can be written as:

$$Sh = - \left. {\frac{\partial \Phi }{{\partial \xi }}} \right|_{\xi = 0}$$
(57)

### Skin friction

Skin friction for CSF is as under:

$$Sf\left( {\xi ,\tau } \right) = \left( {\frac{\partial w}{{\partial \xi }} - \eta \frac{{\partial^{3} w}}{{\partial \xi^{3} }}} \right)$$
(58)

As the given flow model is between two parallel plates. Therefore, the skin friction at the left and right plates is as under:

$$Sf_{lp} \left( {0,\tau } \right) = \left( {\frac{\partial w}{{\partial \xi }} - \eta \frac{{\partial^{3} w}}{{\partial \xi^{3} }}} \right)_{\xi = 0}$$
(59)
$$Sf_{rp} \left( {1,\tau } \right) = \left( {\frac{\partial w}{{\partial \xi }} - \eta \frac{{\partial^{3} w}}{{\partial \xi^{3} }}} \right)_{\xi = 1}$$
(60)

where $$Sf_{lp} (.)$$ and $$Sf_{rp} (.)$$ denotes the skin friction at left and right plates respectively.

## Results and discussion

This section provide fractal-fractional derivative model of unsteady MHD generalized Couette flow of CSF in channel with embedded in porous media with power law kernel. Numerical solutions for the proposed problem are obtained using the implicit finite difference method. In this study we have found the influence of fractal dimension and fractional operator on the fluid motion, fluid temeprature and concentration graphically. Furthermore, for clear understanding the influence of all parameters is shown through graphs which effect the fluid motion, temperature and concentration.

Figure 1 shows the physical sketch of the proposed problem. The effect of fractal dimension $$\beta$$ on velocity profile is highlighted in Fig. 2. From the graph it is quite clear that for greater values of fractal parameter result a decay in the fluid velocity it is due to the power law kernel. The effect of fractional parameter $$\alpha$$ is shown in Fig. 3. From the figure a decreasing in the CSF velocity is noticed. The influence of $$\alpha$$ on CSF velocity is similar to the effect of fractal dimension $$\beta$$ on the velocity field. The comparison between fractal-fractional CSF velocity and fractional velocity is plotted in Fig. 4. From the figure it can be noticed that the magnitude of fractional velocity is greater than fractal-fractional velocity. In this paper we have added a parameter $$\beta$$ known as fractal-fractional dimension. This parameter $$\beta$$ shows the combined effect of fractal-fractional derivative with fractional derivative. By introducing $$\beta$$ in CSF flow can explain better memory effect compared to fractional and classical CSF fluid. The effect of Gr and Gm on CSF velocity is highlighted in Figs. 5 and 6 respectively. From the graphs it is very clear that for greater values of Gr and Gm result an increase in the CSF velocity. This increase in the fluid velocity is due to the fact that Gr and Gm are the ratios of boyancy forces and viscous forces which are responsible to accelerate the fluid motion as a result CSF velocity increases. The effect of magnetic parameter $$M$$ is highlighted if Fig. 7. From the figure it is clear that increasing $$M$$ result a decrease in the fractal-fractional CSF velocity. This is due to the fact that for greater values of $$M$$ Lorentz forces increases in the CSF which control the boundary layer thickness as a result velocity of the fractal-fractional CSF decreases. The influence of porosity $$K$$ on fractal-fractional CSF velocity is highlighted in Fig. 8. From this figure it seems that the CSF velocity increases with the greater values of $$K$$ it is due to the fact that increasing $$K$$ increases the pores in the channel as a result the fluid velocity accelerates.

The effect of couple stress parameter $$\eta_{A}$$ is depicted in Fig. 9. This figure shows the influence of $$\eta_{A}$$ on the fractal-fractional velocity in channel. Increasing the values of couple stress parameter $$\eta_{A}$$ increases the viscosity of the fluid as a result the CSF velocity retards in channel. From this figure it can also be noticed that for $$\eta_{A} = 0$$ shows the comparison of simple Newtonian viscous fluid with fractal-fractional velocity in channel. The influence of external pressure $$P$$ is highlighted in Fig. 10. From the figure it can be noticed that increasing the values of $$P$$ result an increase in the CSF velocity in channel.

The comparison of fractal-fractional temperature with fractal temperature $$\beta$$ and fractional temperature $$\alpha$$ is highlighted in Fig. 11. From the comparison we can see that the classical temperature is higher than fractal-fractional, fractal and fractional temperature. The influence of fractional parameter $$\alpha$$ and fractal parameter $$\beta$$ on temperature is depicted in Figs. 12 and 13 respectively. From both the figures it can be noticed that for the greater values of fractional parameter $$\alpha$$ and fractal parameter $$\beta$$ the temeprature of the CSF in channel reduces. The effect of $$\Pr$$ on temperature profile is highlighted in Fig. 14. From the figure we can observe that for larger values of $$\Pr$$ the CSF temperature decreases it is due to the fact that increasing $$\Pr$$ results a decrease in thermal conductivity of the fluid as a result the temperature of the fluid decreases. The influence of Reynolds umber on temeprature profile is shown in Fig. 15. From the figure it can be observed that greater values of Reynolds number decrease the temperature of CSF in channel.

The comparison of fractal-fractional, fractal and fractional concentration is highlighted in Fig. 16. From the figure one can noticed that the magnitude of classical concentration is higher than the concentration for fractal and fractional concentration. The effect of fractional parameter $$\alpha$$ and fractal parameter $$\beta$$ on concentration profile is highlighted in Figs. 17 and 18 respectively. From both the figures it can be noticed that for the greater values of fractional parameter $$\alpha$$ and fractal parameter $$\beta$$ in both the cases the concentration profile of the CSF in channel reduces. The influence of Reynolds number $${\text{Re}}$$ on concentration profile is highlighted in Fig. 19. From this figure the greater values of Reynolds number decreases the concentration of the CSF in channel it is due to the fact that greater Reynolds number decreases the CSF motion as a result the concentration of the fluid decreases. The effect of Schmidth number Sc on concentration profile is highlighted in Fig. 20. Form the figure it can be notice that increasing the values of Sc result a decrease in the CSF concentration profile. The comparison between the present solution to the already published result is highlighted in Fig. 21. From the figure it can be seen that by putting (Gr = 0), (Gm = 0), P = 0 and $$\frac{1}{K} \to 0$$ our solutions reduced to the solutions recently obtained by Akgul and Siddique12 which verify our obtained solutions.

The skin friction for left and right plate are evaluated and presented in Tables 1 and 2 respectively. From these tables skin friction variation can be noticed for varying different parameters. The bold values in the tables show the change in skin friction in that specific parameter. Similarly, Table 3 shows the Nusselt number variation for different parameters. Table 4 shows the variation in Sherwood number for different parameters. In the Tables 1, 2, 3 and 4 bold values represents the changes in the specific parameter and its effect on the ski friction, Nusselt number and Sherwood number.

## Conclusion

The present paper is focused to study the applications of fractal and fractional derivative on the unsteady MHD CSF in channel with power law kernel. The fractal-fractional CSF is assumed to flow in channel embedded in porous medium. The unsteady CSF with heat and mass transfer passes through the channel in the presence of external pressure. The implicit finite difference method is applied to obtain the numerical solutions of the proposed fractal-fractional CSF model of generalized Couette flow. During the analysis of the present study the following results are obtaiend.

• The velocity, temperature and concentration of the magnetohydrodynamics CSF in porous channel decreases for the greater values of both fractional parameter $$\alpha$$ and fractal parameter $$\beta$$ respectively.

• In the fluid dynamics the fractal-fractional model of CSF explain good memory effect as compared to fractional model of CSF.

• CSF velocity in porous channel decrease with the increase in $$\alpha$$,$$\beta$$,$$M$$ and $$\eta_{A}$$ .

• CSF velocity in porous channel increases with the higher values of $$Gr$$,$$Gm$$ and $$K$$.

• The temperature of the CSF is decreases with the higher values of $$\alpha ,$$$$\beta$$,$$\Pr$$ and $${\text{Re}}$$.

• The concentration of the CSF is decreases with the higher values of $$\alpha ,$$$$\beta$$,$$Sc$$ and $${\text{Re}}$$.

## Data availability

Data is available upon reasonable request to corresponding author.

## References

1. 1.

Tenreiro Machado, J. A. et al. Some applications of fractional calculus in engineering. Math. Probl. Eng https://doi.org/10.1155/2010/639801 (2010).

2. 2.

Dalir, M. & Bashour, M. Applications of fractional calculus. Appl. Math. Sci 4(21), 1021–1032 (2010).

3. 3.

Tavazoei, M. S., Haeri, M., Jafari, S., Bolouki, S. & Siami, M. Some applications of fractional calculus in suppression of chaotic oscillations. IIEEE Trans. Ind. Electron 55(11), 4094–4101. https://doi.org/10.1109/TIE.2008.925774 (2008).

4. 4.

Sabatier, J. A. T. M. J., Agrawal, O. P., & Machado, J. T. Fract. Calc. Appl. Anal, (Vol. 4, No. 9). Dordrecht: Springer, (2007).

5. 5.

Koca, I. & Atangana, A. Solutions of Cattaneo-Hristov model of elastic heat diffusion with Caputo-Fabrizio and Atangana-Baleanu fractional derivatives. J. Therm. Sci 21(6), 2299–2305. https://doi.org/10.2298/TSCI160209103K (2017).

6. 6.

Atangana, A., & Baleanu, D. New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model. arXiv preprint arXiv: 1602.03408, (2016).

7. 7.

Arif, M., Ali, F., Sheikh, N. A., Khan, I. & Nisar, K. S. Fractional model of couple stress fluid for generalized Couette flow: A comparative analysis of Atangana-Baleanu and Caputo-Fabrizio fractional derivatives. IEEE Access 7, 88643–88655 (2019).

8. 8.

Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications (Elsevier, 1998).

9. 9.

Abdeljawad, T. & Baleanu, D. On fractional derivatives with generalized Mittag-Leffler kernels. Adv. Differ. Equ. 2018(1), 1–15 (2018).

10. 10.

Atangana, A. Fractal-fractional differentiation and integration: Connecting fractal calculus and fractional calculus to predict complex system. Chaos Solitons Fractals 102, 396–406 (2017).

11. 11.

Ali, Z., Rabiei, F., Shah, K. & Abdul Majid, Z. Dynamics of SIR mathematical model for COVID-19 outbreak in Pakistan under Fractal-fractional derivative. Fractals 2, 1 (2021).

12. 12.

Akgül, A. & Siddique, I. Novel applications of the magnetohydrodynamics couple stress fluid flows between two plates with fractal-fractional derivatives. Numer. Methods Partial Differ. Equ. 37(3), 2178–2189 (2021).

13. 13.

Esmonde, H. Fractal and fractional derivative modelling of material phase change. Fractal Fract. 4(3), 46 (2020).

14. 14.

Akgül, A. Analysis and new applications of fractal fractional differential equations with power law kernel. Discrete & Continuous Dynamical Systems-S, (2020).

15. 15.

Muzychka, Y. S. & Edge, J. Laminar non-Newtonian fluid flow in noncircular ducts and microchannels. J. Fluids Eng 130, 11 (2008).

16. 16.

Chhabra, R. P. & Richardson, J. F. Non-Newtonian flow in the process industries: fundamentals and engineering applications (Butterworth-Heinemann, 1999).

17. 17.

Pawar, S. S. & Sunnapwar, V. K. Experimental studies on heat transfer to Newtonian and non-Newtonian fluids in helical coils with laminar and turbulent flow. Exp. Therm. Fluid Sci 44, 792–804 (2013).

18. 18.

Hoyt, J. W. Some applications of non-Newtonian fluid flow. In Rheology Series Vol. 8 797–826 (Elsevier, 1999).

19. 19.

Sohail, M., Naz, R., Shah, Z., Kumam, P. & Thounthong, P. Exploration of temperature dependent thermophysical characteristics of yield exhibiting non-Newtonian fluid flow under gyrotactic microorganisms. AIP Adv. 9(12), 125016 (2019).

20. 20.

Bao, K., Lavrov, A., & Nilsen, H. M. Numerical Modelling of Non-newtonian Fluid Flow in Fractures and Porous Media. In ECMOR XV-15th European Conference on the Mathematics of Oil Recovery (pp. cp-494). European Association of Geoscientists & Engineers, (2016, August).

21. 21.

Gowda, R. P. et al. Thermophoretic particle deposition in time-dependent flow of hybrid nanofluid over rotating and vertically upward/downward moving disk. Surf. Interfaces 22, 100864 (2021).

22. 22.

Gowda, R. P., Rauf, A., Kumar, R. N., Prasannakumara, B. C. & Shehzad, S. A. Slip flow of Casson-Maxwell nanofluid confined through stretchable disks. Indian J. Phys 2, 1–9 (2021).

23. 23.

Zhou, S. S. et al. Nonlinear mixed convective Williamson nanofluid flow with the suspension of gyrotactic microorganisms. Int. J. Mod. Phys. B 35(12), 2150145 (2021).

24. 24.

Stokes, V. K. Couple stresses in fluids. In Theories of Fluids with Microstructure 34–80 (Springer, 1984).

25. 25.

Arif, M., Ali, F., Khan, I. & Nisar, K. S. A time fractional model with non-singular kernel the generalized Couette flow of couple stress nanofluid. IEEE Access 8, 77378–77395 (2020).

26. 26.

Adesanya, S. O., Souayeh, B., Rahimi-Gorji, M., Khan, M. N. & Adeyemi, O. G. Heat irreversibiility analysis for a couple stress fluid flow in an inclined channel with isothermal boundaries. J. Taiwan Inst Chem Eng 101, 251–258 (2019).

27. 27.

Krishna, M. V. & Chamkha, A. J. MHD peristaltic rotating flow of a couple stress fluid through a porous medium with wall and slip effects. Spec. Top. Rev. Porous Media 10, 3 (2019).

28. 28.

Khan, N. A., Mahmood, A. & Ara, A. Approximate solution of couple stress fluid with expanding or contracting porous channel. Eng. Comput. 2, 1 (2013).

29. 29.

Akhtar, S. Flows between two parallel plates of couple stress fluids with time-fractional Caputo and Caputo-Fabrizio derivatives. Eur. Phys. J. Plus 131(11), 1–13 (2016).

30. 30.

Kumar, R. N. et al. Impact of magnetic dipole on ferromagnetic hybrid nanofluid flow over a stretching cylinder. Phys. Scr 96(4), 045215 (2021).

31. 31.

Gowda, R. P., Kumar, R. N., Rauf, A., Prasannakumara, B. C. & Shehzad, S. A. Magnetized flow of sutterby nanofluid through cattaneo-christov theory of heat diffusion and stefan blowing condition. Appl. Nanosci 2, 1–10 (2021).

32. 32.

Gowda, R. P., Kumar, R. N., Prasannakumara, B. C., Nagaraja, B. & Gireesha, B. J. Exploring magnetic dipole contribution on ferromagnetic nanofluid flow over a stretching sheet: An application of Stefan blowing. J. Mol. Liq 335, 116215 (2021).

33. 33.

Yusuf, T. A., Mabood, F., Prasannakumara, B. C. & Sarris, I. E. Magneto-bioconvection flow of williamson nanofluid over an inclined plate with gyrotactic microorganisms and entropy generation. Fluids 6(3), 109 (2021).

34. 34.

Kumar, R. N. et al. Impact of magnetic dipole on thermophoretic particle deposition in the flow of Maxwell fluid over a stretching sheet. J. Mol. Liq 334, 116494 (2021).

35. 35.

Naveen Kumar, R., Punith Gowda, R. J., Prasanna, G. D., Prasannakumara, B. C., Nisar, K. S., & Jamshed, W. Comprehensive study of thermophoretic diffusion deposition velocity effect on heat and mass transfer of ferromagnetic fluid flow along a stretching cylinder. P I Mech Eng E-J Pro, 09544089211005291, (2021).

36. 36.

Punith Gowda, R. J., Naveen Kumar, R., Jyothi, A. M., Prasannakumara, B. C. & Sarris, I. E. Impact of binary chemical reaction and activation energy on heat and mass transfer of marangoni driven boundary layer flow of a non-Newtonian nanofluid. Processes 9(4), 702 (2021).

37. 37.

Gowda, R. P. et al. Computational modelling of nanofluid flow over a curved stretching sheet using Koo-Kleinstreuer and Li (KKL) correlation and modified Fourier heat flux model. Chaos, Solitons Fractals 145, 110774 (2021).

38. 38.

Li, Y. X. et al. Dynamics of aluminum oxide and copper hybrid nanofluid in nonlinear mixed Marangoni convective flow with entropy generation: Applications to renewable energy. Chin. J. Phys 2, 1 (2021).

39. 39.

Xiong, P. Y. et al. Dynamics of multiple solutions of Darcy-Forchheimer saturated flow of Cross nanofluid by a vertical thin needle point. Eur. Phys. J. Plus 136(3), 1–22 (2021).

40. 40.

Ali, F., Ahmad, Z., Arif, M., Khan, I. & Nisar, K. S. A time fractional model of generalized Couette flow of couple stress nanofluid with heat and mass transfer: Applications in engine oil. IEEE Access 8, 146944–146966. https://doi.org/10.1109/ACCESS.2020.3013701 (2020).

## Acknowledgements

The authors are thankful to the editors and reviewers for their constructive comments to improve the manuscript. The authors acknowledge the financial support provided by the Center of Excellence in Theoretical and Computational Science (TaCS-CoE), KMUTT. Moreover, this research project is supported by Thailand Science Research and Innovation (TSRI) Basic Research Fund: Fiscal year 2021 under project number 64A306000005.

## Author information

Authors

### Contributions

M.A. and P.K. formulated the problem. M.A. solved the problem and wrote the whole manuscript. P.K. performed the investigation. W.K, A.A. and T.S. review the final draft.

### Corresponding author

Correspondence to Poom Kumam.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

### Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Arif, M., Kumam, P., Kumam, W. et al. Analysis of newly developed fractal-fractional derivative with power law kernel for MHD couple stress fluid in channel embedded in a porous medium. Sci Rep 11, 20858 (2021). https://doi.org/10.1038/s41598-021-00163-3

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s41598-021-00163-3