Numerical Simulation of MHD Couette Flow of a Fuzzy Nanofluid through an Inclined Channel with Thermal Radiation Effect

Siddique, Imran; Zulqarnain, Rana Muhammad; Nadeem, Muhammad; Jarad, Fahd

doi:https://doi.org/10.1155/2021/6608684

Computational Intelligence and Neuroscience

On this page

Abstract Introduction Preliminaries Results and Discussion Conclusion Abbreviations Data Availability Conflicts of Interest References Copyright Related Articles

Special Issue

Artificial Intelligence and Machine Learning-Driven Decision-Making

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 6608684 | https://doi.org/10.1155/2021/6608684

Numerical Simulation of MHD Couette Flow of a Fuzzy Nanofluid through an Inclined Channel with Thermal Radiation Effect

Imran Siddique,¹Rana Muhammad Zulqarnain,²Muhammad Nadeem,¹and Fahd Jarad^3,4

Academic Editor: Ahmed Mostafa Khalil

Received14 Sept 2021

Revised06 Oct 2021

Accepted09 Oct 2021

Published31 Oct 2021

Abstract

The present study especially concerns the investigation of the Couette flow and heat transfer with thermal radiation through an inclined channel. Single-wall carbon nanotube (SWCNT) and multiple-wall carbon nanotube (MWCNT) are nanoparticles embedded in the host fluid. The dimensionless highly nonlinear differential equations (DEs) are solved via numerical scheme bvp4c. The effects of the physical parameters on heat transfer are presented in the form of graphs. The results demonstrate that the heat transfer is enhanced by using solid particle frictions (SWCNT and MWCNT). The large estimation of a magnetic parameter declines the velocity component. The current and existing results with their comparisons are shown in the tabular form for the validation of our code. The current results are in good agreement with their existing results. Generally, fuzziness or uncertainty is inherent in modeling, analysis, and experimentation. Due to the uncertain environmental conditions, fuzziness broadly exists in various engineering heat transfer problems. In this work, the nanoparticles’ volume fraction of the SWCNT and MWCNT is taken as uncertain parameters in terms of triangular fuzzy numbers (TFNs). The TFNs are controlled by the which has less computational effort for analyzing the fuzziness or uncertainties. Also, a comparison between the SWCNT and MWCNT through the membership function and the variability of the uncertainty is studied.

1. Introduction

The study of the flow of third-grade fluids [1–3] over an inclined channel is an important application in engineering, science, and technology. Some of these applications can be found in materials manufactured by the extraction procedure particularly in polymer processing, the flow of synovial fluid in human joints, geological flows inside the Earth’s mantle, microfluids, drilling of oil and gas wells, etc. In fluid dynamics, the study of fundamental flow, namely, Couette flows to attract the researchers by several non-Newtonian fluids because of their uses in the technology and engineering industry. The unidirectional flow is used in polymer engineering such as die flow, injection molding, extrusion, plastic forming, continuous casting, and asthenosphere flows [4–7]. Magnetohydrodynamics (MHD) deals with the study of the motion of electrical fluids in the presence of a magnetic field. MHD flow has significant importance applications in the inclined channel such as geophysical, astrophysical, metallurgical processing, MHD generators, pumps, geothermal reservoirs, polymer technology, and mineral industries. MHD fluid uses as a lubricant to stop the unexpected variation of fluid viscosity with temperature under certain norms. Thermal radiation is a process by which energy, in the form of electromagnetic radiation, is emitted by a heated surface in all directions and travels directly to its point of absorption at the speed of light; thermal radiation does not require an intervening medium to carry it. Thermal radiation plays an important role in science and technology, such as the furnace design, glass production, ship compressors, space vehicles, plasma physics, propulsion systems, spacecraft, flow structure of atomic plants, internal combustion engines, combustion processes, gas turbines, solar radiations, and solar power technology. Kamran and Siddique [8] calculated the analytical solutions for MHD flow between infinite parallel plates. Siddiqui et al. [9] calculated the deliberated flow of a third-grade fluid between two inclined parallel plates with heat transfer using the homotopy perturbation method (HPM). Aiyesimi et al. [10] studied the MHD Couette flow and Poiseuille flow problems using the homotopy perturbation method (HPM) and perturbation method (PM). They examined that both methods have the same solutions. Later on, Aiyesimi et al. [11] studied the solution of MHD Couette flow and Poiseuille flow problems of the temperature and velocity profile utilizing the perturbation method (PM). Yusuf et al. [12] studied the Couette flow with the effect of Soret, Dufour, and thermal radiation using the modified domain decomposition method (MADM) in an inclined channel. Okedayo et al. [13] used the MADM to study the viscous dissipation, Joule heating, and MHD flow of a third-grade fluid in an inclined channel. Farooq et al. [14] used the PM to study the Poiseuille flow of a couple of stress fluids with variable viscosity between two parallel inclined heated plates. Zeeshan et al. [15] studied the impact of variable thermal conductivity, heat generation, heat flux, and inclined uniform magnetic field on Poiseuille flow employing the homotopy analysis method (HAM). In this way, there are large kinds of literature about MHD and thermal radiation such as [16–23].

