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

buoyantSimpleFoam Error

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

Like Tree1Likes
  • 1 Post By HPE

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 22, 2021, 07:48
Default buoyantSimpleFoam Error
  #1
New Member
 
Ataberk
Join Date: Nov 2019
Posts: 14
Rep Power: 4
berkozaydin is on a distinguished road
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!

Code:
/*--------------------------------*- 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;
    }
}


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





Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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)
berkozaydin is offline   Reply With Quote

Old   June 22, 2021, 08:17
Default
  #2
HPE
Senior Member
 
HPE's Avatar
 
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 10
HPE is on a distinguished road
- try another turbulence model. sometimes epsilon of the kEpsilon model destabilises the computation.
saidc. likes this.
HPE is offline   Reply With Quote

Reply

Tags
buoyantsimplefoam, chtmultiregionfoam, heat transfer, openfoam

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
DPM udf error haghshenasfard FLUENT 0 April 13, 2016 06:35
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh gschaider OpenFOAM Community Contributions 300 October 29, 2014 18:00
OpenFOAM without MPI kokizzu OpenFOAM Installation 4 May 26, 2014 09:17
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 09:31
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50


All times are GMT -4. The time now is 01:46.