[Slaninář]

Hello,
I'm trying to simulate two MRFs in OpenFOAM 11. MRFProperties is:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "constant";
    object      MRFProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


MRF1
{
    cellZone    Rotor;
    active      yes;

    // Fixed patches (by default they 'move' with the MRF zone)
    nonRotatingPatches (Rotor_in Rotor_out );

    origin    (0 0 0);
    axis      (0 0 1);
    // omega     -151.843644923507;
    rpm       -1450;
}


MRF2
{
    cellZone    Mrot;
    active      yes;

    // Fixed patches (by default they 'move' with the MRF zone)
    nonRotatingPatches (Mrot_out Mrot_in );

    origin    (0 0 0);
    axis      (0 0 1);
    rpm       -1450;
}




// ************************************************************************* //

However, for incompressibleFluid solver, this error occurs:

/*---------------------------------------------------------------------------*\ 
  =========                 | 
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox 
   \\    /   O peration     | Website:  https://openfoam.org 
    \\  /    A nd           | Version:  11 
     \\/     M anipulation  | 
\*---------------------------------------------------------------------------*/ 
Build  : 11 
Exec   : foamRun 
Date   : Apr 27 2024 
Time   : 13:06:18 
Host   : "slaninar-extreme4" 
PID    : 105828 
I/O    : uncollated 
Case   : /media/slaninar/0f1eafad-3b93-4c73-9544-109935e24a87/Programy/OpenFOAM/slaninar-10/run_org/PumpOFI/OpenFOAM_1 
nProcs : 1 
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). 
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) 
allowSystemOperations : Allowing user-supplied system call operations 

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 
Create time 

Create mesh for time = 0 

patchToPatch: Calculating couplings between 2660 source faces and 10952 target faces 
    Source min/average/max coverage = 0.860403/0.99973/1 
    Target min/average/max coverage = 0.90105/0.999733/1 
    Source average openness/error/depth/angle = 1.34718e-16/2.43049e-15/3.25552e-16/0 
    Source max openness/error/depth/angle = 7.97206e-16/2.43812e-14/2.90293e-15/0 
    51018 couplings calculated in 0.296447s 
patchToPatch: Calculating couplings between 10796 source faces and 20636 target faces 
    Source min/average/max coverage = 0.623984/0.999557/1 
    Target min/average/max coverage = 0.225696/0.999953/1 
    Source average openness/error/depth/angle = 1.21878e-16/0.000701521/0.000723702/0.126176 
    Source max openness/error/depth/angle = 6.64406e-16/0.38887/0.158765/24.1916 
    126220 couplings calculated in 0.654447s 
patchToPatch: Calculating couplings between 33155 source faces and 21924 target faces 
    Source min/average/max coverage = 0.5/0.999924/1 
    Target min/average/max coverage = 0.726963/0.999915/1 
    Source average openness/error/depth/angle = 1.35812e-16/0.00186038/0.0220941/0.75074 
    Source max openness/error/depth/angle = 1.83427e-15/0.385519/0.480225/32.9898 
    214126 couplings calculated in 1.56214s 
fvMeshStitcher: Connecting 
    Cell min/avg/max openness = 0/4.93197e-17/3.68435e-16 
Selecting solver incompressibleFluid 
Selecting viscosity model constant 
Selecting turbulence model type RAS 
Selecting RAS turbulence model realizableKE 
RAS 
{ 
    model           realizableKE; 
    turbulence      on; 
    printCoeffs     on; 
    A0              4; 
    C2              1.9; 
    sigmak          1; 
    sigmaEps        1.2; 
} 

Creating MRF zone list from MRFProperties 
    creating MRF zone: MRF1 
    - selecting cells using cellZone Rotor 
    - selected 2354137 cell(s) with volume 0.00201508 
    creating MRF zone: MRF2 
    - selecting cells using cellZone Mrot 
    - selected 79156 cell(s) with volume 0.0001888 

SIMPLE: Convergence criteria found 
        p: tolerance 1e-06 
        U: tolerance 1e-06 
        "(k|omega)": tolerance 1e-06 


SIMPLE: Operating solver in steady-state mode with 1 outer corrector 
SIMPLE: Operating solver in SIMPLE mode 



Starting time loop 

Constructing face momentum Uf 


--> FOAM FATAL ERROR:  
MRFZoneName not specified 

    From function const Foam::MRFZone& Foam::MRFPatchField::MRFzone(const Foam::objectRegistry&) const 
    in file cfdTools/general/MRF/derivedFvPatchFields/MRFPatchField/MRFPatchField.C at line 78. 

FOAM exiting

If only one MRF zone is in MRFProperties everything works.
I found this thread but it is old and it doesn't solve my problem.
Does anybody know how to solve it?

[Marcelinus]

Hi there,



I got the same error a couple days ago and this was my fix:


In the velocity boundary file you need to specify the name of the MRF zone if there are multiple. Like this:

MRFRegion1
            {
                    type                 MRFnoSlip;
                    MRFZoneName    zone1;
            }

          MRFRegion2
            {
                    type                  MRFnoSlip;
                    MRFZoneName     zone2;
             }

You can get the name of the zone from the MRFProperties file. Hope this helps.

zone1
         {

                   cellZone    cell_1;
                   active       yes;

                   //nonRotatingPatches ();

                   origin    (-5.480-0.3003.000);
                   axis      (-0.346410.2000);

                   //omega     12.566370;
                   rpm            365;
          }

          zone2
          {
                   cellZone        cell_2;
                   active           yes;

                   //nonRotatingPatches ();

                   origin      (4.7752550.7071074.085786);
                   axis        (0.367423-0.2121320.424264);

                   //omega     12.566370;
                   rpm            69;
          }

[Slaninář]

It works. Now incompressible solver calculates two MRFs and results look believable. Thank you very much!