Mesoscopic Analysis of Dynamic Droplet Behavior on Wetted Flat and Grooved Surface for Low Viscosity Ratio

2017 ◽  
Vol 139 (5) ◽  
Author(s):  
Saurabh Bhardwaj ◽  
Amaresh Dalal

In the present study, the interfacial dynamics of displacement of three-dimensional spherical droplet on a rectangular microchannel wall considering wetting effects are studied. The two-phase lattice Boltzmann Shan−Chen model is used to explore the physics. The main focus of this study is to analyze the effect of wettability, low viscosity ratio, and capillary number on the displacement of spherical droplet subjected to gravitational force on flat as well as grooved surface of the channel wall. The hydrophobic and hydrophilic natures of wettabilities on wall surface are considered to study for viscosity ratio, M≤1. The results are presented in the form of temporal evolution of wetted length and wetted area for combined viscosity ratios and wettability scenario. In the present study, it is observed that in dynamic droplet displacement, the viscosity ratio and the capillary number play a significant role. It is found that as the viscosity ratio increases, both the wetted area and the wetted length increase and decrease in the case of hydrophilic and hydrophobic wettable walls, respectively. The groove area on the vertical wall tries to entrap fraction of droplet fluid in case of hydrophilic surface of the vertical wall, whereas in hydrophobic case, droplet moves past the groove without entrapment.

Author(s):  
Saurabh Bhardwaj ◽  
Amaresh Dalal

In the present study, the interfacial dynamics of displacement of three dimensional spherical droplet on a rectangular microchannel wall considering wetting effects are studied. The two-phase lattice Boltzmann Shan-Chen model is used to explore the physics. The main focus of this study is to analyse the effect of wettability, low viscosity ratio and capillary number on the displacement of spherical droplet subjected to gravitational force. The hydrophobic and hydrophilic nature of wettabilities on wall surface are considered to study with capillary number, Ca=0.1, 0.35 and 0.66 and viscosity ratio, M ≤ 1. The results are presented in the form of temporal evolution of wetted length and wetted area for combined viscosity ratios and wettability scenario. In the present study, it is observed that in dynamic droplet displacement, the viscosity ratio and capillary number play a significant role. It is found that as viscosity ratio increases, both the wetted area and the wetted length increase and decrease in the case of hydrophilic and hydrophobic wettable wall respectively.


2020 ◽  
Vol 17 (162) ◽  
pp. 20190792 ◽  
Author(s):  
Candido Diaz ◽  
Daniel Maksuta ◽  
Gaurav Amarpuri ◽  
Akio Tanikawa ◽  
Tadashi Miyashita ◽  
...  

Contaminants decrease adhesive strength by interfering with substrate contact. Spider webs adhering to moths present an ideal model to investigate how natural adhesives overcome contamination because moths' sacrificial layer of scales rubs off on sticky silk, facilitating escape. However, Cyrtarachninae spiders have evolved gluey capture threads that adhere well to moths. Cyrtarachne capture threads contain large glue droplets oversaturated with water, readily flowing but also prone to drying out. Here, we compare the spreading and adhesion of Cyrtarachne akirai glue on intact mothwings, denuded cuticle and glass to the glue of a common orb-weaving spider, Larinioides cornutus, to understand how C. akirai glue overcomes dirty surfaces. Videos show that C. akirai 's glue spreading accelerates along the underlying moth cuticle after the glue seeps beneath the moth scales—not seen on denuded cuticle or hydrophilic glass. Larinioides cornutus glue droplets failed to penetrate the moth scales, their force of adhesion thus limited by the strength of attachment of scales to the cuticle. The large size and low viscosity of C. akirai glue droplets function together to use the three-dimensional topography of the moth's scales against itself via capillary forces. Infrared spectroscopy shows C. akirai glue droplets readily lose free-flowing water. We hypothesize that this loss of water leads to increased viscosity during spreading, increasing cohesive forces during pull-off. This glue's two-phase behaviour shows how natural selection can leverage a defensive specialization of prey against themselves and highlights a new design principle for synthetic adhesives for adhering to troublesome surfaces.


1986 ◽  
Vol 164 ◽  
pp. 305-336 ◽  
Author(s):  
Madalena M. Dias ◽  
Alkiviades C. Payatakes

A theoretical simulator of immiscible displacement of a non-wetting fluid by a wetting one in a random porous medium is developed. The porous medium is modelled as a network of randomly sized unit cells of the constricted-tube type. Under creeping-flow conditions the problem is reduced to a system of linear equations, the solution of which gives the instantaneous pressures at the nodes and the corresponding flowrates through the unit cells. The pattern and rate of the displacement are obtained by assuming quasi-static flow and taking small time increments. The porous medium adopted for the simulations is a sandpack with porosity 0.395 and grain sizes in the range from 74 to 148 μrn. The effects of the capillary number, Ca, and the viscosity ratio, κ = μo/μw, are studied. The results confirm the importance of the capillary number for displacement, but they also show that for moderate and high Ca values the role of κ is pivotal. When the viscosity ratio is favourable (κ < 1), the microdisplacement efficiency begins to increase rapidly with increasing capillary number for Ca > 10−5, and becomes excellent as Ca → 10−3. On the other hand, when the viscosity ratio is unfavourable (κ > 1), the microdisplacement efficiency begins to improve only for Ca values larger than, say, 5 × 10−4, and is substantially inferior to that achieved with κ < 1 and the same Ca value. In addition to the residual saturation of the non-wetting fluid, the simulator predicts the time required for the displacement, the pattern of the transition zone, the size distribution of the entrapped ganglia, and the acceptance fraction as functions of Ca, κ, and the porous-medium geometry.


