|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Alex Pressley
Join Date: Feb 2019
Posts: 3
Rep Power: 7 ![]() |
Hi all,
I have been using buoyantBoussinesqPimpleFoam for some time now quite happily in OpenFoam 6 and earlier. However, I would like to be able to convert to version 7, where the boussinesq solver has been integrated into the general buoyantPimpleFoam solver (I think it should be easier to simulate a heat sink of specific power in this solver). Unfortunately, I have not been able to run a model that even comes close to working other than the hotRoomBoussinesq tutorial. Has anybody else had more success? The model of mine that I am trying to run wildly diverges, especially in T results, despite there being no heat source or sink, just a momentum source to simulate a pump. Originally, it would do one or two time steps then in the next exceed the number of permitted iterations and spit out an error like Code:
--> FOAM FATAL ERROR: Negative initial temperature T0: -[large number] Code:
#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib64/libc.so.6" #3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) at ??:? #4 Foam::DILUPreconditioner::DILUPreconditioner(Foam::lduMatrix::solver const&, Foam::dictionary const&) at ??:? #5 Foam::lduMatrix::preconditioner::addasymMatrixConstructorToTable<Foam::DILUPreconditioner>::New(Foam::lduMatrix::solver const&, Foam::dictionary const&) at ??:? #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatrix::solver const&, Foam::dictionary const&) at ??:? #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #8 ? at ??:? #9 ? at ??:? #10 ? at ??:? #11 ? at ??:? #12 __libc_start_main in "/lib64/libc.so.6" #13 ? at ??:? /opt/site/sge/default/spool/node405/job_scripts/1039376: line 9: 80501 Floating point exceptionbuoyantPimpleFoam I would be very grateful for any help from someone who has had success with Boussinesq calculations in OpenFOAM 7. Thank you, Alex |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
Join Date: Nov 2020
Posts: 53
Rep Power: 5 ![]() |
Hey, maybe try to put an fvOptions limitTemperature and see if the solver works well.
|
|
![]() |
![]() |
![]() |
Tags |
boussinesq, buoyantboussinesqpimple, buoyantpimplefoam, openfoam7 |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Frequently Asked Questions about Installing OpenFOAM | wyldckat | OpenFOAM Installation | 3 | November 14, 2023 11:58 |
boussinesq approximation when temperature difference is large | a.lone | FLUENT | 3 | September 15, 2020 15:07 |
Air without Boussinesq approximation | Sidharthkp | CFX | 2 | February 14, 2020 04:32 |
Boussinesq approximation | Abdul Aziz Jaafar | Main CFD Forum | 19 | May 3, 2019 11:00 |
How to contribute to the community of OpenFOAM users and to the OpenFOAM technology | wyldckat | OpenFOAM | 17 | November 10, 2017 15:54 |