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

MRF solvers give diffenrent results when MRF zone's diameter is different

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 14, 2013, 04:25
Default MRF solvers give diffenrent results when MRF zone's diameter is different
  #1
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 849
Rep Power: 18
sharonyue is on a distinguished road
Hi,

Nowdays I intend to make a test of my mrf case's result with the paper.But I find the result vary with the different diameter of my rotating size.I made a comparision in two situation in which only the MRF zones size is different.As you can see in the image.(its a slice of a 3D case)

The first one's MRF zone is smaller than the other.expect this everything is the same.But after 2000 iterations,the velocity is a little different.I dont know if it gets convergence but in 2500s iterations,the velocity looks dont change anymore.

in this thread linnemann happend to run into this situation. but I didnot try his case.http://www.openfoam.org/mantisbt/view.php?id=409

Code:
Time = 2093

smoothSolver:  Solving for Ux, Initial residual = 0.000157922, Final residual = 1.01677e-05, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.000172252, Final residual = 9.70874e-06, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.000168408, Final residual = 1.10209e-05, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.00286416, Final residual = 6.07327e-05, No Iterations 3
time step continuity errors : sum local = 0.000134333, global = -1.86372e-05, cumulative = 0.176384
smoothSolver:  Solving for epsilon, Initial residual = 5.02636e-05, Final residual = 1.76898e-06, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.000214597, Final residual = 1.16413e-05, No Iterations 2
ExecutionTime = 3003.18 s  ClockTime = 3019 s
Code:
MRF1
{
    type            MRFSource;
    active          true;
    selectionMode   cellZone;
    cellZone        innerCylinder;

    MRFSourceCoeffs
    {
	nonRotatingPatches ();
        origin      (0 0 0);
        axis        (0 1 0);
        omega       31.4159;
    }
}
Attached Images
File Type: jpg paper's mesh.jpg (93.1 KB, 102 views)
File Type: jpg 11.jpg (94.0 KB, 124 views)
sharonyue is online now   Reply With Quote

Old   February 18, 2014, 05:30
Default Are there some news?
  #2
New Member
 
Join Date: Feb 2010
Posts: 7
Rep Power: 0
Dothan is on a distinguished road
Hi,
I found a similar result issue in a different geometry MRF problem.
Are there some news?

Dothan
Dothan is offline   Reply With Quote

Old   March 10, 2014, 12:22
Default
  #3
Member
 
Pengchuan Wang
Join Date: Nov 2012
Location: Michigan USA
Posts: 58
Rep Power: 14
pechwang is on a distinguished road
Hi everyone,

Same thing happens to me as well. Different size of MRFzones gives different results.
And for me, none of these gives me the same results as the analytical solution.
pechwang is offline   Reply With Quote

Old   March 11, 2014, 03:44
Default
  #4
New Member
 
Join Date: Feb 2010
Posts: 7
Rep Power: 0
Dothan is on a distinguished road
Hi pechwang,
Do you have a turbomachinery analytical solution to compare CFD solutiuons?
I know only Eckardt Impeller.

Dothan
Dothan is offline   Reply With Quote

Old   March 11, 2014, 11:06
Default
  #5
Member
 
Pengchuan Wang
Join Date: Nov 2012
Location: Michigan USA
Posts: 58
Rep Power: 14
pechwang is on a distinguished road
Hi Dothan,

No, the analytical solution is for lubrication model. Unfortunately, MRF cannot give me correct results.

Thanks,
Pengchuan
pechwang is offline   Reply With Quote

Old   May 29, 2015, 08:00
Default
  #6
Senior Member
 
M. Montero
Join Date: Mar 2009
Location: Madrid
Posts: 155
Rep Power: 17
be_inspired is on a distinguished road
Same problem here. When MRFZone limit is near walls ( sometimes it is required because stator zone is quite near) the pressure field around the walls is modified. The pressure lines are modified to accommodate its shape to the limit of the MRF Zone so suction lobe is affected and so all the aerodynamics of the wall (wake,....). The solution is wrong.

I do not have AMI interfaces so the solution proposed does not work.
Did you start a bug report?
be_inspired is offline   Reply With Quote

Old   August 9, 2019, 05:25
Default
  #7
Senior Member
 
Kmeti Rao
Join Date: May 2019
Posts: 145
Rep Power: 8
Krao is on a distinguished road
Hi all,

I am also confused about the size of MRF region. With all the parameters constant, each time I change my MRF region size I get different results. I have made several simulations in this regard and found out that, whenever MRF domain is very close to upper part of the propeller, It usually gives wrong results and wrong representation of flow. A domain size of around 1.1D width and 0.3-0.5D height have given good results, which are closer to the experimental values. Very small MRF regions have given poor results.

Because of having two MRF region closer to each other I am unable to find a proper solution to the above problem of domain sizing. It would be great if anyone of you can share your experience with MRF domain size and how you overcome this problem.

Thank you
Krao
Krao is offline   Reply With Quote

Old   February 10, 2021, 11:52
Default
  #8
New Member
 
Join Date: Oct 2019
Posts: 13
Rep Power: 7
anthony761 is on a distinguished road
I have observed the same with centrifugal compressor calculations on a commercial software. I saw that with increasing rotational domain (downstream of the impeller), the compression ratio increases.
anthony761 is offline   Reply With Quote

Old   May 28, 2022, 08:50
Default
  #9
Member
 
Wolfram Schneider
Join Date: Jan 2018
Location: Germany
Posts: 57
Rep Power: 8
Wolfram is on a distinguished road
well, okay ... I am not the only one facing this issue....

Did anyone get to any solution?
Wolfram is offline   Reply With Quote

Old   May 30, 2022, 10:13
Default
  #10
Member
 
Roland
Join Date: Mar 2009
Location: Netherlands
Posts: 93
Rep Power: 17
sylvester is on a distinguished road
See this old thread and the excellent explanation of Eugene here: Should an empty MRF zone really induce this much cross-flow?
sylvester is offline   Reply With Quote

Old   May 31, 2022, 05:46
Default
  #11
Member
 
Wolfram Schneider
Join Date: Jan 2018
Location: Germany
Posts: 57
Rep Power: 8
Wolfram is on a distinguished road
Quote:
Originally Posted by sylvester View Post
See this old thread and the excellent explanation of Eugene here: Should an empty MRF zone really induce this much cross-flow?

This is the answer. Although the MRF-Approach is not really suitable for my case, the post is calming me at least searching for my mistake in the setup.
Thank you!
Wolfram is offline   Reply With Quote

Reply


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
rasInterFoam appears to give erroneous results for “pd”. richard_kenny OpenFOAM Bugs 6 August 2, 2013 18:12
Must periodic oscillation give periodic results? zonexo Main CFD Forum 6 May 13, 2007 16:36


All times are GMT -4. The time now is 18:16.