From surface flow visualization on a MIRA notchback reference model using oil flow technique, the topology of the singular points of the skin friction lines are delineated. Separation and reattachment lines at the front screen, at the A-pillar, at the C-pillar and on the rear side of the car including the trunk have been identified. It is worth to mention that two vortices emerging from the top of the trunk coil in the opposite direction as that vortices starting from the C-pillar edge. The positions of the singular points and the separation and reattachment lines and the foci of the vortices provide a sensitive database for validation of CFD-codes

Optimized high-order compact(OHOC) schemes were proposed, which have high spatial order of truncation and resolution to simulate the aeroacoustic problems due to unsteady compressible flows. Generally, numerical schemes are categorized explicit or implicit by time-marching method. In this research, OHOC differences which were developed with explicit time-marching method is used to have implicit formulation and the implicit OHOC differences result in block hepta-diagonal matrix. This paper presents the comparisons between the explicit and implicit OHOC schemes with a second order accuracy of time in the 1-d linear wave convection problem, and between the explicit OHOC scheme of 4th-order accuracy in time and the implicit OHOC scheme of 1st-order accuracy in tine for the 1-d nonlinear wave convection problem. With these comparisons, the characteristics of implicit OHOC scheme are shown in the point of CFL number.

The practical use of the particle image velocimetry(PIV), a whole-field velocity measurement method, requires the use of fast, reliable, computer-based methods for tracking velocity vectors. The full search block matching, the most widely studied and applied technique both in area of PIV and Image Coding and Compression, is computationally costly. Many less expensive alternatives have been proposed mostly in the area of Image Coding and Compression. Among others, TSS, NTSS, HPM are introduced for the past PIV analysis, and found to be successful. But, these algorithms are based on small dynamic range, 7 pixels/frame in maximum displacement. To analyze the images with large displacement, Even and Odd field image separation and a simple version of multi-resolution hierarchical procedures are introduced in this paper. Comparison with other algorithms are summarized. A Results of application to the turbulent backward step flow shows the improvement of new algorithm.

Automatic mesh generation procedures applied to industrial now problems lead to complex mesh topologies where usually no special considerations to mesh resolution are taken. In the present study a fast and flexible solution algorithm in combination with generalized higher order discretization schemes is presented and its application to intake port calculation is demonstrated.

In this paper, the on-going effort and progress for developing an integrated software for automated CFD application is described. As an outcome of the effort devoted so far, a new system, ICFDIT, is developed and introduced in this paper. The new system can be used to solve fluid dynamics problems in a convenient graphical environment, and it includes a pre-processor, a main-processor, and a post-processor. Usage of the system and examples are demonstrated, and some issues for improvement of the system are discussed.

The cross-flow fan performance and its sound noise characteristics are predicted by computational methods. The unsteady incompressible Navier-Stokes equations in moving coordinates are solved by a SMAC method on unstructured triangular meshes, using a sliding mesh technique at the interface between the domain rotating with blades and the rest stationary part. The computationally predicted fan performance was favorably compared with experiment, and some numerical aspects of simulating the cross-flow fan are discussed. With the computed unsteady flow field, aeroacoustic sound noise of the fan is predicted by the Lighthill-Curie equation. The unsteady surface pressure fluctuations on stabilizer enables a prediction of BPF noise of the uniform pitch blade fan quite accurately. The aeroacoustic sound noise characteristics of both uniform and random pitch blade fans are also examined by SPL spectrum analysis.

The pressure-based methods are very popular in CFD because it requires less computer core memory compared to other coupled or density-based solvers. Currently structured-mesh methodology based on pressure-based algorithm is quite mature to apply to the practical problems. The unstructured mesh method needs much more computer memory than the structured-mesh method. However the pressure-based method utilizing the sequential approach does not require very large memory used for unstructured-mesh density-based solvers. The present study has developed the unstructured grid pressure-based method. Cell-centered finite volume method was selected due to robustness for imposing various boundary conditions and easy implementation of higher-order upwind scheme. The predictive capability of present method has validated against several benchmark problems.

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

