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

non-null velocity at walls despite noSlip condition

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

Like Tree7Likes
  • 1 Post By New_Old
  • 2 Post By Kraneberger
  • 2 Post By snak
  • 2 Post By snak

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 31, 2016, 05:41
Question non-null velocity at walls despite noSlip condition
  #1
New Member
 
Jay B.
Join Date: Nov 2016
Posts: 12
Rep Power: 6
New_Old is on a distinguished road
Hi all!

It's all in the title.

- OF version: v1606+
- solver: simpleFoam
- general problem description: "quick and dirty" laminar flow simulation with stream wise flow periodicity and symmetry BCs on the model boundaries which are not cyclic or walls.

The mesh and other numerical settings might not be optimal, but the simulation runs smoothly and with 20 layers of elements across the domain (10 in the model, x2 due to symmetry BC), I would expect the zero velocity condition at the walls to be enforced without a problem.

I attach here the complete setup, minus the Gmsh-generated mesh file (too big). Also, a screen capture of the velocity at the walls, which will give you a good idea of what the missing mesh looks like.

Basically, the velocity field at the wall in the xy-plane is almost identical to the one on the symmetry plane 10 prism layers away...

Any clue?
Attached Files
File Type: zip periodicity_streamwise_simpleFoam.zip (77.1 KB, 2 views)
kcavatar likes this.
New_Old is offline   Reply With Quote

Old   December 31, 2016, 07:20
Lightbulb Update
  #2
New Member
 
Jay B.
Join Date: Nov 2016
Posts: 12
Rep Power: 6
New_Old is on a distinguished road
I've switched the initial condition at the walls from noSlip to fixedValue, and this enforced zero velocity successfully.

Apparently the noSlip condition is new, and is supposed to be equivalent to fixedValue with a value of (0 0 0):
https://github.com/OpenFOAM/OpenFOAM...98d0b13f22b515

Maybe a bug?
New_Old is offline   Reply With Quote

Old   September 13, 2017, 02:37
Default
  #3
New Member
 
Ben
Join Date: Oct 2016
Posts: 25
Rep Power: 7
bentkj is on a distinguished road
Hi, although late, are you able to confirm this? I've been using the noSlip condition for U at my wing wall I'm still getting nonzero values for the surface velocities.

edit: i'm guessing there's something to do with how the data is viewed from paraview as well? I'm unsure about this. Hoping someone could enlighten!
bentkj is offline   Reply With Quote

Old   December 15, 2017, 09:48
Default
  #4
Member
 
Lennart
Join Date: Feb 2016
Posts: 46
Rep Power: 7
elmo555 is on a distinguished road
I'm also experiencing non-zero velocities at my noSlip boundaries (when viewed in paraview). Might this be due to the mapping from cells to faces? The velocity in the cell center of the boundary cells is not equal to zero, I guess, and mapping or interpolating this value to the face may lead to non-zero wall velocities. Is there a way to access the face values?
elmo555 is offline   Reply With Quote

Old   February 10, 2021, 23:01
Default
  #5
New Member
 
Sumit
Join Date: Jul 2017
Posts: 4
Rep Power: 6
sumitzanje is on a distinguished road
Has anyone found any detailed information on this? I'm also getting non zero velocity at walls with noSlip BC's.
sumitzanje is offline   Reply With Quote

Old   August 20, 2021, 03:52
Default noSlip not zero on wall
  #6
New Member
 
Join Date: Oct 2019
Posts: 10
Rep Power: 3
Kraneberger is on a distinguished road
I can confirm the problem with the noSlip boundary condition on walls. I ran the tutorial airfoil2D and attached the velocity profiles at the wall. y-axis is U_magnitude and x-axis runs from underneath the profile to above. I tested the boundary conditions:


1. type fixedValue;
value uniform (0 0 0);


2. type movingWallVelocity;
value uniform (0 0 0);


3. type noSlip;


the first two are identical, as expected and the third does not reach zero at the wall.


Greetings,
Kraneberger


fixedValue.png

movingWallVel.png

noSlipVel.png
granzer and piu58 like this.
Kraneberger is offline   Reply With Quote

Old   August 21, 2021, 07:40
Default
  #7
Member
 
