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

MRFSimpleFoam: wrong boundary conditions on rotating walls

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

Like Tree1Likes
  • 1 Post By cves

Reply
 
LinkBack Thread Tools Display Modes
Old   April 15, 2010, 11:37
Default MRFSimpleFoam: wrong boundary conditions on rotating walls
  #1
New Member
 
Christian
Join Date: Feb 2010
Location: Payerne, Switzerland
Posts: 9
Rep Power: 7
cves is on a distinguished road
Hello,


I am using OpenFoam1.5-dev, svn version 1677, with SIG Turbomachinery add-ons in order to perform the ERCOFTAC centrifugal pump case-study. I have found that the two faceSet made for the setup of MRFSimpleFoam lead to a wrong velocity value on the rotating walls. Indeed the velocity at the wall is set to zero instead of Omega x Radius.


I have followed the instructions and setups given in the openFoam Wiki :


http://openfoamwiki.net/index.php/Si...vaned_diffuser


However, I had to decrease the k and epsilon relaxation factor value from 0.7 in the original system/fvSolution file to 0.5 in order to achieved 5000 time steps. I give you the images (U.png and linear.png) of the velocity and the convergence that I get after the 5000 time steps.



I have compared these figures with the ones given by :


http://www.tfd.chalmers.se/~hani/pdf_files/OFW4_2009_Petit_slides.pdf


which has the same setups except that the svn version is 1240. I saw that the problem comes from the wrong velocity imposed on the rotating blades as you can see the small blue lines around the rotating blades which means a zero velocity value is imposed (U_zoom.png).


I found that the mixerVessel2D has not this problem and the only difference between these two tutorials was the two faceSet for system/faceSetDict_rotorFaces and system/faceSetDict_noBoundaryFaces, which are not done in the mixerVessel2D tutorial. So I remove this two faceSet for the ERCOFTAC centrifugal pump tutorial and run 5000 times steps (see the convergence and velocity profile after 5000 time steps U_noFaceSet.png and linear_noFaceSet.png).


In this case the velocity is correct on the rotating blade but the computation converge into a non-physical result.


If someone has the solution to correct the boundary value on the rotating blades and get a well physical converged solution, please let me know.


Best regards


Christian
Attached Images
File Type: jpg U.jpg (40.5 KB, 65 views)
File Type: png linear.png (9.3 KB, 55 views)
File Type: jpg U_zoom.jpg (30.6 KB, 60 views)
File Type: jpg U_noFaceSet.jpg (32.6 KB, 56 views)
File Type: png linear_nofaceset.png (6.3 KB, 50 views)
cves is offline   Reply With Quote

Old   April 20, 2010, 03:47
Default Correction of the ercoftacCentrifugalPump tutorial
  #2
New Member
 
Christian
Join Date: Feb 2010
Location: Payerne, Switzerland
Posts: 9
Rep Power: 7
cves is on a distinguished road
Hello,

After a few days of research, I have found a way which corrects the ercoftacCentrifugalPump tutorials and can be used as first setup for a MRFSimpleFoam case study for OpenFOAM-1.5-dev svn version greater than 1242. Indeed the faceSet are useless in these versions, OpenFOAM is able to found the faces for MRFSimpleFoam by itself.

This correction leads to a good convergence without any tuning of the relaxation factors. I have attached some pictures of my results. Moreover the converged result is physical and close to the result of:

http://www.tfd.chalmers.se/~hani/pdf_files/OFW4_2009_Petit_slides.pdf

You just have to replace the makeMesh and constant/MRFZones files of your tutorials by the ones that I have attached to this post. My files are for the stiched mesh but you can also use it for the ggi mesh by adding the ggi boundaries in the nonRotatingPatches of the MRFZones file.

Best regards

Christian
Attached Images
File Type: png linear.png (5.4 KB, 47 views)
File Type: jpg U.jpg (36.9 KB, 55 views)
File Type: jpg p.jpg (38.9 KB, 48 views)
Attached Files
File Type: txt makeMesh.txt (1.5 KB, 94 views)
File Type: txt MRFZones.txt (1.1 KB, 105 views)
cves is offline   Reply With Quote

Old   April 20, 2010, 07:13
Default
  #3
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 307
Rep Power: 9
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618
Also in 1.6
waynezw0618 is offline   Reply With Quote

Old   May 3, 2010, 03:38
Default GGI case
  #4
New Member
 
Christian
Join Date: Feb 2010
Location: Payerne, Switzerland
Posts: 9
Rep Power: 7
cves is on a distinguished road
Hello,

Here are the modified files to run the GGI case correctly.

Christian
Attached Files
File Type: txt makeMesh.txt (1,015 Bytes, 122 views)
File Type: txt MRFZones.txt (1.1 KB, 139 views)
marcelgt87 likes this.
cves is offline   Reply With Quote

Reply

Tags
boundary condition, centrifugal pump, mrfsimplefoam, velocity

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Wind turbine simulation Saturn CFX 34 October 16, 2014 05:27
inlet velocity boundary condition murali CFX 5 August 3, 2012 08:56
Boundary conditions for rotating fan rohit CFX 2 July 22, 2008 09:20
Rotating interpolated boundary condition hani OpenFOAM Running, Solving & CFD 0 July 4, 2006 07:09
Boundary Conditions Jan Ramboer Main CFD Forum 11 August 16, 1999 08:59


All times are GMT -4. The time now is 01:24.