• 한국전산구조공학회논문집
• /
• 제16권1호
• /
• pp.95-104
• /
• 2003

본 논문은 사각형 연료 탱크 내 비점성, 비압축성, 비회전 유동에 대한 슬로싱 주파수 응답의 유한요소 해석을 다룬다. 지배방정식으로 포텐셜 이론을 기반으로 한 라플라스 방정식을 적용한다. 슬로싱 운동이 작다고 가정하여 선형화된 자유표면 조건을 적용하였고, 변수분리기법을 이용하여 이론해를 구하였다. 점성 감쇠에 따른- 에너지 소산의 영향을 구현하기 위해 가상치 점성 계수를 도입하였으며, 이고 인해 공진 주파수에서 응답의 발산을 방지할 수 있나. 슬로싱 응답의 최대 진폭을 예측하기 위해 9절점 요소를 사용한 유한요소법을 이용하여 해석하였다. 슬로싱 높이, 유체 내부 동수압 및 내부 유체력의 수치 결과는 이론해와 잘 일치하였다. 유한요소 시험 프로그램을 검증한 후, 유체높이에 따른 슬로싱 주파수 응답 특성을 분석하였다.

• 대한조선학회논문집
• /
• 제35권4호
• /
• pp.1-10
• /
• 1998

본 논문에서는 비압축성 Newtonian 점성유동에서 초기에 순간 출발하는 2차원 실린더 주위의 유동을 해석하기 위해서, 와도를 기저로 한 수치해석기법을 제안하고 있다. Helmholtz 분리 형태로 표현된 Navier-Stokes방정식에서 유도되는 와도전달방정식과 압력방정식, 그리고 벡터등식에서 유도되는 속도-와도 관계식을 이 문제의 지배방정식으로 택하고, 경계조건으로는 물체표면에서 와도와 압력의 연성관계와 힘의 평형을 고려한 동적와도경계조건과 동적압력조건이 제시된다. 이 지배방정식과 경계조건을 수치적으로 처리하기 위하여, 와도와 압력이 연성되어 있는 경계조건은 Wu등(1994)이 제안한 대로, 연성관계를 유지한 채로 식을 분리하는 방법을 이용하였고, 와도전달 방정식은 유한체적법으로 계산하였으며, 그 식에 포함된 대류항을 처리하는 방법으로 TVD 방법을 이용하였다. 속도는 Biot-Savart적분항이 포함된 벡터등식에서 panel방법으로 구하고, 압력방정식은 형태가 Poisson방정식이므로 역시 panel방법을 이용하였다. 계산에 사용된 격자로 정규격자를 이용하고, 결과를 다른 수치적, 해석적 결과와 비교하여 그 타당성을 검증하였다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.410-415
• /
• 2001

An analytical and numerical examination of second-order fractional-step methods and boundary condition for the incompressible Navier-Stokes equations is presented. In this study, the compatibility condition for pressure Poisson equation and its boundary conditions, stability, and numerical accuracy of canonical fractional-step methods has been investigated. It has been found that satisfaction of compatibility condition depends on tentative velocity and pressure boundary condition, and that the compatible boundary conditions for type D method and approximately compatible boundary conditions for type P method are proper for divergence-free velocity for type D and approximately divergence-free for type P method. Instability of canonical fractional-step methods is induced by approximation of implicit viscous term with explicit terms, and the stability criteria have been founded with simple model problems and numerical experiments of cavity flow and Taylor vortex flow. The numerical accuracy of canonical fractional-step methods with its consistent boundary conditions shows second-order accuracy except $D_{MM}$ condition, which make approximately first-order accuracy due to weak coupling of boundary conditions.

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

The numerical methodology for computing tile impact forces and water entry behaviors of high speed water entry bodies was been developed. Since the present method assumed the impact occurs within a very short time interval. the viscous effects do not have enough time to play a significant role in the impact forces, that is, the flow around a water-entry object was assumed as an incompressible potential flow and is solved by the source panel method. The elements fully submerged into the water are routinely treated, but the elements intersected by the effective planar free surface are redefined and reorganized to be amenable to the source panel method. To validate the present code, it was applied to disk, cone and ogive model and compared with experimental data. Good agreement was obtained. The water entry behavior such as the bouncing phenomena from the free surface was also simulated using the impact forces and two degree of freedom dynamic equation. Physically acceptable results were obtained.

