• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.31-38
• /
• 1999

Numerical analysis is conducted on the deformation of free surface of magnetic fluid. Steady magnetic fields are induced by a circular current loop. Governing equations of magnetic fields are solved by using the concept of vector potential. The free surface of magnetic fluid is formed by the balance of surface force, gravity, pressure difference, magnetic normal pressure and magnetic body force. The deformations of free surface of magnetic fluid are qualitatively clarified. And, the patterns of steady non-uniform magnetic fields induced by a circular current loop are quantitatively presented. The shape of free surface attained by the polar fluid approach is rougher and higher than that attained by the quasi-steady approach.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.1
• /
• pp.23-29
• /
• 2008

Two high resolution compressive schemes, CICSAM(Ubbink, 1997) and HRIC(Muzaferija & Peric, 1999), in interface capturing method are reviewed briefly with respect to the extended forms suitable for unstructured meshes. And then those are applied to three typical test cases of translation test, shearing flow test and collapsing water problem with an obstacle. It is accomplished by implementing the high resolution schemes in the in-house CFD code(PowerCFD) for computing 3-D flow with an unstructured cell-centered method, which is based on the finite-volume technique and fully conservative. The calculated results show that CICSAM is better than HRIC with respect to accuracy and robustness, although either scheme can be used as a good choice for free surface or two-phase flow simulation.

• Journal of the Korean Society of Visualization
• /
• v.16 no.1
• /
• pp.54-60
• /
• 2018

It is important to develop more efficient and productive casting processes for an automated high precision molten-metal casting system. Detailed analysis of molten-metal flow in the casting process by the numerical approach will help to optimize the control of a ladle. In this study, the smoothed particle hydrodynamics method was applied to analyze casting flow characteristics with different tilting angular speed and initial molten-metal level. The smoothed particle hydrodynamics technique has advantages to easily handle non-linear free surface behavior with the absence of a computational mesh. We found that tilting angular speed has relatively greater effect on the casting flowrate and that the effect of the initial molten-metal level is only minor. Further extensive study will be necessary to find an optimal condition for high efficient casting system.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2401-2405
• /
• 2008

In the present study, a three-dimensional least square/level set based two-phase flow code was developed for the simulation of three-dimensional sloshing problems using finite element discretization. The present method can be utilized for the analysis of a free surface flow problem in a complex geometry due to the feature of FEM. Since the finite element method is employed for the spatial discretization of governing equations, an unstructured mesh can be naturally adopted for the level set simulation of a free surface flow without an additional load for the code development except that solution methods of the hyperbolic type redistancing and advection equations of the level set function should be devised in order to give a bounded solution on the unstructured mesh. From the numerical experiments of the present study, it is shown that the proposed method is both robust and accurate for the simulation of three-dimensional sloshing problems.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.247-252
• /
• 2001

The flow around free end of a finite circular cylinder(FC) embedded in an atmospheric boundary layer has been investigated experimentally. The experiments were carried out in a closed-return type subsonic wind tunnel with varying aspect ratio of the finite cylinder mounted vertically on a flat plate. The wake structures behind a 2-D cylinder and a finite cylinder located in a uniform flow were also measured for comparison. Reynolds number based on the cylinder diameter was about Re=20,000. A hot-wire anemometer was employed to measure the wake velocity and the mean pressure distributions on the cylinder surface were also measured. The flow past the FC free end shows a complicated three-dimensional wake structure and flow phenomenon is quite different from that of 2-D cylinder. The three-dimensional flow structure was attributed to the downwashing counter rotating vortices separated from the FC free end. As the FC aspect ratio decreases, the vortex shedding frequency is decreased and the vortex formation length is increased compared to that of 2-D cylinder. Due to the descending counter-rotating twin-vortex, in the region near the FC free end, regular vortex shedding from the cylinder is suppressed and the vortex formation region is hardly established. In the wake center region, the mean velocity for the FC located in atmospheric boundary layer has large velocity deficit, compared to that of uniform flow.

• Journal of Navigation and Port Research
• /
• v.30 no.8 s.114
• /
• pp.631-636
• /
• 2006

In this paper, the wave pattern around the hull with the transom stern advancing on the free surface with a constant speed was taken into consideration. To solve the problem the numerical analysis program was developed using Rankine source panel method based on potential flow analysis technique. The non-linearity of the free surface boundary conditions was fully satisfied. To verify the validity of the developed program the numerical calculations for Athena hull and KCS(KRISO container ship) hull was performed. The results of the numerical computation was compared with the ones of the model test experiment.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.11a
• /
• pp.90-95
• /
• 2004

