simpleFoam+MRF vs pipleFoam+dynamicMesh
2 Attachment(s)
Dear all,
I have a question regarding the difference between the two approaches, MRF and dynamic mesh. I’m referring to the propeller tutorial: Code:
/foam/tutorials/incompressible/pimpleFoam/RAS/propeller However I modified the topoSet dictionary to split the patch of rotating zone into inlet outlet and cylinder. Then I tried to evaluate the Thrust of the propeller using the two approaches: 1) pimpleFoam +dynamicMeshDict: this is an unsteady analysis that imposes the rotation of the innerCylinderSmall cellZone using the dynamicMeshDict file. The solver is pimpleFoam and I used the same setting of the tutorial case. After 0.3s of physical time the Thrust (y direction) is 369.55N Residuals are 6E-3 for pressure and 1E-3 for velocities. 2) simpleFoam + MRFProperties: a steady analyses. I added the MRFProperties dictionary in place of the dynamicMeshDict. Same rotational velocity, same rotation axis. I modified the fvSolution and fvSchemes accordingly with the simpleFoam solver. After 4000 iterations the Thrust (y direction) is 940.27N. Residuals are 6E-4 for pressure and 6E-5 for velocities. Both the analysis seems to converge (to a stable integral value of Thrust); however the results are very different. Where am I wrong? Also the streamlines are different: in pimpleFoam+dynamicMesh the U field represents the absolute velocity, however what is the U field of the MRF zone in the simpleFoam case? These streamlines don’t seem neither the absolute or relative velocities.. Attached is the postprocessing in Paraview, you can notice how both the pressure field and the streamlines are different. Thanks, WhiteW |
2 Attachment(s)
The error seems related to cyclicAMI patches when using MRF.
I run an analysis with simpleFoam without MRF untill it converges (first pictures) then I added the MRFProperties file to constant folder with a low rotational velocity. And the results after 10 iterates in on the second picture. The pressure field completely changes in 10 iterates, even if the rotational velocity of MRF is 0.001 rad/s... How is it possible? Here is my MRFProperties: Code:
MRF1 And this is the boundary file: Code:
13 WhiteW |
I'm trying to do the same. I I have ran the propeller case and now I'm trying to run it with simpleFoam. Have you managed to run it correctly?
|
Hello,
If I am right in your simpleFoam + MRF setup you have to tell to OpenFOAM that your cyclic patch are non rotating. so : Quote:
Best regards |
In which file do I have to set this command?
Quote:
|
All times are GMT -4. The time now is 18:09. |