Simulation of a chimney with heat in open atmosphere

Jurado · December 7, 2020, 09:50

Good afternoon,

I am an phd student in machine learning and CFD. I am trying to model a chimney in an neutral open atmosphere with hot areas on top of the chimney. However, I cannot get my head around it and all my trial have unfortunately failed.

I used the buoyantPimpleFoam from openFoam-6 from openFoam.org to try to simulate it. I think the issue reside in the boundary conditions, especially the ones for pressure and p_rgh. I initiate the pressure on the whole domain with a setField using the equation for isotherm atmosphere.

To present more precisely my case, here are some images to present it in the attached files

My simulation domain is a cube, with a face for :

Inlet
Two sides with symmetry
A groud
An outlet
A roof that I have tried as an outlet and symmetry, both failing

Because of the size of this simulation being too big, I have made a simplified case with a rectangular chimney that shows the same issues that I share in the attached files aswell as images of the geometry. there seem to be a problem at some point a vortex start to create itself behind the chimney and then the whole domain start to diverge

Here is also my 0 files:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include        "include/initialConditions"

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (1 0 0);

boundaryField
{

    // inlet air

     inletAir
    {
        type            atmBoundaryLayerInletVelocity;
        #include        "include/ABLConditionsInlet"
    }

    // symmetry

     symmetryRight
    {
        type            symmetry;
    }

     symmetryLeft
    {
        type            symmetry;
    }

    // solid surface

     wallGround
    {
        type            uniformFixedValue;
        uniformValue    (0 0 0);
        value           uniform (0 0 0);
    }

     chimney
    {
        type            uniformFixedValue;
        uniformValue    (0 0 0);
        value           uniform (0 0 0);
    }

     hot_area
    {
        type            uniformFixedValue;
        uniformValue    (0 0 2);
        value           uniform (0 0 2);
    }

    // outlet

     symmetryRoof
    {
        type            symmetry;
//type                freestreamVelocity;    
//freestreamValue    uniform (3.15251 0 0);
    }

     outletPolAir
    {
        type            freestreamVelocity;
    freestreamValue    uniform (1 0 0);
    }

}


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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 293;

boundaryField
{

    // inlet air
    
  inletAir
  {
        type            fixedValue;
    value            uniform 293;
  }

    // symmetry

  symmetryRight
  {
        type            symmetry;
  }

  symmetryLeft
  {
        type            symmetry;
  }

    // solid surface

  wallGround
  {
        type            fixedValue;
    value            uniform 293;
  }

  chimney
  {
          type            zeroGradient;
  }

  hot_area
  {
        type            fixedValue;
    value            uniform 343;
  }

    // outlet

  symmetryRoof
  {

    type        symmetry;
//        type            fixedValue;
//        value            uniform 293;
  }

  outletPolAir
  {
          type            zeroGradient;
  }
}


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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 101325;

boundaryField
{

    // inlet air

    inletAir
    {
        type            zeroGradient;
    }

    // symmetry

    symmetryRight
    {
        type            symmetry;
    }

    symmetryLeft
    {
        type            symmetry;
    }

    // solid surface

    wallGround
    {
        type            zeroGradient;
    }
     chimney
    {
        type            zeroGradient;
    }
     hot_area
    {
        type            zeroGradient;
    }

    // outlet

    symmetryRoof
    {
    type        symmetry;
//      type            fixedValue;
//    value        uniform 101203;    
    }

    outletPolAir
    {
        type            fixedFluxPressure;      
    }

}


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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 101325;

