Interplay of bulk soluble surfactants and interfacial kinetics governs the stability of two-layer channel flows.

The linear stability of two-layer channel flows in the presence of bulk-soluble surfactants is investigated here, taking into account the rheological properties of the interface. The interfacial stresses are quantified using the Boussinesq-Scriven model, while the surfactant kinetics is assumed to follow the Frumkin isotherm, which accounts for their non-ideal behavior. Our results show that in general, the bulk solubility of surfactants has a stabilizing effect on the interface, both with and without the presence of inertia. On the other hand, the interfacial viscosities play a more complex role, depending on the viscosity ratios of the two fluids, the thickness of the fluid layers, the strength of the surface tension gradients, and the extent of inertia. We show that depending on the strength of inertia and the variability in the surface tension, the interfacial rheology may either stabilize or destabilize the base flow. However, for sufficiently small Reynolds numbers, the surface viscosity always has a stabilizing influence. Our results may be used to better design stable co-flow systems with applications in various processes such as surface coating, preparation of fluid lenses, as well as in a host of multi-purpose microfluidic devices.

The breakdown of Darcy's law in a soft porous material.

We perform direct numerical simulations of the flow through a model of deformable porous medium. Our model is a two-dimensional hexagonal lattice, with defects, of soft elastic cylindrical pillars, with elastic shear modulus G, immersed in a liquid. We use a two-phase approach: the liquid phase is a viscous fluid and the solid phase is modeled as an incompressible viscoelastic material, whose complete nonlinear structural response is considered. We observe that the Darcy flux (q) is a nonlinear function - steeper than linear - of the pressure-difference (ΔP) across the medium. Furthermore, the flux is larger for a softer medium (smaller G). We construct a theory of this super-linear behavior by modelling the channels between the solid cylinders as elastic channels whose walls are made of material with a linear constitutive relation but can undergo large deformation. Our theory further predicts that the flow permeability is an universal function of ΔP/G, which is confirmed by the present simulations.

Confinement effects in premelting dynamics.

We examine the effects of confinement on the dynamics of premelted films driven by thermomolecular pressure gradients. Our approach is to modify a well-studied setting in which the thermomolecular pressure gradient is driven by a temperature gradient parallel to an interfacially premelted elastic wall. The modification treats the increase in viscosity associated with the thinning of films, studied in a wide variety of materials, using a power law and we examine the consequent evolution of the confining elastic wall. We treat (1) a range of interactions that are known to underlie interfacial premelting and (2) a constant temperature gradient wherein the thermomolecular pressure gradient is a constant. The difference between the cases with and without the proximity effect arises in the volume flux of premelted liquid. The proximity effect increases the viscosity as the film thickness decreases thereby requiring the thermomolecular pressure driven flux to be accommodated at higher temperatures where the premelted film thickness is the largest. Implications for experiment and observations of frost heave are discussed.

Fingering instability and mixing of a blob in porous media.

The curvature of the unstable part of the miscible interface between a circular blob and the ambient fluid in two-dimensional homogeneous porous media depends on the viscosity of the fluids. The influence of the interface curvature on the fingering instability and mixing of a miscible blob within a rectilinear displacement is investigated numerically. The fluid velocity in porous media is governed by Darcy's law, coupled with a convection-diffusion equation that determines the evolution of the solute concentration controlling the viscosity of the fluids. Numerical simulations are performed using a Fourier pseudospectral method to determine the dynamics of a miscible blob (circular or square). It is shown that for a less viscous circular blob, there exist three different instability regions without any finite R-window for viscous fingering, unlike the case of a more viscous circular blob. Critical blob radius for the onset of instability is smaller for a less viscous blob as compared to its more viscous counterpart. Fingering enhances spreading and mixing of miscible fluids. Hence a less viscous blob mixes with the ambient fluid quicker than the more viscous one. Furthermore, we show that mixing increases with the viscosity contrast for a less viscous blob, while for a more viscous one mixing depends nonmonotonically on the viscosity contrast. For a more viscous blob mixing depends nonmonotonically on the dispersion anisotropy, while it decreases monotonically with the anisotropic dispersion coefficient for a less viscous blob. We also show that the dynamics of a more viscous square blob is qualitatively similar to that of a circular one, except the existence of the lump-shaped instability region in the R-Pe plane. We have shown that the Rayleigh-Taylor instability in a circular blob (heavier or lighter than the ambient fluid) is independent of the interface curvature.

