Electromagnetohydrodynamic (EMHD) Flow in a Microchannel with Random Surface Roughness

This study investigates the effect of small random transverse wall roughness on electromagnetohydrodynamic (EMHD) flow is in a microchannel, employing the perturbation method based upon stationary random function theory. An exact solution of a random corrugation function ξ, which is a measure of the flow rate deviated from the case without the roughness of two plates, is obtained by integrating the spectral density. After the sinusoidal, triangular, rectangular, and sawtooth functions that satisfy the Dirichlet condition are expanded into the Fourier sine series, the spectral density of the sine function is used to represent the corrugation function. Interestingly, for sinusoidal roughness, the final expression of the corrugation function is in good agreement with our previous work. Results show that no matter the shape of the wall roughness, the flow rate always decreases due to the existence of wall corrugation. Variations of the corrugation function and the flow rate strongly depend on fluid wavenumber λ and Hartmann number Ha. Finally, the flow resistance is found to become small, and the flow rate increases with roughness that is in phase (θ = 0) compared with the one that is out of phase (θ = π).


Introduction
Studies on dynamic micropumps based on electromagnetohydrodynamic (EMHD) flow and electroosmosis have received much attention in recent years.Micropumps are largely used in a diverse range of fields, including biological and chemical analysis [1], drug delivery [2], microelectronic cooling [3], DNA hybridization, and space exploration [4], among others.In particular, EMHD micropumps offer many advantages over other types of non-mechanical micropumps, including a simple fabrication process, bidirectional pumping ability, and the usability of medium conducting liquid.The mechanism of EMHD flow is based on the Lorentz force, which is generated by applying an electric current to the conducting liquid across the channel in the presence of a perpendicular magnetic field.EMHD flow in microchannels was first proposed by Jang and Lee [5].Ko and Dulikravich [6] studied non-reflective boundary conditions for EMHD flow simulations.Chakraborty and Paul [7] established a mathematical model to study the comprehensive influence of electromagnetic hydrodynamics on fluid flow control in parallel flat rectangular microchannels.In the absence of an electric field, Awatif et al. [8] analyzed the effect of magnetic force and moderate Reynolds number on MHD Jeffrey hybrid nanofluid through the peristaltic channel: application of cancer treatment.Qian and Bau [9] investigated the basic theory of EMHD in low conductivity solutions and described the various applications of EMHD.Asterios and Eugen [10] considered the electromagnetohydrodynamic (EMHD) free convection flow of a weakly conductive fluid, such as seawater, in an electromagnetic actuator.Buren et al. [11] employed a perturbation method to investigate the problem of EMHD flow in a microchannel with slightly corrugated walls.Tso and Sundaravadivelu [12] investigated the capillary flow between parallel plates in the presence of an electromagnetic field.Si and Jian [13] studied the EMHD micropump in Jeffrey fluid through two parallel corrugated wall microchannels.Wasim et al. [14] studied the dynamic irreversible process of the flow of magnetohybrid nanofluids of Poiseuille in microchannels.Faisal et al. [15] analyzed the effect of pressure gradient on MHD flow of a tri-hybrid Newtonian nanofluid in a circular channel.Bhatti et al. [16] studied simultaneous EMHD dissipative natural convection in microchannels containing porous media saturated with viscoelastic fluids.However, the above studies are based on smooth microchannel surfaces.In effect, almost all microchannel surfaces have a certain degree of roughness, either incurred during the fabrication process or because of the adsorption adhesion of other species, such as macromolecules [17].Thus, it is necessary to consider the flow in a rough microchannel.Additionally, surface roughness can be artificially created to support axial rotation or mixing [18,19].Thus, many researchers have studied the influence of wall roughness on laminar flow since the 1970s.Wang [20] initially investigated how roughness affected the Stokes flow between corrugated plates.Chu [21] studied the effect of wavy surface roughness on the fluid flow in an annulus under microscopic slip conditions.Tashtoush and Al-Odat [22] analyzed the effect of a magnetic field on heat and fluid flow on a wavy surface with a variable heat flux.Chen and Cho [18] studied the mixing characteristics of two-dimensional electroosmotic flow through cosine wavy microchannels.Chang et al. [23] studied the electroosmotic flow in sinusoidal wall microtubules with the perturbation method.Lei et al. [24] used boundary perturbation analysis to conduct theoretical research of electroosmotic pumping in bumpy microtubules.Su et al. [25] studied electro-osmotic flow through a microchannel with corrugated wavy walls under the Debye-Hückel approximation.Ng and Wang [26] studied the Darcy-Brinkman flow through a corrugated channel.Bujurke and Kudenatti [27] ignored the electric field and analyzed the influence of surface roughness on the extrusion film behavior between two rectangular plates and conductive fluid in the presence of transverse magnetic field.Ashmawy [28] investigated the effect of surface roughness on the rate of flow and mean velocity of two stress fluid flows through a sinusoidal corrugated tube.
In previous theoretical studies, wall roughness is generally represented by sine function or cosine function.However, in actual manufacturing, the roughness of the wall is often random.We hope that we can generalize sine or roughness to a more general case.Phan-Thien [29] first carried out the Stokes flow of incompressible Newtonian fluids between corrugated plates with random roughness, employing the theory of stationary stochastic process.The result was in good agreement with that obtained by Wang [21], when wall roughness reduces to a sinusoidal corrugation.Subsequently, Phan-Thien [30] generalized the above results to channels and pipes with parallel stationary random surface roughness.Faltas et al. [31] extended Phan-Thien's approach [29] to a sparse porous medium to study the effects of Basset-type slip boundary condition and impermeable walls.The results indicated that the slip parameter is a very important factor for almost perfect slip and for high permeability.They also noted that the type of corrugation has an effect on the flow rate and the pressure gradient.Faltas et al. [32] proposed a boundary perturbation method to study an improved micropolar Brinkman model, using a microwave grained cylindrical tube filled with a porous medium through which micropolar fluid flows.At the same time, a stationary stochastic model was used to simulate the surface roughness of the pipe.
However, to our knowledge, no one has studied the effect of random roughness on EMHD flow.Therefore, the main objective of this paper is to generalize the effect of sinusoidal corrugation on EMHD flow to random wall surfaces by using the perturbation method and some properties of stationary random function.The rest of the paper is organized as follows.In Section 2, the perturbation approach is used to construct and solve the EMHD equations defining surface roughness.In Section 3, we introduce the conditions that should be satisfied from for the definitions of random roughness and concrete roughness, and show four different types of roughness.In Section 4, the problem is studied parametrically and the effects of different parameters on the flow are discussed in detail.Section 5 provides a summary of the full text.

Materials and Methods
We consider the 2D EMHD flow of the viscous, incompressible, conductive Newtonian liquid confined to two fixed corrugated walls with an average distance of 2H parts.A Cartesian coordinate system (x, y, z) is adopted with the origin fixed at the middle of the channel.The flow variables are independent of y; therefore, the velocity vector has the form: In our model, the channel length L in the x-direction and width W in the y-direction are assumed to be much larger than the layer thickness, that is W, L >> 2H, and the flow is driven both by the pressure gradient and Lorenz force.The channel is subjected to a uniform magnetic field B * in the z direction and a uniform electric field E * in the y direction, which induce a Lorenz force J * × B * in x-direction, as shown in Figure 1.Where J * is the electrical current density and J * = σ(E * + u * × B * ) [11], σ is the electrical conductivity, u * is the liquid velocity vector with components u * , v * , w * in the x, y, z directions, respectively.
sinusoidal corrugation on EMHD flow to random wall surfaces by using the perturbation method and some properties of stationary random function.The rest of the paper is organized as follows.In Section 2, the perturbation approach is used to construct and solve the EMHD equations defining surface roughness.In Section 3, we introduce the conditions that should be satisfied from for the definitions of random roughness and concrete roughness, and show four different types of roughness.In Section 4, the problem is studied parametrically and the effects of different parameters on the flow are discussed in detail.Section 5 provides a summary of the full text.

Materials and Methods
We consider the 2D EMHD flow of the viscous, incompressible, conductive Newtonian liquid confined to two fixed corrugated walls with an average distance of 2H parts.A Cartesian coordinate system (x, y, z) is adopted with the origin fixed at the middle of the channel.The flow variables are independent of y; therefore, the velocity vector has the form: In our model, the channel length L in the x-direction and width W in the y-direction are assumed to be much larger than the layer thickness, that is W, L >> 2H, and the flow is driven both by the pressure gradient and Lorenz force.The channel is subjected to a uniform magnetic field B * in the z direction and a uniform electric field E * in the y direction, which induce a Lorenz force J * × B * in x-direction, as shown in Figure 1.Where J * is the electrical current density and J * = σ(E * + u * × B * ) [11], σ is the electrical conductivity, u * is the liquid velocity vector with components u * , v * , w * in the x, y, z directions, respectively.The dimensional equations of continuity and momentum can be expressed as: where p * is the pressure of the liquid, µ represents the dynamic viscosity, ρ stands for density, and is the gradient operator.The boundary conditions satisfy no-slip and no-penetration conditions [11], which means that the layer of fluid in The dimensional equations of continuity and momentum can be expressed as: concrete roughness, and show four different types of roughness.In Section 4, the problem is studied parametrically and the effects of different parameters on the flow are discussed in detail.Section 5 provides a summary of the full text.

