U ovom radu povezujemo numeričku metodu s ciljem simulacije protoka i prijenosa topline nanofluida sa stohastičkim modeliranjem ulaznih i materijalnih parametara. Razvijena je brza integralna metoda granične domene za rješavanje upravljačkih jednadžbi i postavljanje determinističkog rješavača protoka i prijenosa topline. Nadalje, metoda stohastičke kolokacije koristi se u kombinaciji s determinističkim kodom za simulaciju toka kako bi se nesigurnost prenijela s ulaznih na izlazne parametre. Razvijeni algoritam koristi se za simulaciju prirodne konvekcije nanofluida u zatvorenoj šupljini. Simulacijom vrlo velikog broja slučajeva i primjenom stohastičke analize, uspjeli smo identificirati relativni utjecaj različitih ulaznih parametara. Rezultati pokazuju da nesigurnost ulaznih parametara rezultira snažnijim utjecajem u režimima strujanja u kojima dominira konvekcija.