Heat transfer is important in industrial areas to launch transportation of energy in the system. Heat transfer is an essential task due to the requirement for energy in the world. In recent years, nanofluids have great attention because of their large-scale usage in industries. There are many applications of nanofluids in science and engineering, such as heat exchangers, solar cells, cooling of electronic equipment, solar water heaters, cooling of diesel-electric generators, nuclear reactors, cooling, and heating in buildings. The range of nanoscale particles is 1–100 nm, which benefits to increase the thermophysical properties of the nanofluid. Due to an increase in energy prices, the management of heat transfer is a vital role in energy systems. The nanofluids are the mixture of solid-liquid which contain nanoparticles for enhancing the heat transfer, investigated by Choi and Eastman [24]. The performance of nanoparticles in a heat transfer mechanism is better than the base fluid; the reason is that the suspended ultrafine particles increase the thermal conductivity of the mixture and raise their competency of energy transfer. Numerous literature studies [25–27] disclose the low volume fractions (1%–20% volume) for better performance of thermal conductivity of the suspensions; we can take more than 20% nanoparticles’ concentrations. Carbon nanotubes (CNTs) are cylindrical shape materials made of graphene. There are two types of CNTs, SWCNT and MWCNT, which are used in this work. As compared to other nanomaterials, SWCNT and MWCNT are of great importance for scholars due to their significant thermal conductivities, bent without any damage, and mechanical power. Saqib et al. [28] studied natural convection flow using carboxymethyl cellulose (CMS) as a base fluid and CNTs as hybrid nanofluids between two vertical parallel plates. For the exact solution, the Caputo–Fabrizio fractional derivative with the Laplace transform method is used. Hatami and Ganji [29] discussed the natural convection flow of sodium alginate (SA) as a base fluid and copper (Cu) and silver (Ag) as a nanofluid between two vertical parallel plates using the differential transform method (DTM). Khan et al. [30] studied the Navier slip boundary conditions for heat transfer using CNTs. Noreen et al. [31] studied the velocity slips and thermal analysis on MHD peristaltic flow in an asymmetric channel using CNTs. Ebaid and Al Sharif [32] deliberated the effect of MHD on nanofluids’ motion using CNTs and enhanced the heat transfer rate. The flow and heat transport of a special kind of second-grade hybrid nanofluid over a permeable stretching/shrinking sheet was addressed by Roy and Pop [33].

The flow of fluids with heat transfer is essential in science and engineering. After controlling the physical quantities like chemical diffusion, magnetic effect, volume friction of nanoparticles, and heat transmission, the fluid models are transformed into linear and nonlinear DEs. After governing these physical problems, they are converted into linear or nonlinear DEs. In general, the physical problems with involved geometry, coefficients, parameters, and initial and boundary conditions greatly affect the solution of DEs. Then, the coefficients, parameters, and initial and boundary conditions are not crisp due to the mechanical defect, experimental error, and measurement error. So, in this situation, fuzzy set theory is a powerful tool for a better understanding of the considered phenomena, and it is more accurate than assuming the crisp or classical physical problems. To be more specific, FDEs play a significant role in reducing uncertainity and providing a correct approach to characterise a physical issue with unknown parameters, initial and boundary conditions.

In 1965, Zadeh [34] presented the fuzzy set theory (FST). FST is a very valuable tool to define the situation in which information is imprecise, vague, or uncertain. FST is completely defined by its membership function or belongingness, and the membership function describes each element of the universe of discourse by a number from the [0, 1] interval. On the contrary, the degree of nonbelongingness is a complement to “one” of the membership degree or belongingness. Fuzzy number (FN) can be expected as a function whose range is specified from zero to one. Every numerical value in the range is allocated a definite grade of the membership function where zero signifies the minimum possible grade and one is the maximum possible grade. Arithmetic operations on FNs were developed by Dubois and Prade [35]. Different types of FNs can be categories in triangular, trapezoidal, and Gaussian fuzzy numbers. Here, we consider TFNs for the sake of completeness. Also, many scholars have applied FST to obtain well-known results in the field of commerce and science [36–46].

The information of dynamical systems modeled by partial or ordinary differential equations is commonly incomplete, vague, or uncertain, while FDEs represent a proper way to model the dynamical systems under vagueness or uncertainty. This impreciseness or vagueness can be defined mathematically using FNs or TFNs. In recent years, there have been many studies revolving around the concept of FDEs. Seikala [47] introduced the fuzzy differentiability concept. Later on, Kaleva [48] presented fuzzy differentiation and integration. Kandel and Byatt [49] introduced the FDEs in 1987. Buckley et al. [50] used two methods’ extension principle and FNs for the solution of FDEs. Gasilov et al. [51] studied the system of FDEs with the TFNs. Salahshour et al. [52] studied the fuzzy logistic equation and alley effect using FDEs with the help of TFNs. Biswal et al. [53] studied the natural convection of nanofluid flow between two parallel plates using the HPM in a fuzzy environment. The volume fraction of nanoparticles is considered as the TFN and also shows that the fuzzy result is better than a crisp result. Borah et al. [54] discussed the MHD flow of second-grade fluids in a fuzzy environment using fractional derivatives Atangana–Baleanu and Caputo–Fabrizio. The nondimensional governing equations are converted into fuzzified governing equations with the help of the Zadeh extension principle and triangular fuzzy number. MHD and ohmic heating on the third-grade fluid in an inclined channel in a fuzzy environment was investigated by Nadeem et al. [55]. To discuss the uncertainty, the triangle membership function was used.

