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

[chtMultiRegionSimpleFoam] Control MRF inside fluid region

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 27, 2016, 12:37
Default [chtMultiRegionSimpleFoam] Control MRF inside fluid region
  #1
Member
 
Huynh Phong Thanh
Join Date: Aug 2013
Location: Ho Chi Minh City
Posts: 64
Rep Power: 4
hiuluom is on a distinguished road
I would like to simulation conjugate heat transfer fluid and solid with a zone MRF in fluid. The mesh was built from STL SHM with 4 zones. This is the of my solve.
  1. The problem:
    When using chtMultiRegionSimpleFoam solver, I must to split mesh into region including MRF zone and other zone. So I need to merge fluid contain MRF zone.
  2. Solution
    My method is to get cell zone using mergeMeshs because splitMeshRegions will divide all zone from cell zone.
    Here is example description:
    Model:

    Front view:



    Side view:



    Iso view:



    The model includes solid1, solid2, fluid and MRF zone
  3. Command:
    Split mesh into regions with the name v1(solid1), v2(solid1), v3(fluid), v4(MRFzone):

    Code:
    splitMeshRegions -cellZonesOnly -overwrite -useFaceZones
    I will combine v4 to v3. In v3 and v4 folder in constant directory, I make sub-directory 0, constant, system inside v3, v4 and then merge mesh. Remember that in system folder must have controlDict, fvSolution and fvScheme file.

    Quote:
    mergeMeshes v3 v4 -overwrite
    The last step of the mesh genesis is renaming of the boundary, changing patch type. Because boundary of v3 will have mappWall between v3 and v4. Changing this wall to cyclicAMI type. This managed by the utility createPatch.

    Code:
    patches
    (
        {
            name            v4_ami1;
            patchInfo       
            { 
                type cyclicAMI; 
                neighbourPatch v4_ami2;
            }
            constructFrom   patches;
            patches         ( default-interior_v3_to_v4 );
        }
        {
            name            v4_ami2;
            patchInfo       
            { 
                type cyclicAMI; 
                neighbourPatch v4_ami1; 
            }
            constructFrom   patches;
            patches         ( default-interior_v4_to_v3 );
        }
    )

The result seemly good, but I am not sure the cyclicAMI from face mappWall type is the same which make from MRF zone.

I will show result of my test model later
Attached Images
File Type: png frontview.png (4.6 KB, 45 views)
File Type: png isoview.png (8.8 KB, 45 views)
File Type: png sideview.png (1.7 KB, 45 views)
hiuluom is offline   Reply With Quote

Old   June 28, 2016, 22:03
Default
  #2
Member
 
Huynh Phong Thanh
Join Date: Aug 2013
Location: Ho Chi Minh City
Posts: 64
Rep Power: 4
hiuluom is on a distinguished road
Hi everybody,

I show my test case with fluid, solid, fan and propeller. The velocity of propeller in my guess.

This is domain:


Inside solid zone has water volume and fan.



Rotating velocity of propeller



water velocity inside solid make by fan



Finally, the path line and temperature of simulation



My question: Do the method to solve MRF with chtMultiregionSimpleFoam and the result test model as I mention above correct or incorrect for this model?
If anybody suggest run MRF with multi region better, you can suggest me other ways.

Best,
Thanh
Attached Images
File Type: png geo1.png (18.2 KB, 36 views)
File Type: png geo3.png (127.4 KB, 36 views)
File Type: png u1.png (82.6 KB, 36 views)
File Type: png u2.png (86.2 KB, 36 views)
File Type: png t1.png (153.7 KB, 36 views)

Last edited by hiuluom; June 29, 2016 at 10:46.
hiuluom is offline   Reply With Quote

Old   June 30, 2016, 05:14
Default
  #3
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,397
Blog Entries: 5
Rep Power: 24
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Dear Thanh, where do you have your MRF ? First I have to mention that I am not familiar with SRF/MRF but in your case you have a propellor that is rotating by a special speed or is there your MRF ? In any case, the streamlines you get are stationary but should be instationary due to the fact that the fan is rotating and "cutting" through the fluid.

I tried to understand the geometry but I am sorry, do not get the point. Can you explain it again?
__________________
Best regards,
Tobias Holzmann

Some interesting OpenFOAM tutorials and videos on www.Holzmann-cfd.de
Tobi is online now   Reply With Quote

