- **OpenFOAM Running, Solving & CFD**
(*http://www.cfd-online.com/Forums/openfoam-solving/*)

- - **SonicTurbFoam crashes when turbulence is added**
(*http://www.cfd-online.com/Forums/openfoam-solving/59213-sonicturbfoam-crashes-when-turbulence-added.html*)

Hi all,
I've been doing comHi 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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) 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 ================================================= |

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