Next Article in Journal
Distributed Optimal Placement Generators in a Medium Voltage Redial Feeder
Next Article in Special Issue
Numerical and Computational Analysis of Magnetohydrodynamics over an Inclined Plate Induced by Nanofluid with Newtonian Heating via Fractional Approach
Previous Article in Journal
Two-Dimensional Turbulent Thermal Free Jet: Conservation Laws, Associated Lie Symmetry and Invariant Solutions
Previous Article in Special Issue
Theoretical Survey of Time-Dependent Micropolar Nanofluid Flow over a Linear Curved Stretching Surface
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Model of Water Treatment by Nanoparticles in a Channel with Adjustable Width under a Magnetic Field

1
Department of Computer Software and Automated Sysytem, Belgorod State Technological University Named after V.G. Shukhov, 308012 Belgorod, Russia
2
Department of Mathematics, COMSATS University Islamabad, Abbottabad Campus, Abbottabad 22060, Pakistan
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(8), 1728; https://doi.org/10.3390/sym14081728
Submission received: 25 July 2022 / Revised: 7 August 2022 / Accepted: 9 August 2022 / Published: 18 August 2022
(This article belongs to the Special Issue Symmetrical Mathematical Computation in Fluid Dynamics)

Abstract

:
The process of water treatment by nanoparticles is one of the most considerable subjects in the cross-field of hydrodynamics, chemistry, and mathematics. This paper is dedicated to the case of the flows that appear when squeezing and stretching a channel with mixing of water, nanoparticles, and contaminants. It is assumed that fluid is homogeneous at the starting moment, the parameters of the nanoparticles and contaminants are known, and there is a constant non-homogeneous magnetic field applied to the system. The flow starts moving when the walls of the channel shift to each other. Exact and numerical solutions of the system of ordinary differential equations are used to obtain the results. The article gives an answer to the question about stability of the flow and proposes the technique to evaluate the essential characteristics of the system to achieve the treatment process efficiency. The main result is that the considered system shows excellent properties concerning purification of water on the selected part of the squeezing stage. This effect does not appear without a magnetic field. The mentioned properties are: decreasing of nanoparticle concentration to zero inside of the unsteady layer under magnetic field close to 1 T and this effect stays until the channel become about 10% of initial width as a result of the squeezing.

1. Introduction

Magneto-Hydro-Dynamics (MHD) is the crossroad of physical and mathematical methods which are used to analyze the motion of fluids in magnetic or electro-magnetic fields. The triple word magneto-hydro-dynamics consists of “magneto” meaning magnetic, “hydro” meaning liquid, and “dynamics” referring to the object movement due to forces. Sometimes MHD is named as “hydro-magnetics” or “magneto-fluid-dynamics”.
Hannes Alfven [1], a famous Swedish research analyst, first explored the MHD fluid. His contribution in the field of plasma physics is remarkable. The magnetic field is an electric current in the flowing conducting fluid. The electric current flowing through the conductive fluid generates a force on the fluid and affects the magnetic field. The electrically conductive flows were analyzed by numerical simulation in various applications, such as in parallel film [2,3], applied magnetic field [4], and treatment of polymers in microgravity [5].
There are two problems among the whole specter of MHD interests that concern the present paper’s subject. The first one is common enough and it is called the stability problem. It is targeted to answer the question “what are the conditions for the flow to become stable?” Of course, the magnetic field influence is taken into account first of all. The second problem is particles’ association both for nanoparticles and for contaminants. It attempts to answer the question “how do magnetic field and flow velocity impact the scalar fields of concentrations?”.
In the present study the main goal is to arrange the zone in the flow where the concentration of nanoparticles decreases extremely. This zone can be treated as the source of clean water, because the contaminants have been already adsorbed by nanoparticles and if they are being removed, then contaminants are being removed too. The stability of the flow is an essential factor in this concern and it takes attention into consideration.
The latest references concerning instability explorations can be illustrated by [2,3,4,5]. Stability of two fluids’ MHD flow inside a flat channel was explored by Z. Hussain et al. in [5] where the corresponding exact solution was presented. Of course, there are a lot of research articles and conferences papers devoted to this problem, but we only mention the work [6] of A.S. Dalkılıç et al., where the heat conduction of hybrid nanofluid CNT-SiO2 with water was analyzed experimentally; these results are closed to those in the present paper, in spite of the difference in physical background. Numerical simulations of the MHD flow stability in a rectangular channel were made in [7]. Stability and some other physical aspects have been discussed in [8] with respect to the problem of water treatment.
The influence of certain functioning parameters is inspected, and notable results in [9] are found that the rate of heat transfer is exaggerated along with the skin friction coefficient while the suction/injection and magnetic parameters are intensified. The results also signified that the rise in the volume fraction of the nanoparticle and the decline in the unsteadiness parameter demonstrates a downward attribution towards the heat transfer performance and skin friction coefficient. The skin friction coefficient intensifies in conjunction with the local Nusselt number by enhancing the suction/injection parameter and the nanoparticle volume fraction past an exponentially stretching/shrinking sheet in [10]. A vertical magnetic field is incorporated in the flow while the energy equation is modeled in [11] taking the Dufour and Soret effects as heat generation through viscous dissipation, and ohmic heating plays a vital role in designing various engineering devices such as electric stoves, heaters, thermistors, food processing, etc. It is observed in [12] that heat and mass transfer considerably enhance along vertical disk movement. Additionally, magnetic field, temperature ratio, and radiation parameter significantly enhance the temperature field, while reaction rate parameter and Schmidt number decrease the concentration profile. The hydrodynamic stability of the Couette flow of an electrically conducting fluid flowing through a porous medium in a parallel-plate channel with a normal magnetic field is investigated in [13]. Our results are devoted mostly to distribution of instable areas in space unlike the results of [13] where the influences of speed, viscosity, and magnetic field on the flow stability have been explored. These results show that the considered problem was explored in different physical cases, but not for the case described above for this study.
The second problem appointed above is under consideration in [14] where the polymerization problem has been explored. The concentrations are closely related to heat transfer. That is why most methods that have been used for heat transfer explorations are applicable to concentration analysis as well. Thus, Jing et al. [15] studied the heat transfer with the help of five distinct sorts of nanoparticles in unsteady flow of nanofluids.
The results of a magnetic field influence on the flow of heat transmission of Fe2O3 with a uniform heat flux was examined by Kandelousi [16]. He examined that the flow of nanofluids and heat transmission across the absorbent medium.
Nanofluids have many applications in different industries, for example, in the automotive industry it is used for enhancing the lubricant and liquid coolants in the systems of cooling-air radiators. That is why a wide range of methods are used for heat-transmission problem solving. In 2016, A. Hussanan et al. [17] solved the particular heat transfer problem by using analytical technique.
The present paper focuses mostly on the case of water treatment, especially for industry applications. Because of this reason, we mostly deal with concentrations instead of heat. However, we use similar methods, of course. The enhancement of heat transmission basically depends on the thermal conductivities of nanoparticles, the volume concentration of particles, and the mass flow rate. Almost the same parameters are under consideration in the case of the present paper.
It is tempting to make water treatment process cycling independent from any filter: the process become more efficient and reliable. Nanoparticles are considered as good adsorbers of the contaminants, but the question about how to remove nanoparticles from the solution, the suitable answer is not found yet. The present study shows the way to remove nanoparticles by means of mechanical and magnetic operations. The squeezing of the channel is considered particularly because it is the stage of nanoparticle removal. Although some results are related to the stretching stage as well, the graphs and conclusions are made concerning squeezing only. We can highlight the next research stages:
(1)
To build the model for the flow of viscid nanofluid with contaminant in the non-homogeneous magnetic field;
(2)
Analytical analysis of the model from the previous stage;
(3)
Numerical experiment for various parameters of the model from the previous stage.

