LBM- Nanofluid

@arun · April 3, 2017, 03:05

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);

medoumza · November 27, 2018, 10:47

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

April 3, 2017, 03:05	LBM- Nanofluid	#1
@arun New Member Arun Sathiyamoorthi Join Date: Apr 2017 Posts: 3 Rep Power: 9	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)+(phirhop); double betanf = (rhobetanf/rhonf); double munf = muf/pow((1-phi),2.5); double ks = ( kp+2kf-2phi(kf-kp))/(kp+2kf+phi(kf-kp)); double knf = kskf; double Prnf = (munfcpnf) /(knf); double Ranf = (rhobetanfkfrhocpnfmufRa)/ (rhofbetafknfrhofCpfmunf); double visco = NYMa(1/sqrt(3)(sqrt(Prnf/Ranf))); double rho0 = 6.0; double alpha = visco/Prnf; double omega = 1.0/((3.0visco)+0.5); double omegat =1.0/((3.0alpha)+0.5); double gbeta = Ranfviscovisco/(float(NXNXNX)Prnf); double rhocpnf = ((1-phi)(rhofCpf))+(phirhopCpp); double cpnf = (rhocpnf/rhonf); double rhobetanf = ((1-phi)(rhofbetaf))+(phirhopbetap); double betanf = (rhobetanf/rhonf); double munf = muf/pow((1-phi),2.5); double ks = ( kp+2kf-2phi(kf-kp))/(kp+2kf+phi(kf-kp)); double knf = kskf; double Prnf = (munfcpnf) /(knf); double Ranf = (rhobetanfkfrhocpnfmufRa)/ (rhofbetafknfrhofCpfmunf); double visco = NYMa(1/sqrt(3)(sqrt(Prnf/Ranf))); double rho0 = 6.0; double alpha = visco/Prnf; double omega = 1.0/((3.0visco)+0.5); double omegat =1.0/((3.0alpha)+0.5); double gbeta = Ranfviscovisco/(float(NXNXNX)Prnf);

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Nanofluid evaporation under radiation	sab10012	FLUENT	0	October 20, 2015 13:38
PowerFlow (LBE) vs. Traditional (Navier Stokes)?	CFD Newbie	Main CFD Forum	44	January 13, 2014 14:56
Practical difference between FVM and LBM?	pharg_yrartibra	Main CFD Forum	2	November 25, 2013 11:52
Nanofluid Modeling at Microchannel Heat-Sink	mech-engineer	CFX	2	June 4, 2012 06:36
Simulate a nanofluid	ingenieroelder	CFX	1	May 11, 2012 00:46