A molecular dynamics model to measure forces between cellulose fibril surfaces: on the effect of non-covalent polyelectrolyte adsorption