2. Materials and Methods

The schematic outlook of the channel is shown in the Figure 1. A semi-infinite rectangular channel is considered as a model of a narrow channel inside a sponge, and it is supposed that the channel is symmetric, and the solution can be reflected to another half-plane. The walls of the channel are partially penetrable, and they are shifting with constant velocity a ˙ = Ω . The shift may be both stressing ( Ω < 0 ) and stretching ( Ω > 0 ). In the experiment the walls squeezed at a uniform rate Ω < 0 , and porosity value is V 0 , 1 . The fluid of the flow is a viscid nanofluid of (Fe3O4/H2O). Although for this study any nanoparticles with nonzero magnetic moment can be considered, the Fe3O4 nanoparticles are chosen due to their accessibility and low price. In the flow there is a contaminant which has known mass and concentration values in standard laboratory IUPAC conditions.
The assumptions of the model are formulated in accordance with the well-known physical chemistry theory (see, for example, [18]):
  • Before the flow starts to move, there is homogeneous fluid inside the channel and nanoparticles are not aggregated and not adsorbed; this can be physically made by a sonication process;
  • The start of the considered process coincides with the moment when the sonicator turns off, the squeezing starts, and the magnetic field is applied;
  • It is supposed that the further mentioned values had been gained in other experiments and they are known:
    Any nanoparticle ( ν ) has the initial aggregation ability 0 A a 0 1 value, that shows that which part of the aggregable particles will be really aggregated, and the initial adsorption ability 0 A s 0 1 value, that shows what part of the adsorbable particles will be really adsorbed; the values are functions of time and position in general, but they do not depend on position at the initial time;
    Adsorption of the nanoparticles ( ν ) means that a particle of contaminant ( c ) being inside the r s -ball around a fixed nanoparticle, will be attached to the nanoparticle during the time τ s , and the mass of the nanoparticle will be increased by the mass of the attached particle; at the same time, the aggregation ability A a decreases by k s a times, adsorption ability A s decreases by k s s times;
    Aggregation of the nanoparticles ( ν ) means that a nanoparticle being inside the r a -ball around a fixed nanoparticle will be attached to the nanoparticle during the time τ a , and the fixed nanoparticle mass will be increased by the mass of the attached particle; at the same time the aggregation ability A a decreases by k a a times, adsorption ability A s decreases by k a s times.
