A fully quantum-mechanical model for carrier scattering transport in semiconductor intersubband devices was applied to modelling of carrier dynamics in quantum cascade lasers. The standard model uses the envelope function and effective mass approximations to solve electron band structure under an applied bias. The k(.)p model has been employed in p-type systems where the more complex band structure requires it. The resulting wave functions are then used to evaluate all relevant carrier-phonon, carrier-carrier and alloy scattering rates from each quantised state to all others within the same and the neighbouring period. This piece of information is then used to construct a rate equation for the equilibrium carrier density in each sub-band and this set of coupled rate equations are solved self-consistently to obtain the carrier density in each eigenstate. The latter is a fundamental description of the device and can be used to calculate the current density and gain as a function of the applied bias and temperature, which in turn yields the threshold current and expected temperature dependence of the device characteristics