2019 ◽  
Vol 29 (5) ◽  
pp. 1563-1589 ◽  
Author(s):  
Navid Ahmadi Cheloii ◽  
Omid Ali Akbari ◽  
Davood Toghraie

Purpose This study aims to numerically investigate the heat transfer and laminar forced and two-phase flow of Water/Cu nanofluid in a rectangular microchannel with oblique ribs with angle of attacks equal to 0-45°. This simulation was conducted in the range of Reynolds numbers of 5-120 in volume fractions of 0, 2 and 4 per cent of solid nanoparticles in three-dimensional space. Design/methodology/approach This study investigates the effect of the changes of angle of attack of rectangular rib on heat transfer and hydrodynamics of two-phase flow. This study was done in three-dimensional space and simulation was done with finite volume method. SIMPLEC algorithm and second-order discretization of equations were used to increase the accuracy of results. The usage of nanofluid, application of rips with different angles of attacks and using the two-phase mixture method is the distinction of this paper compared with other studies. Findings The results of this research revealed that the changing angle of attack of ribs is an effective factor in heat transfer enhancement. On the other hand, the existence of rib on the internal surfaces of a microchannel increases friction coefficient. By increasing the volume fraction of nanoparticles, due to the augmentation of fluid density and viscosity, the pressure drop increases significantly. For all of the angle of attacks studied in this paper, the maximum rate of performance evaluation criterion has been obtained in Reynolds number of 30 and the minimum amount of performance evaluation criterion was been obtained in Reynolds numbers of 5 and 120. Originality/value Many studies have been done in the field of heat transfer in ribbed microchannel. In this paper, the laminar flow in the ribbed microchannel Water/Cu nanofluid in a rectangular microchannel by using two-phase mixture method is numerically investigated with different volume fractions (0-4 per cent), Reynolds numbers (5-120) and angle of attacks of rectangular rib in the indented microchannel (0-45°).


2021 ◽  
Author(s):  
Peter Mora ◽  
Gabriele Morra ◽  
Dave Yuen ◽  
Ruben Juanes

Abstract We present a suite of numerical simulations of two-phase flow through a 2D model of a porous medium using the Rothman-Keller Lattice Boltzmann Method to study the effect of viscous fingering on the recovery factor as a function of viscosity ratio and wetting angle. This suite involves simulations spanning wetting angles from non-wetting to perfectly wetting and viscosity ratios spanning from 0.01 through 100. Each simulation is initialized with a porous model that is fully saturated with a "blue" fluid, and a "red" fluid is then injected from the left. The simulation parameters are set such that the capillary number is 10, well above the threshold for viscous fingering, and with a Reynolds number of 0.2 which is well below the transition to turbulence and small enough such that inertial effects are negligible. Each simulation involves the "red" fluid being injected from the left at a constant rate such in accord with the specified capillary number and Reynolds number until the red fluid breaks through the right side of the model. As expected, the dominant effect is the viscosity ratio, with narrow tendrils (viscous fingering) occurring for small viscosity ratios with M ≪ 1, and an almost linear front occurring for viscosity ratios above unity. The wetting angle is found to have a more subtle and complicated role. For low wetting angles (highly wetting injected fluids), the finger morphology is more rounded whereas for high wetting angles, the fingers become narrow. The effect of wettability on saturation (recovery factor) is more complex than the expected increase in recovery factor as the wetting angle is decreased, with specific wetting angles at certain viscosity ratios that optimize yield. This complex phase space landscape with hills, valleys and ridges suggests the dynamics of flow has a complex relationship with the geometry of the medium and hydrodynamical parameters, and hence recovery factors. This kind of behavior potentially has immense significance to Enhanced Oil Recovery (EOR). For the case of low viscosity ratio, the flow after breakthrough is localized mainly through narrow fingers but these evolve and broaden and the saturation continues to increase albeit at a reduced rate. For this reason, the recovery factor continues to increase after breakthrough and approaches over 90% after 10 times the breakthrough time.


2005 ◽  
Vol 128 (1) ◽  
pp. 131-141 ◽  
Author(s):  
Mahmoud A. Ahmed

A theoretical analysis has been developed to predict the critical height and the location of the onset of gas entrainment during discharge from a stratified two-phase region through two oriented-side branches mounted on a vertical wall. In this analysis, a point sink model was first developed, followed by a more accurate three-dimensional finite branch model. The models are based on a new modified criterion for the onset of gas entrainment. The theoretically predicted critical height and the location of the onset of gas entrainment are found to be a function of the mass rate of each branch (Fr1 and Fr2), the distance between the centerlines of the two branches (L∕d), and the inclination angle (θ). The effects of these variables on the predicted critical height and the onset location were investigated. Furthermore, comparison between the theoretically predicted results and the available experimental data was carried out to verify the developed models. The comparison shows that the predicted results are very close to the measured data within a deviation percentage of 12% at Fr1>10. This small deviation percentage reflects a good agreement between the measured and predicted results.


