Hydrodynamic analysis of the magnetic field dependent viscous fluid flow and thermosolutal convection between rotating channels.

According to research, exposing a person to a magnetic field enhances blood flow and minimizes their risk of suffering a heart attack. Ferrohydrodynamics is the study of fluid motion mechanics that is affected by strong magnetic polarisation forces (FHD). Ferrofluids may transmit heat in a variety of ways by using magnetic fluids. This behaviour is demonstrated by liquid-cooled speakers, which utilise less ferrofluid to prevent heat from reaching the speaker coil. This modification boosts the coil's ability to expand, which enables the loudspeaker to create high-fidelity sound. It is investigated how the fluid dynamics of spinning, squeezing plates are affected by thermosolutal convection and a magnetic field dependent (MFD) viscosity. Standard differential equations are used to represent the equations of the modified form of Navier Stokes, Maxwell's, and thermosolutal convection. The magnetic field, modified velocity field equations, and thermosolutal convection equations all yield suitable answers. Additionally computed and thoroughly detailed are the MHD torque and fluid pressure that are imparted to the top plate. To create a technique with quick and certain convergence, the resulting equations for uniform plates are solved using the Homotopy Analysis Method (HAM) with appropriate starting estimates and auxiliary parameters. The validity and reliability of the HAM outcomes are shown by comparing the HAM solutions with the BVP4c numerical solver programme. It has been found that a magnetic Reynolds number lowers the temperature of the fluid as well as the tangential and axial components of the velocity field. Additionally, when the fluid's MFD viscosity rises, the axial and azimuthal components of the magnetic field behave in opposition to one another. This study has applications in the development of new aircraft take-off gear, magnetorheological airbags for automobiles, heating and cooling systems, bio-prosthetics, and biosensor systems.

Introduction

The study of the mechanics of fluid motion that is impacted by powerful magnetic polarisation effects is known as ferrohydrodynamics (FHD). Ferrofluids may transmit heat in a variety of ways by using magnetic fluids. This behaviour is demonstrated by liquid cooled speakers, which utilise less ferrofluid to prevent heat from reaching the speaker coil[1]. This modification boosts the coil’s ability to expand, which enables the loudspeaker to create high fidelity sound. A drop of ferrofluid going through the body can be used to deliver medications to a precise area in the human body using magnetic fluids[2], which can be propelled by a magnetic field. Lenz’s law states that an electric current develops as it moves through a magnetic field, and the conductor also develops its own magnetic field. Since the currents are observed by the movement of the fluid passing through the magnetic field, the Lorentz force acts on the mass and increases the speed of the. In MHD, the field affects movement and vice versa. The concept is hence non-linear[3,4]. Researchers have shown that being exposed to a magnetic field enhances blood flow to the body and lowers the chance of a heart attack. In an experiment conducted by Rongjia Tao and Ke Huang of the University of Michigan and Temple University, respectively, it was discovered that a minute’s exposure to a magnetic field perpendicular to the direction of blood flow caused a small sample of blood in a tube viscera to compress blood viscosity by 33 to 4.75cp. Within 200 minutes of no more contact to the fluid, the viscosity may have grown noticeably to 5.4cp, which is still well within the acceptable range of 7cp.

