• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.4
• /
• pp.379-388
• /
• 2008

A hybrid Cartesian/immersed boundary code is expanded to simulate flow field around a three-dimensional body which undergoes large dynamic deformation. Immersed boundary nodes are automatically distributed based on the edges crossing triangles on body boundary. Velocity vectors are reconstructed at those immersed boundary nodes along local normal lines to the boundary. The reconstruction of pressure is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other results on laminar flow over a sphere. The code is applied to simulate flow around a foil which is attached to a body of revolution and rotates under periodic deformation. The periodic variation of the tip vortex is observed and the effects of the deformation on hydrodynamic force acting on the body are investigated.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.2
• /
• pp.79-87
• /
• 2013

In this paper, an Euler equation solver based on a Cartesian-grid method and non-uniform staggered grid system is applied to predict the ship motion response and added resistance in waves. Water, air, and solid domains are identified by a volume-fraction function for each phase and in each cell. For capturing the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with a weighed line interface calculation (WLIC) method. The volume fraction of solid body embedded in a Cartesian-grid system is calculated by a level-set based algorithm, and the body boundary condition is imposed by volume weighted formula. Added resistance is calculated by direct pressure integration on the ship surface. Numerical simulations for a Wigley III hull and an S175 containership in regular waves have been carried out to validate the newly developed code, and the ship motion responses and added resistances are compared with experimental data. For S175 containership, grid convergence test has been conducted to investigate the sensitivity of grid spacing on the motion responses and added resistances.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.109-114
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.1
• /
• pp.139-149
• /
• 1994

A numerical simulation of unsteady axisymmetric turbulent flow was performed for a reciprocating engine including port/valve assembly. The governing equations based on a nonorthogonal coordinate formulation with Cartesian velocity components were used and discretised by the finite volume method with non-staggered variable arrangements. The modified $\kappa-\xi$. turbulence model which included the effect of compressibility was used. The results of twodimensional transient calculation for the axisymmetric configuration were compared with the experimental data. Although slightly low rms velocity was predicted compared to the experimental data, predicted velocity distributions at the valve exit and in-cylinder region showed good agreements with the experimental data. The flow at the valve exit was separated at the same valve lift position with the experimental data. Two vortices incylinder region were generated during the initial intake process. The clockwise main vortex became strong and moved upward to the top wall. The counter-clockwise second vortex became weak and stick to the upper left corner of the cylinder. After middle intake process, new vortex adjacent to upper cylinder wall appeared by the piston motion and therefore, the in-cylinder flow was formed into three vortices. The cylinder pressure just before bottom dead center of piston was higher than inlet pressure and then the reverse flow occured at the valve exit. The in-cylinder flow characteristics were strongly dependent on piston motion, but insensitive to valve motion.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.395-406
• /
• 2011

In this research a numerical simulation method is developed for moving body in free surface flows using fixed staggered rectangular grid system. The non-linear free surface near the body is defined by marker-density method. The body boundary is defined by line segment connecting the points where the body surface and grid line meet. Continuity equation and Navier-Stokes equations are used as governing equations and the equations are coupled with two-step projection method. The velocities and pressures of body boundary and free surface cells are calculated with simultaneous iterative method. To treat a body movement in a fixed grid system, the volume displaced by moving body is added to the divergence of the body boundary cell. For the verification of the present numerical method. vortex shedding period of advancing cylinder is calculated and the period is compared with existing experiment results. Moreover, added mass and damping coefficients of a vertically excited box are calculated and the computed results are compared with published experiment results. Impulsive pressure and water level variation due to sloshing phenomenon are simulated and the results are compared with published experiment results. Varying the plunger shape, the waves generated by plunging type wave maker are compared with the 2nd order Stokes wave theory The plunger shape generating the wave that shows the best agreement with the theory is represented.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.6
• /
• pp.619-627
• /
• 2004

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method. The results agree well with previous numerical and experimental results. This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications. The usefulness of this method will be confirmed when we solve the complex geometries and moving bodies.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.89-99
• /
• 1998

한국형 차세대원자로에서는 비상노심 안전주입수가 저온관을 통하지 않고 원자로용기에 직접 주입된다. 원자로용기의 가압열충격과 열수력적 관점에서 최적의 노즐위치를 결정하기 위해서 전산유체역학을 활용하였다. 상용 전산유체코드인 CFX를 이용하여 원자로 하향유로를 모사하는 해석대상 격자를 다중불록으로 형성한 다음 유동장을 비압축성 Navier-Stokes 운동량 방정식, 에너지 방정식과 표준 k-ε 난류모형 등으로 모형화하여 3차원 비정상상태 계산을 수행하였다. CFX에서는 경계 밀착좌표계, 비엇물림격자와 SIMPLE 알고리즘을 사용한다. 본 연구결과 원자로용기의 가압열충격 관점에서 가장 보수적인 사고인 증기관 파단사고시에도 열적혼합이 잘 일어나 가압열충격이 발생할 가능성이 없는 것으로 판단되며 안전주입수 노즐이 저온관 바로 위에 위치할 때 원자로 하향유로 내의 온도 분포가 가장 균일하여 열적 혼합 관점에서는 최적의 위치로 판단된다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.2
• /
• pp.198-209
• /
• 1999

This study reports the selection of dependent variables for momentum equations in general curvilinear coordinates. Catesian, covariant and contravariant velocity components were examined for the dependent variable. The focus of present study is confined to staggered grid system Each dependent variable selected for momentum equations are tested for several flow fields. Results show that the selection of Cartesian and covariant velocity components intrinsically can not satisfy mass conservation of control volume unless additional converting processes ore used. Also, Cartesian component can only be used for the flow field in which main-flow direction does not change significantly. Convergence rate for the selection of covariant velocity component decreases quickly as with the increase of non-orthogonality of grid system. But the selection of contravariant velocity component reduces the total mass residual of discretized equations rapidly to the limit of machine accuracy and the solutions are insensitive to the main-flow direction.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1500-1508
• /
• 2001

This paper represents the importance of dependent variables in non-orthogonal curvilinear coordinates just as the importance of those variables of convective scheme and turbulence model in computational fluid dynamics. Each of Cartesian, physical covariant and physical contravariant velocity components was tested as the dependent variables of momentum equations in the staggered grid system. In the flow past a circular cylinder, the results were computed to use each of three variables and compared to experimental data. In the skewed driven cavity flow, the results were computed to check the grid dependency of the variables. The results used in Cartesian and physical contravariant components of velocity in cylinder flow show the nearly same accuracy. In the case of Cartesian and contravariant component, the same number of vortex was predicted in the skewed driven cavity flow. Vortex strength of Cartesian component case has about 30% lower value than that of the other two cases.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.916-921
• /
• 2003

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method(FVM). This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications.