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

buoyantSimpleFoam crashes at Time 2

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

Like Tree1Likes
  • 1 Post By oclerc

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 20, 2023, 13:58
Default buoyantSimpleFoam crashes at Time 2
  #1
New Member
 
Olivier
Join Date: Oct 2023
Posts: 2
Rep Power: 0
oclerc is on a distinguished road
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:

Code:
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  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
oclerc is offline   Reply With Quote

Old   October 21, 2023, 05:59
Default
  #2
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,001
Rep Power: 25
Yann will become famous soon enough
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
Yann is offline   Reply With Quote

Old   October 21, 2023, 15:12
Default
  #3
New Member
 
Olivier
Join Date: Oct 2023
Posts: 2
Rep Power: 0
oclerc is on a distinguished road
Hey Yann,

Thanks a lot, it did indeed solve the issue.

Cheers,
Olivier
Yann likes this.
oclerc is offline   Reply With Quote

Reply

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
[Other] Contribution a new utility: refine wall layer mesh based on yPlus field lakeat OpenFOAM Community Contributions 58 December 23, 2021 03:36
Slow calculation time on CFD server killian153 Hardware 4 July 27, 2021 08:28
Inconsistencies in reading .dat file during run time in new injection model Scram_1 OpenFOAM 0 March 23, 2018 23:29
pimpleDyMFoam computation randomly stops babapeti OpenFOAM Running, Solving & CFD 5 January 24, 2018 06:28
a problem with convergence in buoyantSimpleFoam skuznet OpenFOAM Running, Solving & CFD 6 November 15, 2017 13:12


All times are GMT -4. The time now is 22:28.