Computational studies on nanofluids composed of metal organic frameworks were performed using molecular modeling techniques. Grand Canonical Monte Carlo simulations were used to study the adsorption behavior of 1,1,1,3,3-pentafluoropropane (R-245fa) in a MIL-101 metal organic frameworks at various temperatures. To understand the stability of the nanofluid composed of MIL-101 particles, we performed molecular dynamics simulations to compute potentials of mean force between hypothetical MIL-101 fragments terminated with two different kinds of modulators in R-245fa and water. Our computed potentials of mean force results indicate that the metal organic frameworks particles tend to disperse better in water than in R-245fa. The reasons for this difference in dispersion were analyzed and are discussed in the paper. Our results agree with experimental results indicating that the potential models employed and modeling approaches provide good descriptions of molecular interactions and the reliabilities.