A time fractional model of Brinkman-type nanofluid with ramped wall temperature and concentration

Nanofluid is an innovative heat transfer fluid with the potential to significantly enhance the heat transfer performance of traditional fluids. By adding various types of nanoparticles to ordinary base fluids, several attempts have been made to boost the rate of heat transfer and thermal conductivity. The unsteady electrically conducting flow of Brinkman-type nanofluid over an infinite vertical plate with ramping wall temperature and concentration is investigated in this article. Water is taken as the base fluid, and multi-walled carbon nanotubes are distributed equally throughout it. The Caputo-Fabrizio fractional derivative, which has a non-singular kernel, is used to generalize the classical model. The Laplace transform technique has been utilized to achieve exact solutions. Furthermore, various graphs for fractional and physical parameters are used to represent the solutions. All figures are drawn for both conditions, that is, ramped and isothermal wall temperature and concentration. The velocity field increases for greater values of thermal and mass Grashof numbers while the reverse effect is observed for Hartman number, Brinkman parameter and volume fraction. Moreover, the obtained results are also reduced to the already published results in order to show the validation of the present results. The results are used to calculate the skin friction, Nusselt number, and Sherwood number. The heat transfer of pure water is increased by 17.03% when 4% of nanoparticles are added to it which will of course increase the efficiency of solar collectors and solar pools. Moreover, the mass transfer decreases by 3.18% when 4% of nanoparticles which are dispersed in it.