Materials and Methods
We consider the 2D EMHD flow of the viscous, incompressible, conductive Newtonian liquid confined to two fixed corrugated walls with an average distance of 2H parts.A Cartesian coordinate system (x, y, z) is adopted with the origin fixed at the middle of the channel.The flow variables are independent of y; therefore, the velocity vector has the form: In our model, the channel length L in the x-direction and width W in the y-direction are assumed to be much larger than the layer thickness, that is W, L >> 2H, and the flow is driven both by the pressure gradient and Lorenz force.The channel is subjected to a uniform magnetic field B * in the z direction and a uniform electric field E * in the y direction, which induce a Lorenz force J * × B * in x-direction, as shown in Figure 1.Where J * is the electrical current density and J * = σ(E * + u * × B * ) [11], σ is the electrical conductivity, u * is the liquid velocity vector with components u * , v * , w * in the x, y, z directions, respectively.The dimensional equations of continuity and momentum can be expressed as: where p * is the pressure of the liquid, μ represents the dynamic viscosity, ρ stands for density, and is the gradient operator.The boundary conditions satisfy no-slip and no-penetration conditions [11], which means that the layer of fluid in * • u * = 0, (2) solve the EMHD equations defining surface roughness.In Section 3, we introduce the conditions that should be satisfied from for the definitions of random roughness and concrete roughness, and show four different types of roughness.In Section 4, the problem is studied parametrically and the effects of different parameters on the flow are discussed in detail.Section 5 provides a summary of the full text.

Materials and Methods
We consider the 2D EMHD flow of the viscous, incompressible, conductive Newtonian liquid confined to two fixed corrugated walls with an average distance of 2H parts.A Cartesian coordinate system (x, y, z) is adopted with the origin fixed at the middle of the channel.The flow variables are independent of y; therefore, the velocity vector has the form: In our model, the channel length L in the x-direction and width W in the y-direction are assumed to be much larger than the layer thickness, that is W, L >> 2H, and the flow is driven both by the pressure gradient and Lorenz force.The channel is subjected to a uniform magnetic field B * in the z direction and a uniform electric field E * in the y direction, which induce a Lorenz force J * × B * in x-direction, as shown in Figure 1.Where J * is the electrical current density and J * = σ(E * + u * × B * ) [11], σ is the electrical conductivity, u * is the liquid velocity vector with components u * , v * , w * in the x, y, z directions, respectively.The dimensional equations of continuity and momentum can be expressed as: where p * is the pressure of the liquid, μ represents the dynamic viscosity, ρ stands for density, and is the gradient operator.The boundary conditions satisfy no-slip and no-penetration conditions [11], which means that the layer of fluid in solve the EMHD equations defining surface roughness.In Section 3, we introduce the conditions that should be satisfied from for the definitions of random roughness and concrete roughness, and show four different types of roughness.In Section 4, the problem is studied parametrically and the effects of different parameters on the flow are discussed in detail.Section 5 provides a summary of the full text.

Materials and Methods
We consider the 2D EMHD flow of the viscous, incompressible, conductive Newtonian liquid confined to two fixed corrugated walls with an average distance of 2H parts.A Cartesian coordinate system (x, y, z) is adopted with the origin fixed at the middle of the channel.The flow variables are independent of y; therefore, the velocity vector has the form: In our model, the channel length L in the x-direction and width W in the y-direction are assumed to be much larger than the layer thickness, that is W, L >> 2H, and the flow is driven both by the pressure gradient and Lorenz force.The channel is subjected to a uniform magnetic field B * in the z direction and a uniform electric field E * in the y direction, which induce a Lorenz force J * × B * in x-direction, as shown in Figure 1.Where J * is the electrical current density and J * = σ(E * + u * × B * ) [11], σ is the electrical conductivity, u * is the liquid velocity vector with components u * , v * , w * in the x, y, z directions, respectively.The dimensional equations of continuity and momentum can be expressed as: where p * is the pressure of the liquid, μ represents the dynamic viscosity, ρ stands for density, and is the gradient operator.The boundary conditions satisfy no-slip and no-penetration conditions [11], which means that the layer of fluid in * p * + µ solve the EMHD equations defining surface roughness.In Section 3, we introd conditions that should be satisfied from for the definitions of random roughn concrete roughness, and show four different types of roughness.In Section 4, the p is studied parametrically and the effects of different parameters on the flow are di in detail.Section 5 provides a summary of the full text.

