• 대한기계학회논문집B
• /
• 제36권7호
• /
• pp.721-730
• /
• 2012

본 연구에서는 벽면부착에 의해 야기되는 다상유동에 대한 수치적 연구를 제시한다. 먼저 다상유동 해석을 위해 표면장력에 대한 CSF(Continuum Surface Force) 모델 및 벽면부착 경계조건 모델을 비정렬격자계에 적합하도록 수치해석방법을 정립시키고, Myong(2009)이 개발한 비정렬격자계와 VOF 방법으로 체적포착법(volume capturing method)을 사용한 수치해석방법(코드)에 삽입하였다. 또한 본 수치해석방법을 사용하여 중력을 포함하여 어떤 외력도 존재하지 않고 오직 벽면부착에 의해 야기되는 유동현상인 원통형 탱크의 바닥에 위치한 얕은 물풀(water pool)에 대해 물이 벽면을 적시는 경우와 적시지 않는 경우에 대해 수치해석 하였다. 연구결과, 본 수치해석방법은 벽면부착에 의해 야기되는 다상유동 문제에 대한 유용성이 입증되었다.

• 한국전산유체공학회지
• /
• 제19권4호
• /
• pp.20-28
• /
• 2014

An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 추계 학술대회논문집
• /
• pp.51-56
• /
• 2006

Several high resolution schemes such as OSHER, MUSCL, SMART, GAMMA, WACEB and CUBISTA are applied to two typical test cases of a translation test and a collapsing water column problem for the accurate capturing of fluid interfaces. It is accomplished by implementing the high resolution schemes in the in-house CFD code(PowerCFD) for computing 3-D flow with an unstructured cell-centered method, which is based on the finite-volume technique and fully conservative. The calculated results are found to show that SMART scheme gives the best performance with respect to accuracy and robustness.

• Nuclear Engineering and Technology
• /
• 제41권8호
• /
• pp.1045-1052
• /
• 2009

A calculation model is developed to predict the transient free surface flow on the containment floor following a loss-of-coolant accident (LOCA) of pressurized water reactors (PWR) for the use of debris transport evaluation. The model solves the two-dimensional Shallow Water Equation (SWE) using a finite volume method (FVM) with unstructured triangular meshes. The numerical scheme is based on a fully explicit predictor-corrector method to achieve a fast-running capability and numerical accuracy. The Harten-Lax-van Leer (HLL) scheme is used to reserve a shock-capturing capability in determining the convective flux term at the cell interface where the dry-to-wet changing proceeds. An experiment simulating a sudden break of a water reservoir with L-shape open channel is calculated for validation of the present model. It is shown that the present model agrees well with the experiment data, thus it can be justified for the free surface flow with accuracy. From the calculation of flow field over the simplified containment floor of APR1400, the important phenomena of free surface flow including propagations and interactions of waves generated by local water level distribution and reflection with a solid wall are found and the transient flow rates entering the Holdup Volume Tank (HVT) are obtained within a practical computational resource.

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

The computation of moving interface by the level set method typically requires reinitializations of level set function. An inaccurate estimation of level set function ${\phi}$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, accurate and robust reinitialization process is essential to the free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates directly level set function ${\phi}$ using a normal vector in the interface without solving the re-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1splitting FEM are adopted to discretize advection equation of the level set function and the Navier-Stokes equation, respectively. Advection equation of free surface and re-initialization process are validated with benchmark problems, i.e., a broken dam flow and time-reversed single vortex flow. The simulation results are in good agreement with the existing results.

• 대한조선학회논문집
• /
• 제47권3호
• /
• pp.279-290
• /
• 2010

In this paper, a VOF method called RHRIC (refined high resolution intertace capturing) is introduced for solving the motion of the free surface and applied to the simulation of the advection of rigid interiaces of different shapes and a 20 dam-break problem, which are typical benchmark test cases. The numerical results for the interface advection cases are compared to the analytic solutions, while the available experimental data and other numerical results of various free surface methods for the dam-break problem are provided for the validation of the proposed VOF method. The same simulations were also carried out using the original HRIC scheme and a modified HRIC scheme called MHRIC for comparison. Although the RHRIC uses a simple order scheme, a basis of the original HRIC scheme, lower than the third-order ULTIMATE-QUICKEST used by the MHRIC, it provides an improved accuracy over the two previous HRIC methods.

• 한국해양공학회지
• /
• 제14권1호
• /
• pp.11-16
• /
• 2000

수중날개를 대상으로 임의로 변형하는 이층류의 수직계산을 수행하였다. 계산은 비구조격자와 접면포착법을 이용하여 자유표면의 형상을 결정하였다 본 연구에 의하면 낮은 물속 깊이와 높은 프루드수에서 쇄파현상이 발생하였다 수중날개에 대하여 쇄파형상의 효과를 연구하였고 실험값과 계산값을 상호비교함으로서 격자의 민감성을 수직적으로 확인하였다.

• Steel and Composite Structures
• /
• 제37권5호
• /
• pp.551-569
• /
• 2020

Vibration of vertically aligned-monolayered-nonuniform nanorods consist of functionally graded materials with elastic supports has not been investigated yet. To fill this gap, the problem is examined using the elasticity theories of Eringen and Gurtin-Murdoch. The geometrical and mechanical properties of the surface layer and the bulk are allowed to vary arbitrarily across the length. The nonlocal-surface energy-based governing equations are established using differential-type and integro-type formulations, and solved by employing the Galerkin method by exploiting admissible modes approach and element-free Galerkin (EFG). Through various comparison studies, the effectiveness of the EFG in capturing both nonlocal-differential/integro-based frequencies is proved. A constructive parametric study is also conducted, and the roles of nanorods' diameter, length, stiffness of both inter-rod's elastic layer and elastic supports, power-law index of both constituent materials and geometry, nonlocal and surface effects on the dominant frequencies are revealed.

• 대한조선학회논문집
• /
• 제46권3호
• /
• pp.219-231
• /
• 2009

Numerical simulations have been carried out for a 138K LNG Carrier (KLNG) model ship with free surface, using WAVIS 2.0 and Fluent 6.3.26 with various $y^+$ values and different grid densities. Level-set method for free surface capturing was adopted in WAVIS, while VOF has been used in Fluent. The calculated results were compared with the experiment data. Resistance coefficient, wave pattern, wave profile along the hull surface, axial velocity contours and transverse vectors had been analyzed. When the first $y^+$ value was fixed at 60, the simulation results from both WAVIS and Fluent were improved as the number of grids increased. The convergence time of WAVIS was much shorter than that of Fluent. Furthermore, WAVIS predicted the velocity field and the wave profile along the hull surface better than Fluent. However, Fluent gave better wave patterns.

• 한국소음진동공학회논문집
• /
• 제15권4호
• /
• pp.395-403
• /
• 2005

The dynamic characteristics of delaminated smart composite laminates are studied using animproved layerwise laminate theory. The theory is capable of capturing interlaminar shear stresses that are critical to delamination. The presence of discrete delamination is modeled through the use of Heaviside unit step functions. Stress free boundary conditions are enforced at all free surfaces. Continuity in displacement field and transverse shear stresses are enforced at each ply level. In modeling piezoelectric composite plates, a coupled piezoelectric-mechanical formulation is used in the development of the constitutive equations. Numerical analysis is conducted to investigate the effect of nonlinearity in the transient vibration of bimodular behavior caused by the contact impact of delaminated interfaces. Composite plates with delamination, subject to external loads and voltage history from surface bonded sensors, are investigated and the results are compared with corresponding experimental results and plates without delamination.