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

Data center cfd /w OpenFOAM boundary conditions?

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 19, 2016, 09:57
Default Data center cfd /w OpenFOAM boundary conditions?
  #1
New Member
 
Yannic E.
Join Date: Nov 2016
Posts: 6
Rep Power: 10
jaroz is on a distinguished road
Hello

I am trying my luck here, so I need to simulate the airflow in a data center and have issues defining boundary conditions.
First of all, floor tiles release air with constant mass flow rate (at pressure p+delta p) and a certain temperature setpoint T into the room. This seems ok

Next, there are servers, which host an internal fan, and are also a volumetric heat source.
I would need for the air to enter the server, which is defined as a Cell zone, and add a velocity source term for the fan, e.g. (50,0.0) , as well as a heat source term in W.
I am using buoyantBoussinesqSimpleFoam solver, which can only define temperature source term, so how to define heat source?

- Do I need to define cyclic boundaries on the inlet and Outlet of the server? But then i think Iit is not possible to model the heat transfer.

Finally, is it necessary to use chtmultiregionfoam for this problem and define the servers as solid regions? There it is possible to model conjugate heat transfer with source term h.
jaroz is offline   Reply With Quote

Old   December 19, 2016, 17:02
Default Use the scalarSemiImplicitSource type in fvOptions
  #2
Member
 
Arvind Jay
Join Date: Sep 2012
Posts: 97
Rep Power: 15
arvindpj is on a distinguished road
In order to specify volumetric heat source you could use the scalarSemiImplicitSource type in fvOptions as below:

Code:
heatSource
{
    type            scalarSemiImplicitSource;
    active          true;
 
    scalarSemiImplicitSourceCoeffs
    {
        selectionMode   all; // all, cellSet, cellZone, points
        cellSet         c1;
        volumeMode      specific; // absolute;
        injectionRateSuSp
        {
            T     (0.1 0);
        }
    }
}
Please refer this blog for more info.

Cheers
arvindpj is offline   Reply With Quote

Old   December 20, 2016, 06:03
Smile
  #3
New Member
 
Yannic E.
Join Date: Nov 2016
Posts: 6
Rep Power: 10
jaroz is on a distinguished road
Hello

Thanks very much for the post, especially the linked blog was very helpful for this!

If you don't mind I have developed some follow up questions:

1. So it seems for the Heat source term SC[K/s]:=q˙v/ρc. This means I can find out S_C by dividing the Total Heat (in W) by the density and capacity of the medium?

2. Is it possible to add Velocity source term to the cell zone, like described in
http://caefn.com/openfoam/fvoptions-meanvelocityforce
In addition to the heat source term?

For the Boussinesq simplification, that is described here, my Temperatures will vary approx. from 20°C-60°C (for hot spots). Is the Boussinesq approach okay here or Delta T too large?

Thanks
jaroz is offline   Reply With Quote

Old   December 20, 2016, 18:05
Default
  #4
Member
 
Arvind Jay
Join Date: Sep 2012
Posts: 97
Rep Power: 15
arvindpj is on a distinguished road
Hello

Thanks very much for the post, especially the linked blog was very helpful for this!

If you don't mind I have developed some follow up questions:

1. So it seems for the Heat source term SC[K/s]:=q˙v/ρc. This means I can find out S_C by dividing the Total Heat (in W) by the density and capacity of the medium?

yes

2. Is it possible to add Velocity source term to the cell zone, like described in
http://caefn.com/openfoam/fvoptions-meanvelocityforce
In addition to the heat source term?

yes, and have all in the same fvoption file.

For the Boussinesq simplification, that is described here, my Temperatures will vary approx. from 20°C-60°C (for hot spots). Is the Boussinesq approach okay here or Delta T too large?

Boussinesq simplification is valid if beta(T-T0) < 1; beta = thermal expansion coefficient; T0 = operating Temp (K)

Last edited by arvindpj; December 21, 2016 at 11:33.
arvindpj is offline   Reply With Quote

Old   January 29, 2017, 12:00
Default
  #5
New Member
 
Yannic E.
Join Date: Nov 2016
Posts: 6
Rep Power: 10
jaroz is on a distinguished road
As a follow up question, my model is developed and solving nicely.
The issue I am facing is that the heat distribution in the medium doesn't seem to be homogeneous, but rather following a gauss-like curve with the following fvOptions code:


Code:
heatSource
    {
        type scalarSemiImplicitSource;
        active true;
        scalarSemiImplicitSourceCoeffs
        {
            selectionMode cellZone;
            cellZone servers;
            volumeMode specific;
            injectionRateSuSp
            {
                T ( 37.0 0);
            }

        }

    }




Higher heat load is clearly in the middle of the geometry. Is it possible to generate a uniform heat source term for the cell zone, where the center and boundaries are equally hot?

Cheers

Last edited by jaroz; January 29, 2017 at 15:37.
jaroz is offline   Reply With Quote

Old   February 21, 2017, 05:35
Default
  #6
New Member
 
Yannic E.
Join Date: Nov 2016
Posts: 6
Rep Power: 10
jaroz is on a distinguished road
Sorry to bring up the post again, but I am still facing some issues. I have modeled the room successfully, but heat distribution remains uneven.

I have since figured out the source for this, air is recirculating from the back of the servers back into the server, this results in heating up.
The internal fan is modeled as a momentum source, is there a way to model the cell zone to deny backflow (with a boundary condition on the outflow side?)
jaroz is offline   Reply With Quote

Old   October 8, 2021, 08:05
Default
  #7
Member
 
Paulo
Join Date: Jun 2011
Posts: 34
Rep Power: 15
strobel is on a distinguished road
Quote:
Originally Posted by jaroz View Post
Sorry to bring up the post again, but I am still facing some issues. I have modeled the room successfully, but heat distribution remains uneven.

I have since figured out the source for this, air is recirculating from the back of the servers back into the server, this results in heating up.
The internal fan is modeled as a momentum source, is there a way to model the cell zone to deny backflow (with a boundary condition on the outflow side?)
Hi, have you solved this issue? I usually try porous media to avoid this recirculation, but I don't know how this would fit on your modelling.
strobel is offline   Reply With Quote

Reply


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
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 08:30
Some Problems about the Boundary Conditions in OpenFoam lzgwhy OpenFOAM Running, Solving & CFD 47 October 10, 2017 09:33
Error - Solar absorber - Solar Thermal Radiation MichaelK CFX 12 September 1, 2016 06:15
Implementation of boundary conditions for FVM Tom Main CFD Forum 7 August 26, 2014 06:58
New topic on same subject - Flow around race car Tudor Miron CFX 15 April 2, 2004 07:18


All times are GMT -4. The time now is 13:59.