Materials and Methods
We consider the 2D EMHD flow of the viscous, incompressible, con Newtonian liquid confined to two fixed corrugated walls with an average distanc parts.A Cartesian coordinate system (x, y, z) is adopted with the origin fixed at the of the channel.The flow variables are independent of y; therefore, the velocity ve the form: In our model, the channel length L in the x-direction and width W in the y-d are assumed to be much larger than the layer thickness, that is W, L >> 2H, and the driven both by the pressure gradient and Lorenz force.The channel is subject uniform magnetic field B * in the z direction and a uniform electric field E * in the y di which induce a Lorenz force J * × B * in x-direction, as shown in Figure 1.Where electrical current density and J * = σ(E * + u * × B * ) [11], σ is the electrical conductivity, liquid velocity vector with components u * , v * , w * in the x, y, z directions, respective The dimensional equations of continuity and momentum can be expressed a where p * is the pressure of the liquid, μ represents the dynamic viscosity, ρ sta density, and is the gradient operator.The boundary con satisfy no-slip and no-penetration conditions [11], which means that the layer of where p * is the pressure of the liquid, µ represents the dynamic viscosity, ρ stands for density, and sinusoidal corrugation on EMHD flow to random wall surfaces by using the perturbation method and some properties of stationary random function.The rest of the paper is organized as follows.In Section 2, the perturbation approach is used to construct and solve the EMHD equations defining surface roughness.In Section 3, we introduce the conditions that should be satisfied from for the definitions of random roughness and concrete roughness, and show four different types of roughness.In Section 4, the problem is studied parametrically and the effects of different parameters on the flow are discussed in detail.Section 5 provides a summary of the full text.

Materials and Methods
We consider the 2D EMHD flow of the viscous, incompressible, conductive Newtonian liquid confined to two fixed corrugated walls with an average distance of 2H parts.A Cartesian coordinate system (x, y, z) is adopted with the origin fixed at the middle of the channel.The flow variables are independent of y; therefore, the velocity vector has the form: In our model, the channel length L in the x-direction and width W in the y-direction are assumed to be much larger than the layer thickness, that is W, L >> 2H, and the flow is driven both by the pressure gradient and Lorenz force.The channel is subjected to a uniform magnetic field B * in the z direction and a uniform electric field E * in the y direction, which induce a Lorenz force J * × B * in x-direction, as shown in Figure 1.Where J * is the electrical current density and J * = σ(E * + u * × B * ) [11], σ is the electrical conductivity, u * is the liquid velocity vector with components u * , v * , w * in the x, y, z directions, respectively.The dimensional equations of continuity and momentum can be expressed as: where p * is the pressure of the liquid, μ represents the dynamic viscosity, ρ stands for density, and is the gradient operator.The boundary conditions satisfy no-slip and no-penetration conditions [11], which means that the layer of fluid in * = ∂ ∂x i + ∂ ∂y j + ∂ ∂z k is the gradient operator.The boundary conditions satisfy no-slip and no-penetration conditions [11], which means that the layer of fluid in contact with a solid body has the same velocity as the body, and the fluid layer molecules do not penetrate the solid layer, are followed as: the upper and lower wavy walls are described by: When m(x) = n(x), we want to study the effect of phase difference on flow, and the plus (mins) sign denotes a situation in which the corrugations are in (half-period out of) phase.The amplitude of the corrugation ε << 1, and will be considered a perturbation parameter for the problems described below.m(x) and n(x) are stationary random functions, which allows all shapes of ripples to be generalized, and can be expressed as Fourier-Stieltjes integrals [33] (spectral representation of m(x) and n(x)): where i 2 = −1, and dZ 1 (λ), dZ 2 (λ) are interval random functions that satisfy [33] dZ 1 (λ) = dZ 2 (λ) = 0, (7a) where the overbar denotes a complex conjugate; δ ij is the Kronecker delta; and F where Ω jk (j,k = 1, 2) are the spectral densities of m(x) and n(x).
The Reynolds number is believed to be low, and the momentum balance equation nonlinear inertia can be ignored.We have introduced the stream function Ψ, which is given in terms of the velocity components, as follows: Supposing the flow is steady, the governing equations of a 2D EMHD flow can be expressed in the dimensional form as: Assuming the mean flow is purely in the x-direction, and the flow rate per unit width is known to be Q.Thus, utilizing H, Q/2, 1/H, and µQ/2H 2 as length, stream function, wave number, and pressure scales, respectively, the non-dimensionalization governing Equations ( 13) and ( 14) become: where Hartmann number Ha is the ratio of Lorentz force to viscous force, and the nondimensional parameter S measures the strength of the external electric field.The dimensionless no-penetration and no-slip boundary conditions are: where n is the normal direction of the channel walls, the normal derivative and n x and n z are the x-and z-components of the normal unit.In view of the randomness of the channel walls, Ψ is defined a stochastic stream function, that is, all dynamical variables now become stochastic processes.
Eliminating pressure gradients from Equations ( 15) and ( 16), we find that the stream function satisfies: If there is no roughness, the flow is along x-direction and the stream function would be a function of z only.However, the presence of surface roughness induces a functional variation in the x-direction as well.Therefore, the stream function Ψ can be represented by the slight perturbation expansion of ε, which takes the following form: Substituting Equation (20) into Equation (19), we obtain differential equations for each power of ε: The boundary conditions Equations ( 17) and ( 18) can be expanded in Taylor series about the mean wall positions z = 1, and z = −1.From which the boundary conditions become: In boundary conditions, the argument "±1" denotes the evaluations at z = ±1, and the prime denotes a derivative with respect to the argument.
The boundary conditions of leading order problem are: The leading order ε 0 -solution yields: where A 1 =

