| A pressure-based computer program for a general Navier-Stokes equations has been developed. Body-fitted coordinate system is employed to handle flows with complex geometry. Non-staggered grid is used while the pressure oscillation is eliminated by a special pressure interpolation scheme. The hybrid algorithm is adopted to discretize the equations and the finite-difference equations are solved by TDMA, while the whole solution is obtained through an under-relaxed iterative process. The pressure field is evaluated using the compressible from of the SIMPLE algorithm.; To test the accuracy and efficiency of the computer program, problems of incompressible and compressible flows are calculated. As examples of inviscid compressible flow problems, flows over a bump with 10% and 4% thickness are computed with the incoming Mach numbers of M{dollar}sb{lcub}infty{rcub}{dollar} = 0.5 (subsonic flow), M{dollar}sb{lcub}infty{rcub}{dollar} = 0.675 (transonic flow and M{dollar}sb{lcub}infty{rcub}{dollar} = 1.65 (supersonic flow). One laminar subsonic flow over a bump with 5% thickness at M{dollar}sb{lcub}infty{rcub}{dollar} = 0.5 is also calculated with the consideration of the full energy equation. With the help of the k-{dollar}epsilon{dollar} model incorporating the wall function, the computations of two turbulent incompressible flows are carried out. One is the flow past a flat plate and the other over a flame holder. As an application to the three-dimensional flow, a laminar flow in a driven cubic cavity is calculated. All the numerical results obtained here are compared with experimental data or other numerical results available in the literature. |