A substance or metal cube made of two or more homogeneous or heterogeneous metals or non-metallic alloys is called an alloy. These alloys mixed mixes with unique characteristics or employed a heat transfer medium. Steel, brass, phosphor bronze, and solder are a few examples of alloys. Aerospace research, powder technology, hip joint replacement, and surgical implantation are just a few of the industries that employ alloys. The three types of alloys are alpha alloys, beta alloys, and alpha beta alloys. An industrially developed titanium alloy combining aluminium and titanium derivatives is called Alpha Alloys. Due to their changed heat transfer bulk, these alloys compromise the quality of the fuse bulk, and as a result, the alpha alloy is used more frequently in temperature application operations. Beta alloy is mostly utilised in cold rolling of sheets and welding systems. Elshekh et al.[5] studied magnetic squeeze film flow from the magnetic field among rotating disk generated effects, considering as,with velocity componentsdecrease the responsiveness of the radial and azimuthal magnetic fields and measure the response to the accompanying disc rotation for a value of the bachelor number. In the flow of the posterior portion of a vertical infinite plate in a rotating structure, Mutua et al.[6] investigated the Stokes challenge. They came to the conclusion that any or all factors can affect the temperature and speed of a liquid. As a result, the rate of heat transmission along the axis and the skin’s resistance are both impacted. For both free cooling and heating in plate variations, the speed profile, magnetic parameters M, and Ec number EC enhance M convection. MHD was investigated by Seth et al.[7] and shown to be amplified with rising values of time and Hall parameters, mass transfer, heat transfer in flow, and mass transfer with porous flat plates, but to be amplified with decreasing values of devolution parameters and magnetic field parameters. In order to maintain heat radiation, Victor[8] employed the finite component approach to examine the volatile MHD free convection cutlet transfer inside the two vertical porous plates. Due to the radiation characteristics and the predental number, it was discovered that the temperature velocity was higher than the velocity. In other words, the fluid particle’s grashof numbers and magnetic properties were unaffected. The volatile MHD QUETT flow in two infinitely parallel porous plates was studied using the coefficient[9] in a bent magnetic field with rate of heat transfer in temperature. The bottom plate was regarded as permanent and porous. He found that when the magnetic field grows, the fluid’s velocity decreases. In the presence of an externally applied magnetic field, Verma[10] examined the lubrication of a magnetic fluid twisting film between two approaching surfaces. He assumed that the applied magnetic field M has the following components. He assumed that the applied magnetic field M has components of the form:trying to squeeze flux with spinning disks is also examined in depth by[11-14] under influences of changing magnetic fields. For the simulation of the magnetic film problem, Rashidi et al.[15] used the Pade DTM hybrid method, which demonstrated good convergence consistency, and versatility. The objective of this research is to examine the existence of exponential heating decay and transversal magnetic field[16-18] for the influence of different electrical conductivity on a free convection flow of a viscous electrically conducive fluid and heat transfer through an isothermal non conductive vertical plat. The heat transfer and aerodynamic forces rate may be conveniently adjusted by applying an appropriate magnetic field. Morley examine the impact of MHD given in[19] known as blanket. The blanket is placed among both magnetic field and plasma’s spindles, absorbing neutrons that transform its energies into heat, which is subsequently taken away by an appropriate cooling, restricting neutrons from entering the magnets and therefore preventing radiation destruction. For drag reduction several fascinating improvements have been presented in the domains of MHD propulsion and remote energy deposition. The flow is represented by the system of nonlinear equations which includes fundamental ohmic law, Maxwell, momentum and energy equations[24-27]. The system of nonlinear equations given in[28] was solved by finite difference technique. Flow of non-Newtonian electrically conducting fluid is an essential phenomena, since we are dealing with fluid flow, which displays various patterns in the effect of magnetic forces in most real circumstances. In such situations, the feature of the flow MHD must be taken into account. Siddiqui et al.[29] utilizing homotopy perturbation method to obtain solutions for squeezed MHD two dimensional flow among parallel plates. Newtonian nonisothermal fluid flow among unsteady squeezed permeable disks for changeable magnetic-field will be investigate. The PDEs well be transform into nonlinear coupled ODEs using similarity transformations. The formed nonlinear DEs describing the flow characteristics within the geometry under observation will also analyze through numerical and analytical approaches. Analyses will also be conducted amongst the solutions. The convergence of results will be explored as well. Flow pattern will be explained graphically[30-32] under the impact of non-dimensional factors. Nazeer et al.[33] stadied the thermal transport of two-phase physiological flow of non Newtonian fluid through an inclined channel with flexible walls. Yassen et al.[34], investigated the theoretical study of transport of MHD peristaltic flow of fluid under the impact of viscous dissipation. Nazeer et al.[35] investigated the heat transmission in a magnetohydrodynamic multiphase flow induced by metachronal propulsion through porous media with thermal radiation. Zubaidi et al.[36] studied the numerical study of squeezing flow past a Riga plate with activation energy and chemical reactions: effects of convective and second-order slip boundary conditions. Quraishi et al.[37] the influence of radially magnetic field properties in a peristaltic flow with internal heat generation: Numerical treatment. Awan et al.[38] examined the numerical treatments to analyze the nonlinear radiative heat transfer in MHD nanofluid flow with solar energy. Raja et al.[39] studied the integrated intelligent computing application for effectiveness of Au nanoparticles coated over MWCNTs with velocity slip in curved channel peristaltic flow. Awan et al.[40] examined the backpropagated intelligent computing networks for 3D nanofluid rheology with generalized heat flux. Majeed et al.[41] studied the mathematical analysis of MHD CNTs of rotating nanofluid flow over a permeable stretching surface. Shahzad et al.[42] investigated the effects of magnetohydrodynamics flow on multilayer coatings of Newtonian and non Newtonian fluids through porous inclined rotating channel. Zeeshan et al., examined the radiative bioconvection nanofluid squeezing flow between rotating circular plates: Semi-numerical study with the DTM Padé approach.

