Numerical simulation of conduit dynamics during paroxysms at Stromboli