Ha
Ha−tanh(Ha) .The boundary conditions of first-order problem are: Based on the expressions of boundary conditions, assuming that the ε 1 -solution is a stationary random function with the following spectral representation: where Φ j (y; λ) (j = 1,2) are functions of y parametrized by λ.Both Φ 1 and Φ 2 hold: where The solutions are: where α j , β j , γ j , δ j (j = 1,2) are constants and given by: where, The mean value of Ψ 1 is zero, due to Equations (7a) and (31).As a result, the mean of stream function depends on the term of order O(ε 2 ).The boundary conditions of secondorder problem are: , Similar to the process of ε 1 -solution representation, the solution to second-order problem can be written as: where ζ jk (z, λ, λ ) are functions of z parameterized by λ and λ .All three functions ζ jk (j, k = 1, 2) satisfy: with the following boundary conditions: Note also that: Although the detailed ε 2 -solution can be calculated, we do not need to know its exact result if one is only interested in the mean quantity.According to Equations (7b-d) and (43), we have: Since the mean of the first-order solutions of the steam function is already zero, in order for the mean of the second-order solutions not being zero, this requires: Therefore, the second-order problem of Equation (44) turns into: The fourth-order ordinary differential equation has the following general solution: where a jk , b jk , c jk , d jk can be determined by boundary conditions.Only the coefficients b jk are needed, and the other coefficients have no contribution to our present problem.The expressions of b jk can be given as: 4(−Ha cosh Ha + sinhHa) 2  (54) Substituting Equation (20) into Equation ( 15), it can be obtained as: Averaging Equation (57) over one wavelength of corrugations, and inserting the relevant results of zero-order and second-order solutions, we have: where: From now, ξ is called as corrugation function used to describe the effect of the random roughness on EMHD flow.It can be seen that from Equation ( 16) over one wavelength of corrugations, the mean pressure gradient in z-direction is zero.Applying the above dimensionless parameters associated with pressure gradient and the definition of S, Equation (58) can be inverted into: The flow rate increases with the corrugations when ξ is negative and decreases with the corrugations when ξ is positive for a certain mean pressure drop and electromagnetic force.When the mean pressure gradient is taken as zero, the flow becomes a pure EMHD flow.
The coefficient of ξ depends on the parameters Ha and λ, which show the impact of corrugations on the driving power required to drive the particular flow rate Q (per unit width of the channel) of flow perpendicular to the corrugations.The driving force is the sum of the mean pressure drop and the electromagnetic force σEB.

Dirichlet's Condition and Applications
It is assumed that m(x) and n(x) satisfy Dirichlet's condition, using harmonic properties of the corrugation function, they can be represented as the Fourier sine series: where λ is wave number, a n and b n are amplitudes, and index n denotes a natural number in the followings for all kinds of corrugations.Under this situation, spectral densities can be expressed as: where δ is the Dirac generalized function.Any shape of the corrugation represented as the Fourier series can be solved in this way.Previous studies have shown that the phase difference between two corrugated walls has significant effects on the flow.This makes us consider two cases where the upper and lower corrugated walls are either in phase or half-period out of phase.Here, 'in phase' corresponds to the exact same direction of the corrugations on the upper and lower walls, while 'half-period out of phase' means the opposite direction of the corrugations on the higher and lower walls.As particular examples of the above study, we consider four types of corrugations: sinusoidal, triangular, rectangular, and sawtooth corrugations.

Sinusoidal Corrugations
Under the sinusoidal corrugations, m(x) and n(x) can be represented as follows: with n = 1 and a n = b n = 1 in Equation (61).In this case, the spectral densities and corresponding corrugation function can be written as: Interestingly, by substituting the relevant coefficients α 1 , β 1 , γ 1 , δ 1 , k 1 , and k 2 in Equations ( 36)-(40) into Equation ( 65), we find that the final analytical expression of the corrugation function ξ 1 is the exactly same as the results obtained by Buren and Jian [34].

Triangular Corrugations
The EMHD flow diagram of the triangular corrugations are shown in Figure 2.

Triangular Corrugations
The EMHD flow diagram of the triangular corrugations are shown in Figure 2.For the triangular corrugations, m(x) and n(x) can be represented as: The Fourier coefficient an and bn can be determined as: Corresponding spectral densities and the corrugation function ξ2 are: For the triangular corrugations, m(x) and n(x) can be represented as: The Fourier coefficient a n and b n can be determined as: Corresponding spectral densities and the corrugation function ξ 2 are: where, (69b)