Old   June 30, 2016, 10:44
Default
  #4
Member
 
Huynh Phong Thanh
Join Date: Aug 2013
Location: Ho Chi Minh City
Posts: 64
Rep Power: 4
hiuluom is on a distinguished road
Dear Tobi,

Thanks your reply.

Sorry my description did not detail. I update my explanation:

This is MRF zone as well as cyclicAMI



After running splitMeshRegions, I have 4 zones: Air, MRF_zone (including propeller), water and solid.
I merge the MRF_zone into Air. So in the Air and MRF_zone sub-directory in constant directory, I make folder 0, constant and system, respectively.
In the Air sub-directory I copy polyMesh into constant and rename polyMesh outside to polyMesh.old. Similarly for the MRF_zone.
And then I start merge 2 meshes by command mergeMeshe. When finish, I copy polyMesh directory in constant/Air/constant/polyMesh to constant/Air/.

Finally, I set createPatchDict file for Air in system/Air to change the boundary mapedwall between Air and MRF_zone to cyclicAMI and run topoSet for MRF_zone also.

I upload my case on Microsoft onedrive, you can download it at:https://1drv.ms/u/s!AkLk2AeMilR1gtZ73idBPAHaSBbRXg

I think the operation is complicated. Hope you understand it.

Best regards,
Thanh
Attached Images
File Type: jpg model.JPG (50.8 KB, 26 views)
hiuluom is offline   Reply With Quote

Old   June 30, 2016, 11:07
Default
  #5
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,397
Blog Entries: 5
Rep Power: 24
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Just a question - do you use AMI or MRF now ?
__________________
Best regards,
Tobias Holzmann

Some interesting OpenFOAM tutorials and videos on www.Holzmann-cfd.de
Tobi is online now   Reply With Quote

Old   June 30, 2016, 11:30
Default
  #6
Member
 
Huynh Phong Thanh
Join Date: Aug 2013
Location: Ho Chi Minh City
Posts: 64
Rep Power: 4
hiuluom is on a distinguished road
Dear Tobi,

What does it mean your question?

In my experience if I want to run a propeller simulation, I make MRF zone as well as cyclicAMI from face MRF. Do my understand miss?
hiuluom is offline   Reply With Quote

Old   June 30, 2016, 12:03
Default
  #7
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,397
Blog Entries: 5
Rep Power: 24
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Yes normally you are wrong.

AMI (arbitrary mesh interface) are the patches that connects a rotating cell zone with a stationary cellZone and hence you are able to rotate the rotating cell zone (where your propeller is in). The AMI just interpolates and maps values from one side of the baffle to the other one.

MRF is a model to avoid moving meshes and make life easier. See here: tutorials/compressible/rhoPimpleFoam/ras/mixerVessel2D

That was the reason why I was confused about your MRF and AMI. So finally you just added the motion libraries to the cht solver (actually you changed the derivative of the fvMesh to become a non-const object).

If you did this, I hope you also have all flux corrections in your equations because due to the rotation, you get mesh-fluxes that have to be corrected.
__________________
Best regards,
Tobias Holzmann

Some interesting OpenFOAM tutorials and videos on www.Holzmann-cfd.de
Tobi is online now   Reply With Quote

Old   June 30, 2016, 12:38
Default
  #8
Member
 
Huynh Phong Thanh
Join Date: Aug 2013
Location: Ho Chi Minh City
Posts: 64
Rep Power: 4
hiuluom is on a distinguished road
Hi Tobi,

I always run propeller with MRFzone. In addtions I also get cyclicAMI and the result has been good. I based on a tutorials/incompressible/pimpleDymFoam/propeller to set up case.

As you told me, if I simulate with one fluid no other zones, do I only need to add MRFzone name in fvOptions file?

Otherwise I make cyclicAMI and do not need use fvOptions file, but how to give the rotate velocity?

Best regards,
Thanh

Last edited by hiuluom; June 30, 2016 at 21:38.
hiuluom is offline   Reply With Quote

Old   June 30, 2016, 21:38
Default
  #9
Member
 
Huynh Phong Thanh
Join Date: Aug 2013
Location: Ho Chi Minh City
Posts: 64
Rep Power: 4
hiuluom is on a distinguished road
I saw difference cyclicAMI and MRF zone in fvOptions file that is zone active or inactive. "cyclicAMI and do not need use fvOptions file" that is my stupid questions
Thanks Tobi.
hiuluom is offline   Reply With Quote

