A study is presented of photothermal effects associated with nanosecond-pulsed laser-illuminated subwavelength metallic nanoparticles in aqueous solutions. Computational electromagnetic and fluid analysis are used to model fundamental aspects of the photothermal process taking into account energy conversion within the nanoparticle at plasmon resonance, heat transfer to the fluid, homogeneous bubble nucleation, and the dynamic behaviour of the bubble and surrounding fluid. Various nanoparticle geometries are modelled including spheres, nanorods and tori. The analysis demonstrates that the laser intensity and pulse duration can be tuned to achieve controllable bubble generation without exceeding the melting temperature of the particle. The analysis also shows that the particle geometry can be tuned to optimize photothermal energy conversion for bubble generation at wavelengths that span the UV to NIR spectrum. Multiparticle systems are studied and a cooperative heating effect is demonstrated for particles that are within a few radii of each other. This provides more robust bubble generation using substantially reduced laser energy as compared to single-particle systems. The modelling approach is discussed in detail and should be of considerable use in the development of new photothermal applications.