Let the fluid have volume density at the point x , y equal to
ρ = m w   n w + m c   n c + m ν   n ν ,
where m is the particle’s mass, n is concentration in the fluid and the indices denote as water (w), contaminant (c) and nanoparticles ( ν ). We denote as δ c = n ¯ w n ¯ c   , δ ν = n ¯ w n ¯ ν the concentration factors for contaminant and nanoparticles, where overlined values are related to the standard laboratory conditions.
In order to make the expressions more readable, we use the following notations for the constants and functions. New constants are defined as follows
α s = 4 π r s 3 3 τ s ,   α a = 4 π r a 3 3 τ a ,
c s s = 1 k s s τ s k s s , c a s = 1 k a s τ s , c a a = 1 k a a τ a k a a , c s a = 1 k s a τ a .
Additionally, new functions will be next
A s = A s n c , A a = A a n ν .
The abovementioned propositions concerning the model imply the following differential equations
d m ν d t = t m ν + u x m ν + v y m ν = α s m c n c 2 A s + α a n ν 2 m ν A a ,
d n c d t = t n c + u x n c + v y n c = α s n c 2 n ν A s ,
d n ν d t = t n ν + u x n ν + v y n ν = α a n ν 3 A a ,
d A s d t = t A s + u x A s + v y A s = c s s A s + c a s A a ,
d A a d t = t A a + u x A a + v y A a = c a a A a + c s a A s .
We have continuity and Navier–Stokes equations additionally (Landau and Lifshitz, [19]):
ρ ˙ + ρ x u + ρ y v + u x ρ + v y ρ = 0 ,
ρ t u + u x u + v y u = x p + η x x u + y y u ζ + η 3 x x u + x y v + n ν f ^ x ,
ρ t v + u x v + v y v = y p + η x x v + y y v ζ + η 3 x y u + y y v + n ν f ^ y ,
where
ρ = ρ w + m c m w δ c n c + m ν m w δ ν n ν .
The auxiliary conditions are as follows
u = 0 , v = V Ω   for   y = a ,   t 0 ,   x ,
u = 0 , v = V Ω   for   y = a ,   t 0 ,   x ,
u = v = 0   for   y = ± a ,   t 0 ,   x .
Additionally, there are initial values of the variables at the start of the process:
n ν 0 = n ν | t = 0 , n c 0 = n c | t = 0 , m ν 0 = m ν | t = 0 , A a 0 = A a | t = 0 , A s 0 = A s | t = 0 .
The above relations are not depending on x , y for the consideration of this paper.
Two Equations (11) and (12) can be formulated as one equation with p function eliminated. Therefore, we have seven equations totally: Equations (5)–(10) plus one Equation made from Equations (11) and (12). There are seven dependent variables ( u , v , n ν , n c , m ν , A a   and   A s ) and three independent variables ( x , y , t ). The system has total order equal to 9 and we have a relevant number of boundary conditions because the auxiliary conditions operate with values which are uniform in space.
According to the Figure 1 and the model assumptions, the magnetic field components which act on the magnetic momentum m and mass m ν can be determined from the known laws of electromagnetism ([20] and other textbooks). Thus, the magnetic field of the thin solenoid at the point r = x , y , z is given by the Biot and Savart law
d B = μ 0 I 4 π d l × r r 3 .
Additionally, (from the same source) the magnetic force that acts on the magnetic momentum m oriented along the external magnetic field B , is determined by
f = m B x x + B y x , B x y + B y y ,
where the third coordinate is dropped. Direct computations from these and a new notation of y y + a (due to the configuration shown in Figure 1) give the following components of the forces:
f ^ x = 3 H m 6 x 3 y + a 2 x 5 5 x y + a 4 ( x 2 + y + a 2 ) 7 2 ( x 2 + 2 y + a 2 ) 2 + x 2 y + a 2 , f ^ y = 12 H m y + a 3 ( x 2 + y + a 2 ) 5 2 ( x 2 + 2 y + a 2 ) 2 + x 2 ( y + a ) 2 ,
where H is the relevant constant, and it is easy to see that
y f ^ x x f ^ y = 0 .
It is clear that constant homogeneous magnetic field does not force the nanoparticles to move. That is why we need a non-homogeneous magnetic field to attract nanoparticles to the bottom of the channel. The corresponding vector field configurations for the initial point and another three positions during the squeezing are shown in Figure 2, where it should be noted that the lengths of arrows are not in the same scale for all shots. Actually, at the last stage the magnetic field is the strongest, but at the same time its heterogeneity is the highest at the moment and it acts farther along the x axis.
The configuration of the magnetic field, shown in Figure 2, looks like a usual magnetic field of the strip magnet. This is confirmation of the formulas (18) because the graphs are made from them.

3. Results

3.1. The Problem Solving

