[hoemmaria]

Hello!

I am using OpenFoam 4.1 and I am trying to implement an actuator disk. The problem I am having is that velocity and pressure is not calculated. At the first time step it says "No iterations = 0". The initial velocity and pressure files 0/U and 0/p are copied to all the following time steps. This is obviously not what I want. Here are my files and process:

I want to use the simpleFoam solver and I introduce the actuation disk in a fvOptions file as follows:

constant/fvOptions:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

disk1
{
    type            actuationDiskSource;
    active          on;

    actuationDiskSourceCoeffs
    {
        fields      (U);

        selectionMode   cellSet;
        cellSet         actuationDisk1;
        diskDir         (1 0 0);    // Orientation of the disk
        Cp              0.4;
        Ct              0.5;
        diskArea        0.16;
        upstreamPoint   (4.9 0.5 0.5);
    }
}

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

0/p :

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
      type          zeroGradient;

    }
    outlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    walls
    {
        type            zeroGradient;
    }
}

0/U:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            pressureInletVelocity;
        value           uniform (0 0 0);
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0);
        value           uniform (0 0 0);
    }
    walls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}

I am using the simpleFoam solver and I have not made any changes to the pEqn.H, UEqn.H, simpleFoam.C or createFields.H.

My blockMeshDict only includes inlet, outlet and walls.

I use topoSet with the following file:

system/topoSetDict

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      topoSetDict;
}

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

actions
(
    // actuationDisk1
    {
        name    actuationDisk1CellSet;
        type    cellSet;
        action  new;
        source  boxToCell;
        sourceInfo
        {
            box (4.9 0.1 0.1) (5.1 0.9 0.9);
        }
    }
    {
        name    actuationDisk1;
        type    cellZoneSet;
        action  new;
        source  setToCellZone;
        sourceInfo
        {
            set actuationDisk1CellSet;
        }
    }
);

Both epsilon and k are computed as normal when run the code, but nothing happens to neither the velocity nor the pressure. I have also tried the turbineSiting tutorial where the actuationDiskSource is being used, and it runs without problem. But I do not seem to find the cause of my case failing compared to the turbineSiting tutorial.

I any of you have tried something similar before or have any suggestions to what I might be doing wrong, I am most grateful for help.


Kind regards,
Maria Hoem