# Lattice Boltzmann method_Thermal Poiseuille flow

 Register Blogs Members List Search Today's Posts Mark Forums Read September 9, 2017, 13:43 Lattice Boltzmann method_Thermal Poiseuille flow #1 Member   shikha Join Date: Aug 2013 Posts: 40 Rep Power: 11 Hello Friends, I am working on a Thermal Poiseuille flow problem using lattice Boltzmann method. The problem statement is the top and bottom boundaries are assigned with hot temperature and left and right open boundaries are set with zero (non-dimensional) BC. Here is a sample code: for(i=0;i<=nx;i++) { g[i][ny]=(t_t*(w_f+w_f))-g[i][ny]; //Dirichlet boundary condition g[i][ny]=(t_t*(w_f+w_f))-g[i][ny]; g[i][ny]=(t_t*(w_f+w_f))-g[i][ny]; g[i]=(t_b*(w_f+w_f))-g[i]; g[i]=(t_b*(w_f+w_f))-g[i]; g[i]=(t_b*(w_f+w_f))-g[i]; } /*---------------------left and right boundary condition----------------*/ for(j=0;j<=ny;j++) { g[j]=-g[j]; g[j]=-g[j]; g[j]=-g[j]; g[nx][j]=2*g[nx-1][j]-g[nx-2][j]; g[nx][j]=2*g[nx-1][j]-g[nx-2][j]; g[nx][j]=2*g[nx-1][j]-g[nx-2][j]; However, I am getting the correct result and unable to understand where I am going wrong. Please if anyone can enlighten me. Thanks in advance. Shikha   September 9, 2017, 13:47 #2 Member   shikha Join Date: Aug 2013 Posts: 40 Rep Power: 11 I am not getting the results   September 9, 2017, 18:39 #3
Senior Member

Lane Carasik
Join Date: Aug 2014
Posts: 692
Rep Power: 13 Quote:
 Originally Posted by SHIKHA BHUYAN Hello Friends, I am working on a Thermal Poiseuille flow problem using lattice Boltzmann method. The problem statement is the top and bottom boundaries are assigned with hot temperature and left and right open boundaries are set with zero (non-dimensional) BC. Here is a sample code: for(i=0;i<=nx;i++) { g[i][ny]=(t_t*(w_f+w_f))-g[i][ny]; //Dirichlet boundary condition g[i][ny]=(t_t*(w_f+w_f))-g[i][ny]; g[i][ny]=(t_t*(w_f+w_f))-g[i][ny]; g[i]=(t_b*(w_f+w_f))-g[i]; g[i]=(t_b*(w_f+w_f))-g[i]; g[i]=(t_b*(w_f+w_f))-g[i]; } /*---------------------left and right boundary condition----------------*/ for(j=0;j<=ny;j++) { g[j]=-g[j]; g[j]=-g[j]; g[j]=-g[j]; g[nx][j]=2*g[nx-1][j]-g[nx-2][j]; g[nx][j]=2*g[nx-1][j]-g[nx-2][j]; g[nx][j]=2*g[nx-1][j]-g[nx-2][j]; However, I am getting the correct result and unable to understand where I am going wrong. Please if anyone can enlighten me. Thanks in advance. Shikha
What exactly is your issue? You said you are getting the correct result.   September 10, 2017, 06:24 #4 Member   shikha Join Date: Aug 2013 Posts: 40 Rep Power: 11 I am not getting the result: The temperature profile is not correct   September 10, 2017, 23:02 #5
Senior Member

Lane Carasik
Join Date: Aug 2014
Posts: 692
Rep Power: 13 Quote:
 Originally Posted by SHIKHA BHUYAN I am not getting the result: The temperature profile is not correct
How is it not correct? You're not providing enough detail.  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 Jing Main CFD Forum 8 October 5, 2018 17:02 ashbab OpenFOAM Pre-Processing 1 June 27, 2015 07:01 DavidOuellette OpenFOAM Running, Solving & CFD 0 February 19, 2013 22:34 Joseph Bergevin Main CFD Forum 1 January 25, 2006 14:56 Rosie FLUENT 1 December 6, 2002 16:52

All times are GMT -4. The time now is 21:01.