• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• 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.

• Ocean Systems Engineering
• /
• 제3권3호
• /
• pp.219-236
• /
• 2013

The accurate prediction of motion in waves of a marine vehicle is essential to assess the maximum sea state vs. operational requirements. This is particularly true for small crafts, such as Autonomous Surface Vessels (ASV). Two different numerical methods to predict motions of a SWATH-ASV are considered: an inviscid strip theory initially developed at MIT for catamarans and then adapted for SWATHs and new a hybrid strip theory, based on the numerical solution of the radiation forces by an unsteady viscous, non-linear free surface flow solver. Motion predictions obtained by the viscous flow method are critically discussed against those obtained by potential flow strip theory. Effects of viscosity are analyzed by comparison of sectional added mass and damping calculated at different frequencies and for different sections, RAOs and motions response in irregular waves at zero speed. Some relevant conclusions can be drawn from this study: influence of viscosity is definitely non negligible for SWATH vessels like the one presented: amplitude of the pitch and heave motions predicted at the resonance frequency differ of 20% respectively and 50%; in this respect, the hybrid method with fully non-linear, viscous free surface calculation of the radiation forces turns out to be a very valuable tool to improve the accuracy of traditional strip theories, without the burden of long computational times requested by fully viscous time domain three dimensional simulations.

• 대한기계학회논문집B
• /
• 제26권12호
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• Structural Engineering and Mechanics
• /
• 제66권3호
• /
• pp.411-421
• /
• 2018

Although a rubber isolation cushion can reduce the dynamic response of a structure itself, it has little influence on the height of a sloshing wave and even may induce magnification action. Vertical baffles are set into a base-isolated Concrete Rectangular Liquid Storage Structure (CRLSS), and baffles are opened as holes to increase the energy dissipation of the damping. Problems of liquid nonlinear motion caused by baffles are described using the Navier-Stokes equation, and the space model of CRLSS is established considering the Fluid-Solid Interaction (FSI) based on the Finite Element Method (FEM). The dynamic response of an isolated CRLSS with various baffles under an earthquake is analyzed, and the results are compared. The results show that when the baffle number is certain, the greater the number of holes in baffles, the worse the damping effects; when a single baffle with holes is set in juxtaposition and double baffles with holes are formed, although some of the dynamic response will slightly increase, the wallboard strain and the height of the sloshing wave evidently decrease. A configuration with fewer holes in the baffles and a greater number of baffles is more helpful to prevent the occurrence of two failure modes: wallboard leakage and excessive sloshing height.

• 한국해양공학회지
• /
• 제33권2호
• /
• pp.131-138
• /
• 2019

This study focused on the roll motion characteristics of a two-dimensional (2D) rectangular floating structure under regular beam sea conditions. An experiment was conducted in a 2D wave tank for a roll free decay test in calm water and the roll motion in a range of regular waves with and without heave motion to investigate the motion response and heave influence on the roll motion. A numerical study was carried out using Reynolds-averaged Navier Stokes (RANS)-based CFD simulations. A grid convergence test was conducted to accurately capture the wave condition on the free surface based on the overset mesh and wave forcing method. It was found in the roll free decay test that the numerical results agreed well with the experimental results for the natural roll period and roll damping coefficient. It was also observed that the heave motion had an impact on the roll motion, and the responses of the heave and roll motion from the CFD simulations were in reasonable agreement with those from the experiment.

• 한국해양공학회지
• /
• 제34권6호
• /
• pp.406-418
• /
• 2020

In this study, a nonlinear wave simulation code is developed using a higher-order spectral (HOS) method. The HOS method is very efficient because it can determine the solution of the boundary value problem using fast Fourier transform (FFT) without matrix operation. Based on the HOS order, the vertical velocity of the free surface boundary was estimated and applied to the nonlinear free surface boundary condition. Time integration was carried out using the fourth order Runge-Kutta method, which is known to be stable for nonlinear free-surface problems. Numerical stability against the aliasing effect was guaranteed by using the zero-padding method. In addition to simulating the initial wave field distribution, a nonlinear adjusted region for wave generation and a damping region for wave absorption were introduced for wave generation simulation. To validate the developed simulation code, the adjusted simulation was carried out and its results were compared to the eighth order Stokes theory. Long-time simulations were carried out on the irregular wave field distribution, and nonlinear wave propagation characteristics were observed from the results of the simulations. Nonlinear adjusted and damping regions were introduced to implement a numerical wave tank that successfully generated nonlinear regular waves. According to the variation in the mean wave steepness, irregular wave simulations were carried out in the numerical wave tank. The simulation results indicated an increase in the nonlinear interaction between the wave components, which was numerically verified as the mean wave steepness. The results of this study demonstrate that the HOS method is an accurate and efficient method for predicting the nonlinear interaction between waves, which increases with wave steepness.

• 한국전산구조공학회논문집
• /
• 제34권1호
• /
• pp.59-69
• /
• 2021