Rectangular Corrugations
The EMHD flow diagram of the rectangular corrugations are shown in Figure 3.
Micromachines 2023, 14, x FOR PEER REVIEW 13 o The EMHD flow diagram of the rectangular corrugations are shown in Figure 3.In case of the rectangular corrugations, m(x) and n(x) can be represented as: The equations of the upper and lower walls are then given, respectively: Corresponding spectral densities and the corrugation function ξ3 are:  In case of the rectangular corrugations, m(x) and n(x) can be represented as: The equations of the upper and lower walls are then given, respectively: where a n = 4 π(2n−1) .Corresponding spectral densities and the corrugation function ξ 3 are: where,

Sawtooth Corrugations
The EMHD flow diagram of the sawtooth corrugations are shown in Figure 4.
Micromachines 2023, 14, x FOR PEER REVIEW 14 o The EMHD flow diagram of the sawtooth corrugations are shown in Figure 4.For the sawtooth corrugations, m(x) and n(x) can be represented as: The equations of the upper and lower walls are then given respectively as: Corresponding spectral densities and the corrugation function ξ4 are: ) where,

Results and Discussion
Taking the same parameter values, through comparison, we found that our pres result is in exact agreement with the findings of Buren and Jian [34].However, compa with those [34], we extended the shape of the wall from sine to random shape, and calculation process was simplified.The influence of Ha on the corrugation function ξ For the sawtooth corrugations, m(x) and n(x) can be represented as: The equations of the upper and lower walls are then given respectively as: where a n = ± 2 π(2n−1) .Corresponding spectral densities and the corrugation function ξ 4 are: where,

