A major challenge in the combustion of biomass fuels is the heterogeneity of ash-forming elements, which may cause a wide range of ash-related problems. Understanding the melting tendency of the coarse ash fractions is necessary to mitigate agglomeration and slagging. This work aims to evaluate the melting tendency of the K-Ca-Mg-Si-P-O system by use of thermodynamic equilibrium calculations. The formation of condensed phases were systematically assessed in a combustion atmosphere, varying temperatures, and composition. Compositional ranges were based on fuel ash data extracted from the Phyllis 2 database. The speciation and degree of polymerization of phosphates, silicates, and melts were evaluated and indicated a systematic variation in composition. The melt fraction was predicted as a function of temperature and composition. The melting tendency was modeled for three systems, i.e., a P-dominated, a Si-dominated, and a mixed Si-P system. Four ratios between K 2 O, CaO, MgO, SiO 2 , and P 2 O 5 were found to have a large effect on the melting tendency of the ash mixtures: the ratio between network formers (SiO 2 , P 2 O 5 ), K 2 O to total network modifiers, CaO to CaO + MgO, and the ratio of network formers to total ash oxides. This modeling approach showed qualitative agreement with ash-related issues seen in previous lab-scale experiments in bubbling fluidized bed and fixed bed combustion. Practical implications of the results are discussed from the perspective of fuel design with the aim of preventing ashrelated problems. This study presents a novel method of applying thermodynamic equilibrium calculations for a broad range of compositions and shows potential for predicting ash-related issues related to the melting of coarse ash fractions.