• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1149-1164
• /
• 1997

In this study, the applicability of the RNG k-.epsilon. model to the analysis of the complex flows is studied. The governing equations based on a non-orthogonal coordinate formulation with Cartesian velocity components are used and discretized by the finite volume method with non-staggered variable arrangements. The predicted results using the RNG k-.epsilon. model of three complex flows, i.e., the flow over a backward-facing step and a blunt flat plate, the flow around a 2D model car are compared to these from the standard k-.epsilon. model and experimental data. That of the unsteady axisymmetric turbulent flow within a cylinder of reciprocating model engine including port/valve assembly and the spray characteristics within a chamber of direct injection model engine are compared to these from the standard k-.epsilon. model and experimental data. The results of reattachment length, separated eddy size, average surface pressure distribution using the RNG k-.epsilon. model show more reasonable trends comparing with the experimental data than those using the modified k-.epsilon. model. Although the predicted rms velocity using the modified k-.epsilon. model is lower considerably than the experimental data in incylinder flow with poppet valve, predicted axial and radial velocity distributions at the valve exit and in-cylinder region show good agreements with the experimental data. The spray tip penetration predicted using the RNG k-.epsilon. model is more close to the experimental data than that using the modified k-.epsilon. model. The application of the RNG k-.epsilon. model seems to have some potential for the simulations of the unsteady turbulent flow within a port/valve-cylinder assembly and the spray characteristics over the modified k-.epsilon. model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.2
• /
• pp.187-196
• /
• 2013

A numerical model was developed to predict the density wave oscillation (DWO) in the furnace wall tubes of a fossil-fired once-through boiler. The transient flow fields in the tubes were obtained using a 1D finite volume method in the time domain. A header model was also implemented to simulate the parallel tube connection of the wall tubes. The inlet and outlet mass flow variation in one of the parallel tubes was examined after a heat perturbation to find the DWO. After successful verification with experimental results reported in literature, the developed model was applied to the wall tubes of a 700-MW boiler furnace. In contrast to the simulation of Takitani's experiment, in which the unstable power thresholds tended to rise in the reduced bypass channel flow, no remarkable changes were observed in the power thresholds in the parallel channel modeling of the wall tubes of the boiler furnace.

• Journal of Biomedical Engineering Research
• /
• v.18 no.2
• /
• pp.187-196
• /
• 1997

The three-dimensional flow analysis using the finite volume method is presented to compare the steady flow characteristics of blood with those of blood substitutes such as water and aqueous polymer solution in an idealized double branching model. The model is used to simlllate the region of the abdominal aorta near the celiac and superior mesenteric branches. Apparent viscosities of blood and the aqueous Separan solution are represented as a function of shear rate by the Carreau model, Water and aqueoiu Separan AP-273 500wppm solution are frequently used as blood substitutes in vitro experiments. Water is a typical Newtonian fluid and blood and Separan solution are non-Newtonian fluids. Flow phenomena such as velocity distribution, pressure variation and wall shear stress distribution of water, blood and polymer solution are quite different due to differences of the rheological characteristics of fluids. Flow phenomena of polymer solution are qualitatively similar to those of blood but the phenomena of water are quite different from those of blood and polymer solution. It is recommended that a lion-Newtonian fluid which exhibits very similar rheological behavior to blood be used in vitro experiments. A non-Newtonian fluid whose rheological characteristics are very similar to those of blood should be used to obtain the meaninylll hemodynamic data for blood flow in vitro experiment and by numerical analysis

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.363-374
• /
• 2009

The 2D flow around 13 similar stay-vane profiles with different trailing edge geometries is investigated to determinate the main characteristics of the excitation forces for each one of them and their respective dynamic behaviors when modeled as a free-oscillating system. The main goal is avoid problems with cracks of hydraulic turbines components. A stay vane profile with a history of cracks was selected as the basis for this work. The commercial finite-volume code $FLUENT^{(R)}$ was employed in the simulations of the stationary profiles and, then, modified to take into account the transversal motion of elastically mounted profiles with equivalent structural stiffness and damping. The k-$\omega$ SST turbulence model is employed in all simulations and a deforming mesh technique used for models with profile motion. The static-model simulations were carried out for each one of the 13 geometries using a constant far field flow velocity value in order to determine the lift force oscillating frequency and amplitude as a function of the geometry. The free-oscillating stay-vane simulations were run with a low mass-damping parameter ($m^*{\xi}=0.0072$) and a single mean flow velocity value (5m/s). The structural bending stiffness of the stay-vane is defined by the Reduced Velocity parameter (Vr). The dynamic analyses were divided into two sets. The first set of simulations was carried out only for one profile with $2{\leq}Vr{\leq}12$. The second set of simulations focused on determining the behavior of each one of the 13 profiles in resonance.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.86-94
• /
• 2010

A thermal-hydraulic code, named CUPID, has been developed for the analysis of transient two-phase flows in nuclear reactor components. A two-fluid three-field model was used for steam-water two-phase flows. To obtain numerical solutions, the finite volume method was applied over unstructured cell-centered meshes. In steam-water two-phase flows, a phase change, i.e., evaporation or condensation, results in a great change in the flow field because of substantial density difference between liquid and vapor phases. Thus, two-phase flows are very sensitive to the local pressure distribution that determines the phase change. This in turn puts emphasis on the accurate evaluation of local pressure gradient. This paper presents a new reconstruction method to evaluate the pressure gradient at cell centers on unstructured meshes. The results of the new scheme for a simple test function, a gravity-driven cavity, and a wall boiling two-phase flow are compared with those of the previous schemes in the CUPID code.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.1-5
• /
• 2010

