Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Soret and Dufour effects on MHD squeezing flow of Jeffrey fluid in horizontal channel with thermal radiation

  • Nur Azlina Mat Noor,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Software, Validation, Writing – original draft, Writing – review & editing

    Affiliation Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, Johor Bahru, Johor, Malaysia

  • Sharidan Shafie,

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision

    Affiliation Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, Johor Bahru, Johor, Malaysia

  • Y. S. Hamed,

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision

    Affiliation Department of Mathematics and Statistics, College of Science, Taif University, Taif, Saudi Arabia

  • Mohd Ariff Admon

    Roles Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing

    ariffadmon@utm.my

    Affiliation Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, Johor Bahru, Johor, Malaysia

Abstract

The fluid flow with chemical reaction is one of well-known research areas in the field of computational fluid dynamic. It is potentially useful in the modelling of flow on a nuclear reactor. Motivated by the implementation of the flow in the industrial application, the aim of this study is to explore the time-dependent squeeze flow of magnetohydrodynamic Jeffrey fluid over permeable medium in the influences of Soret and Dufour, heat source/sink and chemical reaction. The presence of joule heating, joule dissipation and radiative heat transfer are analyzed. The flow is induced due to compress of two surfaces. Conversion of partial differential equations (PDEs) into ordinary differential equations (ODEs) is accomplished by imposing similarity variables. Then, the governing equations are resolved using Keller-box approach. The present outcomes are compared with previously outcomes in the literature to validate the precision of present outcomes. Both outcomes are shown in close agreement. The tabular and graphical results demonstrate that wall shear stress and velocity profile accelerate with the surfaces moving towards one another. Moreover, the concentration, temperature and velocity profiles decreasing for the increment of Hartmann numbers and Jeffrey fluid parameters. The impacts of heat generation/absorption, joule dissipation and Dufour numbers enhance the heat transfer rate and temperature profile. In contrast, the temperature profile drops and the heat transfer rate boosts when thermal radiation increases. The concentration profile decelerates, and the mass transfer rate elevates with raise in Soret number. Also, the mass transfer rate rises for destructive chemical reaction and contrary result is noted for convective chemical reaction.

1 Introduction

The movement of two parallel surfaces approaching one another is caused by external stress. The design of squeezing flow is adapted by researchers in the mechanical appliances such as modelling of oil flow in the bearings, lubrication system, injection moulding and hydraulic lift. The fundamental studies on the behaviour of flow in two surfaces was explored by Stefan [1]. Formulation model of squeeze flow with corresponded boundary conditions is derived via lubrication principle. Further, numerous works are done to analyse the behaviour of squeezing flow in the different geometries. The squeeze flow of elliptic and rectangle geometries was reported by Reynolds [2] and Archibald [3]. The previous works were formulated using Reynolds equation. However, Ishizawa [4] and Jackson [5] stated that Reynolds equation are not suitable in the analysis of squeeze flow in the high velocity and porous thrust bearings. Therefore, the fundamental mathematical model of squeeze flow has been revised and renewed in the various studies [611].

The research on boundary layer flow of non-Newtonian fluid attracted the interest of scientists due to the widespread in engineering application. Several models was proposed to discover the rheological behavior of non-Newtonian fluids. It is discovered that Jeffrey model is a simplest linear model with time derivatives as a substitute to convective derivatives [12]. It is categorized as shear thinning fluid because of high shear viscosity and yield stress [13]. The constitutive equation of Jeffrey fluid model was originally proposed by Pavlovskii [14] to investigate the dynamics of aqueous polymer solution. It is discovered that the addition of small amount of polymers in viscous fluid decrease the friction caused by drag force in the fluid flow [15]. Moreover, the model portrays the viscoelastic characteristics for the polymer industries by considering the relaxation and retardation parameter [16]. Common example of Jeffrey fluid is low-concentrated aqueous polymer solution [17]. Many authors used Jeffrey fluid model for studying the aqueous Polyacrylamide solution [18], blood flow in narrow arteries [19], movement of chyme in small intestine [20] and food bolus through esophagus [21] by applying the Jeffrey fluid model.

The studies of fluid flow over a porous medium has gain considerable attention due to the development of Darcy Law. A solid matrix with interconnected voids is known as permeable medium. It is explored in industrial and natural cases, for instances groundwater flow, petroleum reservoir rocks, engine coolant system and aircraft wings in the permeable cavities [22]. Moreover, the hydrodynamic of magnetic field in the electrical conducted fluid, magnetohydrodynamics (MHD) is widely reviewed because of its applications in MHD pump and generator. The hydromagnetic flow of Jeffrey fluid is discovered in many geometries. Hayat et al. [23] analysed the presence of injection or suction on MHD squeezing flow of Jeffrey fluid in the two permeable surfaces using homotopy analysis method (HAM) analytically. The squeeze flow of Jeffrey fluid with magnetic field on the stretching permeable lower plate with injection or suction was explored by Muhammad et al. [24]. Furthermore, Rao and Sreenadh [25] reviewed the magnetohydrodynamic flow of Jeffrey fluid on porous shrinking and stretching plate. The effect of hydromagnetic on Jeffrey fluid flow over porous medium in a circular tube was discovered by Nallapu and Radhakrishnamacharya [26]. Besides, the flow of Jeffrey fluid in the influence of magnetic field at a stagnation point was reviewed by Ahmad and Ishak [27]. The motion of the flow is due to stretched vertical plate. The problem was solved via Keller-box approach numerically.

