[Sakun]

Hello everyone,


I am trying to simulate my high-speed single blade simulation in rhoCentralFoam. It is running for couple of iterations and suddenly crash. Error that generate doesn't give much information either.
I need to run this case using kOmega SST and i have attached fvSchemes and fvSolution files below.


Please someone help me on this issue?


Error that generated;

deltaT = 4.11195e-09
Mean and max Courant Numbers = 0.00342632 1.6291
Time = 6.09249e-07

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
smoothSolver:  Solving for Ux, Initial residual = 4.96493e-05, Final residual = 5.54071e-07, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 5.67593e-05, Final residual = 7.83343e-07, No Iterations 2
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib/x86_64-linux-gnu/libc.so.6
#3  ? in /lib/x86_64-linux-gnu/libm.so.6
#4  ? at ??:?
#5  ? at ??:?
#6  ? in /usr/lib/openfoam/openfoam2206/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam
#7  ? in /lib/x86_64-linux-gnu/libc.so.6
#8  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#9  ? in /usr/lib/openfoam/openfoam2206/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam
Floating point exception (core dumped)

fvSchemes;

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

fluxScheme          Kurganov;

ddtSchemes
{
    default             Euler;
    
}

gradSchemes
{
    default         Gauss linear;

}

divSchemes
{
    
    default         none;
    
    div(tauMC)      Gauss linear;

    div(phi,U)      Gauss Minmod;
    //div(phi,U)      Gauss linearUpwind limitedU;
    //div(phi,U)      Gauss linearUpwindV limitedU;
    //div(phi,U)      Gauss upwind;

    turbulence1      Gauss linear;
    turbulence2      Gauss linearUpwind limitedT;
    turbulence3      Gauss limitedLinear 1;

    div(phi,k)      $turbulence1;
    div(phi,omega)  $turbulence1;

    energy1         Gauss linearUpwind limitedE;
    energy2         Gauss limitedLinear 1;
    energy3         Gauss LUST limitedE;
    energy4         Gauss linear;

    div(phi,e)        $energy4;

//enthalpy -> h
    div(phi,h)      $energy2;
    div(phi,K)      $energy2;

    //div(phid,p)     Gauss upwind;
    //div(phid,p)     Gauss LUST default;
    div(phid,p)     Gauss limitedLinear 1;
    div(phiv,p)        Gauss limitedLinear 1;
    //div(phid,p)     Gauss MUSCL;

    div(((rho*nuEff)*dev2(T(grad(U)))))    Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
    
}

interpolationSchemes
{
    default         linear;
    
    reconstruct(rho) vanLeer;
    reconstruct(U)  vanLeerV;
    reconstruct(T)  vanLeer;
    
}

snGradSchemes
{
    default         corrected;
}

wallDist
{

    method          meshWave;
}

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

fvSolution;

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "(rho|rhoU|rhoE)"
    {
        solver          diagonal;
    }

    "(rho|U|h|k|e|omega)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-06;
        relTol          0.01;
    }

    
}

relaxationFactors
{
    fields
    {
        p             0.3;
        rho           0.1;
    }
    equations
    {
        U                       0.5;
        h                       0.5;
        "(e|k|omega)"     0.5;
    }
}

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

Thank you very much for your valuable time

[Homogenic]

I have also obtained the same set of errors using reactingFOAM for a combustion simulation. Try reducing the time-steps and the maximum Courant number in your solution, that was the problem in my case.

Regards.

[Sakun]

Originally Posted by Homogenic

I have also obtained the same set of errors using reactingFOAM for a combustion simulation. Try reducing the time-steps and the maximum Courant number in your solution, that was the problem in my case.

Regards.

Hello,


Thank you very much for the reply,


I did actually tried those options but not succeeded

[chc]

I don’t know if you’re still working on this, but I think the use of relaxation factors in your fvSolutions may be impacting your results. Since rCF is a transient solver, using relaxation would skew that and possibly induce some errors. Other than that, since your time step is dropping low, another possibility is it could be a result of poor mesh quality or skewed cells.