A new algorithm for surface tension model in moving particle methods

2007 ◽  
Vol 55 (3) ◽  
pp. 225-240 ◽  
Author(s):  
Shuai Zhang ◽  
Koji Morita ◽  
Kenji Fukuda ◽  
Noriyuki Shirakawa
Author(s):  
Eiji Ishii ◽  
Taisuke Sugii

Predicting the spreading behavior of droplets on a wall is important for designing micro/nano devices used for reagent dispensation in micro-electro-mechanical systems, printing processes of ink-jet printers, and condensation of droplets on a wall during spray forming in atomizers. Particle methods are useful for simulating the behavior of many droplets generated by micro/nano devices in practical computational time; the motion of each droplet is simulated using a group of particles, and no particles are assigned in the gas region if interactions between the droplets and gas are weak. Furthermore, liquid-gas interfaces obtained from the particle method remain sharp by using the Lagrangian description. However, conventional surface tension models used in the particle methods are used for predicting the static contact angle at a three-phase interface, not for predicting the dynamic contact angle. The dynamic contact angle defines the shape of a spreading droplet on a wall. We previously developed a surface tension model using inter-particle force in the particle method; the static contact angle of droplets on the wall was verified at various contact angles, and the heights of droplets agreed well with those obtained theoretically. In this study, we applied our surface tension model to the simulation of a spreading droplet on a wall. The simulated dynamic contact angles for some Weber numbers were compared with those measured by Šikalo et al, and they agreed well. Our surface tension model was useful for simulating droplet motion under static and dynamic conditions.


2008 ◽  
Vol 183 (1) ◽  
pp. 21-26 ◽  
Author(s):  
Guangzheng Zhou ◽  
Wei Ge ◽  
Jinghai Li

Author(s):  
Eiji Ishii ◽  
Taisuke Sugii

The particle method is a useful approach to simulate fluid flows within micro/nano spaces such as micro-electromechanical systems, nano-in-print processes, and head-disk interfaces of hard disk drives. Particle methods are based on continuum dynamics, and some studies have recently extended the scope of these methods to approaches within micro/nano spaces. Surface tension is a dominant force in the fluid flow within micro/nano spaces. However, surface-tension models used in the particle methods need to be improved to achieve more stable and accurate simulation. In the present study, we developed a new surface tension model for the particle method using inter-particle force to improve the stability and accuracy of simulation; the inter-particle force was given by the derivation of potential energy in space. The developed surface tension model was verified using simple benchmark tests: pressure in a round droplet and oscillation period of a square liquid-droplet. The predicted pressure in a round droplet agreed well with that given by the Young-Laplace equation, and the predicted oscillation period of a square droplet agreed well with that given by Lamb’s theory. The wall-adhesion was also verified at various contact angles; heights of droplets on the wall agreed well with those given theoretically. We found that our new surface tension model was useful for simulating fluid flow within micro/nano spaces for particle method.


Author(s):  
Kailun Guo ◽  
Ronghua Chen ◽  
Suizheng Qiu ◽  
Wenxi Tian ◽  
Guanghui Su ◽  
...  

Multiphase flow widely exists in the nature and engineering. The two-phase flow is the highlight of the studies about the flow in the vessel and steam explosion in nuclear severe accidents. The Moving Particle Semi-implicit (MPS) method is a fully-Lagrangian particle method without grid mesh which focuses on tracking the single particle and concerns with its movement. It has advantages in tracking complex multiphase flows compared with gird methods, and thus shows great potential in predicting multiphase flows. The objective of this thesis is to develop a general multiphase particle method based on the original MPS method and thus this work is of great significance for improving the numerical method for simulating the instability in reactor severe accident and two-phase flows in vessel. This research is intended to provide a study of the instability based on the MPS method. Latest achievements of mesh-free particle methods in instability are researched and a new multiphase MPS method, which is based on the original one, for simulating instability has been developed and validated. Based on referring to other researchers’ papers, the Pressure Poisson Equation (PPE), the viscosity term, the free surface particle determination part and the surface tension model are optimized or added. The numerical simulation on stratification behavior of two immiscible flows is carried out and results are analyzed after data processing. It is proved that the improved MPS method is more accurate than the original method in analysis of multiphase flows. In this paper, the main purposes are simulating and discussing Rayleigh-Taylor (R-T) instability and Kelvin-Helmholtz (K-H) instability. R-T and K-H instability play an important role in the mixing process of many layered flows. R-T instability occurs when a lower density fluid is supported by another density higher fluid or higher density fluid is accelerated by lower density fluid, and the resulting small perturbation increases and eventually forms turbulence. K-H instability is a small disturbance for two different densities, such as waves, at the interface of the two-phase fluid after giving a fixed acceleration in the fluid. Turbulence generated by R-T instability and K-H instability has an important effect in applications such as astrophysics, geophysics, and nuclear science.


Author(s):  
Rida S. N. Mahmudah ◽  
Masahiro Kumabe ◽  
Takahito Suzuki ◽  
LianCheng Guo ◽  
Koji Morita ◽  
...  

Understanding the freezing behavior of molten metal in flow channels is of importance for severe accident analysis of liquid metal reactors. In order to simulate its fundamental behavior, a 3D fluid dynamics code was developed using Finite Volume Particle (FVP) method, which is one of the moving particle methods. This method, which is fully Lagrangian particle method, assumes that each moving particle occupies certain volume. The governing equations that determine the phase change process are solved by discretizing its gradient and Laplacian terms with the moving particles. The motions of each particle and heat transfer between particles are calculated through interaction with its neighboring particles. A series of experiments for fundamental freezing behavior of molten metal during penetration on to a metal structure was also performed to provide data for the validation of the developed code. The comparison between simulation and experimental results indicates that the present 3D code using the FVP method can successfully reproduce the observed freezing process such as molten metal temperature profile, frozen molten metal shape and its penetration length on the metal structure.


Sign in / Sign up

Export Citation Format

Share Document