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;
div(phi,U) Gauss vanAlbada;
div(phid,p) Gauss vanAlbada;
div(phi,e) Gauss vanAlbada;
}
...
When I try to run the computation I get the following output:
Code:
...
Create time
Create mesh for time = 0
Reading thermodynamicProperties
Reading transportProperties
Reading field p
Reading field T
Calculating field e from T
Reading field U
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 gradScheme for snGradCorr(U)
Lookup gradScheme for snGradCorr(U)
Lookup gradScheme for snGradCorr(U)
Lookup divScheme for div(phi,U)
Lookup gradScheme for grad(magSqr(U))
Lookup ddtScheme for ddt(rho,U)
Lookup gradScheme for grad(p)
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 gradScheme for grad(U)
Lookup interpolationScheme for interpolate(rho)
Lookup laplacianScheme for laplacian(mu,e)
Lookup fluxRequired for e
Lookup gradScheme for snGradCorr(e)
Lookup divScheme for div(phi,e)
Lookup gradScheme for grad(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?
Thanks for reading my thread,
bye
Julian
|