U ovom radu spajamo numeričku metodu usmjerenu na simulaciju strujanja i prijenosa topline nanofluida sa stohastičkim modeliranjem ulaznih i materijalnih parametara. Kako bi se simulirali nanofluidi, razvijena je interna numerička metoda koja se temelji na rješenju 3D formulacije brzine i vrtloga Navier-Stokesovih jednadžbi. Za rješavanje upravljačkih jednadžbi i postavljanje determinističkog rješavača protoka i prijenosa topline korištena je brza integralna metoda granične domene. Razvijeni algoritam koristi se za simulaciju prirodne konvekcije nanofluida u zatvorenoj šupljini. Nesigurnost prisutna u ulaznim parametrima propagira se na izlaz od interesa putem metode stohastičke kolokacije. Predstavljena je stohastička sredina, varijanca i momenti višeg reda izlaznih vrijednosti. Nenametljiva priroda metode stohastičke kolokacije omogućuje da prethodno potvrđeni deterministički kod ostane nepromijenjen. Stohastička analiza otkriva da nesigurnost ulaznih parametara najviše utječe na izlazne rezultate u područjima gdje se pojavljuju visoki gradijenti polja protoka.