|
[Sponsors] |
Simulating Creeping (Stokes) Flow in OpenFOAM |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 27, 2014, 12:36 |
Simulating Creeping (Stokes) Flow in OpenFOAM
|
#1 |
New Member
Kevin Shanley
Join Date: Jun 2014
Location: New Paltz NY
Posts: 2
Rep Power: 0 |
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 |
|
August 16, 2014, 10:24 |
|
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,982
Blog Entries: 45
Rep Power: 128 |
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" ):
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:
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
__________________
|
|
September 5, 2014, 21:32 |
|
#3 |
New Member
Kevin Shanley
Join Date: Jun 2014
Location: New Paltz NY
Posts: 2
Rep Power: 0 |
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! |
|
July 9, 2015, 16:41 |
|
#4 |
New Member
jafar
Join Date: Jul 2014
Posts: 16
Rep Power: 12 |
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 |
|
Tags |
creeping flow, cubic law, openfoam 2.2, simplefoam, stokes flow |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |