• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.11
• /
• pp.717-726
• /
• 2016

This study evaluated the effect of aspect ratio of an enclosure with a heated inner circular cylinder on three-dimensional natural convection. The immersed boundary method was used to model the inner circular cylinder based on the finite volume method. The Rayleigh number was varied between $10^5$ and $10^6$, and the Prandtl number was maintained at 0.7. The aspect ratio of the three-dimensional enclosure was changed in steps of 1 within a range of 1-4 by increasing the width of the enclosure. In this study, the flow and thermal fields in the enclosure reached the steady state, and showed a mirror-symmetric pattern with respect to the center plane (x=0). In addition, the surface-averaged Nusselt number of the inner circular cylinder increased, while the total surface-averaged Nusselt number of the enclosure walls decreased with increase in the aspect ratio of the enclosure.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.204-213
• /
• 2000

This study deals with jet impingement, which is extensively used in the process industries to achieve intense heating, cooling or drying rates and also widely employed as a test flow for turbulent models due to its complex flow configuration, on a flat plate by numerical methods. In this calculation, the finite volume method was employed to solve the Navier-stokes equation based on the non-orthogonal coordinate with non-staggered variable arrangement. To get a better understanding for the fluid flow and heat transfer characteristics of the turbulent jet impingements, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model was adapted and compared with the experimental data and the result of standard $k-{\varepsilon}$ turbulent model. Numerical calculations were carried out with various flow rates, nozzle to plate distances. In the case of the axisymmetric jet impingement on a flat plate, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model showed better agreement with the experimental data than the standard $k-{\varepsilon}$ turbulent model in the prediction of the mean velocity profiles, the turbulent velocity profiles. the turbulent shear stress and the heat transfer rate. The highest heat transfer rate can be obtained when the impingement occurs within the potential core..

• Journal of Korea Water Resources Association
• /
• v.43 no.1
• /
• pp.105-117
• /
• 2010

This paper presents a new and simple technique to perform MUSCL reconstruction for solving 2D shallow water equations. The modified MUSCL scheme uses weighted area ratio to apply non-uniform grid in stead of the previous method that equally distributed the difference of conservation variables to each interface. The suggested method can physically reconstruct conservation variables in case of uniform grid as well as non-uniform grid. In this study, Unsplit scheme applicable to unstructured grid is used and efficient slope limiter of TVD scheme is used to control numerical oscillation which can be occurred in modified MUSCL scheme. For accurate and efficient treatment of bed slope term, the modified MUSCL scheme is coupled with the surface gradient method. The finite volume model applied to suggested scheme is verified through a comparison between numerical solution and laboratory measurements data such as the simulations of isolated building test case and Bellos's dam break test case.

• Journal of the Society of Naval Architects of Korea
• /
• v.46 no.5
• /
• pp.488-497
• /
• 2009

The freak wave, also known as New-Year-Wave in the north Atlantic, is relatively large and spontaneous ocean surface wave that can sink even large ships and destroy maritime structures. To understand oceanic conditions that develop freak waves, we simulated and generated two versions of scale-downed waves (1:64 and 1:42) in a numerical wave tank and compared the results with the experiment in wave flume. Both of the breaking and non-breaking waves were generated in the simulation. The numerical simulation was implemented based on the finite volume method and a genetic optimization algorithm. Random values were assigned as the initial values for the parameter in the control function, which produced signals representing the motion of wave-maker. The same signal obtained from the optimization process was used for both of the simulation and the experiment. By varying the object function and restrictions of the simulation, a best profile of design wave was selected based on the characteristics, height and period of simulated waves. Results showed that the simulation and experiment with the scale of 1:42 agreed better with freak waves in the natural condition. The presented simulation method will contribute to saving the time and cost for conducting subsequent response analyses of motion under freak waves in the course of the model test for ship and maritime structure.

• Journal of Power System Engineering
• /
• v.9 no.2
• /
• pp.81-87
• /
• 2005

This study is numerical analysis on the increasing temperature characteristics of vaporizer fin for liquefied natural gas with super low temperature. Existing LNG vaporizers use the direct contact heat transfer mode where the extreme super low temperature LNG of $-162^{\circ}C$ flows inside of the tubes and about $20^{\circ}C$air flows on outside of the fin. Recently, the vaporizers with great enhanced performance compared to conventional type have been developed to fulfill these requirements. The vaporizing characteristic of LNG vaporizer with air as heat source has a fixed iced. These characteristic cause a low efficiency in vaporizer, total plant cost and installing space can be increased. The vaporizing characteristics of LNG via heat exchanger with air are analytically studied for an air heating type vaporizer. This study is intended to supply the design data for the domestic fabrication of the thickness and angle vaporizer fin. Governing conservation equations for mass, momentum and energy are solved by STAR-CD based on an finite volume method and SIMPLE algorithm. Calculation parameter is fin thickness, setup angle and LNG temperature. If the vaporization performance of the early stage and late stage of operating is considered, the case of ${\phi}=90^{\circ}$ was very suitable. In this paper was estimated that the heat transfer was most promoted in case of THF=2mm.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.368-378
• /
• 2005

