Modeling non-ideal compressible flows in the context of computational fluid-dynamics (CFD) requires the calculation of thermodynamic state properties at each step of the iterative solution process. To this purpose, the use of a built-in fundamental equation of state (EoS) in entropic form, i.e., s= s(e, ρ), can be particularly cost-effective, as all state properties can be explicitly calculated from the conservative variables of the flow solver. This approach can be especially advantageous for massively parallel computations, in which look-up table (LuT) methods can become prohibitively expensive in terms of memory usage. The goal of this research is to: i) develop a fundamental relation based on the entropy potential; ii) create a data-driven model of entropy and its first and second-order derivatives, expressed as a function of density and internal energy; iii) test the performance of the data-driven thermodynamic model on a CFD case study. Notably, two Multi-Layer Perceptron (MLP) models are trained on a synthetic dataset comprising 500k thermodynamic state points, obtained by means of the Span-Wagner EoS. The thermodynamic properties are calculated by differentiating the fundamental equation, thus ensuring thermodynamic consistency. Conversely, thermodynamic stability is properly enforced during the regression process. Albeit the method is applicable to the development of equation of state models for arbitrary fluids and thermodynamic conditions, the present work only considers siloxane MM in the single phase region. The MLP model is implemented in the open-source SU2 software  and is used for the numerical simulation of non-ideal compressible flows in a planar converging-diverging nozzle. Finally, the accuracy and the computational performance of the data-driven thermodynamic model are assessed by comparing the resulting flow field, the wall time and the memory requirements with those obtained with direct calls to a cubic EoS, and with a LuT method.