This work describes the implementation of an adaptive procedure for viscoelastic flows. Finite element simulations are conducted using a Taylor-Galerkin/pressure correction scheme. The constitutive equations considered are those for an Oldroyd-B and a Phan-ThienlTanner model. The adaptive meshing te
Temporal acceleration techniques for viscoelastic flows
✍ Scribed by Matallah, H. ;Townsend, P. ;Webster, M. F.
- Publisher
- John Wiley and Sons
- Year
- 1998
- Tongue
- English
- Weight
- 251 KB
- Volume
- 14
- Category
- Article
- ISSN
- 1069-8299
No coin nor oath required. For personal study only.
✦ Synopsis
This paper presents a ®nite element study based on a technique associated with time extrapolation to accelerate the convergence rate to the steady state for viscoelastic ¯ows. The approach adopted is a local extrapolation method attributed to Neville. Temporal extrapolation is embedded within a time-marching Taylor±Galerkin/pressure-correction scheme as applied to the solution of model channel ¯ow, 4 : 1 plane contraction ¯ow and ¯ow past a circular cylinder. In particular, consideration is given to obtaining steadystate solutions for an Oldroyd-B model. When extrapolation is performed for stress and velocity or pressure, then stress and velocity overshoot, which consequently leads to divergence. In contrast, a stable numerical scheme emerges when only the stress is extrapolated.
📜 SIMILAR VOLUMES
## Abstract The steady three‐dimensional exterior flow of a viscoelastic non‐Newtonian fluid is approximated by reducing the corresponding nonlinear elliptic–hyperbolic system to a bounded domain. On the truncation surface with a large radius __R__, nonlinear, local second‐order artificial boundary
A multigrid convergence acceleration technique has been developed for solving both the Navier-Stokes and turbulence transport equations. For turbulence closure a low-Reynolds-number q-o turbulence model is employed. To enable convergence, the stiff non-linear turbulent source terms have to be treate
We develop a class of ®lters based upon the numerical solution of high-order elliptic problems in d which allow for independent determination of order and cut-o wave number and which default to classical Fourier-based ®lters in homogeneous domains. However, because they are based on the solution of