Each mixing valve is characterized by its kvs-value and the valve characteristics. In mounted condition, the valve characteristics are deformed and the degree of the deformation is determined by the valve authority. Considering already installed mixing valves, neither the valve authority nor the course of the undeformed valve characteristics are exactly known. The lack of information poses a challenge when physically modelling the relationship between the inlet temperatures, the control signal and the outlet temperature. Based on Bezier curves, a universal non-linear grey-box model to describe the thermal valve behaviour is proposed. It is compared with a state space as well as a Hammerstein model and shows the advantages over black-box approaches. The trained model is very robust, able to predict measurement data very accurately and can be used in control design. The parameter estimation is performed using the Open Source platform JModelica.org.