In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids

No abstract

The present work is concerned with a semi-implicit modification of the GENSMAC method for solving the two-dimensional time-dependent incompressible Navier-Stokes equations in primitive variables formulation with a free surface. A projection method is employed to uncouple the velocity components and pressure, thus allowing the solution of each variable separately (a segregated approach). The viscous terms are treated by the implicit backward method in time and a centered second order method in space, and the nonlinear convection terms are explicitly approximated by the high order upwind variable-order nonoscillatory scheme method in space. The boundary conditions at the free surface couple the otherwise segregated velocity and pressure fields. The present work proposes a method that allows the segregated solution of free surface flow problems to be computed by semi-implicit schemes that preserve the stability conditions of the related coupled semi-implicit scheme. The numerical method is applied to both the simulation of free surface and to confined flows. The numerical results demonstrate that the present technique eliminates the parabolic stability restriction required by the original explicit GENSMAC method, and also found in segregated semi-implicit methods with time-lagged boundary conditions. For low Reynolds number flows, the method is robust and very efficient when compared to the original GENSMAC method.

This work presents a numerical application of a generic conformation tensor transformation for simulating highly elastic flows of non-Newtonian fluids typically observed in computational rheology. In the Kernel-conformation framework [14], the conformation tensor constitutive law for a viscoelastic fluid is transformed introducing a generic tensor transformation function. The numerical stability of the application of the Kernel-conformation for highly elastic flows is ultimately related with the specific kernel function used in the matrix transformation, but also to the existence of singularities introduced either by flow geometry or by the characteristics of the constitutive equation. In this work, we implement this methodology in a free-surface Marker-And-Cell discretization methodology implemented in a finite differences method. The main contributions of this work are two fold: on one hand, we demonstrate the accuracy of this Kernel-conformation formulation using a finite differences method and free surfaces; on the other hand, we assess the numerical efficiency of specific kernel functions at high-Weissenberg number flows. The numerical study considers different viscoelastic fluid flow problems, including the Poiseuille flow in a channel, the lid-driven cavity flow and the die-swell free surface flow. The numerical results demonstrate the adequacy of this methodology for high Weissenberg number flows using the Oldroyd-B model.

A practical high resolution upwind differencing scheme for the numerical solution of convection-dominated transport problems is presented. The scheme is based on TVD and CBC stability criteria and is implemented in the context of the finite difference methodology. The performance of the scheme is investigated by solving the 1D/2D scalar advection equations, 1D inviscid Burgers’ equation, 1D scalar convection–diffusion equation, 1D/2D compressible Euler’s equations, and 2D incompressible Navier–Stokes equations. The numerical results displayed good agreement with other existing numerical and experimental data

This article deals with the study of the development and application of the high-order upwind ADBQUICKEST scheme, an adaptative bounded version of the QUICKEST for unsteady problems (Commun. Numer. Meth. Engng 2007; 23:419-445), employing both linear and nonlinear convection term discretization. This scheme is applicable to a wide range of computational fluid dynamics problems, where transport phenomena are of special importance. In particular, the performance of the scheme is assessed through an extensive numerical simulation study of advection-diffusion problems. The scheme, implemented in the context of finite difference methodology, combines a good approximation of shocks (or discontinuities) with a good approximation of the smooth parts of the solutions. In order to assess the performance of the scheme, seven problems are solved, namely (a) advection of scalars; (b) non-linear viscous Burgers equation; (c) Euler equations of gas dynamics; (d) Newtonian flow in a channel; (e) axisymmetric Newtonian jet flow; (f) axisymmetric non-Newtonian (generalized Newtonian) flow in a pipe; and (g) collapse of a fluid column. The numerical experiments clearly show that the scheme provides more consistent solutions than those found in the literature. From the study, the flexibility and robustness of the ADBQUICKEST scheme is confirmed by demonstrating its capability to solve a variety of linear and nonlinear problems with and without discontinuous solutions

No abstract

SUMMARYThis paper is concerned with the numerical solutions of time dependent two-dimensional incompressible flows. By using the primitive variables of velocity and pressure, the Navier-Stokes and mass conservation equations are solved by a semi-implicit finite difference projection method. A new bounded higher order upwind convection scheme is employed to deal with the non-linear (advective) terms. The procedure is an adaptation of the GENSMAC (J. Comput. Phys. 1994; 110:171-186) methodology for calculating confined and free surface fluid flows at both low and high Reynolds numbers. The calculations were performed by using the 2D version of the Freeflow simulation system (J. Comp. Visual. Science 2000; 2:199-210). In order to demonstrate the capabilities of the numerical method, various test cases are presented. These are the fully developed flow in a channel, the flow over a backward facing step, the die-swell problem, the broken dam flow, and an impinging jet onto a flat plate. The numerical results compare favourably with the experimental data and the analytical solutions.