The investigation on flow with viscous or joule dissipation is a subject of interest due to the several applications such as injection molding, high-rate extrusion and high temperature in polymer processes. The influence of joule dissipation is significant in the high velocity or viscosity fluid [28]. Hayat et al. [29] discussed the heat transfer of Jeffrey fluid on a stretched surface with joule dissipation. The convective heat transfer of Jeffrey fluid flow across a stretched plate with MHD, joule dissipation, joule heating, heat generation/absorption and radiative heat transfer was discovered by Ahmed et al. [30]. The impacts of joule dissipation and heating on the flow of Jeffrey fluid across a stretching surface was reviewed by Ahmad and Ishak [31]. Later, Zokri et al. [32] reported the numerical solution of Jeffrey fluid flow in a horizontal circular cylinder with the influence of viscous dissipation.

The electromagnetic waves emitted caused by the heat of substance is called thermal radiation. It is classified as one of the fundamental mechanisms of heat transfer. The sight of thermal radiation is common in the power generation, space vehicles, gas turbines and nuclear reactor chilling [33]. The effects of radiative heat transfer on Jeffrey fluid flow were reviewed for several geometries. Hayat et al. [34] explored the flow of Jeffrey fluid across a stretched sheet in the influence of thermal radiation at a stagnation point. The impacts of thermal radiation on squeezing flow of Jeffrey fluid in two disks was reported by Hayat et al. [35]. Furthermore, Hayat et al. [36] analysed the mixed convection flow of Jeffrey fluid past an inclined stretching sheet in the presence of thermal radiation. Kavita et al. [37] reported the oscillatory flow of MHD Jeffrey fluid on a vertical channel with radiative heat transfer.

The simultaneous thermal and mass transfer with chemical reaction has arises in many practical process including the flow in a desert and evaporation at the water surface. A chemical reaction among the fluid and the foreign particles takes place in the chemical industrial process. The common examples are the polymer production, the food processing and the manufacture of ceramics or glassware [38]. Alsaedi et al. [39] analysed the convective thermal transfer of Jeffrey fluid on a stretched surface with chemical reaction. The influences of radiative heat transfer and chemical reaction on oscillation flow of MHD Jeffrey fluid in horizontal channel was explored by Idowu et al. [40]. The flow and radiative heat transfer of Jeffrey fluid on a vertical porous surface with MHD and chemical reaction was discussed by Rao et al. [41]. Next, Saleem et al. [42] investigated the impacts of chemical reaction, heat source/sink and thermophoresis on convective thermal transfer of magneto-Jeffrey fluid on a rotated cone. The presence of chemical reaction on squeezing flow of Jeffrey nanofluid with magnetic field and velocity slip was examined by Noor et al. [43].

The phenomenon of thermal and mass transfer or double diffusion in a moving fluid plays a significant role in the field of petroleum reservoirs, nuclear waste disposal and air pollution [44]. It is noteworthy that the double diffusion process becomes more complicated due to the simultaneous occurrence of the driving potentials of heat and mass fluxes. The heat flux due to concentration gradient is indicated as diffusion-thermo or Dufour impact. Meanwhile, the mass flux due to temperature gradient is indicated as thermo-diffusion or Soret impact [45]. The Soret and Dufour term is found in non-dimensional concentration and energy equation, respectively. Generally, the Soret and Dufour impacts are not considered because the magnitude order is smaller than the impact specified by Fourier and Fick′s laws. Nevertheless, the impacts are considered when the presence of species at surface of fluid region have low density than the surrounding fluid [46]. Many researchers have analyzed the impacts of Soret and Dufour on flow with different geometries. Hayat et al. [47] explored the hydromagnetic flow of Casson fluid on a stretched surface with effect of Soret and Dufour. The mixed convection flow of nanofluid past a nonlinear stretching and shrinking surfaces in the influences of MHD, radiative heat transfer, Soret and Dufour was discussed by Pal et al. [48]. Further, Ullah et al. [49] examined the flow of Casson fluid across a nonlinear stretched plate with convective and slip boundaries. The presence of MHD, radiative heat transfer, joule dissipation, joule heating, heat generation or absorption, chemical reaction and Soret and Dufour was studied in the problem. The analysis of Soret and Dufour on squeeze flow of magneto-Casson fluid betweeen two surfaces with heat source/sink, joule heating and dissipation, chemical reaction and radiative heat transfer was investigated by Naduvinamani and Shankar [50].

The above cited papers reveal that the research focusing on squeeze flow of Jeffrey fluid over two surfaces are limited. Moreover, the thermal and mass transfer on squeezing flow of Jeffrey fluid is not yet covered. Thus, the aim of research is to explore Soret and Dufour impacts on time-dependent squeezing flow of MHD Jeffrey fluid embedded in a permeable medium with heat source/sink, joule heating and dissipation, chemical reaction and thermal radiation. The and later resolved through Keller-box technique. The numerical solutions are compared with published outputs in literature and shown in close agreement. Graphical outputs of concentration, temperature, and velocity profiles with correlated parameters are observed.

