Ching-Yao Chen

Miscible displacements in capillary tubes. Part 2. Numerical simulations

Numerical simulations are presented which, in conjunction with the accompanying experimental investigation by Petitjeans & Maxworthy (1996), are intended to elucidate the miscible flow that is generated if a fluid of given viscosity and density displaces a second fluid of different such properties in a capillary tube or plane channel. The global features of the flow, such as the fraction of the displaced fluid left behind on the tube walls, are largely controlled by dimensionless quantities in the form of a Peclet number Pe , an Atwood number At , and a gravity parameter. However, further dimensionless parameters that arise from the dependence on the concentration of various physical properties, such as viscosity and the diffusion coefficient, result in significant effects as well. The simulations identify two distinct Pe regimes, separated by a transitional region. For large values of Pe , typically above O (10), a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the strength of the gravitational forces, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip and others with toroidal recirculation regions inside the finger. Simulations that account for the experimentally observed dependence of the diffusion coefficient on the concentration show the evolution of fingers that combine steep external concentration layers with smooth concentration fields on the inside. In the small- Pe regime, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). An attempt was made to evaluate the importance of the Korteweg stresses and the consequences of assuming a divergence-free velocity field. Scaling arguments indicate that these effects should be strongest when steep concentration fronts exist, i.e. at large values of Pe and At. However, when compared to the viscous stresses, Korteweg stresses may be relatively more important at lower values of these parameters, and we cannot exclude the possibility that minor discrepancies observed between simulations and experiments in these parameter regimes are partially due to these extra stresses.

Miscible porous media displacements in the quarter five-spot configuration. Part 1. The homogeneous case

A detailed two-part computational investigation is conducted into the dynamical evolution of two-dimensional miscible porous media flows in the quarter five-spot arrangement of injection and production wells. High-accuracy direct numerical simulations are performed that reproduce all dynamically relevant length scales in solving the vorticity–streamfunction formulation of Darcys law. The accuracy of the method is demonstrated by a comparison of simulation data with linear stability results for radial source flow. Within this part, Part 1 of the present investigation, a series of simulations is discussed that demonstrate how the mobility ratio and the dimensionless flow rate denoted by the Peclet number Pe affect both local and integral features of homogeneous displacement processes. Mobility ratios up to 150 and Pe -values up to 2000 are investigated. For sufficiently large Pe -values, the flow near the injection well gives rise to a vigorous viscous fingering instability. As the unstable concentration front approaches the central region of the domain, nonlinear interactions between the fingers similar to those known from unidirectional flows are observed, such as merging, partial merging, and shielding, along with secondary tip-splitting and side-branching instabilities. At large Pe -values, several of these fingers compete for long times, before one of them accelerates ahead of the others and leads to the breakthrough of the front. In contrast to unidirectional flows, the quarter five-spot geometry imposes both an external length scale and a time scale on the flow. The resulting spatial non-uniformity of the potential base flow is observed to lead to a clear separation in space and time of large and small scales in the flow. Small scales occur predominantly during the early stages near the injection well, and at late times near the production well. The central domain is dominated by larger scales. Taken together, the results of the simulations demonstrate that both the mobility ratio and Pe strongly affect the dynamics of the flow. While some integral measures, such as the recovery at breakthrough, may show only a weak dependence on Pe for large Pe -values, the local fingering dynamics continue to change with Pe . The increased susceptibility of the flow to perturbations during the early stages provides the motivation to formulate an optimization problem that attempts to maximize recovery, for a constant overall process time, by employing a time-dependent flow rate. Within the present framework, which accounts for molecular diffusion but not for velocity-dependent dispersion, simulation results indeed indicate the potential to increase recovery by reducing the flow rate at early times, and then increasing it during the later stages.

Miscible porous media displacements in the quarter five-spot configuration. Part 2. Effect of heterogeneities

Direct numerical simulations are employed to investigate the coupling between the viscous fingering instability and permeability heterogeneities for miscible displacements in quarter five-spot flows. Even moderate inhomogeneities are seen to have a strong effect on the flow, which can result in a complete bypass of the linear growth phase of the viscous fingering instability. In contrast to their homogeneous counterparts (cf. Part 1, Chen & Meiburg 1998), heterogeneous quarter five-spot flows are seen to exhibit a more uniform dominant length scale throughout the entire flow domain. In line with earlier findings for unidirectional displacements, an optimal interaction of the mobility and permeability related vorticity modes can occur when the viscous length scale is of the same order as the correlation length of the heterogeneities. This resonance mechanism results in a minimal breakthrough recovery for intermediate correlation lengths, at fixed dimensionless flow rates in the form of a Peclet number Pe . However, for a constant correlation length, the recovery does not show a minimum as Pe is varied. Confirming earlier observations, the simulations show a more rapid breakthrough as the variance of the permeability variations increases. However, this tendency is far more noticeable in some parameter regimes than in others. It is furthermore observed that relatively low variances usually cannot change the tendency for a dominant finger to evolve along the inherently preferred diagonal direction, especially for relatively small correlation lengths. Only for higher variances, and for larger correlation lengths, are situations observed in which an off-diagonal finger can become dominant. Due to the nonlinear nature of the selection mechanisms at work, a change in the variance of the heterogeneities can result in the formation of dominant fingers along entirely different channels.

Miscible quarter five-spot displacements in a Hele-Shaw cell and the role of flow-induced dispersion

