CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Help needed with laminar couette flow boundary conditions

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 19, 2009, 09:25
Default Help needed with laminar couette flow boundary conditions
  #1
New Member
 
Join Date: Mar 2009
Location: Germany
Posts: 4
Rep Power: 17
bengt is on a distinguished road
Dear Foamers,

first of all, THANKS A LOT for the great support in the forum and the great job the OpenFoam development team is doing!! Still being a beginner with cfd and OpenFoam, I am currently trying to simulate a microscopic incompressible couette flow for a measured sliding contact geometry, using simpleFoam. The aim is a comparison of the calculated friction torque with the experimental values. The computational domain (170 x 170 µm) is a portion of a fully flooded tribological sliding contact, as similarly described by P. Brajdic in the OpenFoam thesis repository. The fluid film thickness between the fixed upper wall and the moving lower wall (sliding in x-direction) is in the range of approx. 0,1..10 µm; the blockMesh consists of approx. 2,8Mcells of hex type and checkMesh tells me that the mesh quality is fine. The main problem now is the correct definition of the boundary conditions, since it's very difficult to find information on this rather special case. After some trials with exploding models , I am currently using the following, but I am not sure if they are correct:

"inlet" (fluid film cross section normal to x): p = 0, u = zeroGradient
"outlet"(see above): p = 0, u = zeroGradient
side1 (fluid film cross section parallel to x): p = 0, u = zeroGradient
side2 (as above): p = 0, u = zeroGradient
fixed (upper) wall: p = zeroGradient, u = (0, 0, 0)
moving (lower) wall: p = zeroGradient, u = (Ux, 0, 0)

The patch types in the boundary file are all "patch" except for the walls being set to type "wall".
Generally, after approximately 4,5 thousand steps the results for x are looking fine (at least for me ), but there seems to be some convergence issue concerning Uy, Uz and p. However, the fluid film thickness (or gap height) does not vary at all in the sliding direction x (one fluid film cross section simply extruded in x), so I think that the algorithm only deals with the numerical errors regarding Uy, Uz and p . The logfile entries are as follows:

Time = 4795

smoothSolver: Solving for Ux, Initial residual = 8.28444717201096e-07, Final residual = 1.17861528964069e-16, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.000420799816002115, Final residual = 1.15677735105817e-14, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.00045689356317463, Final residual = 1.21403757603192e-14, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00246394853928313, Final residual = 1.12631494042653e-05, No Iterations 3
time step continuity errors : sum local = 6.29921027180131e-09, global = -1.88504368307571e-10, cumulative = -8.84159167321985e-06
ExecutionTime = 168453.17 s ClockTime = 168462 s

Time = 4796

smoothSolver: Solving for Ux, Initial residual = 8.27636882199578e-07, Final residual = 1.17828578411332e-16, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.000420804884209455, Final residual = 1.15685143238247e-14, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.000456888391577294, Final residual = 1.21400758956957e-14, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0024639398031611, Final residual = 1.12633772620776e-05, No Iterations 3
time step continuity errors : sum local = 6.29377884031872e-09, global = -1.88393757795753e-10, cumulative = -8.84178006697764e-06
ExecutionTime = 168488.27 s ClockTime = 168497 s

Time = 4797

smoothSolver: Solving for Ux, Initial residual = 8.26829853129783e-07, Final residual = 1.17743398700947e-16, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.000420809949537425, Final residual = 1.15680804424955e-14, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.000456883224658707, Final residual = 1.21399060586643e-14, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00246393105191049, Final residual = 1.12636035684504e-05, No Iterations 3
time step continuity errors : sum local = 6.2883465087325e-09, global = -1.88280746920101e-10, cumulative = -8.84196834772456e-06
ExecutionTime = 168523.53 s ClockTime = 168532 s

The solver tolerances for both U and p are 1e-7 and 0.01. Ux is in the range of 0..20e-3 m/s.

Could anybody please give me some feedback if my assumptions are correct? I'd appreciate your help very much!!

Have a nice day, kind regards
Bengt
bengt is offline   Reply With Quote

Reply

Tags
boundary conditions, convergence, couette, simplefoam

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
boundary layer specifications for laminar flow ranjith CFX 7 October 16, 2017 09:09
Help needed with boundary conditions JP FLUENT 1 October 2, 2008 03:30
boundary conditions - help needed!!! Elena Saldaeva CFX 1 May 23, 2008 04:05
Boundary Conditions and solver for low Mach laminar air jet fbisetti OpenFOAM Running, Solving & CFD 1 September 13, 2006 13:49
Help needed FLUENT- Flow over rough boundary Trevor F FLUENT 0 October 17, 2002 21:29


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