Introduction
Nanotechnology has grown rapidly in recent decades in a variety of disciplines of engineering and science, including electronics and power generation, where heat transmission is a common occurrence. The additional technological progress need specialized heat transfer management. Many researchers have attempted to study an effective heat transfer medium for this purpose. Many techniques have been explored but failed channel and extending the surface area. For the enhancement of the heat transfer rates, many experiments characterized and examined different nanofluids. Nanofluids are a combination of nanoparticles and traditional base fluids. Water, motor oil, ethylene, and crude oil are examples of base fluids. Molybdenum disulfide, aluminum, graphite, silver, copper, and other nanoparticles dispersed in the basic fluids. Scientists have been attempting to solve this problem since the beginning to improve the characteristics of various kinds of fluids. They are trying to improve the deficiencies of these regular fluids. The nanoparticles dispersion in the base fluid had been discovered for the first time by Choi and Eastman. 1 Various Nanofluids in heat transfer applications have been described, produced, and evaluated by several researchers. 2,3 For the combination of nanoparticles with base fluids, several theoretical and experimental research are conducted. In base fluids including oil, water, and concrete, nanoparticles such as metallic oxides, carbide ceramics, nitride ceramics, and carbon in different forms are widely employed. However, because to their high thermal conductivity, some recent research suggest that carbon nanotubes (CNTs) are better suitable for heat transfer rates. 4,34 The heat conductivity of CNTs nanoparticles is high, while their density is low.
A thermal energy transfer from a greater concentration to a lower concentration is referred to as ''transfer of heat.'' Mass transfer is the net movement of mass from one location to another. Conduction, convection, and radiation are all ways for transferring heat. The combined effect of bouncy thermal diffusion causes the mass and heat transfer are also being studied at the same time. Cooling processes, solar collectors, food processing, and many other lubricants such as engine oil are all examples of this phenomenon in modern technology and industries. An isothermal process is a thermodynamic process that occurs at a constant temperature. Heat transfer can be seen in this step since the temperature remains constant. An isothermal process is one in which the amount of work done is due to the system's net heat content. We can say that an isothermal process is a constant temperature process. Isothermal processes include changes in the state or phase of various fluids due to melting and evaporation. Thermal equilibrium is maintained because heat is transferred in and out of the system at a slow rate. The heat produced by a device is referred to as ''thermal.'' The word ''isothermal'' refers to the concept of ''equal heat,'' which is defined as thermal equilibrium. Isothermal means ''equal heat,'' as ''iso'' means ''equal.'' Different researchers investigated the heat transfer in fluids by considering the isothermal process. Akram and Nadeem 5 examined the effects of a magnetic field and heat transfer on two-dimensional Jaffrey fluid peristaltic motion in an asymmetric channel. The effect of heat radiation along with ramped wall boundary conditions on the fluid flow was studied by Anwar et al. 6 Khan et al. 7 investigated free convection Maxwell fluid flow with ramping wall temperature and concentration. Through an infinite isothermal vertical plate, Soundergekar 8 calculated the exact solution to the viscous flow problem.
When it comes to the flow through a porous material, classical Navier Stokes' equation fails. Several important physical processes in porous media have been explored for transport of mass and heat in the literature. Other variables included the effect of local thermal equilibrium between the solid phases and the fluid, thermal dispersion, porosity variance, and anisotropic porous media. Many scientific investigations have been done through permeable medium utilizing generalized model during the last few decades. The Darcy model is defined in detail by Brinkman-Forch-Homier. Brinkman introduced one of his unique models for flow of fluid among a permeable media. 9,10 As Brinkman explained, when a flow moves through a highly permeable medium, it is referred to as a Brinkman-type fluid flow. Using a magnetic field, Varma and Babu 11 studied the steady flow of an incompressible fluid across a channel with varying permeability medium. Using the Brinkman model, Hsu and Cheng 12 studied fluid flow through a semi-infinite vertical plate immersed in porous material. Rajagopal 13 looked on fluid flow through a porous medium and included the influence of the porous media in the governing equations. Vafai and Tien 14,15 proposed a generalized transport model. In addition, a lot of research is reported on the effects of inertial forces and solid walls through a porous media. Vafai 16 and Hadim and Vafai 17 wrote a book chapter on the importance aspects of porous media transport phenomena. Vafai and Hadim 18 and So¨zen and Vafai 19 done extensive research in porous media on multiphase transport, and models were developed. These investigations initiated and explored a systematic report of several various phase models with the phase change process. Alazmi and Vafai 20 show the differences between porous media's heat transfer and flow transport models. Vafai et al. 21 examined the literature for a number of turbulent models in permeable medium. Porous media theory is also frequently applied in biological applications. Khanafer and Vafai 22 looked into turbulent flows in porous medium with drug delivery and magnetic particles movement. Khaled and Vafai 23 studied the diffusion process in brain tissues, blood movement in cancers, and bioheat transport in tissues using the porous media.
Many physical problems are outside the reach of simple classical models to explain. Fractional calculus was developed to investigate these problems. The generalized form of classical calculus is fractional calculus.
After Leibniz introduced the definition of the nth order derivative, fractional calculus was developed. Del Hospital asked Leibniz what would happen if fractional-order is considered. 24 Following that, several scientists started to think about it, proposing alternative meanings of fractional derivatives. Fractional calculus has practical applications in modern technology electrochemistry, electromagnetism, electric circuits, voltage dividers, neuron models, 3-D chaotic structures, and chaotic circuits, heat and mass transfer phenomenon, mathematical epidemiology. [25][26][27][28][29][30] It was stated by Leibnitz in a 1695 letter to L. Hospital. 31 During the first part of the nineteenth century, systematic investigations on fractional calculus were conducted. 32 Fractional derivatives are defined in many ways in the literature and applied to various sciences, technology, and engineering sectors to explore realworld applications. 33 From open literature, it is cleared that derivatives of integer orders are local in nature while derivatives in the fractional form are non-local, which retain memory property. 34,35 To the best of the author's knowledge, the literature upon the nanofluid flow problems with the effect of the ramped and isothermal wall boundary conditions is very rare. Therefore, based on the above literature, the present study is carried out for the MHD flow of Brinkman type nanofluid with the joint effect ramped and isothermal wall boundary conditions on temperature and concentration. Water is considered as a base fluid while MWCNTs are equally dispersed in the base fluid in order to enhance its performance regarding it's applications. Furthermore, to generalize the classical model, Caputo-Fabrizio (CF) time fractional time has been used to convert the classical Brinkman model with energy and concentration equations into a fractional form. With the governing momentum equation, the generalized energy and concentration equations are coupled in the momentum equation. For the exact solutions of coupled system of PDEs, the Laplace transform technique has been used. Additionally, the present results are also validated and reduced to the already published work. The embedded parameters have been illustrated using a variety of graphs with physical explanation.

