We will consider perturbations of 2D CFTs by multiple relevant operators. The massive phases of such perturbations can be labeled by conformal boundary conditions. Cardy's variational ansatz approximates the vacuum state of the perturbed theory by a smeared conformal boundary state. In this talk we will discuss the limitations and propose generalisations of this ansatz using both analytic and numerical insights based on TCSA. In particular we analyse the stability of Cardy's ansatz states with respect to boundary relevant perturbations using bulk-boundary OPE coefficients. We show that certain transitions between the massive phases arise from a pair of boundary RG flows. The RG flows start from the conformal boundary on the transition surface and end on those that lie on the two sides of it. As an example we work out the details of the phase diagram for the Ising field theory and for the tricritical Ising model perturbed by the leading thermal and magnetic fields. Based on arXiv:2306.13719.