OPEN ACCESS
A coupled Conservative Level Set – Moving Mesh – Immersed Boundary method is formulated and validated against the three-dimensional gravity-driven falling drop problem. First, by employing Conservative Level-Set (CLS) method, the multiphase domain can be successfully handled, while the mass conservation is controlled. Then, by using an Arbitrary Lagrangian-Eulerian formulation (i.e. a moving mesh), the simulation domain can be optimized by reducing the domain size and by allowing an improved mesh, resulting in a computational resources saving. Finally, the use of an Immersed Boundary (IB) method allows to deal with intricate geometries. All these functionalities result in a versatile and robustness method to simulate bubbles/drops problems in complex geometries. The mentioned method was successfully used to thoroughly study the falling of a drop against a plane surface, providing detailed results including velocity evolution, mesh independence study, evolution of the vertical position of the drop, streamlines and vorticity fields, and profiles evolution.
arbitrary Lagrangian-Eulerian, complex geometries, falling drop, finite volume method, fluid-structure interaction, immersed boundary, level set, multiphase flow, open boundary condition, unstructured mesh
[1] Crowe, C., Sommerfeld, M. & Tsuji, Y., Multiphase flows with droplets and particles, CRC Press: Florida, 1997.
[2] Becher, H. & Burns, P.N., Handbook of contrast echocardiography: left ventricular function and myocardial perfusion, Springer: New York, 2000.
[3] Clift, R., Bubbles, drops and particles, Dover Publications: New York, 1978.
[4] O’Reilly Meehan, R., Donnelly, B., Nolan, K., Persoonsa, T. & Murray, D.B., Flow structures and dynamics in the wakes of sliding bubbles. International Journal of Multiphase Flow, 84, pp. 145–154, 2016. https://doi.org/10.1016/j.ijmultiphaseflow.2016.03.010
[5] Podvin, B., Khoja, S., Moraga F. & Attinger, D., Model and experimental visualizations of the interaction of a bubble with an inclined wall. Chemical Engineering Science, 63, pp. 1914–1928, 2008. https://doi.org/10.1016/j.ces.2007.12.023
[6] Han, J. & Tryggvason, G., Secondary breakup of axisymmetric liquid drops. I. Acceleration by a constant body force. Physics of Fluids, 11(12), pp. 3650–3667, 1999. https://doi.org/10.1063/1.870229
[7] Balcázar, N., Jofre, L., Lehmkuhl, O., Castro, J. & Rigola, J., A finite-volume/level-set method for simulating two-phase flows on unstructured grids. International Journal of Multiphase Flow, 64, pp. 55–72, 2014. https://doi.org/10.1016/j.ijmultiphaseflow.2014.04.008
[8] Balcázar, N., Lehmkuhl, O., Jofre, L. & Oliva, A., Level-set simulations of buoyancydriven motion of single and multiple bubbles. International Journal of Heat Fluid Flow, 56, pp. 91–107, 2015. https://doi.org/10.1016/j.ijheatfluidflow.2015.07.004
[9] Balcázar, N., Rigola, J., Castro, J. & Oliva, A., A level-set model for thermocapillary motion of deformable fluid particles. International Journal of Heat Fluid Flow, 62, pp. 324–343, 2016. https://doi.org/10.1016/j.ijheatfluidflow.2016.09.015
[10] Gutiérrez, E., Balcázar N., Bartrons, E. & Rigola, J. Numerical study of Taylor bubbles rising in a stagnant liquid using a level-set/moving-mesh method. Chemical Engineering Science, 164, pp. 158–177, 2017. https://doi.org/10.1016/j.ces.2017.02.018
[11] Peskin, C.S., Numerical analysis of blood flow in the heart. Journal of Computational Physics, 25(3), pp. 220–252, 1977. https://doi.org/10.1016/0021-9991(77)90100-0
[12] Sun, X. & Sakai, M., Numerical simulation of two-phase flows in complex geometries by using the volume-of-fluid/immersed-boundary method. Chemical Engineering Science, 139, pp. 221–240, 2016. https://doi.org/10.1016/j.ces.2015.09.031
[13] Deen, N.G., Annaland, M. van S. & Kuipers, J.A.M., Direct numerical simulation of complex multi-fluid flows using a combined front tracking and immersed boundary method. Chemical Engineering Science, 64(9), pp. 2186–2201, 2009. https://doi.org/10.1016/j.ces.2009.01.029
[14] Harten, A., The artificial compression method for computation of shocks and contact discontinuities. III Self-adjusting hybrid schemes. Mathematics of Computation, 32, pp. 363–389, 1978. https://doi.org/10.2307/2006149
[15] Brackbill, J.U., Kothe, D.B. & Zemach, C., A continuum method for modeling surface tension. Journal of Computarional Physics, 100, pp. 335–354, 1992. https://doi.org/10.1016/0021-9991(92)90240-y
[16] Guermond J.L., Minev, P. & Shen, J., An overview of projection methods for incompressible flows. Computer Methods in Applied Mechanics and Engineering, 195, pp. 6011–6045, 2006. https://doi.org/10.1016/j.cma.2005.10.010
[17] Estruch, O., Lehmkuhl, O., Borrell, R., Pérez Segarra, C.D. & Oliva, A., A parallel radial basis function interpolation method for unstructured dynamic meshes. Computers & Fluids, 80, pp. 44–54, 2013. https://doi.org/10.1016/j.compfluid.2012.06.015
[18] Sussman, M., Smereka, P. & Osher, S., A level set method for computing solutions to incompressible two phase flow. Journal of Computational Physics, 119, pp. 146–159, 1994. https://doi.org/10.1006/jcph.1994.1155
[19] Favre, F., Colomer, G., Lehmkuhl, O. & Oliva, A., Numerical simulations of conjugate convection combined with surface thermal radiation using an Immersed-Boundary Method. Journal of Physics: Conference Series, 745, p. 032017, 2016. https://doi.org/10.1088/1742-6596/745/3/032017
[20] Muradoglu, M. & Kayaalp, A.D., An auxiliary grid method for computations of multiphase flows in complex geometries. Journal of Computational Physics, 214, pp. 858–877, 2006. https://doi.org/10.1016/j.jcp.2005.10.024
[21] Harmathy, T.Z., Velocity of large drops and bubbles in media of infinite or restricted extent. AIChE Journal, 6, pp. 281–288, 1960. https://doi.org/10.1002/aic.690060222