A moving mesh method for non-isothermal multiphase flows
Repository URI
Repository DOI
Change log
Authors
Abstract
In this thesis, a numerical method is developed for simulating non-isothermal multiphase flows, which are important in many technical applications such as crystal growth and welding. The method is based on the arbitrary Lagrangian Eulerian method of Li (2013). The interface is represented explicitly by mesh lines, and is tracked by an adaptive moving unstructured mesh. The
Firstly, a thorough study is presented on the method’s capability in numerically representing the force balance condition on the interface. An inaccurate representation of this condition induces the non-physical spurious currents, which degrade the simulation accuracy especially when the viscous damping is weak (small Ohnesorge number,
Secondly, a numerical treatment of interface topology changes is incorporated into our method for studying problems with interface breakup. Thanks to the adaptive mesh generator, the thin region between the interface boundary and another boundary consists of one layer of elements. The interface topology change is performed once the minimum distance between the two boundaries falls below a pre-set scale
Finally, an FEM solver for temperature is developed and the non-isothermal effects are included in our method for the purpose of simulating non-isothermal multiphase flows. The modified method is validated to be accurate through three benchmark examples: natural convection in a cavity, thermocapillary convection of two layers, and droplet migration subject to a temperature gradient. Our method is then applied to investigate the liquid bridge breakup with thermocapillary effect. The non-isothermal liquid bridge breakup in the viscous and inertial regimes are studied. It has been found that the inertial regime breakup exhibits different pinchoff shapes as the Capillary number increases, and that the viscous regime breakup is accelerated by the thermocapillary motion.