Numerical Analysis of the Blow-Up of One-Dimensional Polymer Fluid Flow with a Front

Abstract

One-dimensional flows of an incompressible viscoelastic polymer fluid that are qualitatively similar to the solutions of Burgers’ equation are described on the basis of mesoscopic approach for the first time. The corresponding initial boundary-value problem is posed for the system of quasilinear differential equations. The numerical algorithm for solving it is designed and verified. The algorithm uses the explicit fifth-order scheme to approximate unknown functions with respect to time variable and the rational barycentric interpolations with respect to space variable. A method for localization of singular points of the solution in the complex plain and for adaptation of the spatial grid to them is implemented using the Chebyshev-Padé approximations. Two regimes of evolution of the solution to the problem are discovered and characterized while using the algorithm: regime 1—a smooth solution exists in a sufficiently large time interval (the singular point moves parallel to the real axis in the complex plane); regime 2—the smooth solution blows up at the beginning of evolution (the singular point reaches the segment of the real axis where the problem is posed). We study the influence of the rheological parameters of fluid on the realizability of these regimes and on the length of time interval where the smooth solution exists. The obtained results are important for the analysis of laminar-turbulent transitions in viscoelastic polymer continua.