Kinetic modelling of parallel plasma and heat flows in the scrape-off-layer (SOL) is presented. Detached [1-4] and attached regimes are studied. A kinetic Edge Localized Mode (ELM) [5] model is simulated by temporally varying the incoming heat pourer from the core. Deviations of the plasma distribution functions from Maxwellian are observed. Our simulations are performed using the collisional PIC code W1 [6], which employs a non-stationary 1D2V fully kinetic model of a high recycling SOL plasma.