In the literature review, it scrutinized that there is no such study that has been accounted in an MHD Couette flow with pure water as the base fluid and CNTs (SWCNTs and MWCNTs) as nanoparticles through the inclined heated channel in the presence of thermal radiation. The governing nonlinear coupled ordinary differential equations are solved by numerical scheme bvp4c. Furthermore, after checking the accuracy of bvp4c, the results of existing works in the literature are compared. The nanoparticles’ volume fraction performs the main role in the enhancement of heat transfer and thermophysical properties. The particular applications of nanoparticles are in solar water heating, engine cooling, cooling of transformer oil, cooling of the radiator, cooling in machining, and defense. Some researchers take nanoparticles’ volume fraction 0 to 4% for the enhancement of heat transfer. However, from different sources [25–27], we can take nanoparticles’ volume fraction more than 20% for better performance of heat transfer. So, the volume fraction of nanoparticles (SWCNTs and MWCNTs) is taken as a fuzzy number or uncertain because it depends on the added nanoparticle volume or dimensions of particles or the shape of particles of the fluid. In this respect, volume fraction has been taken as an uncertain parameter in terms of fuzzy numbers in the present study. To handle this problem in a fuzzy environment, FDEs and the triangular fuzzy plot with the approach were used. Also, the comparison of SWCNTs and MWCNTs is discussed through triangular fuzzy plots.

2. Preliminaries

In this section, some basic notations and definitions are used.

Definition 1 (see [34]). Fuzzy set is defined as a set of ordered pairs such that , where is the universal set and is the membership function of and mapping defined as .

Definition 2 (see [35]). The or of a fuzzy set is a crisp set and defined by , where .

Definition 3 (see [35, 55]). Let with membership function be called a membership function of the TFN ifMembership functions allow us to graphically represent a fuzzy set. The x-axis signifies the universe of discourse, while the y-axis signifies the degrees of membership in the [0, 1] interval. The TFNs with peak (or center) a₂, left width a₂- a₁ > 0, and right width a₃-a₂ > 0 are transformed into the interval numbers through α-cut approach, which is written as , where . A TFN and the -cut of the membership function are seen in Figure 1. An arbitrary TFN satisfies the following conditions: (i) is an increasing function on [0, 1]. (ii) is a decreasing function on [0, 1]. (iii) on [0, 1]. (iv) and are bounded on left continuous and right continuous at [0, 1], respectively. (v) where is a crisp number at .
must be a closed interval for every ; also, is called the level of credibility or presumption. Membership function or grade is also named as a grade of possibility or grade of credibility for a given number. Also, Figure 1 describes the membership function of a triangular fuzzy uncertainty. So, the triangular fuzzy uncertainty is defined as (lower bound), (most belief value), and (upper bound).

Definition 4 (see [47, 50, 55]). Let I be a real interval. A mapping is called a fuzzy process, defined as and . The derivative of a fuzzy process is defined by .

Definition 5 (see [47, 50, 55]). Let be a fuzzy-valued function defined on I. Let for all -cuts. Assume that and have continuous derivatives or are differentiable, for all and ; then, Similarly, we can define higher-order ordinary derivatives in the same way. An FN by an ordered pair of functions satisfies the following conditions: (i) and are continuous on [0, 1]. (ii) is an increasing function on [0, 1]. (iii) is a decreasing function on [0, 1]. (iv) on [0, 1].

3. Problem Formulation

Consider the steady incompressible Couette flow of pure water as a base fluid, and SWCNT and MWCNT nanoparticles are added separately in an inclined channel with heat transfer. The distance between two inclined parallel plates is 2H with the x-axis as the flow direction. The upper plate is moving with constant speed U, while the lower plate is fixed. The plates having constant temperatures with the upper plate instantaneously change in temperature . An unchanging magnetic field is applied in the y-direction and is expected undisturbed as the induced magnetic field is neglected under the assumption of a small magnetic Reynolds number. Radiative heat flux is considered here. The pressure and ambient air are ignored so that the flow is due to the movement of the upper plate and alone with gravity. The properties of SWCNTs, MWCNTs, and water are given in Table 1.

In the equations, the flow of an incompressible, unidirectional, third-grade fluid with the effects of MHD is given [8–11]:where is the constant density, is the material derivative, is the velocity vector, is the pressure, is the stress tensor, is the electric current density, is the total magnetic field, denotes the imposed magnetic field, and represents the induced magnetic field.where is the electric field (), is the magnetic permeability, and is the electrical conductivity.wherein is the coefficient of viscosity and , and are material constants. The Rivlin–Ericksen tensors are defined as is the identity tensor:

The geometry of the problem is presented in Figure 2.

For radiation, Rosseland approximation is used [16–18]:where is the Stefan–Boltzmann constant, is the absorption coefficient, and is the linear temperature function. is expanded by Taylor series expansion about as

Define velocity and temperature profiles for one-dimensional flows as follows [8, 11, 55]:

Using equations (2)–(6) along with equations (7)–(11),

Introduce the following nondimensional parameters. Eqs. (11) and (12), with the boundary conditions become, work done due to deformation and Joule heating in a nondimensional form is given aswith boundary conditionswith boundary conditionswhere is a third-grade fluid parameter, is a gravitational parameter, is a magnetic parameter, is the Brinkman number, and is the thermal radiation.