This demonstrates the aforementioned heavy emphasis on fluid motion in a magnetic field with high rotational velocity in an MFD viscosity. MFD thermosolutal convection also has applications in various fields such as aerospace-technology, powder technologies, liquid cooled speaker, pump building, auto magneto rheological suspension system, new aviation system landing gears, biological prosthetics, hip joints replacement processes and the surgical implantation procedure. Due to its essential scientific and practical importance, it is sought to characterize the impact of rotation and how MFD thermosolutal convection, MFD viscosity impacts magnetization of the fluids. The information available at the theme revealed that the squeezing-flow MFD viscosity between rotating plates for a viscous fluid under the MFD thermosalutal and externally applied magnetic field in Cartesian coordinates was never reported and is the very first study to carried out in the literature. Hence, the suggestework is the best approach toward such problems and is a way of motivation for researchers bringing a new idea of studying the flow between unsteady rotating parallel plates. In the following sections, the problem is studied, examined and explained using tables and graphs.

Formulation of the problem

In this article the viscous fluid is suppose to be axisymmetric and uncompressed, squeezed between parallel plates, distanced by a length of , in which l represents a plate separation at [20,22,23]. The two plates being squeezed unless they meet at for , and they are separated for . Both plates rotate with a distinct angular velocities of and , where , representing the angular speeds of bottom and top plates correspondingly[20-22]. The bottom plate is fixed, while the top plate moves near or apart from the bottom plate. As compared to magnetic forces, Electrical forces are significantly less and thus neglected in the current research problem. Applied magnetic field inside the fluid generates the induced magnetic field which can be written as:Wherever the magnetic permeability of the external and the inner media among these plates is and , , , and , . On the bottom plate[5] the aforementioned magnetic field parameters are zero. Incompressible fluids is considered, with a variable viscosity provided by where , represents the fluid’s viscosity in the absence of a magnetic-field[25]. The alteration viscosity-coefficient are assumed to be isotropic, i.e. = = . Hence, in the component form can be written as , and . Furthermore, M, B, and H may linked by the equation , in which M is the magnetization when the magnetic field is H and the vacuum magnetic-permeability . The influence of shearing dependency on viscosity is not examined because it has a minor influence on a wide spin and high field monodispersive system. A linear change of magnet viscosity was initially employed as a small field variation, thus where is uniform magnetic-field and is the magnetization when the magnetic field is . The equations which governing the flow and heat/mass transfers in viscous fluid are:

Continuity equation: of Navier Stokes equation: of Navier Stokes equation: of Navier Stokes equation: of Magnetic field equation : of Magnetic field equation: of Magnetic field equation:Energy equation:[25]Concentration equation:[25]where is the electrical conductivity, is specific heat at constant volume and magnetic field, T is fluid temperature, C is solute concentration, M is magnetization, thermal conductivity, , are solute conductivity, is reference density p is pressure, is variable dynamic viscosity and U is the velocity of the fluid.

Boundary conditions

The boundary conditions are chosen as:where is the thermal diffusivity, the specific heat at constant pressure is , u, v and w are respectively the radial, azimuthal and axial components of the velocity field, D is diffusion coefficient, , , and are temperature and concentration at lower and upper plate respectively.

Due to the non-linearity factor and selection of the boundary conditions for the set of coupled equations (1–10), it is too complicated to solve it in the form of partial differential equations (PDEs). To solve this issue, the set of PDEs have been converted into the set of ordinary differential equations (ODEs) by introducing the Lie group of similarity transformations[26]. These similarity transformations are chosen as, ,      ,      ,      ,      ,      ,      ,      ,      .

Equation (1) is identically satisfied and Eqs. (2–9) takes the following formand the boundary conditions are reduced towhere is Eckert number, is strength of magnetic field in z direction, is strength of magnetic field in direction, is squeeze Reynolds number, is Pyromagnetic coefficient and is the salinity magnetic coefficient, is rotational Reynolds number, is Bachelor number, is magnetic Reynolds number, is the relative angular velocity of discs, is magnetic field dependent viscosity parameter, is Prandtl number, is Hartmann number, .

Approximate analytical solution

The analytic method HAM is used to solve system of Eqs. (9–15). Due to HAM, the functions , , , , and can be expressed, by a set of base functions as:[43]where , , , , and are the constant coefficients to be determined. Initial approximations are chosen follows:The auxiliary operators are chosen aswith the following propertieswhere , , , , , , , , , , , , and are arbitrary constants.

The Zeroth order deformation problems can be obtained as:The nonlinear operators of Eqs. (9–14) are defined aswhere is an embedding parameter, , , , , and are the nonzero auxiliary parameter and , , , , and are the nonlinear parameters. For and 1, we haveso we can say that as varies from 0 to 1, , , , , , varies from initial guesses , , , , and to exact solution , , , , and respectively.

Taylor’s series expansion of these functions yields:it should be noted that the convergence of above series strongly depends upon , , , , and .

Assuming that these nonzero auxiliary parameters are chosen so that Eqs. (35–40) converges at . Therefore one can obtain

Differentiating the deformation Eqs. (35–40) with respect to and putting , we havesubject to the boundary conditionswhereand .

Finally, the general solution of Eqs. (59–64) can be written asand so the exact solution , , , , and becomes

