CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Problems on implementing actuation disk on piso solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 25, 2013, 00:36
Default Problems on implementing actuation disk on piso solver
  #1
New Member
 
Frank Yu
Join Date: Jun 2011
Location: Toronto, ON
Posts: 15
Rep Power: 5
frankyux is on a distinguished road
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
frankyux is offline   Reply With Quote

Old   September 28, 2013, 08:01
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 7,093
Blog Entries: 32
Rep Power: 70
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Greetings Frank,

The message is somewhat clear about the problem:
Quote:
keyword type is undefined in dictionary "/home/frank/OpenFOAM/frank-2.0.0/run/case8/constant/sourcesProperties::disk1"
The entry "type" is missing from the block "disk1" in the file "constant/sourcesProperties".

In OpenFOAM 2.0, an example can be found at tutorials/incompressible/windSimpleFoam/turbineSiting/constant/sourcesProperties.
In OpenFOAM 2.1, the example is at tutorials/incompressible/simpleFoam/turbineSiting/constant/sourcesProperties.
In OpenFOAM 2.2, the example is at tutorials/incompressible/simpleFoam/turbineSiting/system/fvOptions.

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   October 2, 2013, 23:06
Default
  #3
New Member
 
Frank Yu
Join Date: Jun 2011
Location: Toronto, ON
Posts: 15
Rep Power: 5
frankyux is on a distinguished road
Thanks Bruno. Problem solved.
I copied an example from an older version where 'type' becomes 'typemodel', that's why the error message comes up.
frankyux 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
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 15:45
Problems with unsteady transonic solver Frank Main CFD Forum 0 July 24, 2006 13:48
Solver Software (?) problems JvK CFX 5 August 9, 2002 13:33
CFX 5.5 Roued CFX 1 October 2, 2001 16:49
Solving 2-d problems by 3-d solver eddy Main CFD Forum 3 September 7, 2000 06:15


All times are GMT -4. The time now is 12:08.