CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Error with 'vanAlbada' scheme

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 18, 2009, 09:31
Default Error with 'vanAlbada' scheme
  #1
Member
 
Julian Krick
Join Date: May 2009
Location: Guelph
Posts: 88
Rep Power: 8
Julian K. is on a distinguished road
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
__________________
grid generation: ICEM CFD 13.0
solver: CFX 13.0
Julian K. is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Implementation of QUICK scheme Romuald Skoda Main CFD Forum 10 August 3, 2010 12:23
AUSM scheme ? Central Scheme boling Main CFD Forum 5 September 9, 2009 19:06
Definition of limiter function for central dirrerencing scheme sebastian_vogl OpenFOAM Running, Solving & CFD 0 January 5, 2009 12:08
Best Meshing scheme for Cylinder Nutrex Main CFD Forum 4 July 29, 2008 11:03
extrapolation in MUSCL scheme Chandra Main CFD Forum 6 February 14, 2007 12:21


All times are GMT -4. The time now is 09:15.