Numerical simulation of wavy falling film flow using VOF method.

*(English)*Zbl 1047.76569Summary: Surface wave dynamics of vertical falling films under monochromatic-frequency flowrate-forcing perturbations is computed by the direct simulation of Navier-Stokes equations using the Volume of Fluid (VOF) method to track free surfaces and the Continuum Surface Force (CSF) model to account for dynamic boundary conditions at free surfaces. The numerical VOF–CSF model is completely formulated, and more attention is given to understanding instabilities of thin films. At low frequency and high flowrate, the small inlet disturbance develops into large solitary waves preceded by small capillary bow waves. The circulation flow compatible with the solitary wave size is observed in the solitary peak. On the other hand, at high frequency and low Re, small-amplitude waves in nearly sinusoidal shape without fore-running capillary waves are formed on the surface. The quasi-periodic waveforms are found to occur at the nearly sinusoidal wave regime. The slight increase in wave-amplitude and wavelength, and decrease in residual thickness as waves evolves downstream are observed for both solitary waves and sinusoidal types. The variation of velocity and pressure along a wave are strong at the wave trough and capillary wave region, due to the large surface curvature there. The pressure variation perpendicular to the wall is negligible and only a small variation is observed at the solitary wave trough and capillary region.

##### MSC:

76M25 | Other numerical methods (fluid mechanics) (MSC2010) |

76A20 | Thin fluid films |

76M20 | Finite difference methods applied to problems in fluid mechanics |

PDF
BibTeX
XML
Cite

\textit{D. Gao} et al., J. Comput. Phys. 192, No. 2, 624--642 (2003; Zbl 1047.76569)

Full Text:
DOI

##### References:

[1] | Chang, H.-C., Wave evolution on a falling film, Annu. rev. fluid mech., 26, 103, (1994) |

[2] | Liu, J.; Gollub, J.P., Solitary wave dynamics of film flows, Phys. fluids, 6, 1702, (1994) |

[3] | Yih, C.-S., Stability of liquid flow down an inclined plane, Phys. fluids, 6, 321, (1963) · Zbl 0116.19102 |

[4] | Liu, J.; Paul, J.D.; Gollub, J.P., Measurements of the primary instabilities of film flows, J. fluid mech., 250, 69, (1993) |

[5] | Benney, B.J., Long waves in liquid films, J. math. phys., 45, 150, (1966) · Zbl 0148.23003 |

[6] | Nakaya, C., Long waves on a thin fluid layer flowing down an inclined plane, Phys. fluids, 18, 1407, (1975) · Zbl 0325.76038 |

[7] | Chang, H.-C., Evolution of nonlinear waves on vertical falling films – a normal form analysis, Chem. eng. sci., 42, 515, (1987) |

[8] | P.L. Kapitza, S.P. Kapitza, Wave flow in thin layers of a viscous fluid, in: D. ter Haar (Ed.), Collected Papers of P.L. Kapitza (vol. II), The Macmillan Company, New York, 1964 · Zbl 0029.17603 |

[9] | Krantz, W.B.; Goren, S.L., Stability of thin liquid films flowing down a plane, Ind. eng. chem. fundam., 10, 91, (1971) |

[10] | Alekseenko, S.V.; Nakoryakov, V.Ye; Pokusaev, B.G., Wave formation on a vertical falling liquid film, Aiche j., 31, 1446, (1985) · Zbl 0578.76030 |

[11] | Nosoko, T.; Yoshimura, P.N., Characteristics of two-dimensional waves on a falling liquid film, Chem. eng. sci., 51, 725, (1996) |

[12] | Liu, J.; Schneider, J.B.; Gollub, J.P., Three-dimensional instabilities of film flows, Phys. fluids, 7, 55, (1995) |

[13] | Ho, L.-W.; Patera, A.T., A Legendre spectral element method for simulation of unsteady incompressible viscous free-surface flows, Comput. methods appl. mech. eng., 80, 355, (1990) · Zbl 0722.76052 |

[14] | Malamataris, N.T.; Papanastasiou, T.C., Unsteady free surface flows on truncated domains, Ind. eng. chem. res., 30, 2211, (1991) |

[15] | Salamon, T.R.; Armstrong, R.C.; Brown, R.A., Traveling waves on vertical films: numerical analysis using the finite element method, Phys. fluids, 6, 2202, (1994) · Zbl 0844.76024 |

[16] | Ramaswamy, B.; Chippada, S.; Joo, S.W., A full-scale numerical study of interfacial instabilities in thin-film flows, J. fluid mech., 325, 163, (1996) · Zbl 0889.76021 |

[17] | S. Krishnamoorthy, B. Ramaswamy, S.W. Joo, Three-dimensional instabilities in heated falling films: a full-scale direct numerical simulation [Conference Paper], American Society of Mechanical Engineers, Fluids Engineering Division (Publication) FED. v 238, n 3, 1996, ASME, New York, NY, USA, pp. 415-420 |

[18] | Hirt, C.W.; Nichols, B.D., Volume of fluid (VOF) method for the dynamics of free boundaries, J. comput. phys., 39, 202, (1981) · Zbl 0462.76020 |

[19] | D.B. Kothe, R.C. Mjolsness, M.D. Torrey, RIPPLE: A Computer Program for Incompressible Flows with Free Surfaces, LA-12007-MS, Los Alamos National Laboratory, 1991 |

[20] | Rider, W.J.; Kothe, D.B., Reconstructing volume tracking, J. comput. phys., 141, 112, (1998) · Zbl 0933.76069 |

[21] | Brackbill, J.U.; Kothe, D.B.; Zemach, C., A continnum method for modeling surface tension, J. comput. phys., 100, 335, (1992) · Zbl 0775.76110 |

[22] | Puckett, E.G.; Almgren, A.S.; Bell, J.B., A high-order projection method for tracking fluid interfaces in variable density incompressible flows, J. comput. phys., 130, 269, (1997) · Zbl 0872.76065 |

[23] | D. Gao, Numerical simulation of surface wave dynamics of liquid metal MHD flow on an inclined plane in a magnetic field with spatial variation, Ph.D. Dissertation, Mechanical Engineering Department, University of California Los Angeles, 2003 |

[24] | Popinet, S.; Zaleski, S., A front-tracking algorithm for accurate representation of surface tension, Int. J. numer. mech. fluids, 30, 775, (1999) · Zbl 0940.76047 |

[25] | Scardovelli, R.; Zaleski, S., Direct numerical simulation of free-surface and interfacial flow, Annu. rev. fluid mech., 31, 567, (1999) |

[26] | Miyara, A., Numerical analysis of a falling liquid film with interfacial waves on an inclined plate, Heat transfer - Asian res., 29, 233, (2000) |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.