# Modeling and Analysis of Magnetic Nanoparticles Injection in Water-Oil Two-Phase Flow in Porous Media under Magnetic Field Effect

Authors: Mohamed F. El-Amin; Ahmed M. Saad; Amgad Salama; Shuyu Sun
Journal: Geofluids
 (2017)

---

## Abstract

In this paper, the magnetic nanoparticles are injected into a water-oil, two-phase system under the influence of an external permanent magnetic field. We lay down the mathematical model and provide a set of numerical exercises of hypothetical cases to show how an external magnetic field can influence the transport of nanoparticles in the proposed two-phase system in porous media. We treat the water-nanoparticles suspension as a miscible mixture, whereas it is immiscible with the oil phase. The magnetization properties, the density, and the viscosity of the ferrofluids are obtained based on mixture theory relationships. In the mathematical model, the phase pressure contains additional term to account for the extra pressures due to fluid magnetization effect and the magnetostrictive effect. As a proof of concept, the proposed model is applied on a countercurrent imbibition flow system in which both the displacing and the displaced fluids move in opposite directions. Physical variables, including water-nanoparticles suspension saturation, nanoparticles concentration, and pore wall/throat concentrations of deposited nanoparticles, are investigated under the influence of the magnetic field. Two different locations of the magnet are studied numerically, and variations in permeability and porosity are considered.


---

## Body

## 1. Introduction

Industry is now looking seriously into using nanotechnology as a viable tool to solving new challenges in several fields. In particular, there has been interest among oil and gas production companies to explore using nanotechnology in solving challenges related to unconventional oil and gas reservoirs, such as those found in tight and shale formations [1–7]. Nanotechnology has been used in different areas of the oil and gas industry from exploration, drilling, production to reservoir monitoring, and refining. The conventional Enhanced Oil Recovery (EOR) methods have several problems from high cost to low oil recovery in addition to operations problems especially in thermal and chemical methods. The nature of nanoparticles results in some useful characteristics such as increased surface area, which at the nanoscale size, does matter when it comes to how molecules react to and bond with each other. So, for example, nanoparticles can be used in EOR, because they are small enough to pass through pore throats in typical reservoirs, and they can be retained by the rock. Ju and Fan [8] calibrated a model for nanoparticles transport in two-phase flow in porous media based on the formulation of the colloid model of fine particles transport in two-phase flow in porous media [9]. El-Amin et al. [10–13] have presented modeling and simulations of nanoparticles transport associated with two-phase flow in porous media. On the other hand, experimental studies of using nanoparticles in EOR have been conducted by Suleimanov et al. [14] and Hendraningrat et al. [15].One of the prospective applications of nanotechnology is nanoferrofluids, as the flow of such fluids can be controlled by the introduction of external magnetic field. This opens the way for various applications from directing flow in reservoir monitoring, diverting flow in acid jobs to control and boost of injection fluid advancement during pressure maintenance to increase oil recovery. The idea of using a strong external magnetic field with large magnetic susceptibility fluid is to mobilize ferrofluid through porous media. Oil recovery can be increased by using nanoparticles with electromagnetic properties (such as iron oxide,Fe2O3, and zinc oxide, ZnO) under waves generated from an electromagnetic source. Both direct and alternating magnetic fields are under investigation; in our case here we will focus on direct magnetic field. The magnetization of the particles and their attraction toward the magnet causes flow of the magnetic particles suspension. The movement of the magnetic nanoparticles under the magnetic field effect is independent of the orientation of the magnet. In the last few years, a number of publications have been considered nanoferrofluids in oil and gas recovery or environmental applications (e.g., [16–23]).Borglin et al. [24] conducted experiments to measure the magnetic induction, which converted to magnetic field strength, at various distances in a direction aligned with the poles. McCaig and Clegg [25] presented equations which indicate that ferrofluid magnetization is variable at all locations far away from the magnet due to the decreasing magnetic field strength. The gradient of magnetic field strength varies also with distance from the magnet. Moreover, the force on the ferrofluid decreases with distance from a magnet. Therefore, we may conclude that the distance from the magnet is important and should be reconsidered with taking into consideration its location in particular in the core-scale as well as the direction of flow. On the other hand, it is known that, in the presence of an external magnetic field, the nanoparticles in ferrofluid become magnetized and are pulled toward a magnet. In this work, we attempt to discuss the effect of location of the magnet on the transport of the nanoparticles in porous media.In the current work, we develop a mathematical model to describe the magnetic nanoparticles-water suspension imbibition into an initially oil saturated porous domain under magnetic field effect. The porous medium is considered initially saturated totally with oil except for a residual amount of the other phase. We consider countercurrent imbibition into a small-scale porous core. This countercurrent imbibition refers to the case in which all the porous medium domain boundaries have no flow except one side. Physical variables are investigated under the influences of magnetic field with two different locations of magnet, namely, right and left to the porous core. Numerical experiments for the two cases are performed and results are introduced in graphical representations.


## 2. Modeling and Mathematical Formulation

