compressible solver rhosimplefoam "sigFe" error

ruby_nuaa · July 17, 2015, 22:07

Hi everyone

I'm trying to simulate a flow field around a helicopter blade in hover using rhoSimpleFoam/rhoSimplecFoam. The version is OF-2.2.2 . But I keep getting this error in 2 interations.

Code:

Create time 

Create mesh for time = 0 


SIMPLE: convergence criteria 
    field p     tolerance 1e-05 
    field U     tolerance 1e-06 
    field h     tolerance 0.001 
    field "(k|omega)"     tolerance 1e-05 

Reading thermophysical properties 

Selecting thermodynamics package 
{ 
    type            hePsiThermo; 
    mixture         pureMixture; 
    transport       sutherland; 
    thermo          hConst; 
    equationOfState perfectGas; 
    specie          specie; 
    energy          sensibleEnthalpy; 
} 

Reading field U 

Reading/calculating face flux field phi 

Creating turbulence model 

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

Creating fintite volume options from fvOptions 

Selecting finite volume options model type MRFSource 
    Source: MRF1 
    - applying source for all time 
    - selecting cells using cellZone MRF-cells 
    - selected 17644240 cell(s) with volume 21.2025 


Starting time loop 

forces forces: 

Time = 0.01 

DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.009296, No Iterations 2 
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.012759, No Iterations 2 
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.0126077, No Iterations 2 
DILUPBiCG:  Solving for h, Initial residual = 0.510074, Final residual = 0.00850355, No Iterations 1 
DICPCG:  Solving for p, Initial residual = 1, Final residual = 0.00993348, No Iterations 170 
DICPCG:  Solving for p, Initial residual = 6.77265e-05, Final residual = 9.45422e-07, No Iterations 72 
DICPCG:  Solving for p, Initial residual = 3.33876e-05, Final residual = 9.58751e-07, No Iterations 35 
time step continuity errors : sum local = 4.33681e-06, global = -6.56944e-08, cumulative = -6.56944e-08 
rho max/min : 2 0.5 
DILUPBiCG:  Solving for omega, Initial residual = 0.0144652, Final residual = 0.000953703, No Iterations 1 
DILUPBiCG:  Solving for k, Initial residual = 0.999988, Final residual = 0.0653421, No Iterations 1 
ExecutionTime = 485.77 s  ClockTime = 486 s 

Time = 0.02 

DILUPBiCG:  Solving for Ux, Initial residual = 0.714487, Final residual = 0.00670596, No Iterations 2 
DILUPBiCG:  Solving for Uy, Initial residual = 0.14934, Final residual = 0.00132467, No Iterations 2 
DILUPBiCG:  Solving for Uz, Initial residual = 0.715737, Final residual = 0.071483, No Iterations 1 
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 0.07682, No Iterations 1 
#0  Foam::error::printStack(Foam::Ostream&) at ??:? 
#1  Foam::sigFpe::sigHandler(int) at ??:? 
#2   in "/lib64/libc.so.6" 
#3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:? 
#4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:? 
#5  
 at ??:? 
#6  __libc_start_main in "/lib64/libc.so.6" 
#7  
 at ??:? 
#

Now i will describe my case. The compute domain is cylindrical. There are 4 boundaries defined as “top bottom blade far ”. “Top” stands for the top of the cylinder which meant to be an inlet. “Far & bottom” stand for the sides and bottom of the cylinder as outlet. The “blades” means wall.

Here are my p and U files
p

Code:

dimensions      [1 -1 -2 0 0 0 0]; 

internalField   uniform 1.0e5; 

boundaryField 
{ 
    "blade.*" 
    { 
        type            zeroGradient; 
    } 
  
   
   "top.*" 
    { 
        /*type            totalPressure; 
         p0             uniform 102152; 
         U              U; 
         phi            phi; 
         rho            none; 
         psi            none; 
         gamma           1.4; 
         value          uniform 102152;*/ 
           type          zeroGradient; 
    } 
   "far|bottom.*" 
    { 
        type            fixedValue; 
         value          $internalField; 
         /* type            totalPressure; 
         p0             uniform 98564; 
         U              U; 
         phi            phi; 
         rho            none; 
         psi            none; 
         gamma           1.4; 
         value          uniform 98564;*/ 
    } 


}

U

