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

Simulating Creeping (Stokes) Flow in OpenFOAM

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

Like Tree3Likes
  • 2 Post By wyldckat
  • 1 Post By shanleyk

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 27, 2014, 12:36
Default Simulating Creeping (Stokes) Flow in OpenFOAM
  #1
New Member
 
Kevin Shanley
Join Date: Jun 2014
Location: New Paltz NY
Posts: 2
Rep Power: 0
shanleyk is on a distinguished road
Hello -

I have considerable experience with CFD, but I’m brand new to OpenFOAM. I was able to work my way through the lid driven cavity flow and dam breaking tutorials without really any issues. However, now I am branching off on my own and having some difficulty.


I am trying to start out easy. I am simulating 2D flow between parallel plates. The flow is steady and Re is very low ~0.001 (creeping/Stokes flow). After some research I decided to use simpleFoam. It runs well and convergence looks good. Qualitatively the results look good when I view in paraFoam. However, the pressure drop does not obey the cubic law. The pressure drop across the domain should vary with the cube of the height. It does not. In fact it seems to vary linearly.


I have some more complicated fracture flows to look at, so I want to make sure I am doing this correctly.


Is simpleFoam not the best solver for this application?
Have I set it up incorrectly?
I am using OF 2.2.



I have uploaded a tarball case setup. Included in there is also the log from a run so you can see residuals. You will have to run blockMesh to generate the mesh.



Thank you for any assistance.


Kevin
Attached Files
File Type: zip channelTest.zip (11.8 KB, 63 views)
shanleyk is offline   Reply With Quote

Old   August 16, 2014, 10:24
Default
  #2
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,982
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Greetings Kevin and welcome to the forum!

I've kept a reminder to check this thread when I had the time and today I finally managed to look into this.

So, the issues I can figure out are as follows (I'm referring to a pipe, since it's somewhat similar and easier to write than "parallel plates" ):
  • The pressure drop along the pipe's length is meant to be constant, so there was no problem there.
  • The "fvSolution" was configured to have a very tight "tolerance" values in the "solvers" block. I lighten from 1e-17 to 1e-6 and it would run a lot faster
  • The "divSchemes" in "fvSchemes" are mostly using "upwind". Switching to "linear" should improve the resulting solution.
  • The mesh resolution along the pipe length and diameter is out of proportion. This wouldn't be much of a problem if you were only trying to assess how the flow goes along the pipe. But from what I can see, perhaps you also need more resolution along the diameter, since you're looking for a cubic relation between distance to the walls and pressure variation.
I changed the mesh resolution from "10 1 10000" to "100 1 1000". The results seem mostly the same, but at least the flow speed profile looks more refined.


Problem is that I can't figure out where exactly you're looking for a cubic profile for the pressure. The results I get with these improvements indicates that:
  • the pressure is constant at each cross-section of the pipe, i.e. for a fixed X position I get a constant pressure along the diameter;
  • the pressure drop is constant, which is in accordance with the analytical solution;
  • the flow speed (velocity) profile in the same cross-section does have a parabolic shape which doesn't seem to be quadratic.


Can you provide more details on the simulation you're trying to do and the analytical reference data you're looking at? Because perhaps there was something lost in the translation to the OpenFOAM case, such as whether the two planes should be at the same pressure or perhaps some sort of friction should be present at the walls?

Best regards,
Bruno
rudolf.hellmuth and bagherij like this.
__________________
wyldckat is offline   Reply With Quote

Old   September 5, 2014, 21:32
Default
  #3
New Member
 
Kevin Shanley
Join Date: Jun 2014
Location: New Paltz NY
Posts: 2
Rep Power: 0
shanleyk is on a distinguished road
Wyldckat,
Thank you very much. I had forgotten that I had posted this.
We got our case running. I don't remember all of the details, but mesh was an issue and a few other things.
The cubic law relationship I was looking for is that (aperture height)^3 should be related to pressure drop along the channel. We were able to show this with simpleFoam. We've moved on to more complex geometries and getting nice results.
Thanks!
wyldckat likes this.
shanleyk is offline   Reply With Quote

Old   July 9, 2015, 16:41
Default
  #4
New Member
 
jafar
Join Date: Jul 2014
Posts: 16
Rep Power: 12
bagherij is on a distinguished road
Dear
i use u file for creeping flow around cylnder between to plate.
i want to calculate the magnitude of shear rate tensor(|tau|) which
defined by
tau = -nu()*(fvc::grad(U)+(fvc::grad(U)).T());
|tau|=sqrt(1/2tr(tau^2)

where "nu()" is the dynamic viscosity of the system

where should i add to code ,and how can i get |tau| and compare it with some parameter.
Thanks
bagherij is offline   Reply With Quote

Reply

Tags
creeping flow, cubic law, openfoam 2.2, simplefoam, stokes flow

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
About Some Concepts:Laminar flow, turbulent flow, steady flow and time-dependent flow Jing Main CFD Forum 8 October 5, 2018 18:02
Materials / Tutorials for simulating flow over spillway basud0817 FLOW-3D 4 August 25, 2015 05:39
pipe flow OpenFOAM stix OpenFOAM Running, Solving & CFD 1 March 4, 2013 17:44
Bidomain Stokes flow bluemoth OpenFOAM Running, Solving & CFD 0 December 29, 2005 00:36
steady state creeping flow dominik Main CFD Forum 4 March 29, 2004 10:29


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