Here, denote the density, viscosity, thermal conductivity, specific heat, electrical conductivity, and nanoparticles’ volume fraction of nanofluids, respectively [29, 32].

The properties of the base fluid and nanoparticles are listed in Table 1.

3.1. Formulation of the Crisp Problem into the Fuzzy Problem Using FDEs

The velocity and temperature are affected by a small change in the value of the volume fraction of nanoparticles. Some researchers take the volume fraction of nanoparticles in the range [0.01–0.04], so the point is that the flow of fluid just depends on these values. Then, uncertainty arises due to the fixed crisp values of the volume fraction of nanoparticles. So, it is better to handle a difficult problem in a fuzzy environment by taking volume fraction as the FN.

For the fuzzy form, the governing coupled differential equations (15)–(18) can be written aswith boundary conditions

For the fuzzy solution, equations (15)–(18) can be converted into FDEs using the approach. So, according to Definitions 4 and 5, we havewith boundary conditionswhere “^___” stands for the fuzzy form and fuzzy velocity profile is . Here, is the lower bound and the upper bound of fuzzy velocity profiles. Similarly, the fuzzy temperature profiles are .

The crisp values and TFNs of these FNs are listed in Table 2. The TFN defined the variation of the FN at each . The TFNs are used to describe the triangular membership functions of the FNs which range from 0 to 1; see Figure 1. The investigated ranges are generally used to build up the said problem.

Now we proposed a numerical method bvp4c to solve crisp differential equations and FDEs. It is a Lobatto IIIa formula with three stages based on the finite-difference algorithm. It has a collocation polynomial, and in [a, b], the collocation formula yields a sixth-order accurate uniform C1 continuous solution. For error control and mesh selection, the continuous solution residual is employed. Further solvers, such as NDSOLVE, HPM, HAM, and ADM, and many other similar approaches are less consistent than bvp4c.

We transformed the governing ODEs to the system of the first order as follows.

Let

Boundary conditions are

To find the solution of Eqs. (24) to (28), a code in MATLAB software for bvp4c method is constructed.

4. Results and Discussion

In this section, pure water is chosen as a base fluid, and SWCNT and MWCNT are nanoparticles added into the base fluid to enhance the heat transfer rate between two heated inclined parallel plates. The numerical solutions of governing coupled nonlinear differential equations are achieved through numerical technique bvp4c. The effect of crisp thermophysical parameters, such as nanoparticle volume fraction , gravitational parameter K, magnetic parameter M, Brinkman number Br, and nondimensionless non-Newtonian viscosity , on velocity and temperature profiles is drawn in Figures 3–8.

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)

Tables 3 and 4 show the comparison of velocity and temperature fields for different values of [10, 11] when . The validation of our results of the present study was found to be in excellent agreement.

Figure 3 demonstrates the impact of on velocity and temperature profiles when other parameters are fixed. It is seen that the velocity and temperature of the nanofluid increase with increasing the values of because of the increase in the boundary layer thickness. Figure 4 displays the influence of gravitational parameter K on velocity and temperature profiles when other physical parameters are fixed. It is noticed that the velocity and temperature of the nanofluid rise rapidly in the center of the inclined plates with increasing the values of K. The reason is that when K increases, the upper plate movement expands the nanofluid velocity, which strengthens Joule dissipation; thereby, the rate of heat transfer is enhanced. The effect of M on the velocity and temperature profile is seen in Figure 5. When the value of M is increased, the velocity drops due to the Lorentz forces, but the temperature profile rises. The rate at which the nanofluid flow declines is less noticeable from the lower plate to the higher plate. Physically, when the magnetic parameter rises, the rate of nanofluid flow near the moving upper plate is significantly reduced compared to the stationary lower plate. Also, it is noticed that increases in the magnetic parameter cause an increase in Joule heating, which raises the heat transfer rate. Figure 6 visualizes the effect of N on the nanofluid of the temperature profile. As depicted, it is seen that the temperature profile performs a decreasing function with increasing N, therefore, representing a damping impact on the heat transfer performance and nanofluid flow. The effect of the viscous dissipation parameter on the temperature profile is demonstrated in Figure 7. It is observed that the temperature of the nanofluid increases with increasing the values of . The physical meaning of is when increases, the dissipation of heat in the boundary layer region increases, so the flow of the nanofluid is enhanced, and the heat transfer rate grows.

The effect of via the velocity profile of nanomaterials SWCNT and MWCNT is shown in Figure 8(a). When increases, the velocity of the nanofluid decreases. The velocity of the regular fluid is maximum at , which means that, by increasing the nanofluid volume fraction, it becomes denser. Physically, it means that the boundary layer of nanofluids is thicker than regular fluids; consequently, the velocity shows a decreasing behavior with increasing values of . The effect of on the temperature profile of nanomaterials SWCNT and MWCNT is depicted in Figure 8(b). It can be seen that the temperature profiles increase with an increase in . Physically, this is true because of the thicker thermal boundary layer which increases the heat transfer.

