Next Article in Journal
Improving the Energy Efficiency of Vehicles by Ensuring the Optimal Value of Excess Pressure in the Cabin Depending on the Travel Speed
Next Article in Special Issue
Shedding of Cavitation Clouds in an Orifice Nozzle
Previous Article in Journal
Measuring Turbulent Flows: Analyzing a Stochastic Process with Stochastic Tools
Previous Article in Special Issue
Prediction of Critical Heat Flux during Downflow in Fully Heated Vertical Channels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hartmann Flow of Two-Layered Fluids in Horizontal and Inclined Channels

School of Mechanical Engineering, Faculty of Engineering, Tel-Aviv University, Tel-Aviv 6997801, Israel
*
Author to whom correspondence should be addressed.
Fluids 2024, 9(6), 129; https://doi.org/10.3390/fluids9060129
Submission received: 1 May 2024 / Revised: 21 May 2024 / Accepted: 27 May 2024 / Published: 30 May 2024
(This article belongs to the Special Issue Numerical Modeling and Experimental Studies of Two-Phase Flows)

Abstract

:
The effect of a transverse magnetic field on two-phase stratified flow in horizontal and inclined channels is studied. The lower heavier phase is assumed to be an electrical conductor (e.g., liquid metal), while the upper lighter phase is fully dielectric (e.g., gas). The flow is defined by prescribed flow rates in each phase, so the unknown frictional pressure gradient and location of the interface separating the phases (holdup) are found as part of the whole solution. It is shown that the solution of such a two-phase Hartmann flow is determined by four dimensionless parameters: the phases’ viscosity and flow-rate ratios, the inclination parameter, and the Hartmann number. The changes in velocity profiles, holdups, and pressure gradients with variations in the magnetic field and the phases’ flow-rate ratio are reported. The potential lubrication effect of the gas layer and pumping power reduction are found to be limited to low magnetic field strength. The effect of the magnetic field strength on the possibility of obtaining countercurrent flow and multiple flow states in concurrent upward and downward flows, and the associated flow characteristics, such as velocity profiles, back-flow phenomena, and pressure gradient, are explored. It is shown that increasing the magnetic field strength reduces the flow-rate range for which multiple solutions are obtained in concurrent flows and the flow-rate range where countercurrent flow is feasible.

1. Introduction

Flows of electrically conducting fluids in an electromagnetic field are found across a wide variety of phenomena, which are extensively investigated in the field of magnetohydrodynamics (MHD) and possess important technological applications (e.g., [1,2]). While much of the existing research on MHD focuses on single-phase flows, numerous practical challenges involve multiphase flow systems, particularly within the nuclear and petroleum industries, geophysics, and MHD power generation [3,4]. Among the applications are two-phase liquid metal magnetohydrodynamics (LMMHD) generators, which attracted attention in the industry due to their relatively simple structure, absence of moving parts, high-efficiency conversion, and reduced environmental pollution [5]. Alongside the liquid metal, another fluid phase (gas or liquid) is employed to convert thermal energy into kinematic energy. MHD flows can also transport weakly conducting fluids in microscale systems, e.g., in the microchannel networks, such as the microchannel networks of lab-on-a-chip devices [6,7], where the presence of a second non-conductive fluid enhances the mobility of the conducting fluid. Magnetic field-driven micropumps are in increasing demand due to their long-term reliability in generating flow, absence of moving parts, low power requirements, flow reversibility, and efficient mixing [8,9]. In all those applications, the system performance is dependent on the gas–liquid flow pattern in the channel, and the stability of the interface between the phases, which have been subject to experimental and numerical investigations (e.g., [10,11,12,13,14]).
Since the pioneering works of Hartmann and Lazarus in 1937 [15] the Poiseuille flow of an electrically conducting fluid in a transverse magnetic field (i.e., Hartmann flow) has been thoroughly studied in the literature (e.g., [16,17,18,19]). Hartmann [20] presented an analytical solution for the velocity profile of an isothermal laminar flow in an electrically insulated channel subjected to a wall-normal uniform magnetic field. It was demonstrated that the magnetic field flattens the velocity profile, a phenomenon often termed the ‘Hartmann effect’ in the literature. However, the presence of a second phase can drastically affect the velocity profile and the associated pressure gradient.
Due to the density difference between the phases, a stratified flow configuration is commonly encountered in both horizontal and inclined conduits. The part of the flow cross-section area occupied by the heavier fluid is referred to as the holdup, being an additional characteristic of the flow pattern. Shail [21] investigated the Hartmann flow of a conducting fluid between two horizontal insulating plates, with a layer of non-conducting fluid separating the top wall from the conducting fluid. It was discovered that the flow rate of the conducting fluid can be increased by approximately 30 percent for suitable ratios of the depths and viscosities of the two fluids used in an electromagnetic pump. Owen et al. [22] introduced two mechanistic models: a simple homogeneous model and a more sophisticated film flow model for computing the pressure drop in an MHD two-phase flow at high Hartmann numbers. Lohrasbi and Sahai [23] derived analytical solutions for the velocity and temperature profiles in a two-phase laminar steady MHD flow between two horizontal plates with one conducting phase. Malashetty and Umavathi [24], analytically investigated a similar system in an inclined channel, solving the nonlinear coupled momentum and energy equations in both phases. They found that, for a constant thickness of the layers, increasing the magnetic field has a dampening effect on the velocity of the conducting layer, akin to a single-phase flow. They also noted an increase in the velocity of the conducting fluid with a decrease in the thickness of the upper electrically non-conducting layer and/or with an increase in the channel inclination angle. A subsequent study by Umavathi et al. [25] addressed the magnetohydrodynamic Poiseuille–Couette flow and heat transfer of two immiscible fluids between inclined parallel plates. However, for isothermal flow, the effect of the different gravity driving force in the two layers is not considered in the model equations. Recently, Shah et al. [26] derived approximate analytical solutions for the velocity and temperature fields of unsteady MHD generalized Couette flows of two immiscible and electrically conducting fluids flowing between two horizontal electrically insulated plates subjected to an inclined magnetic field and an axial electric field. In all of these studies, the problem is solved for a pre-defined flow configuration (i.e., in situ holdup of the heavy layer and the pressure gradient).
An important issue in the context of two-phase flows is the occurrence of gravity-driven multiple solutions in inclined channels for a specified two-phase system and fixed operational conditions. It has been established in the two-phase flow literature [27,28,29,30] that there always exist two possible solutions for the holdup of steady countercurrent flow and up to three distinct steady-state solutions can be obtained within a limited range of the flow parameters in concurrent upward- and downward-inclined flows. The existence of multiple solutions for the holdup for some range of the superficial velocities in liquid–liquid flows was verified experimentally [28,29], where experimental verification of the existence of multiple solutions for the holdup within a certain range of superficial velocities in liquid–liquid flows has been conducted. However, the feasibility of achieving countercurrent flow in the presence of a magnetic field, and its impact on the potential for multiple stratified flow configurations in concurrent and countercurrent scenarios has not yet been explored.
In this study, we investigate the Hartmann flow of a conducting fluid between two parallel insulating plates, with a layer of non-conducting fluid separating the top wall from the conducting fluid. The presented analytical solution allows for the first time determination of the in situ holdup and the pressure gradient in horizontal or inclined channels for a specified two-phase system and fixed operational conditions. We identify the corresponding dimensionless input parameters that dictate the in situ stratified flow configuration. Furthermore, we utilize the solution to examine the impact of the presence of a second non-conductive layer and the intensity of the magnetic field on the flow characteristics in concurrent and countercurrent flows of the fluids. Special attention is given to the potential lubrication and pumping power saving achieved by introducing the gas phase, as well as to the operational conditions for which multiple solutions can be obtained for the stratified flow configuration in inclined flows.

2. Two-Phase Horizontal Flow

We consider isothermal, axial, steady and fully developed laminar stratified flow of two immiscible fluids between two infinite electrically insulating plates positioned at z = 0 and z = H (see Figure 1). Under these conditions, the pressure gradient, G = d p / d x , is the same in both layers and is constant (e.g., [31]). For horizontal flow, the channel inclination to the horizontal is β = 0. A plane and smooth interface between the phases is located at z = h . The lower layer (1) is occupied by an electrically conducting fluid (e.g., mercury), which is affected by a constant transverse magnetic field B = B 0 e z , while the upper lighter fluid is assumed to be non-conductive (e.g., air). Then, the x-component of the momentum equations in the two layers are as follows (e.g., Shail [21]):
η 1 d 2 u 1 ( z ) d z 2 σ 1 B 0 2 u 1 z = G 0 z h
  η 2 d 2 u 2 ( z ) d z 2 = G h z H
