CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Error with 'vanAlbada' scheme (http://www.cfd-online.com/Forums/openfoam-solving/65551-error-vanalbada-scheme.html)

 Julian K. June 18, 2009 09:31

Error with 'vanAlbada' scheme

Hallo,

I'm using the sonicFoam solver for my flow problem. For the discretization of the convective terms I'd like to use the the 'vanAlbada' scheme. Here is the excerpt from my /system/fvSchemes file:

Code:

...
divSchemes
{
default        none;
}
...

When I try to run the computation I get the following output:

Code:

...
Create time

Create mesh for time = 0

Calculating field e from T

Reading/calculating face flux field phi

Starting time loop

Time = 1e-08

Lookup interpolationScheme for interpolate(rho)
Courant Number mean: 0 max: 0
Lookup ddtScheme for ddt(rho)
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Lookup laplacianScheme for laplacian(mu,U)
Lookup fluxRequired for U
Lookup divScheme for div(phi,U)
Lookup ddtScheme for ddt(rho,U)
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 3.76603e-12, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 8.97676e-10, No Iterations 1
Lookup interpolationScheme for interpolate(rho)
Lookup laplacianScheme for laplacian(mu,e)
Lookup fluxRequired for e
Lookup divScheme for div(phi,e)
#0  Foam::error::printStack(Foam::Ostream&) in "/home/julian/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/home/julian/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  ?? in "/lib64/libc.so.6"
#3  Foam::LimitedScheme<double, Foam::vanAlbadaLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/julian/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#4  Foam::limitedSurfaceInterpolationScheme<double>::weights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/julian/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#5  Foam::fv::gaussConvectionScheme<double>::fvmDiv(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) const in "/home/julian/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/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>&, Foam::word const&) in "/home/julian/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/sonicFoam"
#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>&) in "/home/julian/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/sonicFoam"
#8  main in "/home/julian/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/sonicFoam"
#9  __libc_start_main in "/lib64/libc.so.6"
#10  Foam::regIOobject::readIfModified() in "/home/julian/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/sonicFoam"
Gleitkomma-Ausnahme

Maybe the error occured due to an incorrect definition of the scheme - I couldn't find any example and, since I'm quiet unfamiliar with OpenFOAM, I couldn't read it from the source code.

Does anyone have a suggestion?