• 설비공학논문집
• /
• 제28권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
• /
• 제4권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.

• 설비공학논문집
• /
• 제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.

• 한국전산유체공학회지
• /
• 제16권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.

• 대한기계학회논문집B
• /
• 제39권4호
• /
• pp.325-333
• /
• 2015

본 연구에서는 층류 점성유동에서 원형 실린더를 기준으로 홈을 이용한 실린더의 형상변화가 유동 특성에 미치는 영향 분석을 2차원 수치 해석으로 수행하였다. 사용한 좌표계는 직교 좌표계로서 가상경계법과 표시 입자 기법을 사용하여 실린더를 구현하였다. 홈의 개수와 위치를 변수로 Re=40, 50에서 계산하였으며, 홈의 수와 위치에 의한 유동 불안정성을 해석하여 임계 레이놀즈 수를 구하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제17권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.

• 한국전산유체공학회지
• /
• 제22권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).

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

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