As is well known, the extrusion rate of polymers from a cylindrical tube or slit (a ``die) is in practice limited by the appearance of ``melt fracture instabilities which give rise to unwanted distortions or even fracture of the extrudate. We present the results of a weakly nonlinear analysis which gives evidence for an intrinsic generic route to melt fracture via a weakly nonlinear subcritical instability of viscoelastic Poiseuille flow. This instability and the onset of associated melt fracture phenomena appear at a fixed ratio of the elastic stresses to viscous stresses of the polymer solutionte
A modal stability analysis shows that plane Poiseuille flow of an Oldroyd-B fluid becomes unstable to a `center mode with phase speed close to the maximum base-flow velocity, $U_{max}$. The governing dimensionless groups are the Reynolds number $Re = rho U_{max} H/eta$, the elasticity number $E = lambda eta/(H^2rho)$, and the ratio of solvent to solution viscosity $eta_s/eta$; here, $lambda$ is the polymer relaxation time, $H$ is the channel half-width, and $rho$ is the fluid density. For experimentally relevant values (e.g., $E sim 0.1$ and $beta sim 0.9$), the predicted critical Reynolds number, $Re_c$, for the center-mode instability is around $200$, with the associated eigenmodes being spread out across the channel. In the asymptotic limit of $E(1 -beta) ll 1$, with $E$ fixed, corresponding to strongly elastic dilute polymer solutions, $Re_c propto (E(1-beta))^{-frac{3}{2}}$ and the critical wavenumber $k_c propto (E(1-beta))^{-frac{1}{2}}$. The unstable eigenmode in this limit is confined in a thin layer near the channel centerline. The above features are largely analogous to the center-mode instability in viscoelastic pipe flow (Garg et al., Phys. Rev. Lett., 121, 024502 (2018)), and suggest a universal linear mechanism underlying the onset of turbulence in both channel and pipe flows of suffciently elastic dilute polymer solutions.
Recently, detailed experiments on visco-elastic channel flow have provided convincing evidence for a nonlinear instability scenario which we had argued for based on calculations for visco-elastic Couette flow. Motivated by these experiments we extend the previous calculations to the case of visco-elastic Poiseuille flow, using the Oldroyd-B constitutive model. Our results confirm that the subcritical instability scenario is similar for both types of flow, and that the nonlinear transition occurs for Weissenberg numbers somewhat larger than one. We provide detailed results for the convergence of our expansion and for the spatial structure of the mode that drives the instability. This also gives insight into possible similarities with the mechanism of the transition to turbulence in Newtonian pipe flow.
A necessary condition for the validity of the linear viscoelastic model for a (passive) polymeric cylinder with an ultrasonic hysteresis-type absorption submerged in a non-viscous fluid requires that the absorption efficiency is positive (Qabs > 0) satisfying the law of the conservation of energy. This condition imposes restrictions on the values attributed to the normalized absorption coefficients for the compressional and shear-wave wavenumbers for each partial-wave mode n. The forbidden values produce negative axial radiation force, absorption and extinction efficiencies, as well as an enhancement of the scattering efficiency, not in agreement with the conservation of energy law. Numerical results for the radiation force, extinction, absorption and scattering efficiencies are performed for three viscoelastic (VE) polymer cylinders immersed in a non-viscous host liquid (i.e. water) with particular emphasis on the shear-wave absorption coefficient of the cylinder, the dimensionless size parameter and the partial-wave mode number n. Mathematical constraints are established for the non-dimensional absorption coefficients of the longitudinal and shear waves for a cylinder (i.e. 2D case) and a sphere (i.e. 3D case) in terms of the sound velocities in the VE material. The analysis suggests that the domain of validity for any viscoelastic model describing acoustic attenuation inside a lossy cylinder (or sphere) in a non-viscous fluid must be verified based upon the optical theorem.
Hydrodynamic flow of charge carriers in graphene is an energy flow unlike the usual mass flow in conventional fluids. In neutral graphene, the energy flow is decoupled from the electric current, making it difficult to observe the hydrodynamic effects and measure the viscosity of the electronic fluid by means of electric current measurements. In particular, we show that the hallmark Poiseuille flow in a narrow channel cannot be driven by the electric field irrespective of boundary conditions at the channel edges. Nevertheless one can observe nonuniform current densities similarly to the case of the well-known ballistic-diffusive crossover. The standard diffusive behavior with the uniform current density across the channel is achieved under the assumptions of specular scattering on the channel boundaries. This flow can also be made nonuniform by applying weak magnetic fields. In this case, the curvature of the current density profile is determined by the quasiparticle recombination processes dominated by the disorder-assisted electron-phonon scattering -- the so-called supercollisions.
The main purpose of this work is to simulate two-phase flow in the form of immiscible displacement through anisotropic, three-dimensional (3D) discrete fracture networks (DFN). The considered DFNs are artificially generated, based on a general distribution function or are conditioned on measured data from deep geological investigations. We introduce several modifications to the invasion percolation (MIP) to incorporate fracture inclinations, intersection lines, as well as the hydraulic path length inside the fractures. Additionally a trapping algorithm is implemented that forbids any advance of the invading fluid into a region, where the defending fluid is completely encircled by the invader and has no escape route. We study invasion, saturation, and flow through artificial fracture networks, with varying anisotropy and size and finally compare our findings to well studied, conditioned fracture networks.