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

Modelling Thin film lubrication in a thrust bearing

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 13, 2015, 22:44
Default Modelling Thin film lubrication in a thrust bearing
  #1
New Member
 
arvind prabhakar
Join Date: Nov 2014
Posts: 2
Rep Power: 0
arvind9008 is on a distinguished road
Hi Everyone,

I am new to this forum. I am trying to model a 2-dimensional incompressible flow in a thin air lubricating film i.e. a in 1d bearing (Re between 100 to 1000) but, I am not getting the resulting correct pressure profile.
The model geometry and flow conditions are essentially a couette flow problem (Only top wall moves causing velocity profile, I have attached image of flow domain for your reference) however, the bottom wall is stationary and non parallel to the top wall ( This I believe requires grid transformation OR non-dimensionalizing the equations as per geometry, in my case I have non dimensionalized as per film thickness to length ratio (epsilon) ).
I have written the code for couette flow for solving 2 dimensional Navier stokes equations (using simple algorithm on staggered grid which was also available on this forum) and it works correctly when epsilon is set equal to 1 for a standard couette flow geometry as per results from John Anderson's book "CFD - Basics with Applications", section 9.4. (I have attached the matlab code as well as a pdf file) However, The section does not discuss anything about the pressure field output from couette flow so I have no idea about what kind of pressure field should be generated.
I would really appreciate help from you guys on this forum to guide me in solving this problem on the geometry I have attached.

Thanks and Regards,
Arvind
Attached Images
File Type: png 1d bearing geometry.png (23.5 KB, 42 views)
Attached Files
File Type: pdf couette flow.pdf (140.5 KB, 38 views)
arvind9008 is offline   Reply With Quote

Old   December 17, 2015, 06:47
Default
  #2
Senior Member
 
Troy Snyder
Join Date: Jul 2009
Location: Akron, OH
Posts: 219
Rep Power: 18
tas38 is on a distinguished road
Quote:
Originally Posted by arvind9008 View Post
Hi Everyone,

I am new to this forum. I am trying to model a 2-dimensional incompressible flow in a thin air lubricating film i.e. a in 1d bearing (Re between 100 to 1000) but, I am not getting the resulting correct pressure profile.
The model geometry and flow conditions are essentially a couette flow problem (Only top wall moves causing velocity profile, I have attached image of flow domain for your reference) however, the bottom wall is stationary and non parallel to the top wall ( This I believe requires grid transformation OR non-dimensionalizing the equations as per geometry, in my case I have non dimensionalized as per film thickness to length ratio (epsilon) ).
I have written the code for couette flow for solving 2 dimensional Navier stokes equations (using simple algorithm on staggered grid which was also available on this forum) and it works correctly when epsilon is set equal to 1 for a standard couette flow geometry as per results from John Anderson's book "CFD - Basics with Applications", section 9.4. (I have attached the matlab code as well as a pdf file) However, The section does not discuss anything about the pressure field output from couette flow so I have no idea about what kind of pressure field should be generated.
I would really appreciate help from you guys on this forum to guide me in solving this problem on the geometry I have attached.

Thanks and Regards,
Arvind
More commonly, the slider bearing is solved via the thin film equation, i.e. Reynolds eqn. (Couette and Poiseuille velocity profiles dumped into continuity eqn.). See the solution for the plane slider in the following link (includes plots of pressure profile):

http://rotorlab.tamu.edu/me626/Notes...D_bearings.pdf

I have not looked at your code in detail but I suspect your issue may lie in the boundary conditions. What are the boundary conditions for the momentum and pressure correction equations at the inflow and outflow?
tas38 is offline   Reply With Quote

Old   December 18, 2015, 20:47
Default
  #3
New Member
 
arvind prabhakar
Join Date: Nov 2014
Posts: 2
Rep Power: 0
arvind9008 is on a distinguished road
Hi,
Thanks for the reply.
I have the codes for reynold's equation applied to thin films. It predicts pressure profile well however, due to neglection of inertia terms the pressure predicted is not as accurate as the one predicted by solution of complete navier stokes.
the case i am analyzing is for Re=1000 and higher. inertia effects become predominant therefore, requiring modelling of complete navier stokes.
although, I have to model the compressible form of navier stokes equations but i want to ensure I have results for the incompressible modelling first so I can be sure and move to modelling the compressible system of equations.
This is kind of stressing me out because, This is to form part of my thesis that I have to defend in spring 16.
Any inputs from you guys will be of help to me.
Thanks,
Arvind
arvind9008 is offline   Reply With Quote

Old   December 19, 2015, 04:29
Default
  #4
New Member
 
Syed Aaquib Hazari
Join Date: Apr 2015
Location: Delft
Posts: 16
Rep Power: 11
hazari.aaquib is on a distinguished road
Hi,
It maybe that turbulence may affect your case.
Are you sure there's no turbulence affecting your case?
Since Re is around 1000 there could possibly be traces of turbulence.
hazari.aaquib is offline   Reply With Quote

Old   December 19, 2015, 05:17
Default
  #5
Senior Member
 
Troy Snyder
Join Date: Jul 2009
Location: Akron, OH
Posts: 219
Rep Power: 18
tas38 is on a distinguished road
Arvind,

Getting back to your original post, I was under the impression that the only issue was that you were unsure if your pressure profile was 'correct'. Please clarify your issues as the Reynolds codes should give you an idea of what your pressure profile should look like.

If your concern is inertia or turbulence muddying the waters, simply run a case with smaller film height, lower sliding velocity, or higher viscosity (take your pick) to drop the Re # and then compare with the Reynolds code results.
tas38 is offline   Reply With Quote

Reply

Tags
2 dimensional, cfd code, couette flow, incompressible, navier stokes equations

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
Evaporation of thin water film Stephen Waite CFX 6 July 3, 2014 22:14
Multiphase Journal Bearing CFD: Modelling Air Entrainment eaxym CFX 1 May 13, 2014 19:29
VOF - Eliminating shear on the surface of a thin liquid film saisanthoshm88 Fluent Multiphase 0 October 27, 2013 12:22
Air bearing modelling Dave Larkham CFX 1 October 19, 2008 19:10
Wall Film modelling Computational Dynamics Siemens 0 May 11, 2001 06:52


All times are GMT -4. The time now is 10:40.