Search for: All records

The flow in a Hele-Shaw cell with a time-increasing gap poses a unique shrinking interface problem. When the upper plate of the cell is lifted perpendicularly at a prescribed speed, the exterior less viscous fluid penetrates the interior more viscous fluid, which generates complex, time-dependent interfacial patterns through the Saffman–Taylor instability. The pattern formation process sensitively depends on the lifting speed and is still not fully understood. For some lifting speeds, such as linear or exponential speed, the instability is transient and the interface eventually shrinks as a circle. However, linear stability analysis suggests there exist shape invariant shrinking patterns if the gap $b(t)$ is increased more rapidly: $b(t)=\left (1-({7}/{2})\tau \mathcal {C} t\right )^{-{2}/{7}}$ , where $\tau$ is the surface tension and $\mathcal {C}$ is a function of the interface perturbation mode $k$ . Here, we use a spectrally accurate boundary integral method together with an efficient time adaptive rescaling scheme, which for the first time makes it possible to explore the nonlinear limiting dynamical behaviour of a vanishing interface. When the gap is increased at a constant rate, our numerical results quantitatively agree with experimental observations (Nase et al. , Phys. Fluids , vol. 23, 2011, 123101). When we use the shape invariant gap $b(t)$ , our nonlinear results reveal the existence of $k$ -fold dominant, one-dimensional, web-like networks, where the fractal dimension is reduced to almost unity at late times. We conclude by constructing a morphology diagram for pattern selection that relates the dominant mode $k$ of the vanishing interface and the control parameter $\mathcal {C}$ . 

We present a quasi-incompressible Navier–Stokes–Cahn–Hilliard (q-NSCH) diffuse interface model for two-phase fluid flows with variable physical properties that maintains thermodynamic consistency. Then, we couple the diffuse domain method with this two-phase fluid model – yielding a new q-NSCH-DD model – to simulate the two-phase flows with moving contact lines in complex geometries. The original complex domain is extended to a larger regular domain, usually a cuboid, and the complex domain boundary is replaced by an interfacial region with finite thickness. A phase-field function is introduced to approximate the characteristic function of the original domain of interest. The original fluid model, q-NSCH, is reformulated on the larger domain with additional source terms that approximate the boundary conditions on the solid surface. We show that the q-NSCH-DD system converges to the q-NSCH system asymptotically as the thickness of the diffuse domain interface introduced by the phase-field function shrinks to zero ( $\epsilon \rightarrow 0$ ) with $\mathcal {O}(\epsilon )$ . Our analytic results are confirmed numerically by measuring the errors in both $L^{2}$ and $L^{\infty }$ norms. In addition, we show that the q-NSCH-DD system not only allows the contact line to move on curved boundaries, but also makes the fluid–fluid interface intersect the solid object at an angle that is consistent with the prescribed contact angle. 

Preparation of thin films by dissolving polymers in a common solvent followed by evaporation of the solvent has become a routine processing procedure. However, modeling of thin film formation in an evaporating solvent has been challenging due to a need to simulate processes at multiple length and time scales. In this work, we present a methodology based on the principles of linear non-equilibrium thermodynamics, which allows systematic study of various effects such as the changes in the solvent properties due to phase transformation from liquid to vapor and polymer thermodynamics resulting from such solvent transformations. The methodology allows for the derivation of evaporative flux and boundary conditions near each surface for simulations of systems close to the equilibrium. We apply it to study thin film microstructural evolution in phase segregating polymer blends dissolved in a common volatile solvent and deposited on a planar substrate. Effects of the evaporation rates, interactions of the polymers with the underlying substrate and concentration dependent mobilities on the kinetics of thin film formation are studied.