Mathematical formulation
The current problem considers the unsteady, laminar, and MHD flow of Brinkman-type nanofluid through an infinite plate. The fluid is thought to be electrically conductive. The impact of the magnetic field is chosen along the y-axis, which is perpendicular to the flow direction, while the fluid flow is examined along the xaxis as shown in Figure 1. Both the plate and the fluid are initially at rest for t ł 0, with a constant temperature T ' and concentration C ' . The plate starts oscillating with H t ð Þ cos vt ð Þ at t = 0 + , and the temperature and concentration of the plate have grown to T ' + T w À T ' ð Þ t=t 0 and C ' + C w À C ' ð Þ t=t 0 , respectively. The governing equations for the Brinkman type fluid flow along with transfer of mass and heat are given by 34 : Here, r, c p , k, D,b,Ĩ o ,Ṽ , T , C and are the notations that are used for the density, heat capacity, thermal conductivity, mass diffusivity, body forces, interactional forces, velocity vector, temperature, concentration, and Cauchy stress tensor respectively. The following are the fields for velocity, temperature, and concentration which can be seen in Figure 1 given above:Ṽ Using equation (4), equations (1)-(3) in components form for Brinkman type nanofluid flow becomes: where b Ã is known as the Brinkman parameter, m nf is the dynamic viscosity, and the electrical conductivity is s nf , u is the velocity, r nf is the density, uniform magnetic field is represented by B 0 , g is the gravitational acceleration, b T ð Þ nf is the thermal expansion coefficient, where the temperature is represented by T 1 , C 1 is the concentration, heat capacitance is symbolized as C P ð Þ nf , and k nf is the thermal conductivity of nanofluid. The required initial and boundary conditions are as follows 7,34 : Where mathematical expressions for r nf , m nf , s nf , (b T ) nf , (bc) nf , (C p ) nf , and k nf , are given by. 34 Nanofluid, base fluid, and solid CNTs nanoparticles are denoted by the subscripts the subscripts nf , f and s in equation (9). The thermophysical properties of the considered base fluid and nanoparticles are given in Table 1 which is given below. Introducing the following dimensionless variables for nondimensionalization: By using the dimensionless variables given in equation (10) and dropping * notation, with initial and boundary conditions, the dimensionless form of the governing equations becomes: Pr with Table 1. Thermo-physical properties of H 2 O and MWCNTs. 34 Material Here, b, M, Gr, Gm, Pr and Sc is the Brinkman fluid parameter, magnetic number, thermal Grashof number, mass Grashof number, Prandtl number, and Schmidt number, respectively. The dimensionless governing equations (11)-(13) are converted to Caputo-Fabrizio time Fractional model as: Where CF D a t : ð Þ is the Caputo-Fabrizio (CF) timefractional operator defined by 35 : Where N a ð Þ is the normalization function with

Solution of the problem
This section provides the exact solutions of the considered fractional model by using the Laplace technique:

Solutions for temperature field
We get the following results by applying the Laplace transform to equation (16) and utilizing equation (14).
The boundary conditions provided in equation (20) are used to find the solutions of equation (13) in the Laplace transform domain as: where l 0 = x 4 Pr x 5 , g 0 = 1 1Àa , g 1 = g 0 a. It's crucial to note that equation (21) is the Laplace transform domain solution of the energy equation with ramped wall temperature. Equation (21) can be further reduced as follows: where: The inverse Laplace transform is used to convert equation (23) back to its original form, yielding: where H t À 1 ð Þis the Heaviside unit step function. The term Y Ramp h, t ð Þ is obtained as: Where f Ramp h, t, l 0 , g 0 , When 0\t\1, the exact energy solutions for ramped wall temperature are found in equation (26). For the isothermal temperature condition, equation (16) is resolved, yielding: such that The inverse Laplace transform is used to equation (28) to obtain the time domain solutions for isothermal temperature, yielding: The exact solutions of equation (16) for ramped and isothermal temperatures are represented by equations (24) and (29), respectively.

