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

SonicTurbFoam crashes when turbulence is added

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 21, 2008, 08:49
Default Hi all, I've been doing com
  #1
New Member
 
Denis Semyonov
Join Date: Mar 2009
Location: Lappeenranta, Finland
Posts: 1
Rep Power: 0
gastovski is on a distinguished road
Hi all,

I've been doing computations of a transonic flow around an airfoil with high Re, so I've been using sonicTurbFoam solver for this case. Everything went well while I was using laminar turbulence model or no turbulence at all, but when I've added kEpsilon turbulence solver start to crash after the very first timestep iteration (or even in the very beginig in case of kOmegaSST model).

For the inlet BC I use fixed velocity value (M_infinity = 0.8 - 0.9), fixed small k and epsilon. For the outlet BC I've tried both fixed velocity outlet and zero gradient.

Could anybody help me to find out, where is the possible error?

Error messages for the k-epsilon case look like:
================================================
Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package hThermo<puremixture<consttransport<speciethermo<hc onstthermo<perfectgas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model kEpsilon
Creating field DpDt


Starting time loop

Time = 1e-07

Courant Number mean: 0.000261797 max: 0.115901
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Ux, Initial residual = 0.998782, Final residual = 0.000437629, No Iterations 1000
GAMG: Solving for Uy, Initial residual = 1, Final residual = 4.432e-09, No Iterations 2
GAMG: Solving for h, Initial residual = 0.538489, Final residual = 0.0813286, No Iterations 1000
GAMG: Solving for p, Initial residual = 1, Final residual = 1.39313e-08, No Iterations 1000
GAMG: Solving for p, Initial residual = 3.71938e-08, Final residual = 1.68472e-12, No Iterations 1
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.7637e-10, global = 8.82502e-12, cumulative = 8.82502e-12
GAMG: Solving for p, Initial residual = 6.59328e-07, Final residual = 3.56788e-11, No Iterations 1
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.11386e-09, global = -1.70428e-10, cumulative = -1.61603e-10
#0 Foam::error::printStack(Foam:stream&) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/tls/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<foam::fvpatchfield,>(Foam::GeometricF ield<double,>&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTurbul enceModels.so"
#5 Foam::tmp<foam::geometricfield<double,> > Foam::operator/<foam::fvpatchfield,>(Foam::tmp<foam::geometricfie ld<double,> > const&, Foam::GeometricField<double,> const&) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTurbul enceModels.so"
#6 Foam::compressible::turbulenceModels::kEpsilon::co rrect() in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTurbul enceModels.so"
#7 main in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/sonicTur bFoam"
#8 __libc_start_main in "/lib64/tls/libc.so.6"
#9 __gxx_personality_v0 in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/sonicTur bFoam"
Floating exception
195.862u 6.226s 3:42.68 90.7% 0+0k 0+0io 2515pf+0w
=================================================
gastovski 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
Possible sonicTurbFoam bug guusama OpenFOAM Bugs 1 December 17, 2007 13:58
SonicTurbFoam Kaushik Balakrishnan (Balakrishnan) OpenFOAM Running, Solving & CFD 5 May 28, 2007 10:20
added mass Kuba Broniszewski FLUENT 2 February 23, 2007 10:39
SonicTurbFoam anja OpenFOAM Pre-Processing 1 May 22, 2006 06:50
Effect of added CPU ? steve podleski Main CFD Forum 8 August 10, 1999 20:26


All times are GMT -4. The time now is 16:07.