Optimal convergence control parameters

It must be remarked that the series solutions (59–64) contain the nonzero auxiliary parameters , ,, , , and which determine the convergence region and also rate of the homotopy series solutions. To obtain the optimal values of , , , , and here the so called average residual error defined by Liao[16] were used as:Due to Liao[16]where is the total squared residual error. Total average squared residual error is minimized by employing Mathematica package BVPh 2.0[16].

Physical quantities of interest

Torque exerted by fluid on the plates

The frictional moment or torque which the fluid exerts on the upper plate is given[26]but so the above equation becomes,For the lower plate the corresponding result is

The pressure or normal force of the fluid on plates

According to HAM za et al.[26], the pressure or the normal force which the fluid exerts on the upper plate is given as:where in for conditions on the side of the plate and p(r, 1, t) denotes the pressure at the edge of the disc at time t[26]. Let us assume that and using Eq. (8), we getwhere . using Eqs. (25), (26) , we havewhere , Which is the dimensionless pressure on the upper plate. The positive or negative numerical values of will be according the force acting by the fluid on the upper plate is in the positive or negative direction of the z-axis respectively.

Error analysis

Error assessments are performed to ensure that the derived system of ODEs is analytically accurate with the lowest residual errors. The validity of HAM techniques utilising the highest residual error is also investigated in this study. Approximations of the order of approximation are used in the study. Error analyses in Fig. 1 and tabulated data in Tables 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 and 12 are offered to validate the legitimacy of the findings for the various included physical parameters. Figure 2 is drawn to show the 3D-view of the velocity field components, magnetic field components and heat/mass distribution.

Figure 1

Error profile of , , , , and for fixed values of , , , , , , , , , , and .

Table 1

Optimal values of convergence control parameters versus different orders of approximation with fixed values of , , , , , , , , , , and .

Order	h_f	h_g	h_h	h_k	h_θ	h_Ψ	εmt
2	-182.02	-2.8096	-1.8780	-1.8760	-1.8468	-1.8428	2.0114×10^-7
3	-190.87	-1.9148	-0.9691	-0.9812	-0.9728	-0.9712	1.1235×10^-10
4	-191.92	-1.8999	-0.9952	-1.0173	-0.9665	-0.9665	4.3529×10^-13
5	-190.97	-1.9025	-0.9882	-0.9810	-0.9530	-0.9498	-2.5293×10^-15
6	-191.79	-1.9184	-0.9828	-0.9456	-0.9498	-0.9412	-6.6548×10^-17

Table 2

Total residual error for different order of approximations taking fixed values of , , , , , , , , , , and .

m	ϵ^f_m	ϵ^g_m	ϵ^h_m	ϵ^k_m	ϵ^θ_m	ϵ^Ψ_m
2	1.231×10^-7	2.864×10^-6	2.245×10^-6	2.752×10^-5	8.534×10^-9	5.234×10^-6
5	3.342×10^-14	4.127×10^-13	3.421×10^-14	2.531×10^-13	3.274×10^-17	3.819×10^-15
8	4.173×10^-21	3.223×10^-20	3.425×10^-21	4.426×10^-20	4.342×10^-24	3.124×10^-21
11	7.421×10^-28	8.658×10^-27	1.238×10^-27	3.153×10^-27	3.446×10^-31	2.456×10^-29c
14	2.642×10^-34	2.451×10^-32	1.395×10^-32	3.327×10^-32	2.148×10^-35	3.469×10^-33
17	2.351×10^-34	1.478×10^-32	1.868×10^-32	2.231×10^-32	2.424×10^-35	4.324×10^-33
20	2.351×10^-34	1.478×10^-32	8.223×10^-33	2.925×10^-32	1.941×10^-35	5.970×10^-33
23	2.351×10^-34	1.478×10^-32	1.068×10^-32	2.741×10^-32	1.760×10^-35	5.691×10^-33
26	2.351×10^-34	1.478×10^-32	1.377×10^-32	3.049×10^-32	1.700×10^-35	5.845×10^-33
30	2.351×10^-34	1.478×10^-32	1.038×10^-32	2.432×10^-32	1.941×10^-35	6.654×10^-33

Table 3

Computations for ,, , , , with , , , , , , , , , , , and various values of .

η	HAM result					BVP4c result
	f(η)	g(η)	h(η)	θ(η)	Ψ(η)	f(η)	g(η)	h(η)	θ(η)	Ψ(η)
