[samiam1000]

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]

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

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]

Hi,

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

Bye
Tian

[samiam1000]

Thanks a lot: Im gonna test it!

[samiam1000]

Hi Tjan,

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

Thanks a lot,
Samuele

[Tobi]

Hi,

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

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]

Originally Posted by Tobi

Hi,

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

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