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;     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?