0	0	1	0	1	1	0	1	0	1	1
0.1001	0.0261	0.9267	0.0899	0.8987	0.8974	0.0261	0.9264	0.0899	0.8989	0.8977
0.2002	0.0668	0.8497	0.1911	0.7987	0.7978	0.0669	0.8499	0.1912	0.7989	0.7978
0.3003	0.1258	0.7703	0.2716	0.6989	0.6954	0.1259	0.7704	0.2717	0.6989	0.6955
0.4004	0.1839	0.7100	0.3654	0.5986	0.5944	0.1839	0.7101	0.3655	0.5987	0.5945
0.5005	0.2582	0.6579	0.4619	0.4983	0.4938	0.2582	0.6579	0.4619	0.4983	0.4938
0.6006	0.3324	0.6245	0.5723	0.3994	0.3946	0.3325	0.6246	0.5724	0.3995	0.3947
0.7007	0.3985	0.5875	0.6765	0.2993	0.2951	0.3986	0.5876	0.6766	0.2994	0.2952
0.8008	0.4622	0.5568	0.7845	0.1994	0.1959	0.4623	0.5569	0.7846	0.1995	0.1959
0.9009	0.4981	0.5315	0.8965	0.0997	0.0976	0.4982	0.5316	0.8966	0.0998	0.0977
1	0.5	0.5	1	0	0	0.5	0.5	1	0	0

Table 4

Convergence of HAM solution for different orders of approximation for , , , , and when , , , , , , , , , , and .

m	f^′′(0)	-g^′(0)	-h^′(0)	-k^′(0)	-θ^′(0)	-Ψ^′(0)
1	3.30248492	0.76935826	-0.76249352	-0.86198432	1.01629436	1.00261843
6	3.30248445	0.76935839	-0.76249348	-0.861984256	1.01629429	1.00261862
12	3.30248472	0.76935856	-0.76249352	-0.861984261	1.01629435	1.00261869
18	3.30248472	0.76935856	-0.76249352	-0.861984261	1.01629435	1.00261869
24	3.30248472	0.76935856	-0.76249352	-0.861984261	1.01629435	1.00261869
30	3.30248472	0.76935856	-0.76249352	-0.861984261	1.01629435	1.00261869

Table 5

Fluid pressure and torques on upper plate with , , , , , , , , , , and various values of .

ξ_sq	HAM result		BVP4c result		%Error = |N.R-H.RN.R|
	F_pres	g^′(1)	F_pres	g^′(1)	F_pres	g^′(1)
0.1	19.2226	-0.4558	19.2225	-0.4562	5.20224×10^-6	2.19346×10^-3
0.5	5.2753	-1.7021	5.2752	-0.4724	1.89566×10^-4	5.87475×10^-4
1	3.4494	-2.6614	3.4493	-0.5128	2.89914×10^-4	3.75728×10^-4
1.5	2.8340	-3.3410	2.8339	-0.3400	3.52871×10^-4	2.99303×10^-4

Table 6

Fluid pressure and torques on upper plate with , , , , , , , , , , and various values of .

ξ_vis	HAM result		BVP4c result		%Error = |N.R-H.RN.R|
	F_pres	g^′(1)	F_pres	g^′(1)	F_pres	g^′(1)
0.1	1.9722	-5.5855	1.9721	-5.5856	5.07074×10^-4	1.79032×10^-4
0.5	3.6515	-2.4989	3.6514	-2.4990	2.73868×10^-4	4.00160×10^-4
1	5.6661	-1.5789	5.6660	-1.5790	1.76491×10^-4	6.33312×10^-4
1.5	7.6712	-1.1632	7.6711	-1.1633	1.30359×10^-4	8.59623×10^-4

Table 7

Fluid pressure and torques on upper plate with , , , , , , , , , , and various values of .

ξ_a	HAM result		BVP4c result		%Error = |N.R-H.RN.R|
	F_pres	g^′(1)	F_pres	g^′(1)	F_pres	g^′(1)
0.1	19.2024	-0.4561	19.2023	-0.4562	5.20771×10^-6	2.19202×10^-3
0.5	6.9653	-0.4723	6.9652	-0.4724	1.43571×10^-4	2.11775×10^-3
1	-34.0176	-0.5127	-34.0175	-0.5128	2.93966×10^-6	1.95008×10^-3
1.5	-268.29	-0.3399	-268.28	-0.3400	3.72745×10^-4	2.94204×10^-3

Table 8

Computations for , , , and with , , , , , , , , , , and various values of .

ξ_vis	HAM result					BVP4c result
	f^′′(0)	-g^′(0)	-h^′(0)	-k^′(0)	-θ^′(0)	f^′′(0)	-g^′(0)	-h^′(0)	-k^′(0)	-θ^′(0)
0.1	3.9532	1.8329	0.9267	0.7936	1.3785	3.9532	1.8329	0.9267	0.7936	1.3785
1	3.7789	0.3857	0.7837	1.037	1.3087	3.7789	0.3857	0.7837	1.037	1.3087
2	3.5985	0.0963	0.5893	1.4687	1.3097	3.5985	0.0963	0.5893	1.4687	1.3097
3	3.3789	0.0065	0.3784	1.7499	1.3001	3.3789	0.0065	0.3784	1.7499	1.3001