During the past 50 years methods for predicting wave overtopping of coastal structures have coastal structures have continuously been developed Wave overtopping is one of the most important processes for the design of seawalls. The term 'wave overtopping' is used here to refer to the processes where waves hit a sloping structure run up the slope and, if the crest level of the slope is lower than the highest run up level, overtop the structure. Wave overtopping is dependent on the processes associated with breaking wave. The Numerical model is based on Navier-Stokes equation and Marker-Density Function of method for nonlinear free-surface flow by Miyata & Park(1995). The influence of how the slopes of seawalls, wave type and crest freeboard affect overtopping discharges has been investigated. The research of study using the new development nonlinear free-surface flow numerical model SOLA-VOF are presented.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1265-1273
• /
• 1999

PIV(Particle Image Velocimetry) velocity field measurements inside the snout of a1/10 scale model of the Zn plating process were carried out at the strip speed $V_s=1.5m/s$. Aluminum powder particles ($1{\mu}m$) and atomized olive oil ($3{\mu}m$) were used as seeding particles to simulate the molten Zinc flow and deoxidization gas flow, respectively. A pulsed Nd:Yag laser and a $2K{\times}2K$ high-resolution CCD camera were synchronized for the PIV velocity field measurement. From flow visualization study, it is found that the liquid flow in the Zn pot is dominantly governed by the uprising flow caused by the rotating sink roll, with its effect on the steel strip inside the snout largely diminished by installing of the snout. The deoxidization gas flow in front of the strip inside the snout can be characterized by a large-scale vortex rotating clockwise direction formed by the moving strip. In the rear side of the strip, a counter-clockwise vortex is formed and some of the flow entrained by the moving strip impinges on the free surface of molten zinc. The liquid flow in front of the strip is governed by the flow entering the snout, caused by the spinning sink roll. Just below the free surface a counter-clockwise vortex is formed near the snout wall. The moving strip affects dominantly the flow behind the strip inside the snout, and large amount of the liquid flow follows the moving strip toward the sink roll. The thickness of the flow following the strip is very thin in the front side due to the uprising flow, however thick boundary layer is formed in the rear side of the strip. Its thickness is increased as moving downstream toward the sink roll. Inside the snout, the deoxidization gas flow above the free surface is much faster than the liquid flow in the zinc pot. Due to the larger influx of the flow following the moving strip in the rear side of the strip, higher percentage of imperfection can be anticipated on the rear surface of the strip.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.3
• /
• pp.9-18
• /
• 1985

The numerical calculation for solving boundary-value problem related to potential flows with a free surface is carried out by application of the localized finite element method. Only forced motion of 2-D body in infinitely deep fluid is considered, although this schemes is equally applicable to any first order time-harmonic problems of similar nature. The infinite domain of the fluid is separated into the inner flow field and the outer flow field with common inter-surface boundary. The finite element method is applied to obtain the solution in the inner flow field and the Green functions are utilized to represent the solution in the outer flow field. At the inter-surface boundary, the continuity of the value of potential and the normal derivative of the potential(i.e. matching condition) is conserved. The present method has better computational efficiency than the previous LFEM and the integral equation method of Frank. This enhanced computational efficiency is presumably due to the fact that the present method gives a symmetric coefficient matrix and requires less computational time in calculating the influence coefficient matrix of Green function than the integral equation method. And the irregular frequency desen't exist because the uniqueness of the solution is assured by the such that the exact free surface condition is satisfied on the boundary of the localized finite element region(i.e. inner region). As an example of the above method, the hydrodynamic forces for the circular cylinder and the rectangular cylinders are calculated. In the computed results, the small number of singularity distribution segments($3{\sim}6$) give good result relative to Ursell's and Vugts'.

• Composites Research
• /
• v.31 no.6
• /
• pp.379-384
• /
• 2018

This paper aims to study flow characteristics of epoxy resin w.r.t. the sizing agents treated on the carbon fibers which have the same surface morphologies before sizing treatment. Dynamic contact angle (DCA) was measured to evaluate wettability of a single carbon fiber. Wicking test and Vacuum Assisted Resin Transfer Molding (VARTM) were performed to find relation between DCA measurement results and impregnation characteristics. In addition, surface properties of the carbon fibers such as surface free energy and chemical compositions were measured and interfacial shear strength (IFSS) between the carbon fiber and the resin were experimentally characterized by using micro-droplet tests. According to these experimental results, the sizing agent for carbon fibers should have appropriate level of surface free energy and good chemical compatibility with the resin to reconcile resin flow characteristics and interfacial strength.