Consider suspension of magnetic nanoparticles injected in an isothermal incompressible water-oil two-phase flow under an external magnetic field. In the following, we describe the mathematical modeling of the problem under consideration.In the following subsections, we firstly introduce the magnetic force and other magnetic modeling. The magnetic body force, which acts as a body force on the nanoparticles suspension per unit volume, appears in the extended Darcy’s law as presented in the second subsection. Other magnetic parameters such as magnetization and magnetic field strength are also presented in the first subsection. Then, the governing flow equations such as momentum (extended Darcy’s law) and mass conservation (saturation equation) are provided in the second subsection. The third subsection is devoted to the modeling of nanoparticles transport which is coupled with the flow equations via velocity and saturation. Finally, the initial and boundary conditions are presented in the last subsection.


### 2.1. Magnetic Force

The magnetization (polarization) of the nanoparticles suspension interacts with the external permanent magnetic field to produce attractive forces on each particle. The external magnetic force acts as a body force on the nanoparticles suspension per unit volume which can be expressed as [16](1)Fmag=μ0M∂H∂z,where μ0 is the magnetic permeability, M is the magnetization, and H is the magnetic field strength. The magnetization M is a function of H, approximated by (2)M=a1tan-1b1H,where the parameters a1 and b1 depend on the particular type of the ferromagnetic material. The values of the initial susceptibility and the saturation magnetization are controlled by the parameters a1 and b1, respectively. The larger b1 leads to a larger initial susceptibility which corresponds to larger particles or agglomeration of particles. The range of a1 may be of order 104–105, while the order of b1 may be of order 10−6–10−5. The magnetic field strength in 1D may be written as follows [16]: (3)Hz=Brπμ0tan-1abza2+b2+z21/2-tan-1abz+La2+b2+z+L21/2,where Br is the residual magnetization (in this work, Br=1.19 [T]) and L is the distance between the poles of the magnet. In Figure 1, Fmag, Hz, and ∂Hz/∂z are plotted against z.Figure 1
F
m
a
g, Hz, and ∂Hz/∂z are plotted against z.Oldenburg et al. [16] reported that magnetic particles (ferrofluids) could have both miscible and immiscible behaviors with water. For example, an oil-based ferrofluid is immiscible with water while aqueous ferrofluids are miscible in water. The immiscible case should be described by two-phase aqueous/nonaqueous mixtures regardless of external magnetic field strength, while the case of miscible in water should be described as single-phase aqueous mixtures. Under no external magnetic field, the two fluids are miscible. When a magnetic field is applied, it affects the ferrofluid part of the aquatic mixture in such a way that it behaves as immiscible. That is, the ferrofluids motion is dictated by the magnetic field. When the magnetic field is relatively strong, its effect on the ferrofluids may be influential. Therefore, a two-phase flow model is suitable for the flow of ferrofluids under strong external magnetic fields, but when the magnet is at larger distances from the flow (thus segregation effects are unimportant), the fully miscible single-phase liquid model is appropriate. Both of these models assume that the mixture is not so dilute. Pure ferrofluid and dilute colloidal suspension have approximately a thousand particles per cubic meter (see [16, 26]). So, one may conclude that the miscible and immiscible models may be applicable for ferrofluid fractions of 10% in water, and they are probably satisfactory in much more dilute mixtures [16].In the current model, we treat the water-nanoparticles suspension as a miscible mixture while it is immiscible with the oil phase. Oldenburg et al. [16] assumed that magnetization increases linearly with the ferrofluid mass fraction c, (4)Mc=Mc=1c.In the case of the immiscible two-phase (water-oil) model, Oldenburg et al. [16] assumed that magnetization increases linearly with the ferrofluid saturation, (5)MSw=MSw=1Sw.For the immiscible conceptualization, density and viscosity are considered properties of the respective phases, and no mixing relations are required. Therefore, in the current model we consider immiscible two-phase (water-oil) model; however, the water-nanoparticles suspension is a miscible mixture. So, one may expresses the the magnetization of this system as, (6)MSw,c=MSw=1,c=1Swc.Following the work of Reeves and Cranwell [27] and Herbert et al. [28] on brine-water mixtures, Oldenburg et al. [16] have assumed that the volumes of pure water and ferrofluid are additive and define for the mixture density as (7)1ρw=1-cρw,p+cρf,where ρw,p [kg·m−3] is the density of water component and ρf [kg·m−3] is the density of ferrofluid component.The viscosity of the nanoparticles-water mixture is calculated by the following linear relationship [16, 28]: (8)μw=μw,p1+1.35c,where μw,p is the viscosity of pure water.


### 2.2. Flow Model