Solution of concentration
By applying the Laplace transform technique in the same way applied in the temperature field, the solution of concentration equation given in equation (17) By inverting the Laplace transform, we get: The term F Ramp (y, t) is obtained as: where When 0\t\1, the exact solutions for ramped wall concentration are found in equation (35). For the isothermal boundary condition, equation (17) is re-solved, yielding: such that F iso h, s ð Þ= g Iso (h, s, g 2 , g 1 ) = 1 s exp Àh ffiffiffiffiffiffiffiffiffiffiffiffiffi g 2 s s + g 1 r : ð37Þ By taking the inverse Laplace transform, we get: Equations (33) and (35) represent the exact solutions of equation (16) for ramped and isothermal concentrations, respectively.

Solution for the velocity field
By applying the Laplace transformation to equation (15) and utilizing equation (14), yields: After further simplification and incorporating equations (21) and (30), we obtain: The transformed boundary conditions are: c h, s ð Þ= s s 2 + v 2 and c(', s) = 0, where l 1 = 1 x 1 Gr. The solution of equation (40) where l 0 g 0 = l 10 , Equation (41) can be rewritten in a more appropriate way as: where Upon inversion of the Laplace transform, equation (42) becomes: where c c h, t ð Þ= e Àh ffiffiffi ffi here l 4 = l 2 À l 1 g 1 l 5 = Gr 0 l 2 + l 3 g 1 l 2 3 , l 7 = l 2 À g 2 l 5 = Gm 0 l 2 + g 2 g 1 l 2 6 , where Ã represents the convolution product. It's worth noting that equation (47) describes velocity field solutions with ramping wall boundary conditions. For isothermal boundary conditions, equation (15) is now re-solved as: To write equation (52) more simplified and convenient form as The terms c c h, s ð Þ, The inverse Laplace transform of equation (46), which has the following form, yields the solution for isothermal temperature: where c 4(iso) h, t ð Þ= e Àh ffiffiffi ffi and

Limiting case
The obtained results given in equations (47) and (55) can be reduced to the results published by Saqib et al. 34 by ignoring mass Grashoff number that is, Gm ! 0. By doing so, it leads us to the following results for both, ramped and isothermal wall boundary conditions: and The results given in equations (57)  Nusselt number, Sherwood number, and skin friction

Nusselt number
The Nusselt number has the following mathematical expression:

Sherwood number
The Sherwood number has the following mathematical expression:

Skin friction
The Skin friction has the mathematical expression as:

Results and discussion
The goal of this study is to see how a Brinkman-type nanofluid flows through a vertical plate with the influence of applied magnetic field. The influence of ramped and isothermal boundary conditions on concentration and temperature is also examined. To generalize the classical model, Caputo-Fabrizio fractional derivative is applied. Additionally, water is considered as a base fluid while MWCNTs are dispersed to enhance the thermal properties of regular water. The Laplace transform technique is used to analyze exact solutions. Temperature, concentration, and velocity field closed form solutions are computed. The physical representation of the present problem can be seen in Figure 1.  Figures 11 and 12 show the effect of concentration distributions. It's worth noting that we use t = 0:5 and t = 1:5 for ramped and isothermal temperatures, respectively. Skin friction, Nusselt number, and Sherwood number are all computed and tabulated. Figure 2 shows the effect of a on the velocity field for both ramping and isothermal wall boundary conditions. Fractional models, in compared to classical models, are more general. Fractional models have many solutions, allowing experimentalists to match their results to theoretical results in excellent agreement, which is not possible in the classical model for a = 1.    The impact of b on the velocity distribution is seen in Figure 3. The drag force is directly proportional to b. The greater the influence of b, the stronger the drag forces become, resulting in a reduction in the velocity field. The effect of f on the velocity field can be seen in Figure 4. When nanoparticles disperse in a conventional fluid, the viscosity of the base fluid increases, slowing fluid velocity. Figure 5 illustrates the influence of M on the velocity field. Fluid velocity drops when M is raised. The physics behind this is that when M increases, Lorentz forces are generated, which oppose fluid motion and create resistive forces on the fluid flow, lowering fluid velocity. Figures 6 and 7 show the effect of Gr and Gm on the velocity field. The velocity of the nanofluid rises in both cases. This velocity pattern is evident because raising Gr and Gm increases buoyancy forces while lowering viscous forces, causing fluid motion to accelerate. Figure 8 shows the effect of a on the temperature field for both ramping and isothermal wall boundary conditions. Fractional models, in compared to classical models, are more general. Fractional models give us several solutions for different values of a, which means that fractional models give a choice to the experimentalists to match their results with best fit solution. The effect of f on the temperature field can be seen in Figure 9. By increasing f, the temperature field clearly increased. As the heat transfer rate of the base fluid increases as the amount of nanoparticles increases, the temperature field rises. Figure 10 depicts the impact of time t on the temperature field. It has been found that in both situations, t raises the fluid's temperature. It is also detected that for ramped boundary condition that is, 0\t\1, the temperature at the boundary changes while for the case of isothermal boundary conditions t ø 1, the temperature reaches to its maximum and remains unchanged.
The effect of a on concentration distribution is shown in Figure 11. It is found that by considering fractional model, it gives us different curves of concentration profile. It also gives us the solutions of concentration profile in classical sense by considering a = 1. Fractional models give a choice to the experimentalists to match their results with best fit solution. The effect    of t on concentration distribution is displayed in Figure 12. It has been observed that for both cases, t increases the concentration field. It is also observed that for ramped boundary conditions that is, 0\t\1, the concentration at the boundary changes, while for the case of isothermal boundary conditions t ø 1, the concentration reaches its maximum and remains unchanged.
For the validation of the present results, we also reduced the present solutions to the solutions of Saqib et al. 34 by taking Gm ! 0 which is given as a limiting case in section 3.4. We also portrayed the comparison of the present and published results by sketching graph. The comparison between the present and published results of Saqib et al. 34 is given in Figure 13. From the figure, it is observed that the present results overlapped with the results of Saqib et al. 34 for both, ramped and isothermal cases when Gm ! 0 which validates the accuracy of our solutions.
The thermophysical properties of the base fluid as well as nanoparticles is given in Table 1. The variation in skin friction against different parameters is tabulated in Table 2. Skin friction is calculated for both the cases, that is, classical as well as the fractional case. Similarly, variation in Nusselt number against f is calculated and displayed in Table 3. It is observed that the dispersion of MWCNTs in regular water can enhance the heat    transfer rate of pure water up to 17.02% by adding 4% of nanoparticles. This means that water can store more heat and hence can enhance the performance of waterbased solar collectors, solar heaters, solar pools, solar panels, etc. A decrease in the mass transfer is portrayed in Table 4. The rate of mass transfer reduces as the number of nanoparticles increases. This is due to the regular fluid's density and viscosity increasing as the number of nanoparticles increases. The table demonstrates that raising 4% of nanoparticles reduces the rate of mass transfer by up to 3.18%.

Conclusion
This research is conducted to determine the closed form solutions for Brinkman-type nanofluid flow over a vertical plate. The effects of ramping and isothermal wall boundary conditions on the fluid's temperature and concentration are examined in the presence of an applied magnetic field. The nanoparticles utilized are MWCNTs, and the base fluid is water. The classical model is then generalized by using the Caputo-Fabrizio fractional derivative, which has recently become the most often used fractional derivative. The solutions of the coupled system is obtained via the Laplace transform technique. The obtained results are also validated by reducing and comparing it with the already published results by making some parameters absent. The graphs also illustrate the results that have been obtained. The following are the main points of the study: A decrease in the nanofluid velocity can be observed for large values of a for the case of ramped wall boundary conditions, while the opposite impact is observed for isothermal wall boundary conditions. The present results are reducable to the classical Brinkman type nanofluid model by taking a ! 1.
The obtained results are also reducable to the already published results of Saqib et al. 34 by taking Gm ! 0. The velocity profile decreases for larger values of M and b.
By raising Gr and Gm, the velocity profile can be enhanced. Increasing f has resulted in a reduction in velocity profile. By raising f, we can increase the temperature profile while reducing the concentration profile.
Increasing the value of f by 0.04 improves heat transfer by up to 17.02%.
Increasing the value of f by 0.04 reduces mass transfer by up to 3.18%.