CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Floating point exception (core dumped) in transonicMRFDyMFoam (https://www.cfd-online.com/Forums/openfoam-solving/134378-floating-point-exception-core-dumped-transonicmrfdymfoam.html)

sam.ho April 29, 2014 01:30

Floating point exception (core dumped) in transonicMRFDyMFoam
 
3 Attachment(s)
Dear Foamers:)
I am simulating a radial inflow turbine bassage with cyclic BC on periodic surfaces and overlapGgi at interface between Upstream and impeller and impeller and downstream.
BC and initial values are attached. Rotational speed is 2303.1833 RPS.
I am using transonicMRFDyMFoam.
When I run the case I am getting the following error.
What could be the Problem.
Please guide me.
Code:

sml5kor@BMH301562:~/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_25.04.2014$ transonicMRFDyMFoam
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM Extend Project: Open source CFD        |
|  \\    /  O peration    | Version:  1.6-ext                              |
|  \\  /    A nd          | Web:      www.extend-project.de                |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 1.6-ext
Exec  : transonicMRFDyMFoam
Date  : Apr 29 2014
Time  : 10:48:49
Host  : BMH301562
PID    : 17621
Case  : /home/sml5kor/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_25.04.2014
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh staticFvMesh
Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
Allocating field rho

Reading field U

Reading/calculating face flux field phi

Creating MRF model

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
    Prt            0.9;
    alphaK1        0.85034;
    alphaK2        1;
    alphaOmega1    0.5;
    alphaOmega2    0.85616;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1          0.075;
    beta2          0.0828;
    betaStar        0.09;
    a1              0.31;
    c1              10;
}

Create Riemann solver

Allocating field rhoU

Allocating field rhoE

Allocating physDeltaT list for RK and Dual-Time Stepping


Starting time loop


 Time = 0.001

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
rho L2 Residual: 5.19706734243579
rho LInf Residual: 6.20573608985886

DILUPBiCG:  Solving for omega, Initial residual = 0.00246886183755971, Final residual = 3.63820329446736e-09, No Iterations 1
bounding omega, min: -10676725.9096259 max: 81480676744.9903 average: 143370435.93326
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 2.09741132236342e-10, No Iterations 3
bounding k, min: 0 max: 8675.38724898184 average: 8659.50215296127

    ExecutionTime = 22.74 s

 performance data:
  Head (m)  = 5.67916442932999e+301
  TOmega (W) = 53.4715203093837
  Eff (%)    = -9.59774223073677
  Forces    = (0.0997502247768037 -0.207201815211534 1.00882792747902)
  Moments    = (-0.00637111953414955 -0.0145051677747836 -0.00369499754906964)

 Fluid power output:
  dEm (W) = -557.126030517272
  Head (m) = 5.67916442932999e+301

forces output:
    forces(pressure, viscous)((0.0224118219595117 -0.0123873558684275 0.0493579078322118) (0.0773384028172919 -0.194814459343106 0.959470019646804))
    moment(pressure, viscous)((-0.00643573254076094 0.000227855362752766 -0.000200221603628329) (6.46130066113849e-05 -0.0147330231375363 -0.00349477594544131))


 Time = 0.002

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding rho, min: -7.23432991465607 max: 87.0577875756529 average: 0.0697171344957788
bounding rhoE, min: -4199206.588189 max: 51019869.3767982 average: 30668.4900360114
bounding U max: 379103.25524316
bounding h, min: -566302.774909988 max: 361487480.588983 average: 571891.979824179
rho L2 Residual: 5.2129917782982
rho LInf Residual: 6.18558520131841

Floating point exception (core dumped)



All times are GMT -4. The time now is 17:21.