Numerical simulation of interface tracking between two immiscible micropolar and dusty fluids

Abstract

The development of the interfacetracking between two immiscible liquids is one of the most challenging problems in fluid dynamics. The simulation of the lubrication process, oil extraction, and the physicochemical mechanism that govern the behaviour of the uncertain conduct of interface is a key element in “gas-oil” or “oil–water” natural gas-oil reservoirs systems and should be tracked for optimal recovery. Motivated by the interface evolution process, in this article, an unsteady flow of two immiscible Eringen micropolar and Dusty (fluid-particle suspension) fluids is considered through a horizontal channel and the functioning of the interface is analysed. Although the channel wall is hyper-stick and no-slip, it is believed that the fluid–fluid interface is unsteady and can also be deformed from one location to another; thus, the single momentum equation is coupled for monitoring the interface profiles using the VOF (Volume of fluid) technique. The modelled coupled- partial differential equations are numerically solved by the modified B-spline differential quadrature method. The interface profile under the influence of different parameters, i.e., Froude, Capillary, wave and Reynolds number, amplitude, the ratio of viscosities, and time is analysed. Three applied pressure gradients namely constant Ge = 10, periodic Ge = 10Sin(t), and decaying Ge = e-t are considered with various constant values of other hydrodynamic parameters. It's been noted that thetopologydevelops with time asthe pulsing pattern repeats quicker for a longer period before becoming stable. The flow's qualitative characteristics remain preserved, the interface begins to move vertically as the amplitude is elevated. A higher interfacial drift is seen in the applied periodic pressure gradient case as compared to the constant and decaying pressure gradient scenario. The tracking does take high time to be stable in the presence of the micropolar parameter with all three applied pressure gradient cases.