• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.3
• /
• pp.1015-1027
• /
• 1996

In this paper, a numerical procedure for the calculation of the incompressible Navier-Stokes equations in a generalized nonorthogonal body fitted coordinate system is proposed and is validated through three test problems. Present numerical procedure derives the pressure equation by using the pressure substitution method on the regular grid system, and discretized momentum equations are based on the covariant velocity components. Cavity flow, backward facing step flow, and two dimensional channel flow with a sinusoidal wavy wall are chosen as three test problems. Numerical solutions obtained by present procedure shows a good agreement with previous numerical and/or experimental results. Convergence rate is also satisfactory.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2009.05a
• /
• pp.133-137
• /
• 2009

Numerical simulations of flows in an axisymmetric supersonic inlet with bleed regions were performed. For the simulations, the existing code which solves the RANS(Reynolds Averaged Navier-Stokes) equations and 2-equation turbulence model equations was transformed to axisymmetric form and bleed boundary condition was applied to the code. In this paper, the modified code was validated by comparing the results against an experimental data and other computational results for flow on a bump and over an oblique shock with bleed region. Using the code, numerical simulations were performed for the flow in the inlet with multiple bleed regions.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.31-35
• /
• 2006

The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.66-66
• /
• 2003

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.8-13
• /
• 2008

We present a numerical simulation technique and some preliminary results of the impact and spreading of a droplet containing particles on the solid substrate in 2D. We used the 2nd-order Adams-Bashforth / Crank-Nicholson method to solve the Navier-Stokes equation and employed the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension. The impact velocity has been generated by the instantaneous gravity. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles and the discontinuous Galerkin method has been used for the stabilization of the interface advection equation. We investigated the droplet spreading by the inertial force and discussed effects of the presence of particles on the spreading behavior using an example problem. We observed reduced oscillation and spread for the particulate droplet.

• 한국가시화정보학회:학술대회논문집
• /
• 2005.12a
• /
• pp.103-106
• /
• 2005

Lattice Boltzmann method is a new numerical method of investigating the fluid flow which have been solved by Navier-Stokes equation recently. It is known that making the single and parallel algorithms of the Lattice Boltzmann equation is easier than those of Navier-Stokes equations. Also, we can simulate the two phase flow using either the 'Interaction Potential model ' introduced by Shan and Chen. In this paper, we first compared the 3D cavity results of Lattice Boltzmann method with other numerical results for validation and showed the 3D phase transition and its simple application by using the ' Interaction Potential model'

• Journal of computational fluids engineering
• /
• v.7 no.2
• /
• pp.25-32
• /
• 2002

The velocity and temperature profiles in the cabin of the automobile affect greatly the comfortableness of passengers. In this paper, the three dimensional flow and temperature analysis in the cabin of real automobile have been peformed. The three dimensional Navier-Stokes equation solver was validated by comparing with the other numerical data of a 1/5 scale model. The temperature field of cavity was also analyzed for the validation of energy equation solver. After the code validation, the numerical analysis of real field of flow and temperature of an automobile was peformed and the present result provides the insight of flow and temperature field of the inside of cabin.

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

Aerodynamic sensitivity analysis is performed for the Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method and a direct differentiation method respectively. Like the mean flow equations, the turbulence model equations are also hand-differentiated to accurately calculate the sensitivity derivatives of flow quantities with respect to design variables in turbulent viscous flows. Both the direct differentiation code and the adjoint variable code adopt the same time integration scheme with the flow solver to efficiently solve the differentiated equations. The sensitivity codes are then compared with the flow solver in terms of solution accuracy, computing time and computer memory requirements. The sensitivity derivatives obtained from the sensitivity codes with different turbulence models are compared with each other. Using two-equation turbulence models, it is observed that a usual assumption of constant turbulent eddy viscosity in adjoint methods may lead to seriously inaccurate results in highly turbulent flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.75-80
• /
• 2005

The present work was undertaken to numerically analyze the defrosting phenomena of windshield glass. In order to analysis the phase change from frost to water on windshield glass by discharging hot air from a defroster nozzle, the flow and the temperature field of the cabin interior, the heat transfer through the windshield glass, and the phase change of frost should be solve simultaneously. In the present work, the flow field was obtained by solving 3-D incompressible Navier-Stokes equations, and the temperature field was computed from the incompressible energy equation. The phase change process was solved by the enthalpy method. For the code validation, the temperature and the phase change of the driven cavity were calculated. The calculation showed a good agreement with other numerical results. Then, the present code was applied to the defrosting problem of a real automobile, and a good agreement with the experimental data was also obtained.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.189-190
• /
• 2003

A New kinetic scheme for the compressible Navier-Stokes equations is developed. While the conventional approach, such as KFVS scheme, employs the splitting algorithm and computes the numerical flux on the basis of the collisionless equation, the present approach employs the splitting algorithm in the evaluation of the numerical flux, where the collision effect is explicitly taken into account. However, the initial condition employed in the computation is slightly different from the conventional Chapman-Enskog NS distribution function. The present study also reveals the background of the existing kinetic schemes. such as the KFVS scheme and Gas-Kinetic BGK scheme.