A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation


Autoria(s): CENICEROS, Hector D.; ROMA, Alexandre M.; SILVEIRA-NETO, Aristeu; VILLAR, Millena M.
Contribuinte(s)

UNIVERSIDADE DE SÃO PAULO

Data(s)

20/10/2012

20/10/2012

2010

Resumo

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

National Science Foundation (NSF)

National Science Foundation (NSF)[DMS 0609996]

Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP)[04/13781-1]

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP)[06/57099-5]

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq)[155491/2006-7]

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Identificador

COMMUNICATIONS IN COMPUTATIONAL PHYSICS, v.8, n.1, p.51-94, 2010

1815-2406

http://producao.usp.br/handle/BDPI/30578

10.4208/cicp.050509.141009a

http://dx.doi.org/10.4208/cicp.050509.141009a

Idioma(s)

eng

Publicador

GLOBAL SCIENCE PRESS

Relação

Communications in Computational Physics

Direitos

restrictedAccess

Copyright GLOBAL SCIENCE PRESS

Palavras-Chave #Front-tracking #immersed boundary method #level set method #adaptive mesh refinements #semi-implicit methods #multilevel multigrid #closest point transform #semi-backward difference formula #RAYLEIGH-TAYLOR INSTABILITY #IMMERSED BOUNDARY METHOD #NAVIER-STOKES EQUATIONS #MULTI-FLUID FLOWS #FRONT-TRACKING #PROJECTION METHOD #NUMERICAL SIMULATIONS #INTERFACE TRACKING #DROP DEFORMATION #3 DIMENSIONS #Physics, Mathematical
Tipo

article

original article

publishedVersion