• Transactions of the Korean Society of Automotive Engineers
• /
• v.11 no.2
• /
• pp.157-163
• /
• 2003

This work was undertaken for the numerical analysis of defrosting phenomena of automobile windshield. To analyze the defrost, the flow and temperature field of cabin interior, heat transfer through the windshield glass, and phase change of the frost should be analyzed simultaneously. The flow field was obtained by solving the 3-D unsteady Navier-Stokes equation and the temperature field was computed by energy equation. The phase-change process of Stefan problem was solved by enthalpy method. For code validation, the temperature field of the driven cavity was calculated. The result of calculation shows a good agreement with the other numerical results. Then, the present code was applied to the defrosting analysis of a real automobile and, also, a good agreement with experiment was obtained.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.53-61
• /
• 1999

Two-dimensional, unsteady, incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. The compressible code involves a conventional upwind-differenced scheme for the convective terms and LU-SGS scheme for temporal integration. The incompressible code with pseudo-compressibility method also adopts the same schemes as the compressible code. Three two-equation turbulence models are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by predicting the flow around the RAE 2822 transonic airfoil and the NACA 4412 airfoil, respectively. In addition, both the incompressible and compressible code are used to compute the flow over the NLR 7301 airfoil with flap to study the compressible effect near the high-loaded leading edge. The grid systems are efficiently generated using Chimera overlapping grid scheme. Overall, the κ-ω SST model shows closer agreement with experiment results, especially in the prediction of adverse pressure gradient region on the suction surfaces of high-lift airfoils.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.1
• /
• pp.35-44
• /
• 2017

Supercavitation is one of the most attractive technologies to achieve high speed for underwater vehicles. However, the multiphase flow with high-speed around the supercavitating vehicle (SCV) is difficult to simulate accurately. In this paper, we use modified the turbulent viscosity formula in the Standard K-Epsilon (SKE) turbulent model to simulate the supercavitating flow. The numerical results of flow over several typical cavitators are in agreement with the experimental data and theoretical prediction. In the last part, a flying SCV was studied by unsteady numerical simulation. The selected computation setup corresponds to an outdoor supercavitating experiment. Only very limited experimental data was recorded due to the difficulties under the circumstance of high-speed underwater condition. However, the numerical simulation recovers the whole scenario, the results are qualitatively reasonable by comparing to the experimental observations. The drag reduction capacity of supercavitation is evaluated by comparing with a moving vehicle launching at the same speed but without supercavitation. The results show that the supercavitation reduces the drag of the vehicle dramatically.

• Journal of Ship and Ocean Technology
• /
• v.10 no.4
• /
• pp.1-11
• /
• 2006

The finite volume based multi-block RANS code, WAVIS developed at MOERI is applied to the numerical self-propulsion test. WAVIS uses the cell-centered finite volume method for discretization of the governing equations. The realizable $k-{\epsilon}$ turbulence model with a wall function is employed for the turbulence closure. The free surface is captured with the two-phase level set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. The propeller forces are obtained using an unsteady lifting surface method based on potential flow theory. The numerical procedure followed the self-propulsion model experiment based on the 1978 ITTC performance prediction method. The self-propulsion point is obtained iteratively through balancing the propeller thrust, the ship hull resistance and towing force that is correction for Reynolds number difference between the model and full scale. The unsteady lifting surface code is also iterated until the propeller induced velocity is converged in order to obtain the propeller force. The self-propulsion characteristics such as thrust deduction, wake fraction, propeller efficiency, and hull efficiency are compared with the experimental data of the practical container ship. The present paper shows that hybrid RANS and potential flow based numerical method is promising to predict the self-propulsion parameters of practical ships as a useful tool for the hull form and propeller design.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.8 s.101
• /
• pp.962-967
• /
• 2005

The twin-impulse wave discharged from two parallel tubes is investigated to see flow patterns, compared with the single impulse wave. In the present study, the merging phenomena and propagation characteristics of the impulse waves are investigated by experiment and numerical computation. The Harten-Yee's total variation diminishing scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The Mach number Ms of incident shock wave is lower than 1.5 and the distance between the tubes is between 1.2 and 4.0. In the shock tube experiment, the twin impulse waves are visualized by a Schlieren optical system in order to validate the computational result. It is shown that on the symmetric axis between two parallel tubes, the peak pressure produced by the twin impulse waves and its location strongly depend upon the tube distance and the incident shock Mach number, Ms. The predicted Schlieren images show a good agreement with the measured twin-impulse wave.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.640-645
• /
• 2001

The present study is to investigate the characteristics of the impulse waves discharged from the exit of the convergent and divergent pipes. An experiment is carried out using a shock tube with an open end and is compared to the computation of the axisymmetric, compressible, unsteady Euler equations, which are solved by the second-order total variation diminishing(TVD) scheme. For the computational work, some initial compression waves are assumed inside the pipe so that those are identical to the experimental ones of the shock tube. The results show that the peak pressures of the impulse waves discharged from the exit of convergent and divergent pipes decrease with an increase in the wavelength of the initial compression wave. All of the impulse waves have a strong directivity toward the pipe axis, regardless of the exit type of the pipe employed. The impulse waves discharged from the divergent pipe are stronger than those from the straight pipe, while the impulse waves of the convergent pipe are weaker than those from the straight pipe. It is believed that the convergent pipe can playa role of a passive control to reduce the peak pressure of the impulse wave. The present computations represent the experimented impulse waves with a good accuracy.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.34 no.1
• /
• pp.87-99
• /
• 1992