where η 1 , η 2 are the fluids’ viscosities and σ 1 is the heavier fluid’s electric conductivity. The solutions for the velocity profiles are as follows:
u 1 z ~ = A 1 cosh H a   z ~ + A 2 sinh H a   z ~ G σ 1 B 0 2 0 z ~ h ~
u 2 z ~ = G H 2   2 η 2 z ~ 2 + B 1 H z ~ + B 2 h ~ z ~ 1
where z ~ = z / H , h ~ = h / H is the conductive (heavier) fluid holdup and Ha is the Hartmann number H a = B 0 H σ 1 η 1 . Applying no-slip boundary conditions at z ~ = 0 and z ~ = 1 , and the continuity of velocities and viscous shear stresses at the fluids’ interface z ~ = h ~ yields the following:
A 1 = G σ 1 B 0 2 ,           B 2 = G H 2   2 η 2 B 1 H
A 2 = G H 2   η 2 1 h ~ 2 1 H a 2 η 12 1 h ~ cosh H a h ~ 1 1 H a sinh H a h ~ H a η 12 cosh H a h ~ + 1 1 h ~ sinh H a h ~
B 1 = G H η 2 h ~ + 1 H a s i n h H a h ~ + H a η 1 A 2 G H 2 cosh H a h ~
The volumetric flow rates (per unit channel width) in the lower and upper layers are as follows:
Q 1 = H 0 h ~ u 1 z ~ d z ~ = G H 3 H a 3 η 1 sinh H a h ~ + A 2 H H a cosh H a h ~ 1 G H 3 h ~ H a 2 η 1
Q 2 = H h ~ 1 u 2 z ~ d z ~ = G   H 3   6 η 2 2 h ~ 3 + 3 h ~ + B 1 H 2 2 1 h ~ 2 + 2 h ~
As the left hand side (l.h.s) of both Equations (8) and (9) are proportional to G, the ratio of Q1/Q2 yields an implicit equation of calculation of the holdup of the conductive layer in the channel, h ~ = h / H :
Q 12 = Q 1 Q 2 = U 1 s U 2 s = 1 H a 3 1 η 12 sinh H a h ~ + H a 2 A 2   η 2 G H 2 cosh H a h ~ 1 H a h ~ η 12 1 6 2 h ~ 3 + 3 h ~ B 1   η 2 2 G H 1 h ~ 2
where U 1,2 s = Q 1,2 / H is the superficial velocity of the fluid in the channel cross-section. In view of Equations (6) and (7), the following expressions are introduced:
A 2   η 2 G H 2 = 1 h ~ 2 1 H a 2 η 12 1 h ~ cosh H a h ~ 1 1 H a sinh H a h ~ H a η 12 cosh H a h ~ + 1 1 h ~ sinh H a h ~
B 1   η 2 G H = h ~ + 1 H a s i n h H a h ~ + H a η 1 A 2 G H 2 cosh H a h ~
and A 2 η 1 G H 2 = η 12 A 2 η 2 G H 2 .
The calculation of the flow-rate ratio for a specified holdup via Equation (10) (with (11) and (12)) is straightforward. However, in practice, the input flow rates of the fluids are known, whereas the location of the interface is unknown. Equation (10) indicates that h ~ = h ~ Q 12 , η 12 , H a . In the limit of Ha → 0, Equation (10) converge to the expression obtained for Q 12 h ,   ~ η 12 for a two-layer Poiseuille flow in the absence of a magnetic field (e.g., [32,33]).
Once Equation (10) is solved for the holdup, the corresponding pressure gradient G can be determined. For example, using the solution for the holdup in Equation (9) yields the following:
G = d p d x = d p f d x = U 2 s   H 2   6 η 2 2 h ~ 3 + 3 h ~ B 1 H 2 G 1 h ~ 2
where B1/G is independent of G (see Equation (7)). Obviously, in horizontal flow, the total pressure gradient is identical to the frictional pressure gradient. The corresponding dimensionless pressure gradient,   P f 2,0   P f 1,0 , can be obtained by normalizing d p / d x with respect to the frictional pressure gradient of the single-phase flow of either the lighter or the heavier fluid obtained in the absence of a magnetic field, ( d p f / d x ) 2 s , 1 s 0 = 12 η 2,1 U 2 s , 1 s /H2, respectively, whereby
P f 2,0 = d p / d x ( d p f / d x ) 2 s 0 = 1 2 1 h ~ 3 3 h ~ + 2 + 3 η 2 B 1 G H ( 1   h ~ ) 2
P f 1,0 = d p / d x ( d p f / d x ) 1 s 0 = P f 2,0 Q 12   η 12
where the first superscript (1 or 2) represents the phase selected to scale the two-phase pressure gradient, and the second superscript (0) indicates that the pressure gradient in the absence of a magnetic field is used for the scaling (see also the nomenclature list). Note that B 1   η 2 G H is also determined by ( Q 12 , η 12 , H a ) (see Equation (7)).
The dimensionless pressure gradient that would be obtained in the case of a single-phase (SP) flow of the conductive (heavier) fluid under the same magnetic field strength (i.e., the same Ha) is as follows:
P f s 1,0 = d p / d x 1 s ( d p f / d x ) 1 s 0 = H a 2 12 1 2 H a t g h ( H a 2 )
Hence, when scaling the two-phase pressure gradient with respect to the SP flow of the heavy phase in the presence of the same magnetic field, the dimensionless pressure gradient is given by the following:
P f 1 = P f 1,0 P f s 1 = 12 G ~ 1,0 H a 2   1 2 H a t g h ( H a 2 )
The corresponding dimensionless pumping power (power factor) required for the given flow rates is as follows:
P o = P f 1 ( U 1 s + U 2 s ) U 1 s = P f 1 ( 1 + Q 21 )
Hence, the dimensionless pressure gradient (either P f 2,0 ,   P f 1,0 or P f 1 ) and the power factor are also determined in terms of ( Q 12 , η 12 , H a ) .
Finally, upon scaling the velocity with respect to U 1 s , the dimensionless velocity profiles are also obtained in terms of those three dimensionless parameters:
u ~ 1 z ~ = u 1 z U 1 s = A ~ 1 cosh H a z ~ 1 + A ~ 2 sinh H a z ~
u ~ 2 z ~ = u 2 z U 1 s = B ~ 2 z ~ 2 1 + B ~ 1 z ~ 1
with
A ~ 1 = 12 P f 1,0 H a 2 ;     B ~ 2 = 6 η 12 P f 1,0
A ~ 2 = 12 η 12 P f 1,0 1 h ~ 2 1 H a 2 η 12 1 h ~ cosh H a h ~ 1 1 H a sinh H a h ~ H a η 12 cosh H a h ~ + 1 1 h ~ sinh H a h ~
B ~ 1 = 12 η 12 P f 1,0 h ~ + 1 H a s i n h H a h ~ + H a η 1 A 2 G H 2 cosh H a h ~

3. Two-Phase Inclined Flow

In the case of inclined flow, the momentum equations in the two layers are
η 1 d 2 u 1 ( z ) d z 2 σ 1 B 0 2 u 1 z = d p d x ρ 1 g s i n β = G 1 0 z h
η 2 d 2 u 2 ( z ) d z 2 = d p d x ρ 2 g s i n β = G 2 h z H
where β is the angle of the channel downward inclination to the horizontal (x is in the downward direction) and ρ 1 , ρ 2 are the fluids’ densities. The model assumptions are the same as those indicated in the case of the horizontal channel (see Section 2). The terms describing the gravity body force, which drives the flow along with the imposed pressure gradient, are added to the equations. The solution of (24) and (25) for the velocity profile that satisfies the no-slip boundary conditions at z ~ = 0 and z ~ = 1 is as follows:
u 1 z ~ = G 1 σ 1 B 0 2 cosh H a   z ~ 1 + A 2 sinh H a   z ~ 0 z ~ h ~
u 2 z ~ = G 2   2 η 2 z ~ 2 1 + B 1 z ~ 1 h ~ z 1
Imposing the boundary conditions of velocities and shear stress continuity at the interface yields the following:
B 1 H = G 2 H 2 h ~ η 2 + G 1 H 2 η 2 1 H a sinh H a h ~ + A 2 Ha η 12 cosh H a h ~
A 2 = G 2 H 2   2 η 2 1 h ~ G 1 H 2 H a 2 η 1 1 h ~ cosh H a h ~ 1 G 1 H 2 H a η 2 sinh H a h ~ η 12 Ha cosh H a h ~ + 1 1 h ~ sinh H a h ~
The volumetric fluxes in the lower and upper layers are as follows:
Q 1 = U 1 s H = H 0 h ~ u 1 z ~ d z ~ = G 1 H 3 H a 3 η 1 sinh H a h ~ + A 2 H H a cosh H a h ~ 1 G 1 H 3 h ~ H a 2 η 1
Q 2 = U 2 s H = H h ~ 1 u 2 z ~ d z ~ = 1 h ~ 2 G 2   H 3   6 η 2 h ~ + 2 + B 1 H 2 2
Substituting B 1 H (Equation (28)) in Equation (31) and dividing Equations (30) and (31) by U 1 s yields the following:
1 = 12 G ~ 1 1,0 H a 3 sinh H a h ~ + A 2 ~ H a cosh H a h ~ 1 + 12 G ~ 1 1,0 h ~ H a 2 ,
U 2 s U 1 s = 2 1 h ~ 2 { G ~ 2 1,0 η 12 h ~ + 2 + [ 3 G ~ 2 1,0 η 12 h ~ + 3 G ~ 1 1,0 η 12 H a sinh H a h ~ A 2 ~ H a η 12 4 cosh H a h ~ ] }
where the following relations were introduced:
G 2 U 1 s η 2 / H 2 = 12 η 12 G 2 ( d p f / d x ) 1 s 0 = 12 η 12 G ~ 2 1,0
G 1 U 1 s η 1 / H 2 = 12 G 1 ( d p f / d x ) 1 s 0 = 12 G ~ 1 1,0
and A ~ 2 = A 2 / U 1 s , which by Equation (29) reads as follows:
A ~ 2 = 12 G ~ 2 1,0 η 12 1 h ~ 2 + G ~ 1 1,0 H a 2 1 h ~ cosh H a h ~ 1 + G ~ 1 1,0 η 12 H a sinh H a h ~ H a η 12 cosh H a h ~ + 1 1 h ~ sinh H a h ~
In principle, once U 1 s ,   U 2 s are known, Equations (32) and (33) can be solved for the unknown holdup,   h ~ , and pressure gradient, d p / d x , embedded both in G ~ 1 1,0 and G ~ 2 1,0 . However, to reveal the dimensionless parameters of the solution in inclined flows, the solution form is further manipulated.
The total pressure gradient is a sum of the frictional and gravitational (hydrostatic) pressure gradients; hence:
d p d x = d p d x f + d p d x g = d p d x f + ρ 1 h ~ + ρ 2 ( 1 h ~ ) g s i n β
Therefore:
G 2 = d p d x ρ 2 g s i n β = d p d x f +   h ~ ( ρ 1 ρ 2 ) g s i n β
and
G ~ 2 1,0 = G 2 ( d p f / d x ) 1 s 0 = P ~ f 1,0 +   h ~ Y 1,0   ;   Y 1,0 = ( ρ 1 ρ 2 ) g s i n β ( d p f / d x ) 1 s 0
where P f 1,0 is the dimensionless frictional pressure gradient in the inclined two-phase flow and Y 1,0 is the (a priori known) inclination parameter. Similarly,
G 1 = d p d x ρ 1 g s i n β = d p d z f 1 h ~ ( ρ 1 ρ 2 ) g s i n β
and
G ~ 1 1,0 = G 1 ( d p f / d x ) 1 s 0 = P ~ f 1,0 1 h ~ Y 1,0
Upon substituting (39) and (41) in Equations (32) and (33) (and in (36)), a solution for the holdup h ~ and the dimensionless frictional gradient P f 1,0 can be obtained in terms of the four (a priori known) dimensionless parameters ( η 12 , Q 12 , H a ,   Y 1,0 ). The values of those parameters are a priori known for a specified two-flow system (i.e., channel geometry and inclination, fluid properties, and magnetic field intensity).
In fact, Equation (32) can be used to obtain an explicit expression for P f 1,0 in terms of ( η 12 , H a , Y 1,0 , h ~ ) , which can then be substituted into Equation (33) to obtain an implicit equation for the holdup. This step is important to explore the possibility of multiple solutions for prescribed values of the dimensionless parameters, since this becomes rather straightforward once a single algebraic equation for the holdup is obtained.
In view of Equation (37), the total pressure gradient in inclined flow (normalized with respect to the frictional pressure gradient of mercury flow in the absence of a magnetic field) is given by the following:
P 1,0 = d p / d x ( d p f / d x ) 1 s 0 = P f 1,0 + P g 1,0 = P f 1,0 + h ~ Y 1,0 + ρ 2 ρ Y 1,0
where the sum of the last two terms on the right hand side (r.h.s) represents the dimensionless hydrostatic pressure gradient, P g 1,0 . When normalized with respect to the frictional pressure gradient obtained for SP flow of the conductive layer under the same magnetic field,   d p f / d x 1 s (see Equation (16)), the dimensionless total pressure gradient is then P 1 = P 1,0 / P f s 1,0 ,     ( P f 1 and P g 1 are the frictional and hydrostatic components). Note that the density ratio,   ρ 12 = ρ 1 / ρ 2 , is an additional dimensionless parameter that has to be prescribed in order to obtain the total pressure and power factors. However, in the case the lighter (non-conductive) fluid is a gas, the last term on the r.h.s of Equation (42) is negligible and practically P 1,0 G ~ 2 1,0 .
To examine the effect of the gas layer on the pressure gradient, its value compared to the total pressure gradient in SP flow of the conductive fluid which develops under the same Ha is of interest. It is given by the following:
P 1 T = P 1 1 + ρ 2 ρ Y 1,0 P f s 1,0
and the corresponding power factor is P o = P 1 T ( 1 + Q 21 ) . Once the holdup and pressure gradient have been obtained, the corresponding dimensionless velocity profile in the two layers are given by the following:
u ~ 1 z ~ = u 1 z ~ U 1 s = A ~ 1 cosh H a z ~ 1 + A ~ 2 sinh H a z ~
u ~ 2 z ~ = u 2 z ~ U 1 s = B ~ 2 z ~ 2 1 + B ~ 1 z ~ 1
The dimensionless coefficients are given by the following:
A ~ 1 = 12 G ~ 1 1,0 H a 2 ;   B ~ 2 = G 2 H 2   2 η 2 U 1 s = 6 η 12 G ~ 2 1,0
A ~ 2 is given in Equation (36) and by using Equation (28):
B ~ 1 = 12 G ~ 2 1,0 η 12 h ~ G ~ 1 1,0 η 12 H a s i n h H a h ~ + A ~ 2 H a η 12 cosh H a h ~