Nonmodal linear stability analysis of miscible viscous fingering in porous media.

The nonmodal linear stability of miscible viscous fingering in a two-dimensional homogeneous porous medium has been investigated. The linearized perturbed equations for Darcy's law coupled with a convection-diffusion equation is discretized using a finite difference method. The resultant initial value problem is solved by a fourth-order Runge-Kutta method, followed by a singular value decomposition of the propagator matrix. Particular attention is given to the transient behavior rather than the long-time behavior of eigenmodes predicted by the traditional modal analysis. The transient behaviors of the response to external excitations and the response to initial conditions are studied by examining the ε-pseudospectra structures and the largest energy growth function, respectively. With the help of nonmodal stability analysis we demonstrate that at early times the displacement flow is dominated by diffusion and the perturbations decay. At later times, when convection dominates diffusion, perturbations grow. Furthermore, we show that the dominant perturbation that experiences the maximum amplification within the linear regime lead to the transient growth. These two important features were previously unattainable in the existing linear stability methods for miscible viscous fingering. To explore the relevance of the optimal perturbation obtained from nonmodal analysis, we performed direct numerical simulations using a highly accurate pseudospectral method. Furthermore, a comparison of the present stability analysis with existing modal and initial value approach is also presented. It is shown that the nonmodal stability results are in better agreement than the other existing stability analyses, with those obtained from direct numerical simulations.

Onset of fingering instability in a finite slice of adsorbed solute.

The effect of a linear adsorption isotherm on the onset of fingering instability in a miscible displacement in the application of liquid chromatography, pollutant contamination in aquifers, etc., is investigated. Such fingering instability on the solute dynamics arise due to the miscible viscus fingering (VF) between the displacing fluid and sample solvent. We use a Fourier pseudo-spectral method to solve the initial value problem that appears in the linear stability analysis. The present linear stability analysis is of generic type and it captures the early-time-diffusion-dominated region which was never expressible through the quasi-steady-state analysis (QSSA). In addition, it measures the onset of instability more accurately than the QSSA methods. It is shown that the onset time depends nonmonotonically on the retention parameter of the solute adsorption. This qualitative influence of the retention parameter on the onset of instability resemblances with the results obtained from direct numerical simulations of the nonlinear equations. Moreover, the present linear stability method helps for an appropriate characterization of the linear and nonlinear regimes of miscible VF instability and also can be useful for the fluid flow problems with the unsteady base state.

Effect of Péclet number on miscible rectilinear displacement in a Hele-Shaw cell.

The influence of fluid dispersion on the Saffman-Taylor instability in miscible fluids has been investigated in both the linear and the nonlinear regimes. The convective characteristic scales are used for the dimensionless formulation that incorporates the Péclet number (Pe) into the governing equations as a measure for the fluid dispersion. A linear stability analysis (LSA) is performed in a similarity transformation domain using the quasi-steady-state approximation. LSA results confirm that a flow with a large Pe has a higher growth rate than a flow with a small Pe. The critical Péclet number (Pec) for the onset of instability for all possible wave numbers and also a power-law relation of the onset time and most unstable wave number with Pe are observed. Unlike the radial source flow, Pec is found to vary with t0. A Fourier spectral method is used for direct numerical simulations (DNS) of the fully nonlinear system. The power-law dependence of the onset of instability ton on Pe is obtained from the DNS and found to be inversely proportional to Pe and it is in good agreement with that obtained from the LSA. The influence of the anisotropic dispersion is analyzed in both the linear and the nonlinear regimes. The results obtained confirm that for a weak transverse dispersion merging happens slowly and hence the small wave perturbations become unstable. We also observ that the onset of instability sets in early when the transverse dispersion is weak and varies with the anisotropic dispersion coefficient, ε, as ∼√[ε], in compliance with the LSA. The combined effect of the Korteweg stress and Pe in the linear regime is pursued. It is observed that, depending on various flow parameters, a fluid system with a larger Pe exhibits a lower instantaneous growth rate than a system with a smaller Pe, which is contrary to the results when such stresses are absent.