Table 9

Computations for , , , and with , , , , , , , , , , and various values of .

ξ_mag	HAM result					BVP4c result
	f^′′(0)	-g^′(0)	-h^′(0)	-k^′(0)	-θ^′(0)	f^′′(0)	-g^′(0)	-h^′(0)	-k^′(0)	-θ^′(0)
0.5	3.1389	0.2598	0.7994	1.4829	3.2687	3.1389	0.2598	0.7994	1.4829	3.2687
1	3.1372	0.2598	0.5839	1.3799	1.3087	3.1372	0.2598	0.5839	1.3799	1.3087
2	3.1351	0.2599	0.4981	1.2911	1.3097	3.1351	0.2599	0.4981	1.2911	1.3097
3	3.1314	0.2599	0.2389	1.1999	1.3001	3.1314	0.2599	0.2389	1.1999	1.3001

Table 10

Computations for , , and with , , , , , , , , , , and various values of .

ξ_ec	HAM Result					BVP4c result
	f^′′(0)	-g^′(0)	-h^′(0)	-θ^′(0)	-Ψ^′(0)	f^′′(0)	-g^′(0)	-h^′(0)	-θ^′(0)	-Ψ^′(0)
1	3.0719	0.2981	0.1819	1.1371	1.0988	3.0719	0.2981	0.1819	1.1371	1.0988
2	3.0719	0.2981	0.1819	1.2811	1.1381	3.0719	0.2981	0.1819	1.2811	1.1381
3	3.0719	0.2981	0.1819	1.3101	1.2847	3.0719	0.2981	0.1819	1.3101	1.2847
4	3.0719	0.2981	0.1819	1.3919	1.3899	3.0719	0.2981	0.1819	1.3919	1.3899

Table 11

Computations for , , , and with , , , , , , , , , , and various values of .

ξ_pm	HAM Result					BVP4c result
	f^′′(0)	-g^′(0)	-h^′(0)	-k^′(0)	-θ^′(0)	f^′′(0)	-g^′(0)	-h^′(0)	-k^′(0)	-θ^′(0)
0.1	3.1301	0.1299	0.8102	0.6989	1.1127	3.1301	0.1299	0.8102	0.6989	1.1127
0.5	3.1301	0.1299	0.8102	0.6989	1.2994	3.1301	0.1299	0.8102	0.6989	1.2994
1	3.1301	0.1299	0.8102	0.6989	1.3999	3.1301	0.1299	0.8102	0.6989	1.3999
2	3.1301	0.1299	0.8102	0.6989	1.5391	3.1301	0.1299	0.8102	0.6989	1.5391

Table 12

Computations for , , , and with , , , , , , , , , , and various values of .

ξ_sm	HAM result					BVP4c result
	f^′′(0)	-g^′(0)	-h^′(0)	-k^′(0)	-θ^′(0)	f^′′(0)	-g^′(0)	-h^′(0)	-k^′(0)	-θ^′(0)
0.1	3.9749	0.3769	0.3996	1.0929	0.9998	3.9749	0.3769	0.3996	1.0929	0.9998
0.5	3.9749	0.3769	0.3996	1.0929	1.0069	3.9749	0.3769	0.3996	1.0929	1.0069
1	3.9749	0.3769	0.3996	1.0929	1.0135	3.9749	0.3769	0.3996	1.0929	1.0135
1.5	3.9749	0.3769	0.3996	1.0929	1.0299	3.9749	0.3769	0.3996	1.0929	1.0299

Figure 2

Error profile of , , , , and for fixed values of , , , , , , , , , , and .

The greatest average residual-error of f(eta), g(eta), h(eta), k(eta), ta(eta), and psi(eta) using various approximating orders is shown in Fig. 3. In sub-figures it is evident that the error is virtually constantly decreased to . Table 1 shows the overall residual error for various approximation orders, with fixed values, Table 2 shows comparison of HAM Results to numerical values for velocity field, heat and mass distributions variables. The BVP4c Results for , , , , , and given in Table 3 provides further conformation of the authenticity for our results. Results converge almost on a approximation order. As the order of approximation is increased, the result converges to an exact solution. The results of HAM and BVP4c are compared in Tables 6, 7, 8, 9, 10, 11 and 12 for varied values of and fixed values of other associated parameters.

Figure 3

Geometry of the problem.

Results and discussions

