We explore heat transport across an ion Coulomb crystal beyond the harmonic regime by tuning it across the structural phase transition between the linear and zigzag configurations. This demonstrates that the control of the spatial ion distribution by varying the trapping frequencies renders ion Coulomb crystals an ideal test-bed to study heat transport properties in finite open system of tunable non-linearities. 64.60.Ht, 05.70.Fh, 37.10.Ty Ultracold ion Coulomb crystals represent one of the most promising platforms for the simulation of many-body physics thanks to the high degree of spatial and temporal control of mesoscopic ion crystals they afford us with [1][2][3]. Recent years have seen a shift away from the study of ground and thermal state properties, towards the exploration of the potential role of ion traps as a test-bed for models of nonequilibrium statistical mechanics.In this context it is important to recognize that in addition to the electronic spin degrees of freedom trapped ions also possess motional degrees of freedom that can exhibit highly non-trivial static and dynamical properties including classical and quantum phase transitions. Indeed, ion Coulomb crystals confined in ion traps may support a wide variety of phases including a linear chain and a doubly-degenerate zigzag phase, extending further to increasingly complex configurations in two and three spatial dimensions [4][5][6]. The associated structural phase transitions between those configurations are generally of first order, with the exception of the linear-to-zigzag phase transition which is known to be of second order [7][8][9]. As a symmetry breaking scenario, it provides a natural testing ground for universal dynamics of phase transitions and topological defect formation [9][10][11][12], recently explored in the laboratory [13][14][15][16].Another fundamental setting in non-equilibrium statistical mechanics considers the thermal transport in low-dimensional systems, which exhibits a rich variety of anomalous features, including the breakdown of Fourier's law of heat conduction, instances in which sub-diffusive and super-diffusive behavior can be observed [17][18][19], as well as the divergence of the thermal conductivity with the system size [20,21]. Most rigorous theoretical results have been obtained for exactly solvable quasi-free models while systems with non-linearities are typically exceedingly difficult to treat. Equally, the controlled generation of non-linear physics in mesoscopic ion crystals is non-trivial and much recent progress has concerned harmonic models of complex networks and trapped-ion chains [22][23][24]. The richest phenomenology however can be expected in nonintegrable models [25] which mandates the development of both theoretical and experimental methods for their examination.In this Letter, we advance further the case for trapped ions as a model system for the examination of challenging problems in mesoscopic physics by considering continuously driven ion chains between two thermal reservoirs as a plat...