1982년 10월 21일 스페인에서 발생한 Tous 댐 붕괴 시 발생한 유출로 인하여 댐에서 5km 하류에 위치한 도시지역을 홍수파가 관통하면서 시간에 따라 변화되는 홍수파의 특성을 수치모형을 이용하여 분석하였다. 분석에 적용된 수치모형은 비구조적 2차원 유한체적법과 불연속 흐름을 모의하기 위해 시간과 공간상으로 1차 정확도를 가진 HLLC 기법을 적용한 모형이다. 본 연구에서는 도시지역의 침수현상을 두 가지 Case에 대해 구성하였으며, 모의결과를 침수흔적 자료 및 기존의 수치모의 결과와 비교하였다. 첫 번째 Case는 Mulet and Alcrudo(2004)에 의한 수치모의결과와의 비교를 통해 모형의 적용성 검증을 하였으며, 두 번째 Case는 도시 지역에서의 건물 투수여부에 따라 변화되는 홍수파의 흐름특성 및 지점별 최대 침수심을 모의하였을 경우이다. Case 1에서 모형의 검증결과 수치모형을 통해 계산된 지점별 수심변화는 비교적 잘 일치하는 결과를 얻었으며, Case 2의 경우 지점별 수심은 Case 1에 비해 다소 과소산정 되는 경향을 나타내었다.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.6
• /
• pp.507-514
• /
• 2019

Computational studies of accelerating flow around 2D Circular Cylinder was performed to investigate characteristics of wake field and drag forces. Previous studies had revealed that drag on the cylindrical body in accelerating flow is much greater than that in the flow with constant velocity; however, the underlying physics on the drag increase has not been clearly investigated. In order to investigate the drag increase and its relationship with wake development, this study employed a finite-volume based CFD code, Fluent 13.0 with k-ω SST model for turbulence effects. Inflows are modeled with varied accelerations from 0.4905 to 9.81m/s². The drag computed in the present study is in good agreement with previous studies, and clearly shows the increase compared to the drag on the body in the flow with constant velocity. The results also show that drag crisis observed at high Reynolds number in the case of the flow with constant velocity is also found in the case of accelerating flow. The analysis for wake and recirculation length shows that conventional vortex shedding does not occur even at high Reynolds number and the drag increase is larger at higher acceleration.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.39 no.3
• /
• pp.283-289
• /
• 2015

In this study, the lightweight design of a brake bracket for a composite bogie was studied by considering two brake bracket models with thicknesses of 12t and 9t, respectively. For achieving this goal, finite element analysis and topology optimization were conducted. Firstly, the largest cross-sectional areas of the vertical and horizontal plates of the brake bracket were selected as the design variables. As the constraint, the Z-axis displacement of the brake bracket was increased by 2.5 units from the initial displacement value. The minimum volume fraction of the design regions was chosen as the objective function. The full model comprised a composite bogie frame and brackets attached together. However, to reduce the analysis time, 1D beam elements were used instead of the composite bogie frame by ensuring its equivalence with the full model. The result revealed that the weights of the 12t and 9t models of the brake bracket were reduced to 60 kg and 31 kg, respectively.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.1
• /
• pp.49-56
• /
• 1997

Heat-transfer enhancement is seeked through modifications of fin surface. Real life plate-fin heat exchangers have complex three-dimensional geometries. Fins can have arrays of dimples and are attached to rows of penetrating tubes. To isolate the effect of surface modification, we model the real flow by a two-dimensional channel flow with a dimple on one side. The flow is analysed by solving the incompressible Navier-Stokes equation by a finite volume method on a generalized boundary-fitted coordinate. Results show a trapped vortex inside the dimple for all cases computed. Local maximum of Nusselt number occurs near the downstream end of the dimple, due to such a vortex. Location of the vortex does not change with respect to the wall temperature change, but moved downstream when Reynolds number increases. This, together with the results that in all cases vortex core is somewhat downstream of the dimple center, suggests that the mean flow above continuously feeds the kinetic energy to the recirculating flow. Heat transfer enhancement and pressure losses are studied through analysing the relevant dimensionless parameters like, Nusselt number and friction factor. In all cases computed, dimpled channel flow experiences less pressure loss than two-dimensional Poiseuille flow.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.86-92
• /
• 2014

In this paper, prediction of separation trajectory for Two-stage-To-Orbit space launch vehicle has been numerically simulated by using an aerodynamic database based on steady state analysis. Aerodynamic database were obtained for matrix of longitudinal and vertical positions. The steady flow simulations around the launch vehicle have been made by using a 3-D RANS flow solver based on unstructured meshes. For this purpose, a vertex-centered finite-volume method was adopted to discretize inviscid and viscous fluxes. Roe's finite difference splitting was utilized to discretize the inviscid fluxes, and the viscous fluxes were computed based on central differencing. To validate this flow solver, calculations were made for the wind-tunnel experiment model of the LGBB TSTO vehicle configuration on steady state conditions. Aerodynamic database was constructed by using flow simulations based on test matrix from the wind-tunnel experiment. ANN(Artificial Neural Network) was applied to construct interpolation function among aerodynamic variables. Separation trajectory for TSTO launch vehicle was predicted from 6-DOF equation of motion based on the interpolated function. The result of present separation trajectory calculation was compared with the trajectory using experimental database. The predicted results for the separation trajectory shows fair agreement with reference[4] solution.