boundaryField
{

    // inlet air

  inletAir
  {
        type            zeroGradient;
  }

    // symmetry

  symmetryRight
  {
      type              symmetry;
  }

  symmetryLeft
  {
      type              symmetry;
  }

    // solid surface

  wallGround
  {
        type            zeroGradient;
  }
  hot_area
  {
        type            zeroGradient;
  }
  chimney
  {
        type            zeroGradient;
  }


    // outlet

  outletPolAir
  {
      type        calculated;
      value        101325;
  }

  symmetryRoof
  {
    type        symmetry;
//      type        calculated;
//      value        95696.491;
  }

}

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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    #include "include/ABLConditionsDomain"

    // symmetry

     symmetryRight
    {
        type            symmetry;
    }
     symmetryLeft
    {
        type            symmetry;
    }
     symmetryRoof
    {
        type            symmetry;
    }
    // solid surface

     wallGround
    {
        type            nutkAtmRoughWallFunction;
        z0              $z0Ground;
        value           uniform 0.0;
    }

     chimney
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

     hot_area
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

    // inlet and outlet

    ".*"
    {
        type            calculated;
        value           uniform 0;
    }
}


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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include        "include/initialConditions"

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform $turbulentKE;

boundaryField
{

    // inlet air

    inletAir
    {
        type            atmBoundaryLayerInletK;
        #include        "include/ABLConditionsInlet"
    }

    // symmetry

     symmetryRight
    {
        type            symmetry;
    }

     symmetryLeft
    {
        type            symmetry;
    }
     symmetryRoof
    {
        type            symmetry;
    }

    // solid surface

     wallGround
    {
        type            kqRWallFunction;
        value           uniform 0.0;
    }

     chimney
    {
        type            kqRWallFunction;
        value           $internalField;
    }

     hot_area
    {
        type            kqRWallFunction;
        value           $internalField;
    }

    // outlet

    ".*"
    {
        type            inletOutlet;
        inletValue      uniform $turbulentKE;
        value           $internalField;
    }
}

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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -3 0 0 0 0];

#include        "include/initialConditions"

internalField   uniform $turbulentEpsilon;

boundaryField
{

    // inlet air

    inletAir
    {
        type            atmBoundaryLayerInletEpsilon;
        #include        "include/ABLConditionsInlet"
    }

    // symmetry

    symmetryRight
    {
        type            symmetry;
    }
    symmetryLeft
    {
        type            symmetry;
    }
     symmetryRoof
    {
        type            symmetry;
    }

    // solid surface

    wallGround
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.4;
        E               9.8;
        value           $internalField;
    }

     chimney
    {
        type            epsilonWallFunction;
        value           $internalField;
    }

     hot_area
    {
        type            epsilonWallFunction;
        value           $internalField;
    }

    // outlet

    ".*"
    {
        type            inletOutlet;
        inletValue      uniform $turbulentEpsilon;
        value           $internalField;
    }
}


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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{

    // inlet air

  inletAir
  {
    type            calculated;
    value           uniform 0;
  }

    // symmetry

  symmetryRight
  {
    type                symmetry;
  }

  symmetryLeft
  {
    type                symmetry;
  }
    // solid surface

  wallGround
  {
    type               compressible::alphatWallFunction;
    value              uniform 0;
  }

  chimney
  {
    type               compressible::alphatWallFunction;
    value              uniform 0;
  }

  hot_area
  {
    type               compressible::alphatWallFunction;
    value              uniform 0;
  }

    // outlet

  symmetryRoof
  {
    type         symmetry;
//    type           calculated;
//    value          uniform 0;
  }

  outletPolAir
  {
    type           calculated;
    value          uniform 0;
  }

}


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

For the includes:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

flowVelocity         (1 0 0);
pressure             0;
turbulentKE          1.3602E-01;
turbulentEpsilon     9.3257E-04;

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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

Uref                 3.15251;
Zref                 300;
zDir                 (0 0 1);
flowDir              (1 0 0);
z0                   0.5;
zGround              0;
value                $internalField;

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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

z0Ground                   uniform 0.1;
z0Block                    uniform 0.01;

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

mAlletto · December 8, 2020, 06:05

can you include an Allrun and Allclean script in your tar file? so it is easier for other people to know what you've done

