Pressure-based residual smoothing operators for multistage pseudocompressibility algorithms

F. Sotiropoulos, G. Constantinescu

Research output: Contribution to journalArticlepeer-review

15 Scopus citations

Abstract

Implicit residual smoothing operators for accelerating the convergence of explicit, multistage, artificial compressibility algorithms are developed using ideas from pressure-based methods. The velocity derivatives in the continuity equation and the pressure gradient terms in the momentum equations are discretized in time implicitly. The discrete system of equations is linearized in time producing a block implicit operator which is approximately factorized and diagonalized via a similarity transformation. The so-derived diagonal operator depends only on the metrics of the geometric transformation and can, thus, be implemented in an efficient and straightforward manner. It is combined with the standard implicit residual smoothing operator and incorporated in a four-stage Runge-Kutta algorithm also enhanced with local time-stepping and multigrid acceleration. Linear stability analysis for the three-dimensional Navier-Stokes equations and calculations for laminar flows through curved square ducts and pipes demonstrate the damping properties and efficiency of the proposed approach particularly on large-aspect ratio, highly skewed meshes.

Original languageEnglish (US)
Pages (from-to)129-145
Number of pages17
JournalJournal of Computational Physics
Volume133
Issue number1
DOIs
StatePublished - May 1 1997
Externally publishedYes

Bibliographical note

Funding Information:
This work was supported by a grant from the Electric Power Research Institute (EPRI), monitored by Mr. J. L. Tsou, and carried out under the general guidance and supervision of Professor V. C. Patel. Some calculations during the development stage of this work were carried out on the Cray C90 at the San Diego Supercomputing Center (SDSC).

Fingerprint

Dive into the research topics of 'Pressure-based residual smoothing operators for multistage pseudocompressibility algorithms'. Together they form a unique fingerprint.

Cite this