Register Blogs Members List Search Today's Posts Mark Forums Read

 June 18, 2009, 09:31 Error with 'vanAlbada' scheme #1 Member   Julian Krick Join Date: May 2009 Location: Guelph Posts: 88 Rep Power: 8 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, Foam::limitFuncs::magSqr>::limiter(Foam::GeometricField const&) const in "/home/julian/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #4 Foam::limitedSurfaceInterpolationScheme::weights(Foam::GeometricField const&) const in "/home/julian/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #5 Foam::fv::gaussConvectionScheme::fvmDiv(Foam::GeometricField const&, Foam::GeometricField&) const in "/home/julian/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #6 Foam::tmp > Foam::fvm::div(Foam::GeometricField const&, Foam::GeometricField&, Foam::word const&) in "/home/julian/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/sonicFoam" #7 Foam::tmp > Foam::fvm::div(Foam::GeometricField const&, Foam::GeometricField&) 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 __________________ grid generation: ICEM CFD 13.0 solver: CFX 13.0