Code:

dimensions      [0 1 -1 0 0 0 0]; 

internalField   uniform (0.0 0 0); 

boundaryField 
{ 
    "blade.*" 
    { 
        type            fixedValue;  
        value           uniform (0 0 0); 
    } 
   
    
   "top.*" 
    { 
        type            fixedValue; 
        value           uniform (0 0 0); 
        /*type            pressureInletOutletVelocity; 
        value           uniform (0 0 0);*/ 
    } 
    "far|bottom.*" 
    { 
        type            inletOutlet; 
        inletValue      uniform (0 0 0); 
        value           uniform (0 0 0); 
        /*type            pressureInletOutletVelocity; 
        value           uniform (0 0 0); */       
    } 

}

and Here are the fvSchemes and fvSolution

Code:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 

ddtSchemes 
{ 
    default         Euler; 
} 

gradSchemes 
{ 
    default         Gauss linear; 
    grad(p)         faceLimited leastSquares 1; //cellLimited Gauss linear 1; //Gauss linear; 
    grad(U)         faceLimited leastSquares 1; 
} 

divSchemes 
{ 
    default         none; 
    div(phi,U)      Gauss upwind; 
    //div(phi,U)      Gauss linearUpwind grad(U); 
    div(phi,k)      Gauss upwind; 
    div(phi,omega) Gauss upwind; 
    div((nuEff*dev(T(grad(U)))))     Gauss upwind; //linear; 
    div((muEff*dev2(T(grad(U)))))     Gauss linear; 
    div(phi,Ekp)      Gauss upwind; 
    div(phi,h)      Gauss upwind;  //limitedLinear 1; 
    div(phi,e)      Gauss upwind;  //limitedLinear 1; 
    div(phi,K)      Gauss linear; 

     /*div(phi,U)      Gauss limitedLinearV 1; 
    div(phi,h)      Gauss limitedLinear 1; 
    div(phi,k)      Gauss limitedLinear 1; 
    div(phi,omega) Gauss limitedLinear 1; 
    div((muEff*dev2(T(grad(U))))) Gauss linear; 
    div(phi,K)      Gauss linear;*/ 
} 

laplacianSchemes 
{ 
    default         Gauss linear limited 0.33; 
} 

interpolationSchemes 
{ 
    default         linear; 
} 

snGradSchemes 
{ 
    default         corrected;         //limited  0.33; 
} 

fluxRequired 
{ 
    default         no; 
    pcorr           ; 
    p               ; 
}

fvSolution

Code:

solvers 
{ 
       rho 
    { 
        solver          PCG; 
        preconditioner  DIC; 
        tolerance       1e-5; 
        relTol          0.1; 
    } 

    rhoFinal 
    { 
        $rho; 
        tolerance       1e-5; 
        relTol          0; 
    } 

   p 
    { 
        solver          PCG; 
        preconditioner  DIC; 
        tolerance       1e-06; 
        relTol          0.01; 
    } 

    U 
    { 
        solver          PBiCG; 
        preconditioner  DILU; 
        tolerance       1e-05; 
        relTol          0.1; 
    } 

    k 
    { 
        solver          PBiCG; 
        preconditioner  DILU; 
        tolerance       1e-05; 
        relTol          0.1; 
    } 

    omega 
    { 
        solver          PBiCG; 
        preconditioner  DILU; 
        tolerance       1e-05; 
        relTol          0.1; 
    } 

    "(e|h|R)" 
    { 
        solver          PBiCG; 
        preconditioner  DILU; 
        tolerance       1e-05; 
        relTol          0.1; 
    } 

    nuTilda 
    { 
        solver          PBiCG; 
        preconditioner  DILU; 
        tolerance       1e-05; 
        relTol          0.1; 
    } 
} 

SIMPLE 
{ 
    nNonOrthogonalCorrectors 2; 
    rhoMin          rhoMin [ 1 -3 0 0 0 ] 0.5; 
    rhoMax          rhoMax [ 1 -3 0 0 0 ] 2.0; 
    residualControl 
    { 
        p               1e-5; 
        U               1e-6; 
        h               1e-3; 
        "(k|omega)" 1e-5; 
    } 
} 