First of all, we should note that the only time factor in the system is the wall moving rate Ω . Hence, it is natural to unify the time with the y coordinate using the new independent variable
ξ = y a ,   a = Ω t + a 0 ,
which is unitless and runs over 1 , 1 . The unification can be made while the rate Ω is constant in the model and we can write as
y = 1 a ξ ,   t = ξ Ω a ξ .
Let us try to eliminate the second variable in order to transform the problem to the dynamical system formulation. If we put
t + u x + v y z .
For some (in general, complex) variable z , then we will have
u x + v ξ Ω a ξ = z ,   ξ = a v ξ Ω z u x .
Additionally, it follows from Equation (5) that:
z m ν = α s m c n c 2 A s + α a n ν 2 m ν A a .
And so on. We have the next transformation rules
x , ξ x , z , x x = 1 , x ξ = 0 , z x = u , z ξ = v ξ Ω a , x x = 1 , x z = 0 , ξ x = a u v ξ Ω , ξ z = a v ξ Ω .
In order to gain the ODE system, we should find expressions in terms of z ,   x variables for all derivatives involved in Equations (11) and (12), and it is necessary to make an identity from the Equation (10).
First of all, let us note that in coordinates x , ξ we obviously have the next:
x u = 0 , ξ u = a v ξ Ω z u , x v = 0 , ξ v = a v ξ Ω z v .
Additionally, it is possible to evaluate v explicitly from Equation (10) that transforms to the following one:
z v v ξ Ω = z ρ ρ .
From here, we have
ρ 0 v 0 ρ e Ω a z + ξ Ω = v ,
With some ρ 0 x v 0 . Taking z from here and substituting it into Equation (25), we find
z = a Ω ln ρ 0 v 0 ρ v ξ Ω ,
ξ z = a v ξ Ω = a Ω ρ v ξ Ω ξ ρ ρ 2 v ξ Ω ξ v Ω ρ v ξ Ω 2 ,
1 = ξ ρ v ξ Ω Ω ρ ξ v Ω + 1 ,
ξ ρ ρ = ξ v v ξ Ω .
Equation (32) combined with Equation (27) shows that ρ 0 v 0 ρ does not depend on x . Integration of the Equation (27) gives the solution for v :
v = ξ Ω + 1 V Ω e i π 1 + ξ 2 ,
It obviously obeys the conditions on v in the real part. From here, we have
ρ 0 v 0 ρ R z = 1 V Ω e Ω a z + i π ξ + 1 2 .
Let us now solve Equations (5)–(9) which are equations on the only independent variable z :
z m ν = α s m c n c 2 A s + α a n ν 2 m ν A a ,
z n c = α s n c 2 n ν A s ,
z n ν = α a n ν 3 A a ,
z A s = c s s A s + c a s A a ,
z A a = c a a A a + c s a A s .
The solution can be found easily, and it has the following form:
n ν z = 2 α a 0 z A a μ d μ + 1 n ν 0 2 1 2 ,
n c z = 1 α s 1 α s n c 0 + 0 z A s μ n ν μ d μ 1 ,
m ν z = m c n c z n c 0 n ν z + n ν 0 m ν 0 n ν z ,
A s = F 1 e b 1 z + F 2 e b 2 z ,
A a = F 1 b 1 c s s c a s e b 1 z + F 2 b 2 c s s c a s e b 2 z ,
where
b 1 , 2 = c a a + c s s ± c a a + c s s 2 4 c a a c s s c a s c s a 2 .
Additionally, in general,
F 1 = A a 0 c a s A s 0 b 2 c s s c a a + c s s 2 4 c a a c s s c a s c s a , F 2 = A a 0 c a s + A s 0 b 1 c s s c a a + c s s 2 4 c a a c s s c a s c s a .
Index 0 corresponds to the z = 0 position.
In consequence of Equations (11) and (12), by equivalence of mixed derivatives, we derive:
y ρ t u + u x u + v y u ρ ( t y u + y u x u + y v y u + u x y u + v y y u t x v x u x v x v y v u x x v v y x v ) + η x x y u + y y y u x x x v x y y v ζ + η 3 x x y u + x y y v x x y u x y y v + x ρ t v + u x v + v y v + f ^ x y n ν f ^ y x n ν = 0 .
This complicated equation can be essentially simplified using our previous considerations. Thus, involving Equations (25)–(27), we have from Equation (47)
z ρ ρ + Ω a z u z z u + η ρ ρ 0 2 v 0 2 e 2 Ω a z ( 4 z ρ ρ Ω a + 2 Ω 2 a 2 + z z ρ ρ + z ρ 2 ρ 2 z u + 3 z ρ ρ + Ω a z z u + z z z u ) + f ^ x ρ z n ν f ^ y ρ 0 v 0 e Ω a z + f ^ x ρ u x n ν = 0 .
Let us look for solutions in the form of
u = U R e Ω a z .
For some constant U and note that
x n ν = 0 .
Then, from Equation (48) we derive
i π R 2 a e Ω a z u u + η ρ ρ 0 2 e 2 Ω a z π 2 4 a 2 R 2 e 2 Ω a z u + 3 i π R 2 a e Ω a z u + u + f ^ x ρ z n ν = 0 .
While we have from Equation (34)
R = R Ω a + i π 2 a R 2 e Ω a z , R = R Ω 2 a 2 + i π a Ω a R 2 e Ω a z π 2 2 a 2 R 3 e 2 Ω a z .
Then,
u = U i π 2 a R 2 e 2 Ω a z ,
u = U π 2 2 a 2 R 3 e 3 Ω a z ,
z ln R = Ω a + i π 2 a R e Ω a z ,
z z ln R = π 2 4 a 2 R 2 e 2 Ω a z .
Substitution into the Equation (51) shows that the function (49) turns to zero for all terms in Equation (51) except the last one. In order to solve the equation, let us use the variation of the constant method. Suppose that U = U z . Then,
u = U i π 2 a R 2 e 2 Ω a z + U R e Ω a z ,
u = U π 2 2 a 2 R 3 e 3 Ω a z + U i π a R 2 e 2 Ω a z + U R e Ω a z .
Additionally, we have from Equation (51) that
U ρ 0 v 0 e Ω a z + i π η 2 a U η R e Ω a z f ^ x z n ν = 0 .
This equation can help us to solve the first problem: to find flow stability conditions.

3.2. Stability of the Flow

In order to find the stability conditions, we need to formulate the very definition of the stability. The stability breaks when velocities have the chance to become infinite and the critical condition is
u = i u ^ i e ω c i t , max i R e   ω c i = 0 ,
where u ^ i is allowed to be at most polynomial in time. The flow is stable when all R e   ω i < 0 and it is unstable when R e   ω i > 0 for at least one value of i . Of course, the same is true for the other velocity component v . Let us recall the first equation of Equation (52) and solve it:
R = e Ω a z 1 v 0 i π 2 a z 1 ,
where v 0 = R 0 . It is easy to derive that
v 0 = i 1 V Ω e i π ξ 0 2 ,     z = 2 a π 1 V Ω e i π ξ 2 e i π ξ 0 2 ,
where ξ 0 is value of ξ that corresponds to z = 0 .
Let us consider the vertical velocity v . Taking into account that z is time equivalent, we have the following stability condition from Equations (28) and (60):
R e   Ω a z < 0 .
Recalling that Ω < 0 when squeezing and Ω > 0 when stretching, we have from Equations (62) and (63) that the flow is stable when:
s q u e e z i n g :   cos π ξ 2 < cos π ξ 0 2 , s t r a t c h i n g :   cos π ξ 2 cos π ξ 0 2 .
In particular, if there is no magnetic field applied, just ξ 0 = 0 has a physical sense and hence we have very simple conditions:
  • The flow has no stability restrictions from the v component when squeezing;
  • The flow is always unstable when stretching.
However, in the presence of the magnetic field, the value of ξ 0 may be non-zero and we will have something like the following (the restrictions derived only from the v component)
  • The flow is stable for ξ > ξ 0 and unstable for ξ < ξ 0 when squeezing;
  • The flow is stable for ξ < ξ 0 and unstable for ξ > ξ 0 when stretching.
