• Three constitutive formulations of the Leonov model are numerically investigated. • It is verified that keeping positive-definiteness of the conformation tensor is crucial for computational evolution while appropriate interpolation of exponentially convecting domain is essential for accuracy. • Swirling asymmetric corner vortices are described in 3D contraction flow modeling. For the purpose of defeating the high Deborah (or Weissenberg) number problems in viscoelastic flow modeling, the matrix-logarithmic formulation of the constitutive equations has been suggested and proven probably the most effective. This transformation of the conformation tensor endows proper resolution of the boundary layer with exponential convection as well as conservation of its positive definiteness. Herein we make an attempt to evaluate numerical consequences of these two features independently for the Leonov viscoelastic model mainly in the contraction flow. To this end, we from the viewpoint of computational stability and accuracy investigate three formulations written with variables such as the original conformation tensor c (c form), its matrix-logarithm (h form) and additive decomposition (b form). 2D planar 4:1 contraction flow in the limit of vanishing inertia in addition to its 1D toy problem is considered in the framework of finite element method. The h and b forms are designed to preserve the positive definiteness of c , and furthermore the h form grants efficient interpolation of exponentially convecting variables. The 1D analogy expresses the characteristics qualitatively equivalent to the original problem in a simpler fashion. The c form modeling of 2D flow exhibits typical nonevolutionary behavior when the Deborah number ( De ) exceeds a certain limit, and this degradation gets worse with finer mesh. While the h and b forms have always been computationally stable with proper description of steady behavior or elastic instability otherwise, the h form manifests faster mesh convergence. As a complementary study, the h form calculation of full 3D flow along an axisymmetric contraction channel has been performed, which also attests well-posedness without upper bound in De . Its steady result is corroborated with the one under axisymmetric flow approximation. Loss of convergence in the axisymmetric flow regime at high De is verified to result from the occurrence of 3D instability, which is then well portrayed in 3D modeling. In this instability, swirling of asymmetric structure for the corner vortices about the geometrical symmetry axis plausibly mimics experimental observation. Consequently, maintaining the positive definiteness in terms of h or b form is crucial for numerically stable evolution, and the h form interpolation delineating appropriately the exponential convection assures the fastest mesh convergence, hence enhancing the accuracy.
Read full abstract