• Journal of computational fluids engineering
• /
• v.18 no.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.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.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
• /
• v.41 no.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.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.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.

• Journal of Industrial Technology
• /
• v.25 no.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.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.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.

• 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.

• Journal of computational fluids engineering
• /
• v.18 no.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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.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.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.3
• /
• pp.248-257
• /
• 2011

In the present study, we proposed a new numerical wave tank model to analyze the vertical tension-leg circular floating bodies, using a 2-D Navier-Stokes solver. An IBM(Immersed Boundary Method) capable of handling interactions between waves and moving structures with complex geometry on a standard regular Cartesian grid system is coupled to the VOF(Volume of Fluid) method for tracking the free surface. Present numerical results for the motions of the floating body were compared with existing experimental data as well as numerical results based on FAVOR(Fractional Area Volume Obstacle Representation) algorithm. For detailed examinations of the present model, the additional hydraulic experiments for floating motions and free surface transformations were conducted. Further, the versatility of the proposed numerical model was verified via the numerical and physical experiments for the general rectangular floating bodies. Numerical results were compared with experiments and good agreement was archived.