Now, let us turn to the velocity component u and consider the flow under the magnetic field influence. From Equations (34), (40), (49) and (59) we can find that all possible exponents in the expression of u are the following:
B Ω a z ,   B z , B + Ω a z , 1 2 b 1 z , b 1 3 2 b 2 z , b 2 3 2 b 1 z , 1 2 b 2 z ,
where B is (real-valued) an exponential factor of u f ^ x in “time” z . All values in Equation (65) must have negative real parts to provide flow stability. It follows from here and from conditions for v , that for squeezing the conditions take the form
z r < 0 ,   B z i < Ω a z r ,   0 > 2 3 b 2 > b 1 > 3 2 b 2 .
In the case of stretching B = max i B i and inequalities are the following
z r > 0 , B z i < Ω a z r ,   0 < 2 3 b 2 < b 1 < 3 2 b 2 .
Using Equation (18), the value of B is evaluated as
B = 5 1 V Ω a ξ 0 + 1 .
Thus, the MHD-flow is stable when the following conditions hold:
  • Squeezing ( Ω < 0 ): ξ > ξ 0 ,   0 > 2 3 b 2 > b 1 > 3 2 b 2 ;
  • Stretching ( Ω > 0 ): ξ < ξ 0 ,   0 < 2 3 b 2 < b 1 < 3 2 b 2 .
With the common condition of
5 1 V ξ 0 + 1 < cos π ξ 2 cos π ξ 0 2 sin π ξ 0 2 sin π ξ 2 .

3.3. The Application of the Model

In order to illustrate our model, the colloid fluid based on water with calcium hydro carbonate of 20 (that is slightly hard water) is considered and Fe3O4 nanoparticle characteristics were used in the model. The following parameters correspond to the chosen fluid:
r s = 2.1013 10 8   m ,   r a = 1.41842 10 8   m ,   τ s = 4.2026 10 2   s ,   τ a = 2.83684 10 2   s ,   k a a = 1.1 ,   k a s = 10.2 ,   k s a = 10.5 ,   k s s = 1.6 .
The values of parameters are evaluated from reference values [21].
We built graphs for n ν ξ for fixed a 0 = 0.1   m and Ω = 0.005   m / s under magnetic field made by a permanent alnico magnet that gives induction of about 0.6 and 0.9 T in the form as shown in Figure 2.
It is possible to determine the distribution of nanoparticle concentrations using the Equations (40)–(44) and definition of z ξ from Equation (62). The formula is too large, and we illustrate it using a graph.
In Figure 3, the situation without magnetic field is shown for the cases of standard nanoparticle concentrations (green color), 20% reduced nanoparticle concentrations (blue color), and 20% increased nanoparticle concentrations (red color). It is clear that concentration oscillates inside the unstable zone and takes a constant value in the stable zone. At the same time, it is notable that the amplitude of the oscillations grows together with nanoparticle concentrations: in particular, the relative concentration decreases outside the turbulent zone. However, the form of the oscillations does not change.
The result for a magnetic field of 0.6 T impact is shown in Figure 4 and Figure 5 showing the case of 0.9 T. It can be noted that the magnetic field shifts down the maximum of the (c) curve. Another notable point is that concentration decreases almost to zero inside the turbulent zone under a strong enough magnetic field. This result is new.
It should be noted that all presented curves were constructed using numeric evaluation, and they are just approximations. However, the underlined effects are confirmed for any large number of steps.
As it is clear from Figure 3, Figure 4 and Figure 5, the presence of a magnetic field does not affect the concentration of nanoparticles in the laminar layers of the flow. However, it has an extremely high influence on the concentration inside the turbulent layer: it decreases almost to zero on the squeezing stroke whenever the magnetic field is strong enough. It is possible to show the influence of the magnetic field on the nanoparticle concentrations inside the turbulent layer with an average concentration versus field magnitude graph, shown in Figure 6.
It may be helpful to see what happens with the nanoparticle concentrations when the squeezing rate changes. Let us show the concentrations versus vertical position for the cases of the rate decreasing and increasing twice. Please consult Figure 7 and Figure 8, where only standard concentrations are considered.
The graphs show that the low squeezing rate gives a better result: clearness in the turbulent layer stays for an increased duration. However, it is clear that this dependency is nonlinear and the best rate should be determined specially.

3.4. The Physical Interpretation and Model Proofs

The concentration decreasing effect, shown above, has a simple physical background. To demonstrate it, let us discuss the mathematical model predictions first.
The prediction, does not depend on the presence of nanoparticles, is instability appearance: the flow can be stable under squeezing, but it always unstable under stretching. Anybody can see it using the ordinary French press (for brewing tea). This phenomenon, as far as the authors know, has not been previously deduced from the laws of hydrodynamics. However, the mathematical model supports it. This is evidence in favor of the correctness of another prediction of the model: nanoparticle concentration decreases in the zone of unstable flow under a magnetic field at the beginning of the flow squeezing.
The physical sense of the shown effects is next. Instability of the flow causes the nanoparticles to move faster, i.e., to increase their energy. This leads to oscillations of the concentration inside the instable zone without an external magnetic field. If the nanoparticles have a magnetic moment (Fe3O4 has it), then they choose the special direction to move under an external non-homogeneous magnetic field and if the magnitude of the field is strong enough, the particles tend to escape the unstable flow zone to the area where they can move slower. This will not occur under stretching because there are no stable zones.

4. Discussion

