Study of Europa's crust via a Stefan model with convection

Here we adopt the stream-function/vorticity formulation of the Navier-Stokes equations for the flow of the mushy ice and a Stefan condition combined with a front-fixing technique for the front evolution. Our numerical discretization is based upon an ENO scheme. Mathematical model and numerical procedure have been thoroughly tested and have the advantage of yielding accurate numerical solutions via relatively coarse space discretization grids. For applications in this field the present one is the first attempt, at our knowledge, to solve a complete Stefan condition with convection flow, obtaining a good match with other numerical solutions in the literature.