Comparison of Surface Tension Models in Smoothed Particles Hydrodynamics Method

2019 ◽  
Author(s):  
Lijing Yang ◽  
Milad Rakhsha ◽  
Dan Negrut
Keyword(s):  
Surface Tension ◽  
Liquid Droplet ◽  
Surface Force ◽  
Accurate Estimation ◽  
Droplet Deformation ◽  
Two Phase ◽  
Mesh Free ◽  
Interface Curvature ◽  
Particle Force ◽  
Smoothed Particles Hydrodynamics

Abstract We compare two surface tension models to solve two-phase fluid interaction problems in the context of the mesh-free Smoothed Particles Hydrodynamics (SPH) method. The Continuum Surface Force (CSF) model (later extended to Continuum Surface Stress, CSS), originally derived from grid-based numerical methods, requires an accurate estimation of the interface curvature to express the surface tension. Unlike CSF, the Inter-Particle Force (IPF) model is more robust in this regard as it draws on a molecular dynamics foundation by considering how the pairwise interaction forces between particles within a cutoff distance act in relation to producing the surface tension. Herein, we rely on second-order consistent gradient and Laplacian operators to improve the accuracy of SPH formulations as well as on a particle shifting technique to “disorder” particles from non-differentiable interface geometries. A 3D liquid droplet deformation test is used to compare CSF and IPF in terms of their pressure field and kinetic energy dissipation accuracy.

A Microfluidic Tensiometer

2004 ◽  
Author(s):  
Hans C. Mayer ◽  
Shelley L. Anna
Keyword(s):  
Surface Tension ◽  
Transport Mechanism ◽  
Microfluidic Devices ◽  
Dynamic Surface Tension ◽  
Two Phase ◽  
Surface Tension Data ◽  
Dynamic Surface ◽  
Theoretical Predictions ◽  
Surfactant Transport ◽  
Interface Curvature

Recent theoretical predictions indicate that a shift in surfactant transport mechanism from diffusion controlled to kinetically controlled occurs at highly curved interfaces where the length scale is on the same order as feature sizes in microfluidic devices. At present, experimental evidence of this shift in transport mechanism is lacking and this is due to the limitations on degree of interface curvature imposed by traditional methods of surface tension measurement. We show that the measurement of dynamic surface tension is possible at highly curved interfaces using a microfluidic tensiometer that utilizes glass micropipettes to control curvature dimension. Comparison of dynamic surface tension data from our microfluidic tensiometer with data obtained from traditional techniques will validate the theoretical arguments reported, and will improve understanding of two phase flows in microfluidic devices.

Surface Tension Model for Particle Method Using Inter-Particle Force Derived From Potential Energy

2012 ◽  
Author(s):  
Eiji Ishii ◽  
Taisuke Sugii
Keyword(s):  
Surface Tension ◽  
Fluid Flow ◽  
Dimensional Space ◽  
Three Dimensional ◽  
Flow Simulation ◽  
Particle Method ◽  
Surface Force ◽  
Continuous Surface ◽  
Particle Force ◽  
Surface Tension Model

Fluid-flow simulation within micro/nano spaces is essential for designing micro/nano devices, such as those in micro-electro-mechanical systems and nanoimprint processes. Surface tension is a dominant force in the fluid flow within micro/nano spaces. Surface-tension models can be classified into two groups: models based on continuous surface force in immiscible phases, and models based on inter-particle force in miscible phases. The surface-tension model based on inter-particle force for modeling interactions between materials would fit fluid-flow simulation within micro/nano spaces better than the surface-tension model based on continuous surface force. We developed a surface tension model using inter-particle force for use with a particle method in a past study. However, workings of inter-particle forces in miscible phases were not verified. Furthermore, accuracy in three-dimensional simulation needed to be verified. These subjects were verified in this study using simple benchmark tests. First, cohesion based on potential energy was validated to qualitatively check the workings of inter-particle force. The phase separation from the mixed two-phase flow due to inter-particle force was simulated. Next, the inter-particle force at the gas-liquid interface was quantitatively verified using the theory of the Young-Laplace equation; the pressure in a droplet was compared in two- and three-dimensional simulations, and the predicted pressures in a droplet agreed well with this theory. The inter-particle force at the gas-liquid-solid interface for the wall adhesion of a droplet was also verified; the results for wall adhesion in three-dimensional space agreed much better than that in two-dimensional space. We found that our surface tension model was useful for simulating the fluid flow within micro/nano spaces.

