CFD Online Discussion Forums - Second order scheme with buoyantSimpleFoam

Hi to all!

I have a problem running a simple problem of mixed natural and forced convection in a square cavity, where the flow enters at the top of the cavity, and exits at the bottom of the opposite side (Blay cavity).

The simulation runs without problem using upwind scheme for convective therms, but If I turn to second order scheme (limitedLinear) I obtain the following error:

Create mesh for time = 2000

Reading environmentalProperties
Reading thermophysical properties

Selecting thermodynamics package hThermo<puremixture<consttransport<speciethermo<hc onstthermo<perfectgas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model kEpsilon
Calculating field g.h

Creating field pd

Starting time loop

Time = 2001

DILUPBiCG: Solving for Ux, Initial residual = 0.0401578, Final residual = 2.36105e-06, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.891016, Final residual = 8.7564e-07, No Iterations 6
DILUPBiCG: Solving for h, Initial residual = 0.00604865, Final residual = 7.18185e-06, No Iterations 3
DICPCG: Solving for pd, Initial residual = 0.911804, Final residual = 7.97762e-09, No Iterations 86
time step continuity errors : sum local = 7.71829e-07, global = 1.66559e-09, cumulative = 1.66559e-09
rho max/min : 1.22225 1.14107
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::Istream&) in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<foam::dilupreco nditioner>::New(Foam::lduMatrix::solver const&, Foam::Istream&) in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::Istream&) in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#8 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#9 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<foam::fvmatrix<doubl e> > const&) in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#10 Foam::compressible::turbulenceModels::kEpsilon::co rrect() in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTurbul enceModels.so"
#11 main in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/buoyantS impleFoam"
#12 __libc_start_main in "/lib/libc.so.6"
#13 Foam::regIOobject::readIfModified() in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/buoyantS impleFoam"
Errore di virgola mobile

Reducing under-relaxation coefficients the error still remain the same.

The scheme used for the convective terms are the following:

divSchemes
{
default none;
div(phi,U) Gauss limitedLinear 1;
div(phi,h) Gauss limitedLinear 1;
div(phi,k) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div(phi,R) Gauss linear;
/* div(phi,U) Gauss upwind;
div(phi,h) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;*/
div(R) Gauss linear;
div((muEff*dev2(grad(U).T()))) Gauss linear;
}

Many thanks

Michele