In this chapter we describe the numerical methods that can be used to solve the ZNG equations in the context of a dynamical simulation. These equations consist of a generalized GP equation (3.21) for the condensate and a Boltzmann equation (3.42) for the thermal component. The fact that the two equations are coupled makes their numerical solution more complex than when either is considered on its own. Indeed, the distinct quantum and classical aspects of the problem require specifically tailored numerical methods. Although most of these methods are well established and described elsewhere (Taha and Ablowitz, 1984; Sanz-Serna and Calvo, 1994), we provide in this chapter a detailed pedagogical discussion that will serve as a guide to those interested in carrying out such calculations for trapped Bose gases. This chapter is based on the papers of Jackson and Zaremba (2002a,b).

There are two main parts to the numerical problem. The first is developing a method for solving the time-dependent GP equation for an arbitrary three-dimensional geometry. This we take up in Section 11.1. Second, a method is needed for solving the Boltzmann equation that accounts for the dynamics of the thermal component. Here one must deal both with the Hamiltonian dynamics of the thermal atoms, as they move in the self-consistent mean field of the condensate and thermal cloud, and with the collisions that take place between the thermal atoms themselves (the C22 collisions) and between the thermal atoms and the condensate (the C12 collisions). The methods used to account for these two distinct collisional processes are taken up in Section 11.3. As we shall see, collisions play an important role and cannot be neglected even when the dynamical behaviour is dominated by mean-field interactions.