An implicit finite-difference procedure for unsteady three-dimensional flow capable of handling arbitrary geometry through the use of general coordinate transformations is described. Viscous effects are optionally incorporated with a "thin-layer" approximation of the Navier-Stokes equations. An implicit approximate factorization technique is employed so that the small grid sizes required for spatial accuracy and viscous resolution do not impose stringent stability limitations. Results obtained from the program include transonic inviscid or viscous solutions about simple body configurations. Comparisons with existing theories and experiments are made. Numerical accuracy and the effect of three-dimensional coordinate singularities are also discussed.

Implicit Finite-Difference Simulations of Three-Dimensional Compressible Flow

Finite-difference procedures are used to solve either the Euler equations or the "thin-layer" Navier-Stokes equations subject to arbitrary boundary conditions. An automatic grid generation program is employed, and because an implicit finite-difference algorithm for the flow equations is used, time steps are not severely limited when grid points are finely distributed. Computational efficiency and compatibility to vectorized computer processors is maintained by use of approximate factorization techniques. Computed results for both inviscid and viscous flow about airfoils are described and compared to viscous known solutions.

/pdf/implicit-finite-difference-simulation-of-flow-about-2sbe4xojqu.pdf

Implicit Finite-Difference Simulation of Flow about Arbitrary Two-Dimensional Geometries

Spectral methods for problems in complex geometries

Boundary-fitted coordinate systems for numerical solution of partial differential equations—A review

The transition of the cylinder wake is investigated experimentally in a water channel and is computed numerically using a finite‐difference scheme. Four physically different instabilities are observed: a local ‘‘vortex‐adhesion mode,’’ and three near‐wake instabilities, which are associated with three different spanwise wavelengths of approximately 1, 2, and 4 diam. All four instability processes can originate in a narrow Reynolds‐number interval between 160 and 230, and may give rise to different transition scenarios. Thus, Williamson’s [Phys. Fluids 31, 3165 (1988)] experimental observation of a hard transition is for the first time numerically reproduced, and is found to be induced by the vortex‐adhesion mode. Without vortex adhesion, a soft onset of three‐dimensionality is numerically and experimentally obtained. A control‐wire technique is proposed, which suppresses transition up to a Reynolds number of 230.

http://berndnoack.com/publications/Zhang1995pf.PDF

On the transition of the cylinder wake

A method is presented for automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multi-connected two-dimensional region containing any number of arbitrarily shaped bodies. No restrictions are placed on the shape of the boundaries, which may even be time-dependent, and the approach is not restricted in principle to two dimensions. With this procedure the numerical solution of a partial differential system may be done on a fixed rectangular field with a square mesh with no interpolation required regardless of the shape of the physical boundaries, regardless of the spacing of the curvilinear coordinate lines in the physical field, and regardless of the movement of the coordinate system in the physical plane. A number of examples of coordinate systems and application thereof to the solution of partial differential equations are given. The FORTRAN computer program and instructions for use are included.

/pdf/boundary-fitted-curvilinear-coordinate-systems-for-solution-14g1ykpk03.pdf

Boundary-fitted curvilinear coordinate systems for solution of partial differential equations on fields containing any number of arbitrary two-dimensional bodies

A procedure for numerical solution of the time-dependent, two-dimensional incompressible Navier-Stokes equations that can treat the unsteady laminar flow about bodies of arbitrary shape, such as two-dimensional airfoils, multiple airfoils, and submerged hydrofoils, as naturally as it can deal with the flow about simple bodies. The solution is based on a method of automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multiconnected region containing any number of arbitrarily shaped bodies. The curvilinear coordinates are generated as the solution of two elliptical partial differential equations with Dirichlet boundary conditions, one coordinate being specified to be constant on each of the boundaries, and a distribution of the other being specified along the boundaries. The solution compares excellently with the Blasius boundary layer solution for the flow past a semiinfinite flat plate.

Use of numerically generated body-fitted coordinate systems for solution of the Navier-Stokes equations

A method for numerical solution of the Navier-Stokes equations for the flow about arbitrary airfoils or other bodies is presented. This method utilizes a numerically generated curvilinear coordinate system having a coordinate line coincident with the body contour. Streamlines, velocity profiles, and pressure and force coefficients for several airfoils and an arbitrary rock are given. Potential flow solutions are also presented. The procedure capable of treating multiple-element airfoils, and potential flow results are presented.

/pdf/numerical-solution-of-the-navier-stokes-equations-for-5wkhu5851a.pdf

Numerical solution of the Navier-Stokes equations for arbitrary two-dimensional airfoils

Solutions of the navier-stokes equations in various flow regimes on fields containing any number of arbitrary bodies using boundary-fitted coordinate systems

C. W. Mastin

Papers

Boundary-fitted curvilinear coordinate systems for solution of partial differential equations on fields containing any number of arbitrary two-dimensional bodies

Use of numerically generated body-fitted coordinate systems for solution of the Navier-Stokes equations

Numerical solution of the Navier-Stokes equations for arbitrary two-dimensional airfoils

Solutions of the navier-stokes equations in various flow regimes on fields containing any number of arbitrary bodies using boundary-fitted coordinate systems