• Journal of Korean Society of Environmental Engineers
• /
• v.33 no.4
• /
• pp.251-257
• /
• 2011

In the present study, a 3-dimensional (3D) numerical model was developed to mimic the micro porous structure of a geological $CO_2$ storage reservoir. Especially, 3D modeling technique assigning random pore size to a 3D micro porous structure was devised. Numerical method using CFD (computational fluid dynamics) was applied for the 3D micro porous structure to calculate supercritical $CO_2$ flow field. The three different configurations of 3D micro porous model were designed and their flow fields were calculated. For the physical conditions of $CO_2$ flow, temperature and pressure were set up equivalent to geological underground condition where $CO_2$ fluid was stored. From the results, the characteristics of the supercritical $CO_2$ flow fields were scrutinized and the influence of the micro pore configuration on the flow field was investigated. In particular, the pressure difference and consequent $CO_2$ permeability were calculated and compared with increasing $CO_2$ flow rate.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.1
• /
• pp.18-25
• /
• 2017

This study deals with numerically comparing performance according to rudder shape called 'Twisted rudder and Wavy twisted rudder'. In comparison with conventional rudder, rudder with wavy shape has showed a better performance at high angles of attack($30^{\circ}{\sim}40^{\circ}$) due to delaying stall. But most of study concerned with wavy shape had been performed in uniform flow condition. In order to identify the characteristics behind a rotating propeller, the present study numerically carries out an analysis of resistance and self-propulsion for KCS with twisted rudder and wavy twisted rudder. The turbulence closure model, Realizable $k-{\epsilon}$, is employed to simulate three-dimensional unsteady incompressible viscous turbulent and separation flow around the rudder. The simulation of self-propulsion analysis is performed in two step, because of finding optimization case of wavy shape. The first step presents there are little difference between twisted rudder and case of H_0.65 wavy twisted rudder in delivered power. So two kind of rudders are employed from first step to compare lift-to-drag ratio and torque at high angles of attack. Consequently, the wavy twisted rudder is presented as a possible way of delaying stall, allowing a rudder to have a better performance containing superior lift-to-drag ratio and torque than twisted rudder at high angles of attack. Also, as we indicate the flow visualization, check the quantity of separation flow around the rudder.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.6
• /
• pp.515-521
• /
• 2017

Large Cavitation Tunnel (LCT) of KRISO enables us to conduct cavitation tests of the propeller attached to a ship model. As the ship model tests are done at rather high Reynolds number of 107~108, flow measurement system such as pitot tube cannot be employed because of structural safety problems in its system and difficulties in installing it within the test section. Thus, KRISO has developed new 3-D LDV system used in large test section of LCT. There are several difficulties in using 3-D LDV, which did not allow efficient operation of it. The first trouble was the calibration using the conventional pin hole. To make the focus with same laser-beam waists at the wanted position, the high spatial resolution CCD is utilized in the calibration procedure for 3-D LDV. The off-axis configuration provides two velocity components in the horizontal plane and on-axis configuration gives third velocity component in the vertical plane. The horizontal velocity components are also obtained in the coincidence mode, which prevents any misleading results in the off-axis configuration. The nominal wake of Aframax tanker model is measured by the developed 3-D LDV system. The measured hull wake showed good agreement with that obtained by CFD calculation.

• Wind and Structures
• /
• v.23 no.4
• /
• pp.313-350
• /
• 2016

Taipei 101 Tower, which has 101 stories with height of 508 m, is located in Taipei where typhoons and earthquakes commonly occur. It is currently the second tallest building in the world. Therefore, the dynamic performance of the super-tall building under strong wind actions requires particular attentions. In this study, Large Eddy Simulation (LES) integrated with a new inflow turbulence generator and a new sub-grid scale (SGS) model was conducted to simulate the wind loads on the super-tall building. Three-dimensional finite element model of Taipei 101 Tower was established and used to evaluate the wind-induced responses of the high-rise structure based on the simulated wind forces. The numerical results were found to be consistent with those measured from a vibration monitoring system installed in the building. Furthermore, the equivalent static wind loads on the building, which were computed by the time-domain and frequency-domain analysis, respectively, were in satisfactory agreement with available wind tunnel testing results. It has been demonstrated through the validation studies that the numerical framework presented in this paper, including the recommended SGS model, the inflow turbulence generation technique and associated numerical treatments, is a useful tool for evaluation of the wind loads and wind-induced responses of tall buildings.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.169-172
• /
• 2002

When a plane shockwave reflects ken a concave wall, it is focused at a certain location, resulting in extremely high local pressure and temperature. This focusing is due to a nonlinear phenomenon of shock wave. The focusing phenomenon has been extensively applied to many diverse folds of engineering and medical treatment as well. In the current study, the focusing of shock wave over a reflector is numerically investigated using a CFD method. The Harten-Yee total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The incident shock wave Mach number $M_{s}\;of\;1.1{\~}l.3$ is applied to the parabolic reflectors with several different depths. Detailed focusing characteristics of the shock wave are investigated in terms of peak pressure, gasdynamic and geometrical foci. The results obtained are compared with the previous experimental results. The results obtained show that the peak pressure of shock wave focusing and its location strongly depend on the magnitude of the incident shock wave and depth of parabolic reflector. It is also found that depending up on the depth of parabolic reflector, the weak shock wave focusing process can classified into three distinct patterns : the reflected shock waves do not intersect each other before and after focusing, the reflected shock waves do not intersect each other before focusing, but intersect after focusing, and the reflected shock waves intersect each other before and after focusing. The predicted Schlieren images represent the measured shock wave focusing with a good accuracy.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.13 no.2
• /
• pp.124-130
• /
• 2014

