• 한국전산유체공학회지
• /
• 제18권4호
• /
• pp.17-24
• /
• 2013

In this study, we consider the characteristics of channel flow in the presence of an infinite streamwise array of equispaced identical rotating circular cylinders. This flow configuration can be regarded as a model representing a micro channel or an internal heat exchanger with cylindrical vortex generators. A numerical parametric study has been carried out by varying Reynolds number based on the bulk mean velocity and the cylinder diameter, and the gap between the cylinders and the channel wall for some selected angular speeds. An immersed boundary method was employed to facilitate implementing the cylinders on a Cartesian grid system. No-slip condition is employed at all solid boundaries including the cylinders, and the flow is assumed to be periodic in the streamwise direction. The presence of the rotating circular cylinders arranged periodically in the streamwise direction causes a significant topological change of the flow, leading to increase of mean friction on the channel walls. More quantitative results as well as qualitative physical explanations are presented to justify the effectiveness of rotating cylinders to modify flow topology, which might be used to enhance heat transfer on the channel walls.

• 대한조선학회논문집
• /
• 제53권1호
• /
• pp.1-9
• /
• 2016

It is important to research and understand the physical phenomenon around a semi-submersible offshore structure on waves and currents because the wave run-up and load occurs owing to the waves and currents. In this study, the numerical simulations are performed about flow around a fixed semi-submersible offshore structure. The Modified Marker-density method is adopted in the present computation procedure, this method is one of the various methods to define the free-surface. The present computation results are compared with existing experimental and numerical simulation(VOF method) results. And, the computation results are relatively coincident with the existing results of model test and numerical simulation by VOF method.

• Nuclear Engineering and Technology
• /
• 제41권10호
• /
• pp.1347-1360
• /
• 2009

A multi-dimensional component for the thermal-hydraulic system analysis code, MARS, was developed for a more realistic three-dimensional analysis of nuclear systems. A three-dimensional and two-fluid model for a two-phase flow in Cartesian and cylindrical coordinates was employed. The governing equations and physical constitutive relationships were extended from those of a one-dimensional version. The numerical solution method adopted a semi-implicit and finite-difference method based on a staggered-grid mesh and a donor-cell scheme. The relevant length scale was very coarse compared to commercial computational fluid dynamics tools. Thus a simple Prandtl's mixing length turbulence model was applied to interpret the turbulent induced momentum and energy diffusivity. Non drag interfacial forces were not considered as in the general nuclear system codes. Several conceptual cases with analytic solutions were chosen and analyzed to assess the fundamental terms. RPI air-water and UPTF 7 tests were simulated and compared to the experimental data. The simulation results for the RPI air-water two-phase flow experiment showed good agreement with the measured void fraction. The simulation results for the UPTF downcomer test 7 were compared to the experiment data and the results from other multi-dimensional system codes for the ECC delivery flow.

• 설비공학논문집
• /
• 제8권3호
• /
• pp.375-385
• /
• 1996

Turbulent flow and heat transfer characteristics around staggered tube banks were studied using the 3-D Navier-Stokes equations and energy equation governing a steady incompressible flow, which were reformulated in a non-orthogonal coordinate system with cartesian velocity components and discretized by the finite volume method with a non-staggered variable arrangement. The predicted turbulent kinetic energy using RNG $k-{\varepsilon}$ model was lower than that of standard $k-{\varepsilon}$ model but showed same result for mean flow field quantities. The prediction of the skin friction coefficient using RNG $k-{\varepsilon}$ model showed better trend with experimental data than standard $k-{\varepsilon}$ model result. The inclined flow showed higher velocity and skin friction coefficient than transverse flow because of extra strain rate ($\frac{{\partial}w}{{\partial}y}$). Also, this was why the inclined flow showed higher local heat transfer coefficient than the transverse flow.

• 산업기술연구
• /
• 제25권A호
• /
• pp.191-201
• /
• 2005

A numerical analysis was done for a tunnel fire in a 1000m road tunnel. A cartesian coordinate was adopted to make a computational grid sytem which has 448,000 computational cells. A transient flow phenomena in the tunnel was simulated by the commercial code of PHEONICS from the ignition of fire to 600 seconds by the interval of 100 seconds. Total computational time of about 44 hours was required to get a convered solution in each time step. The purpose of this research is to analyze of the backlayering pheonomena and recirculation flow in a tunnel. The compuational results say that the backlayering does not happens near the fire of vehicle in this case because the vehicle fire is located at the outside of recirculation zone of flow ocuured near the jet fan. In this research, onset of backlayering pheonomena could be escaped if jet fan is set 95m in front of the the fire of vehicle.

• 대한조선학회논문집
• /
• 제51권4호
• /
• pp.265-273
• /
• 2014

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.

• 대한기계학회논문집B
• /
• 제21권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.

• 한국전산유체공학회지
• /
• 제18권2호
• /
• pp.85-92
• /
• 2013

In this study, we consider heat transfer enhancement in laminar channel flow by means of an infinite streamwise array of equispaced identical circular cylinders. This flow configuration can be regarded as a model representing a micro channel or an internal heat exchanger with cylindrical vortex generators. A numerical parametric study has been carried out by varying Reynolds number based on the bulk mean velocity and the cylinder diameter, and the gap between the cylinders and the channel wall. An immersed boundary method was employed to facilitate to implement the cylinders on a Cartesian grid system. No-slip condition is employed at all solid boundaries including the cylinders, and the flow is assumed to be periodic in the streamwise direction. Also, the Prandtl number is fixed as 0.7. For thermal boundary conditions on the solid surfaces, it is assumed that heat flux is constant on the channel walls, while the cylinder surfaces remain adiabatic. The presence of the circular cylinders arranged periodically in the streamwise direction causes a significant topological change of the flow, leading to heat transfer enhancement on the channel walls. The Nusselt number averaged on the channel wall is presented for the wide ranges of Reynolds number and the gap. A significant heat transfer enhancement is noticed when the gap is larger than 0.8, while the opposite is the case for smaller gaps. More quantitative results as well as qualitative physical explanations are presented to justify the effectiveness of varying the gap to enhance heat transfer from the channel walls.

• 대한기계학회논문집
• /
• 제18권1호
• /
• pp.139-149
• /
• 1994

A numerical simulation of unsteady axisymmetric turbulent flow was performed for a reciprocating engine including port/valve assembly. The governing equations based on a nonorthogonal coordinate formulation with Cartesian velocity components were used and discretised by the finite volume method with non-staggered variable arrangements. The modified $\kappa-\xi$. turbulence model which included the effect of compressibility was used. The results of twodimensional transient calculation for the axisymmetric configuration were compared with the experimental data. Although slightly low rms velocity was predicted compared to the experimental data, predicted velocity distributions at the valve exit and in-cylinder region showed good agreements with the experimental data. The flow at the valve exit was separated at the same valve lift position with the experimental data. Two vortices incylinder region were generated during the initial intake process. The clockwise main vortex became strong and moved upward to the top wall. The counter-clockwise second vortex became weak and stick to the upper left corner of the cylinder. After middle intake process, new vortex adjacent to upper cylinder wall appeared by the piston motion and therefore, the in-cylinder flow was formed into three vortices. The cylinder pressure just before bottom dead center of piston was higher than inlet pressure and then the reverse flow occured at the valve exit. The in-cylinder flow characteristics were strongly dependent on piston motion, but insensitive to valve motion.

• 한국해안·해양공학회논문집
• /
• 제23권3호
• /
• pp.248-257
• /
• 2011

본 연구에서는 규칙파동장에 있어서 인장계류된 원형부유체의 파랑응답해석에 2차원 Navier-Stokes solver에 기초한 새로운 수치파동수조모델을 제안하였다. 본 수치파동수조모델에서는 이동구조물과 유체와의 상호작용을 해석하기 위하여 직각좌표계에서 임의형상의 불투과경계를 갖는 구조물과 유체와의 연성해석이 가능한 IBM(Immersed Boundary Method)과 자유수면 추적을 위한 VOF(Volume of Fluid)법을 결합하였다. 부유체운동에 대한 수치결과를 기존의 FAVOR(Fractional Area Volume Obstacle Representation)법에 의한 계산결과 및 수리실험과 비교하였다. 게다가, 수치모델의 보다 자세한 검증을 위하여 원형부유체의 동요 및 자유수면변동에 관한 수리모형을 추가로 실시하였고, 제안된 수치모델의 범용성과 타당성을 검증하기 위하여 직사각형부유체에 대한 수치 및 수리실험도 병행하였다. 이로부터 추정된 수치계산결과는 실험결과를 잘 재현하고 있는 것으로 판단되었다.