A Lagrangian-Eulerian finite element algorithm for advection-diffusion-reaction problems with phase change

Beñat Oliveira*, Juan Carlos Afonso, Sergio Zlotnik

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

8 Citations (Scopus)


This paper presents a particle-based Lagrangian-Eulerian algorithm for the solution of the unsteady advection-diffusion-reaction heat transfer equation with phase change. The algorithm combines a Lagrangian formulation for the advection + reaction problem with the Eulerian-based heat source method for the diffusion + phase change problem. The coupling between the Lagrangian and Eulerian subproblems is achieved with a phase change detector scheme based on a local latent heat balance and a consistent/conservative interpolation technique between Lagrangian particles and the Eulerian grid. This technique makes use of an auxiliary (finer) Eulerian grid that provides a simple and efficient method of tracking internal heterogeneities (e.g. phase boundaries), allows the use of higher order integration quadratures, and facilitates the implementation of multiscale techniques. The performance of the proposed algorithm is compared against one- and two-dimensional benchmark problems, i.e. pure rigid-body advection, isothermal and non-isothermal phase change, two-phase advective heat transfer and chemical reactions coupled with diffusion and advection. The numerical results confirm that the proposed solution method is accurate, oscillation-free and useful for and applicable to a wide range of fully coupled problems in science and engineering.

Original languageEnglish
Pages (from-to)375-401
Number of pages27
JournalComputer Methods in Applied Mechanics and Engineering
Publication statusPublished - 1 Mar 2016


Dive into the research topics of 'A Lagrangian-Eulerian finite element algorithm for advection-diffusion-reaction problems with phase change'. Together they form a unique fingerprint.

Cite this