To examine and analyse how the magnetic field and squeezing phenomena affect flow rate in the presence of MFD temperature, optimise the system’s heating/cooling process, lessen fluid turbulence, and scale up flow tracers. using the laminar, unstable, and incompressible MFD viscous fluid flow produced by compressing plates under the influence of no-slip velocity and convective surface boundary conditions. to look into the flow for mass and heat transfer and to offer an analytical and numerical study of it. For the cases of rotating and compressing plates, the influence of the different flow parameters is addressed graphically. The influence of the flow parameters is given for the vertical velocity component , horizontal velocity compnent , magnetic field components and temperature and mass transfer variations, respectively. In this section the impact of rotational parameter , squeeze Reynold number , Prandtl number , Hartman number and all other parameters are respectively analyzed and discussed in detail. Figures 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22 and 23 and Tables 5, 6, 7, 8, 9, 10, 11 and 12 depict effect of the natural parameters that appear during mathematical modelling of this problem. Dimensionless parameters are analysed for , , , , , , and on the velocity field components , , , components of the magnetic field , and heat/mass distributions , respectively. To model physically realistic flows by HAM , representative values are used.

Figure 4

Impact of on and for fixed values of , , , , , , , , , and .

Figure 5

3D-view of the impact of on and for fixed values of , , , , , , , , , and .

Figure 6

Impact of on and for fixed values of , , , , , , , , , and .

Figure 7

3D-view of the impact of on and for fixed values of , , , , , , , , , and .

Figure 8

Impact of on and for fixed values of , , , , , , , , , and .

Figure 9

3D-view of the impact of on and for fixed values of , , , , , , , , , and .

Figure 10

Impact of on and for fixed values of , , , , , , , , , and .

Figure 11

3D-view of the impact of on and for fixed values of , , , , , , , , , and .

Figure 12

Impact of on and for fixed values of , , , , , , , , , and .

Figure 13

3D-view of the impact of on and for fixed values of , , , , , , , , , and .

Figure 14

Impact of on and for fixed values of , , , , , , , , , and .

Figure 15

3D-view of the impact on and for fixed values of , , , , , , , , , and .

Figure 16

Impact of on and for fixed values of , , , , , , , , , and .

Figure 17

Impactof on and for fixed values of , , , , , , , , , and .

Figure 18

3D-view of the Impact on , for fixed values of , , , , , , , , , and .

Figure 19

Impact of of and for fixed values of , , , , , , , , , and .

Figure 20

3D-view of the impact and for fixed values of , , , , , , , , , and .

Figure 21

Impact of on for fixed values of , , , , , , , , , and .

Figure 22

Impact of on and for fixed values of , , , , , , , , , and .

Figure 23

3D-view of and for fixed values of , , , , , , , , , and .

Figures 4 and 5 depicts the effect of on and . The ratio of the fluid flux to magnetic diffusivity is described by . As a result, along streamlines this parameter is critical in evaluating magnetic field diffusion[15]. As seen in Fig. 4, a decrease in the magnetic diffusion enhances the axial and radial velocity distributions neat the lower plate. The minimum value of and appears near the lower plate. As fluid travels from the central region toward upper plate, the influence of on and decreases due to the decrease in the magnetic flux. Figure 5 is made to depict the 3D-view of the velocity components. The influence of on and is shown in Figs. 6 and 7 . As seen in Fig. 6a, the minimum value of appears near the upper plate due to the resistance of magnetic diffusion to the horizontal and vertical velocities, similarly in Fig. 6b, an increase in fluid flux or a reduction in magnetic diffusivity has a clear effect on . Figure 7 show the 3D-view of and for fixed values of other parameters. The effect of is also investigated for and concentration distribution. For this purpose, Figs. 8 and 9 are drawn. A decrease in the magnetic diffusion enhances the fluid viscosity will decrease near the upper plate, similarly both decrease strongly. Figure 9 is made to to observed this phenomenon in a 3D-view.

The effects of MFD viscosity on and are shown in Figs. 10 and 11. A higher viscosity coefficient means that a fluid has more flow resistance. Figure 10 shows that the flow of fluid near the plates reduces in the horizontal direction as the fluid’s viscosity increases because of an increase in the fluid’s flow resistance. It’s also important to note that fluid density begins to rise in the fluid domain’s centre area as it moves in the direction of the top plate. An increase in has an impact on azimuthal velocity, as seen in Fig. 11b. The highest value of has the greatest drop in at the upper disc. Figure 11 shows this phenomena in a 3D-dimensional view. Investigations on the magnetic field’s impact of MFD viscosity are also conducted. The and can be impacted by changes in MFD viscosity, as would be predicted. Figure 12 illustrates how and both display the opposing action of . This phenomena is depicted in Fig. 13 in a 3D-dimensional form.

