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

Second order scheme with buoyantSimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 5, 2008, 10:52
Default Hi to all! I have a problem
  #1
Member
 
Michele Vascellari
Join Date: Mar 2009
Posts: 70
Rep Power: 8
mighelone is on a distinguished road
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:stream&) 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
mighelone is offline   Reply With Quote

Old   February 6, 2008, 06:37
Default Solving the same problem with
  #2
Member
 
Michele Vascellari
Join Date: Mar 2009
Posts: 70
Rep Power: 8
mighelone is on a distinguished road
Solving the same problem with buoyantFoam (PISO instead of SIMPLE) the second order discretization schemes works correctly.

Any ideas?

Michele
mighelone is offline   Reply With Quote

Old   September 26, 2012, 10:28
Default Mixed Convection
  #3
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 198
Rep Power: 7
deji is on a distinguished road
Good day Michele. Are you still working on mixed convection boundary layer flow? I ask such a question because it pertains to my research. I am having problems specifying boundary conditions for a vertically oriented flat plate with a incoming freestream velocity from the bottom.
deji 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
Need help with 1st order FV scheme Andrei Chernousov Main CFD Forum 1 March 14, 2013 04:37
2nd Order Scheme for SimpleFoam bastil OpenFOAM Running, Solving & CFD 0 February 22, 2009 17:20
High order FVM scheme .... cfd_newbie Main CFD Forum 2 February 6, 2008 07:46
High order FVM scheme .... cfd_newbie FLUENT 0 February 4, 2008 14:43
FVM-higher order scheme K S Chang Main CFD Forum 6 December 9, 2003 11:26


All times are GMT -4. The time now is 12:50.