In a mass manufacturing system of optical fibers, the sufficient cooling of glass fibers freshly drawn from a draw furnace is essential, asinadequately cooled glass fibers can lead to poor resin coating on the fiber surface and possibly fiber breakage during the process. In order to improve fiber cooling at a high drawing speed, it is common to use a helium injection into a glass fiber cooling unit in spite of the high cost of the helium supply. The present numerical analysis carried out three-dimensional thermo-fluid computations of the cooling gas flow and heat transfer on moving glass fiber to determine the cooling performance of glass fiber cooling depending on the method of helium injection. The results showed that afront injection of helium is most effective compared to a uniform or rear injection for reducing air entrainment into the unit and thus cooling the glass fibers at a high fiber drawing speed. However, above a certain amount of injected helium, there was no more increase of the cooling effect regardless of the helium injection method.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.15-23
• /
• 2016

In the present study, two phase flows around a projectile vertically launched from an underwater platform have been numerically investigated by using a three dimensional multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. The relative motion between the platform and projectile was described by six degrees of freedom equations of motion with Euler angles and a chimera technique. The propulsive power of the projectile was modeled as the fluid force acting on the lower surface of the body by the compressed air emitted from the underwater platform. Various flow conditions were considered to analyze the fluid-dynamics motion parameters of the projectile. The water level of platform and the current speed around the projectile were the main parametric variables. The numerical calculations were conducted up to 0.75sec in physical time scale. The dynamics tendency of the projectile was almost identical with respect to the water level variation due to the constant buoyancy term. The moving speed of the projectile along the vertical axis inside the platform decreased when the current speed increased. This is because the inflow from outside of the platform impeded development of the compressed air emitted from the floor surface of the launch platform. As a result, the fluid force acting on the lower surface of the projectile decreased, and injection time of the projectile from the platform was delayed.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.98-109
• /
• 2016

The flow and heat transfer characteristics of a heat exchanger having rectangular pin-fin in the flow direction have been investigated numerically. On the bottom plate, the convective boundary conditions for the hot side was given, and the fins were arranged in a channel-type geometric model using the periodic boundary condition in the span-wise direction. Three-dimensional numerical calculations for the flow and conjugate heat transfer problem were conducted using SIMPLE algorithm and $k-{\varepsilon}$ turbulence model. For the slanted pin-fin models, it was found that the downward cooling flow is generated due to the downward pressure gradient component, which can enhance the heat transfer performance near the bottom surface and the fin stem region. Four different inclined angles were considered in the Reynolds number range of 13,500-55,000. The aero-thermal performance of the slanted pin-fin heat exchangers, such as the volume and area goodness factors, were summarized and compared with the baseline plate-fin type heat exchanger quantitatively.

• Journal of Advanced Navigation Technology
• /
• v.16 no.2
• /
• pp.179-189
• /
• 2012

The aim of this paper is to research the change in the turbulent flow and the AOA(Angle Of Attack) occurred by $40^{\circ}$ crosswind to the direction of runway through the three-dimensional numerical analysis and to predict the take-off and landing flight stability. As a result, the maximum amplitude of AOA variation on runway reached $2^{\circ}$ within 3 second because of the wake formed by the constructions in the vicinity of the airport, and the overall effects appeared as an irregular aperiodic forms. Additionally, it was observed that the layout and shape of the buildings effected on the strength of turbulence directly, and the rapid flow generated between the buildings changed into stronger wake and eventually expected that the flow raises serious take-off and landing flight instability.

• The KSFM Journal of Fluid Machinery
• /
• v.11 no.1
• /
• pp.9-17
• /
• 2008

A numerical study for three-dimensional laminar flow in the entrance region of helical tubes connected with straight ones is carried out to investigate the effects of Reynolds number, pitch and curvature ratio on the oscillation periods of the flow. The fully elliptic governing equations were solved by means of a finite volume method. The fully developed laminar flow boundary condition was applied at the straight tube inlet. This results cover a curvature ratio range of 1/10${\sim}$1/320, a pitch range of 0.0${\sim}$3.2, and a Reynolds number range of 62.5${\sim}$2000. A comparison is made with previous experimental correlations and numerical data. The developments of velocity, local and average friction factors are discussed. The average friction factors are oscillatory in the entrance region of helical pipes. It has been found that the angle required for the flow to be similarly developed is most affected by the curvature ratio. The pitch and Reynolds number do not have any significant effect on the angle. The characteristic angle ${\phi}_c(={\phi}/sqrt{\delta})$, or the characteristic length to diameter ratio $s_c(=l\sqrt{\delta} cos(atan{\lambda})/d)$, can be useful to represent the development of flow in helical tubes. As the pitch increases and as the curvature ratio and Reynolds number decrease, the amplitude and the number of flow oscillations along the main streamwise direction decrease.