We studied hydration of four types of small nonpeptidic molecule, flavonoid, by molecular dynamics simulations at 300 K with focusing on three physical quantities: solvent density, solvent site-dipole field, and solvent diffusion. The solvent site-dipole field is a quantity recently introduced by us to study directional ordering of water molecules around solute. The spatial patterns of these quantities showed strong site-dependency around the flavonoids. Common to the four flavonoids, high solvent-density sites around hydrophilic solute atoms were characterized by strong directional ordering of water molecule and by depressed solvent diffusive motions. Contrarily, high solvent-density sites around hydrophobic solute surface were characterized by weak directional ordering. The solvent site-dipole field showed specific ordering patterns of water molecules not only in the first solvent layer but also in the second solvent layer. The spatial patterns of the three quantities were conservative among the four flavonoids whether the intra-flavonoid flexibility was large or not. Thus, an adiabatic approximation, which has been assumed in various theoretical hydration studies, was satisfied well. The hydration at a site in the vicinity of solute was determined mainly by the physico-chemical property of the solute atom group nearest to the solvent site, which supports a phenomenological theorem that the solvent accessible surface area of a solute is proportional to the solvation free energy.