Gold nanostructures capped with citrate and solvated in water were simulated using RNEMD methods. Two kinds of simulations were carried out: One in periodic boundary conditions using slabs of fcc gold with the (111) surface exposed to the solvent, and one using spherical gold nanoparticles fully encased by a spherical shell of the solvent. Spherical nanoparticles were simulated using the non-periodic Langevin Hull algorithm in OpenMD to maintain an external pressure and temperature bath. In both cases, RNEMD can effectively maintain a thermal flux either between two slabs (in periodic boundary conditions) or two concentric spheres (for the spherical nanoparticles).