• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.153-158
• /
• 2000

For the analysis of compressible flow with real gas effect, characteristic form of Roe's Riemann solver was derived again using real gas equation of state and variable specific heal and it was extended to multi component reactive system. From this study, it is known that some correction should be made for the use of existing numerical algorithm. 1) Sonic speed and characteristic variable should be corrected with real gas effect. 2) Roe's average was applicable only with the assumption of constant properties. 3) Artificial damping term and characteristic variables should be corrected but their influence may not be significant.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1996.11a
• /
• pp.107-116
• /
• 1996

Numerical simulations have been conducted to characterize unsteady flow structures in an axisymmetric supersonic inlet diffuser with sinusoidal pressure oscillations at the diffuser exit. The formulation is based on the unsteady Navier-Stokes equations and turbulence closure is achieved using a two-layer model with a too-Reynolds-number scheme for the near-wall treatment. The governing equations are formulated in an integral form, and are discretized by the four-stage Runge-Kutta scheme for temporal terms and the Harten-Yee upwind TVD scheme for convective terms. Results indicated that the inlet shock characteristics are significantly modified by acoustic oscillations originating from the combustor. The characteristics of shock/boundarv-layer interactions (such as the size of separation bubble, terminal shock shape, and vorticity intensity) are also greatly iufluenced by the shock oscillation due to acoustic waves.

• Proceedings of the KSR Conference
• /
• 2004.06a
• /
• pp.685-689
• /
• 2004

Numerical simulation is conducted to clarify the heat transfer and fluid flow characteristics of HVAC(Heating, Ventilating and Air-Conditioning) for double' deck train. The HVAC system is installed under the roof of carbody. In the lay-out of HVAC system, air duct must be installed to supply air to 1st and 2nd floor respectively. The standard k-$\epsilon$ and LES models for turbulence and SIMPLE algorithm for pressure equation hased on finite volume method are used to solve the physic a] HVAC model. To assure convergence, QUICK scheme for momentum equation and the 2nd order upwind scheme for turbulent equations arc used. From the results of simulation, the temperature and velocity magnitude are also distributed uniformly in the interior of double-deck passenger car.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.42-47
• /
• 1998

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.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.98-111
• /
• 1996

The compressible laminar and turbulent viscous flows on a slender body in supersonic speed as well as subsonic speed have been numerically simulated at high angle of attack. The steady and time-accurate compressible thin-layer Navier-Stokes code based on an implicit upwind-biased LU-SGS algorithm has been developed and specifically applied at angles of attack of 20, 30 and 40 dog, respectively. The modified eddy-viscosity turbulence model suggested by Degani and Schiff was used to simulate the case of turbulent flow. Any geometric asymmetry and numerical perturbation have not been intentionally or artificially imposed in the process of computation. The purely numerical results for laminar and turbulent cases, however, show clear asymmetric formation of vortices which were observed experimentally. Contrary to the subsonic results, the supersonic case shows the symmetric formation of vortices as indicated by the earlier experiments.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.13-20
• /
• 1997

A fundamental study for the numerical scheme to simulate unsteady nonlinear waves by solving Euler equations is presented. First a conservation form and a non-conservation form of the Euler equations with a free surface fitted coordinate system are compared. Next, a time splitting fractional step method and an alternating direction implicit(ADI) method for the time integration are compared. For the comparative study, flow calculations around a bottom bump in a channel and a NACA 0012 hydrofoil in a flume are performed. The results show that the ADI method with a third order upwind differencing scheme is very efficient in reducing the computing time with keeping the accuracy, And, there is no distinct difference between two expression forms except that the non-conservative form shows faster wave propagating velocity than the conservation form. Some results are compared with experiments and show good agreement.

• Journal of computational fluids engineering
• /
• v.5 no.2
• /
• pp.20-27
• /
• 2000

An unstructured implicit Euler solver is parallelized on a Cray T3E. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by the Gauss-Seidel implicit scheme. Domain decomposition is accomplished by using the k-way n-partitioning method developed by Karypis. In order to analyze the parallel performance of the solver, flows over a 2-D NACA 0012 airfoil and 3-D F-5 wing were investigated.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.41-48
• /
• 2013

Ice accretion on the solid surface is an importance factor in assessing the performance of aircraft and wind turbine blade. Changes in the external shape due to ice accretion can greatly deteriorate the aerodynamic performance. In this study, a three-dimensional upwind-type second-order positivity-preserving finite volume CFD scheme based on the unstructured mesh topology is developed to simulate two-phase flow in atmospheric icing condition. The code is then validated by comparing with NASA IRT experimental data on the sphere. The present results of the collection efficiency are found to be in close agreement with experimental data and show improvement near the stagnation region.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.29-35
• /
• 1999

A two-dimensional Euler code based on flux vector splitting scheme has been developed to simulate the behavior of supersonic shock wave over the cylinder. AF+ADI scheme was used for time integration. The sliding multiblock technique was implemented to handle the relative motion of the moving cylinder and the stationary tunnel. The code is validated with a problem of subsonic flow around a Naca-0012 airfoil. The Computation results show complex phenomena of the propagation of shock waves and the reflection as expansion wave at tunnel exit.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.21-24
• /
• 2005

Numerical analysis were made for the unsteady flow fields of the rotor system of a Tilt-Rotor aircraft in cruise mode. The Reynolds-averaged thin-layer Wavier-Stokes equations were discretized by Roe's upwind differencing scheme and integrated in time by the LU-SGS algorithm. The computational domain of the rotor system was constructed by seven multi-block Chimera grids. Comparison of pressure coefficient on the surface of the main wing and blades were made for 3cases of advance ratio(0.325, 0.350, 0.375) and thrust and power coefficients for the rotor were compared with experimental data.