Research on the Mesoscopic Characteristics of Kelvin-Helmholtz Instability in Polymer Fluids with Dissipative Particle Dynamics

In this paper, the two-dimensional Kelvin-Helmholtz (KH) instability occurring in the shear flow of polymer fluids is modeled by the dissipative particle dynamics (DPD) method at the coarse-grained molecular level. A revised FENE model is proposed to properly describe the polymer chains. In this revised model, the elastic repulsion and tension are both considered between the adjacent beads, the bond length of which is set as one segment's equilibrium length. The entanglements between polymer chains are described with a bead repulsive potential. The characteristics of such a KH instability in polymer fluid shear flow can be successfully captured in the simulations by the use of the modified FENE model. The numerical results show that the waves and vortexes grow more slowly in the shear flow of the polymer fluids than in the Newtonian fluid case, these vortexes become flat, and the polymer impedes the mixing of fluids and inhibits the generation of turbulence. The effects of the polymer concentration, chain length, and extensibility are also investigated regarding the evolution of KH instability. It is shown that the mixing of two polymer fluids reduces, and the KH instability becomes more suppressed as the polymer concentration increases. The vortexes become much longer with the evolution of the elongated interface as the chain length turns longer. As the extensibility increases, the vortexes become more flattened. Moreover, the roll-up process is significantly suppressed if the polymer has sufficiently high extensibility. These observations show that the polymer and its properties significantly influence the formation and evolution of the coherent structures such as the waves and vortexes in the KH instability progress.