|
[Sponsors] |
February 14, 2014, 16:50 |
SRFSimpleFoam vs SimpleFoam (MRF)
|
#1 |
New Member
Oscar Urbina
Join Date: Nov 2011
Location: Montréal
Posts: 16
Rep Power: 14 |
Hello Foamers,
i have a couple of questions about the rotating reference frame solvers in OpenFoam. i'm simulating a blade in half a cylinder, now the blade is static (represented by a force), and the entire domain has to move. i've been reading the threads and OpenFoamWiki, so from my understanding it would be the same if i use SRFSimpleFoam or if i use SimpleFoam. so my questions are: In SRFSimpleFoam i can define the axis of rotation but i'm unable to define the origin of my rotation, so is the origin going to be (0,0,0) all the time? or if there a way to set it? when i see SimpleFoam with the MRF option (MRFSimpleFoam), i can define the axis and the origin. would it be the equivalent of SRF if i define the entire domain as the cellzone that rotates? also when i use simpleFoam (domain with no obstacles in constant rotation), my inlet bc is fixedvalue (5,0,0), and my omega is 700rpm, but when i see my boundary in parafoam after a couple of iterations it becomes (0, something, something). i understand if it adds the rotational forces but why is setting my axial inflow to 0 when it's suppose to be fixed to 5? and in the internal field my axial velocity tends to 0 and gets negative values even. if someone could throw some light into this i would really appreciate it thanks Oscar |
|
March 10, 2014, 12:53 |
|
#2 |
Member
Pengchuan Wang
Join Date: Nov 2012
Location: Michigan USA
Posts: 58
Rep Power: 13 |
Hi Oscar,
If you put all the cells in the MRFzones, theoretically SRF and MRF will give you the same result. I did a very simple test before, that is using simpleFoam, SRFSimpleFoam and MRFSimpleFoam to run the same case. That means three different solvers solving the same problem, though the setup for each solver is different. Unfortunately, MRFSimpleFoam cannot give me the correct result. The results from simpleFoam and SRFimpleFoam are the same, and agree with my analytical solution. But the result from MRFSimpleFoam does not. And right now I still don't know the reason why the results are different. If you have any news, please update me. Thanks, Pengchuan |
|
March 11, 2014, 10:37 |
|
#3 |
New Member
Oscar Urbina
Join Date: Nov 2011
Location: Montréal
Posts: 16
Rep Power: 14 |
Hi Pengchuan,
I found a solution to my problem, first i discarded SRFSimpleFoam, the solution is too unstable at least for my case. i used SimpleFoam with the MRF option. now as i said before if you set all the cells in the rotationg frame it will crush your boundary conditions, i think it's because MRF is made for more than one frame, so if you put only one frame these kind of problems arise. the solution for me was to leave at least a one cell distance between the rotating frame and the BCs, this way your BCs are respected and the overall solution is not affected (it behaves as SRF). hope it helps, Oscar |
|
March 11, 2014, 10:41 |
|
#4 |
Member
Pengchuan Wang
Join Date: Nov 2012
Location: Michigan USA
Posts: 58
Rep Power: 13 |
Hi Oscar,
Thank you for your reply. Then your results are exactly the same? I mean the results with MRF and without MRF. For me, the results are different, and without MRF is better than the result with MRF. Thanks, Pengchuan |
|
March 16, 2014, 14:53 |
|
#5 |
New Member
Oscar Urbina
Join Date: Nov 2011
Location: Montréal
Posts: 16
Rep Power: 14 |
Hi,
i think they are suppose to be different in simpleFoam without MRF you only have the inertial frame, and in simpleFoam with MRF you have the inertial frame plus the rotating frame in other words you add the centrifugal and coriolisis forces. the only case in which both will have the same results is if you have a domain with no obstacle cause then the forces won't have any effect over the flow. as soon as you add an obtacle let say a force, the solutions will be different. you should read this article in openFoamWiki, here they explain how the MRF works. http://openfoamwiki.net/index.php/Se...RF_development hope it clears out some doubts, good luck Oscar |
|
March 16, 2014, 15:21 |
|
#6 |
Member
Pengchuan Wang
Join Date: Nov 2012
Location: Michigan USA
Posts: 58
Rep Power: 13 |
Hi Oscar,
Thank you for your reply. Yes, you are right, they use different method to solve the equation. However, though they use different methods, they are actually solve the same PDE, since they are simulating the same problem. No matter what method you use, the results should be the same. But right now, I got different results from them, that's what makes me confused. Thanks, Pengchuan |
|
May 13, 2015, 11:51 |
|
#7 | |
New Member
Join Date: May 2015
Posts: 7
Rep Power: 11 |
Hi Oscar,
in your first Message you wrote: Quote:
Cheers, Dennis |
||
December 3, 2015, 16:26 |
|
#8 |
Member
Rudolf Hellmuth
Join Date: Sep 2012
Location: Dundee, Scotland
Posts: 40
Rep Power: 14 |
I am also getting different results in SRF and MRF.
I am simulating a 3-D stirred tank, in which the propeller is on the bottom. In the SRF case, the liquid is stirred along the whole tank, whereas in the MRF case, the liquid is not getting stirred not too far away from the propeller. I have tried many ways to solve the MRF, messing up with different solution algorithms, interpolation schemes, and mesh grades, but the result is always very similar. On the other hand, in the SRF case, I did a test where I deleted the rotor, but it resulted in the whole tank equally stirring (besides the top and the bottom, of course). I guess, the rotor is not essential for the rotating the system in SRF, because the source terms are added everywhere. Interestingly, I did a test in MRF, where I elongated the rotating cellZone up to the top of the tank, but the result was the same as in all other MRF tests, instead of being similar to the SRF case. I am not too sure that SRF cannot be used in steering tanks, where the rotor is not as long as the tank itself. It seems to me that it cannot. I can't be sure, because I don't have access neither to a CFD expert nor to experimental data to validate my case. Could anyone clarify that? Attached is the SRF case without propeller. Boundary conditions are, side and bottom walls U (0 0 0), and top wall slip. |
|
February 24, 2017, 14:17 |
|
#9 | |
New Member
Karthikeyan
Join Date: Jan 2016
Posts: 1
Rep Power: 0 |
Quote:
I am currently working on SRFSimpleFoam, The case you have mentioned is right. cozz we have only one cellzone in this case and that whole set of 3D meshes keep rotating. The blade or propeller are kept in the way of the rotating fluid to simulate the results. In realtime, the blade rotates and the fluid is stationary, for simulation purpose, we keep the blade stable and rotate the fluid for calculation of effects. |
||
November 13, 2018, 18:18 |
|
#10 |
New Member
Luis Fernando
Join Date: May 2013
Posts: 14
Rep Power: 13 |
testing simpleFoam using MRF and leaving one cell outside to define the stationary zone, gives me non-reliable values. The y+ increases up to 5000!
|
|
Tags |
mrfsimplefoam, rotating frame, simplefoam, srfsimplefoam |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
interFoam vs. simpleFoam channel flow comparison | DanM | OpenFOAM Running, Solving & CFD | 12 | January 31, 2020 16:26 |
Simulation of Radial Fan with simpleFoam MRF | nash | OpenFOAM Running, Solving & CFD | 2 | November 5, 2015 11:12 |
simpleFoam MRF validation | mpl | OpenFOAM Running, Solving & CFD | 1 | August 19, 2014 13:17 |
Laminar simpleFoam and inviscid simpleFoam | herenger | OpenFOAM Running, Solving & CFD | 7 | July 11, 2013 07:27 |
Trying to run a benchmark case with simpleFoam | spsb | OpenFOAM | 3 | February 24, 2012 10:07 |