10. 4  Reliability of Two-Phase Simulation is Increased!
- Development of a New Method to Calculate Surface Tension in Lattice Gas Model
 


Fig. 10-9 An example of grids in the LGM simulation

Only one artificial particle can stay on one grid point. If two particles come into one grid point the directions of the particles are changed according to a given collision rule.


Fig. 10-10  A liquid drop obtained by the LGM simulation


Fig. 10-11 A new principle of the new method to calculate the surface tension in the LGM simulation and an example of calculations

The surface tension of the liquid drop is calculated by fitting a part of the surface to a circle and calculating it for a combination of two cases, i.e., the case (a) where the tangential direction is not parallel to the grid edge and the parallel case (b). Subfigure (c) shows the comparison with the result obtained by the Laplace method. For large radius the result of the new method agrees well with that of the Laplace method and for small radius region the Laplace method cannot be used.


Motion of fluids as liquid or gas are traditionally analyzed by solving partial differential equations named the Navier-Stokes equations. Recently, however, methods to solve such fluid motions by lattice gas model (LGM) are often studied by many researchers. The LGM method can deal with a two-phase simulation rather easily, where, phase changes between two different phases such as gas and liquid are important. To carry out the simulation of fluid by LGM, at first, a computational domain is divided into regular grids such as two-dimensional triangular grids (Fig. 10-9), and then many artificial particles placed on the grid points are moved according to a given rule for collisions among the particles. In a result of the simulation of the artificial particles we can see fluid behavior which is the same as that of the ordinary fluid. Even in the case with phase changes there is no technological difficulty such as is found when the partial differential equations are solved.
However, when we use this method to carry out a simulation of liquid drop formation in a gas the shape of the liquid drop does not become circular (Fig. 10-10). It is because the surface tension of the liquid drop is not uniform on the surface due to the discrete structure of the space by the grids. Because in the calculation of the problems including phase change the role of the surface tension is very important, it is necessary that surface tension at phase boundaries with various shapes can be always calculated. Though, in the previous method to calculate the surface tension, it was necessary to introduce an arbitrary parameter, we developed a new method by which a unique solution can be obtained without an arbitrary parameter. By using the new method, therefore, the range of applicability of the LGM is made clear and the reliability of the method is increased (Fig. 10-11).



Reference
K. Ebihara et al., A Method for Calculating the Surface Tension of a Droplet in a Lattice-Gas Model with Long-Range Interaction, Eur. Phys. J. B, 18, 319 (2000).

Select a topic in left column


Persistent Quest - Research Activities 2001
Copyright(c) Japan Atomic Energy Research Institute