CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   request for surfaceScalarField corrected from objectRegistry region0 failed (http://www.cfd-online.com/Forums/openfoam-solving/113547-request-surfacescalarfield-corrected-objectregistry-region0-failed.html)

immortality February 21, 2013 13:40

request for surfaceScalarField corrected from objectRegistry region0 failed
 
Code:

Starting time loop

Courant Number mean: 0 max: 0
deltaT = 6.711409396e-09
Time = 6.71141e-09

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1


--> FOAM FATAL ERROR:

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

3
(
weights
phi
deltaCoeffs
)


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

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  at SFCD.C:0
#3  Foam::surfaceInterpolationScheme<double>::addMeshConstructorToTable<Foam::LimitedScheme<double, Foam::SFCDLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr> >::New(Foam::fvMesh const&, Foam::Istream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4  Foam::surfaceInterpolationScheme<double>::New(Foam::fvMesh const&, Foam::Istream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
#5  Foam::fv::laplacianScheme<Foam::Vector<double>, double>::addIstreamConstructorToTable<Foam::fv::gaussLaplacianScheme<Foam::Vector<double>, double> >::New(Foam::fvMesh const&, Foam::Istream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6  Foam::fv::laplacianScheme<Foam::Vector<double>, double>::New(Foam::fvMesh const&, Foam::Istream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#7  Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#8  Foam::compressible::RASModels::kOmegaSST::divDevRhoReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#9
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
Aborted
Killing PID 4931
 PyFoam WARNING on line 232 of file /usr/local/lib/python2.7/dist-packages/PyFoam/Execution/FoamThread.py : Process 4931 was already dead

my fvScheme is:
Code:

laplacianSchemes
{
    default      Gauss SFCD corrected;//limited .5
    /*laplacian(muEff,U) Gauss linear corrected;
    laplacian(mut,U) Gauss linear corrected;
    laplacian(DkEff,k) Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
    laplacian(DREff,R) Gauss linear corrected;
    laplacian(DomegaEff,omega) Gauss linear corrected;
    laplacian((rho*(1|A(U))),p) Gauss linear corrected;
    laplacian(alphaEff,h) Gauss linear corrected;
    laplacian(k,T)  Gauss linear corrected; 
    laplacian(alpha,e) Gauss linear corrected;
    laplacian(alphaEff,e)  Gauss linear corrected;*/
}

interpolationSchemes
{
    default        SFCD;// cubicCorrection
}

snGradSchemes
{
    default      fourth;//corrected
}


immortality February 26, 2013 19:04

can someone tell me the possible cause?

KYPCK444 September 11, 2013 12:52

I have the same kind of problem.

I modified simpleFoam, on the basis of channelFoam, in order to add the term gradP in UEqn, so that I can use cyclic boundary conditions and simulate an infinite pipe flow.

I also changed some entries in fvSchemes (following these recommendations: http://www.cfd-online.com/Forums/ope...implefoam.html):

ddtSchemes
{
default backward;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear; not changed to upwind as suggested because it does not work at all
}

divSchemes
{
default none;
div(phi,U) Gauss upwind; changed from linear
div(phi,k) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div(phi,R) Gauss limitedLinear 1;
div(R) Gauss linear;
div(phi,nuTilda) Gauss limitedLinear 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss upwind corrected; changed from linear
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate(U) upwind; changed from linear
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}


but when I run my new solver, this is what I obtain:

--> FOAM FATAL ERROR:

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

4
(
weights
phi
deltaCoeffs
nonOrthDeltaCoeffs
)


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

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const& Foam::objectRegistry::lookupObject<Foam::Geometric Field<double, Foam::fvsPatchField, Foam::surfaceMesh> >(Foam::word const&) const in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#3 Foam::surfaceInterpolationScheme<double>::addMeshC onstructorToTable<Foam::upwind<double> >::New(Foam::fvMesh const&, Foam::Istream&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::surfaceInterpolationScheme<double>::New(Foam ::fvMesh const&, Foam::Istream&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleTransportModels.so"
#5 Foam::fv::laplacianScheme<Foam::Vector<double>, double>::addIstreamConstructorToTable<Foam::fv::ga ussLaplacianScheme<Foam::Vector<double>, double> >::New(Foam::fvMesh const&, Foam::Istream&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::fv::laplacianScheme<Foam::Vector<double>, double>::New(Foam::fvMesh const&, Foam::Istream&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#7 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#8 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#9 Foam::incompressible::RASModels::kEpsilon::divDevR eff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#10
in "/home/fsauser/OpenFOAM/fsauser-2.1.x/platforms/linux64GccDPOpt/bin/simplePipeFoam"
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12
in "/home/fsauser/OpenFOAM/fsauser-2.1.x/platforms/linux64GccDPOpt/bin/simplePipeFoam"
Aborted



I am new to openFOAM, and I'm not able yet to fully understand what could be happening with numerical schemes.

any hints? :)

thanks

immortality September 11, 2013 16:07

test adding phi in all laplacian schemes like it:
Code:

laplacian(nuEff,U) Gauss upwind corrected phi;

KYPCK444 September 12, 2013 06:07

Quote:

Originally Posted by immortality (Post 451209)
test adding phi in all laplacian schemes like it:
Code:

laplacian(nuEff,U) Gauss upwind corrected phi;

thank you very much for your reply,
unfortunately it does not work, the errors I get are exactly the same.

Anyway, I tried to further handle fvSchemes in such a way:

ddtSchemes
{
default backward;
}

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

divSchemes
{
default none;
div(phi,U) Gauss upwind phi;
div(phi,k) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div(phi,R) Gauss limitedLinear 1;
div(R) Gauss linear;
div(phi,nuTilda) Gauss limitedLinear 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

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

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

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}


and when I run the solver, simulation starts, but after some steps it crashes, giving what follows:

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::LimitedScheme<double, Foam::limitedLinearLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::fv::gaussConvectionScheme<double>::fvmDiv(Fo am::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#8 Foam::incompressible::RASModels::kEpsilon::correct () in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#9
in "/home/fsauser/OpenFOAM/fsauser-2.1.x/platforms/linux64GccDPOpt/bin/simplePipeFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
in "/home/fsauser/OpenFOAM/fsauser-2.1.x/platforms/linux64GccDPOpt/bin/simplePipeFoam"
Floating point exception


the interesting thing I noticed is that when the simulation crashes, pressure gradient values are something like 6.10715e+10 (it swings from positive to negative values step after step, with absolute values greater and greater)
----> I think it's a matter of initial conditions, but I currently I'm not yet aware of what's the point.

immortality September 12, 2013 06:45

I added phi to schemes,it should add so that solver can perform operations in proper direction,but I don't remember where it should be added,maybe only in divSchemes.

songwukong November 22, 2014 09:05

Hello, everyone, do you find an elegant solution to this problem?

I am looking forward to your reply! : )

shipman May 7, 2015 21:00

Which solver do you use? I had same problem in cavitatingFoam. First you should check your fvScheme as second your fields in 0 folder. Because according to version of openfoam it might change.


All times are GMT -4. The time now is 15:48.