CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

 Julian K. June 18, 2009 09:31

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?