In this work, we present a numerical study on the development length (the length from the channel inlet required for the velocity to reach $99\%$ of its fully-developed value) of a pressure-driven viscoelastic fluid flow (between parallel plates) modelled by the generalised Phan–Thien and Tanner (gPTT) constitutive equation. The governing equations are solved using the finite-difference method, and, a thorough analysis on the effect of the model parameters $\alpha$ and $\beta$ is presented. The numerical results showed that in the creeping flow limit ($Re=0$), the development length for the velocity exhibits a non-monotonic behaviour. The development length increases with $Wi$. For low values of $Wi$, the highest value of the development length is obtained for $\alpha =\beta =0.5$; for high values of $Wi$, the highest value of the development length is obtained for $\alpha =\beta =1.5$. This work also considers the influence of the elasticity number.

In this work, we present the implementation and verification of HiGTree-HiGFlow solver (see for numerical simulation of the KBKZ integral constitutive equation. The numerical method proposed herein is a finite difference technique using tree-based grids. The advantage of using hierarchical grids is that they allow us to achieve great accuracy in local mesh refinements. A moving least squares (MLS) interpolation technique is used to adapt the discretization stencil near the interfaces between grid elements of different sizes. The momentum and mass conservation equations are solved by an implicit method and the Chorin projection method is used for decoupling the velocity and pressure. The Finger tensor is calculated using the deformation fields method and a three-node quadrature formula is used to derive an expression for the integral tensor. The results of velocity and stress fields in channel and contraction-flow problems obtained in our simulations show good agreement with numerical and experimental results found in the literature.

This work introduces a novel numerical method designed to address three-dimensional unsteady free surface flows incorporating integral viscoelastic constitutive equations, specifically the K–BKZ–PSM (Kaye–Bernstein, Kearsley, Zapas–Papanastasiou, Scriven, Macosko) model. The new proposed methodology employs a second-order finite difference approach along with the deformation fields method to solve the integral constitutive equation and the marker particle method (known as marker-and-cell) to accurately capture the evolution of the fluid’s free surface. The newly developed numerical method has proven its effectiveness in handling complex fluid flow scenarios, including confined flows and extrudate swell simulations of Boger fluids. Furthermore, a new semi-analytical solution for velocity and stress fields is derived, considering fully developed flows of a K–BKZ–PSM fluid in a pipe.

This work introduces a novel numerical method designed to address three-dimensional unsteady free surface flows incorporating integral viscoelastic constitutive equations, specifically, the K-BKZ-PSM (Kaye–Bernstein, Kearsley, Zapas - Papanastasiou, Scriven, Macosko) model. The new proposed methodology employs a second-order finite difference approach along with the deformation fields method to solve the integral constitutive equation, and the marker-particle method (known as \emph{Marker-and-Cell}) to accurately capture the evolution of the fluid's free surface. The newly developed numerical method has proven its effectiveness in handling complex fluid flow scenarios, including confined flows and extrudate swell simulations of Boger fluids. Furthermore, a new semi-analytical solution for velocity and stress fields is derived, considering fully developed flows of a K-BKZ-PSM fluid in a pipe.

In this paper we introduce a technique to deal with implicit differential equations exhibiting singularities. Our approach is a geometrical one, we use the concept of contact structure on a manifold associated with the differential equation. In this setting we prove an existence and uniqueness theorem. We also show how it relates to known geometric results for this kind of equation. We also indicate how the method can be implemented by using continuation methods techniques and the BDF (Backward Differentiation Formula).