[farzadmech]

Dear friends
I have changed pitzdaily example(pimpleFoam solver), and I want to add momentum source to that.
I use the topoSetDict as below to create the desired zone;

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      topoSetDict.1;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

actions
(

    {
        name    farzaaadSulzerCellSet;
        type    cellSet;
        action  new;     //remove;  // new;
        source  boxToCell;

        sourceInfo
        {
            box (0.096 0.10 0.0357) (0.404 0.20 0.0643);
        }
    }
    
      {
        name    farzaaadSulzer;
        type    cellZoneSet;
        action  new; // remove; // new;
        source  setToCellZone;
        sourceInfo
        {
            set farzaaadSulzerCellSet;
        }
    } 
    
    
    
);

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

and I used fvOption as below;

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

momentumSource
{

    type            meanVelocityForce;

    active          on;            //on/off switch

    //cellSet       sourceSet;

 

    meanVelocityForceCoeffs

    {
      timeStart   0.;
      duration    1000.;

        selectionMode   cellZone;        //cellSet;       //cellZone;
        cellZone         farzaaadSulzer;  //cellZone        farzadSulzer; 

        fields                  (Uc);

        Ubar            ( 5.0 0.0 0.0);               // ( 0. -1.0 0.);

        relaxation  1.0;

    }

}

What I am doing wrong in my simulation(momentum source did not do anything)? I have attached the figure of my domain and also the momentum source inside it to this post.


Thanks,
Farzad

[farzadmech]

in the "incompressible/pimpleFoam/LES/channel395/", we have;

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

momentumSource
{
    type            meanVelocityForce;

    selectionMode   all;

    fields          (U);
    Ubar            (0.1335 0 0);
}


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

But why it does not work in my pimpleFoam example?


Farzad

[mAlletto]

any special reason you wrote fields

 fields (Uc);

There are no fields Uc you solve for as far I know in OpenFoam