shinji nakagawa
Join Date: Mar 2009
Location: Japan
Posts: 97
Blog Entries: 1
Rep Power: 15
snak is on a distinguished road
Hi,


How did you take data and plot ?


With tutorials/simpleFoam/airFoil2D case, velocities at the wall is zero with noSlip condition as shown here,

Screenshot from 2021-08-21 20-35-29.png


The picture above is obtained only with patch/walls (without internalMesh) using ParaView.
snak is offline   Reply With Quote

Old   August 23, 2021, 11:50
Default
  #8
New Member
 
Join Date: Oct 2019
Posts: 10
Rep Power: 3
Kraneberger is on a distinguished road
Hi snak,


I obtained the plots via "Plot over Line". What version of OpenFoam did you run the tutorial with. I used v1912.



In the attached picture you can see that I get the same results when only importing walls in paraview.


walls noSlip.jpg


Greetings,
Kraneberger
Kraneberger is offline   Reply With Quote

Old   August 23, 2021, 20:42
Default
  #9
Member
 
shinji nakagawa
Join Date: Mar 2009
Location: Japan
Posts: 97
Blog Entries: 1
Rep Power: 15
snak is on a distinguished road
Hi Kraneberger,

I used OpenFOAM v2106. With v1912, I got the same result as shown here:
Screenshot from 2021-08-24 09-24-21.png

I used paraFoam command. paraFoam will stat paraview with the OpenFOAM libraries and reader modules.

When I use a native openfoam reader (*.foam), wall velocity is NOT zero as you say.

It will be better to use openfoam utility such as sample to extract data. PlotOverLine interpolates data and not exactly show your result in some case.
granzer and Kraneberger like this.
snak is offline   Reply With Quote

Old   August 24, 2021, 03:47
Default
  #10
New Member
 
Join Date: Oct 2019
Posts: 10
Rep Power: 3
Kraneberger is on a distinguished road
Thankyou for taking the time, snak.


I tried it the way you suggested and I can find the zero velocities now. I was confused, because I did the same thing three times in a row. For movingWallVelocity and fixedValue I got identical results and with noSlip the graph changed. The mistake must be on my end and is a post-processing error, not an OpenFoam calculation error ;-)


Greetings,
Kraneberger
Kraneberger is offline   Reply With Quote

Old   August 24, 2021, 04:21
Default
  #11
Member
 
shinji nakagawa
Join Date: Mar 2009
Location: Japan
Posts: 97
Blog Entries: 1
Rep Power: 15
snak is on a distinguished road
Hi Kraneberger,

I found the open issue at VTK site.
https://gitlab.kitware.com/vtk/vtk/-/issues/18085

Using ParaView without a OpenFOAM Plugin (or not using paraFoam) causes this misleading visualization. Using fixedValue instead of noSlip may be convenient for you if you want visualize your results without reconstruction etc.


Added:
The problem with noSlip b.c. wll be fixed after the release of ParaView 5.9.1. There is no problem with the recent nightly built like ParaView-5.9.1-1583-g94277a4...

https://www.paraview.org/download/?version=nightly
granzer and Kraneberger like this.

Last edited by snak; August 24, 2021 at 05:33.
snak is offline   Reply With Quote

Old   August 24, 2021, 06:40
Default
  #12
New Member
 
Join Date: Oct 2019
Posts: 10
Rep Power: 3
Kraneberger is on a distinguished road
I see. Thanks for finding the bug report. I will use a more up to date paraview version from now on :-)
Kraneberger is offline   Reply With Quote

Reply

Tags
cyclic, noslip, 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
slip & noslip boundary condition for dynamic mesh, how to implement? oschi OpenFOAM Programming & Development 6 May 21, 2019 04:33
Massflow or average velocity boundary condition Sideshore OpenFOAM Pre-Processing 6 September 27, 2016 15:09
Boundary condition of velocity and pressure at interface for air water pipe flow jignesh_thaker2007 OpenFOAM Running, Solving & CFD 0 June 10, 2014 16:42
Condition of Air velocity in Open Channel simulation wes1204 OpenFOAM Running, Solving & CFD 0 May 9, 2014 05:08
Rotating wall velocity boundary condition for axisymmetric flow problem diwakaranant SU2 1 March 23, 2013 23:45


All times are GMT -4. The time now is 19:34.