• 한국복합재료학회:학술대회논문집
• /
• 한국복합재료학회 2000년도 춘계학술발표대회 논문집
• /
• pp.126-130
• /
• 2000

In this study, flutter characteristics of a composite wing have been studied for the variation of laminate angles in the subsonic, transonic and supersonic flow regime. The laminate angles are selected by the aspect of engineering practice such as 0, $\pm$45 and 90 degrees. To calculate the unsteady aerodynamics for flutter analysis, the Doublet Lattice Method(DLM) in subsonic flow and the Doublet Point Method(DPM) in supersonic flow are applied in the frequency domain. In transonic flow, transonic small disturbance(TSD) code is used to calculate the nonlinear unsteady aerodynamics in the time domain. Aeroelastic governing equation has been solved by v-g method in the frequency domain and also by Coupled Time-Integration Method(CTIM) in the time domain. from the results of present study, characteristics of free vibration responses and aeroelastic instabilities of a composite wing are presented for the set of various lamination angles in the all flow range.

• 대한기계학회논문집B
• /
• 제35권6호
• /
• pp.617-624
• /
• 2011

본 연구에서는 사각용기 내에서의 자유표면을 갖는 2차원 슬로싱 문제에 대하여 고찰하였다. 용기 내부의 유동은 수평방향의 조화운동을 통해 형성되며, 외부 가진 속도는 u=Asin($2{\pi}ft$)으로 정의 된다. 여기서 u, A�� 그리고 f는 외부로부터 작용하는 가진 속도, 변위 그리고 주파수를 각각 나타낸다. 큰 변위(A~O(1)) 슬로싱 문제의 해석을 위한 실험설비를 구축하였으며, 광범위한 영역에서의 PIV실험을 수행하였다. 실험을 통해 자유표면의 움직임(motion)을 각각 서로 다른 물리적 특성을 갖는 세 가지; 선형 슬로싱의 자유표면의 움직임과 유사한 standing wave motion, 사각용기의 속도가 0을 지나는 순간(the moment of turn-over) 측벽에서 발생하는 run-up motion 그리고 측벽으로부터 내부유체로 점차적으로 전파되는 bore motion으로 분류하였다.

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.13-20
• /
• 1997

A fundamental study for the numerical scheme to simulate unsteady nonlinear waves by solving Euler equations is presented. First a conservation form and a non-conservation form of the Euler equations with a free surface fitted coordinate system are compared. Next, a time splitting fractional step method and an alternating direction implicit(ADI) method for the time integration are compared. For the comparative study, flow calculations around a bottom bump in a channel and a NACA 0012 hydrofoil in a flume are performed. The results show that the ADI method with a third order upwind differencing scheme is very efficient in reducing the computing time with keeping the accuracy, And, there is no distinct difference between two expression forms except that the non-conservative form shows faster wave propagating velocity than the conservation form. Some results are compared with experiments and show good agreement.

• 대한수학회지
• /
• 제49권3호
• /
• pp.585-604
• /
• 2012

In this paper we consider the evolution of the rolling stone with a rotationally symmetric nonconvex compact initial surface ${\Sigma}_0$ under the Gauss curvature flow. Let $X:S^n{\times}[0,\;{\infty}){\rightarrow}\mathbb{R}^{n+1}$ be the embeddings of the sphere in $\mathbb{R}^{n+1}$ such that $\Sigma(t)=X(S^n,t)$ is the surface at time t and ${\Sigma}(0)={\Sigma}_0$. As a consequence the parabolic equation describing the motion of the hypersurface becomes degenerate on the interface separating the nonconvex part from the strictly convex side, since one of the curvature will be zero on the interface. By expressing the strictly convex part of the surface near the interface as a graph of a function $z=f(r,t)$ and the non-convex part of the surface near the interface as a graph of a function $z={\varphi}(r)$, we show that if at time $t=0$, $g=\frac{1}{n}f^{n-1}_{r}$ vanishes linearly at the interface, the $g(r,t)$ will become smooth up to the interface for long time before focusing.

• International Journal of Ocean System Engineering
• /
• 제3권1호
• /
• pp.1-9
• /
• 2013

A design procedure for a ship with minimum resistance had been developed using a numerical optimization method called SQP (Sequential Quadratic Programming) combined with computational fluid dynamics (CFD) technique. The frictional resistance coefficient was estimated by the ITTC 1957 model-ship correlation line formula and the wave-making resistance coefficient was evaluated by the potential-flow panel method with the nonlinear free surface boundary conditions. The geometry of the hull surface was represented and modified by B-spline surface modeling technique during the optimization process. The Series 60 ($C_B$=0.60) hull was selected as a parent hull to obtain an optimized hull that produces minimum resistance. The models of the parent and optimized hull forms were tested at calm water condition in order to demonstrate the validity of the proposed methodolgy.

