We study the melting process of a solid under microgravity, driven solely by lateral vibrations that are perpendicular to the applied temperature gradient due to the absence of gravity-induced convection. Using direct numerical simulations with the phase-field method, we examine two-dimensional vibration-induced melting in a square cavity over four orders of magnitude of vibrational Rayleigh numbers, $10^5\le Ra_{{vib}}\le 10^9$. Our results show that as melting progresses, the flow structure transitions from a periodic-circulation regime with diffusion-dominated heat transfer to a columnar regime with vibroconvection. The mean height of the liquid–solid interface follows a power-law dependency with time, $\bar {\xi } \sim \tilde t^{1/(2-2\alpha )}$, where $\alpha = 0$ in the periodic-circulation regime and $\alpha = 1/2$ in the columnar regime. We further observe that within the columnar regime, the morphological evolution of the liquid–solid interface is influenced by the interaction of columnar thermal plumes in the central regions and the peripheral flow near the sidewalls. Specifically, we offer a comprehensive analysis of the plume merging behaviour, which is governed by the aspect ratio ($\bar {\xi }$) of the liquid layer and the intensity of vibration, quantified by the effective vibrational Rayleigh number $Ra_{vib}^{eff}$. We identify the relationship between the number of columnar plumes $K_m$ and $Ra_{vib}^{eff}$, finding that $K_m \sim \bar {\xi }^{-1} (Ra_{vib}^{eff})^{\gamma }$ with the fitting scaling exponent $\gamma = 0.150 \pm 0.025$. We subsequently quantify the characteristics of the interface roughness amplitude evolution in microgravity vibroconvection. Our results indicate that the roughness amplitude exhibits a power-law dependence on the mean height of the liquid layer. Drawing from the Stefan boundary condition, we theoretically deduce this dependence under the assumption of a non-uniform heat flux distribution at the interface, where the theory is corroborated by our numerical simulations.

We report the unified constitutive law of vibroconvective turbulence in microgravity, i.e. $Nu \sim a^{-1} Re_{os}^\beta$ where the Nusselt number $Nu$ measures the global heat transport, $a$ is the dimensionless vibration amplitude, $Re_{os}$ is the oscillational Reynolds number and $\beta$ is the universal exponent. We find that the dynamics of boundary layers plays an essential role in vibroconvective heat transport and the $Nu$-scaling exponent $\beta$ is determined by the competition between the thermal boundary layer (TBL) and vibration-induced oscillating boundary layer (OBL). Then a physical model is proposed to explain the change of scaling exponent from $\beta =2$ in the TBL-dominant regime to $\beta = 4/3$ in the OBL-dominant regime. Our finding elucidates the emergence of universal constitutive laws in vibroconvective turbulence, and opens up a new avenue for generating a controllable effective heat transport under microgravity or even microfluidic environment in which the gravity effect is nearly absent.