The fluid flow with chemical reaction is one of the most significant research areas in the field of computational fluid dynamic due to its industrial engineering applications. There are two categories of chemical reaction namely homogeneous and heterogeneous. The reaction is categorized as homogeneous or heterogeneous based on its occurrence in single phase (gaseous, liquid, or solid) or two phases (solid and gas, gas and liquid or solid and liquid), respectively [51]. Homogenous reaction occurs if the reactants and products are in the same phase while heterogeneous reactions have reactants in two or more phases [52].

The present study is mainly applied in the modelling of flow in a nuclear reactor. The presence of chemical reaction in the mathematical model is important to investigate the flow with nuclear reaction in the nuclear reactor. It is discovered that the lack of control of nuclear reaction may lead to the widespread contamination of air and water. Hence, the nuclear reaction flow is instantaneously stopped when the nuclear power plant accidents happen [53]. Moreover, the significance of MHD and permeable medium is analyzed in the fluid flow. Song et al. [54] discovered the capability of power conversion system raise in the presence of magnetohydrodynamics. The power conversion is essential to nuclear electric propulsion (NEP) system. NEP is simply electric propulsion in which the electricity is generated from a nuclear reactor. The safety standard enhance with the implementation of porous media concept as it accelerates the heat dissipation in the nuclear reactor [55]. The present study explores the following research questions:

  1. How do the mathematical models for unsteady MHD squeezing flow in a porous medium with heat source/sink, joule heating and dissipation, chemical reaction and thermal radiation can be formulated?
  2. How does the presence of Soret and Dufour will affect the temperature and concentration of Jeffrey fluid?
  3. What is the variation of wall shear stress, heat and mass transfer with increasing in magnetic field, porosity, viscosity of Jeffrey fluid, heat source, thermal radiation, chemical reaction, Soret and Dufour impacts?

2 Mathematical formulation

The time dependent MHD flow of Jeffrey fluid induced by squeezing of two surfaces through porous medium with heat source/sink, Soret and Dufour and chemical reaction is studied. Also, the influences of joule dissipation, joule heating and thermal radiation are considered. The distance of two surfaces is y = ±h(t) = ±l(1 − αt)1/2. The two surfaces are moving further when α < 0 and the surfaces are moving closer when α > 0 till t = 1/α with velocity . The lower plate is exerted with the magnetic field B(t) vertically [56]. Fig 1 depicts the geometrical model of Jeffrey fluid flow.

thumbnail
Fig 1. Schematic diagram of Jeffrey fluid between two squeezed plates embedded in a porous medium with transverse magnetic field.

https://doi.org/10.1371/journal.pone.0266494.g001

Referring to Jaluria [57], the constitutive relation of momentum equation for an incompressible fluid in two-dimensional form is given by (1) Here, represents the substantial derivative, ρ F denotes the body force and T is the Cauchy stress tensor. Based on Nadeem and Akbar [58], Cauchy stress tensor T of Jeffrey fluid is denoted by (2) with I is identity tensor and p is pressure. The definition of extra stress tensor S is given by (3) where A1 is the Rivilin-Ericksen tensor.

The derivation of energy equation is according to the first law of thermodynamics which stated energy is conserved. For any closed system, the energy remains constant and cannot be created nor destroyed over time. The energy equation in vector form is given by [59] (4) where q is the heat flux due to temperature gradient, jp.T is the heat flux caused by thermophoretic effect, Φ is the viscous dissipation, is Joule heating with the current density, J. The definition of q and jp.T are expressed as (5) Substitute Eq (5) into Eq (4) yields (6) The radiative heat transfer, qr based on Roseland approximation in two-dimensional form is denoted as [60] (7) The small difference of temperature in the fluid is indicated as T4. The term T4 is expanded as linear function of temperature by Taylor’s series for T and the higher order terms are ignored yields (8) Then, substitue Eqs (7) and (8) into energy Eq (6) yields (9) The definition of J according to Ohm’s law is [28] (10) with B is the total magnetic field and E is the electric field. It is assumed that E is ignored because there are no polarization and external applied electric field. Then, Eq (8) becomes (11) Substitute Eq (11) into the energy Eq (9) yields (12) The viscous dissipation term Φ is expressed as [61] (13) Therefore, the energy Eq (12) can be written as follows (14) or (15)

The governing equations of Jeffrey fluid are reduced to the following equations using boundary layer approximation (16) (17) (18) (19) The corresponded boundary conditions (BCs) are (20) (21) The non-dimensional variables are applied to simplify PDEs into ODEs [62]; (22) Substituting dimensionless variables (22) into Eqs (17), (18) and (19) gives the subsequent forms (23) (24) (25) with dimensionless BCs (26) (27) The pertinent terms in the dimensionless equations are described by (28) Physically, the motion of two surfaces is indicated by squeeze term, with S > 0 portrays the surfaces approaching nearer and S < 0 portrays the surfaces separating apart. Moreover, Darcy, Hartmann and Deborah terms are implemented to manage the velocity field. The temperature is regulated by radiative heat transfer, Eckert and heat generation or absorption terms. Further, the impacts of Dufour and Soret are explored in concentration and temperature graphs. The addition of chemical reaction term is examined in concentration profile.

3 Results and discussion

The governing Eqs (23) to (25) in corresponded BCs (26) and (27) are resolved using Keller-box procedure numerically. The four steps to get the numerical results are

  1. The conversion of ODEs to a system of 1st order equations.
  2. The discretization of 1st order equations into the form of finite difference via central difference method.
  3. The linearization of nonlinear equations with Newton’s method and addressed in the form of matrix-vector.
  4. The linear system is solved using block tri-diagonal elimination technique.

