• Journal of the Korea Computer Graphics Society
• /
• v.15 no.1
• /
• pp.1-6
• /
• 2009

This paper presents a new CIP method for unstructured mesh to reduce the numerical dissipation. To reflect precise physical characteristics, CIP method updates both the physical quantity and the derivative information. The proposed method uses the Finite Volume Method(FVM) to solve the non-advection term of CIP equation. And we performed several experiments to improve the accuracy of third-order interpolation. Our result shows that our algorithm has less numerical dissipation than that of linear advection solver.

• Journal of the Korean Society of Clothing and Textiles
• /
• v.19 no.1
• /
• pp.142-148
• /
• 1995

In this study, FVM (Finite Volume Method) which is one of the 2-dimensional numerical approach has been conducted to anticipate the temperature distribution between skin and clothes by the change of air temperature and fabric characteristics including fabric thickness. Several experimental works have been done to understand the thermal insulation effect (If fabrics on a human body by measuring the averaged temperature in the air layer between skin and clothes or by measuring the thermal resistance of fabrics. However, the formal method is inconvenient to measure the temperature distribution in the air layer to evaluate the insulation rate of the clothes on the skin because the real size of the clearance between skin and the clothes is too small to place the temperature sensor, and in the Tatter method the relationship between human body and the fabrics are ignored. However, the numerical method will be very effective and economical way to evaluate the insulation efficiency of clothes when the computational result is in the reliable range. As the result of this study, the temperature change in the sir layer between skin and clothes was linear to the fabric thickness and this result coincides with many previous experimental results. Moreover, it is possible to predict the optimum fabric thickness for the best thermal insulation in the air layer between skin and clothes.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.09a
• /
• pp.800-807
• /
• 2010

In this study, a series of numerical analyses has been performed in order to evaluate the performance of a full-scale closed-loop vertical ground heat exchanger constructed in Wonju. The circulation pipe HDPE, borehole and surrounding ground were modeled using FLUENT, a finite-volume method (FVM) program, for analyzing the heat transfer process of the system. Two user-defined functions (UDFs) accounting for the difference in the temperatures of the circulating inflow and outflow water and the change of the surrounding ground temperature with depth were adopted in the FLUENT model. The thermal properties of materials estimated in laboratory were used in the numerical analyses to compare the thermal efficiency of the cement grout with that of the bentonite grout used in the construction. The results of the simulation provide a verification of the in situ thermal response test data. The numerical model with the ground thermal conductivity of 4W/mK yielded the simulation result closer to the in-situ thermal response test than with the ground thermal conductivity of 3W/mK. From the results of the numerical analyses, the effective thermal conductivities of the cement and bentonite grouts were obtained to be 3.32W/mK and 2.99 W/mK, respectively.

• Journal of the Korean Society for Precision Engineering
• /
• v.19 no.8
• /
• pp.92-99
• /
• 2002

This work investigated the optimal condition for an uniform deposition growth rate in the vertical cylindric CVD chamber. Heat transfer, surface chemical reaction and mass diffusion in the flow field of CVD chamber h,id been computed using Fluent v5.3 code. A SIMPLE based finite Volume Method (FVM) was adopted to solve the fully elliptic equations for momentum, temperature and concentration of a chemical species. The numerical analysis results show good agreements with the measurements obtained by N. Yoshikawa. The results obtained by the numerical analysis showed that the film growth rate in the center of a susceptor is increasing, as the inner flow approaches to the forced convection. To the contrast, as it approaches to the natural convection, that in the outside of a susceptor is increasing. As the Reynolds number increases, the uniformity may not hold due to the larger temperature gradient at a susceptor surface. Therefore, when the temperature gradient on the surface of a susceptor is zero, the film growth rate becomes uniform on most surface.

• Computers and Concrete
• /
• v.18 no.3
• /
• pp.319-336
• /
• 2016

In this study, both two- and three-dimensional (2D and 3D) finite-volume-based models were developed to analyze the heat transfer mechanisms through the porous structures of cellular concretes under steady-state heat transfer conditions and to investigate the differences between the 2D and 3D modeling results. The 2D and 3D reconstructed pore networks were generated from the microstructural information measured by 3D images captured by X-ray computerized tomography (X-CT). The computed effective thermal conductivities based on the 2D and 3D calculations performed on the reconstructed porous structures were found to be nearly identical to those evaluated from the 2D cross-sectional images and the 3D X-CT images, respectively. In addition, the 3D computed effective thermal conductivity was found to agree better with the measured values, in comparison with the 2D reconstruction and real cross-sectional images. Finally, the thermal conductivities computed for different reconstructed porous 3D structures of cellular concretes were compared with those obtained from 2D computations performed on 2D reconstructed structures. This comparison revealed the differences between 2D and 3D image-based modeling. A correlation was thus derived between the results of the 3D and 2D models.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.4
• /
• pp.274-280
• /
• 2017

Generally, global performance analysis in offshore platforms is performed using potential-based numerical tools, which neglect hydrodynamic viscous effects. In comparison with the potential theory, computational fluid dynamics (CFD) methods can take into account the viscous effects by solving the Navier-Stokes equation using the finite-volume method. The open-source field operation and manipulation (OpenFOAM) C++ libraries are employed for a finite volume method (FVM) numerical analysis. In this study, in order to apply CFD to the global performance analysis of a hull-mooring coupled system, we developed a numerical wave basin to analyze the global performance problem of a floating body with a catenary mooring system under regular wave conditions. The mooring system was modeled using a catenary equation and solved in a quasi-static condition, which excluded the dynamics of the mooring lines such as the inertia and drag effects. To demonstrate the capability of the numerical basin, the global performance of a barge with four mooring lines was simulated under regular wave conditions. The simulation results were compared to the analysis results from a commercial mooring analysis program, Orcaflex. The comparison included the motion of the barge, catenary shape, and tension in the mooring lines. The study found good agreement between the results from the developed CFD-based numerical calculation and commercial software.

• Ocean Systems Engineering
• /
• v.8 no.4
• /
• pp.381-407
• /
• 2018

This paper presents a hybrid numerical approach, which combines a two-phase Navier-Stokes model (NS) and the fully nonlinear potential theory (FNPT), for modelling wave-structure interaction. The former governs the computational domain near the structure, where the viscous and turbulent effects are significant, and is solved by OpenFOAM/InterDyMFoam which utilising the finite volume method (FVM) with a Volume of Fluid (VOF) for the phase identification. The latter covers the rest of the domain, where the fluid may be considered as incompressible, inviscid and irrotational, and solved by using the Quasi Arbitrary Lagrangian-Eulerian finite element method (QALE-FEM). These two models are weakly coupled using a zonal (spatially hierarchical) approach. Considering the inconsistence of the solutions at the boundaries between two different sub-domains governed by two fundamentally different models, a relaxation (transitional) zone is introduced, where the velocity, pressure and surface elevations are taken as the weighted summation of the solutions by two models. In order to tackle the challenges associated and maximise the computational efficiency, further developments of the QALE-FEM have been made. These include the derivation of an arbitrary Lagrangian-Eulerian FNPT and application of a robust gradient calculation scheme for estimating the velocity. The present hybrid model is applied to the numerical simulation of a fixed horizontal cylinder subjected to a unidirectional wave with or without following current. The convergence property, the optimisation of the relaxation zone, the accuracy and the computational efficiency are discussed. Although the idea of the weakly coupling using the zonal approach is not new, the present hybrid model is the first one to couple the QALE-FEM with OpenFOAM solver and/or to be applied to numerical simulate the wave-structure interaction with presence of current.

• Journal of Korea Water Resources Association
• /
• v.42 no.4
• /
• pp.271-279
• /
• 2009

The objective of this study is to develop the scheme to apply one-dimensional finite volume method (FVM) to natural river with complex geometry. In the previous study, FVM using the Riemann approximate solver was performed successfully in the various cases of dam-break, flood propagation, etc. with simple and rectangular cross-sections. We introduced the transform the natural into equivalent rectangular cross-sections. As a result of this way, the momentum equation was modified. The accuracy and applicability of newly developed scheme are demonstrated by means of a test example with exact solution, which uses triangular cross-sections. Secondly, this model is applied to natural river with irregular cross-sections and non-uniform lengths between cross-sections. The results shows that the aspect of flood propagation, location and height of hydraulic jump, and numerical solutions of maximum water level are in good agreement with the measured data. Using the developed scheme in this study, existing numerical schemes conducted in simple cross-sections can be directly applied to natural river without complicated numerical treatment.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.222-228
• /
• 2004

The effects of rotation on the unsteady laminar flow past a circular cylinder is numerically investigated in the present study. We obtained the numerical solutions for unsteady two-dimensional governing equation for the flow using two different numerical schemes. One is an accurate spectral method and another is finite volume method. Above all, the flow around a stationary circular cylinder is investigated to understand the basic phenomenon of flow separation, bluff body wake. Also, the validation of our own codes, expecially based on FVM, is carried out by the comparison of results obtained from our simulations using two different schemes and previous numerical and experimental studies. By the effect of rotation, the mean lift increases and drag deceases, which well represent the previous study.

• Journal of computational fluids engineering
• /
• v.6 no.1
• /
• pp.1-13
• /
• 2001

The natural convection heat transfer and solidification in a stratified pool are studied. The flow and heat transfer characteristics in a heat generating pool are compared between single-layered and double-layered pools. And local Nusselt number distributions on outer walls are obtained to consider thermal loads on a vessel wall. The cooling and solidification of Al₂O₃/Fe melt in a hemispherical vessel are simulated to study the mechanism of heat transfer and temperature distribution. A unstructured mesh is chosen for this study because of the non-orthogonality originated from the boundaries of double-layered pool. Interface between the layers is modeled to be fixed. With this assumption mass flux across the interface is neglected, but shear force and heat flux are considered by boundary conditions. The colocated cell-centered finite volume method is used with the Rhie-Chow interpolation to compute cell face velocity. To prevent non-physical solutions near walls in case body force is large the wall pressure is extrapolated by the way to include body force. The numerical solutions calculated by current method show that averaged downward heat flux of the double-layered pool increases compared to single-layered pool and maximum temperature occurs right below the interface of the layers.