최근 친환경 에너지 개발에 대한 관심의 증가로 해상 및 연안 지역에서 대규모 해양구조물들이 건설되고 있다. 해양구조물은 항상 파랑 하중에 노출되어 있으므로 구조적인 안전성을 확보하기 위해서는 파랑에 대한 정확한 이해와 분석이 필수적이다. 실해역에서 수행되는 실험은 해양파를 이해하기 위한 가장 정확한 방법이지만, 변수의 통제가 어렵고 비용과 규모 측면에서 실험이 제한되는 경우가 많다. 본 연구에서는 수치파동수조를 이용하여 다양한 조건의 파를 생성하고 및 이론식과 비교를 통해 파랑 생성 능력을 검증하였다. 입자 기반 수치해석법인 SPH(Smoothed Particle Hydrodynamics) 기법을 이용하여 3차원 수조 및 피스톤 조파기를 모델링하였으며, 반사파에 대한 영향을 최소화하기 위해 수로 끝단에 질량 가중 감쇠 영역을 설정하여 안정적인 파고 및 유속 계산이 수행될 수 있게 하였다. 목표 파랑 조건에서,상대 수심이 2 이하를 만족하는 경우 파형경사에 관계없이 파고와 유속을 계산한 결과가 이론값과 높은 정확도를 보였다. 그러나 상대수심과 파형경사의 목표값이 증가하고, 측정 위치가 멀어짐에 따라서 최대 10% 이상의 오차가 발생하였다. 수치해석을 이용하여 정확한 계산이 가능한 파랑 범위를 무차원 변수를 이용하여 제안하였으며, 차후 수치해석을 이용한 수치파동수조 검증기준과 유체-구조물 상호작용 해석분야 연구에 효과적으로 활용될 수 있다.

• 대한기계학회논문집A
• /
• 제24권1호
• /
• pp.16-28
• /
• 2000

The steady rotating flows in eccentric annulus has been studied by a numerical method based on the spectral collocation method. The inner cylinder has a constant angular velocity while the outer on e is stationary. Flow between eccentric cylinders is of considerable technical importance as it occurs in journal bearings. In the present work, the governing equations for laminar flow are expressed as Navier-Stokes equations, including the non-linear convection terms. The solutions were utilized i, estimate the effects of the nonlinear terms on the load acting on the rotating cylinder. Based on the half and the full Sommerfeld methods, the load on the rotating cylinder is evaluated with eccentricity, by integrating the pressure and skin friction around the cylinder. The attitude angle and Sommerfeld reciprocal are calculated from the load. Also, the torque on the rotating inner cylinder was calculated. considering the skin friction. The attitude angle and Sommerfeld reciprocal are decreased with eccentricity. Viscous damping coefficient due to the skin friction becomes larger with decreasing the annular space. It is found the non-linear effects of the convection terms on the flow and the load are important. especially on the attitude angle, for relatively wide annular configurations however, the effects on those are minor for very narrow annular ones.

• 소음진동
• /
• 제11권3호
• /
• pp.428-436
• /
• 2001

The rotating inner cylinder executes a periodic translational motion in concentric annulus while the outer one is stationary. In the study of flow-induced vibrations and relaxed instabilities, it is of interest to evaluate the fluid-dynamic forces acting on the rotating inner cylinder. In the present work, the governing equations for the confined flow are expressed as Navier-Stokes equations, including the steady and unsteady terms. The fluid parameters for steady flow generated by the rotating cylinder are determined analytically while the unsteady ones by the oscillatory motion are evaluates by a numerical method based on the spectral collocation method. In order to validate the numerical approach, the numerical results are compared wish the analytical ones given by existing theories, for simple cases where the both approaches are applicable. Good agreement was found between the results. It is found the effects of the Reynolds number, defined by rotating velocity, on the fluid-dynamic forces are important for the case of relatively low oscillatory Reynolds number, defined by oscillatory frequency : j.e., in case of $Re_\omega\gg Re_S$.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2007년도 춘계학술대회논문집
• /
• pp.275-282
• /
• 2007

In this study, a fluid-structure interaction (FSI) analysis system has been developed in order to evaluate the turbine cascade performance with blade structural deformation effect. Relative movement of the rotor with respect to stator is reflected by modeling independent two computational domains. To consider the deformed position of rotor airfoil, dynamic moving grid method is applied. Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation SST $k-{\varepsilon}$ turbulence models are solved to predict unsteady fluid dynamic loads. A fully implicit time marching scheme based on the Newmark direct integration method with high artificial damping is used to compute the fluid-structure interaction problem. Cascade performance evaluations for different elastic axis positions are presented and compared each other. It is importantly shown that the predicted aerodynamic performance considering structural deformation effect of blade can show some deviations compared to the data generally computed from rigid blade configurations and the position of elastic axis also tend to give sensitive effect.