|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Olivier
Join Date: Oct 2023
Posts: 2
Rep Power: 0 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,169
Rep Power: 27 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Olivier
Join Date: Oct 2023
Posts: 2
Rep Power: 0 ![]() |
Hey Yann,
Thanks a lot, it did indeed solve the issue. ![]() Cheers, Olivier |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
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 02:36 |
Slow calculation time on CFD server | killian153 | Hardware | 4 | July 27, 2021 07:28 |
Inconsistencies in reading .dat file during run time in new injection model | Scram_1 | OpenFOAM | 0 | March 23, 2018 22:29 |
pimpleDyMFoam computation randomly stops | babapeti | OpenFOAM Running, Solving & CFD | 5 | January 24, 2018 05:28 |
a problem with convergence in buoyantSimpleFoam | skuznet | OpenFOAM Running, Solving & CFD | 6 | November 15, 2017 12:12 |