The algorithm based on Keller-box approach is developed in MATLAB for the iterative computation to obtain numerical and graphical outputs. The proper values of the step size Δη = 0.01 and thickness of boundary layer η = 1 are compulsory to get the accurate outputs. The convergence criterion is referred to the variation in the current and former outputs of concentration, temperature and velocity. Calculation is ended when the numerical outputs converging to 10−5 [63].

The numerical calculation is done to examine the impacts of S, Ha, λ1, De, Da, Pr, Ec, γ, Rd, R, Sc, Du and Sr on concentration, temperature and velocity profiles. The algorithm built in the MATLAB is validated by comparison with the present results of −f′′(1), −θ′(1) and −ϕ′(1) with previous existing results by Naduvinamani and Shankar [50] as limiting cases. Both outputs are shown in the good agreement as displayed in Table 1.

thumbnail
Table 1. Numerical values of −f′′(1), −θ′(1) and −ϕ′(1) for S when λ1 → ∞, Da → ∞, De = 10−10, Rd = Ha = γ = Du = Sr = 0, δ = 0.1 and R = Sc = Pr = Ec = 1.

https://doi.org/10.1371/journal.pone.0266494.t001

Limiting cases analysis is conducted by setting the quantity of interest to a specific extreme value (usually 0 or ∞). For instance, the present work is reduced to the pioneer work by Naduvinamani and Shankar [50] as shown in Table 1. The values of dimensionless parameter is setting as λ1 → ∞, Da → ∞, Rd = Ha = γ = Du = Sr = 0, De = 10−10, δ = 0.1 and Sc = Ec = R = Pr = 1 when input in the MATLAB program. All the results are discovered in excellent agreement. Hence, it is proven that the Keller-box scheme used to get the present outcomes is accurate and acceptable.

The impacts of S, λ1, Ha, Da and De on axial velocity is plotted in Figs 2 to 6. The motion of surfaces closer is described by S > 0 and the movement of surfaces further is described by S < 0. Fig 2 shows the effects of S on axial velocity. It is noticed that boundary range close on the below surface is 0 ≤ η < 0.5 and the boundary range near on the above surface is 0.5 ≤ η ≤ 1. It is discovered that the velocity decreasing as η < 0.5 and it elevates as η ≥ 0.5 with S > 0. In contrary, the velocity enhancing as η < 0.5 and it decelerates as η ≥ 0.5 with S < 0. The flow cross the confined channel accelerates as the surfaces approaching closer, which resulting in the velocity profile accelerates. In contrast, the deceleration of velocity is due to the flow confronts higher resistance in the wider channel. There is cross flow behaviour at the midpoint of channel. The squeezing parameter does not affect the fluid velocity at ηc = 0.5 (critical point). The similar behaviour of velocity profile with the variation of S is shown in the work by Naduvinamani and Shankar [50]. The influence of λ1 on axial velocity is displayed in Fig 3. The velocity slowing down as η ≤ 0.45 and it accelerates as η > 0.45 as λ1 increasing. It is noticed that the flow decelerates due to the high fluid viscosity and the intermolecular forces of fluid molecules enhance as λ1 increases. Fig 4 portrays the effect of Ha on axial velocity. The velocity reduces as η ≤ 0.45 and it enhances as η > 0.45 with Ha rises. Lorentz force is generated by imposing the magnetic field on the electrical conducted fluid. The opposition on the flow elevates with the presence of Lorentz force and consequently, decelerate the velocity in the boundary area. The behaviour of velocity profiles when varying λ1 and Ha parameters is the same with the results of Hayat et al. [23]. The impact of Da on axial velocity is explored in Fig 5. The velocity increases as η ≤ 0.45 and it declining as η > 0.45 with enhance in Da. The increment of Darcy parameter boosts the permeability of medium and thus, accelerates the flow through permeable medium at the centre of boundary area. Fig 6 shows the influences of De on axial velocity. The velocity enhances at the below surface and it decreasing at the above surface with rising in De. The ratio of retardation time and the observation time is denoted by Deborah term. Retardation time is the delayed reaction to the internal stress or delay of elasticity. The fluid has high viscosity because the retardation time is longer as De increases. It implies that the enhancement of intermolecular forces of fluid molecules, resulting in the flow nearer the upper plate slows down. The variation of De on velocity profile is the same as illustrated in Muhammad et al. [24].

