|
[Sponsors] |
July 15, 2014, 06:19 |
Couette flow in fluent
|
#1 |
Member
Ashutosh
Join Date: Jul 2013
Posts: 98
Rep Power: 13 |
I am trying to analyse couette flow in fluent. Couette flow is basically flow fluid between two parallel plates in which one plate is moving. Due to the motion of one plate the motion is imparted to the fluid.
I am analysing in 2D. I have created rectangle. I am considering the top plate as moving with velocity of 1m/s. So, what should be my boundary conditions to left side, right side, top plate and bottom plate? |
|
July 15, 2014, 06:35 |
|
#2 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,428
Rep Power: 49 |
The "bottom" boundary should be a wall with zero (or -1 m/s, that is up to you) tangential velocity.
"Left" and "right" boundaries should be connected through a periodic boundary condition. BTW: the description of this sub-forum clearly states "do not post fluent questions here" |
|
July 15, 2014, 06:41 |
|
#3 |
Member
Ashutosh
Join Date: Jul 2013
Posts: 98
Rep Power: 13 |
Thanks flotus1 for your reply. Further I would like to know for the periodic boundary condition how should I compute the offset. Should I let be computed automatically by fluent. I actually do not understand this concept of offset.
PS:Sorry for posting question in this forum. Can you also tell how can I transfer this question the software/fluent forum.? |
|
July 15, 2014, 10:01 |
|
#4 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,428
Rep Power: 49 |
The offset is the vector from one boundary in the periodic pair to the other boundary.
Lets assume your domain for Couette flow has an extent of 1m in x-direction, then the offset vector would be (1,0)m. Usually fluent computes the vector correctly, but you should check it again. If Fluent for whatever reason does not compute the vector correctly, you can input it manually. Moving the post is not possible for you, only moderators can do that. I just wanted to let you know because there are tons of questions related to fluent in this forum. |
|
July 16, 2014, 07:22 |
|
#5 |
Member
Ashutosh
Join Date: Jul 2013
Posts: 98
Rep Power: 13 |
I applied the periodic boundary condition. But I am confused about the exact method. Is it necessary to match the periodic edges in meshing?
I get a message cannot create surface from sliding interface zone, creating empty surface. |
|
July 16, 2014, 10:45 |
|
#6 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,428
Rep Power: 49 |
As far as I know, the error message you get is not related to periodic boundaries.
Did you follow the steps outlined in the tutorial "Modeling Periodic Flow and Heat Transfer"? As far as the meshing is concerned, it is not necessary to set up periodicity here. You can do it in fluent, following the steps described in the tutorial. Of course you need to make sure that the topologies of the two faces you want to make periodic are identical. And the meshes on both faces (lines in your case, since you are dealing with 2D) should be identical. Otherwise the periodic boundary condition will rather be an interface that needs lots of interpolation. |
|
July 16, 2014, 11:47 |
|
#7 |
Member
Ashutosh
Join Date: Jul 2013
Posts: 98
Rep Power: 13 |
I followed all the steps given in the tutorial of modeling periodic flow and heat transfer. My edges are totally identical. Now i don't get the error. however, solution does not converge. The residuals oscillate and then flatten out. What could be the reason?
|
|
July 16, 2014, 12:33 |
|
#8 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,428
Rep Power: 49 |
Post a figure of the residuals.
|
|
July 16, 2014, 15:47 |
|
#9 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,897
Rep Power: 73 |
first, check the steady solution for small Re number, check also to have a divergence-free velocity field
|
|
July 17, 2014, 04:10 |
|
#10 |
Member
Ashutosh
Join Date: Jul 2013
Posts: 98
Rep Power: 13 |
Well, after numerous tries I did not get convergence.
This is the image of my residuals: https://imageshack.com/i/hlb214dfj |
|
July 17, 2014, 04:30 |
|
#11 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,428
Rep Power: 49 |
Filippo might be on the right track here.
What is the Reynolds number in your simulation? Are you sure it is low enough for steady laminar flow? |
|
July 17, 2014, 09:19 |
|
#12 |
Member
Ashutosh
Join Date: Jul 2013
Posts: 98
Rep Power: 13 |
Oh, yes I overlooked that aspect. The reynolds number in my case is around 1x10^6. So, I think it must be a turbulent flow. I am actually following tutorial by Penn State University where they haven't turned the turbulent models.
Also, I observed the following: Whenever I use second order upwind for discretisation of momentum I get the case of no convergence and flat residuals. However, with first order upwind I get convergence. The velocity profile is fairly unaffected though. I have plotted my velocity plot on a centerline in my domain. Is it according to the expectations? https://imageshack.com/i/ex8c56afj |
|
July 17, 2014, 09:31 |
|
#13 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,428
Rep Power: 49 |
Pointless to discuss the results of spatially under-resolved two-dimensional DNS solved with a steady state solver.
|
|
July 17, 2014, 09:36 |
|
#14 |
Member
Ashutosh
Join Date: Jul 2013
Posts: 98
Rep Power: 13 |
So, what could be the problem? Is it the periodic boundary condition that is creating a problem?
|
|
July 17, 2014, 10:11 |
|
#15 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,428
Rep Power: 49 |
The problem is the Reynolds number.
At least to my understanding, Couette flow implies laminar flow. The flow you are modeling is turbulent. From the residuals you posted I could derive that you are not using any turbulence model. You can either lower the Reynolds number or use some proper turbulence modeling. The latter would require a finer mesh in wall-normal direction, 10 cells are definitely too few to get yplus values in the correct range. |
|
July 18, 2014, 08:03 |
|
#16 |
Member
Ashutosh
Join Date: Jul 2013
Posts: 98
Rep Power: 13 |
I reduced the Reynolds number to 5000 by reducing the velocity of the moving wall i.e. the top plate. However, I havent got any convergence. Residuals oscillate but their average behaviour is flat.
|
|
July 18, 2014, 09:33 |
|
#17 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,897
Rep Power: 73 |
you should set Re=1
|
|
July 18, 2014, 11:42 |
|
#18 |
New Member
Ali
Join Date: Jul 2014
Posts: 8
Rep Power: 12 |
http://www.mne.psu.edu/cimbala/Learn...uette_flow.pdf
http://www.mne.psu.edu/cimbala/Learn...uette_heat.pdf i don't know you may find these ones but you can use these pdf's for your solutions |
|
July 18, 2014, 11:49 |
|
#19 |
Member
Ashutosh
Join Date: Jul 2013
Posts: 98
Rep Power: 13 |
Yes Ali, I am following the same one. My dimensions and all specifications are the same as in the tutorial. However, I am having certain issues about convergence which I am currently trying to sort out.
|
|
July 18, 2014, 12:13 |
|
#20 |
New Member
Ali
Join Date: Jul 2014
Posts: 8
Rep Power: 12 |
convergence problem is limited for a specific parameter or all parameters are non-convergent? if continuity is the problem you can fix pressure relaxiation value from= solve>controls>solution. you may use for other parameters but i'm not sure if it helps cause i never tried them
|
|
Tags |
boundary condition, couette flow, fluent |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
looking for a smart interface matlab fluent | chary | FLUENT | 24 | June 18, 2021 10:07 |
The fluent stopped and errors with "Emergency: received SIGHUP signal" | yuyuxuan | FLUENT | 0 | December 3, 2013 23:56 |
Does fluent calculate reversed flow situations more simply than it should? | heksel8i | FLUENT | 0 | August 8, 2013 09:03 |
Rotating Couette Flow? | courtjester140 | Main CFD Forum | 2 | February 20, 2012 11:59 |
Is Fluent v4.5 efficient for transonic flow | Pef | FLUENT | 0 | March 8, 2002 11:43 |