Jurado · December 8, 2020, 13:48

Quote:

Originally Posted by mAlletto

can you include an Allrun and Allclean script in your tar file? so it is easier for other people to know what you've done

Hi mAlleto,

Thank you for your replay. I added in this message the files with two allrun, one all run to make the mesh and the other allrun to run the computation on 10 cores. I also suppressed a file that were in constant and blocking computation.

It runs on openfoam6.

Jurado · December 11, 2020, 07:08

I have tested a fixed inlet with a speed of 2.3 m/s and another of 1 m/s. Both end up failing with the 1m/s needing way more time to diverge (around 600s for the 1m/s vs 100s for the 2.3 m/s).

Jurado · December 16, 2020, 08:49

Quote:

Originally Posted by Jurado

I have tested a fixed inlet with a speed of 2.3 m/s and another of 1 m/s. Both end up failing with the 1m/s needing way more time to diverge (around 600s for the 1m/s vs 100s for the 2.3 m/s).

I tried to increase the number of loop and to reduce the dt by lowering the CFL to 0.1 but still strange results.

mAlletto · December 16, 2020, 12:50

What OF version are you using. I tried to run it with OF2006 but it failed.

Why did you specify a free stream value of (1 0 0)?

Code:


     outletPolAir
    {
        type            freestreamVelocity;
        freestreamValue uniform (1 0 0);
    }

Code:


    outletPolAir
    {
        type            fixedFluxPressure;
    }

What I see for the pressure boundary condition that you spesify all gradiends (fixedFluxPressure is kind of a specified gradient) did you try to specify a constant value at the outlet?

Jurado · December 17, 2020, 12:19

Quote:

Originally Posted by mAlletto

What OF version are you using. I tried to run it with OF2006 but it failed.

Why did you specify a free stream value of (1 0 0)?

Code:


     outletPolAir
    {
        type            freestreamVelocity;
        freestreamValue uniform (1 0 0);
    }

Code:


    outletPolAir
    {
        type            fixedFluxPressure;
    }

What I see for the pressure boundary condition that you spesify all gradiends (fixedFluxPressure is kind of a specified gradient) did you try to specify a constant value at the outlet?

Hello mAlletto,

Thank you for your reply.

It is strange that it does not work for you, the version with the allrun works perfectly on my PC, probably an issue with openFoam version. The version I use is openfoam 6 from openfoam.org.

I tried several condition for the output pressure like uniformHydrostaticPressure, calculated, totalPressure, freeStreamPressure but none worked. fixedFluxPressure seemed to be the one giving the "best result".

dlahaye · December 18, 2020, 04:16

How does your case differ from the tutorial ./heatTransfer/buoyantBoussinesqPimpleFoam/hotRoom ?

Can you set up a 2D version allowing more comfortable test run in an initial stage of the project?

Does simpleFoam run on your 2D case?

mAlletto · December 18, 2020, 04:28

I saw you had a wall function for espilon at you inlet of the hot air. For k the same.

I suggest you to do the following:

1) Make a 2D flat plate case where you apply youre B.C at the inlet, outlet and symmetry to test if you setup works in this case. If not, than you have something fundamentally wrong.

2) Make a 2D flat plate case with a patch at bottom where you introduce the hot air

3) Make a 2D case with a chimney

4) Now you can go to 3D

Increasing the complexity of your problem step by step helps you to understand the problem better and reduces the amount of errors you make. Furthermore you find errors more quickly. From my experience you get to your final result a lot faster since you avoid time consuming error search and waiting for the simulation to end.

By the way, you do not need a transient solver. Your problem is steady and you can use a steady state solver.

To test your setup use a very small mesh where the simulation runs through in a few seconds/minutes. This speeds up your search for the error a lot. It should run on one cpu not on 10...

Best

Michael

Jurado · December 18, 2020, 08:35

Quote:

Originally Posted by dlahaye