Fig 7 portrays the impacts of Pr on temperature region. The flow temperature elevating when Pr rises. The raise in Pr indicates the larger specific heat capacity of fluid. The temperature increases because of the heat absorption elevates in the fluid flow. The variation of Ec on temperature region is presented in Fig 8. It is discovered that temperature field boosts with the increment of Ec. The viscous dissipation is represented by Ec. The heat generated by the internal friction of fluid particles rises as Ec increases. Hence, it has resulting the increment of the temperature field. Fig 9 illustrates the influences of Rd on temperature region. The transfer of thermal energy caused by the emission of electromagnetic waves from heated substance is described as thermal radiation. The reduction of temperature profile occurs because the heat transfer from flow area to the surfaces elevates for increasing Rd values. The effect of γ on temperature region is indicated in Fig 10. The heat absorption and heat generation case are characterized by γ < 0 and γ > 0, respectively. The temperature profile drops when γ < 0 and it rising when γ > 0. The heat generation increase the thermal energy of fluid, which cause the temperature profile increases. Meanwhile, a contrary behaviour is discovered in the heat absorption case. Fig 11 discovers the impact of Du on temperature region. The temperature profile enhances with raise in Du. This behavior is due to the reason that the kinematic viscosity in the vicinity of boundary flow decreases. The resistance to the fluid flow is measured by kinematic viscosity. The kinetic energy of fluid particles accelerates because the resistance in the flow slowing down. The influences of Sr on temperature region is depicted in Fig 12. It is discovered that the increment of Sr boost the temperature profile. Soret number is inversely proportional to the kinematic viscosity. It implies that the fluid confronts less resistance and the kinetic energy in the boundary region increases. The variation of temperature profiles for Pr, Ec, Rd, γ, Du and Sr is similar as shown in Naduvinamani and Shankar [50] work.

Fig 13 demonstrates the impacts of Du on concentration region. The concentration in the flow drops for increasing Du. Dufour effect is known as the energy flux generated by concentration differences. It has result in the increment of fluid temperature by lowering the concentration in the flow. The variation of Sr on concentration region is portrayed in Fig 14. The concentration profile declines as Sr rises. Soret effect or thermal diffusion is the mass flux generated by temperature differences. The mass transfer from the fluid area to the above plate accelerate due to Soret effect and thus, decreasing the fluid concentration. Fig 15 explores the variation of Sc on concentration region. The concentration reduces with raise in Sc. The mass diffusion slows down when Sc increases. The decrease in mass diffusivity from the surface to the fluid flow result in the concentration drops. The influences of R on concentration region is illustrated in Fig 16. The chemical reaction impacts are characterised by destructive (R > 0) and constructive (R < 0). It is explored that the concentration region increasing with (R < 0) and it decelerates with (R > 0). The increment of constructive chemical reaction enhances the reaction rate in the boundary layer. This phenomenon raises the fluid concentration. The same behavior of concentration profiles is observed with impacts of Du, Sr, Sc and R in the previous results by Naduvinamani and Shankar [50].

4 Physical quantities of fluid flow

Physically, the dimensionless parameters in the fluid are skin friction, Nusselt and Sherwood terms. The friction force at the surface boundary is described by skin friction. Moreover, Nusselt and Sherwood terms represent the rate of thermal and mass transfer in the fluid and surfaces. The terms of Cfx, Nux and Shx are denoted by [64] (29) The non-dimensional form of Cfx, Nux and Shx are (30)

The numerical outputs of skin friction, Nusselt and Sherwood parameters with variation of the dimensionless parameters are shown in Tables 2 to 4. The effects of S, λ1, De, Ha and Da on wall shear stress is displayed in Table 2. It is found that wall shear stress boosts with increment of S and Ha, in contrast it decreases for elevating De, λ1 and Da. The velocity rises when the plates approaching closer, which resulting in the friction force in the boundary region elevates. Besides, the Lorentz force accelerates the flow near the boundary, and consequently enhances the frictional force in the flow. Meanwhile, the raise in viscosity of Jeffrey fluid owing to the higher values of λ1 and De strengthen the intermolecular force of fluid particles. It has caused the fluid velocity slowing down. Next, the flow over porous medium encounter opposition from the drag force near the boundary. Thus, the drop of velocity profile result in the wall shear stress declines. Table 3 presents the impacts of Pr, Ec, Rd, γ and Du on heat transfer rate. It is noticed that Pr, Ec, Rd, γ and Du elevate the Nusselt number. The ratio of convection heat transfer and diffusion heat transfer is described by Nusselt number. The increment of Pr, Ec, Rd, γ and Du accelerate the kinetic energy, which cause the temperature and convective heat transfer enhances in the flow boundary. Meanwhile, the temperature reduces when Rd increases as shown in Fig 9. The presence of thermal radiation raises the heat transfer from the flow area to the boundary and hence, decreasing the temperature profile. This behaviour promotes the heat transfer rate on the boundary region. The variation of Sc, R and Sr on Sherwood number is exhibited in Table 4. The rate of mass transfer enhances for increasing Sc, R and Sr. The ratio of convection mass transfer to the diffusion mass transfer is denoted as Sherwood number. The drop of concentration profile in Figs 14 to 16 implies that the mass transfer by diffusion decelerates with increase in Sc, Sr and R. Thus, it indicates that this phenomenon boosts the Sherwood number and the convective mass transfer.

thumbnail
Table 2. Numerical results of −(1 + 1/λ1)f′′(1) for S, De, Da, λ1 and Ha when Du = δ = Rd = Ec = Sr = 0.1, γ = 0.01, S = 1 and Pr = R = Sc = 1.5.

https://doi.org/10.1371/journal.pone.0266494.t002

thumbnail
Table 3. Numerical results of for Rd, Pr, γ, Ec and Du when Sr = Ha = 0.1, De = 0.01, Da = S = δ = 1 and R = λ1 = Sc = 1.5.

https://doi.org/10.1371/journal.pone.0266494.t003

thumbnail
Table 4. Numerical outputs of ϕ′(1) for Sc, Sr and R as De = γ = 0.01, Rd = Ec = Du = Ha = 0.1, δ = S = Da = 1 and Pr = λ1 = 1.5.