Considering the external magnetic force, which acts as a body force as explained in the above subsection, the vertical 1D (i.e.,z-direction) velocities of oil phase and nanoparticles-water mixture phase may be written as (9)uo=-Kkroμo∂po∂z-ρog,uw=-Kkrwμw∂pw∂z-ρwcg-μ0MSw,c∂H∂z.The mass conservation equation is (10)ϕ∂Sα∂t+∂uα∂z=0,α=w,o,where ϕ [—] is the porosity, ρα [kg·m−3] is the density of phase α, Sα [—] is the saturation of phase α, and uα [m·s−1] is the velocity of the phase α. w stands for the nanoparticles-water suspension phase, and o stands for the oil phase. K [m2] is the permeability, krα [—] is the relative permeability of phase α, pα [Pa] is the pressure of phase α, g [m·s−2] is the gravitational acceleration, and μα [kg m−1 s−1] is the viscosity of phase α. The fluid saturations for the two-phase flow of water and oil are related by(11)Sw+So=1.In countercurrent imbibition, the sum of the velocities of the wetting and nonwetting phases is zero, that is,(12)ut=uw+uo=0The capillary pressure is defined as(13)pc=po-pw.Since, in this case of study, we have two phases, one magnetized and the other nonmagnetized, the magnetostrictive effect should not be neglected. The magnetized phase pressure will have additional pressure term with the conventional thermodynamic pressure. The resulting pressure is called composite pressure [21], which can be given as(14)pw∗=pw+pm+ps+pn,where pw is the ferrofluid phase dynamic pressure, pm is the fluid magnetic pressure, ps is the magnetostrictive pressure, and pn is the magnetic normal pressure, which is neglected [21] in this study. Also, one can write(15)pc∗=pc-pm+ps,ps=μ0∫0HV∂M∂VH,TdH,pm=μ0∫0HMdH.The specific volume is given as (16)V=∂1ρ=-1ρ2∂ρ, So, (17)ps=-μ0∫0Hρ∂M∂ρH,TdH, where (18)∂M∂ρ=∂M∂a1b1,H∂a1∂ρ+∂M∂b1a1,H∂b1∂ρ+∂M∂Ha1,b1∂H∂ρ. We may approximate ∂a1/∂ρ and ∂b1/∂ρ, when a1,o and b1,o are equal to zero in the nonmagnetized phase, as (19)∂a1∂ρ≈a1-a1,oρw-ρo=a1ρw-ρo,∂b1∂ρ≈b1-b1,oρw-ρo=a1ρw-ρo. Therefore,(20)ps=ρwa1ρw-ρoμ0Htan-1b1H.Similarly, therefore,(21)pm=a1Htan-1b1H-12b1ln⁡b12H2+1.Substituting from (9), (13), (20), and (21) into (12) and eliminating ∂po/∂z, we may have(22)∂pw∗∂z=-fw∂pc∗∂z+ρwfw+ρofog+fwμ0MSw,c∂H∂z,So, the water velocity becomes(23)uw=Kλwfo∂pc∗∂z-Δρg+μ0MSw,c∂H∂z.Therefore, the saturation equation for the water phase becomes(24)ϕ∂Sw∂t+∂∂zKλwfo∂pc∗∂z-Δρg+μ0MSw,c∂H∂z=0,where λw=krw/μw and λo=kro/μo are mobility ratios of water and oil phases, respectively. λt=λw+λo is the total mobility. fw=λw/λt and fo=λo/λt are the flow fraction of water and oil phases, respectively. Δρ=ρw-ρo. The capillary pressure is a function of the normalized saturation which can be given as [29](25)pc=pdln⁡S,where pd is the entry pressure for the imbibition. Moreover, the relative permeabilities are defined as(26)krw=krw0Sa2,kro=kro01-Sb2,where S is the normalized water phase saturation, which is given as(27)S=Sw-Siw1-Sro-Siw,0<S<1,and krw0=krwS=1 and kro0=kroS=0 are the endpoint relative permeability of the water and oil phase, respectively. a2 and b2 are positive numbers. Siw is the irreducible water saturation and Sro is the residual oil saturation.


### 2.3. Nanoparticles Transport Model

Assuming that the nanoparticles exist only in the water phase and have only one size interval. The Brownian diffusion is considered for the nanoparticles and the volumes of pure water and the particles suspension are additive. The transport equation of the nanoparticles-water suspension in the water phase can be written as(28)ϕ∂Swc∂t-∂cs1∂t-∂cs2∂t+∂∂zuwc-ϕSwDdiff+Ddisp∂c∂z=0,where c [m3·m−3] is the volume concentration of nanoparticles in the water phase. cs1 [m3·m−3] is the volume of the nanoparticles in contact with the water phase available on the pore surfaces per unit bulk volume of the porous medium. cs2 [m3·m−3] is the volume of the nanoparticles entrapped in pore throats from the water phase per unit bulk volume of porous medium due to plugging and bridging. Ddiff [m2·s−1] is the molecular diffusion coefficient which can be calculated using the Stokes-Einstein equation,(29)Ddiff=τkBT3πμwd,and τ is the tortuosity of the flow. kB is the Boltzmann constant; T [K] absolute temperature. The mechanical dispersion coefficient Ddisp [m2·s−1] is a function of Darcy’s velocity and may be given as(30)ϕSwDdisp=dl,wuw,where dl,w is the longitudinal dispersion coefficient.The modified Gruesbeck-Collins model [30] for the surface deposition is used in this study [8, 12]. In the presence of a critical velocity of the surface deposition only particle retention occurs while above it retention and entrainment of the nanoparticles take place simultaneously, which can be modeled as(31)∂cs1∂t=γduwc,uw≤ucγduwc-γeuw-uccs1,uw>uc.Also, the rate of entrapment of the nanoparticles in the water phase is given by(32)∂cs2∂t=γptuwc,where γd [m−1] is the rate coefficient for surface retention of the nanoparticles in the water phase. γe [m−1] is the rate coefficient for entrainment of the nanoparticles. uc is the critical velocity for the water phase. γpt [m−1] is the pore throat blocking constant.The porosity variation due to nanoparticles deposition is given as [8, 9](33)ϕ=ϕ0-δϕ,δϕ=cs1+cs2,where ϕ0 is the initial porosity. Therefore, permeability can be changed as [8](34)K=K01-fkf+fϕϕ0l,where K0 is the initial permeability and K=kI,  K0=k0I, where k [m2] and k0 [m2] are positive constants and I is a unit matrix. kf is constant for fluid seepage allowed by the plugged pores. The flow efficiency factor expressing the fraction of unplugged pores available for flow is given by(35)f=1-γfcs2,where γf is the coefficient of flow efficiency for the nanoparticles. The exponent l has a value of the range from 2.5 to 3.5.The variation of relative permeabilities due to the nanoparticles reposition is given by El-Amin et al. [12] as follows:(36)krα,p=1+raθα-1krα,α=w,o,where θα is the ratio of the phase α relative permeability due to nanoparticles adhering; that is,(37)krα,c=θαkrα,α=w,o,where krα,c is the relative permeabilities of water/oil phase when the surfaces per unit bulk volume of the porous media is completely occupied by the nanoparticles. ra=atot/asp such that(38)asp=AϕϕK1/2is a specific area of the sand core and A [m2] is the cross-sectional area. Also, the total surface area in contact with fluids per unit bulk volume is defined as [8](39)atot=6βdδϕand d is the diameter of the nanoparticles in a given interval size.


