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

buoyantPimpleFoam: test case

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

Like Tree2Likes
  • 2 Post By tian

Reply
 
LinkBack Thread Tools Display Modes
Old   June 5, 2013, 10:35
Default buoyantPimpleFoam: test case
  #1
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
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
Attached Files
File Type: gz test.tar.gz (2.8 KB, 6 views)
samiam1000 is offline   Reply With Quote

Old   June 6, 2013, 06:49
Default
  #2
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
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?
samiam1000 is offline   Reply With Quote

Old   June 7, 2013, 02:08
Default Other Test Case
  #3
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 103
Rep Power: 8
tian is on a distinguished road
Hi,

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

Bye
Tian
Attached Files
File Type: zip hotRoom-HVACTool.zip (12.3 KB, 12 views)
samiam1000 and Tobi like this.
tian is offline   Reply With Quote

Old   June 9, 2013, 12:28
Default
  #4
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
Thanks a lot: Im gonna test it!
samiam1000 is offline   Reply With Quote

Old   June 11, 2013, 03:49
Default
  #5
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
Hi Tjan,

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

Thanks a lot,
Samuele
samiam1000 is offline   Reply With Quote

Old   June 11, 2013, 06:22
Default
  #6
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,100
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
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
Tobi is offline   Reply With Quote

Old   January 8, 2014, 13:04
Default
  #7
Member
 
Jamal
Join Date: May 2012
Location: Freiburg
Posts: 54
Rep Power: 4
aujamal20 is an unknown quantity at this point
Quote:
Originally Posted by Tobi View Post
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
aujamal20 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
rhoSimpleFoam - test case atareen64 OpenFOAM Running, Solving & CFD 5 May 1, 2011 17:17
Pressure instabilities with interDyMFoam for the floatingObject case nbadano OpenFOAM Running, Solving & CFD 14 March 8, 2011 09:19
Help Required With Simple Test Case steph79 OpenFOAM Pre-Processing 4 August 3, 2010 07:45
3D TRANSITION TEST CASE venkatesh4386@gmail.com FLUENT 0 March 9, 2009 13:04
Porous Media test case Alex FLUENT 0 April 9, 2006 08:23


All times are GMT -4. The time now is 01:29.