• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 춘계학술대회논문집
• /
• pp.559-564
• /
• 2008

A prediction of binary droplets collision is important in the formation of falling drops and the evolution of sprays. The droplet velocity, impact parameter and drop-size ratio have influence on the interaction of the droplets. By the effect of these parameter, the collision processes are generated with the complicated phenomena. The droplet collision can be classified into four interactions such as the bouncing, coalescence, reflexive separation and stretching separation. In this study, the two-phase flow of the droplet collision was simulated numerically by using the Level Set method. 2D axi-symmetric simulations on the head-on collisions in the coalescence and reflexive separation, and 3D simulation on the off-center collisions in the coalescence and stretching separation were performed. These numerical results showed good agreements with the experimental and analytical results. For tracking the identity of droplets after the collision, transport equation for the volume fraction of the each initial droplet were used. From this, the identities of droplets were analyzed on the collision of droplets having different size.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.109-114
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• 한국전산유체공학회지
• /
• 제14권3호
• /
• pp.63-68
• /
• 2009

The motion of a rising liquid droplet is different that of a bubble motion. Treatment of liquid drops is more complex because internal motion must be considered. A 3D unstructured CFD code has been developed to solve incompressible N-S equation for the droplet simulation. This front-tracking consideration which the interface is tracked explicitly is very available to apply for not only exact interface topology but also the high schmidt number issue, such as $CO_2$ dissolution. This paper is forced on the zig-zag motion of the liquid droplet. The simulation shows that if the rising droplet is located at the corner of the zig-zag path, the velocity is low and shape of the droplet is more spherical shape, results in the less drag coefficient. Twin horse shoe vortexes behind the rising droplet are presented and the topology of the droplet is compared with an experimental result during one period of the path.

• Journal of Ship and Ocean Technology
• /
• 제7권4호
• /
• pp.32-40
• /
• 2003

Numerical simulation code is developed to study the formation of liquid drops from a nozzle into another quiescent liquid. The Navier-Stokes equations for two immiscible, incompressible, Newtonian fluids are solved on a fixed, staggered grid of cylindrical axisymmetric coordinates. Interfacial motion is captured using a Front-Tracking Method. The time variation of interfacial shape simulated by the code is in excellent agreement with experiments. Simulation results show that the viscosity ratio affects the size of the satellite drops.

• 한국가시화정보학회지
• /
• 제19권2호
• /
• pp.32-40
• /
• 2021

Vapor-driven solutal Marangoni flow is governed by the concentration distribution of solutes on a liquid-gas interface. Typically, the flow structure is investigated by particle image velocimetry (PIV). However, to develop a theoretical model or to explain the working mechanism, the concentration distribution of solutes at the interface should be known. However, it is difficult to achieve the concentration profile theoretically and experimentally. In this paper, to find the concentration distribution of solutes around 2D droplet, the reverse tracking method with an artificial neural network based on PIV data was performed. Using the method, the concentration distribution of solutes around a 2D droplet was estimated for actual flow data from PIV experiment.

• 한국전산유체공학회지
• /
• 제15권4호
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Korea-Australia Rheology Journal
• /
• 제21권1호
• /
• pp.59-69
• /
• 2009

The prediction of pressure drop for a droplet flow in a confined micro channel is presented using FE-FTM (Finite Element - Front Tracking Method). A single droplet is passing through 5:1:5 contraction - straight narrow channel - expansion flow domain. The pressure drop is investigated especially when the droplet flows in the straight narrow channel. We explore the effects of droplet size, capillary number (Ca), viscosity ratio ($\chi$) between droplet and medium, and fluid elasticity represented by the Oldroyd-B constitutive model on the excess pressure drop (${\Delta}p^+$) against single phase flow. The tightly fitted droplets in the narrow channel are mainly considered in the range of $0.001{\leq}Ca{\leq}1$ and $0.01{\leq}{\chi}{\leq}100$. In Newtonian droplet/Newtonian medium, two characteristic features are observed. First, an approximate relation ${\Delta}p^+{\sim}{\chi}$ observed for ${\chi}{\geq}1$. The excess pressure drop necessary for droplet flow is roughly proportional to $\chi$. Second, ${\Delta}p^+$ seems inversely proportional to Ca, which is represented as ${\Delta}p^+{\sim}Ca^m$ with negative m irrespective of $\chi$. In addition, we observe that the film thickness (${\delta}_f$) between droplet interface and channel wall decreases with decreasing Ca, showing ${\delta}_f{\sim}Ca^n$ Can with positive n independent of $\chi$. Consequently, the excess pressure drop (${\Delta}p^+$) is strongly dependent on the film thickness (${\delta}_f$). The droplets larger than the channel width show enhancement of ${\Delta}p^+$, whereas the smaller droplets show no significant change in ${\Delta}p^+$. Also, the droplet deformation in the narrow channel is affected by the flow history of the contraction flow at the entrance region, but rather surprisingly ${\Delta}p^+$ is not affected by this flow history. Instead, ${\Delta}p^+$ is more dependent on ${\delta}_f$ irrespective of the droplet shape. As for the effect of fluid elasticity, an increase in ${\delta}_f$ induced by the normal stress difference in viscoelastic medium results in a drastic reduction of ${\Delta}p^+$.

• 대한기계학회논문집
• /
• 제17권12호
• /
• pp.3187-3195
• /
• 1993

Combustion of liquid-fuelled combustion in a high-temperature vitiated-air stream was studied. The mathematical formulation comprise the application of Eulerian conservation equation to the gas phase and Lagrangian equation of droplet motion. The latter is coupled with a droplet-tracking technique (PSI-CELL Model) which regard the droplet phase as a source of mass, momentum, and energy to the gaseous phase. Reaction rate is determined by taking into account the Arrhenius reaction rate based on a single-step reaction mechanism. The calculated profiles show somewhat uncertainess at the upstream, but bases data for designing the combustor followed by 2-phase flow were obtained.

• 한국분무공학회지
• /
• 제6권4호
• /
• pp.31-38
• /
• 2001

The Holographic Particle Velocimetry system can be a promising optical tool for the measurements of three dimensional particle velocities. In this study, the holographic panicle velocimetry system was used to measure the sizes and velocities of droplets formed by a commercial full cone spray nozzle. Uncertainty analysis was performed to identify the sources of all relevant errors and to evaluate their magnitude. The droplet velocities ranged from 10.3 to 13.3 m/s with average uncertainty of ${\pm}1.6m/s$, which is ${\pm}14%$ of the mean droplet velocity. Compared with relatively small uncertainties of velocity components in the normal direction to the optical axis, the uncertainty of the optical axis component is ${\pm}3.6m/s$. This is due to the long depth of field of droplet images in the optical axis, which is inherent feature of holographic system using forward-scattering object wave of particles.

• 한국분무공학회지
• /
• 제7권4호
• /
• pp.23-31
• /
• 2002

In this study, the sizes and velocities of droplets in sprays were measured by image processing method from digital images of local region of sprays. The morphological method based on the Euclidean distance transform, Watershed separation, and perimeter image was adopted for the recognition and separation of overlapped particles. The match probability method was used for the particle tracking and pairing. The measurement results show that the present method may be reliable for the analysis of the motion and distribution of droplets produced by spray and atomization devices.