• 대한조선학회논문집
• /
• 제38권1호
• /
• pp.1-8
• /
• 2001

본 논문에서는 비압축성 Reynolds-Averaged Navier-Stokes 방정식을 수치 해석하여 자유표면을 포함한 선체 주위의 난류 유동을 계산하였다. 정규격자 상에서 공간의 이산화는 2차 정도의 유한차분법을, 시간의 적분에는 4단계 Runge-Kutta법을 이용하였고, 난류 닫힘 조건을 만족시키기 위해 Baldwin-Lomax 난류 모형을 사용하였다. 자유표면의 위치는 운동학적 경계조건식을 Lax-Wendroff법으로 풀어서 구하였고, 자유표면과 격자 경계면을 일치시키기 위해 매 시간마다 새로 계산된 자유표면 위치에 맞추어 격자를 새로 구성하였다. 속도와 압력에 대한 경계조건은 자유표면에서 점성을 무시하여 근사한 동역학적 조건을 적용해서 구하였다. 본 연구에서 개발된 수치해법을 검증하기 위하여 실험자료가 많은 Wigley 선형과 Sries 60 $C_B=0.6$ 선형에 대해 수치계산을 수행하였고 계산된 선체 주위의 파형이 실험 결과와 잘 일치하는 것을 확인하였다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제23권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.

• 한국전산유체공학회지
• /
• 제13권3호
• /
• pp.69-78
• /
• 2008

The high-order numerical method based on the adaptive mesh refinement(AMR) on the quadrilateral unstructured grids has been developed in this paper. This adaptive-grid method, originally developed with MUSCL-TVD scheme, is now extended to the WENO (weighted essentially no-oscillatory) scheme with the Runge-Kutta time integration of fifth order in spatial and temporal accuracy. The multidimensional interpolation was studied in the preliminary research, which allows us to maintain the same order of accuracy for the computation of numerical flux between two adjacent cells of different levels. Some standard benchmark tests are done to validate this method for checking the overall capacity and efficiency of the present adaptive-grid technique.

• Journal of Ship and Ocean Technology
• /
• 제5권1호
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.

• 동력기계공학회지
• /
• 제12권2호
• /
• pp.35-40
• /
• 2008

This paper deals with the free vibrations of conical shells with linearly variable thickness by the transfer influence coefficient method. The classical thin shell theory based upon the Flugge theory is assumed and the governing equations of a conical shell are written as a coupled set of first order matrix differential equations using the transfer matrix. The Runge-Kutta-Gill integration method is used to solve the governing differential equation. The natural frequencies and corresponding mode shapes are calculated numerically for the conical shells with linearly variable thickness and various boundary conditions at the edges. The present method is applied to conical shells with linearly varying thickness, and the effects of the semi-vertex angle, the number of circumferential waves and thickness ratio on vibration are studied.

• 항공우주시스템공학회지
• /
• 제10권1호
• /
• pp.1-7
• /
• 2016

During flight of the aircraft, the vortex merging phenomenon appears under the certain condition between co-rotating vortices which were generated at the wing tip and lifting-surface. And then these merged vortices at both sides show counter-rotating pattern to affect on the downstream of the aircraft. In this paper, the numerical simulations are conducted assuming this phenomenon in two-dimensional co-rotating or counter-rotating vortices pairs. Two-dimensional incompressible Navier-Stokes equations were converted into Vorticity-Streamfunction form and the Galerkin spectral method was adopted. The third order Runge-Kutta method was used for time integration. The effects on the vortex merger and degree of vortex merger were investigated according to time, Reynolds number, and changes in the distance between two vortices.

• 소음진동
• /
• 제2권3호
• /
• pp.203-211
• /
• 1992

The main purpose of this paper is to present both the natural frequencies and the first critical loads of tapered columns. The ordinary differential equation governing the free vibration for tapered columns under compressive axial force is derived. Three kinds of cross sectional shape are considered in the governing equation. The Runge-Kutta method and determinant search method are used to perform the integration of the differential equation and to determine the natural frequencies, respectively. Additionally, the bisection method is used to determine the critical loads. In numerical examples, the effects of compressive axial force on the natural frequencies of tapered columns are investigated varying the end conditions. The first critical loads of tapered columns are determined on the basis of dynamic concepts. The first critical loads of tapered columns are determined on the basis of dynamic concept. The effects of cross sectional shapes are shown and some typical mode shapes are also presented.

• 동력기계공학회지
• /
• 제9권3호
• /
• pp.58-65
• /
• 2005

This paper deals with the free vibrations of truncated conical shell with uniform thickness by the transfer influence coefficient method. The classical thin shell theory based upon the $Fl\ddot{u}gge$ theory is assumed and the governing equations of a conical shell are written as a coupled set of first order differential equations using the transfer matrix. The Runge-Kutta-Gill integration and bisection method are used to solve the governing differential equations and to compute the eigenvalues respectively. The natural frequencies and corresponding mode shapes are calculated numerically for the truncated conical shell with any combination of boundary conditions at the edges. And all boundary conditions and the intermediate supports between conical shell and foundation could be treated only by adequately varying the values of the spring constants. Numerical results are compared with existing exact and numerical solutions of other methods.

• 한국윤활학회:학술대회논문집
• /
• 한국윤활학회 1997년도 제25회 춘계학술대회
• /
• pp.183-189
• /
• 1997

This paper presents the motion of dynamically loaded journal in the connecting rod bearing of reciprocating internal combustion engine. Journal motions in engine bearings have been composed of two components, which was rotational and translational motion. Early study of journal locus in engine bearing had been performed on each motion. This paper has been considered two motions simultaneously. Reynolds equation including the squeeze effect has been analyzed using the ADI method, and real engine bearing and crankshaft system has been considered to calculate the cyclic external force. The equations of journal motion have been derived and then the numerical integration of these equations performed by 4th order Runge-Kutta method. This paper gives various journal orbits in connecting rod bearing depending on cyclic external forces, rotating speeds, and bearing parameters.

• 한국해양공학회지
• /
• 제3권2호
• /
• pp.551-551
• /
• 1989

In general, a remotely operated vehicle(ROV) operates at deep sea. The control system of ROV is composed of two local loops; the first loop placed on the surface vessel monitors and manipulates the attitude of the ROV using joystick, and the second part on the ROV automatically controls thrusters and acquires positional data. This paper presents a position control simulation of a ROV using an adaptive controller and discusses the control effects of two different conditions. The design of an adaptive control system is obtained by the application of a self-tuning controller with the minimization of an appropriate cost function. The parameters of the control system are estimated by a recursive least square method(RLS). In the simulation, a Runge-Kutta method is used for the numerical integration and the generated outputs are obtained by adding measurement errors. Additionally, this paper discusses the mathematical modelling of a ROV and make a survey of control systems.