Flow Induced Dynamics of a Pinned Droplet on the Surface of a Channel

2005 ◽  
Author(s):  
Amirreza Golpaygan ◽  
Nasser Ashgriz
Keyword(s):  
Surface Tension ◽  
Stokes Equations ◽  
Liquid Droplet ◽  
Equilibrium Shape ◽  
High Surface ◽  
Droplet Surface ◽  
Navier Stokes ◽  
Fluid Viscosity ◽  
Droplet Deformation ◽  
Newtonian Liquids

Dynamic behavior of a droplet adhering to the surface of a channel has been modeled under the influence of surrounding fluid. The numerical solution is based on solving Navier-Stokes equations for Newtonian liquids. The study includes the effect of interfacial forces with constant surface tension, also effect of adhesion between the wall and droplet accounted by implementing contact angle at the wall. The Volume-Of-Fluid method is used to numerically determine the deformation of free surface. Droplet deformation and final shapes have been predicted. A reduction in the surface tension allows the droplet to deform much easier. However, an increase in the fluid viscosity, although increases the shear force on the droplet, may not result in the deformation at high surface tension. It is shown that deformation of droplet significantly influences structure of channel flow. Effects of liquid droplet and channel fluid properties, namely density and viscosity, inlet velocity, surface tension and channel geometry on dynamics of the problem have been studied. Two different outcomes have been considered: the first one droplet with equilibrium shape and the other one when breakup of the droplet occurs. The border line between the disintegration region and equilibrium region is determined for different droplet surface tensions.

Numerical Analysis of Droplet Dynamics Under Different Temperature and Cross-Flow Velocity Conditions

2012 ◽  
Vol 134 (4) ◽  
Author(s):  
Manabendra Pathak
Keyword(s):  
Surface Tension ◽  
Flow Velocity ◽  
Liquid Droplet ◽  
Cross Flow ◽  
Effect Of Temperature ◽  
Two Phase ◽  
Navier Stokes Equation ◽  
Droplet Dynamics ◽  
Cross Flow Velocity ◽  
Detachment Time

Numerical investigations have been performed to investigate the dynamics of a single liquid droplet formation in another immiscible cross-flow liquid for different values of cross-flow velocities and temperature difference between the two phases. The transient growth of the droplet and detachment time have been calculated by solving Navier-Stokes equation for two-phase flow using coupled level-set and volume of fluid method. The effect of temperature on surface tension has been incorporated in numerical simulation by modeling the surface tension as a linear function of temperature. The effects of cross-flow velocity and temperature in the evolution and detachment of the droplet have been investigated from the balance of different forces acting on the system. With the increase in cross-flow velocity, the diameter of the droplet decreases and droplet detachment time increases. While with the increase in cross-flow temperature, both the diameter and detachment time decreases.

scholarly journals Depth-Integrated Two-Phase Modeling of Two Real Cases: A Comparison between r.avaflow and GeoFlow-SPH Codes

2021 ◽  
Vol 11 (12) ◽  
pp. 5751
Author(s):  
Seyed Ali Mousavi Tayebi ◽  
Saeid Moussavi Tayyebi ◽  
Manuel Pastor
Keyword(s):  
Rock Avalanche ◽  
Two Phase ◽  
Simulation Code ◽  
Simulation Tools ◽  
Mesh Free ◽  
Particle Hydrodynamics ◽  
Landslide Evolution ◽  
Smoothed Particle ◽  
Hazard And Risk ◽  
Comparison Of The Results