• 대한토목학회논문집
• /
• 제13권2호
• /
• pp.267-277
• /
• 1993

본 논문에서는 유한요소법을 2차원 난류경계층에 적용하였으며, 점성유체의 시간의존 비압축성 운동을 시간과 압력장(場)에서 Navier-Stokes방정식과 vorticity방정식을 이용하여 정식 화하였다. 수치계산방법은 Galerkin방법에 기초하였으며, 난류 경계층의 eddy kinematic viscosity에 대해서는 Prandtl의 혼합거리이론을 도입하였다. 난류 경계층에서 파동에 의한 임의 저면에서 저질의 이동을 수치계산하였다. 유한 요소법에 의해 얻어진 결과는 진동흐름에 의한 경계층과 파동에 의한 경계층에서의 특성의 차이를 분명히 하였다.

• 한국전산유체공학회지
• /
• 제19권2호
• /
• pp.72-78
• /
• 2014

A hybrid particle-mesh method as the combination between the Vortex-In-Cell (VIC) method and penalization method has been achieved in recent years. The VIC method, which is based on the vorticity-velocity formulation, offers particle-mesh algorithms to numerically simulate flows past a solid body. The penalization method is used to enforce boundary conditions at a body surface with a decoupling between body boundaries and computational grids. The main advantage of the hybrid particle-mesh method is an efficient implementation for solid boundaries of arbitrary complexity on Cartesian grids. However, a numerical simulation of flows in large domains is still not too easy. In this study, a multi-domain approach is thus proposed to further reduce computation cost and easily implement it. We validate the implementation by numerical simulations of an incompressible viscous flow around an impulsively started circular cylinder.

• 한국전산유체공학회지
• /
• 제22권1호
• /
• pp.51-58
• /
• 2017

A penalized VIC method offers an efficient hybrid particle-mesh algorithm to simulate an incompressible viscous flow passing a solid body in an infinite domain. In this manner, the computational domain needs to be restricted to a relatively small region to reduce computational cost which would be very high in case of using a large domain. In this paper, we present how to dispose of far-field particles to avoid an unnecessarily large computational domain. The present approach constraints expansion of the domain and thus prevents the incremental computational cost. To validate the numerical approach, a flow around an impulsively started sphere was simulated for Reynolds numbers of 100 and 1000.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• Wind and Structures
• /
• 제1권1호
• /
• pp.111-125
• /
• 1998

This paper deals with the development of variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined element and efficiently used for the construction of a refined mesh without generating distorted elements. A modified Guassian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of the independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems in which the locations to be refined are changed in accordance with the dynamic distribution of velocity gradient, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• 대한조선학회논문집
• /
• 제50권6호
• /
• pp.421-428
• /
• 2013

This study aims at numerically investigating the angle of attack characteristics of the rudder behind a rotating propeller. The rotating propeller of 5 blades and the full spade rudder are placed in the numerical water tunnel with a uniform flow condition to consider propeller-rudder interaction. The turbulence closure model is employed to simulate the three-dimensional unsteady incompressible viscous turbulent flow around the propeller and the rudder. The present numerical method are well verified by comparing with the experimental results. In order to identify the dependence of the angle of attack of the rudder on the rudder angle, a wide range of rudder angles is considered. The present study carried out the quantitative and qualitative analysis of the angle of attack in terms of the pressure distribution, streamlines and the evaluation of the flow incidence, resulting in that the angle of attack increases as we move from the root and the tip to the center of the rudder, regardless of the rudder angle. The distribution of the angle-of-attack along the span is strongly affected by rotating propeller flow and rudder angle. Consequently, the distribution of the angle-of-attack of the oncoming flow against the rudder leading edge plays a role in determination of rudder performance.