[berkozaydin]

Hi,


I try to solve a steady, turbulent, incompressible and heat transfer problem. In a rectangular shape, there are 13 batteries. And they have a heat flux value. In all heat transfer models, I have to change the compressible file to incompressible. Therefore, I always make mistakes.

In this study, I used buoyantSimpleFoam. There are some problems while writing the necessary input. And I got this error.

What should I do? Maybe, you can suggest another solver (probably you will say chtMultiRegion) or we can talk about changing compressible file to incompressible (I don't think it's such an easy matter that I can replace density only.)



And also, I want to add thermophysical properties file. Maybe, I have to change something in this file.




Thanks!

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

mixture
{
    specie
    {
        molWeight       28.96;
    }
    thermodynamics
    {
        Cp              1004.4;
        Hf              0;
    }
    transport
    {
        mu              1.831e-05;
        Pr              0.705;
    }
}


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

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

Create mesh for time = 0


SIMPLE: Convergence criteria found
        p_rgh: tolerance 0.001
        U: tolerance 0.0001
        h: tolerance 0.0001
        "(k|epsilon|omega)": tolerance 0.005

Reading thermophysical properties

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

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
    model           kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model eddyDiffusivity

Reading g

Reading hRef
Calculating field g.h


Reading pRef
Reading field p_rgh

No MRF models present

Radiation model not active: radiationProperties not found
Selecting radiationModel none
No finite volume options present

Starting time loop

Time = 1

DILUPBiCGStab:  Solving for h, Initial residual = 1, Final residual = 0.00536049, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 0.999989, Final residual = 0.0037206, No Iterations 7
time step continuity errors : sum local = 0.269408, global = -0.00834916, cumulative = -0.00834916
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.0321705, Final residual = 7.80071e-06, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 1, Final residual = 0.000295899, No Iterations 1
ExecutionTime = 0.46 s  ClockTime = 0 s

Time = 2

DILUPBiCGStab:  Solving for h, Initial residual = 0.026843, Final residual = 0.000173081, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.892057, Final residual = 0.00336125, No Iterations 3
time step continuity errors : sum local = 766.721, global = -4.69339, cumulative = -4.70174
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.740648, Final residual = 0.000894831, No Iterations 1
bounding epsilon, min: -335.513 max: 60702.7 average: 189.468
DILUPBiCGStab:  Solving for k, Initial residual = 0.983617, Final residual = 0.000904046, No Iterations 1
bounding k, min: -6.35689 max: 628.844 average: 10.4528
ExecutionTime = 0.59 s  ClockTime = 0 s

Time = 3

DILUPBiCGStab:  Solving for h, Initial residual = 0.0810137, Final residual = 0.000642639, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.931238, Final residual = 0.00480081, No Iterations 5
time step continuity errors : sum local = 929.728, global = -13.246, cumulative = -17.9478
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.854529, Final residual = 0.000398156, No Iterations 1
bounding epsilon, min: -1834.41 max: 381686 average: 1660.7
DILUPBiCGStab:  Solving for k, Initial residual = 0.63389, Final residual = 0.000277833, No Iterations 1
bounding k, min: -113.523 max: 7561.63 average: 67.972
ExecutionTime = 0.71 s  ClockTime = 1 s

Time = 4

DILUPBiCGStab:  Solving for h, Initial residual = 0.47617, Final residual = 0.00421619, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
#8  Foam::fvMatrix<double>::solve() in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
#9  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
Floating point exception (core dumped)

[HPE]

- try another turbulence model. sometimes epsilon of the kEpsilon model destabilises the computation.