Numerical modeling of crystal growth under strong magnetic fields: An application to the travelling heater method