• Journal of Ocean Engineering and Technology
• /
• v.7 no.2
• /
• pp.150-161
• /
• 1993

원초변수를 이용한 Navier-Stokes 방정식의 수치계산기법을 개발하고, 이를 응용하여 backward facing step의 층류 유동을 계산하였다. 직교좌표계에서의 비압축성 Navier-Stokes방정식을 풀기위해 시간과 공간항을 2차 정도의 유한 차분을 사용하여 이산화하였고 비교차격자계를 사용하여 양해법으로 수치 계산하였다. 운동량방정식과 연속방정식으로 부터 유도된 압력방정식(pressure-poisson equation)을 이용하여 무발산 조건을 만족시켰ㄲ다. Backward facing step의 층류 유동을 100.$\leq$R$_e$$\leq$1000 범위에 대해서 수치 계산하였으며 실험결과와 잘 일치하는 결과를 구할 수 있었다. 특히 step뒤에서 생기는 박리구간의 길이는 다른 계산결과들보다 실험치에 가까운 값을 얻을 수 있었으며, Re가 600보다 클때는 위쪽 벽에 또 다른 박리 유동이 발생되는 현상이 예측되었다.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.2
• /
• pp.150-160
• /
• 2002

The free vibration characteristics of the triple cylindrical shells filled with fluid are investigated. The triple cylindrical shells are filled with compressible fluid. The boundary condition is clamped at both ends. Analytical method is developed to evaluate natural frequencies of triple cylindrical shells using Sanders' shell theory and courier series expansion by Stokes' transformation. Their results are compared with those of finite element method to verify the validation of the method developed. The modal characteristics of shells filled with fluid at region 1, 2 and 3 are evaluated.

• Proceedings of the Korean Society of Computer Information Conference
• /
• 2023.07a
• /
• pp.585-586
• /
• 2023

본 논문에서는 바람의 이류(Advection)를 고려하여 사운드의 전파를 변형하는 방법을 제시한다. 사운드는 공기와 같은 매질의 진동을 통해 전파되는 파동이며, 이런 바람의 이동 방향은 사운드 에너지 전파에 직접적인 영향을 주며, 본 논문에서는 이를 광선추적법(Raytracing) 기반으로 모델링한다. 기존의 사운드 전파는 물리기반, 기하처리(Geometry processing), 혼합기법(Hybrid method) 등의 방법이 제안됐으며, 다양한 장면에서 좋은 결과를 만들어냈다. 하지만 바람의 움직임은 유체역학을 기반으로 한 나비에-스토크스 방정식(Navier-Stokes equation)에 의해 표현되기 때문에 사운드 전파만으로는 바람의 영향을 고려한 전파 형태를 모델링할 수 없다. 본 논문에서는 바람의 유동 중 이류를 고려하여 사운드 맵을 효율적으로 변형할 수 있는 방법을 제시한다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.

• Journal of Navigation and Port Research
• /
• v.27 no.5
• /
• pp.473-478
• /
• 2003

Numerical study was made on the flow characteristics around a circular pipeline between parallel walls. The incompressible Navier-Stokes equations were solved by using a third-order upwind differential scheme. When the distance near a wall is small enough, the vortex shedding is almost completely suppressed because of the interaction with the wall boundary layer separation. This study aims to clarify the characteristics of the vortex shedding regime as the body approaches a wall as Reynolds number varies. The feature of separated vorticity dynamics is analyzed at different conditions with particular attention to the interaction between the pipeline wake and the induced separation on the plane walls.

• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.575-581
• /
• 1988

To integrate the two-dimensional Navier-Stokes equations numerically in multiply-connected flow regions, the vorticity-stream function formulation is used. The steady stream function value at the surface of the multibody, initially unknown, has been determined interactively by introducing a line integral which requires the single-valuedness of pressure at each interaction step. This procedure is relatively simpler and more efficient than the primitive variable formulation which requires much more computing time and shows poor convergence. Three doubly-connected flow problems are defined and numerically analyzed by the present method. The results have been compared either with earlier existing ones or with the experimental interferograms to demon-strate the validity of the presented method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1240-1247
• /
• 1999

Particle motion through a disk type critical orifice placed in a 3.0cm diameter chamber has been studied numerically. In the simulation, the velocity field is solved using Pantankar's SIMPLER algorithm for the compressible flow and convergence of the computation is confirmed if the mass source at each control volume is smaller than $10^{-7}$. The particle motion in the flow field is solved in Lagrangian method. The particle trajectories showed that the particles injected away from the center line are expanded rapidly. At lower pressures, this expansion phenomena are more dominant. At lower pressures, the clear difference in particle and air speed is showed all the way down to the exit plan. It was found that particles with Stokes number of ca.2.5 tend to focus close to the center line very well except the particles travelling near the wall. However, particles with Stokes number greater than ca.2.5 show a tendency to cross the center line.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.93-114
• /
• 2019

In this article, we introduce a finite difference method for solving the Navier-Stokes equations in rectangular domains. The method is proved to be energy stable and shown to be second-order accurate in several benchmark problems. Due to the guaranteed stability and the second order accuracy, the method can be a reliable tool in real-time simulations and physics-based animations with very dynamic fluid motion. We first discuss a simple convection equation, on which many standard explicit methods fail to be energy stable. Our method is an implicit Runge-Kutta method that preserves the energy for inviscid fluid and does not increase the energy for viscous fluid. Integration-by-parts in space is essential to achieve the energy stability, and we could achieve the integration-by-parts in discrete level by using the Marker-And-Cell configuration and central finite differences. The method, which is implicit and second-order accurate, extends our previous method [1] that was explicit and first-order accurate. It satisfies the energy stability and assumes rectangular domains. We acknowledge that the assumption on domains is restrictive, but the method is one of the few methods that are fully stable and second-order accurate.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.575-589
• /
• 2008

In the present work, we investigate the hydrodynamic behavior of a turbulent bore, such as tsunami bore and tidal bore, generated by the removal of a gate with water impounded on one side. The bore generation system is similar to that used in a general dam-break problem. In order to the numerical simulation of the formation and propagation of a bore, we consider the incompressible flows of two immiscible fluids, liquid and gas, governed by the Navier-Stokes equations. The interface tracking between two fluids is achieved by the volume-of-fluid (VOF) technique and the M-type cubic interpolated propagation (MCIP) scheme is used to solve the Navier-Stokes equations. The MCIP method is a low diffusive and stable scheme and is generally extended the original one-dimensional CIP to higher dimensions, using a fractional step technique. Further, large eddy simulation (LES) closure scheme, a cost-effective approach to turbulence simulation, is used to predict the evolution of quantities associated with turbulence. In order to verify the applicability of the developed numerical model to the bore simulation, laboratory experiments are performed in a wave tank. Comparisons are made between the numerical results by the present model and the experimental data and good agreement is achieved.