How does your case differ from the tutorial ./heatTransfer/buoyantBoussinesqPimpleFoam/hotRoom ?

Can you set up a 2D version allowing more comfortable test run in an initial stage of the project?

Does simpleFoam run on your 2D case?

Thanks both for the reply

My first trial was with this case boundary condition and it worked when the roof is not high. When the roof start reaching the hundred of meters the issue start appearing. My guess is that when the height is small the difference in pressure between the roof and ground is small. However when we start reaching high heights it is not longer the case and it makes the computation fail.

I will set up a 2D version as mAlleto said to ease the trials. I will come back to post it once I am done.

December 18, 2020, 04:28		#9
mAlletto Senior Member Michael Alletto Join Date: Jun 2018 Location: Bremen Posts: 616 Rep Power: 16	I saw you had a wall function for espilon at you inlet of the hot air. For k the same. I suggest you to do the following: 1) Make a 2D flat plate case where you apply youre B.C at the inlet, outlet and symmetry to test if you setup works in this case. If not, than you have something fundamentally wrong. 2) Make a 2D flat plate case with a patch at bottom where you introduce the hot air 3) Make a 2D case with a chimney 4) Now you can go to 3D Increasing the complexity of your problem step by step helps you to understand the problem better and reduces the amount of errors you make. Furthermore you find errors more quickly. From my experience you get to your final result a lot faster since you avoid time consuming error search and waiting for the simulation to end. By the way, you do not need a transient solver. Your problem is steady and you can use a steady state solver. To test your setup use a very small mesh where the simulation runs through in a few seconds/minutes. This speeds up your search for the error a lot. It should run on one cpu not on 10... Best Michael dlahaye likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Some questions about flow boiling simulation in Fluent	beastieboys6	FLUENT	8	November 21, 2017 00:47
[OpenFOAM.org] Patches to compile OpenFOAM 2.2 on Mac OS X	gschaider	OpenFOAM Installation	136	October 10, 2017 18:25
[swak4Foam] groovyBC in openFOAM-2.0 for parabolic velocity bc	ofslcm	OpenFOAM Community Contributions	25	March 6, 2017 11:03
Error - Solar absorber - Solar Thermal Radiation	MichaelK	CFX	12	September 1, 2016 06:15
SparceImage v1.7.x Issue on MAC OS X	rcarmi	OpenFOAM Installation	4	August 14, 2014 07:42

December 8, 2020, 06:05		#2
mAlletto Senior Member Michael Alletto Join Date: Jun 2018 Location: Bremen Posts: 616 Rep Power: 16	can you include an Allrun and Allclean script in your tar file? so it is easier for other people to know what you've done

December 11, 2020, 07:08		#4
Jurado New Member Jurado Join Date: Nov 2017 Posts: 22 Rep Power: 9	I have tested a fixed inlet with a speed of 2.3 m/s and another of 1 m/s. Both end up failing with the 1m/s needing way more time to diverge (around 600s for the 1m/s vs 100s for the 2.3 m/s).

December 16, 2020, 12:50		#6
mAlletto Senior Member Michael Alletto Join Date: Jun 2018 Location: Bremen Posts: 616 Rep Power: 16	What OF version are you using. I tried to run it with OF2006 but it failed. Why did you specify a free stream value of (1 0 0)? Code: outletPolAir { type freestreamVelocity; freestreamValue uniform (1 0 0); } Code: outletPolAir { type fixedFluxPressure; } What I see for the pressure boundary condition that you spesify all gradiends (fixedFluxPressure is kind of a specified gradient) did you try to specify a constant value at the outlet?

December 18, 2020, 04:16		#8
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 802 Blog Entries: 1 Rep Power: 18	How does your case differ from the tutorial ./heatTransfer/buoyantBoussinesqPimpleFoam/hotRoom ? Can you set up a 2D version allowing more comfortable test run in an initial stage of the project? Does simpleFoam run on your 2D case?