CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Couette Flow in Matlab using Pressure Correction Method

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree2Likes
  • 2 Post By mulchah

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 19, 2014, 06:22
Default Couette Flow in Matlab using Pressure Correction Method
  #1
New Member
 
Jack
Join Date: Dec 2013
Posts: 6
Rep Power: 12
mulchah is on a distinguished road
Hi!

I am having some trouble completing my Matlab code for Couette flow (parallel flow between two plates, the top one moving at a specified velocity) using the pressure correction method, as described in section 9 of John Anderson's book "CFD - Basics with Applications".

The trouble is that while my code manages to complete the suggested 300 iterations to solve for the velocity and pressure fields, the mass source term I get at the end of the iterations is d = -0.350 x 10^-5 slug/ft^3.s, different from Anderson's value of d = -0.172 x 10^-5 slug/ft^3.s. Related to this, the horizontal component of velocity, u, looks slightly different after 4 iterations when compared with Anderson's values (there is a graph on page 444 of his book). The vertical component of velocity, v, looks similar after 4 iterations though. This makes me think the horizontal component is incorrect. However, I can't find anything wrong with it in my code!

Having said all that, none of this makes a huge difference to the flowfield as the horizontal component of velocity, u, still looks very linear when it has converged but nevertheless, I would like to know what (small) thing I am doing wrong.

Can you please let me know the error of my ways? My code is attached for your perusal. Thanks in advance!
Attached Files
File Type: docx Couette Flow in Matlab.docx (14.5 KB, 395 views)
zapico and Gal!Leo like this.
mulchah is offline   Reply With Quote

Old   February 19, 2014, 06:27
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,842
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I don't think you can get the exactly same numbers... precision, round-off, etc. are some elements that can affect the number. However, I see that your magnitude order is the same...
FMDenaro is offline   Reply With Quote

Old   February 19, 2014, 06:45
Default
  #3
New Member
 
Jack
Join Date: Dec 2013
Posts: 6
Rep Power: 12
mulchah is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I don't think you can get the exactly same numbers... precision, round-off, etc. are some elements that can affect the number. However, I see that your magnitude order is the same...
Really? But my u velocity field doesn't match Anderson's after K = 4 iterations either. For example, at a vertical height of 0.09 ft, my u velocity prediction is approximately 0.4 ft/s while Anderson's is closer to 0.3 ft/s. This makes me think I am doing something wrong.

If what you are saying is true, then surely my velocity prediction should be a lot closer to his values.
mulchah is offline   Reply With Quote

Old   February 19, 2014, 06:55
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,842
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
check the difference in the velocity profile at the steady state
FMDenaro is offline   Reply With Quote

Old   February 25, 2014, 05:41
Default
  #5
New Member
 
Jack
Join Date: Dec 2013
Posts: 6
Rep Power: 12
mulchah is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
check the difference in the velocity profile at the steady state
Velocity profile in steady state is almost perfectly linear. As I have said already, the code converges in steady state to give the final velocity profile but results from the interim steps do not match Anderson's, as given in his textbook. Thus I know I am doing something wrong but the frustrating thing is, I don't know what it could be.
mulchah is offline   Reply With Quote

Reply

Tags
cfd code, couette flow, john anderson, matlab, pressure correction

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Pressure Correction Method for Cylindrical Coordinates Eslam Reda Main CFD Forum 4 July 4, 2009 02:15
flow simulation across a small fan jane luo Main CFD Forum 15 April 12, 2004 17:49
pressure driven flow by pressure correction method justentered Main CFD Forum 0 December 29, 2003 23:52
pressure gradient term in low speed flow Atit Koonsrisuk Main CFD Forum 2 January 10, 2002 10:52
Method for Stokes flow with continuous pressure gradients? Matthew Cordery Main CFD Forum 0 August 31, 1998 17:58


All times are GMT -4. The time now is 08:05.