Due to the growing populations in areas at high risk of natural disasters, hazard and risk assessments of landslides have attracted significant attention from researchers worldwide. In order to assess potential risks and design possible countermeasures, it is necessary to have a better understanding of this phenomenon and its mechanism. As a result, the prediction of landslide evolution using continuum dynamic modeling implemented in advanced simulation tools is becoming more important. We analyzed a depth-integrated, two-phase model implemented in two different sets of code to stimulate rapid landslides, such as debris flows and rock avalanches. The first set of code, r.avaflow, represents a GIS-based computational framework and employs the NOC-TVD numerical scheme. The second set of code, GeoFlow-SPH, is based on the mesh-free numerical method of smoothed particle hydrodynamics (SPH) with the capability of describing pore pressure’s evolution along the vertical distribution of flowing mass. Two real cases of an Acheron rock avalanche and Sham Tseng San Tsuen debris flow were used with the best fit values of geotechnical parameters obtained in the prior modeling to investigate the capabilities of the sets of code. Comparison of the results evidenced that both sets of code were capable of properly reproducing the run-out distance, deposition thickness, and deposition shape in the benchmark exercises. However, the values of maximum propagation velocities and thickness were considerably different, suggesting that using more than one set of simulation code allows us to predict more accurately the possible scenarios and design more effective countermeasures.

scholarly journals Comparison of Surface Tension Models for the Volume of Fluid Method

Processes ◽  
2019 ◽  
Vol 7 (8) ◽  
pp. 542 ◽  
Author(s):  
Kurian J. Vachaparambil ◽  
Kristian Etienne Einarsrud
Keyword(s):  
Surface Tension ◽  
Multiphase Flow ◽  
Capillary Rise ◽  
Surface Tension Force ◽  
Surface Force ◽  
Parallel Plates ◽  
Time Step ◽  
Mesh Resolution ◽  
Active Research ◽  
Spurious Currents

With the increasing use of Computational Fluid Dynamics to investigate multiphase flow scenarios, modelling surface tension effects has been a topic of active research. A well known associated problem is the generation of spurious velocities (or currents), arising due to inaccuracies in calculations of the surface tension force. These spurious currents cause nonphysical flows which can adversely affect the predictive capability of these simulations. In this paper, we implement the Continuum Surface Force (CSF), Smoothed CSF and Sharp Surface Force (SSF) models in OpenFOAM. The models were validated for various multiphase flow scenarios for Capillary numbers of 10 − 3 –10. All the surface tension models provide reasonable agreement with benchmarking data for rising bubble simulations. Both CSF and SSF models successfully predicted the capillary rise between two parallel plates, but Smoothed CSF could not provide reliable results. The evolution of spurious current were studied for millimetre-sized stationary bubbles. The results shows that SSF and CSF models generate the least and most spurious currents, respectively. We also show that maximum time step, mesh resolution and the under-relaxation factor used in the simulations affect the magnitude of spurious currents.

Simulation of Impaction Between Liquid Droplet and Solid Particles Based on SPH Method

2014 ◽  
Author(s):  
Shuai Meng ◽  
Qian Wang ◽  
Rui Yang
Keyword(s):  
Surface Tension ◽  
Solid Particle ◽  
Liquid Droplet ◽  
Particle Interaction ◽  
Solid Particles ◽  
Seed Particle ◽  
Liquid Droplets ◽  
Particle Hydrodynamics ◽  
Smoothed Particle ◽  
Spray Granulation