Old   July 1, 2016, 03:25
Default
  #10
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,397
Blog Entries: 5
Rep Power: 24
Tobi will become famous soon enoughTobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,


  • if you are using a rotating mesh (mesh is moving in OpenFOAM), you have to define a cellzone that is rotating. To couple the rotating and stationary mesh, you have to need AMI boundary conditions.
  • If you have rotation somewhere and don't want to move the mesh or lets say, you can use the SRF or MRF approach, use it because here you do not need any mesh motion. Therefore you do not have to take care about additional mesh fluxes, interpolations between rotating and stationary mesh and so forth ... it is a model (but I never used it, hence I am not so familiar with that).
If you build a cellZone around your rotor and you really move this cellZone with dynamic mesh, you need AMI and hence there is no MRF approach.

If you somehow model the rotation with the MRF/SRF approach, there is no physical mesh motion in your case (you model the motion but the mesh stays the same during the calculation).
__________________
Best regards,
Tobias Holzmann

Some interesting OpenFOAM tutorials and videos on www.Holzmann-cfd.de
Tobi is online now   Reply With Quote

Old   July 2, 2016, 23:47
Default
  #11
Member
 
Huynh Phong Thanh
Join Date: Aug 2013
Location: Ho Chi Minh City
Posts: 64
Rep Power: 4
hiuluom is on a distinguished road
I tested with simpleFoam OF2.2.x. Both of models take from tutorial.
MRF and AMI 2D model ran with simpleFoam. The result pressure, velocity are not difference.
So I think that if I make AMI ratotating for solver not support dynamic mesh, the value do not effect when compares with MRF. In chtMultiRegionSimple does not support dynamic mesh, I can make AMI face as mention above but it doesn't have a mesh flux, interpolation between station and rotation mesh.


Here is my setup case fvOptions file:

MRF:

Code:
MRF1
{
    type            MRFSource;
    active          true;
    selectionMode   cellZone;
    cellZone        rotor;

    MRFSourceCoeffs
    {
        origin      (0 0 0);
        axis        (0 0 1);
        omega       104.72;
    }
}
AMI:

Code:
MRF1
{
    type            MRFSource;
    active          true;
    selectionMode   cellZone;
    cellZone        rotor;
    
    MRFSourceCoeffs
    {
    	nonRotatingPatches (AMI1 AMI2);
        origin      (0 0 0);
        axis        (0 0 1);
        omega       104.72;
    }
}
You can review in the following figures.
Attached Images
File Type: jpg glyph.JPG (73.9 KB, 6 views)
File Type: jpg p.JPG (46.7 KB, 6 views)
File Type: jpg u.JPG (55.0 KB, 6 views)
File Type: jpg wall.JPG (47.0 KB, 6 views)
hiuluom is offline   Reply With Quote

Old   July 3, 2016, 08:08
Default
  #12
Member
 
Huynh Phong Thanh
Join Date: Aug 2013
Location: Ho Chi Minh City
Posts: 64
Rep Power: 4
hiuluom is on a distinguished road
I tested between AMI and MRF 3D model, got moment coefficient with simpleFoam.
The result also is no difference.
Finally, I think if simpleFoam or chtMultiRegionSimpleFoam solver, which do not support dynamic mesh, can run with AMI or not.
Attached Images
File Type: jpg p-3d.JPG (58.4 KB, 5 views)
File Type: jpg u-3d.JPG (59.0 KB, 5 views)
File Type: jpg momentCoeff.JPG (25.9 KB, 3 views)
hiuluom is offline   Reply With Quote

Reply

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
Problem with an old Simulation FrankW CFX 3 February 8, 2016 05:28
Problem in setting Boundary Condition Madhatter92 CFX 12 January 12, 2016 05:39
chtMultiRegionFoam connection between solid and fluid region of heat exchanger ahab OpenFOAM 0 March 2, 2015 13:19
Difficulty in calculating angular velocity of Savonius turbine simulation alfaruk CFX 8 December 3, 2013 17:51
meshing F1 front wing Steve FLUENT 0 April 17, 2003 12:37


All times are GMT -4. The time now is 04:57.