Abstract

A pure and incompressible material is confined between two plates such that it is heated from below and cooled from above. When its melting temperature is comprised between these two imposed temperatures, an interface separating liquid and solid phases appears. Depending on the initial conditions, freezing or melting occurs until the interface eventually converges towards a stationary state. This evolution is studied numerically in a two-dimensional configuration using a phase-field method coupled with the Navier-Stokes equations. Varying the control parameters of the model, we exhibit two types of equilibria: diffusive and convective. In the latter case, Rayleigh-Bénard convection in the liquid phase shapes the solid-liquid front, and a macroscopic topography is observed. A simple way of predicting these equilibrium positions is discussed and then compared with the numerical simulations. In some parameter regimes, we show that multiple equilibria can coexist depending on the initial conditions. We also demonstrate that, in this bi-stable regime, transitioning from the diffusive to the convective equilibrium is inherently a nonlinear mechanism involving finite amplitude perturbations.

We develop and analyse the first second-order phase-field model to combine melting and dissolution in multi-component flows. This provides a simple and accurate way to simulate challenging phase-change problems in existing codes. Phase-field models simplify computation by describing separate regions using a smoothed phase field. The phase field eliminates the need for complicated discretizations that track the moving phase boundary. However, standard phase-field models are only first-order accurate. They often incur an error proportional to the thickness of the diffuse interface. We eliminate this dominant error by developing a general framework for asymptotic analysis of diffuse-interface methods in arbitrary geometries. With this framework, we can consistently unify previous second-order phase-field models of melting and dissolution and the volume-penalty method for fluid–solid interaction. We finally validate second-order convergence of our model in two comprehensive benchmark problems using the open-source spectral code Dedalus.

No abstract

No abstract