Popis: |
We study the flow of a dilute polymer solution in a wavy channel under steady-state flow conditions by employing the non-equilibrium thermodynamics two-fluid model (Mavrantzas-Beris, 1992), allowing for the coupling between polymer concentration and polymer stresses. The resulting highly complex system of partial differential equations describing inhomogeneous transport phenomena in the fluid are solved with an efficient implementation of the mixed finite-element method. We present numerical results for polymer concentration, stress, velocity and fluxes of polymer as a function of the non-dimensional parameters of the problem (the Deborah number , the Peclet number , the Reynolds number , the ratio of the solvent viscosity to the total fluid viscosity , and the constriction ratio of the channel width ). We find that the constricted part of the wall is depleted of polymer, when the polymer diffusion length-scale, expressed by the ratio of / , increases. The migration is more pronounced for macromolecules characterized by longer relaxation times, and takes place towards the expanding part of the channel or towards the centerplane. Migration is also enhanced by the width variability of the channel: the more corrugated the channel, the stronger the transfer of polymer to the centerplane. This increases the spatial extent of polymer depletion near the wall or induces a zone of sharp variation in polymer stress and concentration, which moves away from the channel wall, especially in lower polymer concentration. The development of a polymer-depleted layer smooths out the boundary layer which is known to arise with Boger fluids at the walls of such corrugated channels or tubes and gives rise to an “apparent” slip in the constricted section of the wall and to a very low value of the drag force on the wall. When and where boundary layers arise, they scale as (1/De) for the stresses and as (De⁄Pe)^(1⁄3) for the concentration. |