• Journal of computational fluids engineering
• /
• v.7 no.4
• /
• pp.19-27
• /
• 2002

An efficient implicit multigrid method is presented for the Navier-Stokes and k-ω turbulence equations. Freezing and limiting strategies are applied to improve the robustness and convergence of the multigrid method. The eddy viscosity and strongly nonlinear production terms of turbulence are frozen in the coarser grids by passing down the values without update of them. The turbulence equations together with the Navier-Stokes equations, however, are consecutively solved on the coarser grids in a loosely coupled fashion. A simple limit for k is also introduced to circumvent slow-down of convergence. Numerical results for the unseparated and separated transonic airfoil flows show that all computations converge well without any robustness problem and the computing time is reduced to a factor of about 3 by the present multigrid method.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.229-232
• /
• 2002

For the purpose of Thermal Protection Material design problem, a numerical analysis of axisymmetric high temperature supersonic impinging jet flows of exhaust gas from combustor on curved surfaces has been accomplished. A modifed CSCM Upwind Navier-Stokes method which is able to cure the carbuncle Phenomena has been developed to study strong shock wave structure and thermodynamic wall properties such as pressure and heat transfer rate on various curved surfaces. The results show that the maximum heat transfer rate which is the most important parameter affecting thermo-chemical surface ablation on the plate did not occur at the center of jet impingement, but rather on a circle slightly away from the center of impingement and the shear stress distribution alone the wall is similar to the wall heat transfer late distribution.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.31-34
• /
• 2002

An efficient implicit multigrid method is presented for the Navier-Stokes and $k-{\omega}$ turbulence equations. Freezing and limiting strategies are applied to improve the robustness and convergence of the multigrid method. The eddy viscosity and strongly nonlinear production terms of turbulence are frozen in the coarser grids by passing down the values without update of them. The turbulence equations together with the Navier-Stokes equations, however, are consecutively solved on the coarser grids in a loosely coupled fashion. A simple limit for k is also introduced to circumvent slow-down of convergence. Numerical results for the unseparated and separated transonic airfoil flows show that all computations converge well without any robustness problem and the computing time is reduced to a factor of about 3 by the present multigrid method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.8 s.239
• /
• pp.921-926
• /
• 2005

In the idealized model problem of the interaction between a planar travelling shock and a symmetric vortex, the physics of shock distortion and quadrupole sound generation are well known to many researchers. However, the authors have distinguished the weak waves reflected and transmitted by the complicated photograph images obtained from a shock tube experiment. In this paper, we introduces a parametric study based on Navier-Stokes simulation and Rankin vortex model to see the difference of shock deformation shapes. Four combination of the strength of shock and vortex are respectively selected from a parameter plane of shock and vortex strength extended to the strong vortex region. The result shows clearly discernable wave morphology for the main parameters, which is not yet explicitly mentioned by other researchers.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.11 no.2
• /
• pp.157-163
• /
• 2003

This work was undertaken for the numerical analysis of defrosting phenomena of automobile windshield. To analyze the defrost, the flow and temperature field of cabin interior, heat transfer through the windshield glass, and phase change of the frost should be analyzed simultaneously. The flow field was obtained by solving the 3-D unsteady Navier-Stokes equation and the temperature field was computed by energy equation. The phase-change process of Stefan problem was solved by enthalpy method. For code validation, the temperature field of the driven cavity was calculated. The result of calculation shows a good agreement with the other numerical results. Then, the present code was applied to the defrosting analysis of a real automobile and, also, a good agreement with experiment was obtained.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.93-100
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved successively in the moving least square sense. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as a FVM.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2010.04a
• /
• pp.508-511
• /
• 2010

본 논문에서는 정상상태 유체-구조 연성문제를 연속체 기반으로 정식화하고 유한요소법을 이용하여 완전 연성된 해를 구하였다. 대 변형을 고려하기 위하여 토탈 라그란지안 정식화를 사용하였으며 유체 및 구조의 비선형성이 고려되었다. 유체와 구조 영역의 형상을 NURBS 곡면을 이용하여 매개화하여 표현하였으며, 형상 최적화를 위해 효율적인 설계민감도 해석법인 애조인 기법을 이용하여 압력, 속도, 변위 등에 대한 설계민감도를 구하였다. 이를 이용하여 최소 컴플라이언스를 갖게 하는 구조물 내부의 유체영역의 설계 등의 수치예제를 통하여 개발된 방법론의 타당성을 확인하였다.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.331-334
• /
• 2010

The supersonic speeds slowing down by shock waves is a common problem during the transonic region. So how to study the status of shock on the surface of airplane and wings is crucial adjective during transonic region. However, the theoretical and computational transonic flow problems were very hard. This paper introduced using Navier-Stokes Schemes to study characteristics of AGARD Wing 445.6 by ANSYS CFX in transonic region. From simulations results, as the Mach number increases, shock waves appear in the flowfield, getting stronger as the speed increases, these shock waves will lead to a rapid increase in drag.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.11 no.4
• /
• pp.55-68
• /
• 2007

First-order least-squares method of a distributed optimal control problem for the incompressible Navier-Stokes equations is considered. An optimality system for the optimal solution are reformulated to the equivalent first-order system by introducing velocity-flux variables and then the least-squares functional corresponding to the system is defined in terms of the sum of the squared $L^2$ norm of the residual equations of the system. The optimal error estimates for least-squares finite element approximations are obtained.

• Nuclear Engineering and Technology
• /
• v.41 no.9
• /
• pp.1191-1204
• /
• 2009

In this study, a CFD methodology is employed to address the problem of the prediction of the flow in a T-junction. An Unsteady Reynolds Averaged Navier-Stokes (URANS) approach has been selected for its low computational cost. Moreover, Unsteady Reynolds Navier-Stokes methodologies do not need complex boundary formulations for the inlet and the outlet such as those required when using Large Eddy Simulation (LES) or Direct Numerical Simulation (DNS). The results are compared with experimental data and an LES calculation. In the past, URANS has been tried on T-junctions with mixed results. The biggest limit observed was the underestimation of the oscillatory behavior of the temperature. In the present work, we propose a comprehensive approach able to correctly reproduce the root mean square (RMS) of the temperature directly downstream of the T-junction for cases where buoyancy is not present.