• Journal of Ocean Engineering and Technology
• /
• v.20 no.5 s.72
• /
• pp.15-22
• /
• 2006

A finite-volume method using multi-block grid technique has been developed in order to investigate the flow field around multi-body. Here the matching grids are used at the interface between blocks and the boundary conditions are imposed there to exchange physical values across the interface. For the test problems, the laminar viscous flaw around one-and two-cylinder are simulated and the results are compared with experiments and other numerical results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.12
• /
• pp.3345-3356
• /
• 1994

A general purpose program, TURBO-3D, for the analysis of 3-D fluid flow in complex geometry has been developed, which employs a standard $k-\varepsilon$ turbulence model and a general nonorthogonal grid system. For the purpose of verification of the program and testing the applicability, turbulent flows in an S-shaped diffuser and turbulent flows over an backward facing step are solved and compared with the earlier results. Comparison with the results by the STAR-CD program has been also made for the same flow configuration and grid structure. The agreements are excellent and hence the program has been verified. Since the present program is applicable only on limited flow phenomena and lacks the pre-and post processor, further improvements toward these directions are being made.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.2 s.146
• /
• pp.186-196
• /
• 2006

The turbulent wake flow of an axisymmetric body at incidence of $10.1^{\circ}$ is investigated by commericial CFD code, Fluent 6.2. Reynolds stress turbulence model with wall function is applied for the turbulent flow computation. For the grid generation, the Gridgen V15 is used. Numerical predictions are compared with experimental data for the validation. The computed results show goof agreements with the experimental measurements, implying that the CFD analysis is a useful and efficient tool for predicting turbulent flow characteristics of wake field of an axisymmetric body at incidence.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.4
• /
• pp.232-239
• /
• 2013

This paper provides the structure of a Reynolds Averaged Navier-Stokes(RANS) based simulation method and its validation results for the ship motion problem. The motion information of the hull computed from the equations of motion is considered in the momentum equations as the relative fluid motions with respect to a non-inertial coordinates system. A finite volume method is used to solve the governing equations, while the free surface is captured by using a two-phase level-set method and the realizable k-${\varepsilon}$ model is used for turbulence closure. For the validation of the present numerical approach, the numerical results of the resistance and motion tests for DTMB 5415 at two ship speeds are compared against available experimental data.

• The KSFM Journal of Fluid Machinery
• /
• v.19 no.5
• /
• pp.50-60
• /
• 2016

A design methodology of the modeled test facility to conserve an injection performance of a passive safety injection system is proposed. This safety injection system is composed of a core makeup tank and a safety injection tank. Individual tanks are connected with pressure balance line on the top side and injection line on the bottom side. It is important to conserve the scaled initial injection flow rate and total injection time since this system can be operated by small gravity head without any active pumps. Differential pressure distribution of the injection line induced by the gravity head is determined by the vertical length and elevation of each tank. However, the total injection time is adjustable by the flow resistance coefficient of the injection line. The scaling methodology for the tank and flow resistance coefficient is suggested. A key point of this test facility design is a scaling analysis for the flow resistance coefficient. The scaling analysis proposed on this paper is based on the volume scaling law with the same vertical length to the prototype and can be extended to a model with a reduced vertical length. A set of passive injection test were performed for the tanks with the same volume and the different length. The test results on the initial flow rate and total injection time showed the almost same injection characteristics and they were in good agreement with the design values.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.247-247
• /
• 2019

수력발전 사업에 있어 Desander 구조물은 주로 고산지대 수력발전댐의 Run-of-river 형식의 발전방식에서 유사로 인한 터빈의 손상을 방지하기 위한 목적으로 설치된다. Desander의 적정 규모는 터빈의 손상을 일으킬 수 있는 유사 입경에 대해 안정적으로 침전을 시킬 수 있는 폭/길이/깊이로 평가할 수 있으며 상대적으로 Desander의 규모가 크게 설계된 경우 초기 공사비 증가하고 반대로 규모가 작게 설계된 경우 터빈의 교체 주기 단축으로 인한 유지관리비가 증가된다. 현재까지 일반적인 Desander 구조물의 설계 방식은 제거 입경의 침전 속도, 유입유량 및 깊이를 변수로 사용하여 경험식(L. Sudry method, Guicciardis method 및 Rouse method)을 통해 규모를 결정해 왔다. 하지만, 3-D 전산유체해석을 통해 유속 흐름 분석으로 직 간접적 Desander 규모의 적정성을 평가할 수 있는 현 시점에서 경험식으로부터 도출된 결과의 신뢰성과 객관성을 검증할 필요가 있다고 판단된다. 본 연구에서는 노르웨이 NSTU에서 개발한 유사의 이송 및 확산해석 기능이 내장된 범용 소프트웨어인 SSIIM을 이용하였다. SSIIM(Simulation of Sediment movements In water Intakes with Multiblock)은 개수로 흐름 상태에서 유사 이동 및 하상 변동을 분석할 수 있도록 개발된 3-D 해석 프로그램이다. SSIIM은 수치해석 방법으로 유한체적법(Finite Volume Method)를 채택하였으며 Navier-Stokes equations을 통해 유체의 흐름을 해석한다. 입력 자료는 유입 유량($m^3/sec$), 유입 유사량(kg/sec), 유출부 수위 및 해당 Desander Structure grid 자료가 사용되며 해석 결과로 Desander 내 grid 별 유속, 수위, 유사 농도 변화 등을 제공한다. 본 연구에서는 SSIIM을 이용하여 제거 목표 유사 입경의 차집 효율(Trap efficiency)로 Desander의 적정 규모를 평가 할 수 있는 설계법을 제안하며 설계 단계에서 결정되는 최소 제거유사 입자와 차집 효율에 의한 Desander의 적정 규모 평가 분석을 파키스탄 A 프로젝트를 대상으로 수행하였다. 연구 성과로 (1)SSIIM을 통해 해석된 차집 효율을 기초로 Desander의 적정 규모를 계획할 경우 경험적 방식에 비해 설계의 객관성과 신뢰성을 제고할 수 있으며 (2)3-D 수치해석을 통해 grid 별 유사농도를 확인 할 수 있어 Desander 형상과 규모에 대한 평가가 가능하다.