relaxationFactors 
{ 
    fields 
    { 
        /*p               0.2; 
        rho             0.05;*/ 
    } 
    equations 
    { 
        U               0.7; 
        k               0.7; 
        omega           0.7; 
        R               0.7; 
        nuTilda         0.7; 
        h               0.7; 
    } 
}

I know that the “float exception core dumped “ means there's divided by zero like condition here. But i checked the thermo variables they seemd ok. I've tried with the totalPressure bc with pressure ,and low the relaxation factors of p and rho to 0.01, but it never worked out. The mesh check is ok btw. So something wrong with my bcs? I must have missed somthing.

Can anyone help me please? I've been struggled for days.

ruby_nuaa · July 17, 2015, 22:10

sorry about the code format mess .
something wrong with my ie , i can't get a response on the quote or # botton.
sorry again.

and please help.

wyldckat · July 18, 2015, 18:58

Quick question: Can you please create and share a small test example, based on your case? Otherwise this will be a very big guessing game

BSengupta · January 7, 2019, 05:24

Hi,

I am having few problems. Although residuals are quite low and no out of bounds still:

Code:

Courant Number mean: 3.35064e-13 max: 0.301157
deltaT = 1.90835e-11
Time = 0.546089

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 3.56934e-12, Final residual = 3.56934e-12, No Iterations 0
smoothSolver:  Solving for Uy, Initial residual = 3.52604e-12, Final residual = 3.52604e-12, No Iterations 0
smoothSolver:  Solving for Uz, Initial residual = 4.2873e-12, Final residual = 4.2873e-12, No Iterations 0
smoothSolver:  Solving for h, Initial residual = 9.95909e-06, Final residual = 9.95909e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 5.94342e-12, Final residual = 5.94342e-12, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.25413e-14, global = -4.70952e-15, cumulative = -4.0427e-07
GAMG:  Solving for p, Initial residual = 5.94342e-12, Final residual = 5.94342e-12, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.25414e-14, global = -4.70952e-15, cumulative = -4.0427e-07
smoothSolver:  Solving for epsilon, Initial residual = 2.97988e-09, Final residual = 2.97988e-09, No Iterations 0
smoothSolver:  Solving for k, Initial residual = 5.14549e-07, Final residual = 5.14549e-07, No Iterations 0
ExecutionTime = 22051.5 s  ClockTime = 22078 s

Courant Number mean: 3.33784e-13 max: 0.301157
deltaT = 1.90102e-11
Time = 0.546089

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 3.56857e-12, Final residual = 3.56857e-12, No Iterations 0
smoothSolver:  Solving for Uy, Initial residual = 3.52543e-12, Final residual = 3.52543e-12, No Iterations 0
smoothSolver:  Solving for Uz, Initial residual = 4.2808e-12, Final residual = 4.2808e-12, No Iterations 0
smoothSolver:  Solving for h, Initial residual = 1.0036e-05, Final residual = 5.10147e-08, No Iterations 1
[1] #0  Foam::error::printStack(Foam::Ostream&)[2] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
 at ??:?
[1] #1  [2] #1  Foam::sigFpe::sigHandler(int)Foam::sigFpe::sigHandler(int) at ??:?
[2] #2  ? at ??:?
[1] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3  Foam::hePsiThermo<Foam::psiThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::hePsiThermo<Foam::psiThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
[2] #4  Foam::hePsiThermo<Foam::psiThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
[1] #4  Foam::hePsiThermo<Foam::psiThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
[2] #5  ? at ??:?
[1] #5  ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
[2] #6  __libc_start_main in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
[1] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #7  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #7  ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
[bidesh:03000] *** Process received signal ***
[bidesh:03000] Signal: Floating point exception (8)
[bidesh:03000] Signal code:  (-6)
[bidesh:03000] Failing at address: 0x3e800000bb8
[bidesh:03000] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f1b958d3f20]
[bidesh:03000] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f1b958d3e97]
[bidesh:03000] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f1b958d3f20]
[bidesh:03000] [ 3] /opt/openfoam6/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE9calculateEv+0x2b6)[0x7f1b9c3b4ac6]
[bidesh:03000] [ 4] /opt/openfoam6/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE7correctEv+0x2a)[0x7f1b9c3cfc1a]
[bidesh:03000] [ 5] /opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam(+0x2ef91)[0x55e7b6cb1f91]
[bidesh:03000] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f1b958b6b97]
[bidesh:03000] [ 7] /opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam(+0x36eca)[0x55e7b6cb9eca]
[bidesh:03000] *** End of error message ***
 in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
