• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.1-10
• /
• 2011

A new computational program, which is based on the CIP/CCUP(Constraint Interpolation Profile/CIP Combined Unified Procedure) method, has been developed to numerically analyse sloshing phenomena dealt as multiphase-flow problems. For the convection terms of Navier-Stokes equations, the RCIP(Rational function CIP) method was adopted and the THINC-WLIC(Tangent of Hyperbola for Interface Capturing-Weighted Line Interface Calculation) method was used to capture the air/water interface. To validate the present numerical method, two-dimensional dam-breaking and sloshing problems in a rectangular tank were solved by the developed method in a stationary Cartesian grid system. In the case of sloshing problems, simulations by using a improved MPS(Moving Particle Simulation) method, which is named as PNU-MPS(Pusan National University-MPS), were also carried out. The computational results are compared with those of experiments and most of the comparisons are reasonably good.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.4
• /
• pp.325-333
• /
• 2015

This study was devoted to elucidating the change in the flow characteristics of a laminar flow past a circular cylinder by modifying the cylinder shape with O-grooves. A numerical analysis was performed in a two-dimensional framework. The cylinder was represented using an immersed boundary method and marker particles on a Cartesian grid system. The number and locations of the O-grooves were the key parameters. An analysis of the flow pattern and flow-induced forces was performed at Re = 40 and 50. In addition, we calculated the critical Reynolds number depending on the number of O-grooves and their locations.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.164-169
• /
• 2009

Study on swimming of microorganisms like, sperm motility, cilia beating, bacterial flagellar propulsion has found immense significance in the field of biological fluiddynamics. Because of the complexity involved, it is challenging for the researchers to model such problems. Immersed boundary method has proved its efficacy in the field of biological fluiddynamics, The present work aims at performing a numerical study on the microorganism locomotion using the immersed boundary method proposed by Peskin[1]. A two-dimensional model of the microorganism is modeled as thin elastic filament described as a sine wave. The neutrally buoyant organism undergoing deformations is immersed in a viscous and incompressible fluid. The fluid quantities are described using Eulerian coordinates and the immersed body is represented by Lagrangian coordinates. The Eulerian and Lagrangian variables are connected by the Dirac delta function. The Navier-Stokes equations governing the fluid flow are solved using the fractional step method on a staggered Cartesian grid system. The developed numerical code in FORTRAN will be validated by comparing the numerical results with the available results.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.28 no.7
• /
• pp.293-298
• /
• 2016

A three-dimensional (3D) numerical model of the vertical ground-coupled heat exchanger is useful for analyzing the modern ground source heat pump system. Furthermore, a detailed description of the inner side of the exchanger allows to account for the effects of the thermal capacity. Thus, both methods are included in the proposed numerical model. For the ground portion, a FDM (Finite Difference Method) scheme has been applied using the Cartesian coordinate system. Cylindrical grids are applied for the borehole portion, and the U-tube configuration is adjusted at the grid, keeping the area and distance unchanged. Two sub-models are numerically coupled at each time-step using an iterative method for convergence. The model is validated by a reference 3D model under a continuous heat injection case. The results from a periodic heat injection input show that the proposed thermal capacity model reacts more slowly to the changes, resulting in lower borehole wall temperatures, when compared with a thermal resistance model. This implies that thermal capacity effects may be important factors for system controls.

• KIEAE Journal
• /
• v.4 no.3
• /
• pp.179-186
• /
• 2004

The evaluation of the indoor environment of the Assembly Hall in the University, which is designed to be a large space, requires efficient design of its heating system that takes into consideration natural convection and the characteristics of the occupant's spaces. Indoor thermal environment was measured in the field and simulated with CFD code. The estimations of temperature distribution and indoor airflow distribution must be carried out simultaneously, as the thermal stratification is induced by natural convection flows. In order to simulate the even distribution of factors affecting the indoor environment, including temperature and airflow, Phoenics is used. The turbulent flow model adopted is the RNG k- model. The inlets and outlets of the air-conditioning systems, material and thermal properties, and the size of the test room ($35m{\times}18m{\times}10m$) are used for the simulation. Since the Assembly Hall is symmetric, half of the space is simulated. A Cartesian grid is used for calculation and the number of grids are respectively $60{\times}45{\times}35$. The results of the computer simulation during winter conditions are compared with the measurements at the typical points in the assembly hall with the heating system. After evaluating the results of the computer simulations, the methods of the heating system and layout are suggested.

• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.1-9
• /
• 2015

A developed code is applied to simulate relative motion of floating bodies in a side-by-side arrangement, including effects of a fender and a hawser. The developed code is based on the flux-difference splitting scheme for immiscible incompressible fluids and the hybrid Cartesian/immersed boundary method. To validate the developed code for free surface flows around deforming boundaries, the water wave generation is simulated, which is caused by bed movement. The computed wave profile and time histories of wave elevation are compared with other experimental and computational results. The effects of a fender and a hawser are modeled by asymmetric force acting on the floating bodies according to a relative displacement with the bounds, in which the fender and the hawser exert no force on the bodies. It has been observed that the floating body can be accelerated by a gap flow due to a phase difference caused by the free surface. Grid independency is established for the computed time history of the body velocity, based on three different size grids.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.44-53
• /
• 2012

Immersed boundary method(IBM) is a numerical scheme proposed to simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies, the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. The weight coefficients of the bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies. For an analysis of moving boundary, we smulated an oscillating circular cylinder with Re=100 and KC(Keulegan-Carpenter) number of 5. The predicted flow fields were compared with experimental data and they also showed good agreements.

• Journal of computational fluids engineering
• /
• v.22 no.1
• /
• pp.81-87
• /
• 2017

The characteristics of flow past a square cylinder submerged under the free surface have been numerically studied. An immersed boundary method was adopted for implementation of the cylinder cross-section in a Cartesian grid system. Also, a level-set method was used to capture the interface of the two fluids. The case for Reynolds number 150 was examined. At the specific Reynolds number, by varying the gap ratio(0.25, 0.40, 0.55, 0.70, 1.00, 1.50, 2.50, 5.00) the effects of the free surface on the force coefficients and Strouhal number of vortex shedding were identified. The presence of the free surface very close to the cylinder significantly affects the shedding pattern, resulting in considerable deviation of the force coefficients and Strouhal number from those of the single-phase flow. In addition, the influence of Froude number was considered in this study. By increasing Froude number(0.2-0.4), flow topology change was identified at the specific gap ratios(0.40, 0.70, 1.50, 5.00).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.8
• /
• pp.1046-1055
• /
• 2000

The flow and thermal phenomena in flat tubes of radiator are analyzed numerically. To predict the characteristics of heat transfer and pressure drop, the flow analysis program for three-dimensional complex geometry is developed, which adopted an non-staggered grid system and Cartesian velocities as dependent variables of the momentum equations. Using the developed program, the effect of tube specifications on the heat transfer characteristics is investigated for various flat tubes. From this study, the following results are obtained; (1) For the same hydraulic diameter($D_h{\doteq}5.2$mm), the Nusselt numbers of three basic modeis(D, J, and H-model) are 8.71, 8.92, and 10.58, respectively, and the pressure drops of D-, J-, and H-model are predicted as $-3.08{\times}10^{-2}\;Pa,\;-3.12{\times}10^{-2}\;Pa,\;and\; -3.98{\times}10^{-2}$ Pa, (2) In case of the same flat tube specification, the fins must be brazed at upper tube surface because the heat is more vividly transferred. Therefore, it is found that the H- model is the most effective tube as a heat exchanger and these results are used as a fundamental data for the design of tube.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.5
• /
• pp.349-357
• /
• 2009

Effect of computational domain size in simulating of periodic obstacle flow has been investigated for the flow past tube banks. Reynolds number, defined by freestream velocity ($U_{\infty}$) and cylinder diameter (d), was fixed as 200, and center-to-center distance (P) as 1.5d. In-line square array and staggered square array were considered. Drag coefficient, lift coefficient and Strouhal number were calculated depending on domain size. Circular cylinders were implemented on a Cartesian grid system by using an immersed boundary method. Boundary condition is periodic in both streamwise and lateral directions. Previous studies in literature often use a square domain with a side length of P, which contains only one cylinder. However, this study reveals that the domain size is improper. Especially, RMS values of flow-induced forces are most sensitive to the domain size.