1. Introduction
Flow of electric charges within or with a fluid or alternatively freely moving in free space is being analyzed as part of magneto-fluid dynamics (MFD). The latter appears also under the acronym MHD (magneto-hydrodynamics) although inaccurately linked to water (hydro) flow due to historical reasons. It applies to plasmas, liquid metals as well as beams of charges moving in free space. The present paper deals with the theoretical demonstration of two linked MFD effects that are being shown to be analogous to fluid flow in a rotating frame of reference. The flow of fluids in a rotating frame of reference has been studied extensively and has applications in geophysics, astrophysics, as well as in engineering. The specific effects related to rotating flows are predominantly a result of centripetal and Coriolis accelerations (Greenspan [
1]) as well as possibly centrifugal buoyancy Vadasz [
2]. The Coriolis effect and the resulting vortex formations have been identified theoretically as well as experimentally. Amar et al. [
3] demonstrated the latter numerically as well as analytically and compared their results with experimental data. They focused on the separation between geostrophic flow and Ekman and Stewartson boundary layers. Asymptotic analyses of rotating flows identify Taylor-Proudman columns and two-dimensional flow at the leading order and Ekman as well as Stewartson boundary layers for higher order corrections [
1], results that were confirmed numerically as well as experimentally (Subbotin et al. [
4], and Burmann and Noir [
5]). The explanation of the appearance of von Karman vortex streets around invisible bluff bodies as captured by satellite images over certain islands in the Atlantic and Pacific Oceans was provided by Vadasz [
6] in terms of Taylor-Proudman columns. Sarkar et al. [
7] investigated the effect of a magnetic field on Taylor-Proudman columns in a rotating electrically conducting fluid (MFD). They concluded that the “application of a magnetic field” “suppresses the Taylor column” in certain circumstances. In an electrically conducting fluid the balance between the Coriolis acceleration and the Lorentz force is the mechanism that controls the Taylor-Proudman column. Extensive research results on the problem of natural convection due to centrifugal buoyancy in rotating porous media were presented by Vadasz [
8,
9,
10,
11,
12,
13] and by Vadasz and Govender [
14]. Vadasz [
8] focused on the centrifugal buoyancy in a porous layer distant from the axis of rotation, Vadasz [
9] analyzed the Coriolis effect on a rotating porous layer heated from below via linear as well as weak nonlinear methods. Other investigations of centrifugal buoyancy in rotating porous media were performed by Saravanan and Vigneshwaran [
15] and Kang et al. [
16].
The effect of the magnetic field on the flow of conducting fluids has been studied extensively, e.g., identifying the law of isorotation (Allen et al. [
17]). An analogy between the Taylor-Proudman theorem for rotating fluids and the law of isorotation for MFD might evolve as a consequence of the present paper’s derivations.
The analogy between the Coriolis and Lorentz forces is not new, e.g., the analogy between the gyrocompass and the magnetic compass was demonstrated by Opat [
18]. In the present paper it is demonstrated also for fluids rather than compasses.
The next section describes the analogy between rotating and magnetic (MFD) flows followed by the theoretical demonstration on how electromagnetic columns are being created when a flowing fluid carrying electric charges is exposed to an externally imposed magnetic field. The consequent two-dimensional magnetostrophic flow then emerges. The final section presents the fact that the conditions required for the electromagnetic columns and magnetostrophic flow to emerge are well satisfied in reality and should then be observable in lab experiments.
2. Analogy between Magneto-Fluid Dynamics (MFD/MHD) and Rotating Flows
The equations governing the isothermal compressible flow in a rotating frame of reference are the continuity and momentum equations presented in the form
where
is the velocity vector,
and
are unit vectors in the
, and
directions, respectively,
is the constant angular velocity of rotation,
is pressure,
is the dynamic viscosity, and
is the position vector. By assuming the fluid to be barotropic and using a linear relationship between pressure and density (this assumption applies to isothermal conditions for an ideal gas, and approximately also for isentropic conditions of the latter, and for liquids) in the form
where
or
. Moving the Coriolis and centripetal terms to the right-hand side of the equation, and dividing Equation (2) by
it yields
where
is the kinematic viscosity. The centripetal acceleration term
has a potential and can be therefore moved under the gradient term in the form
It becomes appealing now to define a generalized reduced pressure term (a specific kinetic energy) in the form
Substituting (6) into (5) leads to
The momentum equation for an MFD fluid, i.e., a fluid that carries electric charges (even when as a whole it is neutrally charged) includes the electric and magnetic fields via the Lorentz force, and when neglecting the gravitational field, the latter being much weaker than the electromagnetic fields can be presented in the form
where
is the electric field due to the distributed charges in the fluid,
is the velocity of the charges, and
is the induced magnetic field that results from the electric current due to the moving charges. When the fluid is exposed to an external magnetic field then an additional term in the form
is to be added in Equation (8). Substituting Equation (3) into Equation (8), dividing Equation (8) by
, and using the definition of the generalized reduced pressure (6) for the case without rotation (i.e.,
) produces
where
is assumed to be constant and equal to the ratio between the total mass of the electric charges to the total electric charge, i.e.,
, and
is the electric charge sign function. Comparing Equation (9) for the magnetic-fluid dynamics flow to Equation (7) for non-magnetic fluid dynamics in a rotating frame of reference one observes the analogy where the gravitational field term
was replaced by the electric field term
, and the Coriolis acceleration term
was replaced by the magnetic term
. Consequently it is plausible to anticipate that the magnetic field term
in MFD has an identical effect as the angular velocity of rotation term
in the Coriolis acceleration of the rotating flow. Note that even the units of
are
identical to the units of
. Since the analogy between MFD flow and rotating flows was established by this comparison, one can therefore anticipate similar effects in MFD flow as in the corresponding Coriolis effects for rotating flow.
A similar analogy applies if the fluid is exposed to an externally imposed magnetic field. Then Equation (9) takes the form
where
is the imposed magnetic field, which is constant, i.e.,
, and
is a unit vector in the direction of the imposed magnetic field.
3. Emerging Electromagnetic Columns in an Electrically Charged Fluid Flowing Subject to an Imposed Magnetic Field
When the fluid is incompressible, or when the fluid can be assumed to be weakly compressible (i.e., the density does not change but only when in a body force term in the momentum equation) the continuity Equation (1) transforms into
By considering the flow subject to an imposed magnetic field, then Equation (11) is to be solved together with Maxwell equations presented in the form
Coulomb law in field form
where
is the permittivity of vacuum.
Ampere law
where
is the speed of light in vacuum, and
is the permeability of free space.
Faraday law of induction
Gauss law for the magnetic field
Note that does not appear in Equations (14)–(16) because it is constant and therefore its derivatives vanish.
Converting the equations into a dimensionless form requires the introduction of scales or characteristic values. The only requirement from such scales is that they are non-zero constants. Selecting , , , , , and as the length, mass density, charge density, velocity, electric charge, and mass scales, respectively, one can introduce the following additional scales in terms of the latter , , , , . One can also relate the mass density scale to the mass scale, i.e., and the charge density scale to the charge scale, i.e., . By using these scales the following dimensionless variables are being defined as follows , ,
,
,
,
,
,
,
. Substituting these dimensionless variables into Equations (11) and (12) transforms the latter into the following dimensionless form
where the following dimensionless groups emerged
Dividing Equation (18) by
produces (after dropping the symbol * as now all variables are dimensionless)
where a rescaled pressure was introduced in the form
and where new dimensionless groups emerged as follows
Assuming steady state, and
which implies
. Then if also
(or
) which implies
, and if
(or
), and if also
(or
) then neglecting the terms in Equation (24) corresponding to the small coefficients leads to
since
by definition presented in Equation (10). Selecting the coordinates axes in a such a way that the direction of the imposed magnetic field
is aligned with the negative direction of the
-axis, i.e.,
transforms (27) into
Taking the curl () of Equation (28) produces
However, by using Equation (17) and evaluating the curl in (29) leads to
Evaluating the dot product yields
Equation (31) is the electromagnetic equivalent to the Taylor-Proudman theorem from rotating flows. It implies in particular that
and since
at
due to impermeability boundary conditions it means that
for all values of
. This indicates that a flow over an object aligned with the imposed magnetic field as presented in
Figure 1 is impossible because it introduces a vertical component of velocity,
, and consequently is violating Equation (31). The conclusion is that the flow will adjust around the object as presented in
Figure 2a.
However, this flow pattern is also independent of z because Equation (31) implies
and consequently the flow pattern extends over the whole height creating a fluid column above the object which behaves like a solid body, as presented in
Figure 2b.
4. Magnetostrophic Flow
The lack of a vertical component of velocity leads to two-dimensional flow and then it becomes appealing to introduce a stream function that causes the continuity Equation (17) to be identically satisfied, i.e.,
Substituting (33) in Equation (28) after performing the vector product on its left-hand side produces
where
is the horizontal gradient operator.
As both the generalized reduced pressure, , and the stream function, , can be related to an arbitrary reference value, the conclusion from (34) is that the stream function and the generalized reduced pressure are the same in the limit of small magnetic Ekman and Rossby numbers. This type of magnetostrophic flow (in analogy to the geostrophic type in rotating flows) means thar isobars represent streamlines at the leading order for and , i.e., for .
5. Parameter Estimation and Suggested Experimental Setup
In this section an evaluation of the conditions for electromagnetic columns and magnetostrophic flow to exist are tested via parameter estimation and a suggested experimental setup is proposed. The objective is in establishing the criteria for the theoretical results to be detectable in practice. The main condition for the emergence of electromagnetic columns and magnetostrophic flow was derived as
For the flow of free charges one can estimate the value of
as the ratio between the mass of the elementary charge (the electron) and the absolute value of its electrical charge. For flow of charges via fluids one has to account for the mass of the fluid too when evaluating the mass, however for rarified gases the former estimation still applies as an approximation. Then for an imposed magnetic field of
the ratio and condition (35) implies
. This means that for length scales of
,
, and
, the condition that the average fluid velocity should satisfy is
,
, and
, respectively. The first two values exceed the speed of light in vacuum. These conditions are easily satisfied in any realistic experimental setups. Even if the value of
turns out to be much larger even by a few orders of magnitude it should not be difficult to satisfy this condition, which on extreme cases will require stronger magnetic fields. Estimates for the secondary conditions, i.e.,
(or
),
(or
), (or
) indicate that they can also be satisfied in reality. For example, for plasmas large velocities are required such that for length scales of
,
, and
, the condition that the average fluid velocity should satisfy when combined with the condition for
is
,
, and
, respectively. For liquid metals the conditions are much less constrained as
,
, and
. Therefore experiments to confirm the emergence of electromagnetic columns and magnetostrophic flow are being suggested and planned by using the constraints evaluated in this section and an experimental setup as described in
Figure 1 and
Figure 2. Using earth magnets for creating the imposed magnetic field is also quite practical given the estimated parameters values.
6. Conclusions
The theoretical derivation identifying an analogy between magneto-fluid dynamics (MFD/MHD) and geostrophic flow in a rotating frame of reference was presented. The latter includes the existence of electromagnetic columns identical to Taylor-Proudman columns. The emergence of these columns is conditional upon very small magnetic Rossby numbers or alternatively very large values of its reciprocal, the electromagnetic number. This condition was shown to be possible to satisfy in reality.
As a result, the electromagnetic fluid flow is two dimensional and the streamlines are being shown to be identical to the pressure lines in complete analogy to rotating geostrophic flows. A possible setup is suggested to confirm the theoretical results experimentally in
Figures 1, 2, and in detail in
Figure 3. Using earth magnets for creating the imposed magnetic field is also quite practical given the estimated parameters values. An imposed magnetic field is anticipated to be affecting a stream of beams of charges or MFD type fluids that is flowing through. A specific type of detection via a tracer, PIV, or optical methods such as interferometry or similar can be used for the detection of the flow of charges. When magnetostrophic flow and magnetic virtual columns are detected, it raises the possibility of the creation of von-Kármán vortex streets as presented in
Figure 3 on the “top view”. This is similar to the creation of von-Kármán vortex streets in the wake of Taylor-Proudman columns in rotating flows as captured by satellite images and presented theoretically by Vadasz [
6].
References
- H.P. Greenspan, Theory of Rotating Fluids, Cambridge University Press, (1968).
- Vadasz, P. Centrifugal buoyancy in a rotating fluid layer next to and distant from the rotation axis. Phys. Fluids 2021, 33, 034123. [Google Scholar] [CrossRef]
- Amar, S.D.; Rashkovan, A.; Ziskind, G. Analysis of a stable bathtub vortex in a rotating container. Phys. Fluids 2023, 35. [Google Scholar] [CrossRef]
- Subbotin, S.; Shmakova, N.; Kozlov, V.; Ermanyuk, E. Nonlinear regimes of inertial wave attractors generated by a precessing lid: Zonal flows and Rossby waves. Phys. Fluids 2023, 35. [Google Scholar] [CrossRef]
- Burmann, F.; Noir, J. Effects of bottom topography on the spin-up in a cylinder. [CrossRef]
- P. Vadasz, “von Kármán vortex streets around invisible bluff bodies”, Physics of Fluids, Vol. 35 (9), pp. 096604 (1-5), (2023). [CrossRef]
- Sarkar, S.; Sahoo, B.; Sekhar, T.V.S. Influence of magnetic field in the control of Taylor column phenomenon in the translation of a sphere in a rotating fluid. Phys. Fluids 2021, 33, 073606. [Google Scholar] [CrossRef]
- Vadasz, P. Stability of free convection in a rotating porous layer distant from the axis of rotation. Transp. Porous Media 1996, 23, 153–173. [Google Scholar] [CrossRef]
- Vadasz, P. Coriolis effect on gravity-driven convection in a rotating porous layer heated from below. J. Fluid Mech. 1998, 376, 351–375. [Google Scholar] [CrossRef]
- P. Vadasz, “Fluid Flow and Thermal Convection in Rotating Porous Media”, Handbook of Porous Media, Marcel Dekker, New York, pp. 395-439, (2000).
- P. Vadasz, Fluid Flow and Heat Transfer in Rotating Porous Media, Springer (Springer Briefs in applied Science and Engineering), ISBN: 978-3-319-20055-2 (Print), 978-3-319-20056-9 (Online), (2016).
- P. Vadasz, “Natural Convection in Rotating Flows”, Handbook of Thermal Science and Engineering, F. A. Kulacki (ed.), Chapter 11, Springer International Publishing AG, pp. 691-758, (2018).
- Vadasz, P. Instability and Convection in Rotating Porous Media: A Review. Fluids 2019, 4, 147. [Google Scholar] [CrossRef]
- Vadasz, P.; Govender, S. Stability and stationary convection induced by gravity and centrifugal forces in a rotating porous layer distant from the axis of rotation. Int. J. Eng. Sci. 2001, 39, 715–732. [Google Scholar] [CrossRef]
- S. Saravanan and S. Vigneshwaran “Centrifugal filtration convection in bidisperse media”, Physics of Fluids, Vol. 32 (8), pp. 084109 (1-8), (2020).
- Kang, J.; Chen, X.; Fu, C.; Tan, W. Centrifugally driven thermal convection in a rotating porous cylindrical annulus. Phys. Fluids 2013, 25, 044104. [Google Scholar] [CrossRef]
- E Allen, J.; Auer, P.L.; Endean, V.G. On the law of isorotation and laboratory experiments. Plasma Phys. 1976, 18, 143–154. [Google Scholar] [CrossRef]
- Opat, G.I. Coriolis and magnetic forces: The gyrocompass and magnetic compass as analogs. Am. J. Phys. 1990, 58, 1173–1176. [Google Scholar] [CrossRef]
|
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. |
© 2024 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).