The purpose of this study is to compare the predictive behaviors of the extended $k-{\varepsilon}$ turbulence model and the standard $k-{\varepsilon}$ turbulence model. Grid dependency is tested with the H-type grid and the O-type grid. Computations have been performed for a circular-to-rectangular transition duct. Numerical results for several sections along the streamwise have been obtained to compare with experimental results. The Reynolds number is 390,000 based on the bulk velocity at the inlet. The computed axial velocity contours, transverse velocity profiles, static pressure contours, peripheral skin friction coefficient, and peripheral wall static pressure distributions have been compared with experimental results. The computed results obtained with the extended $k-{\varepsilon}$ turbulence model show better agreement with experimental results than those obtained with the standard $k-{\varepsilon}$ turbulence model. Comparing to the computed results obtained with the H-type grid and O-type grid, those with H-type grid agree well with experimental results.

Spiral tube heat exchangers can find numerous applications in many engineering field. Flow in spital tubes is interest to engineers due to occurrence of secondary flow which enhances the cross-sectional mixing and the heat transfer rate. In the present study, an incompressible viscous 3-D flow in spiral tubes with rectangular cross-section of various torsion rate and Reynolds number is studied by using a finite volume method. It is shown that the axial velocity profile is affected by the secondary flow motion. Because there is some difference from correlation proposed by Hur et al., a lot of analysis and arrangement of experimental results are needed. This study showed the results of variation of hydrodynamic entry length for torsion and Re numbers.

In this study we calculated the flowfields around the axisymmetric bodies in the rarefied gas regime by using the DSMC. A flat-ended cylinder was selected as a representative axisymmetric body and the gas used for all calculations was nitrogen. With zero angle of attack, an increasingly rapid rise in density and the effect of shock waves near the flat-ended face were examined. And on the cylinder surface velocity slips and boundary layers could be observed in the results. Larsen-Borgnakke model was used for the energy redistribution in inelastic collisions. And by considering all internal energy modes, the distributions of translational, rotational and vibrational temperatures were plotted. The calculations were peformed for various Knudsen numbers, Especially the rotational temperatures calculated for a case were compared with the experimental results. And the simulation results show good agreements with the experimental ones.

A finite-difference method based on conservative supra characteristic method(CSCM) type upwind flux difference splitting has been studied on the bluntness effect on the wall heat transfer rate and wall pressure over blunt-body. The results show that the stagnation heating varies inversely with the square root of the nose radius.

This paper describes aerodynamic heating on a hypersonic vehicle. For this purpose, the 2-D, and 3-D equilibrium code are developed. In order to obtain an accurate solution, AUSMPW+ is used for spatial discretization. Curve fitting data in NASA Reference Publication 1181, 1260 are used to calculate equilibrium properties. To observe aerodynamic heating phenomena, Reynolds number parametric study for diamond airfoil is done, 3-D full Navier-Stokes equation is computed and wall temperature distribution data are obtained. Analyzing these results, we conclude that Reynolds number and secondary flow are important factors in aerodynamic heating.

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

Various discretiation methods of Laplacian operator in the Pressure-Poisson equation are investigated for the solution of incompressible Navier-Stokes equations using the non-staggered grid. Laplacian operators previously proposed by other researchers are applied to a Driven-Cavity problem. The computational results are compared with those of Ghia. The results show the characteristics of the discrete Laplacian operators.

Axisymmetric flows in an optical fiber coating applicator were computed by using structured finite volume solver FLUENT. The numerical computations were performed for various draw speeds from 900m/min to 1,200m/min. The physical phenomena in the coating applicator is the main topic of this paper. The results of the simulations covering streamline, pressure, temperature, and fluid-fluid interface are shown in this paper.

A numerical investigation for the flow, heat and mass transfer characteristics of the grooved evaporating tube with the films flowing down on both the inside and outside tube walls has been carried out. The condensation occurs along the outside wall while the evaporation takes place at the free surface of the inside film. The 3-D transport equations for momentum and energy are solved by using the FVM(Finite Volume Method). The free surface shape is tracked by the moving grid technique satisfying the SCL(Space Conservation Rule). Due to the secondary motion of the fluid, the film thins at the crest, while thickens at the valley. The velocity and temperature fields as well as the amounts of the condensed and evaporated mass have been successfully predicted for various operating conditions and groove shapes.