• 대한조선학회논문집
• /
• 제46권6호
• /
• pp.587-594
• /
• 2009

In the paper the effect of a stern-plate attached to a ship was taken into account. The relationship between the trim angle of a ship and the wave-resistance coefficient induced by the a stern-plate was studied using the potential flow analysis method. Numerical algorithm was described using the panel method and the vortex lattice method(VLM) to simulate the flow phenomena around a ship. The non-linearity of the free surface boundary conditions were considered using the iterative method and the IGE-GMRES(Incomplete Gaussian Elimination-The Generalized Minimal RESidual) algorithm was adopted to solve the linear equation at each iterative step. Numerical calculations were carried out to investigate the validity of the adopted algorithm using KCS(KRISO 3600 TEU Container) hull. Possible cases for attachment of the plate were checked. The results showed that the numerical algorithm could be physically appropriate.

• 대한조선학회지
• /
• 제27권3호
• /
• pp.19-30
• /
• 1990

본 논문에서는 비선형 자유표면파의 수치해법 개발을 위한 노력의 일환으로서, 하밀톤 원리에 근거한 비정상 자유표면파 문제에 대한 범함수를 정의하였다. 즉, 고전적인 하밀톤 원리에 자유표면 경계조건에 따른 수정항을 고려해 줌으로써 본 문제에 적합한 두 범함수를 만들 수 있다, 본 논문에서는 첫째, 선형문제에서의 고유치 문제를 해석하여 두 범함수가 선형 분산관계식의 상한과 하한치를 가지는 것을 보여준다. 둘째로, 비선형 정상문제에 두 범함수를 응용하여 고립파(solitary wave)의 파고 및 속도의 근사 관계식을 유도하였다. 이 근사 관계식들의 수치계산결과는 Longuet-Higgins & Fenton(1974)의 해석해와 비교하여 좋은 결과를 보여준다. 세째로, 전 유체영역을 유한요소로 분할하여 유한요소법에 의한 고립파의 파형 및 특성치들의 계산결과들을 보여준다. 특히 본 논문에서는 분할한 유한요소의 수에 따른 수치해의 수렴을 보여준다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제4권1호
• /
• pp.121-131
• /
• 2000

Melt flow, heat and mass transfer of oxygen have been analyzed numerically in the process of Czochralski single crystal growth of silicon under the influence of misaligned cusp magnetic fields. Since the silicon melt in a crucible for crystal growth is of high temperature and of highly electrical-conducting, experimentation method has difficulty in analyzing the behavior of the melt flow. A set of simultaneous nonlinear equations including Navier-Stokes and Maxwell equations has been used for the modelling of the melt flow which can be regarded as a liquid metal. Together with the melt flow which forms the Marangoni convection, a flow circulation is observed near the comer close both to the crucible wall and the free surface. The melt flow tends to follow the magnetic lines instead of traversing the lines. These flow characteristics helps the flow circulation exist. Mass transfer characteristics influenced by the melt flow has been analyzed and the oxygen absorption rate to the crystal has been calculated and turned out to be rather uniform than in the case of an aligned magnetic field.

• 한국해양공학회지
• /
• 제18권3호
• /
• pp.32-39
• /
• 2004

A design procedure for a ship with minimum total resistance has been developed using a numerical optimization method called SQP(sequential quadratic programming) to search for different optimal hull forms. The frictional resistance has been estimated using the ITTC 1957 model-ship correlation line formula, and the wave resistance has been evaluated using a potential-flow panel method that is based on Rankine sources with nonlinear free surface boundary conditions. The geometry of a hull surface has been modified using B-spline surface patches, during the whole optimization process. The numerical analyses have been carried out for the modified Wilgey hull at three different speeds (Fn=0.25, 0.316, 0.408), and the calculation results were compared.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제3권3호
• /
• pp.201-207
• /
• 2011

In this paper, we focused on computing the higher-harmonic components of the transmitted wave passing over a submerged circular cylinder to show that it is causing a horizontal negative drift force. As numerical models, a circular cylinder held fixed under free surface in deep water is adopted. As the submergence of a circular cylinder decreases and the incident wavelength becomes longer, the higher-harmonic components of the transmitted wave starts to increase. An increase of the higher-harmonic components of the transmitted wave makes the horizontal drift force be negative. It is also found that the higher-harmonic amplitudes averaged over the transmitted wave region become larger with the increase of wave steepness and wavelength as well as the decrease of submergence depth.