• 한국전산유체공학회지
• /
• 제10권4호통권31호
• /
• pp.24-31
• /
• 2005

A dynamic simulation on the missile staging system is conducted with numerical techniques. Both Euler equations and Navier-Stokes equations are numerically solved respectively. The dynamic simulation of two moving bodies is fully integrated into the computational fluid dynamics solution procedure. The Chimera grid scheme is applied in this simulation for unsteady supersonic flow analysis with dynamic modeling. The objective of the study is to investigate the problem pertaining to possible unstability in missile staging. In addition, the computational comparison between in viscid and viscid flow solvers is also performed in this study.

• 한국해양공학회지
• /
• 제29권6호
• /
• pp.427-435
• /
• 2015

This paper provides the numerical results of a simulation of the flow around a propeller working beneath the free surface. A finite volume method is used to solve the unsteady Reynolds averaged Navier-Stokes (URANS) equations, where the wave-making problem is solved using a volume-of-fluid (VOF) method. The numerical analysis focuses on the propeller wake structure affected by the free surface, where we consider another free surface boundary condition that treats the free surface as a rigid wall surface. The propeller wake under the effect of these two free surface conditions shows a reduction in the magnitude of the longitudinal and vertical flow velocities, and its vortical structures strongly interact with the free surface. The thrust and torque coefficient under the free surface effect decrease about 3.7% and 3.1%, respectively. Finally, the present numerical results show a reasonable agreement with the available experimental data.

• 한국자동차공학회논문집
• /
• 제11권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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.389-396
• /
• 2010

In this paper, a new computational code was developed using Chorin's artificial compressibility method to solve the two-dimensional incompressible Navier-Stokes equations. In spatial derivatives, Roe's flux difference splitting was used for the inviscid flux, while central differencing was used for the viscous flux. Furthermore, AF-ADI with dual time stepping method was implemented for accurate unsteady computations. Two-equation turbulence models, Menter's $k-{\omega}$ SST model and Coakley's $q-{\omega}$ model, hae been adopted to solve high-Reynolds number flows. A number of numerical simulations were carried out for steady laminar and turbulent flow problems as well as unsteady flow problem. The code was verified and validated by comparing the results with other computational results and experimental results. The results of numerical simulations showed that the present developed code with the artificial compressibility method can be applied to slve steady and unsteady incompressible flows.

• 한국전산유체공학회지
• /
• 제4권1호
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2006년도 추계학술대회논문집
• /
• pp.848-854
• /
• 2006

In this study, a computational analysis system has been developed in order to investigate flow-induced vibration(FIV) phenomenon for general stator-rotor cascade configurations. Relative movement of the rotor with respect to stator is reflected by modeling independent two computational domains. Fluid domains are modeled using the unstructured grid system with dynamic moving and local deforming methods. Unsteady, Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation SST $k-\omega$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used flow computing the coupled governing equations of the fluid-structure interaction problem. Detailed FIV responses for different flow conditions are presented with respect to time and vibration characteristics are also physically investigated in the time domain.

• Korea-Australia Rheology Journal
• /
• 제17권2호
• /
• pp.47-62
• /
• 2005

The present study deals with a mathematical model describing the dynamic response of heat and mass transfer in blood flow through bifurcated arteries under stenotic condition. The geometry of the bifurcated arterial segment possessing constrictions in both the parent and the daughter arterial lumen frequently appearing in the diseased arteries causing malfunction of the cardiovascular system, is formulated mathematically with the introduction of the suitable curvatures at the lateral junction and the flow divider. The blood flowing through the artery is treated to be Newtonian. The nonlinear unsteady flow phenomena is governed by the Navier-Stokes equations while those of heat and mass transfer are controlled by the heat conduction and the convection-diffusion equations respectively. All these equations together with the appropriate boundary conditions describing the present biomechanical problem following the radial coordinate transformation are solved numerically by adopting finite difference technique. The respective profiles of the flow field, the temperature and the concentration and their distributions as well are obtained. The influences of the stenosis, the arterial wall motion and the unsteady behaviour of the system in terms of the heat and mass transfer on the blood stream in the entire arterial segment are high­lighted through several plots presented at the end of the paper in order to illustrate the applicability of the present model under study.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.55-58
• /
• 2002

Parallel finite element code has been recently developed for the analysis of the incompressible Wavier-Stokes equations using domain decomposition method. Metis and MPI libraries are used for the domain partitioning of an unstructured mesh and the data communication between sub-domains, respectively. For unsteady computation of the incompressible Navier-Stokes equations, 4-step splitting method is combined with P1P1 finite element formulation. Smagorinsky and dynamic model are implemented for the simulation of turbulent flows. For the validation performance-estimation of the developed parallel code, three-dimensional Laplace equation has been solved. It has been found that the speed-up of 40 has been obtained from the present parallel code fir the bench mark problem. Lastly, the turbulent flows around the MIRA model and Tiburon model have been solved using 32 processors on IBM SMP cluster and unstructured mesh. The computed drag coefficient agrees better with the existing experiment as the mesh resolution of the region increases, where the variation of pressure is severe.

• 대한조선학회논문집
• /
• 제32권3호
• /
• pp.62-71
• /
• 1995

비압축성 가정하에 Navier-Stokes 방정식을 이용하여 비정상 점성유동을 수치해석하기 위해서는 매시간 단계에서 타원형 압력 Poisson 방정식의 해를 구해야 하며, 이에 많은 계산시간이 소요된다. 본 논문에서는 직접법을 이용하여 압력 Poisson 방정식을 수치해석하였으며, 분할수 증가에 따른 소요시간 문제를 다루었다. Green 정리를 압력 Poisson 방정식에 적용하면 주어진 문제는 경계치문제로 변환되고, convolution 형의 영역적분은 F.F.T.를 이용하여 계산시간을 단축할 수 있어, 직접법 이용시 소요시간은 경계치문제의 해를 구하는 데에 좌우된다. 직접법의 검증을 위하여 해석해를 알고 있는 경우에 대하여 수치해석하였고, 물체경계조건과 정합문제에 관하여 수치해석 하였는데. 분할수가 (n.n) 시 O($n^{3}$) 미만의 계산시간으로 수치해석할 수 있었다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.251-255
• /
• 2008

Flow field around helicopter involves incompressible flow near the blade root and compressible flow at the blade tip. A problem occurs for low Mach number flow due to the stiffness of the governing equations. Time-derivative preconditioning techniques have been incorporated to reduce the stiffness that occurs at low speed region. The preconditioned form of the compressible Navier-Stokes and Euler equations is used. Computations are performed for the Caradonna-Tung's hovering and non-lifting forward flight case. Computational results are in good agreement with the experimental data.