Extended functionality of model: introduced elastic boundary conditions, a box inside the box and sample velocities from random normal distribution