### 2.4. Initial and Boundary Conditions

At the beginning of the flow the following initial conditions are considered:(40)Sw=Sw0att=0,0≤z≤h,c=cs1=cs2=0att=0,0≤z≤h,where h is the rock depth and Sw0 is the initial water saturation.The boundary conditions are(41)Sw=1-So0,c=c0,cs1=cs2=0att>0,z=0,dSwdz=dcdz=dcs1dz=dcs2dz=0att>0,z=h,where c0 is the concentration of nanoparticles-water suspension on the inlet boundary.


## 2.1. Magnetic Force

The magnetization (polarization) of the nanoparticles suspension interacts with the external permanent magnetic field to produce attractive forces on each particle. The external magnetic force acts as a body force on the nanoparticles suspension per unit volume which can be expressed as [16](1)Fmag=μ0M∂H∂z,where μ0 is the magnetic permeability, M is the magnetization, and H is the magnetic field strength. The magnetization M is a function of H, approximated by (2)M=a1tan-1b1H,where the parameters a1 and b1 depend on the particular type of the ferromagnetic material. The values of the initial susceptibility and the saturation magnetization are controlled by the parameters a1 and b1, respectively. The larger b1 leads to a larger initial susceptibility which corresponds to larger particles or agglomeration of particles. The range of a1 may be of order 104–105, while the order of b1 may be of order 10−6–10−5. The magnetic field strength in 1D may be written as follows [16]: (3)Hz=Brπμ0tan-1abza2+b2+z21/2-tan-1abz+La2+b2+z+L21/2,where Br is the residual magnetization (in this work, Br=1.19 [T]) and L is the distance between the poles of the magnet. In Figure 1, Fmag, Hz, and ∂Hz/∂z are plotted against z.Figure 1
F
m
a
g, Hz, and ∂Hz/∂z are plotted against z.Oldenburg et al. [16] reported that magnetic particles (ferrofluids) could have both miscible and immiscible behaviors with water. For example, an oil-based ferrofluid is immiscible with water while aqueous ferrofluids are miscible in water. The immiscible case should be described by two-phase aqueous/nonaqueous mixtures regardless of external magnetic field strength, while the case of miscible in water should be described as single-phase aqueous mixtures. Under no external magnetic field, the two fluids are miscible. When a magnetic field is applied, it affects the ferrofluid part of the aquatic mixture in such a way that it behaves as immiscible. That is, the ferrofluids motion is dictated by the magnetic field. When the magnetic field is relatively strong, its effect on the ferrofluids may be influential. Therefore, a two-phase flow model is suitable for the flow of ferrofluids under strong external magnetic fields, but when the magnet is at larger distances from the flow (thus segregation effects are unimportant), the fully miscible single-phase liquid model is appropriate. Both of these models assume that the mixture is not so dilute. Pure ferrofluid and dilute colloidal suspension have approximately a thousand particles per cubic meter (see [16, 26]). So, one may conclude that the miscible and immiscible models may be applicable for ferrofluid fractions of 10% in water, and they are probably satisfactory in much more dilute mixtures [16].In the current model, we treat the water-nanoparticles suspension as a miscible mixture while it is immiscible with the oil phase. Oldenburg et al. [16] assumed that magnetization increases linearly with the ferrofluid mass fraction c, (4)Mc=Mc=1c.In the case of the immiscible two-phase (water-oil) model, Oldenburg et al. [16] assumed that magnetization increases linearly with the ferrofluid saturation, (5)MSw=MSw=1Sw.For the immiscible conceptualization, density and viscosity are considered properties of the respective phases, and no mixing relations are required. Therefore, in the current model we consider immiscible two-phase (water-oil) model; however, the water-nanoparticles suspension is a miscible mixture. So, one may expresses the the magnetization of this system as, (6)MSw,c=MSw=1,c=1Swc.Following the work of Reeves and Cranwell [27] and Herbert et al. [28] on brine-water mixtures, Oldenburg et al. [16] have assumed that the volumes of pure water and ferrofluid are additive and define for the mixture density as (7)1ρw=1-cρw,p+cρf,where ρw,p [kg·m−3] is the density of water component and ρf [kg·m−3] is the density of ferrofluid component.The viscosity of the nanoparticles-water mixture is calculated by the following linear relationship [16, 28]: (8)μw=μw,p1+1.35c,where μw,p is the viscosity of pure water.


## 2.2. Flow Model

Considering the external magnetic force, which acts as a body force as explained in the above subsection, the vertical 1D (i.e.,z-direction) velocities of oil phase and nanoparticles-water mixture phase may be written as (9)uo=-Kkroμo∂po∂z-ρog,uw=-Kkrwμw∂pw∂z-ρwcg-μ0MSw,c∂H∂z.The mass conservation equation is (10)ϕ∂Sα∂t+∂uα∂z=0,α=w,o,where ϕ [—] is the porosity, ρα [kg·m−3] is the density of phase α, Sα [—] is the saturation of phase α, and uα [m·s−1] is the velocity of the phase α. w stands for the nanoparticles-water suspension phase, and o stands for the oil phase. K [m2] is the permeability, krα [—] is the relative permeability of phase α, pα [Pa] is the pressure of phase α, g [m·s−2] is the gravitational acceleration, and μα [kg m−1 s−1] is the viscosity of phase α. The fluid saturations for the two-phase flow of water and oil are related by(11)Sw+So=1.In countercurrent imbibition, the sum of the velocities of the wetting and nonwetting phases is zero, that is,(12)ut=uw+uo=0The capillary pressure is defined as(13)pc=po-pw.Since, in this case of study, we have two phases, one magnetized and the other nonmagnetized, the magnetostrictive effect should not be neglected. The magnetized phase pressure will have additional pressure term with the conventional thermodynamic pressure. The resulting pressure is called composite pressure [21], which can be given as(14)pw∗=pw+pm+ps+pn,where pw is the ferrofluid phase dynamic pressure, pm is the fluid magnetic pressure, ps is the magnetostrictive pressure, and pn is the magnetic normal pressure, which is neglected [21] in this study. Also, one can write(15)pc∗=pc-pm+ps,ps=μ0∫0HV∂M∂VH,TdH,pm=μ0∫0HMdH.The specific volume is given as (16)V=∂1ρ=-1ρ2∂ρ, So, (17)ps=-μ0∫0Hρ∂M∂ρH,TdH, where (18)∂M∂ρ=∂M∂a1b1,H∂a1∂ρ+∂M∂b1a1,H∂b1∂ρ+∂M∂Ha1,b1∂H∂ρ. We may approximate ∂a1/∂ρ and ∂b1/∂ρ, when a1,o and b1,o are equal to zero in the nonmagnetized phase, as (19)∂a1∂ρ≈a1-a1,oρw-ρo=a1ρw-ρo,∂b1∂ρ≈b1-b1,oρw-ρo=a1ρw-ρo. Therefore,(20)ps=ρwa1ρw-ρoμ0Htan-1b1H.Similarly, therefore,(21)pm=a1Htan-1b1H-12b1ln⁡b12H2+1.Substituting from (9), (13), (20), and (21) into (12) and eliminating ∂po/∂z, we may have(22)∂pw∗∂z=-fw∂pc∗∂z+ρwfw+ρofog+fwμ0MSw,c∂H∂z,So, the water velocity becomes(23)uw=Kλwfo∂pc∗∂z-Δρg+μ0MSw,c∂H∂z.Therefore, the saturation equation for the water phase becomes(24)ϕ∂Sw∂t+∂∂zKλwfo∂pc∗∂z-Δρg+μ0MSw,c∂H∂z=0,where λw=krw/μw and λo=kro/μo are mobility ratios of water and oil phases, respectively. λt=λw+λo is the total mobility. fw=λw/λt and fo=λo/λt are the flow fraction of water and oil phases, respectively. Δρ=ρw-ρo. The capillary pressure is a function of the normalized saturation which can be given as [29](25)pc=pdln⁡S,where pd is the entry pressure for the imbibition. Moreover, the relative permeabilities are defined as(26)krw=krw0Sa2,kro=kro01-Sb2,where S is the normalized water phase saturation, which is given as(27)S=Sw-Siw1-Sro-Siw,0<S<1,and krw0=krwS=1 and kro0=kroS=0 are the endpoint relative permeability of the water and oil phase, respectively. a2 and b2 are positive numbers. Siw is the irreducible water saturation and Sro is the residual oil saturation.


## 2.3. Nanoparticles Transport Model

Assuming that the nanoparticles exist only in the water phase and have only one size interval. The Brownian diffusion is considered for the nanoparticles and the volumes of pure water and the particles suspension are additive. The transport equation of the nanoparticles-water suspension in the water phase can be written as(28)ϕ∂Swc∂t-∂cs1∂t-∂cs2∂t+∂∂zuwc-ϕSwDdiff+Ddisp∂c∂z=0,where c [m3·m−3] is the volume concentration of nanoparticles in the water phase. cs1 [m3·m−3] is the volume of the nanoparticles in contact with the water phase available on the pore surfaces per unit bulk volume of the porous medium. cs2 [m3·m−3] is the volume of the nanoparticles entrapped in pore throats from the water phase per unit bulk volume of porous medium due to plugging and bridging. Ddiff [m2·s−1] is the molecular diffusion coefficient which can be calculated using the Stokes-Einstein equation,(29)Ddiff=τkBT3πμwd,and τ is the tortuosity of the flow. kB is the Boltzmann constant; T [K] absolute temperature. The mechanical dispersion coefficient Ddisp [m2·s−1] is a function of Darcy’s velocity and may be given as(30)ϕSwDdisp=dl,wuw,where dl,w is the longitudinal dispersion coefficient.The modified Gruesbeck-Collins model [30] for the surface deposition is used in this study [8, 12]. In the presence of a critical velocity of the surface deposition only particle retention occurs while above it retention and entrainment of the nanoparticles take place simultaneously, which can be modeled as(31)∂cs1∂t=γduwc,uw≤ucγduwc-γeuw-uccs1,uw>uc.Also, the rate of entrapment of the nanoparticles in the water phase is given by(32)∂cs2∂t=γptuwc,where γd [m−1] is the rate coefficient for surface retention of the nanoparticles in the water phase. γe [m−1] is the rate coefficient for entrainment of the nanoparticles. uc is the critical velocity for the water phase. γpt [m−1] is the pore throat blocking constant.The porosity variation due to nanoparticles deposition is given as [8, 9](33)ϕ=ϕ0-δϕ,δϕ=cs1+cs2,where ϕ0 is the initial porosity. Therefore, permeability can be changed as [8](34)K=K01-fkf+fϕϕ0l,where K0 is the initial permeability and K=kI,  K0=k0I, where k [m2] and k0 [m2] are positive constants and I is a unit matrix. kf is constant for fluid seepage allowed by the plugged pores. The flow efficiency factor expressing the fraction of unplugged pores available for flow is given by(35)f=1-γfcs2,where γf is the coefficient of flow efficiency for the nanoparticles. The exponent l has a value of the range from 2.5 to 3.5.The variation of relative permeabilities due to the nanoparticles reposition is given by El-Amin et al. [12] as follows:(36)krα,p=1+raθα-1krα,α=w,o,where θα is the ratio of the phase α relative permeability due to nanoparticles adhering; that is,(37)krα,c=θαkrα,α=w,o,where krα,c is the relative permeabilities of water/oil phase when the surfaces per unit bulk volume of the porous media is completely occupied by the nanoparticles. ra=atot/asp such that(38)asp=AϕϕK1/2is a specific area of the sand core and A [m2] is the cross-sectional area. Also, the total surface area in contact with fluids per unit bulk volume is defined as [8](39)atot=6βdδϕand d is the diameter of the nanoparticles in a given interval size.


## 2.4. Initial and Boundary Conditions

At the beginning of the flow the following initial conditions are considered:(40)Sw=Sw0att=0,0≤z≤h,c=cs1=cs2=0att=0,0≤z≤h,where h is the rock depth and Sw0 is the initial water saturation.The boundary conditions are(41)Sw=1-So0,c=c0,cs1=cs2=0att>0,z=0,dSwdz=dcdz=dcs1dz=dcs2dz=0att>0,z=h,where c0 is the concentration of nanoparticles-water suspension on the inlet boundary.


## 3. Results and Discussion

The above highly nonlinear parabolic partial differential equation is solved numerically using an efficient algorithm [31]. The spatial discretization is handled by Galerkin method, while an adaptive time step is used with the time integration. The above governing equations (24), (28), (31), and (32) are solved along with their initial and boundary conditions (40)-(41). The following parameters values are used in the computations; namely [8], γd=16 [m-1], γpt=1.28 [m-1], γe=30 [m-1], uc=4.6×10-6 [m·s−1], and D=5.6×10-8 [m2·s−1]. The nanoparticles diameter is taken as 40 nm and the inlet concentration c0=0.0009. Other different values of c0 were taken into account in our previous study [12]. The remaining model parameters are [12, 32]; h=0.2 m, Swr=Sor=0.1, aw=ao=0.5, μw=μo=0.001 [kg m-1s-1], ρw=1000 [kg·m−3], ρo=660 [kg·m], ϕ0=0.3, l=3, kf=0.6, γf=0.01, krw0=kro0=1, a2=b2=4, β=0.8, k0=20×10-15 [m2], θo=3, θw=0.6. Moreover, the magnetic field parameters are a1=1.5×104, b1=2.4×10-5 (when the magnet is on the right side of the core). But when the magnet is on the left side of the core, they are taken as a1=3×102, b1=3×10-5.In the following we consider two different locations of the magnet, namely, right and left side, with respect to the core and its boundary conditions (see Figures2 and 3). On other words, the magnet location and flow directions are very important factors in this problem. If the magnet location is on the right side of the no-flow boundary of the core, the magnetic field will be assisting the flow, while the opposite is true when the magnet is placed on the left side of the inflow boundary of the rock (opposing flow).Figure 2
Magnet is placed on the right side of the core, with a core length,h=0.2 m, and magnet parameters, a1=1.5×104, b1=2.4×10-5, Br=1.19, a=b=0.0254, L=0.127 m, and μ0=4π10-7, and magnet positioned at x=0.2 m, y=x-0.2 m.Figure 3
Magnet is placed on the left side of the core, with a core length,h=0.2 m, and magnet parameters, a1=3×102, b1=3×10-5, Br=1.19, a=b=0.0254, L=0.127 m, and μ0=4π10-7, and magnet positioned at x=0.2 m, y=x-0.2 m.


### 3.1. Magnet at the Right Side of the Core

Firstly, we present the first case in Figures4–9. Figure 4 illustrates the nanoparticles-water saturation profiles against the rock length with and without magnetic field effect for various values of the time of imbibition. It is interesting to note from this figure the effect of placing the magnet on the right side of the rock as it increases the saturation of nanoparticles-water suspension on this side while it decreases the saturation on the left side of the rock near to the flow inlet. This can be explained through an important fact that is the magnetic force is proportional to magnetic field strength, so fluid that is closer to the magnet is more strongly magnetized and pulled more strongly toward the magnet. For instance, this is in contrast to gravity driven flow in which gravitational force is independent of position and fluid is pulled uniformly under gravity.Figure 4
Saturations profiles against the core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 5
Nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 6
Pore wall deposited nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 7
Pore throat deposited nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 8
Permeability variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 9
Porosity variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.The nanoparticles concentration is plotted against the rock length in Figure5, with and without magnetic field effect for various values of the imbibition time. From this figure, it can be seen that after long time of imbibition, the nanoparticles concentration decreases slightly under the effect of the magnetic field. Figure 6 shows the concentration of deposited nanoparticles on the pore wall as plotted against the core length with and without magnetic field effect for various values imbibition times. Also, Figure 7 shows the concentration of deposited nanoparticles on the pore throat which is plotted against the core length with and without magnetic field effect for various values imbibition times. As expected the deposited nanoparticles have an opposite behavior of the nanoparticles in the water; that is, when the nanoparticles concentration in the water decreases, the deposited nanoparticles concentration increases.In Figure8, the permeability variation is plotted against the distance, with and without the magnetic field effect for various values of imbibition time when the magnet is located on the right side of the core. It is clear from this figure that the reduction in permeability is related to the quantity of the precipitation of nanoparticles on the wall and the imbibition time. Also, in Figure 9, the porosity variation is plotted against the distance with and without the magnetic field effect for various values of imbibition time when the magnet is located on the right side of the core. It can be seen that the change in porosity is related to the quantity of nanoparticles plugging the throat of the pore and imbibition time.


### 3.2. Magnet at the Left Side of the Core

Now, we consider the second case, when the magnet is to the left of the rock (which opposes the flow), in Figures10–13. Saturations profiles are plotted in Figure 10, against the core distance with and without magnetic field effect for various values of imbibition time when the magnet is placed on left side of the core. An interesting behavior can be seen in this figure. Unlike the case of the right magnet, the left magnet resists the flow and reduces water invasion. Moreover, Figure 11 shows the profiles of nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time, in the case of left location of the magnet. It is clear from this figure that the concentration of nanoparticles in water increases especially after long imbibition time. Figure 12 illustrates deposited nanoparticles concentration on the pore wall against the core length, with and without magnetic field effect for various values of imbibition time when the magnet is located on the left side of the core. This figure shows that the magnetic field reduces significantly the deposited nanoparticles on the pore surface. Finally, Figure 13 illustrates the nanoparticles concentration plugging pore throat against core length with and without magnetic field effect for various values of imbibition time when the magnet is located on the left side of the core. One may notice that the magnetic field reduces the nanoparticles plugging pore throat.Figure 10
Saturations profiles against the core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure 11
Nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure 12
Pore wall deposited nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure 13
Pore throat deposited nanoparticles concentration against core distance with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure14 illustrates the permeability profiles against the distance with various values of imbibition time, with and without magnetic field effect when the magnet is on the left side. Similarly, the permeability decreases as a consequence of the precipitation of nanoparticles. The permeability has a minor reduction compared to the initial permeability. Figure 15 shows the porosity variation against the distance with various values of imbibition time, with and without magnetic field effect when the magnet is located on the left side. It is interesting to note reduction in the porosity which is a consequence of the precipitation of nanoparticles on the pore walls. The porosity variation is negligibly small.Figure 14
Permeability variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure 15
Porosity variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.


## 3.1. Magnet at the Right Side of the Core

