A Cartesian grid method for solving the two-dimensional streamfunction-vorticity equations in irregular regions

Research output: Contribution to journalArticlepeer-review

287 Scopus citations

Abstract

We describe a method for solving the two-dimensional Navier-Stokes equations in irregular physical domains. Our method is based on an underlying uniform Cartesian grid and second-order finite-difference/finite-volume discretizations of the streamfunction-vorticity equations. Geometry representing stationary solid obstacles in the flow domain is embedded in the Cartesian grid and special discretizations near the embedded boundary ensure the accuracy of the solution in the cut cells. Along the embedded boundary, we determine a distribution of vorticity sources needed to impose the no-slip flow conditions. This distribution appears as a right-hand-side term in the discretized fluid equations, and so we can use fast solvers to solve the linear systems that arise. To handle the advective terms, we use the high-resolution algorithms in CLAWPACK. We show that our Stokes solver is second-order accurate for steady state solutions and that our full Navier-Stokes solver is between first- and second-order accurate and reproduces results from well-studied benchmark problems in viscous fluid flow. Finally, we demonstrate the robustness of our code on flow in a complex domain.

Original languageEnglish
Pages (from-to)231-275
Number of pages45
JournalJournal of Computational Physics
Volume176
Issue number2
DOIs
StatePublished - 1 Mar 2002

Keywords

  • Cartesian grid
  • Computational fluid dynamics
  • Embedded boundary
  • Finite-difference
  • Incompressible flow
  • Streamfunction-vorticity finite-volume

Fingerprint

Dive into the research topics of 'A Cartesian grid method for solving the two-dimensional streamfunction-vorticity equations in irregular regions'. Together they form a unique fingerprint.

Cite this