Results and Discussion
Taking the same parameter values, through comparison, we found that our present result is in exact agreement with the findings of Buren and Jian [34].However, compared with those [34], we extended the shape of the wall from sine to random shape, and the calculation process was simplified.The influence of Ha on the corrugation function ξ for different values of wave number λ is illustrated in Figure 5, when the corrugations are in the phase θ = 0. Increases in Ha always reduce the impact of corrugations on the EMHD flow for given λ.When Ha is fixed, it was found that the larger λ leads to larger roughness function, which means more driving force is needed to actuate the flow.Within the low range 0 < Ha < 2, the corrugation function declines sharply, and then becomes gentle gradually for a higher Ha number.In addition, it can be seen from Figure 5 that, when wave number is small (such as λ = 0.4), the wall roughness has almost no effect on the flow.function, which means more driving force is needed to actuate the flow.Within the low range 0 < Ha < 2, the corrugation function declines sharply, and then becomes gentl gradually for a higher Ha number.In addition, it can be seen from Figure 5 that, when wave number is small (such as λ = 0.4), the wall roughness has almost no effect on th flow.In order to generalize the result of Buren and Jian [34], we used the same parametrica values to study the triangular shape of the wall roughness in Figure 6. Figure 6a show that, when λ is fixed, the roughness function decreases with the Ha number.When 0.1 < λ < 0.2, the roughness function has a little discrepancy for different Ha.However, when λ 0.2, ξ increases with λ, and a lower value of Ha results in larger roughness function.Figur 6b shows that, with the increase of Ha, the roughness function ξ constantly decreases and gradually attains a stable status.The reason for this is that the electromagnetic force can suppress the effect o roughness and reduce the loss caused by the resistance, which is related to the second order problem.
Figures 7 and 8 show the variations of the corrugation function ξ with Ha and λ fo rectangular and sawtooth shape roughness, respectively.the corrugation function i always positive, which shows that the flow rate always reduces with the increase o roughness when the average pressure drop and electromagnetic fields are prescribed.In order to generalize the result of Buren and Jian [34], we used the same parametrical values to study the triangular shape of the wall roughness in Figure 6. Figure 6a shows that, when λ is fixed, the roughness function decreases with the Ha number.When 0.1 < λ < 0.2, the roughness function has a little discrepancy for different Ha.However, when λ > 0.2, ξ increases with λ, and a lower value of Ha results in larger roughness function.Figure 6b shows that, with the increase of Ha, the roughness function ξ constantly decreases and gradually attains a stable status.
different values of wave number λ is illustrated in Figure 5, when the corrugations are in the phase θ = 0. Increases in Ha always reduce the impact of corrugations on the EMHD flow for given λ.When Ha is fixed, it was found that the larger λ leads to larger roughness function, which means more driving force is needed to actuate the flow.Within the low range 0 < Ha < 2, the corrugation function declines sharply, and then becomes gentle gradually for a higher Ha number.In addition, it can be seen from Figure 5 that, when wave number is small (such as λ = 0.4), the wall roughness has almost no effect on the flow.In order to generalize the result of Buren and Jian [34], we used the same parametrical values to study the triangular shape of the wall roughness in Figure 6. Figure 6a shows that, when λ is fixed, the roughness function decreases with the Ha number.When 0.1 < λ < 0.2, the roughness function has a little discrepancy for different Ha.However, when λ > 0.2, ξ increases with λ, and a lower value of Ha results in larger roughness function.Figure 6b shows that, with the increase of Ha, the roughness function ξ constantly decreases and gradually attains a stable status.The reason for this is that the electromagnetic force can suppress the effect of roughness and reduce the loss caused by the resistance, which is related to the secondorder problem.
Figures 7 and 8 show the variations of the corrugation function ξ with Ha and λ for rectangular and sawtooth shape roughness, respectively.the corrugation function is always positive, which shows that the flow rate always reduces with the increase of roughness when the average pressure drop and electromagnetic fields are prescribed.The reason for this is that the electromagnetic force can suppress the effect of roughness and reduce the loss caused by the resistance, which is related to the second-order problem.
Figures 7 and 8 show the variations of the corrugation function ξ with Ha and λ for rectangular and sawtooth shape roughness, respectively.the corrugation function is always positive, which shows that the flow rate always reduces with the increase of roughness when the average pressure drop and electromagnetic fields are prescribed.
Similarly, we can study the effect of phase difference between upper and lower plates on EMHD flow.For our present four kinds of roughness, only sinusoidal corrugation strictly has a conception of phase difference, and this problem has been studied by Buren and Jian [34].For other kinds of roughness, we only consider two special cases of θ = 0 (in phase) and θ = π (out of phase).
The variations of triangular, rectangular, and sawtooth corrugation function with Ha and λ for phase difference θ = 0 (in phase) and θ = π (out of phase) are illustrated in Figures 9-11.However, in the present derivation, we found that only the Fourier coefficients for different shapes of ripples are different, and these coefficients determine the magnitude of the corrugation function.However, the physical meaning behind Figures 10 and 11 is the same as Figure 9, which reveals that the corrugation function in phase of two plates is smaller than that when the roughness is out of phase.Furthermore, it means that the flow resistance becomes small, and the flow rate increases for in-phase roughness compared with the one out of phase.Similarly, we can study the effect of phase difference between upper and lowe on EMHD flow.For our present four kinds of roughness, only sinusoidal corr strictly has a conception of phase difference, and this problem has been studied b and Jian [34].For other kinds of roughness, we only consider two special cases of phase) and θ = π (out of phase).
The variations of triangular, rectangular, and sawtooth corrugation function and λ for phase difference θ = 0 (in phase) and θ = π (out of phase) are illustrated in 9-11.However, in the present derivation, we found that only the Fourier coeffici different shapes of ripples are different, and these coefficients determine the ma of the corrugation function.However, the physical meaning behind Figures 10 a the same as Figure 9, which reveals that the corrugation function in phase of two smaller than that when the roughness is out of phase.Furthermore, it means that t resistance becomes small, and the flow rate increases for in-phase roughness co with the one out of phase.Similarly, we can study the effect of phase difference between upper and low on EMHD flow.For our present four kinds of roughness, only sinusoidal cor strictly has a conception of phase difference, and this problem has been studied and Jian [34].For other kinds of roughness, we only consider two special cases of phase) and θ = π (out of phase).
The variations of triangular, rectangular, and sawtooth corrugation function and λ for phase difference θ = 0 (in phase) and θ = π (out of phase) are illustrated in 9-11.However, in the present derivation, we found that only the Fourier coeffic different shapes of ripples are different, and these coefficients determine the m of the corrugation function.However, the physical meaning behind Figures 10 the same as Figure 9, which reveals that the corrugation function in phase of two smaller than that when the roughness is out of phase.Furthermore, it means that resistance becomes small, and the flow rate increases for in-phase roughness c with the one out of phase.Form Figure 12, we can note that, although the three corrugations have the same trend, the roughness of rectangular is significantly greater than that of triangular roughness and sawtooth roughness.The reason for this difference, mathematically speaking, is that their Fourier coefficients are different.This inspired us to make the microchannel rough shape triangular when we want to increase flow, and rectangular when we want to promote mixing.Form Figure 12, we can note that, although the three corrugations have the same trend, the roughness of rectangular is significantly greater than that of triangular roughness and sawtooth roughness.The reason for this difference, mathematically speaking, is that their Fourier coefficients are different.This inspired us to make the microchannel rough shape triangular when we want to increase flow, and rectangular when we want to promote mixing.Form Figure 12, we can note that, although the three corrugations have the same trend, the roughness of rectangular is significantly greater than that of triangular roughness and sawtooth roughness.The reason for this difference, mathematically speaking, is that their Fourier coefficients are different.This inspired us to make the microchannel rough shape triangular when we want to increase flow, and rectangular when we want to promote mixing.