The results received in the present study are in agreement with other research and provide additional impact in the specific scientific area.
The study of squeezing flow of Cu water and Cu kerosene under the influence of a magnetic field with a pressure gradient using the regular perturbation method [22] showed that nanoparticle concentration has a direct impact on the flow velocity profile during the squeezing stage (see Figure 6a,b in [22]). This is in agreement with our result: the flow velocity has a direct impact on the nanoparticle concentration (according to Figure 8, the half-rated velocity leads to a much lower concentration compared to that in Figure 7 with double-rated velocity).
In [23], the squeezed MHD flow of water-based metallic nanoparticles over a porous sensor surface in the presence of a heat source was investigated. Although we considered no heat source in our study, the behavior of the characteristic value (in [23] it was temperature) depends on the compression stage in the same way as concentration in our experiment depends on the stage of squeezing.
The study [24] was dedicated to unsteady squeezing flow of nanofluids. Among the results of the paper is the following: the velocity profile provides the oscillatory behavior with the enhancement of the velocity and magnetic magnitude. This is in agreement with our results, although our magnitudes of the field are much greater.
Our study provides the following consequences based on the graphs above. In Figure 3, the flow without a magnetic field is presented. When the channel is squeezed, the nanoparticle concentrations become oscillating inside the unsteady layer and these oscillations are attenuated when the absolute width of the channel decreases. Figure 4 shows the same process, but under a non-homogeneous magnetic field of 0.6 T magnitude and it is clear that the oscillations become irregular. Figure 5 shows the process for a 0.9 T magnetic field and this case is extremely different from those before: the concentration inside the unsteady layer is dropped to zero and stays until the channel becomes narrow; Figure 5c corresponds to the channel of 0.01 m that is 10% of the initial width.
Figure 6 illustrates how the average nanoparticle concentration inside the unsteady layer at the moment related to the half-squeezed channel depends on the magnetic field strength. The strength of the field leads to a decrease in average nanoparticle concentration; hence, it enforces purification ability of the system (up to 90% with a 0.9 T field).
Figure 7 provides a description for the same processes as shown in Figure 3, Figure 4 and Figure 5 for a standard concentration (green-colored line), but with a double squeezing rate. In Figure 7, the lines of different colors correspond to different widths of the channel, however, the same width as before (95%, 50%, and 10% of initial width). The effect observed before becomes weaker essentially: increasing the squeezing rate is unjustifiable.
Figure 8 corresponds to the half squeezing rate (compared to Figure 3, Figure 4 and Figure 5). The required effect is observed and looks the same as it was with the standard rate. Thus, too fast squeezing destroys the purification ability, but slower squeezing preserves the purification properties, but decreases the efficiency (due to time lost).

5. Conclusions

The following highlights are related to this work:
1.
It was a new concept proposed to consider the flows of Poiseuille and Couette with multicomponent liquids: using the complex independent variable (with a time-like real part) that considers the inner time of the flow;
2.
Stability conditions of the nano-liquid flow were formulated both for squeezing and for stretching of the flow;
3.
It was shown that there are oscillations of nanoparticle concentrations across the channel when the magnetic field is absent. These oscillations attenuate with time;
4.
It was shown that nanoparticle concentrations become irregular oscillation under a weak magnetic field. The concentration in the turbulent part of the flow becomes almost zero in a strong magnetic field.
There are several physical consequences of the present study:
-
Any flow configuration shown in Figure 1 becomes partially unsteady when the channel is squeezed: the unstable flow occupies the central part of the channel, and stable layers are close to the walls.
-
When the channel is stretching, the whole flow is unsteady and there is no stable layer.
-
If there is nanofluid in the flow, then the nanoparticle concentration inside the unsteady zone can be extremely (practically down to zero) decreased if a magnetic field of enough magnitude (about 1 T) is applied. This is the purification effect that was the goal of the paper.
-
The value of the squeezing rate is essential for the presence of the purification effect: too fast squeezing destroys the purification abilities even if the magnetic field is strong enough.

Author Contributions

Conceptualization, Z.H. and S.Z.; methodology, S.Z.; software, S.Z.; validation, Z.H., S.Z. and P.K.; formal analysis, S.Z.; investigation, S.Z.; resources, Z.H.; data curation, P.K.; writing—original draft preparation, S.Z.; writing—review and editing, Z.H.; visualization, S.Z.; supervision, S.Z.; project administration, S.Z.; funding acquisition, S.Z. and P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Symbols

NotationMeaningUnit of Measurement
Ω Velocity of walls movingm/s
a ,   a 0 Instant distant between walls, initial distancem
V Porositydimensionless
x Horizontal coordinate along channelm
y Vertical coordinate across channelm
t Times
u Horizontal velocity of the flowm/s
v ,   v l , v u Vertical velocity of the flow inside the channel, lower than channel and upper than channelm/s
α a , α s The values of enhancement of active sphere for aggregation and adsorption, respectivelym3/s
c s s , c s a , c a s , c a a The characteristics of the changes: adsorption by adsorption, aggregation by adsorption, adsorption by aggregation and aggregation by aggregation, respectivelys−1
A s , A a The functions (on t ) of inverse concentration of contaminants multiplied by adsorption ability and inverse concentration of nanoparticle multiplied by aggregation abilitym3
m w   , m c , m ν The masses of the particles of: water, contaminants, and nanoparticles, respectivelykg
n w ,   n c ,   n ν The concentrations of the: water, contaminants, and nanoparticles, respectivelym−3
ρ , ρ w The density of the mixture of water, contaminants, and nanoparticles (instant local value), and the density of water (constant value), respectivelykg/m3
η ,   ζ The dynamical viscosity and the second viscosity (instant local values) of the mixture of water, contaminants, and nanoparticleskg/m3,N·s/m2
f ^ x , f ^ y The components of the magnetic field (local values)A/m
ξ , ξ 0 The dimensionless vertical coordinate, ξ 0 corresponds to z = 0 .dimensionless
z The combined coordinates–time complex-valued variables (both in real and in imaginary parts)
R The function on z equal to ρ 0 v 0 ρ where index 0 related to values for z = 0 m/s
U The constant (function after constant variation) in the u function masks−1
ω i , ω c i The harmonics of u function, critical harmonics, respectivelys−1
The referenced and new numerical values are presented in the next Table.

Values

