A convergent and universally bounded interpolation scheme for the treatment of advection

• Published in: International journal for numerical methods in fluids 2003-01, Vol.41 (1), p.47-75
• Main Authors: Alves, M. A., Oliveira, P. J., Pinho, F. T.
• Format: Article
• Language: English
• Subjects: boundedness; Computational methods in fluid dynamics; CUBISTA; Exact sciences and technology; Flows in ducts, channels, nozzles, and conduits; Fluid dynamics; Fundamental areas of phenomenology (including applications); high resolution scheme (HRS); iterative convergence; Non-newtonian fluid flows; Physics; TVD property; viscoelastic fluids
• ISSN: 0271-2091; 1097-0363
• DOI: 10.1002/fld.428

A high resolution scheme with improved iterative convergence properties was devised by incorporating total‐variation diminishing constraints, appropriate for unsteady problems, into an implicit time‐marching method used for steady flow problems. The new scheme, referred to as Convergent and Universally Bounded Interpolation Scheme for the Treatment of Advection (CUBISTA), has similar accuracy to the well‐known SMART scheme, both being formally third‐order accurate on uniform meshes for smooth flows. Three demonstration problems are considered: (1) advection of three scalar profiles, a step, a sine‐squared, and a semi‐ellipse; (2) Newtonian flow over a backward‐facing step; and (3) viscoelastic flow through a planar contraction and around a cylinder. For the case of the viscoelastic flows, in which the high resolution schemes are also used to represent the advective terms in the constitutive equation, it is shown that only the new scheme is able to provide a converged solution to the prescribed tolerance. Copyright © 2003 John Wiley & Sons, Ltd.