Miscible quarter five-spot displacements in a Hele-Shaw cell are investigated by means of experimental measurements and numerical simulations. The experiments record both the volumetric as well as the surface efficiency at breakthrough as a function of the dimensionless flow rate in the form of a Peclet number and the viscosity contrast. For small flow rates, both of these efficiency measures decrease uniformly with increasing Peclet numbers. At large flow rates, an asymptotic state is reached where the efficiencies no longer depend on the Peclet number. Up to Atwood numbers of approximately 0.5, the less viscous fluid occupies close to 23 of the gap width, which indicates a near-parabolic velocity profile across the gap. Consequently, in this parameter range a Taylor dispersion approach should be well suited to account for flow-induced dispersion effects. For larger viscosity contrasts, accompanying two-dimensional numerical simulations based on Taylor dispersion predict an increased stabilization for hi...

Miscible displacements in capillary tubes: Influence of Korteweg stresses and divergence effects

The question is addressed as to whether Korteweg stresses and/or divergence effects can potentially account for discrepancies observed between conventional Stokes flow simulations (Chen and Meiburg) and experiments (Petitjeans and Maxworthy) for miscible flows in capillary tubes. An estimate of the vorticity and stream function fields induced by the Korteweg stresses is presented, which shows these stresses to result in the formation of a vortex ring structure near the tip of the concentration front. Through this mechanism the propagation velocity of the concentration front is reduced, in agreement with the experimental observations. Divergence effects, on the other hand, are seen to be very small, and they have a negligible influence on the tip velocity. As a result, we can conclude that they are not responsible for the discrepancies between experiments and conventional Stokes simulations.

Miscible droplets in a porous medium and the effects of Korteweg stresses

Numerical simulation results are presented for the displacement of a drop in a porous medium. The drop is surrounded by a more viscous fluid with which it is fully miscible. The simulations are based on a set of augmented Hele–Shaw equations that account for nonconventional, so-called Korteweg stresses resulting from locally steep concentration gradients. Globally, these stresses tend to stabilize the displacement. However, there are important distinctions between their action and the effects of surface tension in an immiscible flow. Since the Korteweg stresses depend on the concentration gradient field, the effective net force across the miscible interface region is not just a function of the drop’s geometry, but also of the velocity gradient tensor. Locally high strain at the leading edge of the drop generates steep concentration gradients and large Korteweg stresses. Around the rear of the drop, the diffusion layer is much thicker and the related stresses smaller. The drop is seen to form a tail, which...

An experimental study on Rosensweig instability of a ferrofluid droplet

We experimentally investigate the interfacial morphologies of Rosensweig instability on an extremely thin layer of ferrofluid droplets under a constant perpendicular magnetic field. Striking patterns consisting of numerous subscale droplets that developed from Rosensweig instability are observed. For a dry plate, on which surface tension dominates, the breaking pattern of subscale droplets can be characterized by a dimensionless magnetic Bond number Bom. In general, a more pronounced instability, which is evident by a greater number of breaking subscale droplets N, arises with a higher Bom. For a magnetic Bond number that is larger than a critical value, we identify a new mode of interfacial breakup pattern, where the central droplet is torn apart with major mass loss. In addition, we found that the volume fractions of breaking subscale droplets are strongly affected by the height variation of the initial fluid surface and appear unevenly distributed with dominance of a central droplet. On the other hand,...

Flow visualization of natural convection of magnetic fluid in a rectangular Hele-Shaw cell

Abstract The nature convection of a magnetic fluid in a Hele-Shaw cell with aspect ratio of one is studied experimentally. Results obtained from heat transfer measurements and shadowgraphs revealed that the vertically imposed magnetic field has a destabilizing influence. The flow instability mode becomes different from that without the magnetic field.

Numerical simulations of fingering instabilities in miscible magnetic fluids in a Hele-Shaw cell and the effects of Korteweg stresses

The first high-accuracy direct numerical simulations of labyrinthine instabilities of miscible magnetic fluids are presented. Great similarities of fingering patterns in a miscible magnetic droplet with the effects of Korteweg stresses to their immiscible equivalents give strong qualitative validations to the numerical results. For conventional modeling that excludes the effects of Korteweg stresses, more vigorous fingerings are observed at stronger magnetic fields and higher viscosity contrast, which are in line with common expectations.

Single-layered multi-color electrowetting display by using ink-jet-printing technology and fluid-motion prediction with simulation

— This paper describes a single-layered multi-color electrowetting display (EWD) by using ink-jet-printing (IJP) technology and comparing different pattern electrodes with the use of the numerical investigations of ANSYS FLUENT®. This work consists of two parts: the first describes the design of implementing a single-layered multi-color EWD and the second demonstrates the application of ANSYS FLUENT® simulation in different pattern electrodes settings on the proposed EWD. The single-layered multi-color EW device was evaluated by using various colored oils without adopting a color filter. The single-layered multi-color EWD at a driving voltage of 25 V can achieve a maximum aperture ratio and reflectivity of 80% and 38.5%, respectively. The colored saturation of R, G, B oils can increase to 50% (NTSC: 13.3–27.8%). In addition, a radiate electrode at the required viewable area condition of 85% and force 5 * Fk, which results in ink stable contraction and a shorter response time of 50% (radiate vs. square), was proposed. The experimental results and simulation demonstrate that ink-jet-printing (IJP) technology along with the use of radiate electrodes can result in a single-layered multi-color EWD with a shorter response time.