One of the most significant physical parameter is the rotational Reynolds parameter. The effect of the rotational Reynolds number could be observed in Figs. 13 and 14. A rise in indicates either an increase in the lower plate angular velocity or a reduction in fluid viscosity. Expanding the value of , in vicinity of plates, the axial velocity reduces. This is because of the rise of angular velocity of the lower plate (both plates rotate in the same direction as pushes fluid in a radial direction, reducing fluid movement in the axial direction near the plates. As the fluid flows into the central region, the rotational strength of the plates decreases, and the velocity begins to increase in the axial direction after the central region, as seen in Fig. 14b. Figure 15 is made to observed this phenomenon in a 3D-view. Figure 16 depicts the effect of on the axial as well as azimuthal portion of the induced magnetic field, . Clearly increase in strengthens .

The impact of the dimensionless strength of the axial magnetic field is depicted in Figs 17 and 18. The influence of on and near plates has been seen to be rather small, but it becomes noticeable when fluid moves into the higher plate. starts to rise at the fluid domain’s centre, then falls as it gets closer to the top disc. Similar to this, drops as grows. The greatest decrease is visible towards the centre of the fluid domain. Figure 18 is designed to show the effect of from a three-dimensional perspective. Figures 19 and 21 show how has an effect on , , and . The applied magnetic field’s dimensionless azimuthal amplitude is . It is important to keep in mind that rising raises both and . This is due to the fact that at the top disc, where the angular velocities of the lower plate’s and are larger than one another, has a bigger impact. Figure 20 is made to display the impact of in a three dimensional view. The impact of Eckert number is examined in Figs. 22 and 23 on and . The Eckert number is the kinetic energy ratio to the fluctuation of the boundary layer in temperature used to characterize the heat dispersion. The figure shows that decreases each and and virtually provide the identical impact due to the increase in the kinetic energy of the fluid. Tables 5, 6, 7, 8, 9, 10, 11, 12 and 13 are made to numerically analyze the impact of various physical parameters on skin friction, magnetic flux, heat and mass fluxes. From these tables it may be observed that all findings correspond quite well with the findings of the BVP4c and HAM. The formula is used to calculate the error between HAM findings (H.F) and the findings that obtained numerically (N.F) by BVP4c. Tables 5, 6 and 7 illustrate the torque and pressure that the fluid exerts upon on top plate. It is noticed that increasing the distance between plates i.e reduces pressure and torque on the top plate. Elshekh et al. and Hughes et al.[5,24] respectively found a rise in fluid pressure on the top plate, which corresponds with their experimental and theoretical findings. Similarly Tables 6 and 7 investigate the effect of and . Also Table 8 depict that an increase in MFD viscosity decrease skin friction and slightly increase heat flux . Similarly effect of and is shown in Tables 9 and 10. Tables 11 and 12 are made to show effects of and on , , , , and .

Table 13

Torques on lower and upper plate with , , , , , and various values of .

ξ_sq	Rashidi et al.[15]		HAM result
	g^′(0)	g^′(1)	g^′(0)	g^′(1)
0.1	-1.08963506	-0.95987349	-1.08963672	-0.95989843
0.2	-1.17203765	-0.93844830	-1.17208879	-0.93846874
0.3	-1.25013649	-0.92615609	-1.25010991	-0.92612331
0.5	-1.39797797	-0.91295593	-1.39792988	-0.91296243
1	-1.73306821	-0.89280536	-1.73309746	-0.89282365
2	-2.28925762	-0.83902117	-2.28929872	-0.83903546

Concluding remarks

The computational model for the constitutive expressions of an unsteady Newtonian fluid is used to simulate the flow between the circular space of porous and squeezing discs in the form of equations. Using the mathematical model for the constitutive expressions of unstable Newtonian fluid, the flow between the circular space of porous and squeezing discs is described in the form of Eqs. (9–14) is subject to the boundary requirements given in Eq. (15). The series solution of the following equations is found using HAM: , , , , , , and for magnetic field components. The extraordinary stability and convergence properties of the HAM have been demonstrated. Through the use of HAM and BVP4c, these equations are compared for numerical investigations. Parametric analysis are carried out for the dimensionless parameters such as MFD viscosity parameter , magnetic Reynolds number , rotational Reynolds number , squeezing Reynolds number , axial magnetic force parameter , tangential magnetic force parameter and Eckert number . In the future, this problem could be investigated in the form of PDEs. The solution of the problem in the PDEs form will explore the physics of the problem in details.

Main upshots of this paper are presented as below:

As the fluid’s MFD viscosity rises, the axial and azimuthal components of the magnetic field behave in opposition to one another.

It has been discovered that a magnetic Reynolds number lowers the temperature of the fluid as well as the tangential and axial components of the velocity field.

With an increase in the bottom plate’s rotating speed, the fluid’s torque and pressure on the plates decrease. The findings of Hughes et al.[27] and Elshekh et al.[23] as well as theoretical and experimental findings are in agreement with this outcome.

For the magnetic field components, the effect of the MFD viscosity has been seen to follow a different pattern.

The Eckert number’s effects on and have been seen to be identical.

Heat and mass transfer coefficients, as well as the distribution of the generated magnetic field, are rising functions of the rotational Reynolds number .

It is also observed from Table 7 that increase in fluid viscosity decrease and increase .