[Fanfei]

Originally Posted by lucchini

Hi Susan,
you get the error because you try to assign to a scalar variable the product between a scalarField and a volScalarField.
You can get it working if you write:

volScalarField z = x*y*mag(U)

I hope it has been useful,
bye

Tommaso

Hi Tommaso:
I want to calculated volScalarField Shildes with volVectorField wallShearStress. The expression is Shildes=mag(wallShearStress)/(rho*g).

      volScalarField  rho
     (
      "rho",
      1000.0
     ); 

     volScalarField shileds
    (
       IOobject
       (
           "shileds",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::NO_WRITE
        ),
        mesh,
       mag(wallShearStress)/(rho*g)      
      );

as compiled, some errors occured: could you give me some advices.

 pimFoam.C:129:7: error: no matching function for call to ‘Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(const char [4], double)’
pimFoam.C:129:7: note: candidates are:
In file included from /home/administrator/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/GeometricField.H:601:0,
                 from /home/administrator/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/GeometricScalarField.H:38,
                 from /home/administrator/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/GeometricFields.H:34,
                 from /home/administrator/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/volFields.H:37,
                 from /home/administrator/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/surfaceInterpolationScheme.C:30,
                 from /home/administrator/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/surfaceInterpolationScheme.H:234,

[Fanfei]

anyone can help me?

[Fanfei]

Hi, Bruno Santos: I wantto add the wallshearstress to pimpleFoam solver.the code is:

            while (pimple.correct())
            {
                #include "pEqn.H"
            }

            if (pimple.turbCorr())
            {
                turbulence->correct();
            }
        }

      Info<<"1  step"<<endl;
      autoPtr<incompressible::RASModel> model;
      Info<<"2 step"<<endl;
      const volSymmTensorField Reff(model->devReff());
      Info<<"3 step"<<endl;

 //   #include "createPhi.H"
  //    autoPtr<incompressible::RASModel> model;
  //  singlePhaseTransportModel laminarTransport(U, phi);
  //    Info<<"3 step"<<endl;
  //  autoPtr<incompressible::RASModel> model
  //  (
  //      incompressible::RASModel::New(U, phi, laminarTransport)
  //  );
 
     label patchWallID=mesh.boundaryMesh().findPatchID("bottomWall");

       volVectorField wallShearStress     //定义剪切应力
       (
           IOobject
            (
              "wallShearStress",
              runTime.timeName(),
              mesh,
              IOobject::NO_READ,
              IOobject::AUTO_WRITE
            ),
            mesh,
            dimensionedVector
            (
              "wallShearStress",
              sqr(dimLength)/sqr(dimTime),
              vector::zero
            )
       );
      {  
        wallShearStress.boundaryField()[patchWallID] =
        (
           -mesh.Sf().boundaryField()[patchWallID]
           /mesh.magSf().boundaryField()[patchWallID]
        ) & Reff.boundaryField()[patchWallID];
       }

when compiled, there is no problem. but as run a test, there're a problem

 GAMG:  Solving for p, Initial residual = 0.0390428, Final residual = 6.86508e-08, No Iterations 25
time step continuity errors : sum local = 1.61558e-07, global = -3.77072e-08, cumulative = -2.71215e-05
smoothSolver:  Solving for epsilon, Initial residual = 0.187626, Final residual = 8.51315e-06, No Iterations 32
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 8.18223e-06, No Iterations 45
1  step
2 step


--> FOAM FATAL ERROR: 
object of type N4Foam14incompressible8RASModelE is not allocated

    From function Foam::autoPtr<T>::operator->()
    in file /home/administrator/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/autoPtrI.H at line 176.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  
 at ??:?
#3  
 at ??:?
#4  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#5  
 at ??:?
Aborted (core dumped)

I think it's may be same parameters undefined, could you give me some advise? Thanks

[wyldckat]

Hi Fei Fan,

I've moved your related posts above from the following threads:

• http://www.cfd-online.com/Forums/ope...ctorfield.html
• http://www.cfd-online.com/Forums/ope...ar-stress.html

I created this new thread in the programming sub-forum, since this is a programming related series of questions


OK, regarding the latest post, what is happening is simple:

1. The solver is telling you that you're trying to use a turbulence model that has not yet been allocated.
2. This was because of you using this code:
   

   autoPtr<incompressible::RASModel> model;
   const volSymmTensorField Reff(model->devReff());

   From a simple point of view, you only defined that the variable "model" is of type "autoPtr<incompressible::RASModel>". But you did not initialize the variable "model" with any values, therefore it's undefined, and in this case, not allocated.
3. If you look at the code commented out a few lines below, you'll see how the turbulence model is being initialized, namely these lines are all related:
   

   #include "createPhi.H"
   singlePhaseTransportModel laminarTransport(U, phi);
   autoPtr<incompressible::RASModel> model
   (
       incompressible::RASModel::New(U, phi, laminarTransport)
   );

4. My advice is that you do a bit more careful studying and training in C++ as a programming language, before stepping directly into OpenFOAM's source code.

Best regards,
Bruno

[Fanfei]

Originally Posted by wyldckat

Hi Fei Fan,

I've moved your related posts above from the following threads:

• http://www.cfd-online.com/Forums/ope...ctorfield.html
• http://www.cfd-online.com/Forums/ope...ar-stress.html

I created this new thread in the programming sub-forum, since this is a programming related series of questions


OK, regarding the latest post, what is happening is simple:

1. The solver is telling you that you're trying to use a turbulence model that has not yet been allocated.
2. This was because of you using this code:
   

   autoPtr<incompressible::RASModel> model;
   const volSymmTensorField Reff(model->devReff());

   From a simple point of view, you only defined that the variable "model" is of type "autoPtr<incompressible::RASModel>". But you did not initialize the variable "model" with any values, therefore it's undefined, and in this case, not allocated.
3. If you look at the code commented out a few lines below, you'll see how the turbulence model is being initialized, namely these lines are all related:
   

   #include "createPhi.H"
   singlePhaseTransportModel laminarTransport(U, phi);
   autoPtr<incompressible::RASModel> model
   (
       incompressible::RASModel::New(U, phi, laminarTransport)
   );

4. My advice is that you do a bit more careful studying and training in C++ as a programming language, before stepping directly into OpenFOAM's source code.

Best regards,
Bruno

Thanks Bruno:
You are right, I have solved it. regards.