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

Compare 2D axi-symmetric and 3D. Is OpenFOAM calculate the velocity wrong?

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

Like Tree4Likes
  • 4 Post By novyno

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 11, 2010, 06:02
Default Compare 2D axi-symmetric and 3D. Is OpenFOAM calculate the velocity wrong?
  #1
Member
 
Join Date: Jun 2010
Posts: 33
Rep Power: 16
RalphS is on a distinguished road
Hi,

I simulate a swirl flow. A big cylinder (in which I took a displacer) following a small cylinder.
I try 2 dimensional axi-symmetric and 3 dimensional. All settings are equal. I used pisoFoam, Rankine vortex at inlet

I think OpenFOAM have problems with calculating the velocity in a 2D case. Enclose two picture.

U_x at the inlet is 10 m/s.
The continuity equation says: u1*A1=u2*A2

In my case, the velocity in the small passage in the middle should be 16 m/s.

u2 = u1 * (r1²/(r1²-rd²)) (With rd, I mean the radius of the displacer)
3D its right!

But in the 2D case the velocity is to high. Does OpenFOAM calculate the area ratio wrong?
If i would calculate like u2 = u1 *(r1/(r1-rd)). I get the velocity in the 2D axi-symmetric case.

The second problem is, that the flow at the axis makes no sense. Here I think OpenFOAM have problems with radial exchange in a 2D axi-symmetric case.

Have you any idea why the velocity is such different between these two cases?
Could it be, that I have forgot something to set in the 2d axi-symmetric case, or is there an additional setting to control something like area ration etc.

Regards
Ralph
Attached Images
File Type: jpg 2D.jpg (12.6 KB, 322 views)
File Type: jpg 3D.jpg (15.0 KB, 289 views)
RalphS is offline   Reply With Quote

Old   June 14, 2010, 12:09
Default 2d case/wedge, is radial area ratio respected?
  #2
New Member
 
Johannes Kneer
Join Date: Mar 2009
Location: Germany, Karlsruhe
Posts: 13
Rep Power: 17
johannesk is on a distinguished road
Hi Ralph,

nice to meet you here

If I understand you correctly, your problem with the 2d wedge case is as follows: the inlet is at a smaller radius than the outlet and has a smaller area as the outlet. Yet the results for U from OF look like they were the same area?

Quote:
Originally Posted by RalphS View Post
Have you any idea why the velocity is such different between these two cases?
Could it be, that I have forgot something to set in the 2d axi-symmetric case, or is there an additional setting to control something like area ration etc.
Could you post your polyMesh/boundary file of the 2d case, or maybe even the complete 2d case?

cheers,
Johannes
johannesk is offline   Reply With Quote

Old   June 14, 2010, 12:39
Default
  #3
Senior Member
 
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 23
l_r_mcglashan will become famous soon enough
Based on what the continuity equation gives you, it sounds to me as if you're not simulating a wedge, but an infinite plane. Are you following this:

http://openfoamwiki.net/index.php/Ma...s/AxiSymmetric
__________________
Laurence R. McGlashan :: Website
l_r_mcglashan is offline   Reply With Quote

Old   June 14, 2010, 13:08
Default
  #4
New Member
 
ernest
Join Date: Jun 2010
Posts: 21
Rep Power: 15
erncyc is an unknown quantity at this point
Hi Johannes...I am a new user to OpenFOAM and I am trying to simulate flow in a cyclone. I imported the mesh from a program caedium which created the .ofm. I get the following errors when i try to run please help thank you. (Ia m using simplefoam). I have also enclosed the offending file (ie RASproperties). thanks!

Errors below:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : c:\cfd\OpenFOAM-1.5\bin\simpleFoam
Date : Jun 14 2010
Time : 20:07:31
Host : hostname
PID : 3416
Case : C:\Users\user\Documents\PhD\CFD\Openfoam\My-tutorials\cyclone\cyclone-flow
nProcs : 1

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
Programme terminated with errors: exit code 1, status 0.
Error messages:


keyword kEpsilonCoeffs is undefined in dictionary "C:\Users\user\Documents\PhD\CFD\Openfoam\My-tutorials\cyclone\cyclone-flow\constant\RASProperties"

file: C:\Users\user\Documents\PhD\CFD\Openfoam\My-tutorials\cyclone\cyclone-flow\constant\RASProperties from line 14 to line 15.

From function dictionary::subDict(const word& keyword)
in file C:\tmp\OpenFOAM-1.5\src\OpenFOAM\db\dictionary\dictionary.C at line 288.

FOAM exiting



Offending file below:

/*--------------------------------*- C++ -*----------------------------------*/
// OpenFOAM 1.6.x
// Generated by Symscape's Caedium on: 14 Jun 2010 16:20:41, www.symscape.com

FoamFile
{
version 2.0;
format binary;
class dictionary;
object RASProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

RASModel kEpsilon;
turbulence on;




erncyc is offline   Reply With Quote

Old   November 12, 2010, 11:17
Default Split wedge patch...
  #5
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 17
novyno is on a distinguished road
Hi RalphS,

You probably have defined a single wedge patch (called for example front&Back) in wich you have put every wedge surfaces.

Try to generate two different patchs for every "wedge side" (for example called front and back). This two patch of course will be "specular".

Now when you run your simulations the continuity equation will give you the "naturals" expected results. :-)

Regards!
__________________
...The best way to acquire new knowledge is to share it...
novyno is offline   Reply With Quote

Old   November 13, 2010, 07:00
Default
  #6
Member
 
Join Date: Jun 2010
Posts: 33
Rep Power: 16
RalphS is on a distinguished road
Close the thread please.

Why you answered a question from june?
RalphS is offline   Reply With Quote

Old   November 13, 2010, 20:51
Default
  #7
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 17
novyno is on a distinguished road
I have answered becouse the thread original question wasn't solved (or the solution was not clearly described). I had the same problem some times ago and I hoped that my answer could help someone who encounter the same problems with wedge patchs.
krene, hfsf, Rokerda and 1 others like this.
__________________
...The best way to acquire new knowledge is to share it...
novyno is offline   Reply With Quote

Reply

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



All times are GMT -4. The time now is 07:27.