CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   buoyantPimpleFoam: test case (http://www.cfd-online.com/Forums/openfoam/118845-buoyantpimplefoam-test-case.html)

samiam1000 June 5, 2013 10:35

buoyantPimpleFoam: test case
 
1 Attachment(s)
Dear all,

I am trying to set up a case with the solver buoyantPimpleFoam, but I can't make the simulation start.

Attached a test case. If you launch it (execute the blockMesh, first), you get an error. Immediately.

Could you have a look?

Thanks a lot,
Samuele

samiam1000 June 6, 2013 06:49

This is the error that you get, when launching the solver:
Code:

lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/NOPCM/testTransient$ buoyantPimpleFoam
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.1.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 2.1.0-0bc225064152
Exec  : buoyantPimpleFoam
Date  : Jun 06 2013
Time  : 12:45:38
Host  : "lab-laptop"
PID    : 4188
Case  : /home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/NOPCM/testTransient
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Selecting thermodynamics package hRhoThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model realizableKE
realizableKECoeffs
{
    Cmu            0.09;
    A0              4;
    C2              1.9;
    sigmak          1;
    sigmaEps        1.2;
    Prt            1;
}

Calculating field g.h

Reading field p_rgh

Creating field dpdt

Creating field kinetic energy K

Courant Number mean: 0.0340283 max: 0.371779

PIMPLE: Operating solver in PISO mode


Starting time loop

Courant Number mean: 0.0340283 max: 0.371779
Time = 0.0001

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  in "/lib/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  Foam::operator/(Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  Foam::diagonalSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7  Foam::fvMatrix<double>::solve() in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/buoyantPimpleFoam"
#8 
 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/buoyantPimpleFoam"
#9  __libc_start_main in "/lib/libc.so.6"
#10 
 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/buoyantPimpleFoam"
Floating point exception
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/NOPCM/testTransient$

It seems that there is a division by zero. But where is it? And why this happens? Because of the wrong BC?

tian June 7, 2013 02:08

Other Test Case
 
1 Attachment(s)
Hi,

here is an other test case for you. Maybe you compare with your case...

Bye
Tian

samiam1000 June 9, 2013 12:28

Thanks a lot: Im gonna test it!

samiam1000 June 11, 2013 03:49

Hi Tjan,

I can't open the archive. It seems corrupted. Could you check it, please?

Thanks a lot,
Samuele

Tobi June 11, 2013 06:22

Hi,

it s not possible for me to test your case but set a value in the file p:

Code:

inlet
{
type buoyantPressure;
value $internalField;
}

And set a value for every "zeroGradient" BC (just to be safe)

That is the only hint at the moment

aujamal20 January 8, 2014 12:04

Quote:

Originally Posted by Tobi (Post 433327)
Hi,

it s not possible for me to test your case but set a value in the file p:

Code:

inlet
{
type buoyantPressure;
value $internalField;
}

And set a value for every "zeroGradient" BC (just to be safe)

That is the only hint at the moment

Hi,
I am working on buoyantPimpleFoam case. I prepared a wedge shaped geometry and using almost similar setting as of hotRoom tutorial. As the simulation is started the deltaT starts decreasing and keep continue to decrease and reaches upto the power of 1e-5/6 and then it stop with an error message of "iterations exceeded".
Is there any possible way to run a stable simulation with buoyantPimpleFoam. working on OF 2.2.x

Regards,

Jamal


All times are GMT -4. The time now is 20:57.