The proposed method enables automatic generation of a volumetric mesh. It can be applied to arbitrary surfaces representations of genus 0. Isogeometric and classical finite elements are combined. The application to thermal problems leads to a substantial accuracy gain.