• Journal of Biomedical Engineering Research
• /
• v.21 no.4
• /
• pp.373-384
• /
• 2000

The goal of this paper is to develop a computational method to predict cancellous bone density distributions based upon continuum levels of volumetric strain. Volumetric strain is defined as the summation of normal strains, excluding shear strains, within an elastic range of loadings. Volumetric strain at a particular location in a cancellous structure changes with changes of the boundary conditions (prescribed displacements, tractions, and pressure). This change in the volumetric strain is postulated to predict the adaptive change in the bone apparent density. This bone remodeling theory based on volumetric strain is then used with the finite element method to compute the apparent density distribution for cancellous bone in both lumbar spine and proximal femur using an iterative algorithm, considering the dead zone of strain stimuli. The apparent density distribution of cancellous bone predicted by this method has the same pattern as experimental data reported in the literature (Wolff 1892, Keller et al. 1989, Cody et al. 1992). The resulting bone apparent density distributions predict Young's modulus and strength distributions throughout cancellous bone in agreement with the literature (Keller et al. 1989, Carter and Hayes 1977). The method was convergent and sensitive to changes in boundary conditions. Therefore, the computational algorithm of the present study appears to be a useful approach to predict the apparent density distribution of cancellous bone (i.e. a numerical approximation for Wolff's Law)

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.6
• /
• pp.429-434
• /
• 2017

A molecular-continuum hybrid method was developed to simulate microscale and nanoscale fluids where continuum fluidics cannot be used to predict Couette flow. Molecular dynamics simulation is used near the solid surface where the flow cannot be predicted by continuum fluidics, and Navier-Stokes equations are used in the other regions. Numerical simulation of Couette flow was performed using the hybrid method to investigate the effect of solid-liquid interaction and surface roughness in a nanochannel. It was found that the solid-liquid interaction and surface roughness influence the boundary condition. When the surface energy is low, slippage occurs near the solid surface, and the magnitude of slippage decreases with increase in surface energy. When the surface energy is high, a locking boundary condition is formed. The roughness disturbs slippage near the solid surface and promotes the locking boundary condition.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.3
• /
• pp.61-67
• /
• 2017

In this study, Burning Area Analysis Program (BAAP) was developed by using VOF method to estimate the burning area of 3D shaped grain. The parametric study of mesh size, burning rate and time interval for numerical calculation was conducted. The result of BAAP is compared with the one from commercial 3D modeling software. Also the internal ballistic analysis was performed using the result of BAAP. In order to estimate the burning area and internal pressure with time, Chemical Equilibrium Analysis (CEA) was conducted with a composition of reduced smoke propellant. As a result, the web-averaged pressure was 5.34 MPa which is similar to the published research result.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.2
• /
• pp.52-58
• /
• 2018

Hysteresis phenomena are often encountered in a wide variety of fluid flow systems used in industrial and engineering applications. Hence, in recent years, a significant amount of research been focusing on clarifying the physics of the flow hysteresis appearing during the transient change of the pressure ratios and influencing the performance of the supersonic wind tunnel. However, investigations on the hysteresis phenomenon, particularly when it occurs inside the supersonic wind tunnel, are rare. In this study, numerical simulations were carried out to investigate the hysteresis phenomena of the shock waves encountered in a supersonic wind tunnel. The unsteady and compressible flow was analyzed with an axisymmetric model, and the N-S equations were solved by using a fully implicit finite volume scheme. The optimal pressure ratio was determined from the hysteresis curves, and the results can be utilized to operate the wind tunnel efficiently.