Numerical Treatment for 3D Squeezed Flow in a Rotating Channel With Soret and Dufour Effects

This article examines magnetohydrodynamic three-dimensional (3D) squeezed ﬂow by a rotating permeable channel subject to Dufour and Soret impacts. Impact of viscous dissipation is also considered. An applied magnetic ﬁeld is considered subject to electrically conducting viscous ﬂuid. The change from the non-linear partial differential framework to the non-linear ordinary differential framework is assumed into position by utilizing appropriate variables. Governing differential frameworks are computed numerically by shooting method. Numerical results have been achieved by considering numerous values of emerging ﬂow parameters. Contributions of inﬂuential parameters on physical quantities are studied thoroughly. Surface drag coefﬁcients and mass and heat transport rates are also processed and examined. Furthermore, the concentration and temperature distributions are reduced for larger values of Soret number. The prime interest of presented study is to model and examine the Dufour and Soret aspects in concentration and energy expressions. To our knowledge, no such analysis has been addressed in the literature yet.


INTRODUCTION
Flow squeezed by two parallel surfaces is an attractive region of research embraced by recent researchers.The presence of squeezed flow in the designing applications having liquid metal grease, polymer, and sustenance ventures is in charge of this premium.Squeezing flow can be utilized in displaying the grease framework.Stefan [1] has started the exploration on lubrication approximation.Squeezed flow of power-law liquid between parallel disks is examined by Leider and Bird [2].Hamza and MacDonald [3] discussed the effect of suction/blowing in the squeezed flow.MHD unsteady squeezed flow between two parallel surfaces is discussed by Bhattacharyya and Pal [4].Fully developed free-convection micropolar fluid flow in a vertical channel is provided by Chamkha et al. [5].Rashidi et al. [6] introduced the investigation of axisymmetric and two-dimensional squeezed flows between parallel walls.Homotopic perturbation solution for MHD squeezed flow between parallel disks is addressed by Domairry and Aziz [7].Hayat et al. [8] explored squeezed flow of second-grade liquid between parallel disks.Three-dimensional squeezed flow in a rotating channel with lower stretchable permeable plate is examined by Munawar et al. [9].Freidoonimehr et al. [10] investigated solution in a rotating channel by taking three-dimensional squeezed nanofluid flow.Few relevant examinations on squeezed flows can be seen through attempts [11][12][13][14][15].

R E T R A C T E D
Simultaneous presence of mass and heat transfer in a moving liquid gives more intricate nature that the fluxes and driving potentials convey between them.It has been seen that temperature gradients and also concentration gradients can produce energy flux.The diffusion-thermo (Dufour) impact is characterized as the heat transport because of concentration gradient while the thermal-diffusion (Soret) impact is the mass transport because of temperature gradient.Mass and heat transport related examinations uncovered the smaller order of magnitude of the Dufour and Soret impacts when contrasted with the impacts of Fourier's and Fick's law and are neglected much of the time.These impacts are critical in nuclear waste disposal, hydrology, geothermal energy, petrology, etc. Soret impact is utilized for partition of isotope and in mixture between gases with almost small and medium sub-atomic weights (H 2 , He) and (N 2 , air) separately.Dufour impact can't be ignored in view of its considerable magnitude [16].Rashidi et al. [17] explained convective MHD flow by a rotating disk subject to diffusion-thermo and thermal-diffusion impacts.Dufour and Soret features in magnetohydrodynamic Casson liquid flow over an extending surface is proposed by Hayat et al. [18].Turkyilmazoglu and Pop [19] examined the Soret impact in natural convection unsteady flow subject to thermal radiation and heat generation.Properties of Dufour and Soret in buoyancydriven MHD flow by a stretching surface is addressed by Pal and Mondal [20].Hayat et al. [21] explored Dufour and Soret impacts in mixed convection peristaltic transport containing nanoliquid by considering slip and Joule heating.Some relevant examinations on Dufour and Soret effects can be seen through attempts [22][23][24][25].
Keeping the above discussion in mind, the present article is organized for magnetohydrodynamic three-dimensional (3D) squeezed flow of viscous liquid in a rotating permeable channel subject to Dufour and Soret effects.Viscous dissipation is also considered.Non-linear partial differential systems are simplified via appropriate transformations to the non-linear ordinary differential systems.Shooting technique is used in order to construct the numerical solution of non-linear flow problem.Salient features of fluid flow and mass and heat transfer are further examined.

PROBLEM FORMULATION
We examine unsteady three dimensional squeezing flow of viscous liquid between two parallel plates which are separated by a distance √ ν (1 − γ t) /a.Upper plate at y = h (t) = √ ν (1 − γ t) /a is moving with velocity − γ 2 ν a(1−γ t) and the lower permeable plate at y = 0 is stretched with velocity (ax/1 − γ t) in which t is always less than 1/γ (see Figure 1).

An angular velocity
= ωj/1 − γ t has been utilized by the fluid and the channel to rotate about y−axis while the lower plate sucks the flow with velocity −V 0 /1 − γ t.Magnetic field of strength (B 0 /1 − γ t) is employed in y−direction [26][27][28][29][30][31][32].The thermophysical characteristics of under discussion fluid are taken to be constant.The governing equations in rotating frame of  reference are defined by in which τ stands for Cauchy stress tensor, V for velocity field, J for magnetic flux and B for current density.
In component form, the resulting expressions of mass, momentum, energy, and concentration in the absence of thermal  radiation are [9,25]: ∂w ∂t     boundary conditions are [9,12]:

∂C ∂t
in which stretching rate of the lower plate, suction, injection, temperature and concentration at the lower plate are symbolized  by a, V 0 > 0, V 0 < 0, T 0 and C 0 respectively.Selecting [9,12]:

R E T R A C T E D
Pressure gradient is eliminated from Equations 4) and (5 and Equation ( 3) is now verified while Equations (4) − (9) have been reduced to

R E T R A C T E D
Here Ec stands for Eckert number, for rotation parameter, S for suction/blowing parameter, M for magnetic number, Pr for Prandtl number, S q for squeezing number, Sr for Soret number, Sc for Schmidt number and Df for Dufour number.These parameters are stated by Expressions of friction drag coefficients are given by and Friction drag coefficients in non-dimensional scale are and The non-dimensional forms of local Nusselt and Sherwood numbers are stated below: in which local Reynolds number is signified by Re x = U 0 x/ν.

SOLUTION METHODOLOGY
By utilizing acceptable boundary conditions on system of the equations, a numerical solution is established employing NDSolve in Mathematica.Shooting method is employed via NDSolve.This method is very friendly in case of small step-size featuring negligible error.As a consequence, both the y and x varied uniformly by a step-size of 0.01 [26,[33][34][35].

NUMERICAL RESULTS AND DISCUSSION
This section elaborates impacts of different flow variables like magnetic number M, Eckert number Ec, squeezing number S q , Soret number Sr and Dufour number Df on velocities F ′ (η) and G (η), temperature θ (η) and concentration φ (η). Figure 2 depicts that how squeezing number S q affects velocity distribution F ′ (η).It is seen that the velocity distribution F ′ (η) is depreciated near the permeable plate where suction impacts are superior.Pressure which increases the flow has been developed due to the mobility of the upper plate toward a stretching permeable plate.Mass conservation constrain has been mollified by incrementing velocity distribution F ′ (η) near the upper plate.Figure 3 demonstrates that how the velocity distribution is get effected by magnetic parameter M.Here the velocity distribution is decayed by enhancing M for (0 ≤ η ≤ 0.4) while opposite trend is seen when (0.4 ≤ η ≤ 1).Physically by increasing magnetic parameter M, the velocity and its gradient are decreased.Therefore the mass conservation constraint is satisfied by introducing the same mass flow rate.In MHD flow, we consider that a cross flow behavior is generated by increasing fluid velocity in the central region which results in balancing of fluid velocity decrement in the wall regions.Figure 4 depicts variation in velocity distribution G (η) for varying squeezing number S q .It is noted that by enhancing S q , an enhancement is appeared in velocity distribution G (η) and this increment is more prominent at central zone of channel.Figure 5 is displayed to depict the influence of magnetic number M on velocity distribution G (η). Greater values of magnetic parameter M constitutes a lower velocity distribution G (η). Figure 6 depicts the impact of squeezing number S q on temperature distribution θ (η).An increment in squeezing number S q leads to weaker temperature θ (η).

R E T R A C T E D
we examined that friction drags are reduced for the greater values of M and S q while it enhances by incrementing S. Table 2 presents numerical estimations of mass and heat transport rates for varying S, Ec, , M, S q , Sr, Pr, Df , and Sc.Here we concluded that mass and heat transport rates are higher when larger estimations of S q are considered.

FIGURE 1 |
FIGURE 1 | Schematic diagram of the problem.

FIGURE 2 |
FIGURE 2 | Sketch of velocity distribution F ′ (η) for squeezing number S q .

FIGURE 4 |
FIGURE 4 | Sketch of velocity distribution G (η) for squeezing number S q .
) Here u, v and w represent velocities in x−, y− and z−directions respectively whereas µ, ρ, ν (= µ/ρ), p and σ stand for dynamic viscosity, density, kinematic viscosity, pressure and electrical conductivity respectively, D for mass diffusion coefficient, α * m for thermal diffusivity, T for temperature, k T for thermal-diffusion, c s for concentration susceptibility, c p for specific heat, C for concentration and T m for fluid mean temperature.Subjected

TABLE 2 |
Numerical data for local Nusselt number and local Sherwood number for varying S, , S q , M, Ec, Sr, Df, Sc, and Pr.

TABLE 1 |
Numerical data for friction drag coefficients for varying S, , S q and M.
Figures 7, 8are sketched to examine that how temperature distribution is get effected by Soret Sr and Dufour Df numbers.From these figures, temperature distribution is weaker for larger Sr while higher trend is seen for larger Df .Soret and Dufour numbers respectively.Here we seen that an increase in Sr and Df show decreasing behavior for concentration distribution φ (η).Table1is tabulated in order to analyze the numerical computations of friction drag coefficients F ′′ (1) and G ′ (1) for varying , S, M, and S q .Here

Table 3
By increasing the Eckert parameterEc, an enhancement is observed in temperature distribution θ (η) .Opposite trend is seen in temperature distribution θ (η) for larger estimations of Dufour and Soret numbers.Concentration distribution φ (η) is decreasing functions Dufour and Soret numbers.The Dufour impact is characterized as the heat transport because of concentration gradient while the Soret impact is the mass transport because of temperature gradient.Heat and mass transport related examinations uncovered the smaller order of magnitude of the Dufour and Soret impacts when contrasted with the impacts of Fourier's and Fick's law and are neglected much of the time.These impacts are critical in nuclear waste disposal, hydrology, geothermal energy, petrology, etc.The present work provides an inspiration for future developments on topic in the regimes of melting heat transfer, variable sheet thickness, Cattaneo-Christov heat flux and Joule heating.