https://doi.org/10.1371/journal.pone.0266494.t004

5 Conclusion

The present work examines the influences of Soret and Dufour on unsteady hydromagnetic flow of Jeffrey fluid through permeable medium with radiative heat transfer, chemical reaction and heat generation or absorption. The presence of joule dissipation and heating was studied. The flow is caused by squeeze within two surfaces. The conversion of PDEs to ODEs via similarity transformation is conducted. Keller-box technique is applied to resolve the governing equations. The impacts of S, λ1, De, Ha, Da, Pr, Rd, Ec, γ, Sc, R, Du and Sr on velocity, temperature and nanoparticles concentration are studied. The importance findings of the discussion are deduced as:

  1. The velocity accelerating when the surfaces approach nearer (S > 0) and it decreasing when the surfaces separate further (S < 0) at the above surface.
  2. The wall shear stress enhances for increasing in S and Ha, in contrast it reduces when De, Da and λ1 increases.
  3. The velocity, concentration, and temperature drop with increment of λ1 and Ha.
  4. The velocity slows down near the upper boundary for enhancing De and Da.
  5. The presence of Ec, γ and Du raise the heat transfer rate and temperature.
  6. The enhancement of Rd elevate the rate of heat transfer and it decreasing the temperature region.
  7. The rate of mass transfer rising and the concentration dropping when Sr enhances.
  8. The concentration boosts in the constructive chemical reaction (R < 0) and it declines in the destructive chemical reaction (R > 0).
  9. The rate of mass transfer decelerates when R < 0 and it increases when R > 0.

