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

SonicFoam forwardStepTutorial doesnbt complete the run

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 10, 2007, 14:44
Default Description: The forwardStep t
  #1
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
Description: The forwardStep tutorial of sonicFoam doesn't complete the simulation and aborts with the error:

Courant Number mean: -8.55342e+25 max: 1.66084e+28
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
#0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xffffe420]
#3 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::Istream&)
#4 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<foam::dilupreco nditioner>::New(Foam::lduMatrix::solver const&, Foam::Istream&)
#5 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::Istream&)
#6 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const
#7 Foam::fvMatrix<foam::vector<double> >::solve(Foam::Istream&)
#8 Foam::lduMatrix::solverPerformance Foam::solve<foam::vector<double> >(Foam::tmp<foam::fvmatrix<foam::vector<double> > > const&)
#9 main
#10 __libc_start_main
#11 __gxx_personality_v0 at ../sysdeps/i386/elf/start.S:122

Increasing the time step to 0.003s circumvents the problem.

Solver/Application: sonicFoam

Source file:

Testcase: forwardStep tutorial for sonicFoam

Platform: All

Version: 1.4

Notes: The issue was found by Florian Fruth. See http://www.cfd-online.com/cgi-bin/Op...how.cgi?1/4657

Regards,
Alberto
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541)
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto is offline   Reply With Quote

Old   June 10, 2007, 15:35
Default As suggested by Hrvoje, the so
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
As suggested by Hrvoje, the solution consists in setting the outlet BC to zeroGradient for T and U.

I attach the corrected tutorial case.

forwardStep.zip

Regards,
A.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541)
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
FOAM The Complete Guide nico765 OpenFOAM Running, Solving & CFD 4 December 6, 2013 10:26
mesh a complete plane horsee FLUENT 13 April 17, 2010 20:20
how to view the complete results shokry FLOW-3D 5 February 3, 2009 13:56
Complete mixing Jamie Main CFD Forum 0 October 6, 2008 10:47
Simple Q. How to complete a Case? W.A. Siemens 3 April 23, 2006 19:41


All times are GMT -4. The time now is 14:06.