The development of efficient industrial oxy-coal boilers can be significantly aided by Computational Fluid Dynamics (CFD) tools, as far as fidelity in modeling coal combustion is also complemented by feasible computational costs. Reduced and predictive models are the most suitable for this application scale. Reduced models feature predictivity when they are validated against a broad range of experiments and targeted by Uncertainty Quantification (UQ) procedures. This work proposes a numerical procedure that uses Bound-to-Bound Data Collaboration (B2B-DC) to derive a reduced char combustion model describing transport phenomena and reactions between char carbon and O-2, CO2 and H2O, in both conventional and oxy-conditions. The approach determines the consistency between a numerical model and an experimental dataset. The latter is made up of the experiments carried out in an optically accessible laminar entrained flow reactor, operated by Sandia National Laboratories. The procedure follows five steps towards predictive modeling capability, namely: quantification of the uncertainty in the experiments, via instrument verification and modeling; development of a physics model and continuous improvement of its fidelity, via model-form uncertainty; identification of the uncertain and most sensitive parameters and of their prior bounds; sampling of the initial uncertain parameter space and training of a surrogate model; validation of the physics model via inference from the data. The last step, also known as inverse problem, is performed by applying the Bound-to-Bound Data Collaboration approach. A char combustion model is found consistent with the experimental data and its validity stands for conventional and oxy-combustion conditions. It accounts for heterogeneous reactions at the particle surface, mass transport of species in the particle boundary layer, pore diffusion and surface area changes. The consistent reduced model overcomes the differences in mass transport and kinetics observed in the experimental campaign. A reduction of the initial degree of uncertainty in both model and experiments is achieved.