[oclerc]

Hi all,

I am currently stuck on the same issue for a long time. I have used the circuitBoardColling as an initial template to solve a simple heat transfer problem where I have an inlet and outlet flow and a box in the middle generating heat at a constant temperature. However, the solving always crashes at time=2:

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 9-a0f1846504f2
Exec   : buoyantSimpleFoam
Date   : Oct 20 2023
Time   : 16:46:10
Host   : "oclerc-heatsink-0-0"
PID    : 2490
I/O    : uncollated
Case   : /cvlabdata2/home/oclerc/iglooWithFridges
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 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 unityLewisEddyDiffusivity

Reading g

Reading hRef
Calculating field g.h


Reading pRef
Reading field p_rgh

No MRF models present

No fvModels present
No fvConstraints present

Starting time loop

Time = 1

DILUPBiCGStab:  Solving for h, Initial residual = 0.000295657, Final residual = 2.74438e-06, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 0.999975, Final residual = 0.00201398, No Iterations 1
time step continuity errors : sum local = 5.78159, global = 3.59638, cumulative = 3.59638
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.98999, Final residual = 0.00561489, No Iterations 1
bounding epsilon, min: -111.168 max: 2354.29 average: 9.78446
DILUPBiCGStab:  Solving for k, Initial residual = 1, Final residual = 1.75555e-05, No Iterations 1
bounding k, min: -4.26375 max: 108.425 average: 0.702237
ExecutionTime = 0.27 s  ClockTime = 0 s

Time = 2

DILUPBiCGStab:  Solving for h, Initial residual = 1, Final residual = 0.00111084, No Iterations 2
#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/openfoam9/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
#8  Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
#9  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"

Does anyone have an insight on what could be causing this crash? I was guessing it could be boundary conditions. You can find the folder to run attached.

example_cfd.zip

[Yann]

Hello Olivier,

This is indeed a matter of boundary conditions: in p_rgh, you set the pressure to 0. Switch it to 1e5 and it runs.

Cheers,
Yann

[oclerc]

Hey Yann,

Thanks a lot, it did indeed solve the issue.

Cheers,
Olivier