• 대한기계학회논문집B
• /
• 제23권8호
• /
• pp.945-956
• /
• 1999

This work reports a set of approximate analytical solutions describing the initial transient process of close-contact melting between a rectangular parallelepiped solid and a flat plate on which either constant temperature or constant heat flux is imposed. Not only relative motion of the solid block tangential to the heating plate, but also the density difference between the solid and liquid phase is incorporated in the model. The thin film approximation reduces the force balance between the solid weight and liquid pressure, and the energy balance at the melting front into a simultaneous ordinary differential equation system. The normalized model equations admit compactly expressed analytical solutions which include the already approved two-dimensional solutions as a subset. In particular, the normalized liquid film thickness is independent of all pertinent parameters, thereby facilitating to define the transition period of close-contact melting. A unique behavior of the solid descending velocity due to the density difference is also resolved by the present solution. A new geometric function which alone represents the three-dimensional effect is introduced, and its properties are clarified. One of the representative results is that heat transfer is at least enhanced at the expense of the increase in friction as the cross-sectional shape deviates from the square under the same contact area.

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

Two-dimensional Navier-Stokes codes are used to simulate the shock tunnel in Konkuk university. In order to design experiments in impulse facilities properly and to interpret data from such facilities, it is necessary to understand how the flow approaches steady state. This is done by determining the transient flow field and flow establishment time around a given model. This will be accomplished by developing appropriate CFD codes which solve the Navier-Stokes equations, and simulating the starting process and resulting unsteady viscous flow phenomena. The starting process in a shock tunnel consists of multiple shock interactions and contact discontinuities, which are difficult to solve with the classical shock capturing schemes. A recently developed high resolution scheme is adapted for resolving the unsteady phenomena of those multiple shock interactions and contact surfaces during the starting process. The bifurcation phenomenon due to the interactions of the reflected shock from the end of the shock tube with the boundary layer generated by the incident shock becomes of particular interest. By comparing with the experiment results, the accuracy of the numerical analysis is validated and it is demonstrated that the properties which can hardly be obtained through the experiment can be estimated.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제8권6호
• /
• pp.589-601
• /
• 2016

The hybrid grid was adopted and numerical prediction analysis of propeller unsteady bearing force considering free surface was performed for mode and full-scale KCS hull-propeller-rudder system by employing RANS method and VOF model. In order to obtain the propeller velocity under self-propulsion point, firstly, the numerical simulation for self-propulsion test of full-scale ship is carried out. The results show that the scale effect of velocity at self-propulsion point and wake fraction is obvious. Then, the transient two-phase flow calculations are performed for model and full-scale KCS hull-propeller-rudder systems. According to the monitoring data, it is found that the propeller unsteady bearing force is fluctuating periodically over time and full-scale propeller's time-average value is smaller than model-scale's. The frequency spectrum curves are also provided after fast Fourier transform. By analyzing the frequency spectrum data, it is easy to summarize that each component of the propeller bearing force have the same fluctuation frequency and the peak in BFP is maximum. What's more, each component of full-scale bearing force's fluctuation value is bigger than model-scale's except the bending moment coefficient about the Y-axis.

• Nuclear Engineering and Technology
• /
• 제18권2호
• /
• pp.107-116
• /
• 1986

새로이 개발된 망상구조 충전탑에서의 중수증류성능해석 및 설계를 위한 정상 및 비정상상태 모델을 수립하였다. 정상상태 모델은 MESH 방정식에 각단의 압력강하를 고려하여 설정하고, Equation Tearing방법으로 그 해를 구하여 중수중류탄내의 농도, 온도 및 압력구배를 얻었다. 비정상상태 거동해석을 위하여 Cohen의 이상단 모델을 수정한 평형단 전이모델을 세웠으며, 그 모델식의 해석적 해를 구함으로써 시간에 따른 중수 농축과정을 예측할 수 있게 되었다. 설정된 모델들은 전환류 중수증류탑에서의 실험결과와 매우 잘 일치하였으며 개발된 충전물임이 높은 이론단수에 낮은 압력강하, 저체류량을 갖는 우수한 중수농축용 충전물임이 확인되었다.

• Journal of Mechanical Science and Technology
• /
• 제20권4호
• /
• pp.533-544
• /
• 2006