Now, we discuss the nanoparticles’ volume fraction of the SWCNT and MWCNT in a fuzzy environment. The governing equations of momentum and energy are converted into FDEs; then, the bvp4c scheme is employed for the numerical solution. The nanoparticle volume fraction ψ is taken as TFN (see Table 2). Moreover, we utilized the α-cut approach (0≤α≤1) to analyze the uncertainty. Then, the velocity and temperature profiles are said to be fuzzy. The controls the fuzzy terms, for instance, if , it will cover the whole interval of nanoparticles’ volume fraction; that is, . α-cut increases as (0.05 to 0.95) the width of lower and upper bounds of fuzzy velocity or temperature profiles decreases. When α-cut = 1, then the lower and upper bounds of fuzzy velocity or temperature profiles are coherent with each other, so they provide a crisp result. It is important to note that if the width between lower and upper bounds of the velocity or temperature profile is less, the uncertainty is less. The fuzzy velocity and temperature profiles are plotted in Figures 9 and 10 for some particular values of , . The triangular membership functions are depicted in Figures 11 and 12 for different values of y.

(a)

(b)

(a)

(b)

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

Figures 9 and 10 show the nanoparticles’ volume fraction of SWCNTs and MWCNTs taken as the TFN (see in Table 1); then, the fuzzy velocity and temperature profiles are controlled by . Figures 9(a) and 10(a) describe the effect of the on fuzzy velocity profile for the volume fraction of SWCNTs and MWCNTs which are TFNs, respectively. It is seen that significantly affects the fuzzy velocity profile of the regular fluid and nanofluid. The lower bound of the velocity profile shows the regular fluid, and the upper bound of the velocity profile displays the nanofluid at . Physically, it means that, at , the thermal boundary layer of is thinner, and the thermal boundary layer of is denser. The velocity of a regular fluid is maximum compared to the velocity of the nanofluid for different values of . When increases, the width between and decreases, and at , they coherent with one another which is a real flow of the nanofluid. It is noted that the width between and is very less, so the uncertainty is less.

The effect of the on temperature profile is represented in Figures 9(b) and 10(b) for the fuzzy volume fraction of SWCNTs and MWCNTs, respectively. The lower bound of the temperature profile shows the regular fluid, and the upper bound of the temperature profile shows the nanofluid at . Physically, it means that, at , the thermal boundary layer of is denser, and the thermal boundary layer of is thinner. Momentum and thermal boundary layers show opposite behavior at , which shows that the heat transfer is maximum. When increases, the width between and decreases, and at , they coherent with one another which is the real flow of the nanofluid. It is noted that the width between and is very less, so the uncertainty is less.

The comparison of SWCNTs and MWCNTs through fuzzy plots for various values of y is depicted in Figures 11 and 12. The volume fraction of SWCNTs and MWCNTs is the TFN. Figures 11 and 12 represent the fuzzy plots of fuzzy velocity and temperature profiles for different values of y. The fuzzy velocity of SWCNTs is greater than the fuzzy velocity of MWCNTs. Physically, it is correct because the density of SWCNTs is less than the density of MWCNTs. Comparing for the fuzzy temperature distribution, SWCNTs have greater thermal conductivity than MWCNTs but less density. So, SWCNTs bear more heat than MWCNTs, and in this study, we suggest that SWCNTs are better for the enhanced heat transfer as compared to MWCNTs. Also, the SWCNTs show better behavior as compared to MWCNTs because of their less width according to the membership function.

5. Conclusion

In this study, the effect of CNTs on MHD Couette flow via an inclined channel in a fuzzy environment is reported. The SWCNTs and MCNTs are nanoparticles, while pure water is the base fluid. The effects of Brinkman number (), magnetic parameter (M), the volume fraction of nanoparticles viscosity (), thermal radiation (N), and gravitational parameter (K) are considered. The impact of different parameters on heat transfer, velocity, and temperature profiles is analyzed. For mathematical computation, bvp4c is used, and the present results are found to be in excellent agreement as compared to existing results. The governing nonlinear DEs are converted into FDEs; then, the numerical technique bvp4c employed. The volume fraction of nanoparticles is considered as TFNs through , and the fuzziness is controlled. Some of the important and convenient achieved results are as follows:(i)The velocity and temperature profiles increase when the gravitation parameter (K) and fluid parameter increase.(ii)The velocity profile decreases when a magnetic parameter (M) increases, and the temperature profile increases near the boundary layer.(iii)The temperature profile decreases when the thermal radiation parameter (N) increases due to heat transfer. The temperature profile increases with an increase in the Brinkman number (iv)The temperature profile increases and the velocity profile decreases when the volume fractions of SWCNTs and MWCNTs increase.(v)The comparison of SWCNTs and MWCNTs shows that SWCNTs have a better heat transfer rate as compared to MWCNTs. Also, the comparison of SWCNTs and MWCNTs is examined through triangular membership plots of fuzzy velocity and temperature profiles, which shows that the width of MWCNTs is greater than that of SWCNTs.(vi)The results indicate that the crisp solution is always in-between the upper and lower solutions when increases from 0 to 1. The lower solution shows the minimum flow of the fluid (regular fluid), and the upper solution shows the maximum flow of the nanofluid when .(vii)In future works, one can use other fuzzy numbers to solve the heat transfer problems.

Abbreviations