1968 ◽  
Vol 8 (04) ◽  
pp. 325-330
Author(s):  
Alan D. Modine ◽  
Keith H. Coats

Abstract A mathematical model has been formulated for simulating three-dimensional displacement of a viscous fluid by a displacing fluid of zero viscosity. The model has been incorporated into a FORTRAN IV computer program for application in low-rate, high-permeability systems. Where applicable, the zero-viscosity program reduces computer time by a factor of 5 to 10 relative to conventional two- and three-dimensional programs. To determine the area of applicability, a gas-oil cross-section model representation of a high-dip, high-permeability reservoir was simulated with the zero-viscosity and conventional two-dimension programs for a range of flow rates up to 80 percent programs for a range of flow rates up to 80 percent of the critical rate. In comparing the two solutions, the conventional one was assumed to be the correct one because its program is based upon a more physically realistic model than that of the physically realistic model than that of the zero-viscosity solution. The two solutions agreed at rates up to 50 percent of the critical; at 80 percent they disagreed significantly. This indicates percent they disagreed significantly. This indicates that the zero-viscosity model, which is quite simple and inexpensive to apply, can be used with accuracy at rates up to at least 50 percent of the critical. This area of applicability is important in improving computational capability, for it is at these lower rates that the conventional programs are excessively costly. At the higher rates, where the zero-viscosity solution is not accurate, the conventional programs are easy and economical to apply. The zero-viscosity model accounts for capillary and gravitational forces, effects of viscosity and relative permeability for the displaced phase, and arbitrary reservoir heterogeneity. The program handles up to 1,800 blocks on an in-core basis. Introduction Computational difficulties caused by slow or metastable convergence in gas-oil calculations using conventional two-phase reservoir simulation programs have been correlatable with the effects of programs have been correlatable with the effects of low viscosity in the gas phase. In many such problems, a very small deviation in the calculated problems, a very small deviation in the calculated flow potentials causes a large deviation in the calculated gas flow due to the low viscosity. Thus, the program is trying to converge on a small variation in potential, which makes the computations difficult. A previous method of overcoming this difficulty has been to introduce in the conventional two-phase calculations an artificial resistance to gas flow; this method causes a more significant variation in the calculated flow potential. This paper describes a new method for treatment of gas-oil problems in which a zero-viscosity gas phase is used. Both methods are based on the assumption that oil mobility is the controlling factor in the displacement and that the behavior is insensitive to gas mobility over a relatively wide range. We show that the two methods give identical results, and since the correct gas mobility is bracketed by the two methods, we may conclude that either method gives valid results for low rate displacements. The chief advantage of the zero-viscosity program is lower computing costs. This report presents a mathematical description of the zero-viscosity model and compares saturation distributions calculated for several typical problems using the zero-viscosity and conventional two-phase programs. ZERO-VISCOSITY MODEL The zero-viscosity model simulates the immiscible displacement of a viscous fluid by a displacing fluid of zero viscosity. The method includes the effects of capillary and gravitational forces, relative permeability and viscosity in the displaced phase, permeability and viscosity in the displaced phase, and arbitrary reservoir heterogeneity. SPEJ P. 325


2014 ◽  
Vol 752 ◽  
pp. 383-409 ◽  
Author(s):  
Mitsuhiro Ohta ◽  
Yu Akama ◽  
Yutaka Yoshida ◽  
Mark Sussman

AbstractIn a low Morton number ($\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}M$) regime, the stability of a single drop rising in an immiscible viscous liquid is experimentally and computationally examined for varying viscosity ratio $\eta $ (the viscosity of the drop divided by that of the suspending fluid) and varying Eötvös number ($\mathit{Eo}$). Three-dimensional computations, rather than three-dimensional axisymmetric computations, are necessary since non-axisymmetric unstable drop behaviour is studied. The computations are performed using the sharp-interface coupled level-set and volume-of-fluid (CLSVOF) method in order to capture the deforming drop boundary. In the lower $\eta $ regimes, $\eta = 0.02 $ or 0.1, and when $\mathit{Eo}$ exceeds a critical threshold, it is observed that a rising drop exhibits nonlinear lateral/tilting motion. In the higher $\eta $ regimes, $\eta = 0.1$, 1.94, 10 or 100, and when $\mathit{Eo}$ exceeds another critical threshold, it is found that a rising drop becomes unstable and breaks up into multiple drops. The type of breakup, either ‘dumbbell’, ‘intermediate’ or ‘toroidal’, depends intimately on $\eta $ and $\mathit{Eo}$.


2010 ◽  
Vol 86 (1) ◽  
pp. 243-259 ◽  
Author(s):  
M. Ferer ◽  
Shelley L. Anna ◽  
Paul Tortora ◽  
J. R. Kadambi ◽  
M. Oliver ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document