The present numerical study investigates flow characters and heat transfer enhancement by the viscoelastic-driven secondary flow and buoyancy effect in a 2:1 rectangular duct. Three versions of thermal boundary conditions involving difference combination of heated walls and adiabatic walls are analyzed in this study. The Reiner-Rivlin model is adopted as a viscoelastic fluid model to simulate the secondary flow and temperature-dependent viscosity model is used. Calculated Nusselt numbers are very good agreement with experimental results for reported viscoelastic fluids. It is found that the heat transfer enhancement is mainly caused by the viscoelastic-driven secondary flow and buoyancy-induced secondary flow play a role of promoting this effect.

A numerical optimization procedure is developed to find the inlet velocity profile that yields the most uniform epitaxial layer in a vertical MOCVD reactor. It involves the solution of fully elliptic equations of motion, temperature, and concentration; the finite volume method based on SIMPLE algorithm has been adopted to solve the Navier-Stokes equations. The overall optimization process is highly nonlinear and has been efficiently treated by the sequential linear programming technique that breaks the non-linear problem into a series of linear ones. The optimal profile approximated by a 6th-degree Chebyshev polynomial is very successful in reducing the spatial non-uniformity of the growth rate. The optimization is particularly effective to the high Reynolds number flow. It is also found that a properly constructed inlet velocity profile can suppress the buoyancy driven secondary flow and improve the growth-rate uniformity.

The present study is aimed to investigate flow characteristics of Two dimensional backward-facing step by numerical approach. A convection conservative difference scheme based upon SOLA algorithm is used for the solution of the two-dimensional incompressible Navier-Stokes equations to simulate the laminar, transitional and turbulent flow conditions at which the experimental data can be available for the backward-facing step. The twenty kinds of Reynolds number are used for the calculations. In an effort to demonstrate that the reported solutions are dependent on the mesh refinement, computations are performed on seven different meshes of uniformly increasing refinement. Also to investigate the result of inflow dependence, two kinds of the inflow profile are chosen for the laminar flow. As criterion of benchmarking the result of numerical simulation, reattachment length is used for the selected Reynolds numbers.

A Navier-Stokes code based on a unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number $k-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for a domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. Unsteadiness inside boundary layer is entrained when a unsteady vortex impinge on the blade surface. It shoves that local peak value inside the boundary layer and also local minimum near the edge of boundary layer as it developes along the blade surface. The unsteadiness inside the boundary layer is almost isolated from the free stream unsteadiness and being convected at local boundary layer speed, less than the free stream value.

An Euler/Navier Stokes solver has been developed for the analysis of steady and unsteady flows. The $q-{\omega}$ turbulent model has been incorporated into the solver in strongly coupled manner for stability and robustness. A new Chimera hole cutting algorithm, Cut-paste algorithm, has been devised for automatic Chimera hole cutting. Number of viscous/inviscid numerical computations demonstrate the accuracy and the versatility of the solver.

2-D flow fields are studied by using a shared memory parallel computer with a parallel flow analysis program which uses domain decomposition method and MPI library for data exchange at overlapped interface. Especially, effects of directional domain decomposition on parallel efficiency are studied for 2-D Lid-Driven cavity flow and flow through square cavity. It is known from the present study that domain decomposition along the main flow direction gives better parallel efficiency in 1-D partitioning than along the other direction. 2-D partitioning, however, is less sensitive to flow directions and gives good parallel efficiency for most of the cases considered.

The Newton-Krylov method on the unstructured grid flow solver using the cell-centered spatial discretization oi compressible Euler equations is presented. This flow solver uses the reconstructed primitive variables to get the higher order solutions. To get the quadratic convergence of Newton method with this solver, the careful linearization of face flux is performed with the reconstructed flow variables. The GMRES method is used to solve large sparse matrix and to improve the performance ILU preconditioner is adopted and vectorized with level scheduling algorithm. To get the quadratic convergence with the higher order schemes and to reduce the memory storage. the matrix-free implementation and Barth's matrix-vector method are implemented and compared with the traditional matrix-vector method. The convergence and computing times are compared with each other.

In this paper, incompressible two-dimensional Navier-Stokes equations are numerically solved for the study of steady laminar flow around a body with the wall effect. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. To investigate the wall effect, numerical computations are carried out for the NACA 0012 section at the various blockage ratios. The pressure and skin friction on the foil surface, velocity pronto in its wake and drag coefficient are investigated as functions of the blockage ratio.

The numerical methodology for simulating water entry behaviors of the high-speed bodies has been developed. Since the present method assumed the impact occurs within a very short time interval, the viscous effects do not have enough time to play a significant role in the impact forces, that is, the flow around a water-entry object was assumed as an incompressible potential flow and is solved by the source panel method. The elements fully submerged into the water are routinely teated, but the elements intersected by the effective planar free surface are redefined and reorganized to be amenable to the source panel method. To validate the present code, it has been applied to the ogive model and compared with experimental data. Good agreement has been obtained. The water entry behavior of the bouncing phenomena from the free surface has been also simulated using the impact forces and two degree of freedom dynamic equation. Physically, acceptable results have been obtained.

본 연구에서는 제한된 수로를 운항하는 선체주위의 유동특성에 관한 연구를 수행하고자 하였다. 일반적으로 운하 또는 하천을 운항하는 선박의 경우 제한된 수심의 영향으로 천수효과가 발생하게 된다. 이러한 천수의 효과와 제한된 폭의 영향으로 선수부분에서 선체보다 선행하는 파도가 전파되어 나아가기도 하며, 이로 인하여 선박은 보다 많은 조파저항을 받게 된다. 본 연구에서는 임의의 형상을 갖는 선체가 폭과 수심이 제한된 운하를 임계속도 근처에서 운항하는 경우에 관하여 폭과 수심을 변화시켜가며 수치계산을 수행하여, 제한수로에서의 임의의 선체주위의 유동특성을 관찰하고자 하였다. 수치계산은 MAC(Marker And Cell)법을 기초로 한 유한차분법(Finite Difference Method)을 사용하였으며, 계산에 사용된 격자계는 임의의 형상에 관하여 격자생성이 용이한 직사각형 격자계(Rectangular Grid System)를 사용하였다.

A Numerical simulation on the flow around a Rudder with blowing is performed by Finite Volume Method. The governing equations are three dimensional incompressible Navier-Stokes equation and Continuity equation, Flow field around a finite Rudder including tip vortex is simulated and the change of the lift force by blowing is analyzed.

A two-dimensional airfoil flying over a wavy wall is calculated by solving the unsteady Euler equation. Unsteady Transonic flow over an NACA00012 airfoil in pitching motion has been computed for code validation. Some numerical results for NACA6409 airfoil under different wave number, wave length, fly height are presented. The numerical results show the variation of lift and pitching moment coefficients are increased as wave length decrease.

The discrete element method is a numerical model capable of describing the mechanical behaviour of assemblies of discs and spheres. The method is based on the use of an explicit numerical scheme in which the interaction of the particles is monitored contact by contact and the motion of the particles modelled particle by particle. In this paper, A two-dimensional model for computing contacts and motions of granular particles of unform, inelasticity is presented. And, code is developed. The primary aim of this paper is to approv computational result of continuum alaysis which is on processing. The end of this paper, that code is tested with several examples.

A method is developed to include the effect of volume expansion in the description of the flame dynamics using G-equation. Line volume-source is used to represent the effect of the exothermic process of combustion with source strength determined by the density difference between the burned and the unburned region. Volume expansion adjusts the flow field to accommodate the increased volume flow rate crossing the flame front. Test result predicted the measured velocity field qualitatively. The method was applied to study the interaction of vortex and premixed flame. Increased volume expansion did not change the initial growth rate of flame area. However, the residence time and flame surface area increased with higher expansion ratios.

Propellants pressurized and fed into the combustion chamber undergoes the mechanical, chemical combustion processes. Along with their distinctive physical characteristics, propellant combustion is typically divided into the processes; injection, atomization, mixing, vaporization and chemical reaction. These processes assumed to happen in a serial manner are strongly coupled, thereby involves formidable physical complexities. In this study a numerical experiment is attempted to simulate the burning sprays due to OFO, FOF triplet / FOOF split doublet injectors. Based on Eulerian-Lagrangian frame, Navier-Stokes equation system for compressible flows is preconditioned with low Reynolds number $k-{\varepsilon}$ turbulent model and time-integrated by LU-SGS, and the sprays are described by DSF model with the characteristics initialized by experimentally determined spray characteristics. Simplified single global reaction model approximates heptane-air reaction. It was observed that FOOF split doublet injector shows better atmization with shortest residence and the FOF triplet injector produces better combustion performance.