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

Floating Point Exception Error for p_rgh in chtMRF

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 8, 2021, 14:39
Default Floating Point Exception Error for p_rgh in chtMRF
  #1
Member
 
Claudia
Join Date: Mar 2021
Posts: 43
Rep Power: 3
Clau.77 is on a distinguished road
Hey guys!


i have a multi region case in OF8 and i am having trouble with my fluid.

To keep it simple i am using steadyState and laminar flow.
The solver log looks like this:
Code:
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 8
Exec   : chtMultiRegionFoam
Date   : May 08 2021
Time   : 19:20:03
Host   : "login02.css.lan"
PID    : 51958
I/O    : uncollated
Case   : /bigwork/nhkfclan/Test_Merci
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region Luft for time = 0

Create solid mesh for region Form for time = 0

Create solid mesh for region Schokolade for time = 0

*** Reading fluid mesh thermophysical properties for region Luft

    Adding to thermoFluid

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

    Adding to rhoFluid

    Adding to UFluid

    Adding to phiFluid

    Adding to gFluid

    Adding to hRefFluid

    Adding to pRefFluid

    Adding to ghFluid

    Adding to ghfFluid

    Adding to turbulenceFluid

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
    Adding to thermophysicalTransport

Selecting thermophysical transport type laminar
Selecting default laminar thermophysical transport model Fourier
    Adding to reactionFluid

Combustion model not active: combustionProperties not found
Selecting combustion model none
    Adding to radiationFluid

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding to KFluid

    Adding to dpdtFluid

    Adding to fieldsFluid

    Adding MRF

No MRF models present

    Adding fvOptions

*** Reading solid mesh thermophysical properties for region Form

    Adding to thermos

Selecting thermodynamics package 
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          eConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

*** Reading solid mesh thermophysical properties for region Schokolade

    Adding to thermos

Selecting thermodynamics package 
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          eConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions


PIMPLE: Region Luft
PIMPLE: No convergence criteria found


PIMPLE: Region Form
PIMPLE: No convergence criteria found


PIMPLE: Region Schokolade
PIMPLE: No convergence criteria found


PIMPLE: Operating solver in steady-state mode with 1 outer corrector
PIMPLE: Operating solver in SIMPLE mode


Region: Luft Courant Number mean: 314.33656 max: 7014.9304
Region: Form Diffusion Number mean: 0.10103954 max: 261.63834
Region: Schokolade Diffusion Number mean: 101.25719 max: 251280.78
deltaT = 3.979466e-05
Region: Luft Courant Number mean: 0.012508916 max: 0.27915677
Region: Form Diffusion Number mean: 4.0208341e-06 max: 0.010411809
Region: Schokolade Diffusion Number mean: 0.0040294954 max: 9.9996333
deltaT = 3.979466e-05
Time = 3.97947e-05


Solving for fluid region Luft
DILUPBiCGStab:  Solving for Ux, Initial residual = 1, Final residual = 4.1108316e-09, No Iterations 2
DILUPBiCGStab:  Solving for Uy, Initial residual = 1, Final residual = 2.4053866e-08, No Iterations 2
DILUPBiCGStab:  Solving for Uz, Initial residual = 1, Final residual = 7.1445854e-08, No Iterations 2
DILUPBiCGStab:  Solving for h, Initial residual = 1, Final residual = 5.3058713e-08, No Iterations 3
Min/max T:284.67219 291.08218
I know that after that step something with p_rgh should be next, but then i am just getting a long error message which says floating point exception at the end. So I guess there is a problem with p_rgh or could it be something else? I know that it means, that something is devided by zero.



In the 0 directory for the fluid (air) i use fixedFluxPressure (and fixedValue at the outlet) for p_rgh and calculate for p. Both internalFields are 1e5.



the fluid is standard air, here the thermophysical properties:
Code:
thermoType
{
    type              heRhoThermo;  
    mixture         pureMixture;
    transport       const; 
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}


mixture 
{
    specie
    {
        //nMoles        1;
        molWeight   28.9;
    }
    thermodynamics
    {
        Cp          1000;
        Hf          0;
    }
    transport
    {
        mu          1.8e-05;
        Pr          0.7;   
    }
}
fvSchemes:
Code:
ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;

    div(phi,U)      Gauss upwind;
    div(phi,K)      Gauss linear;
    div(phi,h)      Gauss upwind;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon)      Gauss upwind;
    div(phi,omega)  Gauss upwind;    
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}
wallDist
{
method        meshWave;
}
Code:
solvers
{

    rho
    {
        solver          diagonal;
    }

    rhoFinal
    {
        $rho;
    }

    p_rgh
    {
        solver           GAMG;
        smoother         symGaussSeidel;
        tolerance        1e-7;
        relTol           0.01;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol           0;
    }

    
    p
    {
        solver           GAMG;
        smoother         symGaussSeidel;
        tolerance        1e-3;
        relTol           0.01;
    }

    p_Final
    {
        $p;
        relTol           0;
    }


    "(U|h|epsilon|omega)"
    {
        solver           PBiCGStab;
        preconditioner   DILU;
        tolerance        1e-7;
        relTol           0.1;
    }

    "(U|epsilon|omega|h)Final"
    {
        $U;
        relTol           0;
    }
}

PIMPLE
{
    momentumPredictor   yes;
}

relaxationFactors
{
    fields
    {
    p_rgh          0.7;
    p_rghFinal  0.3;
    rho         1;
    rhoFinal    0.1;
    }       

 equations
    {
        h              0.3;
    hFinal           0.3;
        U              0.1;
    UFinal           0.1;
    k            0.1;
    //omega        0.1;
    //epsilon        0.1;
    }
}
I played around with the fvSolution and fvSchemes but the error doesn't go away.


Does somebody have an idea where this error might come from and how to fix it? I would appreciate some help, I have been stuck with this error for days.



Also, let me know if I forgot to give important information.


Thanks!
Clau.77 is offline   Reply With Quote

Old   May 9, 2021, 10:41
Default
  #2
Member
 
Claudia
Join Date: Mar 2021
Posts: 43
Rep Power: 3
Clau.77 is on a distinguished road
Update:
it's a problem with the steadyState, different ddtSchemes work. But i still don't understand why
Clau.77 is offline   Reply With Quote

Reply

Tags
chtmultiregion, floatingpointexception, multiregion, openfoam 8, pressure

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
icoFoam floating point exception (8) leizhao512 OpenFOAM Running, Solving & CFD 7 November 1, 2018 12:43
A floating point exception has occurred: floating point exception [Overflow]. starlight STAR-CCM+ 4 May 4, 2016 10:08
A floating point exception - SEM Model yansheng STAR-CCM+ 1 April 4, 2016 05:57
Floating point exception from twoPhaseEulerFoam openfoammaofnepo OpenFOAM Running, Solving & CFD 1 March 19, 2016 14:56
Floating point exception (core dumped) for GAMG solver yuhou1989 OpenFOAM Running, Solving & CFD 2 March 24, 2015 20:28


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