Large eddy simulation (LES) was performed for a spatially developing round jet and its scalar transport at four steps of Reynolds number set between 1200 and 1,000,000. A simulated domain, which extends 30 times the nozzle diameter, includes initial, transitional, and established stage of jet. A modified version of convection outflow condition was proposed in order to diminish the effect of a downstream boundary. Tested were two kinds of subgrid scale (SOS) models: a Smagorinsky model (SM) and a dynamic Smagorinsky model (DSM). In the former model, parameters are kept at empirically deduced constants, while in the latter, they are calculated using different levels of space filtering. Data analysis based on the decay law of jet clearly presented the performance of SGS models. Simulated results by SM and DSM compared favorably with existing measurements of jet and its scalar transport. However, the quantitative accuracy of DSM was better than that of SM at a transitional stage of flow field. Computed parameters by DSM, coefficient for SGS stresses, C R , and SGS eddy diffusivity ratio, Γ SGS , were not far from empirical constants of SM.Optimization of the model coefficient was suggested in DSM so that coefficient C R was nearly equal in the established stage of jet but it was reduced in low turbulence close to the jet nozzle.