Abstract:We establish a compact finite difference algorithm in general curvilinear coordinates with fourth‐order spatial accuracy and second‐order temporal accuracy for the pure streamfunction formulation of the unsteady incompressible Navier‐Stokes equations. This method extends the high‐order pure streamfunction method to general unsteady flow problems with complex geometry and non‐conformal grids. The stability analysis is validated by von‐Neumann linear stability analysis. The accuracy and robustness of the newly proposed method are verified by five numerical experiments. In this paper, a high‐order compact finite difference method in general curvilinear coordinates is proposed for solving unsteady incompressible Navier‐Stokes equations. By constructing the fourth‐order spatial discretization schemes for all partial derivative terms of the pure streamfunction formulation in general curvilinear coordinates, especially for the fourth‐order mixed derivative terms, and applying a Crank‐Nicolson scheme for the second‐order temporal discretization, we extend the unsteady high‐order pure streamfunction algorithm to flow problems with more general non‐conformal grids. Furthermore, the stability of the newly proposed method for the linear model is validated by von‐Neumann linear stability analysis. Five numerical experiments are conducted to verify the accuracy and robustness of the proposed method. The results show that our method not only effectively solves problems with non‐conformal grids, but also allows grid generation and local refinement using commercial software. The solutions are in good agreement with the established numerical and experimental results.

Recent Developments in the Pure Streamfunction Formulation of the Navier-Stokes System

A quasi-Lagrangian finite element method for the Navier-Stokes equations in a time-dependent domain

A Fourth‐order Compact Finite Difference Scheme for the Steady Stream Function–vorticity Formulation of the Navier–Stokes/Boussinesq Equations

A high‐order pure streamfunction method in general curvilinear coordinates for unsteady incompressible viscous flow with complex geometry

A physics-preserving pure streamfunction formulation and high-order compact solver with high-resolution for three-dimensional steady incompressible flows

A high-order upwind compact difference scheme for solving the streamfunction-velocity formulation of the unsteady incompressible Navier–Stokes equations

An efficient compact difference scheme for solving the streamfunction formulation of the incompressible Navier-Stokes equations

A high-order compact scheme for the pure streamfunction (vector potential) formulation of the 3D steady incompressible Navier–Stokes equations

An Upwind Compact Difference Scheme for Solving the Streamfunction-Velocity Formulation of the Unsteady Incompressible Navier-Stokes Equation

On the solutions of Navier-Stokes equations for turbulent channel flows in a particular function class

A high‐order compact scheme for solving the 2D steady incompressible Navier‐Stokes equations in general curvilinear coordinates

Essentially compact schemes for unsteady viscous incompressible flows

A class of exact solutions of the Navier-Stokes equations in three and four dimensions

Convergence of a Finite Difference Scheme for the Navier–Stokes Equations Using Vorticity Boundary Conditions

A Compact Streamfunction-Velocity Scheme for the 2-D Unsteady Incompressible Navier-Stokes Equations in Arbitrary Curvilinear Coordinates

A Higher Order Compact Finite Difference Algorithm for Solving the Incompressible Navier–Stokes Equations

Finite Difference Schemes for Incompressible Flows in Vorticity Formulations

Compact computations based on a stream-function-velocity formulation of two-dimensional steady laminar natural convection in a square cavity.

Improved Crouzeix-Raviart scheme for the Stokes and Navier-Stokes problem

Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations

Navier-Stokes solver using Green's functions I: channel flow and plane Couette flow