:	Cartesian coordinates
u:	Normal component of the flow
T:	Transpose of the matrix
θ (y):	Dimensionless temperature
:	Reference and ambient temperatures
M:	Magnetic parameter
:	Brinkman number
:	Third-grade fluid parameter
P:	Pressure
:	Cauchy stress tensor
:	Material constants
:	Magnetic permeability
:	Coefficient of viscosity
:	Fuzzy temperature profile
:	Stefan–Boltzmann constant
:	Absorption coefficient
N:	Thermal radiation
:	Density of the nanofluid
:	Viscosity of the nanofluid
:	Specific heat of the nanofluid
:	Nanoparticles’ volume fraction
:	Lower fuzzy temperature profile
:	Upper fuzzy temperature profile
d/dt:	Material derivative
:	Electrical conductivity
K:	Gravitational parameter
:	Density of the fluid
:	Level or cut technique
FDE:	Fuzzy differential equation
:	Membership function
:	Fuzzy velocity profile
:	Lower fuzzy velocity profile
:	Upper fuzzy velocity profile
TFN:	Triangular fuzzy number
B:	Total magnetic field
f:	External body force
:	Rivlin–Ericksen tensors
:	Identity tensor
H:	Distance between two plates
:	Thermal conductivity of the nanofluid
:	Electrical conductivity of the nanofluid.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