NotationMeaningUnit of MeasurementValue
r s The characteristic radius of the adsorption for the Fe3O4 nanoparticle m × 10 8 2.1013
r a The characteristic radius of the aggregation for the Fe3O4 nanoparticle m × 10 8 1.41842
τ s The characteristic time of the adsorption for the Fe3O4 nanoparticle s × 10 2 4.2026
τ a The characteristic time of the aggregation for the Fe3O4 nanoparticle s × 10 2 2.83684
k a a Coefficient of aggregation ability decreasing after act of aggregation for Fe3O4 nanoparticle and Ca(HCO3)2 as contaminantn/a 1.1
k a s Coefficient of aggregation ability decreasing after act of adsorption for Fe3O4 nanoparticle and Ca(HCO3)2 as contaminantn/a 10.2
k s a Coefficient of adsorption ability decreasing after act of aggregation for Fe3O4 nanoparticle and Ca(HCO3)2 as contaminantn/a 10.5
k s s Coefficient of adsorption ability decreasing after act of adsorption for Fe3O4 nanoparticle and Ca(HCO3)2 as contaminantn/a 1.6
a 0 Initial value of the channel widthm 0.1
Ω The standard (constant) rate of channel squeezing m / s 0.005
n/aThe final value of the channel width, registered in numeric experiment with strong magnetic field m 0.05
n/aThe width of the unsteady layer %~30

References

  1. Alfven, H. Existence of electromagnetic-hydrodynamic waves. Nature 1942, 150, 405–406. [Google Scholar] [CrossRef]
  2. Hussain, Z.; Zeesahan, R.; Shahzad, M.; Ali, M.; Sultan, F.; Anter, A.M.; Zhang, H.; Khan, N. An optimised stability model for the magnetohydrodynamic fluid. Pramana 2021, 95, 27. [Google Scholar] [CrossRef]
  3. Hussain, Z.; Ali, M.; Shahzad, M.; Sultan, F. Optimized wave perturbation for the linear instability of magnetohydrodynamics in plane Poiseuille flow. Pramana 2020, 94, 49. [Google Scholar] [CrossRef]
  4. Hussain, Z.; Hussain, S.; Kong, T. Instability of MHD Couette flow of an electrically conducting fluid. AIP Adv. 2018, 8, 105209. [Google Scholar] [CrossRef]
  5. Hussain, Z.; Zuev, S.; Kabobel, A.; Ali, M.; Sultan, F.; Shahzad, M. MHD instability of two fluids between parallel plates. Appl. Nanosci. 2020, 10, 5211–5218. [Google Scholar] [CrossRef]
  6. Dalkılıç, A.S.; Yalçın, G.; Küçükyıldırım, B.O.; Öztuna, S.; Eker, A.A.; Jumpholkul, C.; Nakkaew, S.; Wongwises, S. Experimental study on the thermal conductivity of water-based CNT-SiO2 hybrid nanofluids. Int. Commun. Heat Mass Transf. 2018, 99, 18–25. [Google Scholar] [CrossRef]
  7. Anisur, R.; Xu, W.; Li, K.; Dou, H.-S.; Khoo, B.C.; Mao, J. Influence of Magnetic Force on the Flow Stability in a Rectangular duct. Feb. Adv. Appl. Math. Mech. 2019, 11, 24–37. [Google Scholar] [CrossRef]
  8. Hussain, Z.; Abbasi, A.Z.; Ahmad, R.; Bukhari, H.; Shahzad, M.; Sultan, F.; Ali, M. Vibrio cholerae dynamics in drinking water: Mathematical and statistical analysis. Appl. Nanosci. 2020, 10, 4519–4522. [Google Scholar] [CrossRef]
  9. Zainal, N.; Nazar, R.; Naganthran, K.; Pop, I. Unsteady MHD mixed convection flow in hybrid nanofluid at three-dimensional stagnation point. Mathematics 2021, 9, 549. [Google Scholar] [CrossRef]
  10. Zainal, N.A.; Nazar, R.; Naganthran, K.; Pop, I. Unsteady MHD stagnation point flow induced by exponentially permeable stretching/shrinking sheet of hybrid nanofluid. Eng. Sci. Technol. 2021, 24, 1201–1210. [Google Scholar] [CrossRef]
  11. Sharma, K.; Kumar, S.; Narwal, A.; Mebarek-Oudina, F.; Animasaun, I.L. Convective MHD Fluid flow over Stretchable Rotating Disks with Dufour and Soret Effects. Int. J. Appl. Comput. Math. 2022, 8, 159. [Google Scholar] [CrossRef]
  12. Kumar, S.; Sharma, K. Mathematical modeling of MHD flow and radiative heat transfer past a moving porous rotating disk with Hall effect. Multidiscip. Model. Mater. Struct. 2022, 18, 445–458. [Google Scholar] [CrossRef]
  13. Hussain, Z.; Khan, N.; Gul, T.; Ali, M.; Shahzad, M.; Sultan, F. Instability of magneto hydro dynamics Couette flow for electrically conducting fluid through porous media. Appl. Nanosci. 2020, 10, 5125–5134. [Google Scholar] [CrossRef]
  14. Downey, J.P.; Pojman, J.A. Polymer Research in Microgravity: Polymerization and Processing; American Chemical Society: Washington, DC, USA, 2001. [Google Scholar]
  15. Jing, D.; Hu, Y.; Liu, M.; Wei, J.; Guo, L. Preparation of highly dispersed nanofluid and CFD study of its utilization in a concentrating PV/T system. Sol. Energy 2015, 112, 30–40. [Google Scholar] [CrossRef]
  16. Kandelousi, M.S. Effect of spatially variable magnetic field on ferrofluid flow and heat transfer considering constant heat flux boundary condition. Eur. Phys. J. Plus 2014, 129, 248. [Google Scholar] [CrossRef]
  17. Hussanan, A.; Khan, I.; Hashim, H.; Mohamed, M.K.A.; Ishak, N.; Sarif, N.M.; Salleh, M.Z. Unsteady MHD flow of some nanofluids past an accelerated vertical plate embedded in a porous medium. J. Teknol. 2016, 78, 121–126. [Google Scholar] [CrossRef]
  18. Atkins, P.W.; De Paula, J.; Keeler, J. Atkins’ Physical Chemistry, 11th ed.; Oxford University Press: Oxford, UK, 2018; ISBN 978-0-19-876986-6. [Google Scholar]
  19. Landau, L.D.; Lifshitz, E.M. Fluid Mechanics, 2nd ed.; Pergamon Press: Oxford, UK, 1987; 539p. [Google Scholar]
  20. Griffiths, D.J. Introduction to Electrodynamics, 4th ed.; Cambridge University Press: Cambridge, UK, 2017; ISBN 978-1-108-42041-9. [Google Scholar] [CrossRef]
  21. Babichev, A.P.; Babushkina, N.A.; Bratkovskii, A.M. Physical Values; Energoatomizdat: Moscow, Russia, 1991; pp. 123, 124, 370, 376. (In Russian) [Google Scholar]
  22. Sobamowo, M.G.; Akinshilo, A.T. On the analysis of squeezing flow of nanofluid between two parallel plates under the influence of magnetic field. Alex. Eng. J. 2018, 57, 1413–1423. [Google Scholar] [CrossRef]
  23. Kandasamy, R.; Zailani, N.A.B.M.; Jaafar, F.N.B. Impact of nanoparticle volume fraction on squeezed MHD water based Cu, Al2O3 and SWCNTs flow over a porous sensor surface. St. Petersb. Polytech. Univ. J. Phys. Math. 2017, 3, 308–321. [Google Scholar]
  24. Dawar, A.; Shah, Z.; Khan, W.; Idrees, M.; Islam, S. Unsteady squeezing flow of magnetohydrodynamic carbon nanotube nanofluid in rotating channels with entropy generation and viscous dissipation. Adv. Mech. Eng. 2019, 11, 1687814018823100. [Google Scholar] [CrossRef]