[bidesh:03004] *** Process received signal ***
[bidesh:03004] Signal: Floating point exception (8)
[bidesh:03004] Signal code:  (-6)
[bidesh:03004] Failing at address: 0x3e800000bbc
[bidesh:03004] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f1cce2cbf20]
[bidesh:03004] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f1cce2cbe97]
[bidesh:03004] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f1cce2cbf20]
[bidesh:03004] [ 3] /opt/openfoam6/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE9calculateEv+0x2b6)[0x7f1cd4dacac6]
[bidesh:03004] [ 4] /opt/openfoam6/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE7correctEv+0x2a)[0x7f1cd4dc7c1a]
[bidesh:03004] [ 5] /opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam(+0x2ef91)[0x56316ef7af91]
[bidesh:03004] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f1cce2aeb97]
[bidesh:03004] [ 7] /opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam(+0x36eca)[0x56316ef82eca]
[bidesh:03004] *** End of error message ***
^Cmpirun: abort is already in progress...hit ctrl-c again to forcibly terminate

Although I tried to solve it by studying from CFDonline but still problem persists. It is basically 3D O-Grid and at the center the object is placed. So, I have 2 BC namely: farfield and wall.

I have initialized k and epsilon according to:
k/(U*U)=1e-6, epsilon*charecteristic length of object/(U*U*U)=4.5e-7

[Ref: https://www.simscale.com/forum/t/ini...cients/25950/7

I shall be grateful if anyone can help me out. Thank you.

wyldckat · January 8, 2019, 18:12

Quick answer: Partially quoting myself from a recent post of mine:

Quote:

Originally Posted by wyldckat

has "sigFpe" in there, which means that there was a floating point error. This means that either a division by zero or some other illegal mathematical operation was attempted.

The illegal mathematical operation occurred at:

in your case, here:

Code:

[2] #3  Foam::hePsiThermo<Foam::psiThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::hePsiThermo<Foam::psiThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?

Quote:

Originally Posted by wyldckat

So in other words, it occurred at a thermodynamic calculation. Usually this is associated with an absolute pressure value of 0 Pa or a absolute temperature of 0 K.

This likely occurred when you switched from pimpleFoam to rhoPimpleFoam, given that rhoPimpleFoam works with absolute pressures and not relative pressures as in pimpleFoam.

Furthermore, the very small time steps and the very small residual values lead me to suspect that something very wrong is happening in your case, so I really suspect that you are using a pressure value of zero somewhere. It would be as if you suddenly placed your mesh in outer space and the air around escaped at a tremendous speed.

BSengupta · January 8, 2019, 20:36

Thank you very much for your suggestion. I shall look into it the way you said.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Compiling dynamicTopoFvMesh for OpenFOAM 2.1.x	Saxwax	OpenFOAM Installation	25	November 29, 2013 06:34
[swak4Foam] installing funkySetFields	igo	OpenFOAM Community Contributions	1	November 20, 2012 21:16
Compile problem	ivanyao	OpenFOAM Running, Solving & CFD	1	October 12, 2012 10:31
ParaView for OF-1.6-ext	Chrisi1984	OpenFOAM Installation	0	December 31, 2010 07:42
user subroutine error	CFDUSER	CFX	2	December 9, 2006 07:31

July 17, 2015, 22:10		#2
ruby_nuaa New Member Ruby Qian Join Date: Aug 2013 Location: Nanjing,Jiangsu,PRC Posts: 13 Rep Power: 13	sorry about the code format mess . something wrong with my ie , i can't get a response on the quote or # botton. sorry again. and please help.

July 18, 2015, 18:58		#3
wyldckat Retired Super Moderator Bruno Santos Join Date: Mar 2009 Location: Lisbon, Portugal Posts: 10,981 Blog Entries: 45 Rep Power: 128	Quick question: Can you please create and share a small test example, based on your case? Otherwise this will be a very big guessing game

January 8, 2019, 20:36		#6
BSengupta Member Bidesh Sengupta Join Date: Sep 2018 Location: Sngapore Posts: 76 Rep Power: 8	Thank you very much for your suggestion. I shall look into it the way you said.