CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Potentialfoam requesting surfaceScalarField phi

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 14, 2013, 17:53
Default Potentialfoam requesting surfaceScalarField phi
  #1
Member
 
David
Join Date: Aug 2012
Posts: 48
Rep Power: 5
GM_XIII is on a distinguished road
Hi mates,

I am trying to initalize my simplefoam problem with potentialfoam and now i am getting this error:

Code:
--> FOAM FATAL ERROR: 

    request for surfaceScalarField phi from objectRegistry region0 failed
    available objects of type surfaceScalarField are

0
(
)


    From function objectRegistry::lookupObject<Type>(const word&) const
    in file /opt/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 131.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/potentialFoam"
#3  Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const& Foam::objectRegistry::lookupObject<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >(Foam::word const&) const in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4  Foam::surfaceInterpolationScheme<Foam::Vector<double> >::addMeshConstructorToTable<Foam::upwind<Foam::Vector<double> > >::New(Foam::fvMesh const&, Foam::Istream&) in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5  Foam::surfaceInterpolationScheme<Foam::Vector<double> >::New(Foam::fvMesh const&, Foam::Istream&) in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/potentialFoam"
#6  
 at potentialFoam.C:0
#7  
 in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/potentialFoam"
#8  __libc_start_main in "/lib/libc.so.6"
#9  
 in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/potentialFoam"
Aborted
Seems it is looking for the surfaceField phi. Could someone help me with this? Thanks in advance.

I post my fvsolution:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.com               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
    solver          GAMG;
    tolerance       1e-06;
    relTol          0.01;
    smoother        DIC;
    cacheAgglomeration true;
    nCellsInCoarsestLevel 500;
    agglomerator    faceAreaPair;
    mergeLevels     1;
}

U
{
    solver          PBiCG;
    preconditioner  DILU;
    tolerance       1e-06;
    relTol          0.0;
}

k
{
    solver          PBiCG;
    preconditioner  DILU;
    tolerance       1e-06;
    relTol          0.0;
}

epsilon
{
    solver          PBiCG;
    preconditioner  DILU;
    tolerance       1e-06;
    relTol          0.0;
}

}
//version OF 1.7.1:
//SIMPLE
//{
//    nNonOrthogonalCorrectors 6;
//    convergence 5e-3;
//}
//version OF 2.1.1:
SIMPLE
{
    nNonOrthogonalCorrectors 8;

    residualControl
    {
        p               5e-3;
        U               1e-3;
        "(k|epsilon)"   1e-3;
    }
}
//version OF 1.7.1:
//relaxationFactors
//{
//    p               0.05;
//    U               0.05;
//    k               0.01;
//    epsilon         0.01;
//}
//version 2.1.1:
relaxationFactors
{
    fields
    {
        p               0.3;
    }
    equations
    {
        U               0.7;
        k               0.1;
        epsilon         0.1;
    }

}
potentialFlow
{
    nNonOrthogonalCorrectors 10;
}
// ********************1***************************************************** //
and fvschemes
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.com               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default steadyState;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss upwind phi;
    grad(U)         Gauss upwind phi;
}

divSchemes
{
    default         Gauss linear;
    div(phi,U)      Gauss upwind;
    div(phi,k)      Gauss upwind;//Gauss vanLeer;
    div(phi,epsilon) Gauss upwind;//Gauss vanLeer;
// version 1.7.1:
//   div((nuEff*dev(grad(U).T()))) Gauss linear;
// version 2.1.1:
    div((nuEff*dev(T(grad(U)))))  Gauss upwind phi;
}

laplacianSchemes
{
    default         Gauss upwind phi corrected;
    //laplacian(nuEff,U) Gauss linear corrected;
    //laplacian(1,p) Gauss linear corrected;
    //laplacian((1|A(U)),p) Gauss linear corrected;
    //laplacian(DkEff,k) Gauss linear corrected;
    //laplacian(DepsilonEff,epsilon) Gauss linear corrected;
}

interpolationSchemes
{
    default         upwind phi;
    interpolate(U)  upwind phi;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p;
}



// ************************************************************************* //
GM_XIII is offline   Reply With Quote

Old   November 15, 2013, 03:15
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 13
Bernhard is on a distinguished road
I think you should change this
Code:
 interpolate(U)  upwind phi;
To something like this
Code:
 interpolate(U)  linear;
The problem is that simpleFoam will try to create phi from U by an interpolation. Because you want to discretize it using an upwind scheme, you need phi again. By a linear scheme, that is not necessary. I don't think the interpolate(U) is used anywhere else in the code.
Bernhard is offline   Reply With Quote

Old   November 15, 2013, 16:57
Default
  #3
Member
 
David
Join Date: Aug 2012
Posts: 48
Rep Power: 5
GM_XIII is on a distinguished road
Quote:
Originally Posted by Bernhard View Post
I think you should change this
Code:
 interpolate(U)  upwind phi;
To something like this
Code:
 interpolate(U)  linear;
The problem is that simpleFoam will try to create phi from U by an interpolation. Because you want to discretize it using an upwind scheme, you need phi again. By a linear scheme, that is not necessary. I don't think the interpolate(U) is used anywhere else in the code.
Seems logic, thank you mate. I didn't realize that simplefoam take the same shemes in case you dont specify them
GM_XIII 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
MSHArequest for surfaceScalarField phi from objectRegistry msha OpenFOAM Running, Solving & CFD 11 March 21, 2014 04:35
mixerVesselAMI2D's mass is not balancing sharonyue OpenFOAM Running, Solving & CFD 6 June 10, 2013 09:34
Summing phi over all faces in patch NewtoFOAM OpenFOAM Programming & Development 2 May 8, 2013 05:37
Sample phi potentialFoam doubtsincfd OpenFOAM 0 May 16, 2011 03:26
Turbulence Model phi vs phi_ doug OpenFOAM Running, Solving & CFD 4 November 10, 2009 05:33


All times are GMT -4. The time now is 21:33.