Conclusions
We have investigated the impact of small random transverse wall roughness on EMHD flow in a microchannel based on random function theory.The stream function is expanded with small parameters (amplitude of the random roughness), and some properties of stationary random process are employed to gain analytical expression of the random corrugation function ξ, which is a measure of flow rate deviated from the case with no roughness on the two plates.Using the method presented in this paper, we can study the corrugation function of random shapes.When the ripple is sinusoidal, we find the final expression of the corrugation function is the same as that obtained by Buren and Jian [34].No matter what kind of corrugation, identical conclusions can be drawn that increasing wavenumber λ and decreasing Ha number result in a growth of the corrugation

Conclusions
We have investigated the impact of small random transverse wall roughness on EMHD flow in a microchannel based on random function theory.The stream function is expanded with small parameters (amplitude of the random roughness), and some properties of stationary random process are employed to gain analytical expression of the random corrugation function ξ, which is a measure of flow rate deviated from the case with no roughness on the two plates.Using the method presented in this paper, we can study the corrugation function of random shapes.When the ripple is sinusoidal, we find the final expression of the corrugation function is the same as that obtained by Buren and Jian [34].No matter what kind of corrugation, identical conclusions can be drawn that increasing wavenumber λ and decreasing Ha number result in a growth of the corrugation function ξ.These results mean that the electromagnetic force can suppress the effects of roughness, increase flow, and promote production.In addition, we studied the influences of phase differences θ = 0 and θ = π on the corrugation function for all kinds of roughness.We found that the flow resistance becomes small and the flow rate increases for in-phase roughness compared with the one out of phase.In the future, we can use this method to study electroosmosis flow in microchannel with random roughness.

Figure 1 .
Figure 1.Schematic illustration of EMHD flow through a microchannel with a random roughness wall.

Figure 1 .
Figure 1.Schematic illustration of EMHD flow through a microchannel with a random roughness wall.

Figure 1 .
Figure 1.Schematic illustration of EMHD flow through a microchannel with a random roughness wall.

Figure 1 .
Figure 1.Schematic illustration of EMHD flow through a microchannel with a random roughness wall.

Figure 1 .
Figure 1.Schematic illustration of EMHD flow through a microchannel with a random roughness wall.

Figure 1 .
Figure 1.Schematic illustration of EMHD flow through a microchannel with a random ro wall.

Figure 1 .
Figure 1.Schematic illustration of EMHD flow through a microchannel with a random roughness wall.

Figure 6 .
Figure 6.The variations of triangular corrugation function ξ with (a) for different Ha, (b) Ha fo different λ.

Figure 6 .
Figure 6.The variations of triangular corrugation function ξ with (a) λ for different Ha, (b) Ha for different λ.

Figure 6 .
Figure 6.The variations of triangular corrugation function ξ with (a) λ for different Ha, (b) Ha for different λ.

Figure 8 .
Figure 8.The variations of sawtooth corrugation function ξ with (a) λ for different Ha, (b different λ.

Figure 8 .
Figure 8.The variations of sawtooth corrugation function ξ with (a) λ for different Ha, different λ.

Figure 8 .
Figure 8.The variations of sawtooth corrugation function ξ with (a) λ for different Ha, (b) Ha for different λ.Micromachines 2023, 14, x FOR PEER REVIEW

Author
Contributions: N.M.: investigation, methodology, validation, formal analysis, writingoriginal draft, writing-review and editing.Y.S.: visualization, investigation, methodology.Y.J.: supervision, writing-review and editing.All authors have read and agreed to the published version of the manuscript.Funding: The authors acknowledge the financial support provided by the National Natural Science Foundation of China (Grant No. 12262026), the Natural Science Foundation of Inner Mongolia Autonomous Region of China (Grant No. 2021MS01007), the Inner Mongolia Grassland Talent (No. 12000-12102013), and the Research Program of science and technology at Universities of Inner Mongolia Autonomous Region (Grant No. NJZY23054).Data Availability Statement: Data sharing not applicable.Conflicts of Interest:The authors declare no conflict of interest.NomenclatureThe bold terms represent vectors and the others represent scalars.(x, y, z) corrugation m(x), n(x) stationary random functionsdZ i (λ)(i = 1, 2)interval random functionδ ij (j, k = 1, 2) Kronecker delta F ij (j, k = 1, 2) real non-decreasing bounded functions |R ij (s)|(j,k = 1, 2) correlation function Ω jk (j, k = 1, 2)spectral densities of m(x) and n(x).Ψ stream function Q flow rate per unit width Ha Hartmann number S dimensionless electric field strength n normal vector n x , n z the x-and z-components of the unit normal A 1 coefficient of the zero-order solution α j , β j , γ j , δ j (j = 1, 2) coefficients of first-order solutions a jk , b jk , c jk , d jk (j, k = 1, 2) coefficients of second-order solutions