El desarrollo de la tecnología electrónica hace necesario disponer de modelos fiables de los dispositivos y convierte la simulación de los mismos en un elemento estratégico para su desarrollo presente y futuro. El transistor bipolar de heterounión (HBT) es uno de los componentes más prometedores dentro del campo de los dispositivos de alta velocidad, ya que es capaz de controlar densidades de corriente elevadas a frecuencias muy elevadas. Dentro de esta familia destacan los HBTs de InP/InGaAs, que registran actualmente la frecuencia de operación más alta. Estos dispositivos exhiben una discontinuidad abrupta en las bandas de energía (en la interfaz base-emisor), cuya presencia invalida los modelos convencionales de transporte, basados en los mecanismos de arrastre-difusión (DD) o en el modelo hidrodinámico (HD) y ampliados con la emisión termoiónica y la transmisión túnel en la heterounión. Tanto el modelo DD como el HD derivan de aproximaciones de la ecuación de transporte de Boltzmann (BTE) que, alternativamente, puede ser resuelta de forma exacta por el método de Monte Carlo (MC). De esta forma, se superan las limitaciones intrínsecas de los simuladores convencionales y es posible determinar su margen de validez. El inconveniente principal de los simuladores MC es su elevado requerimiento de tiempo de cálculo. El objetivo de la tesis es desarrollar un simulador MC para HBTs abruptos que sea computacionalmente eficiente, por aplicación de técnicas de supercomputación paralela. La memoria de la tesis consta de seis capítulos. En el primero se introducen la problemática de la simulación de dispositivos, los HBTs y el método de MC aplicado a este campo. En el segundo se presenta el método convencional de simulación de HBTs, puesto que su resultado será utilizado como aproximación inicial en el simulador MC desarrollado y como referencia comparativa. El tercer capítulo se dedica a presentar el esquema básico de simulación MC de transistores bipolares, mientras que en el cuarto se dan detalles de la implementación eficiente del simulador MCHBT, con la aplicación de técnicas de supercomputación paralela y la utilización de técnicas específicas como el MC ponderado. El capítulo quinto contiene los resultados de simulación de BJTs de GaAs, para resaltar las diferencias del nuevo método de simulación con el convencional sobre un dispositivo y un material conocidos. Finalmente, se amplía el simulador MC con la resolución de la ecuación de Schrodinger en la interfaz abrupta, para incorporar los fenómenos cuánticos de transporte en el HBT. Asimismo, se utiliza la ecuación de balance del momento para interpretar los resultados y establecer un margen de validez de los simuladores convencionales ampliados para el estudio de HBTs abruptos.