The amplitude and phase lag of surface deformation were determined for a compliant coating under the action of turbulent pressure fluctuations. For this purpose, pressure fluctuations were measured experimentally. The amplitude and duration of coherent wave train of pressure fluctuations were investigated using digital filtration. The transient response was calculated for stabilization of forced oscillations of the coating in approximation of local deformation. The response of coating was analyzed with considerations of its inertial properties and limited duration of coherent harmonics action of pressure fluctuations. It is shown that a compliant coating interacts not with the whole spectrum of pressure fluctuations, but only with a frequency range near the first resonance. According to the analysis, with increasing elasticity modulus of the coating material E, deformation amplitude decreases as 1/E, and dimensionless velocity of the coating surface decreases as $1/\sqrt{E}$. For sufficiently hard coatings, deformation amplitude becomes smaller than the thickness of viscous sublayer, while surface velocity remains comparable to vertical velocity fluctuations of the flow.

• 대한기계학회논문집B
• /
• 제28권4호
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• Korea-Australia Rheology Journal
• /
• 제18권2호
• /
• pp.109-118
• /
• 2006

A finite-element method was employed to analyze axisymmetric unsteady motion of a deformable bubble near the wall. In the present study a deformable bubble in a Newtonian medium near the wall was considered. In solving the governing equations a structured mesh generator was used to describe the collapse of highly deformed bubbles with the Arbitrary Lagrangian Eulerian (ALE) method being employed in order to capture the transient bubble boundary effectively. In order to check the accuracy of the present FE analysis we compared the results of our FE solutions with the result of the collapse of spherical bubbles in a large body of fluid in which solutions can be obtained using a 1D FE analysis. It has been found that 1D and 2D bubble deformations are in good agreement for spherically symmetric problems confirming the validity of the numerical code. Non-spherically symmetric problems were also solved for the collapse of bubble located near a plane solid wall. We have shown that a microjet develops at the bubble boundary away from the wall as already observed experimentally. We have discussed the effect of Reynolds number and distance of the bubble center from the wall on the transient collapse pattern of bubble.

• 대한기계학회논문집B
• /
• 제28권11호
• /
• pp.1348-1358
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation fur flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among seven filling patterns at each tetrahedral control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. The collapse of a water dam and the filling of a fluidity spiral have been analyzed. The numerical results have been in good agreement with the experimental results and the efficiency of the adaptive grid refinement and filling pattern techniques have been verified.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제11권1호
• /
• pp.154-164
• /
• 2019

A rotor dynamic analysis is mandatory for stability and design optimization of submerged propellers and turbines. An accurate simulation requires a proper consideration of fluid-induced reaction forces. This paper presents a bi-directional coupling of a bond graph method solver and an unsteady vortex lattice method solver where the former is used to model the rotor dynamics of the power train and the latter is used to predict transient hydrodynamic forces. Due to solver coupling, determination of hydrodynamic coefficients is obsolete and added mass effects are considered automatically. Additionally, power grid and structural faults like grid fluctuations, eccentricity or failure could be investigated using the same model. In this research work a fast, time resolved dynamic simulation of the complete power train is conducted. As an example, the rotor dynamics of a tidal stream turbine is investigated under two inflow conditions: I - shear flow, II - shear flow + water waves.

• 한국전산구조공학회논문집
• /
• 제19권3호
• /
• pp.295-302
• /
• 2006

본 논문은 자이로콥터의 유한요소 모델링과 로터의 회전에 의한 동하중을 고려한 전산진동해석을 수행하였다. 이를 위해 자이로콥터의 최종 조립된 3차원 CATIA 모델을 구축하였으며, 3차원 데이터를 바탕으로 비구조 질량을 포함한 구조진동해석을 위한 유한요소모델을 생성하였다. 실용적인 전산구조동역학 해석을 위해 상용 유한요소 해석프로그램인 MSC/NASTRAN과 자체 개발한 프로그램을 병용하여 사용하였다. 비행 중 회전하는 로터에 의해 발생하는 동하중은 상용 CFD 프로그램인 FLUENT를 이용하였다. 유체해석과 구조진동해석의 결합을 위해 자체적으로 통합 연계 시스템을 구축하였다. 3차원 구조의 효율적인 진동특성을 고찰하기 위해 모달영역에서 천이응답해석과 주파수응답해석을 수행하였다. 실제 자이로콥터의 연료조건과 비행조건을 고려하였으며, 전산해석을 통하여 고유진동, 주파수 응답 및 천이응답 특성을 고찰하였다.