Sara I Abdelsalam1,2, M M Bhatti3. 1. Basic Science, Faculty of Engineering, The British University in Egypt Al-Shorouk City Cairo 11837 Egypt sara.abdelsalam@bue.edu.eg siabdelsalam@caltech.edu siabdelsalam@yahoo.com. 2. Division of Chemistry and Chemical Engineering, California Institute of Technology Pasadena CA 91125 USA. 3. Shanghai Institute of Applied Mathematics and Mechanics, Shanghai University Shanghai 200072 China.
During the recent years, the study of peristaltic flow has become an increasing interest for various researchers due to its efficient phenomena for the transport of fluids in different biological systems. It is a phenomenon in which a sinusoidal wave arises due to the proportional shrinkage and relaxation of smooth muscles in a human body. In particular, the peristaltic flow is involved in the transport of urine through the kidney to the bladder, transport of cilia, locomotion of spermatozoa (in the male reproductive tract), motion of chyme (in the gastrointestinal tract), motion of ova (in the fallopian tubes), and in the vasomotion of tiny blood vessels. In industry, peristaltic mechanism is very beneficial in transporting different biological fluids such as sanitary and corrosive fluids. For this purpose, many devices, such as heat lung machines, roller pumps, cell separators, and finger pumps, have been introduced in biomedical engineering that follow the fundamentals of peristaltic mechanism. Due to the numerous applications of the peristaltic flow, several researchers investigated the mechanism of peristalsis in different media. For instance, Mekheimer[1] studied the motion of couple stress fluid due to the peristaltic waves through a non-uniform channel. Later, Mekheimer[2] extended the previous problem influenced by the magnetic field considering blood as a couple stress fluid and obtained the exact solutions. The nonlinear peristaltic flow under the influence of a magnetic field through a uniform planar conduit was discussed by Hayat et al.[3] Ellahi et al.[4] investigated the peristalsis of three-dimensional motion of a non-Newtonian fluid in a rectangular canal. The peristaltic flow of the non-Newtonian Williamson fluid with compliant walls was investigated by Ellahi et al.[5] He further analysed the impacts of the wall tension and damping and obtained the series solution with the help of the perturbation technique. Nadeem et al.[6] studied the three-dimensional peristaltic motion of a Jeffrey fluid in a duct having flexible walls and obtained the exact solutions. Mekheimer et al.[7] studied the influence of the relaxation time of a Maxwell fluid together with the MHD peristaltic transport in a microchannel. Some further similar investigations on this topic can be found in the literature.[8-10]In fluid dynamics, a new branch has been introduced namely nanofluid dynamics, which has many applications in biology, medical science, energetics, and engineering processes. Nanotechnology introduces the creation and usage of numerous substances having the nanoscale dimensions from 1 to 100 nm. Basically, a nanofluid (NF) is a fluid that is amalgamated by scattering the nanoparticle (NP) in the base fluid such as body fluids, natural/artificial lubes, and water. Choi[11] was the first who initially concluded that the impact of the nanofluid phenomena was to enhance the energy performance. Although the basic concept of NFs was introduced in the 19th century by a well-known scientist James Clark (a Scottish theoretical physician), later the term nanofluid was introduced officially by Choi. Lee et al.[12] investigated the room-temperature thermal conductivity of NFs as a new class of fluids that was organized by dispersing NPs in water and ethylene glycol (EG). The non-Newtonian NF was investigated by Ellahi et al.[13] with Reynolds' model and Vogel's model using the homotopy method. Nanofluid particles are made up of metals, oxides, nitrides, or carbides having very small diameters (<100 nm). The base fluids can be the EG, lubricants, propylene glycol (PG), biofluids, coolants, emulsion, water, or silk fibroin. Fig. 1 shows the combination of different NPs with base fluids. Nanofluids are applicable and helpful in the understanding of various phenomena such as the enhancement/minimization of the magnitude of heat transfer systems, minimal clogging, miniaturization of the systems, and microchannel clogging.
Fig. 1
Combination of NPs with different base fluids.
Various researchers studied the combined mechanism of peristalsis with NF through different geometrical aspects. For instance, Akbar et al.[14] studied the peristaltic NF flow in an irregular tube. Akbar et al.[15] examined numerically the peristalsis of Williamson fluid in an asymmetric channel. Nadeem et al.[16] discussed the NF peristalsis in an eccentric conduit with heat and mass transfer. Ellahi et al.[17] studied theoretically the peristaltic mechanism of Prandtl nanofluid through a rectangular duct. Nadeem et al.[18] presented a mathematical formulation for the peristaltic motion of a non-Newtonian fluid with NPs. The peristaltic flow of NF having carbon NPs through a permeable channel under the influence of an induced magnetic field was investigated by Akbar et al.[19] Few more relevant studies can be found in the literature.[20-22]The existing fashion in the applications of MHD is towards the strong magnetic fields to take the effect of electromagnetic force into consideration. Consequently, the Hall and ion slip effects are crucial since they have exceptional influence on the current density. The Hall effect and ion slip effect have numerous applications, especially if incorporated with heat transfer, such as in Hall accelerators, refrigeration coils, heating elements, MHD accelerators, and power generators. Moreover, the study of the influence of the magnetic field along with the Hall and ion slip effects on the blood flow in an artery has been found to be very helpful and applicable in magnetic resonance angiography (MRA). This helps to create the images of arteries to explore the existence of stenosis or any other conditions in the arteries of the brain, abdomen, thorax, and kidneys. Magnetic resonance imaging (MRI) is also involved in other applications that involve pumping of blood, hyperthermia, cancer therapy, and magnetic drug targeting. Since the implementation of these applications provides exclusive capabilities to improve the mechanism of peristalsis uses, researchers have devoted much effort towards studying the peristaltic nanofluid with the magnetic field in different conduits. For instance, El Koumy et al.[23] studied the peristaltic motion of a Maxwell fluid under the influence of a strong magnetic field and accordingly the Hall effect through a conduit. Asghar et al.[24] investigated the simultaneous effects of Hall and ion slip along with the ohmic and viscous heating on the peristaltic motion through different ducts. Hayat et al.[25] considered the Hall and ion slip effects on the peristaltic phenomenon of a non-Newtonian Carreau–Yasuda fluid model. Abbasi et al.[26] considered the peristaltic motion of a silver–water nanofluid with the Hall and ion slip effects. They, in addition, considered the ohmic heating and wall characteristics such as the tension of elasticity and damping phenomenon. More comprehensive treatments about magnetic field models can be found in the outlined ref. 27–33.Considering the abovementioned discussion, the primary motivation of the present study was to extend our interest in studying the peristaltic motion of a hyperbolic tangent fluid with the effects of Hall and ion slip through a non-uniform channel taking the chemical reaction into consideration. To the best of our knowledge, this model has not been investigated in any of the referenced state-of-the-art reviews before. The system of equations describing the problem is formed by following the approach of the long wavelength and creeping flow regime. The resulting governing nonlinear partial differential equations have been solved by means of the HPM (homotopy perturbation method). The impacts of all the emerging parameters have been discussed in details with the help of the graphs.
Mathematical formulation
We considered the peristaltic flow of a blood-like incompressible, hyperbolic tangent, and electrically conducting NF under the effect of an externally applied magnetic field. A flow through the 2D non-uniform channel is induced due to the propagation of sinusoidal waves along its walls. The hydromagnetic flow of the nanofluid is considered unsteady and irrotational. We choose the Cartesian coordinate system in a way that x̃-axis is taken along the channel length and ỹ-axis is normal to it, as shown in Fig. 2. The geometry of the peristaltic walls can be described bywhere
Fig. 2
Schematic of the flow problem.
The generalized form of Ohm's law taking the Hall and ion slip effects into consideration can be written asSolving eqn (2), we obtainandwhere βe = ωeτe.The equations of motion governing the flow along with the thermal energy, continuity, and nanoparticle fraction for the blood NF can be written as[34]The stress tensor for hyperbolic tangent fluid is defined asIn the abovementioned equation, we have considered η∞ = 0 and Accordingly, the stress tensor can be rewritten asDefining the dimensionless quantities asWe followed the creeping flow proposition such that the half-width of the conduit was taken small as compared to the peristaltic wavelength. We further speculated that the Reynolds number is low. These assumptions are extensively used in many peristalsis analyses.[30,32,35-38] These approximations are considered in many biological tracts such as in the transport of enzymes to the duodenum. Using eqn (12) in eqn (3)–(10), we obtained the reduced system of equations in the following form:with corresponding boundary conditions aswhere
Method of solution
In this section, we attempted to solve the aforementioned non-linear couple of partial differential equations by means of the HPM. The homotopy for eqn (13)–(15) can be written asThe linear operators L1, L2, and L3 are suggested in the next formsMoreover, we defined the initial guess for the abovementioned linear operators asDefining the following expansionsBy substituting eqn (27)−(29) into eqn (18)−(20) and matching the like powers of q̃, a linear system of differential equations, along with their corresponding boundary conditions, was obtained. With reference to the scheme of HPM, we deduced the solution as q̃ → 1, and we obtainedThe solutions of the temperature profile, nanoparticle concentration, and velocity profile can simply be written aswhere.The instantaneous volume flow rate can be determined through the expressionThus, we can obtain the expression of the pressure gradient, dp/dx, after solving the latter equation. Hence, the dimensionless forms of pressure increase, Δp, and friction force, Δf, by the wall are given bywhere L is the non-uniform channel length.
Numerical results and discussion
In this part, we have discussed the theoretical significance of the developed physical expressions that are involved in the problem based on the current study. A Mathematica toolbox has been used to explore the outcomes arising due to the existence of the Brownian parameter Nb, thermophoresis parameter N, chemical reaction parameter γ, ion slip parameter βi, Hall parameter βe, basic density Grashof number GrF, thermal Grashof number GrT, power law index n, magnetic parameter M, Weissenberg number We, and the average time flow Q̄ into the flow field. More specifically, we investigated their influence on the distributions of temperature θ, concentration Φ, and velocity u, as well as the pressure increase Δp and friction force Δf. In the subsequent figures, the red, blue, and green coloured curves represent the variations of the given variable with the indicated parameter in an ascending order. Same is applicable to the variations of the solid and dashed lines where the solid line indicates a smaller value of the parameter under consideration. We consider that Q(x,t) is the instantaneous volume flow rate, which is cyclic having (x − t) cycle and hence can be written aswhere Q̄ expresses the average of the time flow over one wave cycle.Fig. 3−6 provide insight into the changes in the behaviour of the temperature and concentration distributions on the nanofluid that occur due to changes in the values of the chemical reaction parameter γ, the thermophoresis parameter Nt, and the Brownian parameter Nb. Fig. 3 examines the dependence of θ that is plotted with y for different values of Nb and Nt. The figure shows that the temperature profile θ is semi-parabolic, and there is a considerable increase in θ upon increasing both Nb and Nt. The possible reason is that the Brownian motion causes the nanoparticles to rearrange, forming a blend, which increases the thermal conductivity. Further, eqn (14) and (15) show that the temperature profile is commensurate with the thermophoresis parameter. The latter result may be important in the event that a treatment requires an increase in the temperature of the tissues such as in the case of magnetic hyperthermia treatment where the major aim of hyperthermia is to increase the temperature of malignant tissues above 42 °C. On the other hand, it is shown in Fig. 4 that θ is reduced with an increase in γ for distinct values of Nb. Further, it is shown that as Nb increases, the impact of γ becomes more significant. In this case, it must be taken into account that an increase in the Brownian parameter may result in an increase in the chemical reaction of the nanofluid temperature profile.
Fig. 3
Temperature distribution for various values of Nb and Nt at We = 0.01, M = 1, GrT = 0.5, GrF = 0.6, γ = 0.1, and n = 2.
Fig. 4
Temperature distribution for various values of Nb and γ at Nt = 0.5, We = 0.01, M = 1, GrT = 0.5, GrF = 0.6, and n = 2.
Fig. 5 depicts the influence of Nb and Nt on the concentration of nanofluid, Φ. It is noticed that Nb has a decreasing effect on Φ for different values for Nt, whereas Nt shows a quite opposite effect on Φ for various values of Nb. It is further shown that the Brownian parameter increases the impact of the thermophoresis in the nanofluid concentration profile. Fig. 6 illustrates the behaviour of Φ with various values of the chemical reaction γ at different values of the Brownian parameter. In this figure, it is shown that for a small value of Nb (=1), the chemical reaction seems to weakly affect Φ in the narrow part of the channel where y ∈ [0, 0.48], whereas γ tends to reduce the concentration profile afterwards. Conversely, for a larger value of Nb (=1.5), there is an obvious dual behaviour of γ on the concentration profile. That is, the concentration is enhanced in the region of y > 1.3, where γ enhances the fluid density, whereas it is reduced in the region of y < 1.3 due to the reduction in viscosity.
Fig. 5
Concentration distribution for various values of Nt and Nb at We = 0.01, M = 1, GrT = 0.5, GrF = 0.6, γ = 0.1, and n = 2.
Fig. 6
Concentration distribution for various values of Nb and γ at Nt = 0.5, We = 0.01, M = 1, GrT = 0.5, GrF = 0.6, and n = 2.
The two-dimensional and three-dimensional behaviour of the velocity profile are displayed in Fig. 7−10 for distinct values of the emerged parameters. Fig. 7 helps to elaborate the effects of the slip parameter, βi, and the Hall parameter, βe, on the velocity distribution. This figure reveals that the flow is accelerated for progressive values of βi till a certain turning point at y = 0.8 of negligible slip effect from which the flow decreases substantially afterwards. Similarly, the impact of the Hall parameter is seen to accelerate the flow till the same point before it begins to lag. Fig. 8 demonstrates the influence of the basic density Grashof number GrF and thermal Grashof number GrT on the velocity profile. Physical interpretation of the behaviour of GrF suggests that beyond a certain critical point (y = 0.78), the flow of the nanofluid decelerates with an increase in GrF. This has been predictable since a reduction in the Grashof number in the narrow part of the channel implies an increase in the viscosity causing deceleration in the velocity profile and vice versa. Conversely, the influence of GrT is observed to reduce the velocity distribution till y = 0.78 from which the velocity increases markedly afterwards.
Fig. 7
2D (a) and 3D (b) velocity distribution as a function of y for various values of βi and βe at We = 0.01, M = 1, GrT = 0.5, GrF = 0.6, γ = 0.1, and n = 2. For the 3D velocity profile: (i) yellow shading: βe = 0.5, (ii) green shading: βe = 2.
Fig. 8
2D (a) and 3D (b) velocity distribution as a function of y for various values of GrT and GrF at βe = 0.5, βi = 0.6, We = 0.01, M = 1, γ = 0.1, and n = 2. For the 3D velocity profile: (i) yellow shading: GrF = 1.5, (ii) green shading: GrF = 2.
Fig. 9 illustrates the impact of the power law index n and the magnetic parameter M on the velocity profile. The examination shows that the flow of nanofluid decelerates when the value of M increases in the narrow part of the channel where y ∈ [0, 0.82]. The effect is quite opposite afterwards where the fluid flow is seen to be substantially increasing with an increase of M. It is also noticed that in the narrow part of the channel, the magnitude of velocity is higher in the absence of the magnetic parameter. Thus, the nanofluid velocity can be reduced in this part by the application of a strong magnetic field on the flow. Contrariwise, the influence of n on the velocity distribution is seen to be increasing in the narrow part and decreasing afterwards. Fig. 10 helps to demonstrate the influence of We and Q̄ on the velocity profile. The inspection of this graph reveals that the flow decelerates for progressive values of We in the narrow part of the channel until y = 0.82, where its effect on u is negligible. Later on, the behaviour of the flow is reversed with an increase in We. However, the velocity profile is seen to be remarkably accelerated with an increment in Q̄.
Fig. 9
2D (a) and 3D (b) velocity distribution as a function of y for various values of M and n at βe = 0.5, βi = 0.6, We = 0.01, GrT = 0.5, GrF = 0.6, and γ = 0.1. For the 3D velocity profile: (i) yellow shading: n = 1.7, (ii) green shading: n = 1.8.
Fig. 10
2D (a) and 3D (b) velocity distribution as a function of y for various values of We and Q̅ at βe = 0.5, βi = 0.6, We = 0.01, GrT = 0.5, GrF = 0.6, M = 1, n = 2, and γ = 0.1. For the 3D velocity profile: (i) yellow shading (below): Q̅ = 0.1, (ii) green shading Q̅ = 0.2.
To study the influence of the pertinent parameters on the pressure increase Δp, Fig. 11−14 have been plotted. Fig. 11 demonstrates the effect of the slip parameter βi and the Hall parameter βe on the pressure increase. It is deduced from this graph that Δp is reinforced with escalating both βi and βe. It is further noticed that Δp is almost unperturbed by the variations in βi and βe in the interval t ∈ [0.22, 0.3], whereas the pressure increase attains its maximum value subsequently at t = 0.34. The fact that pressure increase is small in some intervals can be interpreted as the flow can facilely pass without imposition of large pressure, whereas to retain the same flux, large pressure is required. Fig. 12 and 13 are shown to investigate the effects of GrF, GrT, M, and n on the pressure increase Δp. A close look to the graphs reveals that Δp increases with the progressive values of GrT, whereas it decreases with an increase in GrF and M with a maximum value occurring at t = 0.32. It is also shown that M = 0 causes n to have a decreasing effect on the pressure increase till it reaches a critical value of M (= 0.5) where n weakly affects Δp. Afterwards, the behaviour of Δp is totally reversed to be increasing with n at a higher value of M. It is recognized that the pressure increase is higher in the absence of M. This phenomenon symbolizes the fact that the pressure can be controlled by suitably applying the magnetic field on the flow. This is an important factor in the use of magnetic field in physiology as any abrupt change in the intensity of the applied magnetic field can cause severe changes in the systolic/diastolic readings of the patient exposed to the magnetic field. Fig. 14 shows the influence of We and Q̄ on Δp. Evidently, We has an increasing effect on Δp, whereas Q̄ has a decreasing effect on it. Further, the pressure increase attains its maximum value at t = 0.35.
Fig. 11
Pressure increase for various values of βe and βi at We = 0.01, M = 1, γ = 0.1, GrT = 0.5, GrF = 0.6, and n = 2.
Fig. 12
Pressure increase for various values of GrF and GrT at We = 0.01, M = 1, γ = 0.1, βi = 0.6, βe = 0.5, and n = 2.
Fig. 13
Pressure increase for various values of n and M at We = 0.01, γ = 0.1, GrT = 0.5, GrF = 0.6, βi = 0.6, βe = 0.5, and n = 2.
Fig. 14
Pressure increase for various values of We and Q̅ at We = 0.01, γ = 0.1, GrT = 0.5, GrF = 0.6, βi = 0.6, βe = 0.5, M = 1, and n = 2.
Fig. 15−18 provide a perspective of the influence of the parameters under consideration on the friction force Δf. Fig. 15 provides insight into the effect of the slip parameter βi and the Hall parameter βe on Δf. It is seen from this figure that Δf decays prominently with the progressive values of βi and βe. It is further observed that Δf is minimum at t = 0.34. Different variations are observed for the effects of the basic Grashof number GrF and thermal Grashof number GrT on Δf throughout the domain t ∈ [0, 1], as seen in Fig. 16. It is observed that GrF has an increasing effect on the friction force, whereas GrT has a decreasing effect on it. It is also noticed that Δf attains its minimum at t = 0.33. Fig. 17 helps to demonstrate the influence of the Hartmann number and power law index on Δf. It is concluded that Δf increases with an increase in M. Moreover, it is shown that M = 0 causes n to have a decreasing effect on the friction force till it reaches a critical value of M (= 0.5) where n weakly affects Δf. Afterwards, the behaviour of Δf is totally reversed to be increasing with n at a higher value of M. In addition, it is seen that the friction force is lower in the absence of M, and the minimum value of Δf takes place at t = 0.33. Fig. 18 depicts the variations in Δf due to the changes in We and Q̄. Obviously, the friction force is seen to be enhanced noticeably for the higher values of Q̄, whereas it decays with an increase in We. It is also shown that Δf attains its minimum at t = 0.33.
Fig. 15
Friction forces for various values of βe and βi at We = 0.01, M = 1, γ = 0.1, GrT = 0.5, GrF = 0.6, and n = 2.
Fig. 16
Friction forces for various values of GrT and GrF at We = 0.01, M = 1, γ = 0.1, βi = 0.6, βe = 0.5, and n = 2.
Fig. 17
Friction forces for various values of M and n We = 0.01, γ = 0.1, GrT = 0.5, GrF = 0.6, βi = 0.6, and βe = 0.5.
Fig. 18
Friction forces for various values of We and Q̅ at We = 0.01, γ = 0.1, GrT = 0.5, GrF = 0.6, βi = 0.6, and βe = 0.5.
Conclusions
In this study, the peristaltic flow of a blood-like non-Newtonian (hyperbolic tangent) NF is investigated in a non-uniform channel to study the mathematical results under an external magnetic field, nanoparticle concentration, chemical reaction, Hall current, and ion slip conditions. The governing equations along with the boundary conditions are modelled under the long wavelength assumption. The solutions are obtained analytically using the homotopy perturbation technique, and the physical interpretation of the pertinent parameters is discussed. Diagrammatic sketches are given for the physical expressions with the relevant parameters considered in the flow field. The primary findings can be outlined as follows:(i) An increment in the Brownian parameter Nb causes an increase in the thermal conductivity θ, but causes a decrease in the concentration profile Φ.(ii) For a relatively large value of Nb, there is a dual behaviour of the chemical reaction γ on the concentration profile.(iii) The effect of the thermophoresis Nt is to enhance θ and Φ.(iv) The chemical reaction γ has a decreasing effect on both θ and Φ.(v) Unlike the effect of the basic density Grashof number GrT, Hartmann number M, and Weissenberg number We, the slip parameter βi, Hall parameter βe, thermal Grashof number GrF, and power law index n serve to boost the velocity distribution markedly before a certain critical point.(vi) The pressure increase Δp is enhanced with an increase in βi, βe, We, and GrT, whereas it decreases with an increase in GrF, M, and Q̄.(vii) In the absence of the magnetic field, the pressure increase attains the highest value, and the friction force attains the lowest value.(viii) There is a critical point of M at which the behaviours of the pressure increase and friction force are reversed for the progressive values of n.(ix) Contrary to the influence of βi, βe, We, and GrT on the friction force Δf, GrF, M, and Q̄ are seen to enhance Δf prominently.(x) Setting n = 0, We = 0, M = 0, and GrT = GrF = 0 in our analysis, the intrinsic equations governing the flow of Gupta[39] are recovered.(xi) Our study agrees with that reported for Newtonian fluid by Srivastava and Srivastava[36] in the absence of the magnetic field, basic and thermal Grashof numbers, and Weissenberg number.(xii) Upon solving our model in plane or axisymmetric geometries for n = 0, We = 0, M = 0, and GrT = GrF = 0 implies to a consistent physical situation as discussed by Shapiro and Jaffrin.[35](xiii) Upon adopting the analysis of Mekheimer[2] for the Newtonian fluids, our fundamental equations will be in a perfect match with his after setting n = 0, We = 0, M = 0, and GrT = GrF = 0 in our analysis.(xiv) Choosing n = 0, M = 0, and γ = 0 in the present investigation, our system of equations coincides with that of Abbas et al.[40] for the Newtonian fluids.(xv) If n, βi, and βe vanish in the leading equations that govern the current flow field, the system is reduced to that of Rashidi et al.[41] if performed in a non-porous medium in the absence of the radiation parameter, heat source/sink, and Weissenberg number.
Conflicts of interest
There are no conflicts to declare.Velocity components (m s−1)Cartesian coordinate (m)Pressure in fixed frame (N m−2)Wave amplitude (m)Width of the channel (m)Wave velocity (m s−1)Half width at the inletPrandtl numberReynolds numberTime (s)Basic density Grashof numberThermal Grashof numberBrownian motion parameterThermophoresis parameterConstantPower law indexMagnetic field (in Tesla)Weissenberg numberVolume flow rate (m3 s−1)Temperature (K) and concentrationTemperature at the center and at the wallNanoparticle fraction at the center and at the wallPerturbation parameterHartmann numberAcceleration due to gravity (m s−2)Brownian diffusion coefficient (m2 s−1)Thermophoretic diffusion coefficient (m2 s−1)Electric field (V m−1)Fluid velocity (m s−1)Current density (A m2).Chemical reaction parameterNanofluid thermal conductivity (W m K−1)Heat source/sink parameterViscosity of the fluid (N s m−2)Dimensionless temperature profileNanoparticle concentrationElectrical conductivity (S m−1)Second invariant tensor.Wavenumber (m−1)Effective heat capacity of nanoparticle (J/K)Nanofluid kinematic viscosity(m2 s−1)Nanoparticle mass density (kg m−3)Fluid density (kg m−3)Fluid density at the reference temperature (T0) (kg m−3)Volumetric expansion coefficient of the fluidHeat capacity of fluid (J/K)Wavelength (m)Amplitude ratioInfinite shear rate viscosityZero shear rate viscosityTime constantCyclotron frequencyElectron collision timeIon slip parameterHall parameter.