# LBM- Nanofluid

 Register Blogs Members List Search Today's Posts Mark Forums Read April 3, 2017, 02:05 LBM- Nanofluid #1 New Member   Arun Sathiyamoorthi Join Date: Apr 2017 Posts: 3 Rep Power: 8 Hi everyone, I am working Lattice boltzmann method. Recently i faced some problems in that. currently im working in nanofluids. i have calculated all the properties using the respected equations and calculated Rayleigh number and Prandtl number for nanofluids. I have used these two in my coding but im not getting the correct ans. Someone kindly help me. double rhonf = ((1-phi)*rhof)+(phi*rhop); double betanf = (rhobetanf/rhonf); double munf = muf/pow((1-phi),2.5); double ks = ( kp+2*kf-2*phi*(kf-kp))/(kp+2*kf+phi*(kf-kp)); double knf = ks*kf; double Prnf = (munf*cpnf) /(knf); double Ranf = (rhobetanf*kf*rhocpnf*muf*Ra)/ (rhof*betaf*knf*rhof*Cpf*munf); double visco = NY*Ma*(1/sqrt(3)*(sqrt(Prnf/Ranf))); double rho0 = 6.0; double alpha = visco/Prnf; double omega = 1.0/((3.0*visco)+0.5); double omegat =1.0/((3.0*alpha)+0.5); double gbeta = Ranf*visco*visco/(float(NX*NX*NX)*Prnf); double rhocpnf = ((1-phi)*(rhof*Cpf))+(phi*rhop*Cpp); double cpnf = (rhocpnf/rhonf); double rhobetanf = ((1-phi)*(rhof*betaf))+(phi*rhop*betap); double betanf = (rhobetanf/rhonf); double munf = muf/pow((1-phi),2.5); double ks = ( kp+2*kf-2*phi*(kf-kp))/(kp+2*kf+phi*(kf-kp)); double knf = ks*kf; double Prnf = (munf*cpnf) /(knf); double Ranf = (rhobetanf*kf*rhocpnf*muf*Ra)/ (rhof*betaf*knf*rhof*Cpf*munf); double visco = NY*Ma*(1/sqrt(3)*(sqrt(Prnf/Ranf))); double rho0 = 6.0; double alpha = visco/Prnf; double omega = 1.0/((3.0*visco)+0.5); double omegat =1.0/((3.0*alpha)+0.5); double gbeta = Ranf*visco*visco/(float(NX*NX*NX)*Prnf);   November 27, 2018, 09:47 #2
New Member

Med oumza
Join Date: Nov 2015
Posts: 4
Rep Power: 9 Quote:
 Originally Posted by @arun Hi everyone, I am working Lattice boltzmann method. Recently i faced some problems in that. currently im working in nanofluids. i have calculated all the properties using the respected equations and calculated Rayleigh number and Prandtl number for nanofluids. I have used these two in my coding but im not getting the correct ans. Someone kindly help me. double rhonf = ((1-phi)*rhof)+(phi*rhop); double betanf = (rhobetanf/rhonf); double munf = muf/pow((1-phi),2.5); double ks = ( kp+2*kf-2*phi*(kf-kp))/(kp+2*kf+phi*(kf-kp)); double knf = ks*kf; double Prnf = (munf*cpnf) /(knf); double Ranf = (rhobetanf*kf*rhocpnf*muf*Ra)/ (rhof*betaf*knf*rhof*Cpf*munf); double visco = NY*Ma*(1/sqrt(3)*(sqrt(Prnf/Ranf))); double rho0 = 6.0; double alpha = visco/Prnf; double omega = 1.0/((3.0*visco)+0.5); double omegat =1.0/((3.0*alpha)+0.5); double gbeta = Ranf*visco*visco/(float(NX*NX*NX)*Prnf); double rhocpnf = ((1-phi)*(rhof*Cpf))+(phi*rhop*Cpp); double cpnf = (rhocpnf/rhonf); double rhobetanf = ((1-phi)*(rhof*betaf))+(phi*rhop*betap); double betanf = (rhobetanf/rhonf); double munf = muf/pow((1-phi),2.5); double ks = ( kp+2*kf-2*phi*(kf-kp))/(kp+2*kf+phi*(kf-kp)); double knf = ks*kf; double Prnf = (munf*cpnf) /(knf); double Ranf = (rhobetanf*kf*rhocpnf*muf*Ra)/ (rhof*betaf*knf*rhof*Cpf*munf); double visco = NY*Ma*(1/sqrt(3)*(sqrt(Prnf/Ranf))); double rho0 = 6.0; double alpha = visco/Prnf; double omega = 1.0/((3.0*visco)+0.5); double omegat =1.0/((3.0*alpha)+0.5); double gbeta = Ranf*visco*visco/(float(NX*NX*NX)*Prnf);
Dear @arun

I am using alos LBM to simulate the nanofluids behaviour. But, I have a problem with in physical units conversion to the lattice ones. Thank you a lot if you can help me  Thread Tools Search this Thread Show Printable Version Email this Page Search this Thread: Advanced Search Display Modes Linear Mode Switch to Hybrid Mode Switch to Threaded Mode Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules Similar Threads Thread Thread Starter Forum Replies Last Post sab10012 FLUENT 0 October 20, 2015 12:38 CFD Newbie Main CFD Forum 44 January 13, 2014 13:56 pharg_yrartibra Main CFD Forum 2 November 25, 2013 10:52 mech-engineer CFX 2 June 4, 2012 05:36 ingenieroelder CFX 1 May 10, 2012 23:46

All times are GMT -4. The time now is 09:30. 