• Journal of the Korean Society of Marine Environment & Safety
• /
• v.20 no.3
• /
• pp.312-317
• /
• 2014

In this paper, the wave run-up height and depression depth around air-water interface-piercing circular cylinder have been numerically studied. The Reynolds Averaged Navier-Stokes equations (RANS) and continuity equations are solved with Reynolds Stress model (RSM) and volume of fluid (VOF) method as turbulence model and free surface modeling, respectively. A commercial Computational Fluid Dynamics (CFD) software "Star-CCM+" has been used for the current simulations. Various Froude numbers ranged from 0.2 to 1.6 are used to investigate the change of air-water interface structures around the cylinder and experimental data and theoretical values by Bernoulli are compared. The present results showed a good agreement with other studies. Kelvin waves behind the cylinder were generated and its wave lengths are longer as Froude numbers increase and they have good agreement with theoretical values. And its angles are smaller with the increase of Froude numbers.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.10 s.253
• /
• pp.1003-1011
• /
• 2006

Droplets are ejected onto a substrate through a nozzle by pushing liquids in flow channels of drop-on-demand devices. The behavior of ejection and formation of droplets is investigated to enhance the physical understanding of the hydrodynamics involved in inkjet printing. The free surface phenomenon of a droplet is described using $CFD-ACE^{TM}$ which employs the volume-of-fluid (VOF) method with the piecewise linear interface construction (PLIC). Droplet formation characteristics are analyzed in various flow regimes with different Ohnesorge numbers. The computational results show that the droplet formations are strongly dependent on the physical properties of working fluids and the inlet flow conditions. In addition, the wetting characteristics of working fluids on a nozzle influence the volume and velocity of a droplet produced in the device. This study may provide an insight into how a liquid droplet is formed and ejected in a piezoelectric inkjet printing device.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.568-576
• /
• 2005

A 3-D finite element model combined with a volume tracking method is presented in this work to simulate the mold filling for casting processes. Especially, the analysis involves an adaptive grid method that is created under a criterion of element categorization of filling states and locations in the total region at each time step. By using an adaptive grid wherein the elements, finer than those in internal and external regions, are distributed at the surface region through refinement and coarsening procedures, a more efficient analysis of transient fluid flow with free surface is achieved. Adaptive grid based on VOF method is developed in tetrahedral element system. Through a 3-D analysis of the benchmark test of the casting process, the efficiency of the proposed adaptive grid method is verified. Developed FE code is applied to a typical industrial part of the casting process such as aluminum road wheel.

• 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.37 no.6
• /
• pp.461-466
• /
• 2004

This study presents a combined experimental and numerical effort to investigate wave breaking of sinusoidal waves in a surf zone. Numerical predictions are verified by comparing to laboratory measurements. The model solves the Reynolds equations and$textsc{k}$-$\varepsilon$ models for the turbulence analysis. To track the free surface displacement, the volume of fluid method is employed. As the height of incident wave increases, the wave breaking occurs at a closer point of the slope in the numerical model and laboratory experiments with the same depth and period. When a wave breaking occurs, the ratio of wave height becomes larger, with the same wave height and depth, as the period increases.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1037-1048
• /
• 2023

Fishways are constructed to ensure the fish migration because river-crossing structures such as dams and weirs cut off the stream longitudinal connectivity and influence on aquatic ecosystems. However, the passage efficiency of fishes varies depending on flow characteristics in the fishway and fish species. In this study, three-dimensional numerical simulations are carried out using a RANS model and the volume of fluid method for resolving free surface fluctuations to calculate the turbulent flow in the pool and weir type fishway. The Flow velocity and turbulent kinetic energy in the pool of fishway are analyzed according to variation of the upstream water level and the length of pool. The present numerical simulations reasonably well reproduce the stream flow and plunging flow characteristics in the pool. The simulation results show that the stream flow changes to the plunging flow as the length of the pool increases. When the upstream level increases, the stream flow becomes more evident. Key parameters related to the fish migration within the fishway such as the flow velocity and the turbulent kinetic energy are examined to assess the ascending possibility of fishes.

• The Journal of the Acoustical Society of Korea
• /
• v.42 no.3
• /
• pp.250-261
• /
• 2023