4. Results and Discussion

The analytical solution presented in Section 2 and Section 3 enables determining the conductive liquid holdup and the dimensionless pressure gradient in MHD two-phase flow in terms of the set of dimensionless ( η 12 , Q 21 , H a , Y 1,0 ) parameters. The values of these parameters are determined by the channel geometry, fluid properties, and flow rates, and are independent of the in situ flow characteristics. Worth noting is that using the above analytical solution for the calculation of the holdups and velocity profiles appears to be non-trivial. Due to the hyperbolic functions involved, it is necessary to work with very large and very small numbers, which poses an accuracy problem for the usual calculations with 16, and even 32, decimal places. To obtain correct results we had to require 128 or even 256 decimal places, which was possible when the calculations were performed using Maple (https://www.maplesoft.com).
The analytical solution obtained is used to explore the effect of the presence of a second non-conductive (gas) layer and the intensity of the magnetic field on the flow characteristics in horizontal flows and in concurrent and countercurrent inclined flows of the fluids. With this aim, we selected two representative two-phase systems—mercury–air and sodium–argon. Their physical properties are provided in Table 1 and Table 2, respectively. The viscosity ratio of the sodium–argon and the mercury–air systems are η 12 = 11.83 and η 12 = 81.87 , respectively. The higher value of the latter results from the lower viscosity of sodium compared to mercury (and somewhat higher viscosity of argon compared to air). Worth noting is that, for the same H a number, the magnetic field intensity applied in the case of mercury flow is about five times larger than that in argon (e.g., B0 = 0.01 T corresponds to Ha = 5.18 in mercury and to Ha = 25.23 in sodium).

4.1. Horizontal Channel

In a horizontal channel, Y 1,0 = 0 , whereby, for specified fluids (i.e., specified η 12 ), the solutions for the holdup h ~ and the dimensionless pressure gradient P f 1,0 are a function of the phases’ flow-rate ratio, Q 12 , and the Hartmann number, H a .
Figure 2 shows the effect of the magnetic field intensity variation of the holdup with the gas-to-liquid flow-rate ratio Q 21 = U 2 s / U 1 s for mercury–air and sodium–argon two-phase flows. Obviously, the liquid metal holdup decreases with the increase in the gas flow rate, and is higher in the mercury–air system of the larger viscosity ratio. The effect of the magnetic field on the holdup becomes significant for Ha > ~1. Due to the slowing down of the flow by the magnetic damping force, increasing the magnetic field intensity results in a higher holdup of the liquid metal and attenuates its decline with the increase in the gas flow rate. These trends are more pronounced for the same Ha in the more viscous mercury–air system.
Figure 3 shows the effect of the magnetic field intensity on the dimensionless frictional pressure gradient P f 1,0 (normalized with respect to the frictional pressure gradient of a single-phase liquid metal flow with the same U 1 s , but without a magnetic field, see Equation (15)). As expected, increasing Ha results in a higher frictional pressure gradient, whereby the P f 1,0 factor reaches a value of ~1000 for Ha = 103.7 (corresponding to B0 = 0.2 T and 0.041 T for mercury and sodium, respectively). It is well known that, in the absence of a magnetic field, the introduction of a gas flow to the flow of a viscous liquid can result in a lubrication effect (e.g., [31]), where P f 1,0 < 1 values are obtained for some range of low Q 21 . The lubrication effect increases with increasing η 12 (approaching a value of 0.25 in the TP geometry for large η 12 ). For example, in mercury–air flow and Ha = 0, P f 1,0 can reach a value of ~0.35 (65% reduction of the pressure gradient) for Q 21 = 0.05, and P f 1,0 < 1 up to Q 21 < ~10 (see also Figure 4b below for low Ha, which shows the power reduction). Indeed, a reduction of the pressure gradient factor P f 1,0 by the gas flow at low Q 21 values is also obtained for low values of Ha < ~1. However, the figure shows that, in the examined range of Ha > 5, the presence of a gas layer further augments the pressure gradient factor,   P f 1,0 , in particular, when Q 21 > 1 and Ha is relatively low. The sensitivity of the P f 1,0 factor to the gas flow is higher for the less viscous sodium–argon system. In both systems and high Ha, the gas flow rate has a rather small effect on the P f 1,0 factor over a wide range of gas flow rates.
The lubrication effect of the gas flow is noticed over a wider range of Ha when examining the P f 1 factor (Equation (17)), which is the dimensionless frictional pressure gradient normalized with respect to a single-phase liquid metal flow under the same magnetic field intensity. As shown in Figure 4a, the lubrication effect is more pronounced in the case of the more viscous mercury, where values of P f 1 < 1 are obtained up to Ha~25 for Q 21 < 1. The adverse effect of the gas flow on the pressure gradient at Q 21 > ~10 is also reflected in the P f 1 factor, which, for the same Ha, is higher for the sodium–argon system. The lubrication region is more visible in Figure 4b, where the pumping power factor, Po, is also shown for the mercury–air system for lower Ha (B0 = 0.001 T, Ha = 0.5183). Values of Po = P f 1 ( 1 + Q 21 ) < 1 indicate that pumping power saving can be achieved by introducing air flow to the conductive liquid flow. Note that for low Q 21   P o ~ P f 1 (both approach 1 at Q 21 = 0 ) . As shown in Figure 4b, significant power reduction can be achieved for low Ha. For example, for Ha = 0.5183, the minimal value Po = 0.3586 is obtained by adding 3.5% air to the mercury flow. For such low magnetic field strength, the P f 1 ~ P f 1,0 and the value of both is almost the same as that obtained for Ha = 0. The potential for power saving by the air flow is obviously reduced with increasing the magnetic field strength (negligible for Ha > 25), and the minimum Po is shifted to lower Q 21 values.
Further insight into the effect of the gas flow on the holdup and pressure gradient can be obtained by examining the dimensionless velocity profiles (normalized by the conductive liquid superficial velocity, U 1 s ). Figure 5a demonstrates the change in the velocity profile in the mercury–air system when the gas flow rate is increased from Q 21 = 1 to 10 for a constant Ha (=25.92, B0 = 0.05 T). As shown, in both cases, the air moves much faster than the mercury and its velocity profile is practically parabolic. The velocity gradients at the upper wall (hence the shear stress) are sharp and obviously increase with the gas flow rate. The velocity profiles in the mercury are depicted in Figure 5b for a wider range of Q 21 . The magnetic field results in a flat velocity profile in the bulk of the mercury layer, with sharp velocity gradients at the wall and at the interface. As shown, except for very low Q 21 , the maximum velocity is in the air layer, whereby, near the interface, the mercury is dragged by the air flow. The sharp velocity gradients at the interface indicate that high interfacial shear stresses are involved. Increasing the gas flow rate also results in a higher near-wall velocity gradient in the mercury. The increase in the wall shear stresses on both channel walls with the gas flow is reflected in the increased pressure drop factors,   P f 1,0 and P f 1 (in particular for Q 21 > ~ 10 ,  Figure 3 and Figure 4).
Figure 6 shows that, with increasing the intensity of the magnetic field (for a fixed value of Q 21 = 1 ), even though the mercury holdup increases, the mercury velocity gradient near the wall and at the interface becomes steeper, resulting in an increased pressure gradient at higher Ha. The mercury velocity rapidly grows from the bottom wall and quickly attains an almost flat profile, like in the single-phase Hartmann flow. However, contrary to the latter, in the upper part of the mercury layer the velocity does not decay, but steeply grows towards the air velocity at the interface. The thickness of the gas layer decreases at increased applied magnetic fields, while its maximum velocity and interfacial shear stress increase. In the experimental study by Lu et al. [10], such a flow configuration was characterized as a gas jet. This configuration is likely to enhance mixing between the two phases, potentially causing the flow regime to gradually transition to a mixed flow.
The effect of Q 21 and Ha on the shape of the velocity profile can be deduced by examining the variation in the maximal velocity and the interfacial velocity. Figure 7a shows the variation in the maximal velocity in the velocity profile vs. Q 21 for various Ha (B0 ranging from 0.001 to 0.1). The range of Q 21 where the maximal velocity is in the mercury layer is indicated by a dashed line, corresponding to the range of flow rates where the mercury flow drags the air at the interface. In this range of flow rates, the interfacial velocity is obviously lower than the maximal mercury velocity. However, as shown, this flow range is limited to low Q 21 and diminishes with increasing Ha. Examining the interfacial velocity in mercury–air flow vs. Q 21 for various values of Ha (see Figure 7b) shows that it increases with increasing Q 21 , with a rather low sensitivity to the magnetic field strength.

4.2. Inclined Channel

In inclined channels, the holdup and (dimensionless) frictional pressure gradient ( P f 1,0 ,   P f 1 ) are determined by four dimensionless parameters ( η 12 , Q 21 ,   H a , Y 1,0 ) . To obtain the total (dimensionless) pressure gradient factors ( P 1,0 , P 1 , or P 1 T ; see Equations (42) and (43)) the hydrostatic pressure gradient contribution should be accounted for. It can be calculated based on the holdup; however, the density ratio, ρ 12 = ρ 1 / ρ 2 , is then an additional required parameter. Depending on the values set for Q 21 and Y 1,0 , concurrent up-flow, concurrent down-flow or countercurrent flow of the liquid metal and gas are considered. In the following, the effects of the system parameters on the flow characteristics are demonstrated by referring to the mercury–air system, thus fixing the values of η 12 = ( 81.87 ) and ρ 12 (=13,534).

4.2.1. Concurrent Upward Flow

According to the selected flow configuration and coordinates (Figure 1), concurrent upward flow corresponds to negative values of the superficial velocities, U 2 s , U 1 s < 0; hence, Q 21 > 0 and Y 1,0 > 0 .
Figure 8a shows the effect of the magnetic field intensity on the variation in the holdup with the air-to-mercury flow-rate ratio for a fixed inclination parameter value of Y 1,0 = 2075.25 . The latter corresponds to a constant flow rate of the mercury ( U 1 s = 0.005 m/s) and a slight upward inclination of the channel, β = 0.2° (the channel height is 0.02 m). The variation in Q 21 , in this case, actually corresponds to the change in the air flow rate. A comparison of Figure 8 with Figure 2 shows the drastic effect of the slight inclination on the holdup. Due to the retarding gravity force, the metal flow is slowed down. Consequently, even for very shallow channel inclination and low mercury flow rate, the mercury occupies most of the flow cross-section over a wide range of air flow rates ( h ~ > 0.9 up to Q 21 = ~ 200 with low sensitivity to Ha). The holdup (and the hydrostatic pressure drop) steeply declines with the air flow rate for Q 21 > 2000. Then, for some range of high air flow rates, three different solutions for the holdup are obtained for the same air (and mercury) flow rates: the high holdup solution and two additional solutions of lower holdup values. An example of the triple solution is indicated by the three circles in Figure 8b. With increasing the magnetic field strength (i.e., higher Ha), the range of air flow rates where a triple solution is obtained diminishes and is shifted to higher air flow rates (see Figure 8a). For the tested parameter set, the triple solution is feasible up to Ha ~ 52 (Bo ~ 0.1 T). For even higher air flow rates, beyond the triple-solution region, only a single (low) holdup solution is obtained, and the sensitivity of its value to the magnetic field intensity is rather low. Figure 8b shows that, with increasing the mercury flow rate (i.e., reducing Y 1,0 ), the range of gas flow rates where triple solutions are obtained is shifted to a lower Q 21 , but becomes narrower. For U 1 s 0.1   m / s , only a single solution for the holdup is obtained even for low Ha (=0.518). For the same flow-rate ratio ( Q 21 ) the holdup of the mercury layer increases with reducing the mercury flow rate (i.e., increasing Y 1,0 ).
The possibility of obtaining multiple holdup solutions is typical to gravity-dominated two-phase systems, where the body force is of the order of the frictional pressure gradient (or higher), and the forces acting on the fluids can be balanced in more than one flow configuration. A demonstration of the velocity profiles associated with the three different solutions for the holdup is shown in Figure 9. The velocity profile in the air flow cross-section is almost parabolic, whereby the air flow practically determines the values of the pressure gradient and the interfacial shear stress. In the low holdup solution, the gravity force acting on the thin mercury layer is low, and the pressure gradient and air shear at the interface are capable of carrying upward the flow over the entire mercury layer. In the high holdup solution, the mercury experiences the highest counter-flow body force, whereby the pressure gradient and the air interfacial shear are insufficient to carry the entire mercury layer upward. Consequently, a back-flow of the mercury is observed near the lower channel wall, and the wall shear direction is reversed. As shown in the figure, the back-flow is reduced in the middle holdup solution, which results in a lower reversed wall shear.
Back-flow of the heavier fluid is a source of instability of the stratified flow configuration, as it introduces a large disturbance at the inlet, where the fluids are introduced into the channel. The effect of the magnetic field strength on the velocity profiles of the high and middle holdup solutions is demonstrated in Figure 10 (for Y 1,0 = 2075.25 ) . As shown, due to the magnetic field damping force, increasing Ha reduces the back-flow in both flow configurations, although the mercury holdup corresponding to the high holdup solution increases with Ha, and, therefore, is associated with a higher backward gravity force. For the examined parameters, triple-holdup solutions are obtained up to Ha~18 (B0 = 0.035 T). For higher Ha, a single (relatively high) holdup solution is obtained, yet with a significant region of mercury back-flow (e.g., Ha = 51.83 in Figure 10b). The effect of the air flow rate on the velocity profile in the mercury layer for Ha = 51.84 is demonstrated in Figure 10c,d. For this magnetic field intensity, a single solution for the holdup is obtained for any Q 21 . As shown, although the mercury holdup is reduced with increasing the air flow rate, the back-flow intensity increases. Note that the same dimensionless velocity profiles would be obtained for lower superficial velocities of the conductive layer and the gas, provided the value of Y 1,0 is the same (e.g., by referring to a smaller channel size and/or inclination, and/or density difference), and Q 21 and η 12 are unchanged.
The effect of the magnetic field strength and the air flow rate on the pressure gradient is examined in light of Figure 11. As expected, the frictional pressure gradient factor P f 1,0 (scaled with the frictional pressure gradient of single-phase mercury flow for Ha = 0) increases with Ha (see Figure 11a). Introducing the air flow further increases the frictional pressure gradient. In fact, at high Q 21 (beyond the triple-solution region, corresponding to the loop in the pressure gradient curves) the increase in the frictional pressure gradient is dominated by the gas flow and is insensitive to the magnetic field intensity. Examining the frictional pressure gradient factor P f 1 (i.e., scaled with the frictional pressure gradient of mercury flow under the same Ha, Figure 11b) reveals that the lubrication effect, obtained in the horizontal channel at low air flow rates (Figure 7b), is lost at even such a slight upward inclination of the channel. The variation in the hydrostatic pressure gradient factor, P g 1 (see Equation (42)), is shown in Figure 11c. Obviously, for a constant Ha, it follows the trend of the holdup variation with Q 21 (shown in Figure 7). With increasing Ha, the hydrostatic pressure gradient factor decreases, and values < 1 are reached at high Ha and Q 21 .
From the perspective of reducing the pressure gradient by the air flow the values of the total pressure gradient factor, P 1 T (i.e., d p / d x normalized with respect to the total pressure gradient of single-phase mercury flow with the same magnetic field, Equation (43)), should be examined. Figure 11d shows that, at low Q 21 , introducing the air flow practically does not affect the total pressure gradient. However, in the triple-solution region and in its vicinity, the steep reduction in the mercury holdup (and thus the hydrostatic pressure gradient), results in a pronounced reduction in the total pressure gradient compared to the mercury single-phase up-flow (e.g., by a factor of 0.122 at low Ha, and by a factor of 0.285 for B0 = 0.1 T, Ha = 51.83). The potential for pressure drop reduction in this region is reduced with the increase in the magnetic field strength.
The trends of the variation in the holdup with the channel upward inclination are similar to those obtained upon changing the mercury flow rate (Figure 8b). Obviously, an increase in Y 1,0 can be affected by reducing the flow rate of the conductive phase (mercury), by increasing the channel upward inclination, or by increasing the channel size. Figure 12a shows that, for the same Ha and maintaining the same flow rates of both phases (i.e., same Q 21 ) , the holdup of the conductive layer is higher at a steeper upward channel inclination. Also, with increasing β, the triple-solution region is shifted to higher air flow rates. In this region, owing to the lower values of the low and middle holdup solutions at steeper inclinations, smaller values of the total pressure-gradient factor, P 1 T , can be obtained (e.g., 0.00425 for β = 5° in Figure 12b). Obviously, the same range of Y 1,0 examined in Figure 12a can be obtained for a lower superficial velocity of the conductive layer in a smaller channel size, in which case the Q 21 range of the triple-solution region corresponds to lower gas dimensional velocity.
The concurrent down-flow corresponds to Y 1,0 < 0   U 1 s , U 2 s > 0 ( Q 21 > 0). In countercurrent flow, a net downward flow of conductive heavier fluid is considered, namely U 1 s > 0 and Y 1,0 < 0 , while the light fluid flows upward ( U 2 s < 0); hence, Q 21 < 0. Countercurrent flow is a basic configuration in many heat and mass transfer systems.

4.2.2. Countercurrent and Concurrent Downward Flows

Figure 13 demonstrates the effect of the magnetic field intensity on the countercurrent flow characteristics of mercury and air for a given inclination parameter value of Y 1,0 = 259.1 . The latter corresponds to a constant flow rate of the mercury ( U 1 s = 1 m/s) and a channel inclination of β = 5° (the channel height is 0.02 m). The variation in Q 21 can then be attributed to variation in U 2 s . As shown in the figure, countercurrent flow can be established for a limited range of sufficiently low (upward) air flow ( Q 21 < 0 ) and magnetic field strength. In this range, two distinct solutions for the holdup (and the corresponding flow characteristics) are obtained for fixed flow rates of the air and mercury, which merge to a single solution at the flooding point, beyond which countercurrent flow is not feasible. The countercurrent flow region diminishes with increasing Ha. For Ha > 53.7, countercurrent flow of mercury and air is not feasible (see also Figure 15b below). It is worth noting that the possibility of establishing the two high and low holdup configurations in the countercurrent region was demonstrated experimentally by Ullmann et al. [28] for Ha = 0. The flow configuration realized actually depends on the resistance at the heavy phase outlet. The high holdup configuration is obtained by increasing the resistance at the outlet while maintaining the same flow rates.
Figure 13b,c show the variation in frictional pressure gradient factors P f 1,0 , P f 1 with Q 21 and Ha. In the countercurrent flow region, the positive values of the frictional pressure gradient factors are associated with the low-holdup configuration, indicating that the frictional pressure gradient is dominated by the downward mercury flow. Negative values of these factors are obtained (usually, but not exclusively, in the high holdup configuration) for low Ha (e.g., Ha = 0.518 and 5.18 in Figure 13b), which indicates that the frictional pressure gradient is dominated by the upward gas flow. With increasing Ha, the difference between the P f 1,0 and P f 1 values associated with the low and high holdup configurations diminishes (both P f 1,0 values are positive). With high Ha values, the frictional pressure gradients become close to that obtained in single-phase mercury flow with the same Ha (i.e., the P f 1 values are close to 1). The corresponding P f 1,0 values at high Ha are in the concurrent flow outside the range of Figure 13b. The downward flow of the mercury is obviously assisted by the hydrostatic pressure, d p / d x g > 0 ( P g 1 < 0, Figure 13c, and, when P g 1 > P f 1 , no pump is needed to drive the mercury flow.
The difference between the flow characteristics associated with the two-holdup solutions in countercurrent flow can be elucidated by examining the velocity profiles. Figure 14 shows the effect of Ha on the velocity profiles of the two holdup solutions obtained for a downward flow of mercury ( U 1 s = 1 m/s) and a low upward gas flow ( U 2 s = 0.05   m / s , Q 21 = 0.05 ). The velocity profiles in the lower holdup solution (Figure 14a) show that the downward flow of the mercury drags the air downward near the interface (i.e., air back-flow region). With increasing Ha, the mercury holdup increases and its (downward) velocity is reduced, resulting in lower air back-flow near the interface. Figure 14b shows the velocity profiles in the corresponding upper holdup solution. For the low Q 21 considered, the mercury occupies most of the channel, and the air flows in a very thin layer. Here, the air flow drags the mercury upward near the interface (i.e., mercury back-flow region). Such velocity profiles are associated with negative values of P f 1,0 , P f 1 (see Figure 13b,c). As shown, the mercury back-flow in the upper holdup solution diminishes with increasing Ha (for Ha = 25.92, the mercury flows downward over the entire layer). As shown in Figure 13a, for higher upward air flow rates, the thickness of the air layer in the high holdup solution is much larger. Yet, for low Ha, the mercury back-flow phenomenon near the interface is sustained over a wide range of negative Q 21 corresponding to P f 1,0 , P f 1 < 0 (Figure 13b,c).
The variation in the holdup with the magnetic field intensity in concurrent down-flow of mercury and air is elucidated in view of Figure 15 ( Y 1,0 = 259.1 , corresponding to U 1 s = 1 m/s, β=5°, H = 0.02 m). Figure 15a implies that, apparently, a single solution for the holdup is obtained for specified Ha and Q 21 . As expected, the holdup increases with Ha and is reduced by increasing the air flow rate. However, the enlargement of the region of high holdups and low Q 21 (Figure 15b) shows that two additional solutions of higher holdups can be obtained in the range of Ha < 53.92. This range of Ha is slightly higher than the maximal Ha value for obtaining countercurrent flow (Ha~53.7). However, independently of the magnetic field intensity, for the tested inclination parameter value ( Y 1,0 = 259.1 ) , already, only the low holdup solution branch persists in the concurrent flow for Q 21 > 0.006 . Note that the same inclination parameter can be obtained for lower mercury superficial velocity, lower channel inclination, or smaller channel height.
The frictional pressure gradient factor P f 1 is shown in Figure 15c. For low Ha, and (even low) Q 21 values, which are already outside the triple-solution region (see Figure 15a), the low mercury holdup solution is associated with higher frictional pressure gradients compared to single-phase mercury flow under the same Ha. The effect of the airflow on the frictional pressure gradient diminishes with increasing Ha. Enlargement of the low Q 21 (triple solution) region is shown in Figure 15c, where P f 1 values associated with the two additional high holdup solutions at low Q 21 are shown (e.g., for Ha = 25.92 and 51.84). The two high holdup solutions in the triple-solution regions correspond to P f 1 ~ 1   (only slightly lower than the frictional pressure gradients of single-phase mercury under the same Ha) and are lower than that obtained for the low holdup solution. For sufficiently high Ha, only a single holdup solution is obtained in the entire concurrent down-flow region, for which the dimensionless frictional pressure gradient values are only slightly lower than 1 at low Q 21 , due to the air lubrication effect. Practically, for large Ha, the frictional pressure gradient factor,   P f 1 , in the concurrent region, is of the order of 1 in a wide range of Q 21 . For example, for Q 21 = 200, the frictional pressure gradient is only about 5% larger than the single-phase mercury flow. Obviously, it then keeps increasing with further increase in the air flow rate.
The difference between the velocity profiles of the three holdup solutions in concurrent down-flow is demonstrated in Figure 16. As shown, the velocity profile of the lower holdup solution (Figure 16a) is similar to that of the lower solution in countercurrent flow (Figure 14a). The difference is that here the net flow of air is downward, while there is a region of upward (back-flow) of the air near the upper channel surface. The steep mercury velocity gradient at the wall at low Ha is responsible for the large frictional pressure gradient associated with the lower holdup solution discussed above (with reference to Figure 15c). With increasing Ha, the holdup of the mercury increases, its velocity is reduced, and the velocity profile is flatter and becomes more similar to that observed in single-phase mercury flow under the same Ha. The back-flow in the thinner air layer diminishes, and eventually, for sufficiently large Ha, the flow in the entire gas layer is downward, and only one solution exists. A similar effect of Ha on the mercury velocity profile is observed in the middle and high holdup solutions (the latter is not shown), as both correspond to high mercury holdup and the velocity profiles in the mercury layer are similar. The main difference in the velocity profile of the middle and high holdup configurations is in the air layer (shown in Figure 16c,d). While in the upper solution, the flow in the entire air layer is downward for all Ha; in the middle holdup solution, back-flow of the air is still obtained. The increase in the mercury holdup with Ha results in a higher (positive) hydrostatic pressure gradient, which facilitates the downward flow of the heavier conductive liquid. Due to the opposite sign of the hydrostatic and frictional pressure gradient in concurrent down-flow, and also in countercurrent flow, the interpretation of the total pressure gradient value is more complicated than in concurrent up-flow. For example, for a specified single-phase flow rate of mercury in the particular channel size considered, the hydrostatic (dimensional) pressure gradient can be larger than the (negative) frictional pressure gradient. Hence, the total pressure gradient can attain positive values, indicating that the specified mercury flowrate is entirely driven by gravity. In fact, a restriction (valve) at the channel outlet is required to maintain the specified flow rate.
To analyze the effect of the air flow and the magnetic field intensity on the total pressure gradient, in concurrent down-flow and countercurrent flow, the values of d p / d x are normalized with respect to the hydrostatic pressure gradient in single-phase mercury flow, i.e., π T = ( d p d x ) / ρ 1 g s i n β . The obtained π T values are shown in Figure 17. The effect of the air flow rate on the total pressure gradient is demonstrated in Figure 17a, where the π T values are depicted vs. Q 21 . Positive values are obtained in the countercurrent flow region (as long as the magnitude of the magnetic field is low enough to enable countercurrent flow), indicating that the specified mercury downward flowrate is also entirely driven by gravity in the presence of air upward gas flow. Figure 17b shows the π T values when presented vs. the mercury holdup. The range of holdups that correspond to countercurrent flow extends from a holdup of almost 1 (marked by a square) to the holdup for which π T = 0   (marked by a diamond). For this holdup, the frictional pressure gradient for the air–mercury system with the specified mercury flow rate (i.e., specified   Y 1,0 ) is exactly balanced by the hydrostatic pressure gradient. This point corresponds to the lowest holdup solution obtained for Q 21 = 0   (see Figure 13a), where the air is circulating in the channel with a zero net flow. The high holdup configuration (indicated by the squares) at Q 21 = 0 also corresponds to gas circulation in the channel, however, with π T > 0   (see Figure 17a). As shown in Figure 17b, the holdup range of countercurrent flow diminishes with increasing Ha. Countercurrent flow is feasible for magnetic field intensities below a critical value, for which the π T value for single-phase mercury flow (i.e., holdup = 1) is 0 (i.e., the frictional pressure gradient for the single-phase mercury flow is balanced by the hydrostatic pressure gradient). For a higher Ha, the single-phase dimensionless pressure gradient is <1, and a pump is needed to drive the mercury flow, whereby countercurrent flow of air becomes unfeasible.
The negative values of π T in Figure 17 correspond to concurrent down-flow, and in the triple-solution region they correspond to the lower holdup solution. The figure shows that these π T values become more negative as the magnitude of the magnetic field strength and/or the downward air flow are increased, indicating that higher pumping power is needed to drive the concurrent downward mercury–air flow. For Ha values for which triple solutions in concurrent down-flow are feasible, the π T values of the two high holdup solutions (which are close to 1, and cannot be distinguished in Figure 17b) are positive, and are only slightly lower than the π T value of single-phase mercury flow for the same Ha. This indicates that in these configurations the mercury–air down-flow is driven by gravity. It is worth noting that the correspondence of π T = 0 and the low holdup solution of countercurrent flow at Q 21 = 0 is valid in cases of ρ 2 / ρ 1 . Otherwise, values π T > 0 can also be obtained in the low holdup solution at Q 21 = 0 as well as in concurrent flow ( Q 21 > 0 ).
The effect of channel inclination on the holdup and the total pressure gradient factor, π T , in both countercurrent and concurrent downward flow is illustrated in Figure 18. The range of β = 0.5° to 10° corresponds to Y 1,0 = −25.94 to −516.2. Decreasing the channel inclination reduces the range of airflow rates for countercurrent flow (Figure 18a). With decreasing channel inclination, the thickness of the mercury layer associated with the high holdup solution decreases, while thicker mercury layers are associated with the low holdup solution. Conversely, at shallower inclinations, the range where a triple-holdup solution is obtained in concurrent downward flow extends over a wider range of downward airflow rates (refer to Figure 18b, which depicts an enlargement of the high holdup region at low Q 21 ; the corresponding low holdup solution is shown in Figure 18a). While the total pressure gradient factor ( π T ) of the high and middle solutions is positive (values close to 1 for all the examined inclinations, Figure 18d), the π T value of the low holdup solution becomes more negative at shallower channel inclinations (see Figure 18c). This indicates that a pump is required to propel the concurrent downward mercury–air flow over most of the range of Q 21 > 0 , where the low holdup is the only solution of the flow equations.
The influence of channel inclination on the flow characteristics illustrated in Figure 18 can be used to anticipate the effect of the mercury flow rate or the density of the conductive fluid. This is evident from the fact that the same variation in Y 1,0 can be achieved by adjusting the mercury flow rate while keeping β constant, or by altering the density of the conductive layer ( ρ 2 is significantly less than ρ 1 and has minimal impact on the density difference). It is noteworthy that increasing the viscosity of the lighter layer (i.e., reducing η 12 ) has a similar effect to channel inclination (though not depicted), as it diminishes the countercurrent flow region and marginally expands the range of Q 21 for triple solutions in concurrent down-flow. However, for higher η 2 , the value of the single (low) holdup solution in concurrent down-flow (at a specified Q 21 ) is reduced. On the other hand, the impact of altering η 1 is more complex, as it also affects the values of Ha and Y 1,0 , and thus necessitates specific examination of the particular two-phase system of interest.

4.3. Numerical Calculations

As already mentioned above, the hyperbolic functions involved in the analytical solutions pose an accuracy problem for usual calculations with 16, and even 32, decimal places. To obtain correct results, the calculations were carried out in Maple with 128 or 256 decimal places. For example, consider the mercury–air system with H = 0.2   m ,   B 0 = 0.1   T . The holdup noticeably decreases below 1 for Q 21 = 100 and larger. The correct holdup values are calculated using 256 floating point digits, which yields h = 0.962 ,   0.919 ,   0.829 , and 0.645 for Q 121 = 10 2 ,   10 3 ,   10 4 , and 10 5 , respectively. The inaccurate results calculated by Maple using 128 digits are 0.594 ,   0.590 ,   0.585 , and 0.579 , respectively.
An alternative method of calculation of the holdups and the associated velocity profiles and pressure gradients is a numerical solution of the flow, Equations (1), (2), (24), and (25). The numerical schemes do not involve evaluation of the hyperbolic functions, and therefore can be carried out with a usual floating-point precision. We obtained numerical solutions using the finite difference and Chebyshev collocation methods. The Chebyshev collocation method is the same as that presented in [33], and the finite difference method applies central differences on an arbitrary stretched grid. The holdups and pressure drops are calculated using the superposition principle and the secant method as described in [34,35].
For the above example, using the finite difference method for the above parameters and at all values of Q 21 considered, the holdup converges to the value obtained by the analytical solution within the fourth decimal place with 200 uniformly distributed grid points, and within the fifth place with 700 points. Trying to improve the convergence, we applied the t a n h stretching,
x 0.5 + 0.5 t a n h s x 0.5 t a n h 0.5 s ,
near the boundaries and the interface. However, with the stretching parameter gradually increased up to s = 3 , we did not observe any noticeable improvement in the convergence.
The Chebyshev collocation method exhibits much faster, spectral convergence. Thus, for the same test cases, convergence within four decimal places is reached with the truncated series of 50 polynomials, while beyond 90 polynomials eight decimal places are readily converged. Obviously, the analytical solution is important for the validation of the numerical results and for assuring that all possible solutions are considered in inclined flows, where multiple steady-state configurations are feasible. An example of the holdup convergence study and comparison with the analytic solution is presented in Appendix A.
It is emphasized that, for example, the numerical study of stability of these flows is preferably conducted by using the Chebyshev collocation method, as demonstrated in [32,33]. At the same time, consideration of two-phase flows in bounded rectangular ducts or circular pipes (see [34,35]) will require lower-order methods. In these cases, the present results on convergence of the finite difference method will help estimate the computational resources required.

5. Concluding Remarks

We analyzed the characteristics of two-phase laminar stratified flows in horizontal and inclined channels under the effect of a transverse magnetic field, when the lower, heavier liquid is electrically conducting while the upper fluid is an electric insulator. The simplified channel geometry of two infinite parallel plates is considered, allowing for analytical solutions, and thus facilitating the analysis of system parameters’ effects on flow characteristics.
The flow is defined by prescribed flow rates of the fluids so that the pressure gradient needed to drive the flow, as well as the conducting layer holdup, is obtained together with the solution of the whole problem. We focused on two main issues:
(a)
The effect of the magnetic field on the liquid metal holdup and the pressure gradient driving the flow for prescribed flow rates of the fluids.
(b)
The effect of the magnetic field on the multiple states (i.e., different holdups, pressure gradients, etc., for the same operational conditions), which, under certain conditions, may exist in inclined stratified two-phase flows.
The main conclusions are as follows:
(1)
In horizontal flow, there exists a single state whose holdup increases with the increase in magnetic field strength. This increase is explained by the slowing down of the flow by the magnetic damping force. This effect is also valid in concurrent upward and downward flows for operational conditions that correspond to a single holdup solution of the flow governing equations.
(2)
In inclined concurrent flows, the high-density difference between the liquid metal and the gas results in high sensitivity of the flow configuration to the channel inclination, which allows for multiple states (three) for fixed flow rates under low-to-moderate magnetic field intensities.
(3)
In the triple-solution region, growth of the magnetic field strength leads to an increase in the lower and upper holdup values, but to a decrease in the holdup of the middle solution.
(4)
The triple-solution regions are associated with the possibility of a back-flow region (downward) of the heavy phase near the lower wall in upward flow, and a back-flow region (upward) of the gas near the upper wall. However, back-flow is not exclusively associated with multiple solutions and can also be obtained under conditions where a single solution exists.
(5)
Back-flow is a source of instability of the stratified flow configuration, as it introduces a large disturbance at the inlet, where the fluids are introduced into the channel. It can result in the formation of liquid metal slugs in an upward flow or entrainment of the gas into the liquid in a downward flow.
(6)
The magnetic field weakens the back-flow, and the velocity profile of the liquid metal flattens. Consequently, the ranges of flow rates where multiple solutions exist narrow, and then disappear in strong magnetic fields. In this respect, a magnetic field stabilizes the single holdup flow configuration.
(7)
In countercurrent flow, two distinct states exist up to a certain value of the phase’s flow-rate ratio, beyond which countercurrent flow is impossible. In most cases, in the state of the smaller holdup, the heavier liquid metal drags the lighter one near the interface, while in the other state of a larger holdup, the thinner gas layer drags the liquid metal near the interface.
(8)
With the increase in magnetic field strength, the countercurrent region diminishes, and beyond a certain Hartmann number, countercurrent flow is not feasible.
(9)
Compared to the classical single-phase Poiseuille flow, the magnetic field damping force increases the pressure gradient needed to reach the same flow rate. In liquid metal–gas horizontal stratified flows, the addition of a gas layer may yield a lubrication effect so that, compared to the single-phase Hartmann flow, the same flow rate of the metal liquid can be reached with a smaller pressure gradient.
(10)
In horizontal flows, the lubrication effect increases with increasing the liquid–metal/gas viscosity ratio and can reach about 75% pressure gradient reduction with high viscosity liquid metal under low magnetic field strength. The maximal lubrication effect is reached by adding a few percent gas flow rate to that of the liquid metal flow rate, and therefore can also be associated with pumping power reduction. However, the lubrication effect diminishes with increasing the Hartmann number (e.g., for mercury–air flow, it becomes insignificant for Ha > 25).
(11)
In upward concurrent flow, the introduction of gas flow does not result in a lubrication effect and, in fact, increases the frictional pressure gradient. However, as the hydrostatic pressure gradient is reduced in the presence of a gas layer, the total pressure gradient can be significantly reduced. Indeed, a reduction in the total pressure gradient compared to the single-phase flow of the liquid metal is obtained in the triple-solution region and its vicinity. As this region is associated with high gas/liquid flow-rate ratios, power reduction may not be obtained. In any case, the potential for a pressure drop reduction in this region is reduced with the increase in the magnetic field strength.
(12)
In concurrent downward flow, the addition of the gas layer also has a small effect on the frictional pressure gradient in the triple-solution region. However, as the driving force of the hydrostatic pressure is reduced in the presence of a gas layer, there is no benefit to adding gas to the flow of the liquid metal.
(13)
In countercurrent flows, the frictional pressure gradient is higher in the state of the smaller holdup. In both states, it is larger than in single-phase mercury flow under the same magnetic field strength.
The analytical solution has been used to validate the results obtained by the numerical solution of the flow equations, which is intended to be used for studying the effect of the magnetic field on the stability of the laminar stratified flow pattern in horizontal and inclined flows.

Author Contributions

Conceptualization, A.P., A.G., N.B. and A.U.; Software, A.P. and A.G.; Validation, A.P. and A.G.; Formal analysis, A.P., A.G. and N.B.; Resources, A.G. and N.B.; Data curation, A.P. and A.U.; Writing—original draft, N.B. and A.U.; Writing—review & editing, A.P.; Visualization, A.U.; Supervision, A.G.; Funding acquisition, N.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Israel Science Foundation (ISF) grant No 1363/23 and by the Israel Ministry of Aliyah and Integration (for A. Parfenov).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Latin Symbols
B 0 magnetic field fluxkg/(s2−A)
g gravitational accelerationm/s2
G pressure   gradient , = d p / d x Pa/m
h height of interface planem/s
h ~ Holdup , = h / H -
H channel height-
H a Hartmann   number , = B 0 H   σ 1 / η 1 -
p pressurePa
P dimensionless pressure gradient -
P o dimensionless pumping power-
Q volumetric flow ratem3/(m−s)
u velocitym/s
u ~ dimensionless velocity, = u / U 1 s -
U 1 s , 2 s superficial   velocity , = Q 1,2 / H m/s
x horizontal coordinatem
Y inclination parameter, defined in Equation (39)-
z vertical coordinatem
z ~ dimensionless   vertical   coordinate , = z / H -
Greek letters
β channel inclination angle to the horizontal
η dynamic viscositykg/(m-s)
π dimensionless pressure gradient normalized by the hydrostatic pressure gradient of single-phase flow -
ρ densitykg/m3
σ electric conductivity1/(Ω-m)
Subscripts
1 lower phase (conductive fluid)
2upper phase (gas)
12, 21 ratio   ( e . g . ,   η 12 = η 1 /   η 2 )
f frictional pressure
g gravitational (hydrostatic) pressure
s single phase
Superscripts
0value without magnetic field is used for normalization
1,2value in single-phase flow of fluid 1,2 is used for normalization
T normalized by the total pressure gradient in single-phase flow

Appendix A

Comparison between analytical and numerical solution for B 0 = 0.1 T ,   s = 3 .
Finite DifferencesChebyshev CollocationAnalytical
Q12Number of Grid NodesHoldupNumber of Collocation PointsHoldupHoldup
10−6200.2953100.2067
500.2944200.2885
1000.2942500.2940
2000.29411000.2940
5000.29411500.2940
10000.29402000.29400.2940
10−4200.8303100.7694
500.8296200.8253
1000.8294500.8293
2000.82931000.8293
5000.82931500.8293
10000.82932000.82930.8293
10−2200.9626100.9582
500.9624200.9621
1000.9624500.9624
2000.96241000.9624
5000.96241500.9624
10000.96242000.96240.9624
1200.9919100.9918
500.9919200.9919
1000.9919500.9919
2000.99191000.9919
5000.99191500.9919
10000.99192000.99190.9919

References

  1. Branover, H. Magnetohydrodynamic Flow in Ducts; John Wiley: New York, NY, USA, 1978. [Google Scholar]
  2. Kadid, F.Z.; Drid, S.; Abdessemed, R. Simulation of magnetohydrodynamic and thermal coupling in the linear induction MHD pump. J. Appl. Fluid Mech. 2011, 4, 51–57. [Google Scholar]
  3. Nikodijević, D.; Stamenković, Ž.; Milenković, D.; Blagojević, B.; Nikodijevic, J. Flow and heat transfer of two immiscible fluids in the presence of uniform inclined magnetic field. Math. Probl. Eng. 2011, 2011, 132302. [Google Scholar] [CrossRef]
  4. Kim, S.-H.; Park, J.-H.; Choi, H.-S.; Lee, S.-H. Power Generation Properties of Flow Nanogenerator with Mixture of Magnetic Nanofluid and Bubbles in Circulating System. IEEE Trans. Magn. 2017, 53, 4600904. [Google Scholar] [CrossRef]
  5. Wang, Y.; Cheng, K.; Xu, J.; Jing, W.; Huang, H.; Qin, J. A rapid performance prediction method for Two-Phase liquid metal MHD generators based on Quasi-One-Dimensional model. Therm. Sci. Eng. Prog. 2024, 47, 102258. [Google Scholar] [CrossRef]
  6. Haim, H.B.; Jianzhong, Z.; Shizhi, Q.; Yu, X. A magneto-hydrodynamically controlled fluidic network. Sens. Actuators B Chem. 2003, 88, 205–216. [Google Scholar]
  7. Hussameddine, S.K.; Martin, J.M.; Sang, W.J. Analytical prediction of flow field in magnetohydrodynamic-based microfluidic devices. J. Fluids Eng. 2008, 130, 091204. [Google Scholar]
  8. Yi, M.; Qian, S.; Bau, H. A magnetohydrodynamic chaotic stirrer. J. Fluid Mech. 2002, 468, 153–177. [Google Scholar] [CrossRef]
  9. Weston, M.C.; Gerner, M.D.; Fritsch, I. Magnetic fields for fluid motion. Anal. Chem. 2010, 82, 3411–3418. [Google Scholar] [CrossRef] [PubMed]
  10. Lu, P.; Ye, Q.; Yan, X.; Wei, J.; Huang, H. Study on the gas-liquid two-phase flow regime in the power generation channel of a liquid metal MHD system. Int. Commun. Heat Mass Transf. 2022, 137, 106217. [Google Scholar] [CrossRef]
  11. Wang, Y.; Huang, H.; Lu, P.; Liu, Z. Numerical Investigation of Gas–Liquid Metal Two-Phase Flow in a Multiple-Entrance Magnetohydrodynamic Generator. Ind. Eng. Chem. Res. 2022, 61, 4980–4995. [Google Scholar] [CrossRef]
  12. Lee, C.; Ho Kim, H. Velocity measurement of magnetic particles simultaneously affected by two-phase flow and an external magnetic field using dual-sided SPIM-mPIV. Chem. Eng. Sci. 2022, 252, 11727. [Google Scholar] [CrossRef]
  13. Khan, H.A.; Goharzadeh, A.; Jarrar, F. Effect of a uniform magnetic field on a two-phase air-ferrofluid slug flow. J. Magn. Magn. Mater. 2023, 580, 170944. [Google Scholar] [CrossRef]
  14. He, Y.; Wen, G.; Li, Q.; Jiao, F. Field-induced interfacial instabilities in a two-phase ferrofluid flow. Chem. Eng. J. 2024, 485, 14995. [Google Scholar] [CrossRef]
  15. Hartmann, J.; Lazarus, F. Experimental investigations on the flow of mercury in a homogeneous magnetic field. K. Dan. Vidensk. Selsk. Math. Fys. Medd. 1937, 15, 1–45. [Google Scholar]
  16. Davidson, P.A. An Introduction to Magnetohydrodynamics; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  17. Molokov, S.; Moreau, R.; Moffatt, K. Magnetohydrodynamics Historical Evolution and Trends Fluid Mechanics and Its Applications; Springer: Berlin/Heidelberg, Germany, 2007; FMIA; Volume 80. [Google Scholar]
  18. Qian, S.; Bau, H.H. Magneto-hydrodynamics based microfluidics. Mech. Res. Commun. 2009, 36, 10–21. [Google Scholar] [CrossRef] [PubMed]
  19. Kundu, B.; Saha, S. Review and Analysis of Electro-Magnetohydrodynamic Flow and Heat Transport in Microchannels. Energies 2022, 15, 7017. [Google Scholar] [CrossRef]
  20. Hartmann, J. Theory of the laminar flow of an electrically conductive liquid in a homogeneous magnetic field. K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 1937, 15, 1–28. [Google Scholar]
  21. Shail, R. On laminar two-phase flow in magnetohydrodynamics. Int. J. Eng. Sci. 1973, 11, 1103–1108. [Google Scholar] [CrossRef]
  22. Owen, R.G.; Hunt, J.C.R.; Collier, J.G. Magnetohydrodynamic pressure drop in ducted two-phase flows. Int. J. Multiph. Flow 1976, 3, 23–33. [Google Scholar] [CrossRef]
  23. Lohrasbi, J.; Sahai, V. Magnetohydrodynamic heat transfer in two-phase flow between parallel plates. Appl. Sci. Res. 1988, 45, 53–66. [Google Scholar] [CrossRef]
  24. Malashetty, M.S.; Umavathi, J.C. Two-phase magnetohydrodynamic flow and heat transfer in an inclined channel. Int. J. Multiph. Flow 1997, 23, 545–560. [Google Scholar] [CrossRef]
  25. Umavathi, J.C.; Liu, I.-C.; Prathap Kumar, J. Magnetohydrodynamic Poiseuille-Couette flow and heat transfer in an inclined channel. J. Mech. 2010, 26, 525–532. [Google Scholar] [CrossRef]
  26. Shah, N.A.; Alrabaiah, H.; Vieru, D.; Yook, S.-J. Induced magnetic field and viscous dissipation on flows of two immiscible fluids in a rectangular channel. Sci. Rep. 2022, 12, 39. [Google Scholar] [CrossRef] [PubMed]
  27. Landman, M.J. Non-unique holdup and pressure drop in two-phase stratified inclined pipe flow. Int. J. Multiph. Flow 1991, 17, 377–394. [Google Scholar] [CrossRef]
  28. Ullmann, A.; Zamir, M.; Ludmer, Z.; Brauner, N. Stratified laminar countercurrent flow of two liquid phases in inclined tubes. Int. J. Multiph. Flow 2003, 29, 1583–1604. [Google Scholar] [CrossRef]
  29. Ullmann, A.; Zamir, M.; Gat, S.; Brauner, N. Multi-holdups in co-current stratified flow in inclined tubes. Int. J. Multiph. Flow 2003, 29, 1565–1581. [Google Scholar] [CrossRef]
  30. Ullmann, A.; Goldstien, A.; Zamir, M.; Brauner, N. Closure relations for the shear stresses in two-fluid models for stratified laminar flows. Int. J. Multiph. Flow 2004, 30, 877–900. [Google Scholar] [CrossRef]
  31. Goldstein, A.; Ullmann, A.; Brauner, N. Characteristics of stratified laminar flows in inclined pipes. Int. J. Multiph. Flow 2015, 75, 267–287. [Google Scholar] [CrossRef]
  32. Barmak, I.; Gelfgat, A.; Vitoshkin, H.; Ullmann, A.; Brauner, N. Stability of stratified two-phase flows in horizontal channels. Phys. Fluids 2016, 28, 044101. [Google Scholar] [CrossRef]
  33. Barmak, I.; Gelfgat, A.Y.; Ullmann, A.; Brauner, N. Stability of stratified two-phase flows in inclined channels. Phys. Fluids 2016, 28, 084101. [Google Scholar] [CrossRef]
  34. Gelfgat, A.; Brauner, N. Instability of stratified two-phase flows in rectangular ducts. Int. J. Multiph. Flow 2020, 131, 103395. [Google Scholar] [CrossRef]
  35. Gelfgat, A.; Barmak, I.; Brauner, N. Instability of stratified two-phase flows in inclined rectangular ducts. Int. J. Multiph. Flow 2021, 138, 103586. [Google Scholar] [CrossRef]
Figure 1. Sketch of the flow configuration and coordinates.
Figure 1. Sketch of the flow configuration and coordinates.
Fluids 09 00129 g001
Figure 2. The effect of the magnetic field intensity on the liquid metal holdup, h ~ . Comparison of holdup vs. the gas-to-liquid flow-rate ratio in the mercury–air (solid line) and sodium–argon (dashed line) systems for the same Ha. The range of Ha considered corresponds to B0 = 0.01 T to 0.2 T and B0 = 0.0020544 T to 0.041088 T in mercury and sodium, respectively.
Figure 2. The effect of the magnetic field intensity on the liquid metal holdup, h ~ . Comparison of holdup vs. the gas-to-liquid flow-rate ratio in the mercury–air (solid line) and sodium–argon (dashed line) systems for the same Ha. The range of Ha considered corresponds to B0 = 0.01 T to 0.2 T and B0 = 0.0020544 T to 0.041088 T in mercury and sodium, respectively.
Fluids 09 00129 g002
Figure 3. The effect of the magnetic field intensity on the P f 1 , 0 factor (with respect to the pressure gradient in single-phase liquid metal flow without a magnetic field)—comparison of P f 1,0 vs. the gas-to-liquid flow-rate ratio in the mercury–air (solid line) and sodium–argon (dashed line) systems for the same Ha.
Figure 3. The effect of the magnetic field intensity on the P f 1 , 0 factor (with respect to the pressure gradient in single-phase liquid metal flow without a magnetic field)—comparison of P f 1,0 vs. the gas-to-liquid flow-rate ratio in the mercury–air (solid line) and sodium–argon (dashed line) systems for the same Ha.
Fluids 09 00129 g003
Figure 4. The effect of the magnetic field strength on (a) the P f 1 factor (with respect to the pressure gradient in single-phase liquid metal flow under the same magnetic field, see Equation (17))—comparison of P f 1 vs. the gas-to-liquid flow-rate ratio in the mercury–air (solid line) and sodium–argon (dashed line) systems for the same Ha; (b) the power factor, Po (Equation (18)), vs. the gas-to-liquid flow-rate ratio in the mercury–air system.
Figure 4. The effect of the magnetic field strength on (a) the P f 1 factor (with respect to the pressure gradient in single-phase liquid metal flow under the same magnetic field, see Equation (17))—comparison of P f 1 vs. the gas-to-liquid flow-rate ratio in the mercury–air (solid line) and sodium–argon (dashed line) systems for the same Ha; (b) the power factor, Po (Equation (18)), vs. the gas-to-liquid flow-rate ratio in the mercury–air system.
Fluids 09 00129 g004
Figure 5. (a) The effect of the air-to-mercury flow-rate ratio on the dimensionless velocity profile (scaled with respect to mercury superficial velocity). (b) Enlargement of the mercury domain. Ha = 25.92. Mercury—red curves, air—blue curves.
Figure 5. (a) The effect of the air-to-mercury flow-rate ratio on the dimensionless velocity profile (scaled with respect to mercury superficial velocity). (b) Enlargement of the mercury domain. Ha = 25.92. Mercury—red curves, air—blue curves.
Fluids 09 00129 g005
Figure 6. The effect of the magnetic field intensity, Ha, on the mercury dimensionless velocity profile (scaled with respect to mercury superficial velocity), Q 21 = 1 . Mercury—red curves, air—blue curves.
Figure 6. The effect of the magnetic field intensity, Ha, on the mercury dimensionless velocity profile (scaled with respect to mercury superficial velocity), Q 21 = 1 . Mercury—red curves, air—blue curves.
Fluids 09 00129 g006
Figure 7. (a) The dimensionless maximal velocity vs. Q 21 for various Ha. The range of Q 21 where the maximal velocity is in the mercury layer is indicated by the dashed line. For large Q 21 values, the maximal velocity is always in the air and increases proportionally to Q 21 0.6737 . (b) The interfacial velocity in mercury–air flow vs. Q 21 for various values of Ha.
Figure 7. (a) The dimensionless maximal velocity vs. Q 21 for various Ha. The range of Q 21 where the maximal velocity is in the mercury layer is indicated by the dashed line. For large Q 21 values, the maximal velocity is always in the air and increases proportionally to Q 21 0.6737 . (b) The interfacial velocity in mercury–air flow vs. Q 21 for various values of Ha.
Fluids 09 00129 g007
Figure 8. Variation in the mercury holdup in concurrent up-flow of mercury and air: (a) Effect of the magnetic field strength, Ha ( Y 1,0 = 2075.25 , corresponding to U 1 s = 0.005 m/s β = 0.2° and H = 0.02 m). (b) Effect of the mercury superficial velocity for Ha = 0.518 (Y1,0 values in the range of 20.75 to 10376). An example of triple solutions is shown by the 3 circles on the holdup curve corresponding to U 1 s = 0.005 m/s.
Figure 8. Variation in the mercury holdup in concurrent up-flow of mercury and air: (a) Effect of the magnetic field strength, Ha ( Y 1,0 = 2075.25 , corresponding to U 1 s = 0.005 m/s β = 0.2° and H = 0.02 m). (b) Effect of the mercury superficial velocity for Ha = 0.518 (Y1,0 values in the range of 20.75 to 10376). An example of triple solutions is shown by the 3 circles on the holdup curve corresponding to U 1 s = 0.005 m/s.
Fluids 09 00129 g008
Figure 9. Velocity profiles (dimensionless) in the 3-holdup solutions of mercury–air up-flow. (a) (Ha = 0.518, Y 1,0 = 2075.25 ,   Q 21 = 2 × 10 4   U 1 s = 0.005 m/s, U 2 s = 100 m/s. (b) Enlargement of the velocity profiles in the mercury layer. Mercury—red curves, air—blue curves.
Figure 9. Velocity profiles (dimensionless) in the 3-holdup solutions of mercury–air up-flow. (a) (Ha = 0.518, Y 1,0 = 2075.25 ,   Q 21 = 2 × 10 4   U 1 s = 0.005 m/s, U 2 s = 100 m/s. (b) Enlargement of the velocity profiles in the mercury layer. Mercury—red curves, air—blue curves.
Fluids 09 00129 g009
Figure 10. Velocity profiles (dimensionless) in mercury–air up-flow ( Y 1,0 = 2075.25 ,   U 1 s = 0.005 m/s). (a) Effect of Ha on the mercury velocity profiles in the high and middle holdup solutions (for U 2 s = 100 m/s) . Triple solutions are obtained up to Ha~18 (b) Velocity profiles in the mercury layer for higher Ha for U 2 s = 100 m/s, where a single holdup solution is obtained. (c) Effect of the air superficial velocity on the mercury velocity profiles at high magnetic field strength, Ha = 51.84 (B0 = 0.1T). (d) Velocity profile in the air layer. Mercury—red curves, air—blue curves.
Figure 10. Velocity profiles (dimensionless) in mercury–air up-flow ( Y 1,0 = 2075.25 ,   U 1 s = 0.005 m/s). (a) Effect of Ha on the mercury velocity profiles in the high and middle holdup solutions (for U 2 s = 100 m/s) . Triple solutions are obtained up to Ha~18 (b) Velocity profiles in the mercury layer for higher Ha for U 2 s = 100 m/s, where a single holdup solution is obtained. (c) Effect of the air superficial velocity on the mercury velocity profiles at high magnetic field strength, Ha = 51.84 (B0 = 0.1T). (d) Velocity profile in the air layer. Mercury—red curves, air—blue curves.
Fluids 09 00129 g010
Figure 11. Effect of the magnetic field strength, Ha, on the pressure gradient in concurrent up-flow of mercury and air ( Y 1,0 = 2075.25 , corresponding to U 1 s = 0.005 m/s, β = 0.2° and H = 0.02 m). (a) Frictional pressure gradient factor P f 1,0 (with respect to the frictional pressure gradient of SP mercury flow and Ha = 0). (b,c) Frictional and hydrostatic pressure gradient factor, P f 1 and P g 1 , respectively (with respect to the frictional pressure gradient of SP mercury flow with the same Ha). (d) Total pressure gradient normalized with respect to the total pressure gradient of SP mercury flow with the same Ha, P 1 T = P 1 / ( 1 + ρ 2 ρ Y 1,0 P f s 1,0 ) .
Figure 11. Effect of the magnetic field strength, Ha, on the pressure gradient in concurrent up-flow of mercury and air ( Y 1,0 = 2075.25 , corresponding to U 1 s = 0.005 m/s, β = 0.2° and H = 0.02 m). (a) Frictional pressure gradient factor P f 1,0 (with respect to the frictional pressure gradient of SP mercury flow and Ha = 0). (b,c) Frictional and hydrostatic pressure gradient factor, P f 1 and P g 1 , respectively (with respect to the frictional pressure gradient of SP mercury flow with the same Ha). (d) Total pressure gradient normalized with respect to the total pressure gradient of SP mercury flow with the same Ha, P 1 T = P 1 / ( 1 + ρ 2 ρ Y 1,0 P f s 1,0 ) .
Fluids 09 00129 g011
Figure 12. Effect of the channel inclination on the holdup (a) and on the total pressure gradient factor, P 1 T (b), for mercury–air flow, Ha = 5.18, U 1 s = 0.005 m/s, and H = 0.02 m ( Y 1,0 values in the range of 2075.25 to 51815.6, corresponding to β ranging from 0.2° to 5°).
Figure 12. Effect of the channel inclination on the holdup (a) and on the total pressure gradient factor, P 1 T (b), for mercury–air flow, Ha = 5.18, U 1 s = 0.005 m/s, and H = 0.02 m ( Y 1,0 values in the range of 2075.25 to 51815.6, corresponding to β ranging from 0.2° to 5°).
Fluids 09 00129 g012
Figure 13. Countercurrent mercury–air flow—effect of the magnetic field intensity (Ha) and air flow rate on the countercurrent flow characteristics: (a) Mercury holdup (where the two solutions and the flooding point for Ha = 25.92 are indicated by circles). (bd) Frictional and hydrostatic pressure gradient factors, P f 1,0 , P f 1 , and P g 1 , respectively. ( Y 1,0 = 259.1 , corresponding to U 1 s = 1 m/s, β = 5°, H = 0.02 m).
Figure 13. Countercurrent mercury–air flow—effect of the magnetic field intensity (Ha) and air flow rate on the countercurrent flow characteristics: (a) Mercury holdup (where the two solutions and the flooding point for Ha = 25.92 are indicated by circles). (bd) Frictional and hydrostatic pressure gradient factors, P f 1,0 , P f 1 , and P g 1 , respectively. ( Y 1,0 = 259.1 , corresponding to U 1 s = 1 m/s, β = 5°, H = 0.02 m).
Fluids 09 00129 g013
Figure 14. Effect of the magnetic field intensity on the (dimensionless) velocity profiles in the 2-holdup solutions of mercury–air countercurrent flow: (a) lower holdup solution; (b) upper holdup solution ( Y 1,0 = 259.1 , corresponding to U 1 s = 1 m/s, U 2 s = 0.05   m / s   β = 5° H = 0.02 m). Mercury- red curves, air-blue curves.
Figure 14. Effect of the magnetic field intensity on the (dimensionless) velocity profiles in the 2-holdup solutions of mercury–air countercurrent flow: (a) lower holdup solution; (b) upper holdup solution ( Y 1,0 = 259.1 , corresponding to U 1 s = 1 m/s, U 2 s = 0.05   m / s   β = 5° H = 0.02 m). Mercury- red curves, air-blue curves.
Fluids 09 00129 g014
Figure 15. Concurrent mercury–air down-flow—effect of the magnetic field intensity (Ha) and air flow rate on the concurrent flow characteristics: (a,b) Mercury holdup. (c,d) Frictional pressure gradient factor, P f 1 . The r.h.s. figures show enlargement of the low Q 21 region. ( Y 1,0 = 259.1 , corresponding to U 1 s = 1 m/s, β = 5°, H = 0.02 m.)
Figure 15. Concurrent mercury–air down-flow—effect of the magnetic field intensity (Ha) and air flow rate on the concurrent flow characteristics: (a,b) Mercury holdup. (c,d) Frictional pressure gradient factor, P f 1 . The r.h.s. figures show enlargement of the low Q 21 region. ( Y 1,0 = 259.1 , corresponding to U 1 s = 1 m/s, β = 5°, H = 0.02 m.)
Fluids 09 00129 g015
Figure 16. Effect of the magnetic field intensity on the (dimensionless) velocity profiles in the 3-holdup solutions of mercury–air concurrent down-flow: (a) lower holdup solution (for Ha = 17.97 only one solution exists); (b) middle holdup solution; (c,d) enlargement of the air velocity profile in the middle and upper holdup solutions. ( U 1 s = 0.005   m / s   a n d     Y 1,0 = 25.91 , corresponding to U 1 s = 10 m/s, β = 5°, H = 0.02 m.) Mercury—red curves, air—blue curves.
Figure 16. Effect of the magnetic field intensity on the (dimensionless) velocity profiles in the 3-holdup solutions of mercury–air concurrent down-flow: (a) lower holdup solution (for Ha = 17.97 only one solution exists); (b) middle holdup solution; (c,d) enlargement of the air velocity profile in the middle and upper holdup solutions. ( U 1 s = 0.005   m / s   a n d     Y 1,0 = 25.91 , corresponding to U 1 s = 10 m/s, β = 5°, H = 0.02 m.) Mercury—red curves, air—blue curves.
Fluids 09 00129 g016
Figure 17. The effect of the magnetic field intensity (Ha) and the air flow rate ( Q 21 ) on the dimensionless total pressure gradient factor π T (normalized by the hydrostatic pressure gradient in single-phase mercury flow) in concurrent down-flow and countercurrent flow ( Y 1,0 = 259.1 , corresponding to U 1 s = 1 m/s, β = 5°, H = 0.02 m). (a) π T vs.   Q 21 for various Ha. (b) π T vs. the holdup. For each Ha the countercurrent flow extends from a holdup of almost 1 (marked by squares) to the holdup for which π T = 0   (marked by diamonds).
Figure 17. The effect of the magnetic field intensity (Ha) and the air flow rate ( Q 21 ) on the dimensionless total pressure gradient factor π T (normalized by the hydrostatic pressure gradient in single-phase mercury flow) in concurrent down-flow and countercurrent flow ( Y 1,0 = 259.1 , corresponding to U 1 s = 1 m/s, β = 5°, H = 0.02 m). (a) π T vs.   Q 21 for various Ha. (b) π T vs. the holdup. For each Ha the countercurrent flow extends from a holdup of almost 1 (marked by squares) to the holdup for which π T = 0   (marked by diamonds).
Fluids 09 00129 g017
Figure 18. The effect of the channel inclination and the air flow rate ( Q 21 ) on the holdup (a,b) and dimensionless total pressure gradient factor π T (c,d) in countercurrent and concurrent down-flow. The r.h.s. figures are enlargement of the low Q 21 to show the holdup and π T values of the high and middle holdup solutions in the triple-solution region of concurrent down-flow (Ha = 5.184, U 1 s = 0.005   m / s ; the values of Y 1,0 correspond to β = 0.5°, 1°, 2°, 5°, 10°).
Figure 18. The effect of the channel inclination and the air flow rate ( Q 21 ) on the holdup (a,b) and dimensionless total pressure gradient factor π T (c,d) in countercurrent and concurrent down-flow. The r.h.s. figures are enlargement of the low Q 21 to show the holdup and π T values of the high and middle holdup solutions in the triple-solution region of concurrent down-flow (Ha = 5.184, U 1 s = 0.005   m / s ; the values of Y 1,0 correspond to β = 0.5°, 1°, 2°, 5°, 10°).
Fluids 09 00129 g018
Table 1. Properties of mercury (Hg) and air at room temperature.
Table 1. Properties of mercury (Hg) and air at room temperature.
NameNotationMercury (Hg)Air
Densityρ13,534 kg/m31.2 kg/m3
Dynamic viscosityη0.00149 kg/(m∙s)1.8 × 10−5 kg/(m∙s)
Kinematic viscosityν = η/ρ1.1 × 10−7 m2/s1.8 × 10−5 m2/s
Electric conductivityσ106 1/(Ω∙m)
Surface tensionγ0.4589 N/m
Table 2. Properties of liquid sodium (Na) and argon (Ar) at 400 °C.
Table 2. Properties of liquid sodium (Na) and argon (Ar) at 400 °C.
NameNotationSodium (Na)Argon (Ar)
Densityρ856 kg/m30.713
Dynamic viscosityη0.000284 kg/(m∙s)2.4 × 10−5
Kinematic viscosityν = η/ρ3.32 × 10−7 m2/s3.37 × 10−5
Electric conductivityσ4.52 × 106 1/(Ω∙m)
Surface tensionγ0.161 N/m
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Parfenov, A.; Gelfgat, A.; Ullmann, A.; Brauner, N. Hartmann Flow of Two-Layered Fluids in Horizontal and Inclined Channels. Fluids 2024, 9, 129. https://doi.org/10.3390/fluids9060129

AMA Style

Parfenov A, Gelfgat A, Ullmann A, Brauner N. Hartmann Flow of Two-Layered Fluids in Horizontal and Inclined Channels. Fluids. 2024; 9(6):129. https://doi.org/10.3390/fluids9060129

Chicago/Turabian Style

Parfenov, Arseniy, Alexander Gelfgat, Amos Ullmann, and Neima Brauner. 2024. "Hartmann Flow of Two-Layered Fluids in Horizontal and Inclined Channels" Fluids 9, no. 6: 129. https://doi.org/10.3390/fluids9060129

APA Style

Parfenov, A., Gelfgat, A., Ullmann, A., & Brauner, N. (2024). Hartmann Flow of Two-Layered Fluids in Horizontal and Inclined Channels. Fluids, 9(6), 129. https://doi.org/10.3390/fluids9060129

Article Metrics

Back to TopTop