A mathematical model, UNSATR which predicts the seepage flow through the body of dike especially under the tidal fluctuation has been developed. This model has been revised from UNSAT2 model which was developed on the basis of the saturated-unsaturated theory by Neuman. UNSATR has been verified and applied to the hydraulic model in order to estimated the seepage quantity, the formation of free water surface etc. The results lead to the following conclusions : 1. Seepage rates between the mathematical model and hydraulic model experiment are very similar to each other both in constant and transient water level conditions. 2. The lapsed time to be steady state of the free water surface becomes late as the tidal levels are relatively low mainly due to the seepage flow from the unsaturated zone of the body of dike. 3. Under the transient state of water levels, owing to the flow from the unsaturated domain, streamlines crossing to the free water surface are found and time lag during a falling tide may allow the free water surface inside the body of dike to stand at a high level than the outside water level. 4. The utility and validity of UNSATR model are convinced when the analyses on seepage problems through the porous embankment of the soil structures on the conditions of the steady and unsteady states are carried out.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.37 no.5
• /
• pp.53-61
• /
• 1995

Analyses of subsurface drainage effects of farmland with respect to pipe and envelop material are made by the laboratory experiments using soil box to give basic information for the subsurface drainage system planning and design. Three different diameter PVC perforated pipes and a mesh pipe are used with envelop materials such as sand, rice bran, and crushed stone. Steady state subsurface drainage flow rate increased as envelop material changed from sand to rice bran and crushed stone. This indicates that as the hydraulic conductivity of the envelop material increases, the drainage flow rate increases. On the other hand, for a given envelop material, the mesh pipe which has the largest openning area shows the largest flow rate while small diameter PVC pipes show small flow rates. This tells that as the openning area and pipe diameter increase, the flow rate increases, too. Therefore, selection of pipe and envelop material should be made in accordance with the design drainage flow rate. Unsteady state subsurface drainage flow rate with respect to time differs for different envelop material. In case the sand was used as an envelop material, the small diameter PVC pipes show larger flow rates than the large diameter PVC pipe and mesh pipe. When the rice bran was used, the mesh pipe shows the largest flow rate, while small diameter pipes show smaller flow rates. In case the crushed stone was used as an envelop material, the large diameter PVC pipe and mesh pipe show larger flow rates, while small diameter pipes show a little bit smaller flow rates. However, the variation of flow rates among different pipes is the smallest when the crushed stone is used. The flow rate curve with respect to the pipe changes little for the crushed stone envelop which has a large hydraulic conductivity, while that changes much for the sand and rice bran envelops. However, it is difficult to draw a consistent relationship between the drainage flow rate and pipe for all the envelop materials. Since the subsurface drainage experiments are made only under the restricted laboratory condition in this study, further study including field experiment is required.

• International Journal of Automotive Technology
• /
• v.6 no.1
• /
• pp.1-8
• /
• 2005

Recently, the HSDI (High Speed Direct Injection) diesel engine has been spotlighted as a next generation engine because it has a good potential for high thermal efficiency and fuel economy. This study was carried out to investigate the in-cylinder flow characteristics generated in a HSDI diesel engine with a 4-valve type cylinder head. The four kinds of cylinder head were manufactured to elucidate the effect of intake port geometry on the in-cylinder flow characteristics. The steady flow characteristics such as coefficient of flow rate $(C_{f})$, swirl ratio (Rs), and mass flow rate (m,) were measured by the steady flow test rig and the unsteady flow velocity within a cylinder was measured by PIV. In addition, the in-cylinder flow patterns were visualized by the visualization experiment and these results were compared with simulation results calculated by the commercial CFD code. The steady flow test results indicated that the mass flow rate of the cylinder head with a short distance between the two intake ports is $13\%$ more than that of the other head. However, the non-dimensional swirl ratio is decreased by approximately $15\%$. As a result of in-cylinder flow characteristics obtained by PIV and CFD calculation, we found that the swirl center was eccentric from the cylinder center and the position of swirl center was changed with crank angle. As the piston moves to near the TDC, the swirl center corresponded to the cylinder center and the velocity distribution became uniform. In addition, the results of the calculation are in good agreement with the experimental results.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.561-567
• /
• 2000

This paper is concerned with the viscous interaction between rotor and stator. The viscous interaction is caused by wakes from upstream blades. The rotor cascade in the experiment was composed with five blades, and cylinders were placed to make the stator wakes and their locations were about 50 percent upstream of blade chord. The locations of cylinders were varied in the direction of cascade axis with 0, 12.5, 25, 50, and 75 percent of pitch length. The static pressure distributions on the blade surfaces and the velocity distributions in the cascade flow were measured. From the experimental result it was found that the value of velocity defect by a cylinder wake might vary depending on the wake position within the cascade but the value at the cascade exit approached to some constant value regardless of the difference of wake locus. The momentum defect at the downstream from the cascade and the pressure distribution on the blade surfaces showed that the wake flowing near the blade surfaces caused the decrease of lift and the increase of drag regardless of the disappearance of flow separation.