Aerodynamics in a naturally ventilated multi-span greenhouse with plants was analyzed numerically by the computational fluid dynamics (CFD) simulation. To investigate the potential application of CFD techniques to greenhouse design and analysis, the numerical results of the CFD model were compared with the results of a steady-state mass and energy balance numerical model. Assuming the results of the mass and energy balance model as the standard, reasonably good agreement was obtained between the natural ventilation rates computed by the CFD numerical model and the mass and energy balance model. The steady-state CFD model during a sunny day showed negative errors as high as 15% in the morning and comparable positive errors in the afternoon. Such errors assumed to be due to heat storage in the floor, benches, and greenhouse structure. For a west wind of 2.5 m s$^{-1}$ , the internal nonporous shading screens that opened to the east were predicted to have a 15.6% better air exchange rate than opened to the west. It was generally predicted that the presence of nonporous internal shading screens significantly reduced natural ventilation if the horizontal opening of the screen for each span was smaller that the effective roof vent opening.