References

  1. 1. Stefan M. Experiments on apparent adhesion. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science. 1874; 47: 465–466.
  2. 2. Reynolds O. On the theory of lubrication and its application to Mr. Beauchamp tower’s experiments, including an experimental determination of the viscosity of olive oil. Philosophical Transactions of the Royal Society of London. 1886; 177: 157–234.
  3. 3. Archibald FR. Load capacity and time relations for squeeze films. Transactions of the ASME. 1956; 78: 231–245.
  4. 4. Ishizawa S. The unsteady flow between two parallel discs with arbitary varying gap width. Bulletin of the Japan Society of Mechanical Engineers. 1966; 9: 533–550.
  5. 5. Jackson JD. A study of squeezing flow. Applied Scientific Research. 1963; 11: 148–152.
  6. 6. Leider PJ, Bird RB. Squeezing Flow between Parallel Disks: Theoretical Analysis. Industrial & Engineering Chemistry Fundamentals. 1974; 13: 336–341.
  7. 7. Grimm RJ. Squeezing flows of Newtonian liquid films: An analysis including fluid inertia. Applied Scientific Research. 1976; 32: 149–166.
  8. 8. Wang CY. The squeezing of a fluid between two plates. Journal of Applied Mechanics. 1976; 43: 579–583.
  9. 9. Bujurke NM, Achar PK, Pai NP. Computer extended series for squeezing flow between plates. Fluid Dynamics Research. 1995; 16: 173–187.
  10. 10. Rashidi M, Shahmohamadi H, Dinarvand S. Analytic approximate solutions for unsteady two-dimensional and axisymmetric squeezing flows between parallel plates. Mathematical Problems in Engineering. 2008.
  11. 11. Khan U, Ahmed N, Khan SI, Zaidi ZA, Xiao-Jun Y, Mohyud-din ST. On unsteady two-dimensional and axisymmetric squeezing flow between parallel plates. Alexandria Engineering Journal. 2014; 53: 463–468.
  12. 12. Hayat T, Qayyum S, Imtiaz M, Alsaedi A. Three-dimensional Rotating Flow of Jeffrey Fluid for Cattaneo-Christov Heat Flux. AIP Advances. 2016; 6: 1–11.
  13. 13. Noor NAM, Shafie S, Admon MA. Unsteady MHD flow of Casson nanofluid with chemical reaction, thermal radiation and heat generation/absorption. MATEMATIKA. 2019; 35: 33–52.
  14. 14. Pavlovskii VA. On theoretical description of weak aqueous solutions of polymers. Doklady Akademii Nauk SSSR. 1971; 200: 809–812.
  15. 15. Frolovskaya OA, Pukhnachev VV. Analysis of the Models of Motion of Aqueous Solutions of Polymers on the Basis of Their Exact Solutions. Polymers. 2018; 10(6): 684. pmid:30966718
  16. 16. Sandeep N, Sulochana C. Momentum and heat transfer behaviour of Jeffrey, Maxwell and Oldroyd-B nanofluids past a stretching surface with non-uniform heat source/sink. Ain Shams Engineering Journal. 2018; 9(4): 517–524.
  17. 17. Baranovskii ES. Global solutions for a model of polymeric flows with wall slip. Mathematical Methods in the Applied Sciences. 2017; 40: 5035–5043.
  18. 18. Hamza S. Modelling the Effect of Concentration on Non-Newtonian Apparent Viscosity of an Aqueous Polyacrylamide Solution. Global Journal of Physics. 2016; 5: 505–517.
  19. 19. Akbar NS, Nadeem S, Ali M. Jeffrey fluid model for blood flow through a tapered artery with a stenosis. Journal of Mechanics in Medicine and Biology. 2011; 11: 529–545.
  20. 20. Akbar NS, Nadeem S, Lee C. Characteristics of Jeffrey fluid model for peristaltic flow of chyme in small intestine with magnetic field. Results in Physics. 2013; 3: 152–160.
  21. 21. Pandey SK, Tripathi D. Unsteady model of transportation of Jeffrey fluid by peristalsis. International Journal of Biomathematics. 2010; 3: 473–491.
  22. 22. Ahmad S, Farooq M, Anjum A, Mir NA. Squeezing flow of convectively heated fluid in porous medium with binary chemical reaction and activation energy. Advances in Mechanical Engineering. 2019; 11: 1–12.
  23. 23. Hayat T, Sajjad R, Asghar S. Series solution for MHD channel flow of a Jeffrey fluid. Communications in Nonlinear Science and Numerical Simulation. 2010; 15: 2400–2406.
  24. 24. Muhammad T, Hayat T, Alsaedi A, Qayyum A. Hydromagnetic unsteady squeezing flow of Jeffrey fluid between two parallel plates. Chinese Journal of Physics. 2017; 55: 1511–1522.
  25. 25. Rao ME, Sreenadh S. MHD boundary layer flow of Jeffrey fluid over a stretching/shrinking sheet through porous medium. Global Journal of Pure and Applied Mathematics. 2017; 13: 3985–4001.
  26. 26. Nallapu S, Radhakrishnamacharya G. Jeffrey fluid flow through porous medium in the presence of magnetic field in narrow tubes. International Journal of Engineering Mathematics. 2014.
  27. 27. Ahmad K, Ishak A. Magnetohydrodynamic flow and heat transfer of a Jeffrey fluid towards a stretching vertical surface. Thermal Science. 2015.
  28. 28. Ramesh K. Effects of viscous dissipation and Joule heating on the Couette and Poiseuille flows of a Jeffrey fluid with slip boundary conditions. Propulsion and Power Research. 2018; 7: 329–341.
  29. 29. Hayat T, Iqbal Z, Mustafa M, Alsaedi A. Unsteady flow and heat transfer of Jeffrey fluid over a stretching sheet. Thermal Science. 2014; 18: 1069–1078.
  30. 30. Ahmed J, Shahzad A, Khan M, Ali R. A note on convective heat transfer of an MHD Jeffrey fluid over a stretching sheet. AIP Advances. 2015; 5: 1–11.
  31. 31. Ahmad K, Ishak A. MHD flow and heat transfer of a Jeffrey fluid over a stretching sheet with viscous dissipation. Malaysian Journal of Mathematical Sciences. 2016; 10: 311–323.
  32. 32. Zokri SM, Arifin NS, Mohamed MKA, Kasim ARM, Mohammad NF, Salleh MZ. Influence of viscous dissipation on the flow and heat transfer of a Jeffrey fluid towards horizontal circular cylinder with free convection: A numerical study. Malaysian Journal of Fundamental and Applied Sciences. 2018; 14: 40–47.
  33. 33. Hashim Hamid A, Khan M, Khan U. Thermal radiation effects on Williamson fluid flow due to an expanding or contracting cylinder with nanomaterials: Dual solutions. Physics Letters A. 2018; 382: 1982–1991.
  34. 34. Hayat T, Shehzad SA, Qasim M, Obaidat S. Thermal radiation effects on the mixed convection stagnation-point flow in a Jeffrey fluid. Zeitschrift fur Naturforschung. 2011; 66: 606–614.
  35. 35. Hayat T, Qayyum A, Alsaadi F, Awais M, Dobaie AM. Thermal radiation effects in squeezing flow of a Jeffrey fluid. European Phys. J. Plus. 2013; 128: 1–7.
  36. 36. Hayat T, Ashraf MB, Alsulami HH. On mixed convection flow of Jeffrey fluid over an inclined stretching surface with thermal radiation. Heat Transfer Research. 2015; 46: 515–539.
  37. 37. Kavita K, Prasad KR., Kumari BA. Influence of heat transfer on MHD oscillatory flow of Jeffrey fluid in a channel. Advances in Applied Science Research. 2012; 3: 2312–2325.
  38. 38. Sulochana C, Aparna SR, Sandeep N. Heat and mass transfer of magnetohydrodynamic Casson fluid flow over a wedge with thermal radiation and chemical reaction. Heat transfer. 2021; 50(4): 3704–3721.
  39. 39. Alsaedi A, Iqbal Z, Mustafa M, Hayat T. Exact solutions for the MHD flow of a Jeffrey fluid with convective boundary conditions and chemical reaction. Zeitschrift fur Naturforschung. 2012; 67: 517–524.
  40. 40. Idowu AS, Joseph KM, Daniel S. Effect of heat and mass transfer on unsteady MHD oscillatory flow of Jeffrey fluid in a horizontal channel with chemical reaction. Journal of Mathematics. 2013; 8: 74–87.
  41. 41. Rao ME, Sreenadh S, Sumalatha B. Effects of thermal radiation and chemical reaction on an unsteady MHD flow of a Jeffrey fluid past a vertical porous plate with suction. Inter. Journal of Pharmacy & Technology. 2017; 9: 31059–31078.
  42. 42. Saleem S, Al-Qarni MM, Nadeem S, Sandeep N. Convective heat and mass transfer in magneto Jeffrey fluid flow on a rotating cone with heat source and chemical reaction. Communications in Theoretical Physics. 2018; 70: 534–540.
  43. 43. Noor NAM, Shafie S, Admon MA. Slip Effects on MHD Squeezing Flow of Jeffrey Nanofluid in Horizontal Channel with Chemical Reaction. Mathematics. 2021; 9(11): 1215.
  44. 44. Noor NAM, Shafie S, Admon MA. Heat and mass transfer on MHD squeezing flow of Jeffrey nanofluid in horizontal channel through permeable medium. PLoS ONE. 2021; 16(5): e0250402.
  45. 45. Shah NA, Animasaun IL, Chung JD, Wakif A, Alao FI, Raju CSK. Significance of nanoparticle’s radius, heat flux due to concentration gradient and mass flux due to temperature gradient: The case of water conveying copper nanoparticles. Scientific Reports. 2021; 11: 1882. pmid:33479309
  46. 46. Alam MK, Bibi K, Khan A, Noeiaghdam S. Dufour and Soret Effect on Viscous Fluid Flow between Squeezing Plates under the Influence of Variable Magnetic Field. Mathematics. 2021; 9(19): 2404.
  47. 47. Hayat T, Shehzad SA, Alsaedi A. Soret and Dufour effects on MHD flow of Casson fluid. Applied Mathematics and Mechanics. 2012; 33: 1301–1312.
  48. 48. Pal D, Mandal G, Vajravalu K. Soret and Dufour effects on MHD convective–radiative heat and mass transfer of nanofluids over a vertical non-linear stretching/shrinking sheet. Applied Mathematics and Computation. 2016; 287–288.
  49. 49. Ullah I, Khan I, Shafie S. Soret and Dufour effects on unsteady mixed convection slip flow of Casson fluid over a nonlinearly stretching sheet with convective boundary condition. Scientific Reports. 2017; 7: 1–19.
  50. 50. Naduvinamani NB, Shankar U. Thermal-diffusion and thermo-diffusion effects on squeezing flow of unsteady MHD Casson fluid between two parallel plates with thermal radiation. Sadhana. 2019; 44: 1–16.
  51. 51. Bakr AA. Effects of chemical reaction on MHD free convection and mass transfer flow of a micropolar fluid with oscillatory plate velocity and constant heat source in a rotating frame of reference. Communications in Nonlinear Science and Numerical Simulation. 2011; 16(2): 698–710.
  52. 52. Ibrahim SM. Radiation effects on mass transfer flow through a highlyporous medium with heat generation and chemical reaction. SRN Computational Mathematics. 2013; 765408: 1–9.
  53. 53. Ullah I, Bhattacharyya K, Shafie S, Khan I. Unsteady MHD mixed convection slip flow of Casson fluid over nonlinearly stretching sheet embedded in a porous medium with chemical reaction, thermal radiation, heat generation/absorption and convective boundary conditions. PLOSE One. 2016; 11: 1–35.
  54. 54. Song J, An W, Wu Y, Tian W. Neutronics and Thermal Hydraulics Analysis of a Conceptual Ultra-High Temperature MHD Cermet Fuel Core for Nuclear Electric Propulsion. Frontiers in Energy Research. 2018; 6: 1–10.
  55. 55. Saadati H, Hadad K, Rabiee A. Safety margin and fuel cycle period enhancements of VVER-1000 nuclear reactor using water/silver nanofluid. Nuclear Engineering and Technology. 2018; 50(5): 639–647.
  56. 56. Noor NAM, Shafie S, Admon MA. Effects of viscous dissipation and chemical reaction on MHD squeezing flow of Casson nanofluid between parallel plates in a porous medium with Slip Boundary Condition. Eur. Phys. J. Plus. 2020; 135: 1–24.
  57. 57. Jaluria Y. Natural Convection Heat and Mass Transfer. Michigan: Elsevier Science & Technology; 1980.
  58. 58. Nadeem S, Akbar NS. Peristaltic flow of a Jeffrey fluid with variable viscosity in an asymmetric channel. Zeitschrift fur Naturforschung. 2009; 64(11): 713–722.
  59. 59. Anderson JD. Computational Fluid Dynamics. New York: McGraw-Hill Education; 1995.
  60. 60. Prasad KV, Datti PS, Vajravelu K. MHD mixed convection flow over a permeable non-isothermal wedge. Journal of King Saud University—Science. 2013; 25(4): 313–324.
  61. 61. Lund LA, Omar Z, Khan I, Kadry S, Rho S, Mari IA, et al. Effect of Viscous Dissipation in Heat Transfer of MHD Flow of Micropolar Fluid Partial Slip Conditions: Dual Solutions and Stability Analysis. Energies. 2019; 12(24): 1–17.
  62. 62. Noor NAM, Shafie S, Admon MA. MHD Squeezing Flow of Casson Nanofluid with Chemical Reaction, Thermal Radiation and Heat Generation/Absorption. Journal of Advanced Research in Fluid Mechanics and Thermal Sciences. 2020; 68: 94–111.
  63. 63. Noor NAM, Shafie S, Admon MA. Unsteady MHD squeezing flow of Jeffrey fluid in a porous medium with thermal radiation, heat generation/absorption and chemical reaction. Physica Scripta. 2020; 95: 1–13.
  64. 64. Noor NAM, Shafie S, Admon MA. Impacts of chemical reaction on squeeze flow of MHD Jeffrey fluid in horizontal porous channel with slip condition. Physica Scripta. 2021; 96: 1–18.