I. Introduction Due to their high torque density and an excellent efficiency, axial flux PM machines are favorable for vehicle propulsion and wind energy conversion. Although many effort was already done to simulate the electromagnetic properties of the machine using full 3D or multilayer-2D simulations, the work towards modelling of the thermal behavior is still very limited. As the disc-shaped rotors will operate as a fan during rotation, convective heat transfer in the air gap will have a major influence on the thermal design of the machine. Rather than performing full 3D combined fluid and thermal analysis, empiric analytical expressions based on CFD are used in this work to define the convective heat coefficient in different parts of the machine. The final thermal simulations are carried out for the stator and rotor individually, where only a segment needs to be modelled due to thermal periodicity. As a result, the use of this coupled modelling technique results in a significant decrease of the overall simulation time. The coupled electromagnetic and thermal modelling techniques are illustrated on a 4kW axial flux PM machine having the yokeless and segmented armature (YASA) topology. Finally, both the electromagnetic and thermal simulation results were validated on a preliminary prototype. II. Coupled Electromagnetic and Thermal Modelling In the first step of the analysis, the electromagnetic behavior of the machine is modelled. The results from the electromagnetic field computations are used to calculate the corresponding losses in post-processing. In a second step, these losses are used as the heat sources for the thermal simulations. As the axial flux PM machine has an inherent 3D structure, multilayer-2D simulations [1] are used to obtain the electromagnetic properties such as flux-linkage, back-emf and (cogging) torque. In postprocessing, the magnetic flux density pattern of the different layers in the core is used to calculated the core losses. The solutions for the air gap magnetic flux density pattern of each layer are combined and used for a time harmonic calculation of the eddy current losses in the permanent magnets [1]. Together with the Joule losses in the machine winding, the iron losses in the stator cores and eddy current losses in the PM's are used as the source terms in the thermal simulations. As the thermal process is particularly slow with respect to the electromagnetic one, the time average values of the losses are chosen as a source in the thermal model. Instead of using full 3D combined fluid and thermal analysis, empirical equations derived in [2] are used to model the convective heat transfer from each of the surfaces near the air gap region. Next to an empirical formula for the convective heat coefficients of each surface, an expression for the reference temperatures is also derived using dimensionless numbers. The reference temperature is the bulk fluid temperature of the domain near each surface which can be expressed as a function of the average stator, rotor and ambient temper...