Figure 1. Schematic vision of the channel and flow.
Figure 1. Schematic vision of the channel and flow.
Symmetry 14 01728 g001
Figure 2. Magnetic field configurations for start squeezing (a), ¼ squeezing stroke (b), ½ squeezing stroke (c), and ¾ squeezing stroke (d).
Figure 2. Magnetic field configurations for start squeezing (a), ¼ squeezing stroke (b), ½ squeezing stroke (c), and ¾ squeezing stroke (d).
Symmetry 14 01728 g002
Figure 3. Relative concentration of nanoparticles versus vertical position without a magnetic field, at t = 1 s (a), 10 s (b), and 18 s (c): reduced conc. (blue), standard conc. (green), increased conc. (red).
Figure 3. Relative concentration of nanoparticles versus vertical position without a magnetic field, at t = 1 s (a), 10 s (b), and 18 s (c): reduced conc. (blue), standard conc. (green), increased conc. (red).
Symmetry 14 01728 g003
Figure 4. Relative concentration of nanoparticles versus vertical position under the 0.6 T magnetic field, at t = 1 s (a), 10 s (b), and 18 s (c): reduced conc. (blue), standard conc. (green), increased conc. (red).
Figure 4. Relative concentration of nanoparticles versus vertical position under the 0.6 T magnetic field, at t = 1 s (a), 10 s (b), and 18 s (c): reduced conc. (blue), standard conc. (green), increased conc. (red).
Symmetry 14 01728 g004
Figure 5. Relative concentration of nanoparticles versus vertical position under the 0.9 T magnetic field, at t = 1 s (a), 10 s (b), and 18 s (c): reduced conc. (blue), standard conc. (green), increased conc. (red).
Figure 5. Relative concentration of nanoparticles versus vertical position under the 0.9 T magnetic field, at t = 1 s (a), 10 s (b), and 18 s (c): reduced conc. (blue), standard conc. (green), increased conc. (red).
Symmetry 14 01728 g005
Figure 6. Relative concentration of nanoparticles versus magnetic field magnitude at t = 1 s: reduced conc. (blue), standard conc. (green), increased conc. (red).
Figure 6. Relative concentration of nanoparticles versus magnetic field magnitude at t = 1 s: reduced conc. (blue), standard conc. (green), increased conc. (red).
Symmetry 14 01728 g006
Figure 7. Relative concentrations of nanoparticles versus vertical position under 0 T (a), 0.6 T (b), and 0.9 T (c) magnetic fields with a double squeezing rate, at t = 0.5 s (blue), 5 s (green), and 9 s (red).
Figure 7. Relative concentrations of nanoparticles versus vertical position under 0 T (a), 0.6 T (b), and 0.9 T (c) magnetic fields with a double squeezing rate, at t = 0.5 s (blue), 5 s (green), and 9 s (red).
Symmetry 14 01728 g007
Figure 8. Relative concentrations of nanoparticles versus vertical position under 0 T (a), 0.6 T (b), and 0.9 T (c) magnetic fields with a half squeezing rate, at t = 2 s (blue), 20 s (green), and 36 s (red).
Figure 8. Relative concentrations of nanoparticles versus vertical position under 0 T (a), 0.6 T (b), and 0.9 T (c) magnetic fields with a half squeezing rate, at t = 2 s (blue), 20 s (green), and 36 s (red).
Symmetry 14 01728 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zuev, S.; Kabalyants, P.; Hussain, Z. A Model of Water Treatment by Nanoparticles in a Channel with Adjustable Width under a Magnetic Field. Symmetry 2022, 14, 1728. https://doi.org/10.3390/sym14081728

AMA Style

Zuev S, Kabalyants P, Hussain Z. A Model of Water Treatment by Nanoparticles in a Channel with Adjustable Width under a Magnetic Field. Symmetry. 2022; 14(8):1728. https://doi.org/10.3390/sym14081728

Chicago/Turabian Style

Zuev, Sergei, Petr Kabalyants, and Zakir Hussain. 2022. "A Model of Water Treatment by Nanoparticles in a Channel with Adjustable Width under a Magnetic Field" Symmetry 14, no. 8: 1728. https://doi.org/10.3390/sym14081728

APA Style

Zuev, S., Kabalyants, P., & Hussain, Z. (2022). A Model of Water Treatment by Nanoparticles in a Channel with Adjustable Width under a Magnetic Field. Symmetry, 14(8), 1728. https://doi.org/10.3390/sym14081728

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop