Two‐phase fluid motion through porous medium with volume fraction: An application of MATLAB bvp4c solver technique

Heat Transfer ◽  
2021 ◽  
Author(s):  
Debasish Dey ◽  
Barbie Chutia
Ultrasonics ◽  
2017 ◽  
Vol 81 ◽  
pp. 73-80
Author(s):  
Qingbang Han ◽  
Lihua Qi ◽  
Minglei Shan ◽  
Cheng Yin ◽  
Xueping Jiang ◽  
...  

Author(s):  
Aureliano Sancho S. Paiva ◽  
Rafael S. Oliveira ◽  
Roberto F. S. Andrade

We investigate how a plug of obstacles inside a two-dimensional channel affects the drainage of high viscous fluid (oil) when the channel is invaded by a less viscous fluid (water). The plug consists of an Apollonian packing with, at most, 17 circles of different sizes, which is intended to model an inhomogeneous porous region. The work aims to quantify the amount of retained oil in the region where the flow is influenced by the packing. The investigation, carried out with the help of the computational fluid dynamics package ANSYS-FLUENT , is based on the integration of the complete set of equations of motion. The study considers the effect of both the injection speed and the number and size of obstacles, which directly affects the porosity of the system. The results indicate a complex dependence in the fraction of retained oil on the velocity and geometric parameters. The regions where the oil remains trapped is very sensitive to the number of circles and their size, which influence in different ways the porosity of the system. Nevertheless, at low values of Reynolds and capillary numbers Re <4 and n c ≃10 −5 , the overall expected result that the volume fraction of oil retained decreases with increasing porosity is recovered. A direct relationship between the injection speed and the fraction of oil is also obtained.


2020 ◽  
Vol 24 (2 Part B) ◽  
pp. 1045-1054 ◽  
Author(s):  
Mehdi Ahmadi ◽  
Farsani Khosravi

In this paper, the numerical solution of non-Newtonian two-phase fluid-flow through a channel with a cavity was studied. Carreau-Yasuda non-Newtonian model which represents well the dependence of stress on shear rate was used and the effect of n index of the model and the effect of input Reynolds on the attribution of flow were considered. Governing equations were discretized using the finite volume method on staggered grid and the form of allocating flow parameters on staggered grid is based on marker and cell method. The QUICK scheme is employed for the convection terms in the momentum equations, also the convection term is discretized by using the hybrid upwind-central scheme. In order to increase the accuracy of making discrete, second order Van Leer accuracy method was used. For mixed solution of velocity-pressure field SIMPLEC algorithm was used and for pressure correction equation iteratively line-by-line TDMA solution procedure and the strongly implicit procedure was used. As the results show, by increasing Reynolds number, the time of sweeping the non-Newtonian fluid inside the cavity decreases, the velocity of Newtonian fluid increases and the pressure decreases. In the second section, by increasing n index, the velocity increases and the volume fraction of non-Newtonian fluid after cavity increases and by increasing velocity, the pressure decreases. Also changes in the velocity, pressure and volume fraction of fluids inside the channel and cavity are more sensible to changing the Reynolds number instead of changing n index.


Sign in / Sign up

Export Citation Format

Share Document