Second order scheme with buoyantSimpleFoam

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

 February 5, 2008, 10:52 Hi to all! I have a problem #1 Member   Michele Vascellari Join Date: Mar 2009 Posts: 70 Rep Power: 10 Sponsored Links 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>>>> 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&, 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::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&, Foam::Field const&, unsigned char) const in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so" #8 Foam::fvMatrix::solve(Foam::Istream&) in "/home/michele/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so" #9 Foam::lduMatrix::solverPerformance Foam::solve(Foam::tmp > 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

 February 6, 2008, 06:37 Solving the same problem with #2 Member   Michele Vascellari Join Date: Mar 2009 Posts: 70 Rep Power: 10 Solving the same problem with buoyantFoam (PISO instead of SIMPLE) the second order discretization schemes works correctly. Any ideas? Michele

 September 26, 2012, 10:28 Mixed Convection #3 Senior Member   n/a Join Date: Sep 2009 Posts: 198 Rep Power: 9 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Andrei Chernousov Main CFD Forum 1 March 14, 2013 04:37 bastil OpenFOAM Running, Solving & CFD 0 February 22, 2009 17:20 cfd_newbie Main CFD Forum 2 February 6, 2008 07:46 cfd_newbie FLUENT 0 February 4, 2008 14:43 K S Chang Main CFD Forum 6 December 9, 2003 11:26