Communications in Applied Mathematics and Computational Science
- Commun. Appl. Math. Comput. Sci.
- Volume 4, Number 1 (2009), 57-83.
A higher-order upwind method for viscoelastic flow
We present a conservative finite difference method designed to capture elastic wave propagation in viscoelastic fluids in two dimensions. We model the incompressible Navier–Stokes equations with an extra viscoelastic stress described by the Oldroyd-B constitutive equations. The equations are cast into a hybrid conservation form which is amenable to the use of a second-order Godunov method for the hyperbolic part of the equations, including a new exact Riemann solver. A numerical stress splitting technique provides a well-posed discretization for the entire range of Newtonian and elastic fluids. Incompressibility is enforced through a projection method and a partitioning of variables that suppresses compressive waves. Irregular geometry is treated with an embedded boundary/volume-of-fluid approach. The method is stable for time steps governed by the advective Courant–Friedrichs–Lewy (CFL) condition. We present second-order convergence results in for a range of Oldroyd-B fluids.
Commun. Appl. Math. Comput. Sci., Volume 4, Number 1 (2009), 57-83.
Received: 7 August 2008
Revised: 22 May 2009
Accepted: 25 May 2009
First available in Project Euclid: 20 December 2017
Permanent link to this document
Digital Object Identifier
Mathematical Reviews number (MathSciNet)
Zentralblatt MATH identifier
Nonaka, Andrew; Trebotich, David; Miller, Gregory; Graves, Daniel; Colella, Phillip. A higher-order upwind method for viscoelastic flow. Commun. Appl. Math. Comput. Sci. 4 (2009), no. 1, 57--83. doi:10.2140/camcos.2009.4.57. https://projecteuclid.org/euclid.camcos/1513798577