CFD Online Discussion Forums - Problems on implementing actuation disk on piso solver

Hi fellows,

I am tring to implement the actuation disk to the pisoFoam solver. Currently all I know about the way to do it is learned from the tutorial package “turbine siting”. However, i just want to implement it under a simple channel flow, and got the following problem on reading the disk properties.

The modified piso solver is like following:

Code:



  #include "fvCFD.H"

    #include "singlePhaseTransportModel.H"

    #include "turbulenceModel.H"

    #include "IObasicSourceList.H"

     

   

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

     

   

  int main(int argc, char *argv[])

    {

        #include "setRootCase.H"

        #include "createTime.H"

        #include "createMesh.H"

        #include "createFields.H"

        #include "initContinuityErrs.H"

     

   

      IObasicSourceList actuationDisks(mesh);

     

   

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

     

   

      Info<< "\nStarting time loop\n" << endl;

     

   

      while (runTime.loop())

        {

            Info<< "Time = " << runTime.timeName() << nl << endl;

     

   

          #include "readPISOControls.H"

            #include "CourantNo.H"

     

   

          // Pressure-velocity PISO corrector

            {

                // Momentum predictor

     

   

              fvVectorMatrix UEqn

                (

                    fvm::ddt(U.oldTime(), U)

                  + fvm::div(phi, U) 

                + turbulence->divDevReff(U)

                );

     

   

              // Add resistance on the actuation disks

                actuationDisks.addSu(UEqn);

     

   

              UEqn.relax();

     

   

              if (momentumPredictor)

                {

                    solve(UEqn == -fvc::grad(p));

                }

     

   

              // --- PISO loop

                 ....

and the “sourceProperties” file and “makeZones” file is the same as the one in tutorial (changed size and location). When I run the new piso solver on my case, I got the following message:

Code:





  --> FOAM FATAL IO ERROR: 

    keyword type is undefined in dictionary "/home/frank/OpenFOAM/frank-2.0.0/run/case8/constant/sourcesProperties::disk1"

       

  file: /home/frank/OpenFOAM/frank-2.0.0/run/case8_ActuationDisk/cpiso/constant/sourcesProperties::disk1 from line 20 to line 33.

     

   

      From function dictionary::lookupEntry(const word&, bool, bool) const

        in file db/dictionary/dictionary.C at line 400.

       

  FOAM exiting

Am I missing some settings in my case? I cross checked with the tutorial case but couldn't find one.

Thanks a lot.

Regards
Frank