Firstly, we present the first case in Figures4–9. Figure 4 illustrates the nanoparticles-water saturation profiles against the rock length with and without magnetic field effect for various values of the time of imbibition. It is interesting to note from this figure the effect of placing the magnet on the right side of the rock as it increases the saturation of nanoparticles-water suspension on this side while it decreases the saturation on the left side of the rock near to the flow inlet. This can be explained through an important fact that is the magnetic force is proportional to magnetic field strength, so fluid that is closer to the magnet is more strongly magnetized and pulled more strongly toward the magnet. For instance, this is in contrast to gravity driven flow in which gravitational force is independent of position and fluid is pulled uniformly under gravity.Figure 4
Saturations profiles against the core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 5
Nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 6
Pore wall deposited nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 7
Pore throat deposited nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 8
Permeability variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.Figure 9
Porosity variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.The nanoparticles concentration is plotted against the rock length in Figure5, with and without magnetic field effect for various values of the imbibition time. From this figure, it can be seen that after long time of imbibition, the nanoparticles concentration decreases slightly under the effect of the magnetic field. Figure 6 shows the concentration of deposited nanoparticles on the pore wall as plotted against the core length with and without magnetic field effect for various values imbibition times. Also, Figure 7 shows the concentration of deposited nanoparticles on the pore throat which is plotted against the core length with and without magnetic field effect for various values imbibition times. As expected the deposited nanoparticles have an opposite behavior of the nanoparticles in the water; that is, when the nanoparticles concentration in the water decreases, the deposited nanoparticles concentration increases.In Figure8, the permeability variation is plotted against the distance, with and without the magnetic field effect for various values of imbibition time when the magnet is located on the right side of the core. It is clear from this figure that the reduction in permeability is related to the quantity of the precipitation of nanoparticles on the wall and the imbibition time. Also, in Figure 9, the porosity variation is plotted against the distance with and without the magnetic field effect for various values of imbibition time when the magnet is located on the right side of the core. It can be seen that the change in porosity is related to the quantity of nanoparticles plugging the throat of the pore and imbibition time.


## 3.2. Magnet at the Left Side of the Core

Now, we consider the second case, when the magnet is to the left of the rock (which opposes the flow), in Figures10–13. Saturations profiles are plotted in Figure 10, against the core distance with and without magnetic field effect for various values of imbibition time when the magnet is placed on left side of the core. An interesting behavior can be seen in this figure. Unlike the case of the right magnet, the left magnet resists the flow and reduces water invasion. Moreover, Figure 11 shows the profiles of nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time, in the case of left location of the magnet. It is clear from this figure that the concentration of nanoparticles in water increases especially after long imbibition time. Figure 12 illustrates deposited nanoparticles concentration on the pore wall against the core length, with and without magnetic field effect for various values of imbibition time when the magnet is located on the left side of the core. This figure shows that the magnetic field reduces significantly the deposited nanoparticles on the pore surface. Finally, Figure 13 illustrates the nanoparticles concentration plugging pore throat against core length with and without magnetic field effect for various values of imbibition time when the magnet is located on the left side of the core. One may notice that the magnetic field reduces the nanoparticles plugging pore throat.Figure 10
Saturations profiles against the core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure 11
Nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure 12
Pore wall deposited nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure 13
Pore throat deposited nanoparticles concentration against core distance with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure14 illustrates the permeability profiles against the distance with various values of imbibition time, with and without magnetic field effect when the magnet is on the left side. Similarly, the permeability decreases as a consequence of the precipitation of nanoparticles. The permeability has a minor reduction compared to the initial permeability. Figure 15 shows the porosity variation against the distance with various values of imbibition time, with and without magnetic field effect when the magnet is located on the left side. It is interesting to note reduction in the porosity which is a consequence of the precipitation of nanoparticles on the pore walls. The porosity variation is negligibly small.Figure 14
Permeability variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.Figure 15
Porosity variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.


## 4. Conclusions

This paper was devoted to study the magnetic field effects on the the transport of magnetic nanoparticles injected into a two-phase water-oil system in porous media. Firstly, we developed the mathematical model which governs the transport of magnetic nanoparticles in a two-phase, water-oil system in porous media under the effect of an external magnetic field. The countercurrent imbibition in a small-scale porous medium core is considered as an example. Both Brownian diffusion and mechanical dispersion are taken into consideration, with the assumption that only one size interval of the nanoparticles exists. Variation of porosity, permeability, and relative permeabilities due to particles deposition are also considered. Immiscible mixture relationships have been used to determine the magnetization properties, density, and viscosity of this ferrofluid. In such flow field, the pressure includes additional terms, namely, the ferrofluid dynamic pressure, the fluid magnetic pressure, and the magnetostrictive pressure. Two different locations of the magnetic field with respect to the core inflow and no-flow boundaries are considered. Effects of the magnetic field on the physical variables such as saturation, nanoparticles concentrations, porosity, and permeability have been examined. The important result from this simulation is that the magnetic source location has a significant influence on the physical variables of the model. Based on the flow direction and the location of the magnet, the magnetic field can assist or oppose the flow of this two-phase system. Such observed behavior can be used for different applications depending on whether the goal is to aid or delay the injected fluid advancement. According to this investigation, if the magnet is placed next to the no-flow boundary of the studied countercurrent numerical experiment, the saturation of nanoparticles-water suspension increases. Moreover, the concentration of nanoparticles is observed to decrease slightly, which has been attributed to the slight increase of the deposition of nanoparticles. On the other hand, when the magnet is placed next to the inflow/outflow boundary, the magnet resists the flow of the ferrofluid suspension and decreases water invasion capacity, thereby. Furthermore, the nanoparticles concentration seems to increase under the effect of the magnetic field, and the deposited nanoparticles concentration decreases. Both porosity and permeability are reduced due to the nanoparticles adherence to the walls of the porous media.



---

*Source: 3602593-2017-08-28.xml*