The effects of rotary oscillation on the unsteady laminar flow past a circular cylinder. are numerically investigated in the present study. The numerical solutions for the 20 Wavier-Stokes equation are obtained using a finite volume method Tn the framework of an overlapping grid system. The vortex formation behind a circular cylinder and the hydrodynamics of wake flows for different rotary oscillation conditions are analyzed from the results of numerical simulation. The lock-on region is defined as the region that the natural shedding frequency due to the Karmann Vortex shedding and the forcing frequency due to the forced oscillating a cylinder are nearly same, and the quasi-periodic states are observed around that region. At the intersection between lock-on and non-lock-on region the shedding frequency is bifurcated. After the bifurcation, one frequency fellows the forcing frequency($S_f$) and the other returns to the natural shedding frequency($St_0$). in the quasi-periodic states, the variation of magnitudes and relevant phase changes of $C_L$ with forcing phase are examined.

• Journal of the Korean Geotechnical Society
• /
• v.15 no.6
• /
• pp.219-238
• /
• 1999

In the present study, laboratory tests including the seepage-induced consolidation test, suction test, and desiccation shrinkage test are performed to investigate characteristics of the desiccation shrinkage in reclaimed hydraulic fills. Soil samples for laboratory tests are obtained from three sites (districts of Haenam, Kogeum and Koheung in Chunnam area). Desiccation shrinkage settlement caused by three dimensional volume change is numerically evaluated using finite difference technique based on the governing equation proposed by Abu-Hejleh & Znidarcic. Also characteristics of the desiccation shrinkage analyzed from the test results are used as input data for numerical evaluations. Further predicted total settlements including the self-weight consolidation settlement are compared with values measured at the site of Haenam district. Finally, effects of parameters related to the desiccation shrinkage on settlements are examined.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.11 no.7
• /
• pp.1218-1224
• /
• 2000

A half-height waveguide power combiner is designed and analyzed for Ka-band satellite application. The branch line directional coupler is utilized as a power combiner to achieve high port-to-port isolation and low insertion loss. The half height waveguide is adopted to reduce the volume and mass of a power combiner. In this paper a half height waveguide power combiner is designed and analyzed by FDTD and its performance is compared with that of a full-height waveguide power combiner. The designed half-height combiner having optimum order is manufactured and tested. The measurement shows that the designed half-height power combiner satisfies all the performance requirements (insertion loss less than 0.3 dB, reflection loss more than 20dB, port to port isolation more than 20 dB, and port to port phase difference within 5$^{\circ}$) in the satellite communication frequency band of 20.255 GHz to 21.255 GHz.

• Journal of Advanced Marine Engineering and Technology
• /
• v.36 no.8
• /
• pp.1091-1096
• /
• 2012

The aerodynamic characteristics of circular cylinder in a channel are studied to make clear the flow feature by solving the Navier-Stokes equation based on the finite volume method with unstructured grids. Reviews are made on with the vorticity, velocity, dynamic pressure, residual and drag, where the Reynolds numbers are 50 and 100. The flows for $Re{\succeq}50$ shows the vortex shedding in the wake, and the result is the same as the case of moving cylinder. The ground effect of flat bottom results in the growth of vortex, being generated in the upper side of the cylinder and elongated in the rear. As the cylinder approaches to wall, for example 0.6, the cylinder plays as a role of blockage to obstruct the flow between the cylinder and wall. The drag coefficients are compared with others' results to confirm the validity of the present numerical simulation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.6
• /
• pp.1886-1896
• /
• 1991

This paper presents a method to predict the flying height of the head slider in a hard disk drive. Quantitative predicts of the flying height according to the variations of the external load and the disk velocity have been done by numerical computation. In addition, the magnitude of the external load to keep flying height constant were also suggested. The Modified Reynolds' equation driven from hydrodynamic lubrication theory under slip flow condition was used to describe air-bearing system under the slider. To solve the equation, a Finite Volume Method (FVM) has been applied. To determine the final minimum flying height and pitch angle of the head slider, the Secant iteration method is used which update initial guess of the minimum flying height and pitch angle of the slider. In this study, the model head slider has been selected from a real hard disk drive which is equipped in many commercial personal computers. As a result, as the disk velocity increases at constant external load, the minimum flying height and the pitch angle increase due to the in crease of the air-bearing force at the bottom of slider.