In this study, the noise contributions by the duct, stator and rotor, which are the propulsor components, are evaluated to identify the flow noise source in cavitation and non-cavitation conditions on pumpjet propulsion and the noise levels in both conditions are compared. The unsteady incompressible Reynolds averaged Navier-Stokes (RANS) equation based on the homogeneous mixture assumption is applied on the suboff submarine hull and pumpjet propeller in the cavitation tunnel, and the Volume of Fluid (VOF) method and Schnerr-Sauer cavitation model are used to describe the two-phase flow. Based on the flow simulation results, the acoustic analogy formulated by Ffowcs Williams and Hawkings (FW-H) equation is applied to predict the underwater radiated noise. The noise contributions are evaluated by using the three types of impermeable integral surface on the duct, stator and rotor, and the two types of permeable integral surface surrounding the propulsor. As a result of noise prediction, the contribution by the stator is insignificant, but it affects the generation of flow noise source due to flow separation in the duct and rotor, and the noise is predominantly radiated into the upward and right where the flow separations are. Also, the noise is radiated into the thrust direction due to pressure fluctuation between suction and pressure sides on the rotor blades, and the it can be seen that the cavitation effect into the noise can be considered through the permeable integral surface.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.4
• /
• pp.508-516
• /
• 2010

This study aims at investigating the effect of the chamfer on the liquid sloshing in the three-dimensional (3D) rectangular tank. In order to simulate the 3D incompressible viscous two-phase flow in the 3D tank with partially filled liquid, the present study has adopted the volume of fluid (VOF) method based on the finitevolume method which has been well verified by comparing with the results of the relevant previous researches. The effects of the chamfering top corners of the tank on the liquid sloshing characteristics have been investigated. The angle of the chamfering top corners (${\theta}$) has been changed in the range of $0^{\circ}{\leq}{\theta}{\leq}60^{\circ}$(${\Delta}{\theta}=15^{\circ}$) to observe the free surface behavior, and the effect on wall impact load. Generally, as the angle of the chamfering top corners increases, the impact pressure on the upper knuckle point decreases. However it seemed that a critical angle of the chamfering top corners exists to reveal the lowest impact pressure on the wall.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.252-255
• /
• 2008

In this study, the characteristics of crossing a river of Ground Armored Vehicle (GAV) were evaluated by numerical method and real size tests. 3-D hybrid mesh systems were constructed by 3-D models of the GAV, and a commercial software, FLUENT, was used in numerical analysis. In order to deal with multi-phase problem (air and water), Volume Of Fluid (VOF) method was used, and Moving and Deforming Mesh (MDM) was adapted for unsteady motion of GAV. There were two steps in this research. Firstly, stability of the GAV which cruised a river was evaluated by changing several shapes of water-proof-front-wing of the GAV in steady state, and compared results (free surface shape and drag value in 10km/h) with those of real size tests. Secondly, results of unsteady analysis considering weight and moment of inertia of the GAV were presented. There were showed a maximum velocity with a designed water jet and dynamic stability including pitch, roll, and yaw moment. Based on these results, the optimal shape of water-proof-front-wing of the GAV was determined for a proto-type of the GAV.

• Journal of the Semiconductor & Display Technology
• /
• v.9 no.4
• /
• pp.19-23
• /
• 2010

Nanoimprint lithography (NIL) is one of the most versatile and promising technology for micro/nano-patterning due to its simplicity, high throughput and low cost. Recently, one of the major trends of NIL is large-area patterning. Especially, the research of the application of NIL to TFT-LCD field has been increasing. Technical difficulties to keep the uniformity of the residual layer, however, become severer as the imprinting area increases. In this paper we performed a numerical study for a large area NIL (the $2^nd$ generation TFT-LCD glass substrate ($370{\times}470$ mm)) by using finite element method. First, a simple model considering the surrounding wall was established in order to simulate effectively and reduce the computing time. Then, the volume of fluid (VOF) and grid deformation method were utilized to calculate the free surfaces of the resist flow based on an Eulerian grid system. From the simulation, the velocity fields and the imprinting pressure during the filling process in the NIL were analyzed, and the effect of the surrounding wall and the uniformity of residual layer were investigated.