The phenomenon of impaction between liquid droplets and solid particles is involved in many scientific problems and engineering applications, such as impaction between sprayed droplet and solid particles in limestone injection desulfurization system and the collision between a droplet of the liquid to be granulated and a seed particle in fluidized bed spray granulation process. There are a lot of factors affected this phenomenon: droplet and particle size, momentum of both liquid droplet and solid particles, materials, surface conditions of the solid particles and so on. However the experimental or numerical researches have been done mostly pay attention to Specific application or process, so the impaction phenomenon has not been through studied, for example how different factors affected the impaction process with its effect on different applications. This paper focuses on the basic issue of interaction between droplet and solid particles. Three main factors were considered: ratio of diameter between the droplet and solid particle, relative velocity and the surface tension (including the contact angle between droplet and solid particle). All the study is based on simulation using SPH (smoothed particle hydrodynamics) method, and the surface tension is simulated by particle-particle interaction.

Study on C–S and P–R EOS in pseudo-potential lattice Boltzmann model for two-phase flows

2017 ◽  
Vol 28 (09) ◽  
pp. 1750120 ◽  
Author(s):  
Yong Peng ◽  
Yun Fei Mao ◽  
Bo Wang ◽  
Bo Xie
Keyword(s):  
Surface Tension ◽  
Lattice Boltzmann ◽  
Equations Of State ◽  
Density Ratio ◽  
Maximum Density ◽  
Lattice Boltzmann Model ◽  
Two Phase Flows ◽  
Two Phase ◽  
Spurious Currents ◽  
Pseudo Potential

Equations of State (EOS) is crucial in simulating multiphase flows by the pseudo-potential lattice Boltzmann method (LBM). In the present study, the Peng and Robinson (P–R) and Carnahan and Starling (C–S) EOS in the pseudo-potential LBM with Exact Difference Method (EDM) scheme for two-phase flows have been compared. Both of P–R and C–S EOS have been used to study the two-phase separation, surface tension, the maximum two-phase density ratio and spurious currents. The study shows that both of P–R and C–S EOS agree with the analytical solutions although P–R EOS may perform better. The prediction of liquid phase by P–R EOS is more accurate than that of air phase and the contrary is true for C–S EOS. Predictions by both of EOS conform with the Laplace’s law. Besides, adjustment of surface tension is achieved by adjusting [Formula: see text]. The P–R EOS can achieve larger maximum density ratio than C–S EOS under the same [Formula: see text]. Besides, no matter the C–S EOS or the P–R EOS, if [Formula: see text] tends to 0.5, the computation is prone to numerical instability. The maximum spurious current for P–R is larger than that of C–S. The multiple-relaxation-time LBM still can improve obviously the numerical stability and can achieve larger maximum density ratio.

scholarly journals Flow topologies in bubble-induced turbulence: a direct numerical simulation analysis

2018 ◽  
Vol 857 ◽  
pp. 270-290 ◽  
Author(s):  
Josef Hasslberger ◽  
Markus Klein ◽  
Nilanjan Chakraborty
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Turbulent Flows ◽  
Volume Fraction ◽  
Fluid Velocity ◽  
Small Scale ◽  
Flow Topology ◽  
Two Phase ◽  
Local Flow ◽  
Interface Curvature

This paper presents a detailed investigation of flow topologies in bubble-induced two-phase turbulence. Two freely moving and deforming air bubbles that have been suspended in liquid water under counterflow conditions have been considered for this analysis. The direct numerical simulation data considered here are based on the one-fluid formulation of the two-phase flow governing equations. To study the development of coherent structures, a local flow topology analysis is performed. Using the invariants of the velocity gradient tensor, all possible small-scale flow structures can be categorized into two nodal and two focal topologies for incompressible turbulent flows. The volume fraction of focal topologies in the gaseous phase is consistently higher than in the surrounding liquid phase. This observation has been argued to be linked to a strong vorticity production at the regions of simultaneous high fluid velocity and high interface curvature. Depending on the regime (steady/laminar or unsteady/turbulent), additional effects related to the density and viscosity jump at the interface influence the behaviour. The analysis also points to a specific term of the vorticity transport equation as being responsible for the induction of vortical motion at the interface. Besides the known mechanisms, this term, related to surface tension and gradients of interface curvature, represents another potential source of turbulence production that lends itself to further investigation.

Sign in / Sign up

Export Citation Format

Share Document