C. Truesdell and W. Noll, The Non-Linear Field’s Theories of Mechanics, Springer, Berlin, Germany, 3rd edition, 2004.
K. R. Rajagopal, “On the stability of third-grade fluids,” Archives of Mechanics, vol. 32, no. 6, pp. 867–875, 1980.
View at: Google Scholar
K. R. Rajagopal, “Thermodynamics and stability of fluids of third grade,” vol. 339, pp. 351–377, 1980.
View at: Google Scholar
S. I. Natarov and C. P. Conrad, “The role of poiseuille flow in creating depth-variation of asthenospheric shear,” Geophysical Journal International, vol. 190, no. 3, pp. 1297–1310, 2012.
View at: Publisher Site | Google Scholar
T. Hayat, M. Khan, and M. Ayub, “Couette and poiseuille flows of an oldroyd 6-constant fluid with magnetic field,” Journal of Mathematical Analysis and Applications, vol. 298, no. 1, pp. 225–244, 2004.
View at: Publisher Site | Google Scholar
T. Hayat, R. Naz, and M. Sajid, “On the homotopy solution for poiseuille flow of a fourth grade fluid,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 3, pp. 581–589, 2010.
View at: Publisher Site | Google Scholar
T. Chinyoka and O. D. Makinde, “Analysis of transient generalized couette flow of a reactive variable viscosity third-grade liquid with asymmetric convective cooling,” Mathematical and Computer Modelling, vol. 54, no. 1, pp. 160–174, 2011.
View at: Publisher Site | Google Scholar
M. Kamran and I. Siddique, “MHD couette and poiseuille flow of a third grade fluid,” Open Journal of Mathematical Analysis, vol. 1, no. 2, pp. 1–19, 2017.
View at: Publisher Site | Google Scholar
A. M. Siddiqui, A. Zeb, Q. K. Ghori, and A. M. Benharbit, “Homotopy perturbation method for heat transfer flow of a third grade fluid between parallel plates,” Chaos, Solitons & Fractals, vol. 36, no. 1, pp. 182–192, 2008.
View at: Publisher Site | Google Scholar
Y. M. Aiyesimi, G. T. Okedayo, and O. W. Lawal, “Viscous dissipation effect on the MHD flow of a third grade fluid down an inclined plane with ohmic heating,” Mathematical Theory and Modeling, vol. 2, no. 9, pp. 2224–5804, 2012.
View at: Google Scholar
Y. M. Aiyesimi, G. T. Okedayo, and O. W. Lawal, “Effects of magnetic field on the MHD flow of a third grade fluid through inclined channel with ohmic heating,” Journal of Applied Mathematics and Computational Mechanics, vol. 3, no. 2, pp. 1–6, 2014.
View at: Google Scholar
A. Yusuf, Y. M. Aiyesimi, M. Jiya, and G. T. Okedayo, “Analysis of couette flow of a nanofluid in an inclined channel with soret and dufour effects,” American Journal of Computational and Applied Mathematics, vol. 6, no. 2, pp. 57–64, 2016.
View at: Google Scholar
G. T. Okedayo, O. E. Abolarin, and E. O. Titiloye, “A modified adomian decomposition method for a MHD third grade flow with ohmic heating,” Ilorin Journal of Science, vol. 2, no. 2, pp. 455–463, 2015.
View at: Google Scholar
M. Farooq, M. T. Rahim, S. Islam, and A. M. Siddiqui, “Steady poiseuille flow and heat transfer of couple stress fluids between two parallel inclined plates with variable viscosity,” Journal of the Association of Arab Universities for Basic and Applied Sciences, vol. 14, no. 1, pp. 9–18, 2013.
View at: Publisher Site | Google Scholar
A. Zeeshan, N. Shehzad, and R. Ellahi, “Analysis of activation energy in couette-poiseuille flow of nanofluid in the presence of chemical reaction and convective boundary conditions,” Results in Physics, vol. 8, pp. 502–512, 2018.
View at: Publisher Site | Google Scholar
M. S. Ullah, A. Tarammim, and M. J. Uddin, “A study of two dimensional unsteady MHD free convection flow over a vertical plate in the presence of radiation,” Open Journal of Fluid Dynamics, vol. 11, no. 1, pp. 20–33, 2021.
View at: Publisher Site | Google Scholar
S. Afzal, I. Siddique, F. Jarad, R. Ali, S. Abdal, and S. Hussain, “Significance of double diffusion for unsteady carreau micropolar nanofluid transportation across an extending sheet with thermo-radiation and uniform heat source,” Case Studies in Thermal Engineering, vol. 28, Article ID 101397, 2021.
View at: Google Scholar
M. Q. Brewster, Thermal Radiative Transfer and Properties, Wiley, Hoboken, NJ, USA, 1992.
T. Hayat, K. Hutter, S. Asghar, and A. M. Siddiqui, “MHD flows of an oldroyd-B fluid,” Mathematical and Computer Modelling, vol. 36, no. 9-10, pp. 987–995, 2002.
View at: Publisher Site | Google Scholar
S. Islam, “Homotopy perturbations analysis of couette and poiseuille flows of a third-grade fluid with magnetic field,” Science International, vol. 22, no. 3, 2010.
View at: Google Scholar
A. Khan, D. Khan, I. Khan, F. Ali, F. U. Karim, and M. Imran, “MHD flow of sodium alginate-based casson type nanofluid passing through a porous medium with Newtonian heating,” Scientific Reports “Pros R Scotland A”, vol. 8, no. 1, p. 8645, 2018.
View at: Publisher Site | Google Scholar
D. Khan, A. U. Rahman, G. Ali, P. Kumam, A. Kaewkhao, and I. Khan, “The effect of wall shear stress on two phase fluctuating flow of dusty fluids by using light hill technique,” Water, vol. 13, no. 11, p. 1587, 2021.
View at: Publisher Site | Google Scholar
A. Khan, D. Khan, I. Khan, F. Ali, F. U. Karim, and K. S. Nisar, “MHD flow of brinkman type H₂O-Cu, Ag, TiO₂ and Al₂O₃ nanofluids with chemical reaction and heat generation effects in a porous medium,” Journal of Magnetics, vol. 24, no. 2, pp. 262–270, 2019.
View at: Publisher Site | Google Scholar
S. U. S. Choi and J. A. Eastman, “Enhancing thermal conductivity of fluids with nanoparticles,” in Proceedings of the ASME International Mechanical Engineering Congress and Exposition, pp. 99–105, San Francisco, CA, USA, 1995.
View at: Google Scholar
S. Lee, S. U.-S. Choi, S. Li, and J. A. Eastman, “Measuring thermal conductivity of fluids containing oxide nanoparticles,” Journal of Heat Transfer, vol. 121, no. 2, pp. 280–289, 1999.
View at: Publisher Site | Google Scholar
H. Masuda, A. Ebata, K. Teramae, and N. Hishinuma, “Alteration of thermal conductivity and viscosity of liquid by dispersing ultra-fine particles. dispersion of Al₂O₃, SiO₂ and TiO₂ ultra-fine particles,” Netsu Bussei, vol. 7, no. 4, pp. 227–233, 1993.
View at: Publisher Site | Google Scholar
Y. Xuan and Q. Li, “Heat transfer enhancement of nanofluids,” International Journal Heat Fluid Flow, vol. 21, pp. 158–164, 2000.
View at: Publisher Site | Google Scholar
M. Saqib, I. Khan, and S. Shafie, “Natural convection channel flow of CMC-based CNTs nanofluid,” European Physics Journal Plus, vol. 133, no. 549, 2018.
View at: Publisher Site | Google Scholar
M. Hatami and D. D. Ganji, “Natural convection of sodium alginate (SA) non-Newtonian nanofluid flow between two vertical flat plates by analytical and numerical methods,” Case Studies in Thermal Engineering, vol. 2, pp. 14–22, 2014.
View at: Publisher Site | Google Scholar
W. A. Khan, Z. H. Khan, and M. Rahi, “Fluid flow and heat transfer of carbon nanotubes along a flat plate with navier slip boundary,” Applied Nanoscience, vol. 4, pp. 633–641, 2013.
View at: Google Scholar
S. A. Noreen, S. Nadeem, and Z. H. Khan, “Thermal and velocity slip effects on the MHD peristaltic flow with carbon nanotubes in an asymmetric channel: application of radiation therapy,” Applied Nanoscience, vol. 4, pp. 849–857, 2013.
View at: Google Scholar
A. Ebaid and M. A. Al Sharif, “Application of laplace transform for the exact effect of a magnetic field on heat transfer of carbon nanotubes-suspended nanofluids,” Zeitschrift für Naturforschung A, vol. 70, no. 6, pp. 471–475, 2015.
View at: Publisher Site | Google Scholar
N. C. Roy and I. Pop, “Flow and heat transfer of a second-grade hybrid nanofluid over a permeable stretching/shrinking sheet,” European Physics Journal Plus, vol. 135, no. 768, 2020.
View at: Publisher Site | Google Scholar
L. A. Zadeh, “Fuzzy sets,” Information Control, vol. 8, no. 3, pp. 338–353, 1965.
View at: Publisher Site | Google Scholar
D. Dubois and H. Prade, “Operations on fuzzy numbers,” International Journal of Systems Science, vol. 9, no. 6, pp. 613–626, 1978.
View at: Publisher Site | Google Scholar
S. P. Mondal and T. K. Roy, “First order linear homogeneous ordinary differential equation in fuzzy environment based on laplace transform,” Journal of Fuzzy Set Valued Analysis, vol. 2013, pp. 1–18, 2013.
View at: Publisher Site | Google Scholar
H. Zarei, A. V. Kamyad, and A. A. Heydari, “Fuzzy modeling and control of HIV infection,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 893474, 17 pages, 2012.
View at: Publisher Site | Google Scholar
T. Allahviranloo and S. Salahshour, “Applications of fuzzy laplace transform,” Soft Computing, vol. 17, no. 1, pp. 145–158, 2013.
View at: Google Scholar
A. El Allaoui, S. Melliani, and L. S. Chadli, “A mathematical fuzzy model to giving up smoking,” in Proceedings of the IEEE 6th International Conference on Optimization and Appication, (ICOA), pp. 1–6, Beni Mellal, Morocco, April 2020.
View at: Publisher Site | Google Scholar
H. C. Bhandari and K. Jha, “An analysis of microbial population of chemostat model in fuzzy environment,” The Nepali Mathematical Sciences Report, vol. 36, no. 1, 2019.
View at: Publisher Site | Google Scholar
A. Rajkumar and C. Jesuraj, “Mathematical model for dengue virus infected populations with fuzzy differential equations,” Communications in Computer and Information Science, vol. 955, pp. 206–217, 2018.
View at: Publisher Site | Google Scholar
R. M. Zulqarnain, X. L. Xin, I. Siddique, W. Asghar Khan, and M. A. Yousif, “TOPSIS method based on correlation coefficient under pythagorean fuzzy soft environment and its application towards green supply chain management,” Sustainability, vol. 13, no. 4, p. 1642, 2021.
View at: Publisher Site | Google Scholar
R. M. Zulqarnain, I. Siddique, R. Ali, F. Jarad, A. Samad, and T. Abdeljawad, “Neutrosophic hypersoft matrices with application to solve multiattributive decision-making problems,” Complexity, vol. 2021, Article ID 5589874, 17 pages, 2021.
View at: Publisher Site | Google Scholar
R. M. Zulqarnain, I. Siddique, R. Ali, D. Pamucar, D. Marinkovic, and D. Bozanic, “Robust aggregation operators for intuitionistic fuzzy hypersoft set with their application to solve MCDM problem,” Entropy, vol. 23, no. 6, p. 688, 2021.
View at: Publisher Site | Google Scholar
R. M. Zulqarnain, I. Saddique, F. Jarad, R. Ali, and T. Abdeljawad, “Development of TOPSIS technique under pythagorean fuzzy hypersoft environment based on correlation coefficient and its application towards the selection of antivirus mask in COVID-19 pandemic,” Complexity, vol. 2021, Article ID 6634991, 27 pages, 2021.
View at: Publisher Site | Google Scholar
I. Siddique, R. M. Zulqarnain, R. Ali, F. Jarad, and A. Iampan, “Multicriteria decision-making approach for aggregation operators of pythagorean fuzzy hypersoft sets,” Computational Intelligence and Neuroscience, vol. 2021, Article ID 2036506, 19 pages, 2021.
View at: Publisher Site | Google Scholar
S. Seikala, “On the fuzzy initial value problem,” Fuzzy Sets and Systems, vol. 24, no. 3, pp. 319–330, 1987.
View at: Google Scholar
O. Kaleva, “Fuzzy differential equations,” Fuzzy Sets and Systems, vol. 24, no. 3, pp. 301–317, 1987.
View at: Publisher Site | Google Scholar
A. Kandel and W. J. Byatt, “Fuzzy differential equations,” in Proceedings of the International Conference Cybernetics and Society, pp. 1213–1216, Tokyo, Japan, 1978.
View at: Google Scholar
J. J. Buckley, T. Feuring, and Y. Hayashi, “Linear systems of first order ordinary differential equations: fuzzy initial conditions,” Soft Computing—A Fusion of Foundations, Methodologies and Applications, vol. 6, no. 6, pp. 415–421, 2002.
View at: Publisher Site | Google Scholar
N. A. Gasilov, A. G. Fatullayev, and Ş. Emrah Amrahov, “Solution method for a non-homogeneous fuzzy linear system of differential equations,” Applied Soft Computing, vol. 70, pp. 225–237, 2018.
View at: Publisher Site | Google Scholar
S. Salahshour, A. Ahmadian, and A. Mahata, “The behavior of logistic equation with alley effect in fuzzy environment: FDE approach,” International Journal of Applied and Computational Mathematics, vol. 4, no. 62, 2018.
View at: Publisher Site | Google Scholar
U. Biswal, S. Chakraverty, and B. K. Ojha, “Natural convection of nanofluid flow between two vertical flat plates with imprecise parameter,” Coupled Systems Mechanics, vol. 9, no. 3, pp. 219–235, 2020.
View at: Google Scholar
G. Borah, P. Dutta, and G. C. Hazarika, “Numerical study on second-grade fluid flow problems using analysis of fractional derivatives under fuzzy environment,” Soft Computing Techniques and Applications. Advances in Intelligent Systems and Computing, vol. 1248, 2021.
View at: Publisher Site | Google Scholar
M. Nadeem, I. Siddique, F. Jarad, and R. N. Jamil, “Numerical study of MHD third-grade fluid flow through an inclined channel with ohmic heating under fuzzy environment,” Mathematical Problems in Engineering, vol. 2021, Article ID 9137479, 17 pages, 2021.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Imran Siddique et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

572

Downloads

839

Citations