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

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,rid and 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 How 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 alogorithms 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 reproduce,,, results from well-studied benchmark problems in viscous fluid Finally, we demonstrate the robustness of our code on flow in a complex domain, (C) 2002 Elsevier Science (USA).