Organische Halbleiterbauelemente sind eine vielversprechende Technologie, die das Spektrum der optoelektronischen Halbleiterbauelemente erweitert und etablierte Technologien basierend auf anorganischen Halbleitermaterialien ersetzen kann. Für Display- und Beleuchtungsanwendungen werden sie z. B. als organische Leuchtdioden oder Transistoren verwendet. Eine entscheidende Eigenschaft organischer Halbleitermaterialien ist, dass die Ladungstransporteigenschaften stark von der Temperatur im Bauelement beeinflusst werden. Insbesondere nimmt die elektrische Leitfähigkeit mit der Temperatur zu, so dass Selbsterhitzungseffekte, einen großen Einfluss auf die Leistung der Bauelemente haben. Mit steigender Temperatur nimmt die elektrische Leitfähigkeit zu, was wiederum zu größeren Strömen führt. Dies führt jedoch zu noch höheren Temperaturen aufgrund von Joulescher Wärme oder Rekombinationswärme. Eine positive Rückkopplung liegt vor. Im schlimmsten Fall führt dieses Verhalten zum thermischen Durchgehen und zur Zerstörung des Bauteils. Aber auch ohne thermisches Durchgehen führen Selbsterhitzungseffekte zu interessanten nichtlinearen Phänomenen in organischen Bauelementen, wie z. B. die S-förmige Beziehung zwischen Strom und Spannung. In Regionen mit negativem differentiellen Widerstand führt eine Verringerung der Spannung über dem Bauelement zu einem Anstieg des Stroms durch das Bauelement. Diese Arbeit soll einen Beitrag zur mathematischen Modellierung, Analysis und numerischen Simulation von organischen Bauteilen leisten. Insbesondere wird das komplizierte Zusammenspiel zwischen dem Fluss von Ladungsträgern (Elektronen und Löchern) und Wärme diskutiert. Die zugrundeliegenden Modellgleichungen sind Thermistor- und Energie-Drift-Diffusion-Systeme. Die numerische Diskretisierung mit robusten hybriden Finite-Elemente-/Finite-Volumen-Methoden und Pfadverfolgungstechniken zur Erfassung der in Experimenten beobachteten S-förmigen Strom-Spannungs-Charakteristiken wird vorgestellt. / Organic semiconductor devices are a promising technology to extend the range of optoelectronic semiconductor devices and to some extent replace established technologies based on inorganic semiconductor materials. For display and lighting applications, they are used as organic light-emitting diodes (OLEDs) or transistors. One crucial property of organic semiconductor materials is that charge-transport properties are heavily influenced by the temperature in the device. In particular, the electrical conductivity increases with temperature, such that self-heating effects caused by the high electric fields and strong recombination have a potent impact on the performance of devices. With increasing temperature, the electrical conductivity rises, which in turn leads to larger currents. This, however, results in even higher temperatures due to Joule or recombination heat, leading to a feedback loop. In the worst case, this loop leads to thermal runaway and the complete destruction of the device. However, even without thermal runaway, self-heating effects give rise to interesting nonlinear phenomena in organic devices, like the S-shaped relation between current and voltage resulting in regions where a decrease in voltage across the device results in an increase in current through it, commonly denoted as regions of negative differential resistance. This thesis aims to contribute to the mathematical modeling, analysis, and numerical simulation of organic semiconductor devices. In particular, the complicated interplay between the flow of charge carriers (electrons and holes) and heat is discussed. The underlying model equations are of thermistor and energy-drift-diffusion type. Moreover, the numerical approximation using robust hybrid finite-element/finite-volume methods and path-following techniques for capturing the S-shaped current-voltage characteristics observed in experiments are discussed.