Numerical Methods for Protocell Simulations